The toxicity data landscape for environmental chemicals

An automated computational data pipeline to rapidly acquire, score, and rank toxicological data for ecological hazard assessment

Biological Evaluations support Endangered Species Act (ESA) consultation with the US Fish and Wildlife Service and National Marine Fisheries Service by federal action agencies, such as the USEPA, regarding impacts of federal activities on threatened or endangered species. However, they are often time-consuming and challenging to conduct. The identification of pollutant benchmarks or guidance to protect taxa for states and tribes when USEPA has not yet developed criteria recommendations is also of importance to ensure a streamlined approach to Clean Water Act program implementation. Due to substantial workloads, tight regulatory timelines, and the often-protracted length of ESA consultations, there is a need to streamline the development of biological evaluation toxicity assessments for determining the impact of chemical pollutants on ESA-listed species. Moreover, there is limited availability of species-specific toxicity data for many contaminants, further complicating the consultation process. New approach methodologies are being increasingly used in toxicology and chemical safety assessment to rapidly and cost-effectively provide data that can fill gaps in hazard and/or exposure characterization. Here, we present the development of an automated computational pipeline-RASRTox (Rapidly Acquire, Score, and Rank Toxicological data)-to rapidly extract and categorize ecological toxicity benchmark values from curated data sources (ECOTOX, ToxCast) and well-established quantitative structure-activity relationships (TEST, ECOSAR). As a proof of concept, points-of-departure (PODs) generated in RASRTox for 13 chemicals were compared against benchmark values derived using traditional methods-toxicity reference values (TRVs) and water quality criteria (WQC). The RASRTox PODs were generally within an order of magnitude of corresponding TRVs, though less concordant compared with WQC. The greatest utility of RASRTox, however, lies in its ability to quickly and systematically identify critical studies that may serve as a basis for screening value derivation by toxicologists as part of an ecological hazard assessment. As such, the strategy described in this case study can potentially be adapted for other risk assessment contexts and stakeholder needs. Integr Environ Assess Manag 2024;20:2203-2217. © 2024 Society of Environmental Toxicology & Chemistry (SETAC). This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

Incorporating new approach methods (NAMs) data in dose-response assessments: The future is now!

Regulatory dose-response assessments traditionally rely on in vivo data and default assumptions. New Approach Methods (NAMs) present considerable opportunities to both augment traditional dose-response assessments and accelerate the evaluation of new/data-poor chemicals. This review aimed to determine the potential utilization of NAMs through a unified conceptual framework that compartmentalizes derivation of toxicity values into five sequential Key Dose-response Modules (KDMs): (1) point-of-departure (POD) determination, (2) test system-to-human (e.g. inter-species) toxicokinetics and (3) toxicodynamics, (4) human population (intra-species) variability in toxicodynamics, and (5) toxicokinetics. After using several "traditional" dose-response assessments to illustrate this framework, a review is presented where existing NAMs, including in silico, in vitro, and in vivo approaches, might be applied across KDMs. Further, the false dichotomy between "traditional" and NAMs-derived data sources is broken down by organizing dose-response assessments into a matrix where each KDM has Tiers of increasing precision and confidence: Tier 0: Default/generic values, Tier 1: Computational predictions, Tier 2: Surrogate measurements, and Tier 3: Direct measurements. These findings demonstrated that although many publications promote the use of NAMs in KDMs (1) for POD determination and (5) for human population toxicokinetics, the proposed matrix of KDMs and Tiers reveals additional immediate opportunities for NAMs to be integrated across other KDMs. Further, critical needs were identified for developing NAMs to improve in vitro dosimetry and quantify test system and human population toxicodynamics. Overall, broadening the integration of NAMs across the steps of dose-response assessment promises to yield higher throughput, less animal-dependent, and more science-based toxicity values for protecting human health.

High-Throughput Transcriptomics Screen of ToxCast Chemicals in U-2 OS Cells

New approach methodologies (NAMs) aim to accelerate the pace of chemical risk assessment while simultaneously reducing cost and dependency on animal studies. High Throughput Transcriptomics (HTTr) is an emerging NAM in the field of chemical hazard evaluation for establishing in vitro points-of-departure and providing mechanistic insight. In the current study, 1201 test chemicals were screened for bioactivity at eight concentrations using a 24-h exposure duration in the human- derived U-2 OS osteosarcoma cell line with HTTr. Assay reproducibility was assessed using three reference chemicals that were screened on every assay plate. The resulting transcriptomics data were analyzed by aggregating signal from genes into signature scores using gene set enrichment analysis, followed by concentration-response modeling of signatures scores. Signature scores were used to predict putative mechanisms of action, and to identify biological pathway altering concentrations (BPACs). BPACs were consistent across replicates for each reference chemical, with replicate BPAC standard deviations as low as 5.6 × 10-3 μM, demonstrating the internal reproducibility of HTTr-derived potency estimates. BPACs of test chemicals showed modest agreement (R2 = 0.55) with existing phenotype altering concentrations from high throughput phenotypic profiling using Cell Painting of the same chemicals in the same cell line. Altogether, this HTTr based chemical screen contributes to an accumulating pool of publicly available transcriptomic data relevant for chemical hazard evaluation and reinforces the utility of cell based molecular profiling methods in estimating chemical potency and predicting mechanism of action across a diverse set of chemicals.

Food for thought- Paving the way for a UK Roadmap towards optimum consumer safety: Development, Endorsement and Regulatory Acceptance of New Approach Methodologies (NAMs) in Chemical Risk Assessment and Beyond

Advances in biosciences, chemistry, technology, and computer sciences have resulted in the unparalleled development of candidate New Approach Methodologies over the last few years. Many of these are potentially invaluable in the safety assessment of chemicals, but very few have been adopted for regulatory decision making. There is an immediate opportunity to use NAMs in safety assessment where the vision is to be able to predict risk more rapidly, accurately, and efficiently to further assure consumer safety. In order to achieve this, the UK Food Standards Agency (FSA) and the Committee on Toxicity of Chemicals in Food, Consumer Products and the Environment (COT) have developed a roadmap towards acceptance and integration of these new approach methodologies into safety and risk assessments for regulatory decision making. The roadmap provides a UK blueprint for the transition of NAMs from the research laboratory to their use in regulatory decision making. This will require close collaboration across disciplines (chemists, toxicologists, informaticians, risk assessors and others), and across chemical sectors, to develop, verify and utilise appropriate models. Linking up internationally, and harmonization will be fundamental.

In vitro to in vivo extrapolation from 3D hiPSC-derived cardiac microtissues and physiologically based pharmacokinetic modeling to inform next-generation arrhythmia risk assessment

Proarrhythmic cardiotoxicity remains a substantial barrier to drug development as well as a major global health challenge. In vitro human pluripotent stem cell-based new approach methodologies have been increasingly proposed and employed as alternatives to existing in vitro and in vivo models that do not accurately recapitulate human cardiac electrophysiology or cardiotoxicity risk. In this study, we expanded the capacity of our previously established 3D human cardiac microtissue model to perform quantitative risk assessment by combining it with a physiologically based pharmacokinetic model, allowing a direct comparison of potentially harmful concentrations predicted in vitro to in vivo therapeutic levels. This approach enabled the measurement of concentration responses and margins of exposure for 2 physiologically relevant metrics of proarrhythmic risk (i.e. action potential duration and triangulation assessed by optical mapping) across concentrations spanning 3 orders of magnitude. The combination of both metrics enabled accurate proarrhythmic risk assessment of 4 compounds with a range of known proarrhythmic risk profiles (i.e. quinidine, cisapride, ranolazine, and verapamil) and demonstrated close agreement with their known clinical effects. Action potential triangulation was found to be a more sensitive metric for predicting proarrhythmic risk associated with the primary mechanism of concern for pharmaceutical-induced fatal ventricular arrhythmias, delayed cardiac repolarization due to inhibition of the rapid delayed rectifier potassium channel, or hERG channel. This study advances human-induced pluripotent stem cell-based 3D cardiac tissue models as new approach methodologies that enable in vitro proarrhythmic risk assessment with high precision of quantitative metrics for understanding clinically relevant cardiotoxicity.

Identification of neural-relevant toxcast high-throughput assay intended gene targets: Applicability to neurotoxicity and neurotoxicant putative molecular initiating events

The US EPA's Toxicity Forecaster (ToxCast) is a suite of high-throughput in vitro assays to screen environmental toxicants and predict potential toxicity of uncharacterized chemicals. This work examines the relevance of ToxCast assay intended gene targets to putative molecular initiating events (MIEs) of neurotoxicants. This effort is needed as there is growing interest in the regulatory and scientific communities about developing new approach methodologies (NAMs) to screen large numbers of chemicals for neurotoxicity and developmental neurotoxicity. Assay gene function (GeneCards, NCBI-PUBMED) was used to categorize gene target neural relevance (1 = neural, 2 = neural development, 3 = general cellular process, 3 A = cellular process critical during neural development, 4 = unlikely significance). Of 481 unique gene targets, 80 = category 1 (16.6 %); 16 = category 2 (3.3 %); 303 = category 3 (63.0 %); 97 = category 3 A (20.2 %); 82 = category 4 (17.0 %). A representative list of neurotoxicants (548) was researched (ex. PUBMED, PubChem) for neurotoxicity associated MIEs/Key Events (KEs). MIEs were identified for 375 compounds, whereas only KEs for 173. ToxCast gene targets associated with MIEs were primarily neurotransmitter (ex. dopaminergic, GABA)receptors and ion channels (calcium, sodium, potassium). Conversely, numerous MIEs associated with neurotoxicity were absent. Oxidative stress (OS) mechanisms were 79.1 % of KEs. In summary, 40 % of ToxCast assay gene targets are relevant to neurotoxicity mechanisms. Additional receptor and ion channel subtypes and increased OS pathway coverage are identified for potential future assay inclusion to provide more complete coverage of neural and developmental neural targets in assessing neurotoxicity.

Integrated Approach for Testing and Assessment for Developmental Neurotoxicity (DNT) to Prioritize Aromatic Organophosphorus Flame Retardants

Organophosphorus flame retardants (OPFRs) are abundant and persistent in the environment but have limited toxicity information. Their similarity in structure to organophosphate pesticides presents great concern for developmental neurotoxicity (DNT). However, current in vivo testing is not suitable to provide DNT information on the amount of OPFRs that lack data. Over the past decade, an in vitro battery was developed to enhance DNT assessment, consisting of assays that evaluate cellular processes in neurodevelopment and function. In this study, behavioral data of small model organisms were also included. To assess if these assays provide sufficient mechanistic coverage to prioritize chemicals for further testing and/or identify hazards, an integrated approach to testing and assessment (IATA) was developed with additional information from the Integrated Chemical Environment (ICE) and the literature. Human biomonitoring and exposure data were identified and physiologically-based toxicokinetic models were applied to relate in vitro toxicity data to human exposure based on maximum plasma concentration. Eight OPFRs were evaluated, including aromatic OPFRs (triphenyl phosphate (TPHP), isopropylated phenyl phosphate (IPP), 2-ethylhexyl diphenyl phosphate (EHDP), tricresyl phosphate (TMPP), isodecyl diphenyl phosphate (IDDP), tert-butylphenyl diphenyl phosphate (BPDP)) and halogenated FRs ((Tris(1,3-dichloro-2-propyl) phosphate (TDCIPP), tris(2-chloroethyl) phosphate (TCEP)). Two representative brominated flame retardants (BFRs) (2,2'4,4'-tetrabromodiphenyl ether (BDE-47) and 3,3',5,5'-tetrabromobisphenol A (TBBPA)) with known DNT potential were selected for toxicity benchmarking. Data from the DNT battery indicate that the aromatic OPFRs have activity at similar concentrations as the BFRs and should therefore be evaluated further. However, these assays provide limited information on the mechanism of the compounds. By integrating information from ICE and the literature, endocrine disruption was identified as a potential mechanism. This IATA case study indicates that human exposure to some OPFRs could lead to a plasma concentration similar to those exerting in vitro activities, indicating potential concern for human health.

Identifying uncertainty in physical-chemical property estimation with IFSQSAR

This study describes the development and evaluation of six new models for predicting physical-chemical (PC) properties that are highly relevant for chemical hazard, exposure, and risk estimation: solubility (in water SW and octanol SO), vapor pressure (VP), and the octanol-water (KOW), octanol-air (KOA), and air-water (KAW) partition ratios. The models are implemented in the Iterative Fragment Selection Quantitative Structure-Activity Relationship (IFSQSAR) python package, Version 1.1.0. These models are implemented as Poly-Parameter Linear Free Energy Relationship (PPLFER) equations which combine experimentally calibrated system parameters and solute descriptors predicted with QSPRs. Two other ancillary models have been developed and implemented, a QSPR for Molar Volume (MV) and a classifier for the physical state of chemicals at room temperature. The IFSQSAR methods for characterizing applicability domain (AD) and calculating uncertainty estimates expressed as 95% prediction intervals (PI) for predicted properties are described and tested on 9,000 measured partition ratios and 4,000 VP and SW values. The measured data are external to IFSQSAR training and validation datasets and are used to assess the predictivity of the models for "novel chemicals" in an unbiased manner. The 95% PI intervals calculated from validation datasets for partition ratios needed to be scaled by a factor of 1.25 to capture 95% of the external data. Predictions for VP and SW are more uncertain, primarily due to the challenges in differentiating their physical state (i.e., liquids or solids) at room temperature. The prediction accuracy of the models for log KOW, log KAW and log KOA of novel, data-poor chemicals is estimated to be in the range of 0.7 to 1.4 root mean squared error of prediction (RMSEP), with RMSEP in the range 1.7-1.8 for log VP and log SW. Scientific contributionNew partitioning models integrate empirical PPLFER equations and QSARs, allowing for seamless integration of experimental data and model predictions. This work tests the real predictivity of the models for novel chemicals which are not in the model training or external validation datasets.

Probabilistic Reference and 10% Effect Concentrations for Characterizing Inhalation Non-cancer and Developmental/Reproductive Effects for 2,160 Substances

Chemicals assessment and management frameworks rely on regulatory toxicity values, which are based on points of departure (POD) identified following rigorous dose-response assessments. Yet, regulatory PODs and toxicity values for inhalation exposure (i.e., reference concentrations [RfCs]) are available for only ∼200 chemicals. To address this gap, we applied a workflow to determine surrogate inhalation route PODs and corresponding toxicity values, where regulatory assessments are lacking. We curated and selected inhalation in vivo data from the U.S. EPA's ToxValDB and adjusted reported effect values to chronic human equivalent benchmark concentrations (BMCh) following the WHO/IPCS framework. Using ToxValDB chemicals with existing PODs associated with regulatory toxicity values, we found that the 25th %-ile of a chemical's BMCh distribution (POD p 25 BMC h ) could serve as a suitable surrogate for regulatory PODs (Q2 ≥ 0.76, RSE ≤ 0.82 log10 units). We applied this approach to derive POD p 25 BMC h for 2,095 substances with general non-cancer toxicity effects and 638 substances with reproductive/developmental toxicity effects, yielding a total coverage of 2,160 substances. From these POD p 25 BMC h , we derived probabilistic RfCs and human population effect concentrations. With this work, we have expanded the number of chemicals with toxicity values available, thereby enabling a much broader coverage for inhalation risk and impact assessment.

Mouse models to explore the biological and organismic role of DNA polymerase beta

Gene knock-out (KO) mouse models for DNA polymerase beta (Polβ) revealed that loss of Polβ leads to neonatal lethality, highlighting the critical organismic role for this DNA polymerase. While biochemical analysis and gene KO cell lines have confirmed its biochemical role in base excision repair and in TET-mediated demethylation, more long-lived mouse models continue to be developed to further define its organismic role. The Polb-KO mouse was the first of the Cre-mediated tissue-specific KO mouse models. This technology was exploited to investigate roles for Polβ in V(D)J recombination (variable-diversity-joining rearrangement), DNA demethylation, gene complementation, SPO11-induced DNA double-strand break repair, germ cell genome stability, as well as neuronal differentiation, susceptibility to genotoxin-induced DNA damage, and cancer onset. The revolution in knock-in (KI) mouse models was made possible by CRISPR/cas9-mediated gene editing directly in C57BL/6 zygotes. This technology has helped identify phenotypes associated with germline or somatic mutants of Polβ. Such KI mouse models have helped uncover the importance of key Polβ active site residues or specific Polβ enzyme activities, such as the PolbY265C mouse that develops lupus symptoms. More recently, we have used this KI technology to mutate the Polb gene with two codon changes, yielding the PolbL301R/V303R mouse. In this KI mouse model, the expressed Polβ protein cannot bind to its obligate heterodimer partner, Xrcc1. Although the expressed mutant Polβ protein is proteolytically unstable and defective in recruitment to sites of DNA damage, the homozygous PolbL301R/V303R mouse is viable and fertile, yet small in stature. We expect that this and additional targeted mouse models under development are poised to reveal new biological and organismic roles for Polβ.

Combination of computational new approach methodologies for enhancing evidence of biological pathway conservation across species

The ability to predict which chemicals are of concern for environmental safety is dependent, in part, on the ability to extrapolate chemical effects across many species. This work investigated the complementary use of two computational new approach methodologies to support cross-species predictions of chemical susceptibility: the US Environmental Protection Agency Sequence Alignment to Predict Across Species Susceptibility (SeqAPASS) tool and Unilever's recently developed Genes to Pathways - Species Conservation Analysis (G2P-SCAN) tool. These stand-alone tools rely on existing biological knowledge to help understand chemical susceptibility and biological pathway conservation across species. The utility and challenges of these combined computational approaches were demonstrated using case examples focused on chemical interactions with peroxisome proliferator activated receptor alpha (PPARα), estrogen receptor 1 (ESR1), and gamma-aminobutyric acid type A receptor subunit alpha (GABRA1). Overall, the biological pathway information enhanced the weight of evidence to support cross-species susceptibility predictions. Through comparisons of relevant molecular and functional data gleaned from adverse outcome pathways (AOPs) to mapped biological pathways, it was possible to gain a toxicological context for various chemical-protein interactions. The information gained through this computational approach could ultimately inform chemical safety assessments by enhancing cross-species predictions of chemical susceptibility. It could also help fulfill a core objective of the AOP framework by potentially expanding the biologically plausible taxonomic domain of applicability of relevant AOPs.

Interpretable attention-based multi-encoder transformer based QSPR model for assessing toxicity and environmental impact of chemicals

The rising demand from consumer goods and pharmaceutical industry is driving a fast expansion of newly developed chemicals. The conventional toxicity testing of unknown chemicals is expensive, time-consuming, and raises ethical concerns. The quantitative structure-property relationship (QSPR) is an efficient computational method because it saves time, resources, and animal experimentation. Advances in machine learning have improved chemical analysis in QSPR studies, but the real-world application of machine learning-based QSPR studies was limited by the unexplainable 'black box' feature of the machine learnings. In this study, multi-encoder structure-to-toxicity (S2T)-transformer based QSPR model was developed to estimate the properties of polychlorinated biphenyls (PCBs) and endocrine disrupting chemicals (EDCs). Simplified molecular input line entry systems (SMILES) and molecular descriptors calculated by the Dragon 6 software, were simultaneously considered as input of QSPR model. Furthermore, an attention-based framework is proposed to describe the relationship between the molecular structure and toxicity of hazardous chemicals. The S2T-transformer model achieved the highest R2 scores of 0.918, 0.856, and 0.907 for logarithm of octanol-water partition coefficient (Log KOW), octanol-air partition coefficient (Log KOA), and bioconcentration factor (Log BCF) estimation of PCBs, respectively. Moreover, the attention weights were able to properly interpret the lateral (meta, para) chlorination associated with PCBs toxicity and environmental impact.

Changes in chemical occurrence, concentration, and bioactivity in the Colorado River before and after replacement of the Moab, Utah wastewater treatment plant

Long-term (2010-19) water-quality monitoring on the Colorado River downstream from Moab Utah indicated the persistent presence of Bioactive Chemicals (BC), such as pesticides and pharmaceuticals. This stream reach near Canyonlands National Park provides critical habitat for federally endangered species. The Moab wastewater treatment plant (WWTP) outfall discharges to the Colorado River and is the nearest potential point-source to this reach. The original WWTP was replaced in 2018. In 2016-19, a study was completed to determine if the new plant reduced BC input to the Colorado River at, and downstream from, the outfall. Water samples were collected before and after the plant replacement at sites upstream and downstream from the outfall. Samples were analyzed for as many as 243 pesticides, 109 pharmaceuticals, 20 hormones, 51 wastewater indicator chemicals, 20 metals, and 8 nutrients. BC concentrations, hazard quotients (HQs), and exposure activity ratios (EARs) were used to identify and prioritize contaminants for their potential to have adverse biological effects on the health of native and endangered wildlife. There were 22 BC with HQs >1, mostly metals and hormones; and 23 BC with EARs >0.1, mostly hormones and pharmaceuticals. Most high HQs or EARs were associated with samples collected at the WWTP outfall site prior to its replacement. Discharge from the new plant had reduced concentrations of nutrients, hormones, pharmaceuticals, and other BC. For example, all 16 of the hormones detected at the WWTP outfall site had maximum concentrations in samples collected prior to the WWTP replacement. The WWTP replacement had less effect on instream concentrations of metals and pesticides, BC whose sources are less directly tied to domestic wastewater. Study results indicate that improved WWTP technology can create substantial reductions in concentrations of non-regulated BC such as pharmaceuticals, in addition to regulated contaminants such as nutrients.

Identification of neurotoxicology (NT)/developmental neurotoxicology (DNT) adverse outcome pathways and key event linkages with in vitro DNT screening assays

There is a need to assess compounds reliably and quickly for neurotoxicity (NT) and developmental neurotoxicity (DNT). Adverse outcome pathways (AOPs) enable the mapping of molecular events to an apical endpoint in a chemical agnostic manner and have begun to be applied in NT and DNT testing frameworks. We assessed the status of NT/DNT AOPs in the AOP-Wiki (ca. 2/1/23; https://aopwiki.org/), to characterize the state of AOP development, identify strengths and knowledge gaps, elucidate areas for improvement, and describe areas for future focus. AOPs in the Wiki database were assessed for inclusion of NT/DNT molecular events and endpoints, AOP development and endorsement, as well as the linkages of key neurodevelopmental processes with in vitro new approach methods (NAMs). This review found that 41 AOPs have been proposed detailing NT/DNT, of which eight were endorsed by working parties in OECD. Further, this review determined that learning and memory is included as an adverse outcome in eight NT/DNT AOPS, often without distinction regarding the varying forms of learning and memory, regional specification, temporal dynamics, or acquisition mechanisms involved. There is also an overlap with key events (KEs) and in vitro NAMs, which synaptogenesis appeared as a common process. Overall, progress on NT/DNT AOPs could be expanded, adding in modes of action that are missing, improvement in defining apical endpoints, as well as utilizing NAMs further to develop AOPs and identify gaps in current knowledge.

Assessing Contaminants of Emerging Concern in the Great Lakes Ecosystem: A Decade of Method Development and Practical Application

Assessing the ecological risk of contaminants in the field typically involves consideration of a complex mixture of compounds which may or may not be detected via instrumental analyses. Further, there are insufficient data to predict the potential biological effects of many detected compounds, leading to their being characterized as contaminants of emerging concern (CECs). Over the past several years, advances in chemistry, toxicology, and bioinformatics have resulted in a variety of concepts and tools that can enhance the pragmatic assessment of the ecological risk of CECs. The present Focus article describes a 10+- year multiagency effort supported through the U.S. Great Lakes Restoration Initiative to assess the occurrence and implications of CECs in the North American Great Lakes. State-of-the-science methods and models were used to evaluate more than 700 sites in about approximately 200 tributaries across lakes Ontario, Erie, Huron, Michigan, and Superior, sometimes on multiple occasions. Studies featured measurement of up to 500 different target analytes in different environmental matrices, coupled with evaluation of biological effects in resident species, animals from in situ and laboratory exposures, and in vitro systems. Experimental taxa included birds, fish, and a variety of invertebrates, and measured endpoints ranged from molecular to apical responses. Data were integrated and evaluated using a diversity of curated knowledgebases and models with the goal of producing actionable insights for risk assessors and managers charged with evaluating and mitigating the effects of CECs in the Great Lakes. This overview is based on research and data captured in approximately about 90 peer-reviewed journal articles and reports, including approximately about 30 appearing in a virtual issue comprised of highlighted papers published in Environmental Toxicology and Chemistry or Integrated Environmental Assessment and Management. Environ Toxicol Chem 2023;42:2506-2518. © 2023 SETAC. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

Advancing Computational Toxicology by Interpretable Machine Learning

Chemical toxicity evaluations for drugs, consumer products, and environmental chemicals have a critical impact on human health. Traditional animal models to evaluate chemical toxicity are expensive, time-consuming, and often fail to detect toxicants in humans. Computational toxicology is a promising alternative approach that utilizes machine learning (ML) and deep learning (DL) techniques to predict the toxicity potentials of chemicals. Although the applications of ML- and DL-based computational models in chemical toxicity predictions are attractive, many toxicity models are "black boxes" in nature and difficult to interpret by toxicologists, which hampers the chemical risk assessments using these models. The recent progress of interpretable ML (IML) in the computer science field meets this urgent need to unveil the underlying toxicity mechanisms and elucidate the domain knowledge of toxicity models. In this review, we focused on the applications of IML in computational toxicology, including toxicity feature data, model interpretation methods, use of knowledge base frameworks in IML development, and recent applications. The challenges and future directions of IML modeling in toxicology are also discussed. We hope this review can encourage efforts in developing interpretable models with new IML algorithms that can assist new chemical assessments by illustrating toxicity mechanisms in humans.

Reproducibility of organ-level effects in repeat dose animal studies

This work estimates benchmarks for new approach method (NAM) performance in predicting organ-level effects in repeat dose studies of adult animals based on variability in replicate animal studies. Treatment-related effect values from the Toxicity Reference database (v2.1) for weight, gross, or histopathological changes in the adrenal gland, liver, kidney, spleen, stomach, and thyroid were used. Rates of chemical concordance among organ-level findings in replicate studies, defined by repeated chemical only, chemical and species, or chemical and study type, were calculated. Concordance was 39 - 88%, depending on organ, and was highest within species. Variance in treatment-related effect values, including lowest effect level (LEL) values and benchmark dose (BMD) values when available, was calculated by organ. Multilinear regression modeling, using study descriptors of organ-level effect values as covariates, was used to estimate total variance, mean square error (MSE), and root residual mean square error (RMSE). MSE values, interpreted as estimates of unexplained variance, suggest study descriptors accounted for 52-69% of total variance in organ-level LELs. RMSE ranged from 0.41 - 0.68 log10-mg/kg/day. Differences between organ-level effects from chronic (CHR) and subchronic (SUB) dosing regimens were also quantified. Odds ratios indicated CHR organ effects were unlikely if the SUB study was negative. Mean differences of CHR - SUB organ-level LELs ranged from -0.38 to -0.19 log10 mg/kg/day; the magnitudes of these mean differences were less than RMSE for replicate studies. Finally, in vitro to in vivo extrapolation (IVIVE) was employed to compare bioactive concentrations from in vitro NAMs for kidney and liver to LELs. The observed mean difference between LELs and mean IVIVE dose predictions approached 0.5 log10-mg/kg/day, but differences by chemical ranged widely. Overall, variability in repeat dose organ-level effects suggests expectations for quantitative accuracy of NAM prediction of LELs should be at least ± 1 log10-mg/kg/day, with qualitative accuracy not exceeding 70%.

Anti-androgenic activity of novel flame retardants in mixtures: Newly identified contribution from tris(2,3-dibromopropyl) isocyanurate (TDBP-TAZTO)

In recent decades, male infertility has been on the rise, largely attributed to exposure to chemicals with endocrine-disrupting properties. The adverse effects of disrupting androgen actions on the development and reproductive health of children and adolescents have been extensively studied. Flame retardants (FRs), used in consumer products to delay flammability, have been identified as antagonists of the androgen receptor (AR), potentially leading to adverse outcomes in male reproductive health later in life. This study examined the interaction of eight novel FRs with the AR, employing an in vitro AR-dependent luciferase reporter gene assay utilizing MDA-kb2 cells. The investigation revealed the anti-androgenic activity of tris(2,3-dibromopropyl) isocyanurate (TDBP-TAZTO), a frequently detected FR in the environment. Furthermore, TDBP-TAZTO contributed to anti-androgenic activity when combined with six other anti-androgenic FRs. The mixture effects were predicted by three commonly employed models: concentration addition (CA), generalized CA, and independent action, with the CA model showcasing the highest accuracy. This suggests that all FRs act through a similar mechanism, as further confirmed by in silico molecular docking, indicating limited synergy or antagonism. Importantly, in the mixtures, each FR contributed to the induction of anti-androgenic effects at concentrations below their individual effective concentrations in single exposures. This raises concern for public health, especially considering the co-detection of these FRs and their potential co-occurrence with other anti-androgenic chemicals like bisphenols. Therefore, our findings, along with previous research, strongly support the incorporation of combined effects of mixtures in risk assessment to efficiently safeguard population health.

Evaluation of a high-throughput H295R homogenous time resolved fluorescence assay for androgen and estrogen steroidogenesis screening

The H295R test guideline assay evaluates the effect of test substances on synthesis of 17β-estradiol (E2) and testosterone (T). The objective of this study was to leverage commercial immunoassay technology to develop a more efficient H295R assay to measure E2 and T levels in 384-well format. The resulting Homogenous Time Resolved Fluorescence assay platform (H295R-HTRF) was evaluated against a training set of 36 chemicals derived from the OECD inter-laboratory validation study, EPA guideline 890.1200 aromatase assay, and azole fungicides active in the HT-H295R assay. Quality control performance criteria were met for all conditions except E2 synthesis inhibition where low basal hormone synthesis was observed. Five proficiency chemicals were active for both the E2 and T endpoints, consistent with guideline classifications. Of the nine OECD core reference chemicals, 9/9 were concordant with outcomes for E2 and 7/9 for T. Likewise, 9/13 and 11/13 OECD supplemental chemicals were concordant with anticipated effects for E2 and T, respectively. Of the 10 azole fungicides screened, 7/10 for E2 and 8/10 for T exhibited concordant outcomes for inhibition. Generally, all active chemicals in the training set demonstrated equivalent or greater potency in the H295R-HTRF assay, supporting the sensitivity of the platform. The adaptation of HTRF technology to the H295R model provides an efficient way to evaluate E2 and T modulators in accordance with guideline specifications.

Environmental metabolomics uncovers oxidative stress, amino acid dysregulation, and energy impairment in Daphnia magna with exposure to industrial effluents

Anthropogenic activities are regarded as point sources of pollution entering freshwater bodies worldwide. With over 350,000 chemicals used in manufacturing, wastewater treatment and industrial effluents are comprised of complex mixtures of organic and inorganic pollutants of known and unknown origins. Consequently, their combined toxicity and mode of action are not well understood in aquatic organisms such as Daphnia magna. In this study, effluent samples from wastewater treatment and industrial sectors were used to examine molecular-level perturbations to the polar metabolic profile of D. magna. To determine if the industrial sector and/or the effluent chemistries played a role in the observed biochemical responses, Daphnia were acutely (48 h) exposed to undiluted (100%) and diluted (10, 25, and 50%) effluent samples. Endogenous metabolites were extracted from single daphnids and analyzed using targeted mass spectrometry-based metabolomics. The metabolic profile of Daphnia exposed to effluent samples resulted in significant separation compared to the unexposed controls. Linear regression analysis determined that no single pollutant detected in the effluents was significantly correlated with the responses of metabolites. Significant perturbations were uncovered across many classes of metabolites (amino acids, nucleosides, nucleotides, polyamines, and their derivatives) which serve as intermediates in keystone biochemical processes. The combined metabolic responses are consistent with oxidative stress, disruptions to energy metabolism, and protein dysregulation which were identified through biochemical pathway analysis. These results provide insight into the molecular processes driving stress responses in D. magna. Overall, we determined that the metabolic profile of Daphnia could not be predicted by the chemical composition of environmentally relevant mixtures. The findings of this study demonstrate the advantage of metabolomics in conjunction with chemical analyses to assess the interactions of industrial effluents. This work further demonstrates the ability of environmental metabolomics to characterize molecular-level perturbations in aquatic organisms exposed to complex chemical mixtures directly.

Comparison of in silico, in vitro, and in vivo toxicity benchmarks suggests a role for ToxCast data in ecological hazard assessment

Large repositories of in vitro bioactivity data such as US EPA's Toxicity Forecaster (ToxCast) provide a wealth of publicly accessible toxicity information for thousands of chemicals. These data can be used to calculate point-of-departure (POD) estimates via concentration-response modeling that may serve as lower bound, protective estimates of in vivo effects. However, the data are predominantly based on mammalian models and discussions to date about their utility have largely focused on potential integration into human hazard assessment, rather than application to ecological risk assessment. The goal of the present study was to compare PODs based on (1) quantitative structure-activity relationships (QSARs), (2) the 5th centile of the activity concentration at cutoff (ACC), and (3) lower-bound cytotoxic burst (LCB) from ToxCast, with the distribution of in vivo PODs compiled in the Ecotoxicology Knowledgebase (ECOTOX). While overall correlation between ToxCast ACC5 and ECOTOX PODs for 649 chemicals was weak, there were significant associations among PODs based on LCB and ECOTOX, LCB and QSARs, and ECOTOX and QSARs. Certain classes of compounds showed moderate correlation across datasets (eg, antimicrobials/disinfectants), while others, such as organophosphate insecticides, did not. Unsurprisingly, more precise classifications of the data based on ECOTOX effect and endpoint type (eg, apical vs biochemical; acute vs chronic) had a significant effect on overall relationships. Results of this research help to define appropriate roles for data from new approach methodologies in chemical prioritization and screening of ecological hazards.

The ToxCast pipeline: updates to curve-fitting approaches and database structure

Introduction: The US Environmental Protection Agency Toxicity Forecaster (ToxCast) program makes in vitro medium- and high-throughput screening assay data publicly available for prioritization and hazard characterization of thousands of chemicals. The assays employ a variety of technologies to evaluate the effects of chemical exposure on diverse biological targets, from distinct proteins to more complex cellular processes like mitochondrial toxicity, nuclear receptor signaling, immune responses, and developmental toxicity. The ToxCast data pipeline (tcpl) is an open-source R package that stores, manages, curve-fits, and visualizes ToxCast data and populates the linked MySQL Database, invitrodb. Methods: Herein we describe major updates to tcpl and invitrodb to accommodate a new curve-fitting approach. The original tcpl curve-fitting models (constant, Hill, and gain-loss models) have been expanded to include Polynomial 1 (Linear), Polynomial 2 (Quadratic), Power, Exponential 2, Exponential 3, Exponential 4, and Exponential 5 based on BMDExpress and encoded by the R package dependency, tcplfit2. Inclusion of these models impacted invitrodb (beta version v4.0) and tcpl v3 in several ways: (1) long-format storage of generic modeling parameters to permit additional curve-fitting models; (2) updated logic for winning model selection; (3) continuous hit calling logic; and (4) removal of redundant endpoints as a result of bidirectional fitting. Results and discussion: Overall, the hit call and potency estimates were largely consistent between invitrodb v3.5 and 4.0. Tcpl and invitrodb provide a standard for consistent and reproducible curve-fitting and data management for diverse, targeted in vitro assay data with readily available documentation, thus enabling sharing and use of these data in myriad toxicology applications. The software and database updates described herein promote comparability across multiple tiers of data within the US Environmental Protection Agency CompTox Blueprint.

Computational Construction of Toxicant Signaling Networks

Humans and animals are regularly exposed to compounds that may have adverse effects on health. The Toxicity Forecaster (ToxCast) program was developed to use high throughput screening assays to quickly screen chemicals by measuring their effects on many biological end points. Many of these assays test for effects on cellular receptors and transcription factors (TFs), under the assumption that a toxicant may perturb normal signaling pathways in the cell. We hypothesized that we could reconstruct the intermediate proteins in these pathways that may be directly or indirectly affected by the toxicant, potentially revealing important physiological processes not yet tested for many chemicals. We integrate data from ToxCast with a human protein interactome to build toxicant signaling networks that contain physical and signaling protein interactions that may be affected as a result of toxicant exposure. To build these networks, we developed the EdgeLinker algorithm, which efficiently finds short paths in the interactome that connect the receptors to TFs for each toxicant. We performed multiple evaluations and found evidence suggesting that these signaling networks capture biologically relevant effects of toxicants. To aid in dissemination and interpretation, interactive visualizations of these networks are available at http://graphspace.org.

An unbiased, automated platform for scoring dopaminergic neurodegeneration in C. elegans

Caenorhabditis elegans (C. elegans) has served as a simple model organism to study dopaminergic neurodegeneration, as it enables quantitative analysis of cellular and sub-cellular morphologies in live animals. These isogenic nematodes have a rapid life cycle and transparent body, making high-throughput imaging and evaluation of fluorescently tagged neurons possible. However, the current state-of-the-art method for quantifying dopaminergic degeneration requires researchers to manually examine images and score dendrites into groups of varying levels of neurodegeneration severity, which is time consuming, subject to bias, and limited in data sensitivity. We aim to overcome the pitfalls of manual neuron scoring by developing an automated, unbiased image processing algorithm to quantify dopaminergic neurodegeneration in C. elegans. The algorithm can be used on images acquired with different microscopy setups and only requires two inputs: a maximum projection image of the four cephalic neurons in the C. elegans head and the pixel size of the user's camera. We validate the platform by detecting and quantifying neurodegeneration in nematodes exposed to rotenone, cold shock, and 6-hydroxydopamine using 63x epifluorescence, 63x confocal, and 40x epifluorescence microscopy, respectively. Analysis of tubby mutant worms with altered fat storage showed that, contrary to our hypothesis, increased adiposity did not sensitize to stressor-induced neurodegeneration. We further verify the accuracy of the algorithm by comparing code-generated, categorical degeneration results with manually scored dendrites of the same experiments. The platform, which detects 20 different metrics of neurodegeneration, can provide comparative insight into how each exposure affects dopaminergic neurodegeneration patterns.

Characterizing surface water concentrations of hundreds of organic chemicals in United States for environmental risk prioritization

BACKGROUND

Thousands of chemicals are observed in freshwater, typically at trace levels. Measurements are collected for different purposes, so sample characteristics vary. Due to inconsistent data availability for exposure and hazard, it is complex to prioritize which chemicals may pose risks.

OBJECTIVE

We evaluated the influence of data curation and statistical practices aggregating surface water measurements of organic chemicals into exposure distributions intended for prioritizing based on nation-scale potential risk.

METHODS

The Water Quality Portal includes millions of observations describing over 1700 chemicals in 93% of hydrologic subbasins across the United States. After filtering to maintain quality and applicability while including all possible samples, we compared concentrations across sample types. We evaluated statistical methods to estimate per-chemical distributions for chosen samples. Overlaps between resulting exposure ranges and distributions representing no-effect concentrations for multiple freshwater species were used to rank estimated chemical risks for further assessment.

RESULTS

When we apply explicit data quality and statistical assumptions, we find that there are 186 organic chemicals for which we can make screening-level estimates of surface water chemical concentration. Of the original 1700 observed chemicals, this number decreased primarily due to a predominance of censored values (that is, observations indicating concentrations too low to be measured). We further identify 423 chemicals where all measurements were censored but, through consideration of detection limits, risk might still be prioritized based on the detection limits themselves. In the final set of 1.5 million samples, the median environmental concentration of one chemical (acetic acid) exceeded the 5th percentile of no-effect concentrations for the most delicate freshwater species (the highest priority risk condition identified here), and a further 29 chemicals were identified for possible further evaluation based on a small margin between occurrence and toxicity values.

SIGNIFICANCE

This method shows the broad range of chemical concentrations seen for organic chemicals across the country and identifies methods of determining their central tendency, allowing for researchers to characterize higher-than-normal or lower-than-normal surface water conditions as well as providing an overall indication of the presence of organic chemicals in the United States. The highest chemical concentrations did not always indicate the highest-risk conditions. Even when accounting for the high level of uncertainty in these data due to differences in data collection and reporting across the set, some chemicals may still be categorized as higher environmental risk than others using this method, providing value to chemical safety decision makers and researchers by suggesting avenues for more focused investigation.

Evaluating pharmaceuticals and other organic contaminants in the Lac du Flambeau Chain of Lakes using risk-based screening techniques

In an investigation of pharmaceutical contamination in the Lac du Flambeau Chain of Lakes (hereafter referred to as "the Chain"), few contaminants were detected; only eight pharmaceuticals and one pesticide were identified among the 110 pharmaceuticals and other organic contaminants monitored in surface water samples. This study, conducted in cooperation with the Lac du Flambeau Tribe's Water Resource Program, investigated these organic contaminants and potential biological effects in channels connecting lakes throughout the Chain, including the Moss Lake Outlet site, adjacent to the wastewater treatment plant lagoon. Of the 6 sites monitored and 24 samples analyzed, sample concentrations and contaminant detection frequencies were greatest at the Moss Lake Outlet site; however, the concentrations and detection frequencies of this study were comparable to other pharmaceutical investigations in basins with similar characteristics. Because established water-quality benchmarks do not exist for the pharmaceuticals detected in this study, alternative screening-level water-quality benchmarks, developed using two U.S. Environmental Protection Agency toxicological resources (ToxCast database and ECOTOX knowledgebase), were used to estimate potential biological effects associated with the observed contaminant concentrations. Two contaminants (caffeine and thiabendazole) exceeded the prioritization threshold according to ToxCast alternative benchmarks, and four contaminants (acetaminophen, atrazine, caffeine, and carbamazepine) exceeded the prioritization threshold according to ECOTOX alternative benchmarks. Atrazine, an herbicide, was the most frequently detected contaminant (79% of samples), and it exhibited the strongest potential for biological effects due to its high estimated potency. Insufficient toxicological information within ToxCast and ECOTOX for gabapentin and methocarbamol (which had the two greatest concentrations in this study) precluded alternative benchmark development. This data gap presents unknown potential environmental impacts. Future research examining the biological effects elicited by these two contaminants as well as the others detected in this study would further elucidate the ecological relevance of the water chemistry results generated though this investigation.

Evaluation of Complex Mixture Toxicity in the Milwaukee Estuary (WI, USA) Using Whole-Mixture and Component-Based Evaluation Methods

Anthropogenic activities introduce complex mixtures into aquatic environments, necessitating mixture toxicity evaluation during risk assessment. There are many alternative approaches that can be used to complement traditional techniques for mixture assessment. Our study aimed to demonstrate how these approaches could be employed for mixture evaluation in a target watershed. Evaluations were carried out over 2 years (2017-2018) across 8-11 study sites in the Milwaukee Estuary (WI, USA). Whole mixtures were evaluated on a site-specific basis by deploying caged fathead minnows (Pimephales promelas) alongside composite samplers for 96 h and characterizing chemical composition, in vitro bioactivity of collected water samples, and in vivo effects in whole organisms. Chemicals were grouped based on structure/mode of action, bioactivity, and pharmacological activity. Priority chemicals and mixtures were identified based on their relative contributions to estimated mixture pressure (based on cumulative toxic units) and via predictive assessments (random forest regression). Whole mixture assessments identified target sites for further evaluation including two sites targeted for industrial/urban chemical mixture effects assessment; three target sites for pharmaceutical mixture effects assessment; three target sites for further mixture characterization; and three low-priority sites. Analyses identified 14 mixtures and 16 chemicals that significantly contributed to cumulative effects, representing high or medium priority targets for further ecotoxicological evaluation, monitoring, or regulatory assessment. Overall, our study represents an important complement to single-chemical prioritizations, providing a comprehensive evaluation of the cumulative effects of mixtures detected in a target watershed. Furthermore, it demonstrates how different tools and techniques can be used to identify diverse facets of mixture risk and highlights strategies that can be considered in future complex mixture assessments. Environ Toxicol Chem 2023;42:1229-1256. © 2023 SETAC.

Data-Driven Quantitative Structure-Activity Relationship Modeling for Human Carcinogenicity by Chronic Oral Exposure

Traditional methodologies for assessing chemical toxicity are expensive and time-consuming. Computational modeling approaches have emerged as low-cost alternatives, especially those used to develop quantitative structure-activity relationship (QSAR) models. However, conventional QSAR models have limited training data, leading to low predictivity for new compounds. We developed a data-driven modeling approach for constructing carcinogenicity-related models and used these models to identify potential new human carcinogens. To this goal, we used a probe carcinogen dataset from the US Environmental Protection Agency's Integrated Risk Information System (IRIS) to identify relevant PubChem bioassays. Responses of 25 PubChem assays were significantly relevant to carcinogenicity. Eight assays inferred carcinogenicity predictivity and were selected for QSAR model training. Using 5 machine learning algorithms and 3 types of chemical fingerprints, 15 QSAR models were developed for each PubChem assay dataset. These models showed acceptable predictivity during 5-fold cross-validation (average CCR = 0.71). Using our QSAR models, we can correctly predict and rank 342 IRIS compounds' carcinogenic potentials (PPV = 0.72). The models predicted potential new carcinogens, which were validated by a literature search. This study portends an automated technique that can be applied to prioritize potential toxicants using validated QSAR models based on extensive training sets from public data resources.

Prediction of zebrafish embryonic developmental toxicity by integrating omics with adverse outcome pathway

New approach methodologies (NAMs), especially omics-based high-throughput bioassays have been developed rapidly, providing rich mechanistic information such as molecular initiation events (MIEs) and (sub)cellular key events (KEs) in adverse outcome pathways (AOPs). However, how to apply the knowledge of MIEs/KEs to predict adverse outcomes (AOs) induced by chemicals represents a new challenge for computational toxicology. Here, an integrated method named ScoreAOP was developed and evaluated to predict chemicals' developmental toxicity for zebrafish embryos by integrating four related AOPs and dose-dependent reduced zebrafish transcriptome (RZT). The rules of ScoreAOP included 1) sensitivity of responsive KEs demonstrated by point of departure of KEs (POD_KE), 2) evidence reliability and 3) distance between KEs and AOs. Moreover, eleven chemicals with different modes of action (MoAs) were tested to evaluate ScoreAOP. Results showed that eight of the eleven chemicals caused developmental toxicity at tested concentration in apical tests. All the tested chemicals' developmental defects were predicted using ScoreAOP, whereas eight out of the eleven chemicals predicted by ScoreMIE which was developed to score MIEs disturbed by chemicals based on in vitro bioassays data. Finally, in terms of mechanism explanation, ScoreAOP clustered chemicals with different MoAs while ScoreMIE failed, and ScoreAOP revealed the activation of aryl hydrocarbon receptor (AhR) plays a significant role in dysfunction of cardiovascular system, resulting in zebrafish developmental defects and mortality. In conclusion, ScoreAOP represents a promising approach to apply mechanism information obtained from omics to predict AOs induced by chemicals.

A Dynamic In vitro developing testis model reflects structures and functions of testicular development in vivo

To better define appropriate applications of our 3-dimensional testicular co-culture as a model for reproductive toxicology, we evaluated the ability of the model to capture structural and functional elements that can be targeted by reproductive toxicants. Testicular co-cultures were prepared from postnatal day 5 male rats and cultured with a Matrigel overlay. Following a 2-day acclimation period, we characterized functional pathway dynamics by evaluating morphology, protein expression, testosterone concentrations, and global gene expression at a range of timepoints from experimental days 0 to 21. Western blotting confirmed expression of Sertoli cell, Leydig cell, and spermatogonial cell-specific protein markers. Testosterone detected in cell culture media indicates active testosterone production. Quantitative pathway analysis identified Gene Ontology biological processes enriched among genes significantly changing over the course of 21 days. Processes enriched among genes significantly increasing through time include general developmental processes (morphogenesis, tissue remodeling, etc.), steroid regulation, Sertoli cell development, immune response, and stress and apoptosis. Processes enriched among genes significantly decreasing over time include several related to male reproductive development (seminiferous tubule development, male gonad development, Leydig cell differentiation, Sertoli cell differentiation), all of which appear to peak in expression between days 1 and 5 before decreasing at later timepoints. This analysis provides a temporal roadmap for specific biological process of interest for reproductive toxicology in the model and anchors the model to sensitive phases of in vivo development, helping to define the relevance of the model for in vivo processes.

Quantitative analyses of misclassification rates in the hazard classification of environmental chemicals: Evaluation of procedures for deriving predicted no-effect concentrations in the Chemical Substance Control Law in Japan

The quality of chemical management depends more or less on practical procedures used to assess chemicals. This study quantitatively assessed the efficacy of a derivation procedure for calculating no-effect concentrations for screening assessment of environmental hazards under the Chemical Substance Control Law in Japan. We first evaluated the derivation procedure by applying a series of test ecotoxicity datasets to the procedure and calculating the resulting misclassification rates of the hazardous class of chemicals. In this study, a chemical was deemed to have been misclassified if its classification differed from its classification based on the full dataset (chronic toxicity data for three trophic levels), which was defined as the correct assignment. We also calculated the effects of additional uncertainty factors to decrease the variance (i.e., to improve the consistency) of the misclassification rates among cases with different data availability in the derivation procedure. The results showed that the derivation procedure resulted in very high rates of misclassification when only particular sets of ecotoxicity data were available (e.g., only chronic toxicity data of algae were available). Our analyses also showed that the use of additional uncertainty factors improved the consistency of the misclassification rates within the derivation procedure. Our study presents a broadly applicable calculation framework for quantifying error rates in assessment procedures and serves as a case study for future development and reforms of chemical assessment processes and policies, while additional analyses using more extensive ecotoxicity data with various modes of actions are needed in the future.

Probabilistic Points of Departure and Reference Doses for Characterizing Human Noncancer and Developmental/Reproductive Effects for 10,145 Chemicals

BACKGROUND

Regulatory toxicity values used to assess and manage chemical risks rely on the determination of the point of departure (POD) for a critical effect, which results from a comprehensive and systematic assessment of available toxicity studies. However, regulatory assessments are only available for a small fraction of chemicals.

OBJECTIVES

Using in vivo experimental animal data from the U.S. Environmental Protection Agency's Toxicity Value Database, we developed a semiautomated approach to determine surrogate oral route PODs, and corresponding toxicity values where regulatory assessments are unavailable.

METHODS

We developed a curated data set restricted to effect levels, exposure routes, study designs, and species relevant for deriving toxicity values. Effect levels were adjusted to chronic human equivalent benchmark doses (BMDh). We hypothesized that a quantile of the BMDh distribution could serve as a surrogate POD and determined the appropriate quantile by calibration to regulatory PODs. Finally, we characterized uncertainties around the surrogate PODs from intra- and interstudy variability and derived probabilistic toxicity values using a standardized workflow.

RESULTS

The BMDh distribution for each chemical was adequately fit by a lognormal distribution, and the 25th percentile best predicted the available regulatory PODs [R2≥0.78, residual standard error (RSE)≤0.53 log10 units]. We derived surrogate PODs for 10,145 chemicals from the curated data set, differentiating between general noncancer and reproductive/developmental effects, with typical uncertainties (at 95% confidence) of a factor of 10 and 12, respectively. From these PODs, probabilistic reference doses (1% incidence at 95% confidence), as well as human population effect doses (10% incidence), were derived.

DISCUSSION

In providing surrogate PODs calibrated to regulatory values and deriving corresponding toxicity values, we have substantially expanded the coverage of chemicals from 744 to 8,023 for general noncancer effects, and from 41 to 6,697 for reproductive/developmental effects. These results can be used across various risk assessment and risk management contexts, from hazardous site and life cycle impact assessments to chemical prioritization and substitution. https://doi.org/10.1289/EHP11524.

A value of information framework for assessing the trade-offs associated with uncertainty, duration, and cost of chemical toxicity testing

A number of investigators have explored the use of value of information (VOI) analysis to evaluate alternative information collection procedures in diverse decision-making contexts. This paper presents an analytic framework for determining the value of toxicity information used in risk-based decision making. The framework is specifically designed to explore the trade-offs between cost, timeliness, and uncertainty reduction associated with different toxicity-testing methodologies. The use of the proposed framework is demonstrated by two illustrative applications which, although based on simplified assumptions, show the insights that can be obtained through the use of VOI analysis. Specifically, these results suggest that timeliness of information collection has a significant impact on estimates of the VOI of chemical toxicity tests, even in the presence of smaller reductions in uncertainty. The framework introduces the concept of the expected value of delayed sample information, as an extension to the usual expected value of sample information, to accommodate the reductions in value resulting from delayed decision making. Our analysis also suggests that lower cost and higher throughput testing also may be beneficial in terms of public health benefits by increasing the number of substances that can be evaluated within a given budget. When the relative value is expressed in terms of return-on-investment per testing strategy, the differences can be substantial.

Petroleum Hydrocarbon Catabolic Pathways as Targets for Metabolic Engineering Strategies for Enhanced Bioremediation of Crude-Oil-Contaminated Environments

Anthropogenic activities and industrial effluents are the major sources of petroleum hydrocarbon contamination in different environments. Microbe-based remediation techniques are known to be effective, inexpensive, and environmentally safe. In this review, the metabolic-target-specific pathway engineering processes used for improving the bioremediation of hydrocarbon-contaminated environments have been described. The microbiomes are characterised using environmental genomics approaches that can provide a means to determine the unique structural, functional, and metabolic pathways used by the microbial community for the degradation of contaminants. The bacterial metabolism of aromatic hydrocarbons has been explained via peripheral pathways by the catabolic actions of enzymes, such as dehydrogenases, hydrolases, oxygenases, and isomerases. We proposed that by using microbiome engineering techniques, specific pathways in an environment can be detected and manipulated as targets. Using the combination of metabolic engineering with synthetic biology, systemic biology, and evolutionary engineering approaches, highly efficient microbial strains may be utilised to facilitate the target-dependent bioprocessing and degradation of petroleum hydrocarbons. Moreover, the use of CRISPR-cas and genetic engineering methods for editing metabolic genes and modifying degradation pathways leads to the selection of recombinants that have improved degradation abilities. The idea of growing metabolically engineered microbial communities, which play a crucial role in breaking down a range of pollutants, has also been explained. However, the limitations of the in-situ implementation of genetically modified organisms pose a challenge that needs to be addressed in future research.

An unbiased, automated platform for scoring dopaminergic neurodegeneration in C. elegans

Caenorhabditis elegans ( C. elegans ) has served as a simple model organism to study dopaminergic neurodegeneration, as it enables quantitative analysis of cellular and sub-cellular morphologies in live animals. These isogenic nematodes have a rapid life cycle and transparent body, making high-throughput imaging and evaluation of fluorescently tagged neurons possible. However, the current state-of-the-art method for quantifying dopaminergic degeneration requires researchers to manually examine images and score dendrites into groups of varying levels of neurodegeneration severity, which is time consuming, subject to bias, and limited in data sensitivity. We aim to overcome the pitfalls of manual neuron scoring by developing an automated, unbiased image processing algorithm to quantify dopaminergic neurodegeneration in C. elegans . The algorithm can be used on images acquired with different microscopy setups and only requires two inputs: a maximum projection image of the four cephalic neurons in the C. elegans head and the pixel size of the user’s camera. We validate the platform by detecting and quantifying neurodegeneration in nematodes exposed to rotenone, cold shock, and 6-hydroxydopamine using 63x epifluorescence, 63x confocal, and 40x epifluorescence microscopy, respectively. Analysis of tubby mutant worms with altered fat storage showed that, contrary to our hypothesis, increased adiposity did not sensitize to stressor-induced neurodegeneration. We further verify the accuracy of the algorithm by comparing code-generated, categorical degeneration results with manually scored dendrites of the same experiments. The platform, which detects 19 different metrics of neurodegeneration, can provide comparative insight into how each exposure affects dopaminergic neurodegeneration patterns.

Pesticide Prioritization by Potential Biological Effects in Tributaries of the Laurentian Great Lakes

Watersheds of the Great Lakes Basin (USA/Canada) are highly modified and impacted by human activities including pesticide use. Despite labeling restrictions intended to minimize risks to nontarget organisms, concerns remain that environmental exposures to pesticides may be occurring at levels negatively impacting nontarget organisms. We used a combination of organismal-level toxicity estimates (in vivo aquatic life benchmarks) and data from high-throughput screening (HTS) assays (in vitro benchmarks) to prioritize pesticides and sites of concern in streams at 16 tributaries to the Great Lakes Basin. In vivo or in vitro benchmark values were exceeded at 15 sites, 10 of which had exceedances throughout the year. Pesticides had the greatest potential biological impact at the site with the greatest proportion of agricultural land use in its basin (the Maumee River, Toledo, OH, USA), with 72 parent compounds or transformation products being detected, 47 of which exceeded at least one benchmark value. Our risk-based screening approach identified multiple pesticide parent compounds of concern in tributaries of the Great Lakes; these compounds included: eight herbicides (metolachlor, acetochlor, 2,4-dichlorophenoxyacetic acid, diuron, atrazine, alachlor, triclopyr, and simazine), three fungicides (chlorothalonil, propiconazole, and carbendazim), and four insecticides (diazinon, fipronil, imidacloprid, and clothianidin). We present methods for reducing the volume and complexity of potential biological effects data that result from combining contaminant surveillance with HTS (in vitro) and traditional (in vivo) toxicity estimates. Environ Toxicol Chem 2023;42:367-384. Published 2022. This article is a U.S. Government work and is in the public domain in the USA. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

High-content analysis shows synergistic effects of low perfluorooctanoic acid (PFOS) and perfluorooctane sulfonic acid (PFOA) mixture concentrations on human breast epithelial cell carcinogenesis

Perfluoroalkyl substances (PFAS) have been associated with cancer, but the potential underlying mechanisms need to be further elucidated and include studies of PFAS mixtures. This mechanistic study revealed that very low concentrations (500 pM) of the binary PFOS and PFOA mixture induced synergistic effects on human epithelial breast cell (MCF-10A) proliferation. The cell proliferation was mediated by pregnane X receptor (PXR) activation, an increase in cyclin D1 and CDK6/4 levels, decrease in p21 and p53 levels, and by regulation of phosphor-Akt and β-catenin. The PFAS mixture also altered histone modifications, epigenetic mechanisms implicated in tumorigenesis, and promoted cell migration and invasion by reducing the levels of occludin. High-content screening using the cell painting assay, revealed that hundreds of cell features were affected by the PFAS mixture even at the lowest concentration tested (100 pM). The detailed phenotype profiling further demonstrated that the PFAS mixture altered cell morphology, mostly in parameters related to intensity and texture associated with mitochondria, endoplasmic reticulum, and nucleoli. Exposure to higher concentrations (≥50 µM) of the PFOS and PFOA mixture caused cell death through synergistic interactions that induced oxidative stress, DNA/RNA damage, and lipid peroxidation, illustrating the complexity of mixture toxicology. Increased knowledge about mixture-induced effects is important for better understanding of PFAS' possible role in cancer etiology, and may impact the risk assessment of these and other compounds. This study shows the potential of image-based multiplexed fluorescence assays and high-content screening for development of new approach methodologies in toxicology.

Formation and Inventory of Polychlorinated Dibenzo-p-dioxins and Dibenzofurans and Other Byproducts along Manufacturing Processes of Chlorobenzene and Chloroethylene

Chlorinated organic chemicals are produced and used extensively worldwide, and their risks to the biology and environment are of increasing concern. However, chlorinated byproducts [e.g., polychlorinated dibenzo-p-dioxins and dibenzofurans (PCDD/Fs)] formed during the commercial manufacturing processes and present in organochlorine products are rarely reported. The knowledge on the occurrences and fate of unintentional persistent organic chemicals in the manufacturing of organochlorine chemical is necessary for accurate assessment of the risks of commercial chemicals and their production. Here, PCDD/Fs were tracked throughout chlorobenzene and chloroethylene production processes (from raw materials to final products) by target analysis. Other byproducts that can further transform into PCDD/Fs were also identified by performing non-target screening. As a result, the PCDD/F concentrations were mostly the highest in bottom residues, and the octachlorinated congeners were dominant. Alkali/water washing stages may cause the formation of oxygen-containing byproducts including PCDD/Fs and acyl-containing compounds, so more attention should be paid to these stages. PCDD/Fs were of 0.17 and 0.21-1.2 ng/mL in monochlorobenzene and chloroethylene products, respectively. Annual PCDD/F emissions (17 g toxic equivalent in 2018) during chlorobenzene and chloroethylene production were estimated using PCDD/F emission factors. The results can contribute to the improvement of PCDD/F inventories for the analyzed commercial chemicals.

Targeted Per- and Polyfluoroalkyl substances (PFAS) assessments for high throughput screening: Analytical and testing considerations to inform a PFAS stock quality evaluation framework

Per- and polyfluoroalkyl substances (PFAS) represent a large chemical class lacking hazard, toxicokinetic, and exposure information. To accelerate PFAS hazard evaluation, new approach methodologies (NAMs) comprised of in vitro high-throughput toxicity screening, toxicokinetic data, and computational modeling are being employed in read across strategies to evaluate the larger PFAS landscape. A critical consideration to ensure robust evaluations is a parallel assessment of the quality of the screening stock solutions, where dimethyl sulfoxide (DMSO) is often the diluent of choice. Challenged by the lack of commercially available reference standards for many of the selected PFAS and reliance on mass spectrometry approaches for such an evaluation, we developed a high-throughput framework to evaluate the quality of screening stocks for 205 PFAS selected for these NAM efforts. Using mass spectrometry coupled with either liquid or gas chromatography, a quality scoring system was developed that incorporated observations during mass spectral examination to provide a simple pass or fail notation. Informational flags were used to further describe findings regarding parent analyte presence through accurate mass identification, evidence of contaminants and/or degradation, or further describe characteristics such as isomer presence. Across the PFAS-DMSO stocks tested, 148 unique PFAS received passing quality scores to allow for further in vitro testing whereas 57 received a failing score primarily due to detection issues or confounding effects of DMSO. Principle component analysis indicated vapor pressure and Henry's Law Constant as top indicators for a failed quality score for those analyzed by gas chromatography. Three PFAS in the hexafluoropropylene oxide family failed due to degradation in DMSO. As the PFAS evaluated spanned over 20 different structural categories, additional commentary describes analytical observations across specific groups related to PFAS stock composition, detection, stability, and methodologic considerations that will be useful for informing future analytical assessment and downstream HTS efforts. The high-throughput stock quality scoring workflow presented holds value as a tool to evaluate chemical presence and quality efficiently and for informing data inclusion in PFAS or other NAM screening efforts.

A geospatial modeling approach to quantifying the risk of exposure to environmental chemical mixtures via a common molecular target

In the real world, individuals are exposed to chemicals from sources that vary over space and time. However, traditional risk assessments based on in vivo animal studies typically use a chemical-by-chemical approach and apical disease endpoints. New approach methodologies (NAMs) in toxicology, such as in vitro high-throughput (HTS) assays generated in Tox21 and ToxCast, can more readily provide mechanistic chemical hazard information for chemicals with no existing data than in vivo methods. In this paper, we establish a workflow to assess the joint action of 41 modeled ambient chemical exposures in the air from the USA-wide National Air Toxics Assessment by integrating human exposures with hazard data from curated HTS (cHTS) assays to identify counties where exposure to the local chemical mixture may perturb a common biological target. We exemplify this proof-of-concept using CYP1A1 mRNA up-regulation. We first estimate internal exposure and then convert the inhaled concentration to a steady state plasma concentration using physiologically based toxicokinetic modeling parameterized with county-specific information on ages and body weights. We then use the estimated blood plasma concentration and the concentration-response curve from the in vitro cHTS assay to determine the chemical-specific effects of the mixture components. Three mixture modeling methods were used to estimate the joint effect from exposure to the chemical mixture on the activity levels, which were geospatially mapped. Finally, a Monte Carlo uncertainty analysis was performed to quantify the influence of each parameter on the combined effects. This workflow demonstrates how NAMs can be used to predict early-stage biological perturbations that can lead to adverse health outcomes that result from exposure to chemical mixtures. As a result, this work will advance mixture risk assessment and other early events in the effects of chemicals.

Construction of an In Silico Structural Profiling Tool Facilitating Mechanistically Grounded Classification of Aquatic Toxicants

The performance of chemical safety assessment within the domain of environmental toxicology is often impeded by a shortfall of appropriate experimental data describing potential hazards across the many compounds in regular industrial use. In silico schemes for assigning aquatic-relevant modes or mechanisms of toxic action to substances, based solely on consideration of chemical structure, have seen widespread employment─including those of Verhaar, Russom, and later Bauer (MechoA). Recently, development of a further system was reported by Sapounidou, which, in common with MechoA, seeks to ground its classifications in understanding and appreciation of molecular initiating events. Until now, this Sapounidou scheme has not seen implementation as a tool for practical screening use. Accordingly, the primary purpose of this study was to create such a resource─in the form of a computational workflow. This exercise was facilitated through the formulation of 183 structural alerts/rules describing molecular features associated with narcosis, chemical reactivity, and specific mechanisms of action. Output was subsequently compared relative to that of the three aforementioned alternative systems to identify strengths and shortcomings as regards coverage of chemical space.

Simulating toxicokinetic variability to identify susceptible and highly exposed populations

BACKGROUND

Toxicokinetic (TK) data needed for chemical risk assessment are not available for most chemicals. To support a greater number of chemicals, the U.S. Environmental Protection Agency (EPA) created the open-source R package "httk" (High Throughput ToxicoKinetics). The "httk" package provides functions and data tables for simulation and statistical analysis of chemical TK, including a population variability simulator that uses biometrics data from the National Health and Nutrition Examination Survey (NHANES).

OBJECTIVE

Here we modernize the "HTTK-Pop" population variability simulator based on the currently available data and literature. We provide explanations of the algorithms used by "httk" for variability simulation and uncertainty propagation.

METHODS

We updated and revised the population variability simulator in the "httk" package with the most recent NHANES biometrics (up to the 2017-18 NHANES cohort). Model equations describing glomerular filtration rate (GFR) were revised to more accurately represent physiology and population variability. The model output from the updated "httk" package was compared with the current version.

RESULTS

The revised population variability simulator in the "httk" package now provides refined, more relevant, and better justified estimations.

SIGNIFICANCE

Fulfilling the U.S. EPA's mission to provide open-source data and models for evaluations and applications by the broader scientific community, and continuously improving the accuracy of the "httk" package based on the currently available data and literature.

Exposure forecasting - ExpoCast - for data-poor chemicals in commerce and the environment

Estimates of exposure are critical to prioritize and assess chemicals based on risk posed to public health and the environment. The U.S. Environmental Protection Agency (EPA) is responsible for regulating thousands of chemicals in commerce and the environment for which exposure data are limited. Since 2009 the EPA's ExpoCast ("Exposure Forecasting") project has sought to develop the data, tools, and evaluation approaches required to generate rapid and scientifically defensible exposure predictions for the full universe of existing and proposed commercial chemicals. This review article aims to summarize issues in exposure science that have been addressed through initiatives affiliated with ExpoCast. ExpoCast research has generally focused on chemical exposure as a statistical systems problem intended to inform thousands of chemicals. The project exists as a companion to EPA's ToxCast ("Toxicity Forecasting") project which has used in vitro high-throughput screening technologies to characterize potential hazard posed by thousands of chemicals for which there are limited toxicity data. Rapid prediction of chemical exposures and in vitro-in vivo extrapolation (IVIVE) of ToxCast data allow for prioritization based upon risk of adverse outcomes due to environmental chemical exposure. ExpoCast has developed (1) integrated modeling approaches to reliably predict exposure and IVIVE dose, (2) highly efficient screening tools for chemical prioritization, (3) efficient and affordable tools for generating new exposure and dose data, and (4) easily accessible exposure databases. The development of new exposure models and databases along with the application of technologies like non-targeted analysis and machine learning have transformed exposure science for data-poor chemicals. By developing high-throughput tools for chemical exposure analytics and translating those tools into public health decisions ExpoCast research has served as a crucible for identifying and addressing exposure science knowledge gaps.

The chemical landscape of high-throughput new approach methodologies for exposure

The rapid characterization of risk to humans and ecosystems from exogenous chemicals requires information on both hazard and exposure. The U.S. Environmental Protection Agency's ToxCast program and the interagency Tox21 initiative have screened thousands of chemicals in various high-throughput (HT) assay systems for in vitro bioactivity. EPA's ExpoCast program is developing complementary HT methods for characterizing the human and ecological exposures necessary to interpret HT hazard data in a real-world risk context. These new approach methodologies (NAMs) for exposure include computational and analytical tools for characterizing multiple components of the complex pathways chemicals take from their source to human and ecological receptors. Here, we analyze the landscape of exposure NAMs developed in ExpoCast in the context of various chemical lists of scientific and regulatory interest, including the ToxCast and Tox21 libraries and the Toxic Substances Control Act (TSCA) inventory. We examine the landscape of traditional and exposure NAM data covering chemical use, emission, environmental fate, toxicokinetics, and ultimately external and internal exposure. We consider new chemical descriptors, machine learning models that draw inferences from existing data, high-throughput exposure models, statistical frameworks that integrate multiple model predictions, and non-targeted analytical screening methods that generate new HT monitoring information. We demonstrate that exposure NAMs drastically improve the coverage of the chemical landscape compared to traditional approaches and recommend a set of research activities to further expand the development of HT exposure data for application to risk characterization. Continuing to develop exposure NAMs to fill priority data gaps identified here will improve the availability and defensibility of risk-based metrics for use in chemical prioritization and screening. IMPACT: This analysis describes the current state of exposure assessment-based new approach methodologies across varied chemical landscapes and provides recommendations for filling key data gaps.

Effects of dioxins on animal spermatogenesis: A state-of-the-art review

The male reproductive system is especially affected by dioxins, a group of persistent environmental pollutants, resulting in irreversible abnormalities including effects on sexual function and fertility in adult males and possibly on the development of male offspring. The reproductive toxicity caused by dioxins is mostly mediated by an aryl hydrocarbon receptor (AhR). In animals, spermatogenesis is a highly sensitive and dynamic process that includes proliferation and maturation of germ cells. Spermatogenesis is subject to multiple endogenous and exogenous regulatory factors, including a wide range of environmental toxicants such as dioxins. This review discusses the toxicological effects of dioxins on spermatogenesis and their relevance to male infertility. After a detailed categorization of the environmental contaminants affecting the spermatogenesis, the exposure pathways and bioavailability of dioxins in animals was briefly reviewed. The effects of dioxins on spermatogenesis are then outlined in detail. The endocrine-disrupting effects of dioxins in animals and humans are discussed with a particular focus on their effects on the expression of spermatogenesis-related genes. Finally, the impacts of dioxins on the ratio of X and Y chromosomes, the status of serum sex hormones, the quality and fertility of sperm, and the transgenerational effects of dioxins on male reproduction are reviewed.

The development and function of the brain barriers - an overlooked consideration for chemical toxicity

It is well known that the adult brain is protected from some infections and toxic molecules by the blood-brain and the blood-cerebrospinal fluid barriers. Contrary to the immense data collected in other fields, it is deeply entrenched in environmental toxicology that xenobiotics easily permeate the developing brain because these barriers are either absent or non-functional in the fetus and newborn. Here we review the cellular and physiological makeup of the brain barrier systems in multiple species, and discuss decades of experiments that show they possess functionality during embryogenesis. We next present case studies of two chemical classes, perfluoroalkyl substances (PFAS) and bisphenols, and discuss their potential to bypass the brain barriers. While there is evidence to suggest these pollutants may enter the developing and/or adult brain parenchyma, many studies suffer from confounding technical variables which complicates data interpretation. In the future, a more formal consideration of brain barrier biology could not only improve understanding of chemical toxicokinetics but could assist in prioritizing environmental xenobiotics for their neurotoxicity risk.

Evaluation of a rapid, generic human gestational dose model

Chemical risk assessment considers potentially susceptible populations including pregnant women and developing fetuses. Humans encounter thousands of chemicals in their environments, few of which have been fully characterized. Toxicokinetic (TK) information is needed to relate chemical exposure to potentially bioactive tissue concentrations. Observational data describing human gestational exposures are unavailable for most chemicals, but physiologically based TK (PBTK) models estimate such exposures. Development of chemical-specific PBTK models requires considerable time and resources. As an alternative, generic PBTK approaches describe a standardized physiology and characterize chemicals with a set of standard physical and TK descriptors - primarily plasma protein binding and hepatic clearance. Here we report and evaluate a generic PBTK model of a human mother and developing fetus. We used a published set of formulas describing the major anatomical and physiological changes that occur during pregnancy to augment the High-Throughput Toxicokinetics (httk) software package. We simulated the ratio of concentrations in maternal and fetal plasma and compared to literature in vivo measurements. We evaluated the model with literature in vivo time-course measurements of maternal plasma concentrations in pregnant and non-pregnant women. Finally, we prioritized chemicals measured in maternal serum based on predicted fetal brain concentrations. This new model can be used for TK simulations of 859 chemicals with existing human-specific in vitro TK data as well as any new chemicals for which such data become available. This gestational model may allow for in vitro to in vivo extrapolation of point of departure doses relevant to reproductive and developmental toxicity.

Exposure Prioritization (Ex Priori): A Screening-Level High-Throughput Chemical Prioritization Tool

To estimate potential chemical risk, tools are needed to prioritize potential exposures for chemicals with minimal data. Consumer product exposures are a key pathway, and variability in consumer use patterns is an important factor. We designed Ex Priori, a flexible dashboard-type screening-level exposure model, to rapidly visualize exposure rankings from consumer product use. Ex Priori is Excel-based. Currently, it is parameterized for seven routes of exposure for 1108 chemicals present in 228 consumer product types. It includes toxicokinetics considerations to estimate body burden. It includes a simple framework for rapid modeling of broad changes in consumer use patterns by product category. Ex Priori rapidly models changes in consumer user patterns during the COVID-19 pandemic and instantly shows resulting changes in chemical exposure rankings by body burden. Sensitivity analysis indicates that the model is sensitive to the air emissions rate of chemicals from products. Ex Priori's simple dashboard facilitates dynamic exploration of the effects of varying consumer product use patterns on prioritization of chemicals based on potential exposures. Ex Priori can be a useful modeling and visualization tool to both novice and experienced exposure modelers and complement more computationally intensive population-based exposure models.