Utility of In Vitro Bioactivity as a Lower Bound Estimate of In Vivo Adverse Effect Levels and in Risk-Based Prioritization

Advances in computational methods along the exposure to toxicological response paradigm

Human health risk assessment is a function of chemical toxicity, bioavailability to reach target biological tissues, and potential environmental exposure. These factors are complicated by many physiological, biochemical, physical and lifestyle factors. Furthermore, chemical health risk assessment is challenging in view of the large, and continually increasing, number of chemicals found in the environment. These challenges highlight the need to prioritize resources for the efficient and timely assessment of those environmental chemicals that pose greatest health risks. Computational methods, either predictive or investigative, are designed to assist in this prioritization in view of the lack of cost prohibitive in vivo experimental data. Computational methods provide specific and focused toxicity information using in vitro high throughput screening (HTS) assays. Information from the HTS assays can be converted to in vivo estimates of chemical levels in blood or target tissue, which in turn are converted to in vivo dose estimates that can be compared to exposure levels of the screened chemicals. This manuscript provides a review for the landscape of computational methods developed and used at the U.S. Environmental Protection Agency (EPA) highlighting their potentials and challenges.

Use of Zebrafish as a Model Organism to Study Oxidative Stress: A Review

Dioxygen is an integral part of every living organism, but its concentration varies from organ to organ. Production of metabolites from dioxygen may result in oxidative stress. Since oxidative stress has the potential to damage various biomolecules in the cell, therefore, it has presently become an active field of research. Oxidative stress has been studied in a wide range of model organisms from vertebrates to invertebrates, from rodents to piscine organisms, and from in vivo to in vitro models. But zebrafish (adults, larvae, or embryonic stage) emerged out to be the most promising vertebrate model organism to study oxidative stress because of its vast advantages (transparent embryo, cost-effectiveness, similarity to human genome, easy developmental processes, numerous offspring per spawning, and many more). This is evidenced by voluminous number of researches on oxidative stress in zebrafish exposed to chemicals, radiations, nanoparticles, pesticides, heavy metals, etc. On these backgrounds, this review attempts to highlight the potentiality of zebrafish as model of oxidative stress compared with other companion models. Several areas, from biomedical to environmental research, have been covered to explain it as a more convenient and reliable animal model for experimental research on oxidative mechanisms.

Use of new approach methodologies (NAMs) to meet regulatory requirements for the assessment of industrial chemicals and pesticides for effects on human health

New approach methodologies (NAMs) are increasingly being used for regulatory decision making by agencies worldwide because of their potential to reliably and efficiently produce information that is fit for purpose while reducing animal use. This article summarizes the ability to use NAMs for the assessment of human health effects of industrial chemicals and pesticides within the United States, Canada, and European Union regulatory frameworks. While all regulations include some flexibility to allow for the use of NAMs, the implementation of this flexibility varies across product type and regulatory scheme. This article provides an overview of various agencies' guidelines and strategic plans on the use of NAMs, and specific examples of the successful application of NAMs to meet regulatory requirements. It also summarizes intra- and inter-agency collaborations that strengthen scientific, regulatory, and public confidence in NAMs, thereby fostering their global use as reliable and relevant tools for toxicological evaluations. Ultimately, understanding the current regulatory landscape helps inform the scientific community on the steps needed to further advance timely uptake of approaches that best protect human health and the environment.

Are Non-animal Systemic Safety Assessments Protective? A Toolbox and Workflow

An important question in toxicological risk assessment is whether non-animal new approach methodologies (NAMs) can be used to make safety decisions that are protective of human health, without being overly conservative. In this work, we propose a core NAM toolbox and workflow for conducting systemic safety assessments for adult consumers. We also present an approach for evaluating how protective and useful the toolbox and workflow are by benchmarking against historical safety decisions. The toolbox includes physiologically based kinetic (PBK) models to estimate systemic Cmax levels in humans, and 3 bioactivity platforms, comprising high-throughput transcriptomics, a cell stress panel, and in vitro pharmacological profiling, from which points of departure are estimated. A Bayesian model was developed to quantify the uncertainty in the Cmax estimates depending on how the PBK models were parameterized. The feasibility of the evaluation approach was tested using 24 exposure scenarios from 10 chemicals, some of which would be considered high risk from a consumer goods perspective (eg, drugs that are systemically bioactive) and some low risk (eg, existing food or cosmetic ingredients). Using novel protectiveness and utility metrics, it was shown that up to 69% (9/13) of the low risk scenarios could be identified as such using the toolbox, whilst being protective against all (5/5) the high-risk ones. The results demonstrated how robust safety decisions could be made without using animal data. This work will enable a full evaluation to assess how protective and useful the toolbox and workflow are across a broader range of chemical-exposure scenarios.

Grouping of UVCB substances with dose-response transcriptomics data from human cell-based assays

The application of in vitro biological assays as new approach methodologies (NAMs) to support grouping of UVCB (unknown or variable composition, complex reaction products, and biological materials) substances has recently been demonstrated. In addition to cell-based phenotyping as NAMs, in vitro transcriptomic profiling is used to gain deeper mechanistic understanding of biological responses to chemicals and to support grouping and read-across. However, the value of gene expression profiling for characterizing complex substances like UVCBs has not been explored. Using 141 petroleum substance extracts, we performed dose-response transcriptomic profiling in human induced pluripotent stem cell (iPSC)-derived hepatocytes, cardiomyocytes, neurons, and endothelial cells, as well as cell lines MCF7 and A375. The goal was to determine whether transcriptomic data can be used to group these UVCBs and to further characterize the molecular basis for in vitro biological responses. We found distinct transcriptional responses for petroleum substances by manufacturing class. Pathway enrichment informed interpretation of effects of substances and UVCB petroleum-class. Transcriptional activity was strongly correlated with concentration of polycyclic aromatic compounds (PAC), especially in iPSC-derived hepatocytes. Supervised analysis using transcriptomics, alone or in combination with bioactivity data collected on these same substances/cells, suggest that transcriptomics data provide useful mechanistic information, but only modest additional value for grouping. Overall, these results further demonstrate the value of NAMs for grouping of UVCBs, identify informative cell lines, and provide data that could be used for justifying selection of substances for further testing that may be required for registration.

A comparison of the lowest effective concentration in culture media for detection of chromosomal damage in vitro and in blood or plasma for detection of micronuclei in vivo

It is often assumed that genotoxic substances will be detected more easily by using in vitro rather than in vivo genotoxicity tests since higher concentrations, more cytotoxicity and static exposures can be achieved. However, there is a paucity of data demonstrating whether genotoxic substances are detected at lower concentrations in cell culture in vitro than can be reached in the blood of animals treated in vivo. To investigate this issue, we compared the lowest concentration required for induction of chromosomal damage in vitro (lowest observed effective concentration, or LOEC) with the concentration of the test substance in blood at the lowest dose required for biologically relevant induction of micronuclei in vivo (lowest observed effective dose, or LOED). In total, 83 substances were found for which the LOED could be identified or estimated, where concentrations in blood and micronucleus data were available via the same route of administration in the same species, and in vitro chromosomal damage data were available. 39.8 % of substances were positive in vivo at blood concentrations that were lower than the LOEC in vitro, 22.9 % were positive at similar concentrations, and 37.3 % of substances were positive in vivo at higher concentrations. Distribution analysis showed a very wide scatter of > 6 orders of magnitude across these 3 categories. When mode of action was evaluated, the distribution of clastogens and aneugens across the 3 categories was very similar. Thus, the ability to detect induction of micronuclei in bone marrow in vivo regardless of the mechanism for micronucleus induction, is clearly not solely determined by the concentration of test substance which induced chromosomal damage in vitro.

Comprehensive interpretation of in vitro micronucleus test results for 292 chemicals: from hazard identification to risk assessment application

Risk assessments are increasingly reliant on information from in vitro assays. The in vitro micronucleus test (MNvit) is a genotoxicity test that detects chromosomal abnormalities, including chromosome breakage (clastogenicity) and/or whole chromosome loss (aneugenicity). In this study, MNvit datasets for 292 chemicals, generated by the US EPA's ToxCast program, were evaluated using a decision tree-based pipeline for hazard identification. Chemicals were tested with 19 concentrations (n = 1) up to 200 µM, in the presence and absence of Aroclor 1254-induced rat liver S9. To identify clastogenic chemicals, %MN values at each concentration were compared to a distribution of batch-specific solvent controls; this was followed by cytotoxicity assessment and benchmark concentration (BMC) analyses. The approach classified 157 substances as positives, 25 as negatives, and 110 as inconclusive. Using the approach described in Bryce et al. (Environ Mol Mutagen 52:280-286, 2011), we identified 15 (5%) aneugens. IVIVE (in vitro to in vivo extrapolation) was employed to convert BMCs into administered equivalent doses (AEDs). Where possible, AEDs were compared to points of departure (PODs) for traditional genotoxicity endpoints; AEDs were generally lower than PODs based on in vivo endpoints. To facilitate interpretation of in vitro MN assay concentration-response data for risk assessment, exposure estimates were utilized to calculate bioactivity exposure ratio (BER) values. BERs for 50 clastogens and two aneugens had AEDs that approached exposure estimates (i.e., BER < 100); these chemicals might be considered priorities for additional testing. This work provides a framework for the use of high-throughput in vitro genotoxicity testing for priority setting and chemical risk assessment.

The Evolution of Regulatory Toxicology: Where is the Gardener?

There is a need for paradigm change in the methodology employed for toxicological testing and assessment. It could be said that this change is well on its way, through an evolutionary progress analogous to that of natural selection. Darwin's Theory of Evolution has defined the idea of evolution and descendancy since the last third of the 19th century. Increasingly, this concept of 'evolution' is being applied beyond the field of biology. This Comment article discusses the progress of toxicological testing in the context of 'evolutionary pressure' and deliberates how this process can help foster the development, implementation and acceptance of mechanistic and human-relevant methods in this field. By comparing the current regulatory landscape in toxicity testing and assessment to specific elements in Charles Darwin's evolutionary theory, we aim to better understand the needs and requirements for the future.

Next Generation Risk Assessment of the Anti-Androgen Flutamide Including the Contribution of Its Active Metabolite Hydroxyflutamide

In next generation risk assessment (NGRA), non-animal approaches are used to quantify the chemical concentrations required to trigger bioactivity responses, in order to assure safe levels of human exposure. A limitation of many in vitro bioactivity assays, which are used in an NGRA context as new approach methodologies (NAMs), is that toxicokinetics, including biotransformation, are not adequately captured. The present study aimed to include, as a proof of principle, the bioactivity of the metabolite hydroxyflutamide (HF) in an NGRA approach to evaluate the safety of the anti-androgen flutamide (FLU), using the AR-CALUX assay to derive the NAM point of departure (PoD). The NGRA approach applied also included PBK modelling-facilitated quantitative in vitro to in vivo extrapolation (QIVIVE). The PBK model describing FLU and HF kinetics in humans was developed using GastroPlus™ and validated against human pharmacokinetic data. PBK model-facilitated QIVIVE was performed to translate the in vitro AR-CALUX derived concentration-response data to a corresponding in vivo dose-response curve for the anti-androgenicity of FLU, excluding and including the activity of HF (-HF and +HF, respectively). The in vivo benchmark dose 5% lower confidence limits (BMDL05) derived from the predicted in vivo dose-response curves for FLU, revealed a 440-fold lower BMDL05 when taking the bioactivity of HF into account. Subsequent comparison of the predicted BMDL05 values to the human therapeutic doses and historical animal derived PoDs, revealed that PBK modelling-facilitated QIVIVE that includes the bioactivity of the active metabolite is protective and provides a more appropriate PoD to assure human safety via NGRA, whereas excluding this would potentially result in an underestimation of the risk of FLU exposure in humans.

Impact of High-Throughput Model Parameterization and Data Uncertainty on Thyroid-Based Toxicological Estimates for Pesticide Chemicals

Chemical-induced alteration of maternal thyroid hormone levels may increase the risk of adverse neurodevelopmental outcomes in offspring. US federal risk assessments rely almost exclusively on apical endpoints in animal models for deriving points of departure (PODs). New approach methodologies (NAMs) such as high-throughput screening (HTS) and mechanistically informative in vitro human cell-based systems, combined with in vitro to in vivo extrapolation (IVIVE), supplement in vivo studies and provide an alternative approach to calculate/determine PODs. We examine how parameterization of IVIVE models impacts the comparison between IVIVE-derived equivalent administered doses (EADs) from thyroid-relevant in vitro assays and the POD values that serve as the basis for risk assessments. Pesticide chemicals with thyroid-based in vitro bioactivity data from the US Tox21 HTS program were included (n = 45). Depending on the model structure used for IVIVE analysis, up to 35 chemicals produced EAD values lower than the POD. A total of 10 chemicals produced EAD values higher than the POD regardless of the model structure. The relationship between IVIVE-derived EAD values and the in vivo-derived POD values is highly dependent on model parameterization. Here, we derive a range of potentially thyroid-relevant doses that incorporate uncertainty in modeling choices and in vitro assay data.

IVIVE: Facilitating the Use of In Vitro Toxicity Data in Risk Assessment and Decision Making

During the past few decades, the science of toxicology has been undergoing a transformation from observational to predictive science. New approach methodologies (NAMs), including in vitro assays, in silico models, read-across, and in vitro to in vivo extrapolation (IVIVE), are being developed to reduce, refine, or replace whole animal testing, encouraging the judicious use of time and resources. Some of these methods have advanced past the exploratory research stage and are beginning to gain acceptance for the risk assessment of chemicals. A review of the recent literature reveals a burst of IVIVE publications over the past decade. In this review, we propose operational definitions for IVIVE, present literature examples for several common toxicity endpoints, and highlight their implications in decision-making processes across various federal agencies, as well as international organizations, including those in the European Union (EU). The current challenges and future needs are also summarized for IVIVE. In addition to refining and reducing the number of animals in traditional toxicity testing protocols and being used for prioritizing chemical testing, the goal to use IVIVE to facilitate the replacement of animal models can be achieved through their continued evolution and development, including a strategic plan to qualify IVIVE methods for regulatory acceptance.

Latent Variables Capture Pathway-Level Points of Departure in High-Throughput Toxicogenomic Data

Estimation of points of departure (PoDs) from high-throughput transcriptomic data (HTTr) represents a key step in the development of next-generation risk assessment (NGRA). Current approaches mainly rely on single key gene targets, which are constrained by the information currently available in the knowledge base and make interpretation challenging as scientists need to interpret PoDs for thousands of genes or hundreds of pathways. In this work, we aimed to address these issues by developing a computational workflow to investigate the pathway concentration-response relationships in a way that is not fully constrained by known biology and also facilitates interpretation. We employed the Pathway-Level Information ExtractoR (PLIER) to identify latent variables (LVs) describing biological activity and then investigated in vitro LVs' concentration-response relationships using the ToxCast pipeline. We applied this methodology to a published transcriptomic concentration-response data set for 44 chemicals in MCF-7 cells and showed that our workflow can capture known biological activity and discriminate between estrogenic and antiestrogenic compounds as well as activity not aligning with the existing knowledge base, which may be relevant in a risk assessment scenario. Moreover, we were able to identify the known estrogen activity in compounds that are not well-established ER agonists/antagonists supporting the use of the workflow in read-across. Next, we transferred its application to chemical compounds tested in HepG2, HepaRG, and MCF-7 cells and showed that PoD estimates are in strong agreement with those estimated using a recently developed Bayesian approach (cor = 0.89) and in weak agreement with those estimated using a well-established approach such as BMDExpress2 (cor = 0.57). These results demonstrate the effectiveness of using PLIER in a concentration-response scenario to investigate pathway activity in a way that is not fully constrained by the knowledge base and to ease the biological interpretation and support the development of an NGRA framework with the ability to improve current risk assessment strategies for chemicals using new approach methodologies.

Benchmark Dose Modeling Approaches for Volatile Organic Chemicals using a Novel Air-Liquid Interface In Vitro Exposure System

Inhalation is the most relevant route of volatile organic chemical (VOC) exposure; however, due to unique challenges posed by their chemical properties and poor solubility in aqueous solutions, in vitro chemical safety testing is predominantly performed using direct application dosing/submerged exposures. To address the difficulties in screening toxic effects of VOCs, our cell culture exposure system permits cells to be exposed to multiple concentrations at air-liquid interface (ALI) in a 24-well format. ALI exposure methods permit direct chemical-to-cell interaction with the test article at physiological conditions. In the present study, BEAS-2B and primary normal human bronchial epithelial cells (pHBEC) are used to assess gene expression, cytotoxicity, and cell viability responses to a variety of volatile chemicals including acrolein, formaldehyde, 1,3-butadiene, acetaldehyde, 1-bromopropane, carbon tetrachloride, dichloromethane, and trichloroethylene. BEAS-2B cells were exposed to all the test agents, while pHBECs were only exposed to the latter four listed above. The VOC concentrations tested elicited only slight cell viability changes in both cell types. Gene expression changes were analyzed using benchmark dose (BMD) modeling. The BMD for the most sensitive gene set was within one order of magnitude of the threshold-limit value reported by the American Conference of Governmental Industrial Hygienists, and the most sensitive gene sets impacted by exposure correlate to known adverse health effects recorded in epidemiologic and in vivo exposure studies. Overall, our study outlines a novel in vitro approach for evaluating molecular-based points-of-departure in human airway epithelial cell exposure to volatile chemicals.

Application of an Accessible Interface for Pharmacokinetic Modeling and In Vitro to In Vivo Extrapolation

Regulatory toxicology testing has traditionally relied on in vivo methods to inform decision-making. However, scientific, practical, and ethical considerations have led to an increased interest in the use of in vitro and in silico methods to fill data gaps. While in vitro experiments have the advantage of rapid application across large chemical sets, interpretation of data coming from these non-animal methods can be challenging due to the mechanistic nature of many assays. In vitro to in vivo extrapolation (IVIVE) has emerged as a computational tool to help facilitate this task. Specifically, IVIVE uses physiologically based pharmacokinetic (PBPK) models to estimate tissue-level chemical concentrations based on various dosing parameters. This approach is used to estimate the administered dose needed to achieve in vitro bioactivity concentrations within the body. IVIVE results can be useful to inform on metrics such as margin of exposure or to prioritize potential chemicals of concern, but the PBPK models used in this approach have extensive data requirements. Thus, access to input parameters, as well as the technical requirements of applying and interpreting models, has limited the use of IVIVE as a routine part of in vitro testing. As interest in using non-animal methods for regulatory and research contexts continues to grow, our perspective is that access to computational support tools for PBPK modeling and IVIVE will be essential for facilitating broader application and acceptance of these techniques, as well as for encouraging the most scientifically sound interpretation of in vitro results. We highlight recent developments in two open-access computational support tools for PBPK modeling and IVIVE accessible via the Integrated Chemical Environment (https://ice.ntp.niehs.nih.gov/), demonstrate the types of insights these tools can provide, and discuss how these analyses may inform in vitro-based decision making.

Risk-Based Prioritization of Organic Chemicals and Locations of Ecological Concern in Sediment From Great Lakes Tributaries

With improved analytical techniques, environmental monitoring studies are increasingly able to report the occurrence of tens or hundreds of chemicals per site, making it difficult to identify the most relevant chemicals from a biological standpoint. For the present study, organic chemical occurrence was examined, individually and as mixtures, in the context of potential biological effects. Sediment was collected at 71 Great Lakes (USA/Canada) tributary sites and analyzed for 87 chemicals. Multiple risk-based lines of evidence were used to prioritize chemicals and locations, including comparing sediment concentrations and estimated porewater concentrations with established whole-organism benchmarks (i.e., sediment and water quality criteria and screening values) and with high-throughput toxicity screening data from the US Environmental Protection Agency's ToxCast database, estimating additive effects of chemical mixtures on common ToxCast endpoints, and estimating toxic equivalencies for mixtures of alkylphenols and polycyclic aromatic hydrocarbons (PAHs). This multiple-lines-of-evidence approach enabled the screening of more chemicals, mitigated the uncertainties of individual approaches, and strengthened common conclusions. Collectively, at least one benchmark/screening value was exceeded for 54 of the 87 chemicals, with exceedances observed at all 71 of the monitoring sites. Chemicals with the greatest potential for biological effects, both individually and as mixture components, were bisphenol A, 4-nonylphenol, indole, carbazole, and several PAHs. Potential adverse outcomes based on ToxCast gene targets and putative adverse outcome pathways relevant to individual chemicals and chemical mixtures included tumors, skewed sex ratios, reproductive dysfunction, hepatic steatosis, and early mortality, among others. The results provide a screening-level prioritization of chemicals with the greatest potential for adverse biological effects and an indication of sites where they are most likely to occur. Environ Toxicol Chem 2022;41:1016-1041. 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.

Canadian Regulatory Perspective on Next Generation Risk Assessments for Pest Control Products and Industrial Chemicals

In 2012, the Council of Canadian Academies published the expert panel on integrated testing of pesticide’s report titled: Integrating emerging technologies into chemical safety assessment. This report was prepared for the Government of Canada in response to a request from the Minister of Health and on behalf of the Pest Management Regulatory Agency. It examined the scientific status of the use of integrated testing strategies for the regulatory health risk assessment of pesticides while noting the data-rich/poor dichotomy that exists when comparing pesticide formulations to most industrial chemicals. It also noted that the adoption of integrated approaches to testing and assessment (IATA) strategies may refine and streamline testing of chemicals, as well as improve results in the future. Moreover, the experts expected to see an increase in the use of integrated testing strategies over the next decade, resulting in improved evidence-based decision-making. Subsequent to this report, there has been great advancements in IATA strategies, which includes the incorporation of adverse outcome pathways (AOPs) and new approach methodologies (NAMs). This perspective provides the first Canadian regulatory update on how Health Canada is also advancing the incorporation of alternative, non-animal strategies, using a weight of evidence approach, for the evaluation of pest control products and industrial chemicals. It will include specific initiatives and describe how this work is leading to the creation of next generation risk assessments. It also reflects Health Canada’s commitment towards implementing the 3Rs of animal testing: reduce, refine and replace the need for animal studies, whenever possible.

Beyond AOPs: A Mechanistic Evaluation of NAMs in DART Testing

New Approach Methodologies (NAMs) promise to offer a unique opportunity to enable human-relevant safety decisions to be made without the need for animal testing in the context of exposure-driven Next Generation Risk Assessment (NGRA). Protecting human health against the potential effects a chemical may have on embryo-foetal development and/or aspects of reproductive biology using NGRA is particularly challenging. These are not single endpoint or health effects and risk assessments have traditionally relied on data from Developmental and Reproductive Toxicity (DART) tests in animals. There are numerous Adverse Outcome Pathways (AOPs) that can lead to DART, which means defining and developing strict testing strategies for every AOP, to predict apical outcomes, is neither a tenable goal nor a necessity to ensure NAM-based safety assessments are fit-for-purpose. Instead, a pragmatic approach is needed that uses the available knowledge and data to ensure NAM-based exposure-led safety assessments are sufficiently protective. To this end, the mechanistic and biological coverage of existing NAMs for DART were assessed and gaps to be addressed were identified, allowing the development of an approach that relies on generating data relevant to the overall mechanisms involved in human reproduction and embryo-foetal development. Using the knowledge of cellular processes and signalling pathways underlying the key stages in reproduction and development, we have developed a broad outline of endpoints informative of DART. When the existing NAMs were compared against this outline to determine whether they provide comprehensive coverage when integrated in a framework, we found them to generally cover the reproductive and developmental processes underlying the traditionally evaluated apical endpoint studies. The application of this safety assessment framework is illustrated using an exposure-led case study.

Innovation in regulatory approaches for endocrine disrupting chemicals: The journey to risk assessment modernization in Canada

Globally, regulatory authorities grapple with the challenge of assessing the hazards and risks to human and ecosystem health that may result from exposure to chemicals that disrupt the normal functioning of endocrine systems. Rapidly increasing number of chemicals in commerce, coupled with the reliance on traditional, costly animal experiments for hazard characterization - often with limited sensitivity to many important mechanisms of endocrine disruption -, presents ongoing challenges for chemical regulation. The consequence is a limited number of chemicals for which there is sufficient data to assess if there is endocrine toxicity and hence few chemicals with thorough hazard characterization. To address this challenge, regulatory assessment of endocrine disrupting chemicals (EDCs) is benefiting from a revolution in toxicology that focuses on New Approach Methodologies (NAMs) to more rapidly identify, prioritize, and assess the potential risks from exposure to chemicals using novel, more efficient, and more mechanistically driven methodologies and tools. Incorporated into Integrated Approaches to Testing and Assessment (IATA) and guided by conceptual frameworks such as Adverse Outcome Pathways (AOPs), emerging approaches focus initially on molecular interactions between the test chemical and potentially vulnerable biological systems instead of the need for animal toxicity data. These new toxicity testing methods can be complemented with in silico and computational toxicology approaches, including those that predict chemical kinetics. Coupled with exposure data, these will inform risk-based decision-making approaches. Canada is part of a global network collaborating on building confidence in the use of NAMs for regulatory assessment of EDCs. Herein, we review the current approaches to EDC regulation globally (mainly from the perspective of human health), and provide a perspective on how the advances for regulatory testing and assessment can be applied and discuss the promises and challenges faced in adopting these novel approaches to minimize risks due to EDC exposure in Canada, and our world.

Advancing the Adverse Outcome Pathway for PPARγ Inactivation Leading to Pulmonary Fibrosis Using Bradford-Hill Consideration and the Comparative Toxicogenomics Database

Pulmonary fibrosis is regulated by transforming growth factor-β (TGF-β) and peroxisome proliferator-activated receptor-gamma (PPARγ). An adverse outcome pathway (AOP) for PPARγ inactivation leading to pulmonary fibrosis has been previously developed. To advance the development of this AOP, the confidence of the overall AOP was assessed using the Bradford-Hill considerations as per the recommendations from the Organisation for Economic Co-operation and Development (OECD) Users' Handbook. Overall, the essentiality of key events (KEs) and the biological plausibility of key event relationships (KERs) were rated high. In contrast, the empirical support of KERs was found to be moderate. To experimentally evaluate the KERs from the molecular initiating event (MIE) and KE1, PPARγ (MIE) and TGF-β (KE1) inhibitors were used to examine the effects of downstream events following inhibition of their upstream events. PPARγ inhibition (MIE) led to TGF-β activation (KE1), upregulation in vimentin expression (KE3), and an increase in the fibronectin level (KE4). Similarly, activated TGF-β (KE1) led to an increase in vimentin (KE3) and fibronectin expression (KE4). In the database analysis using the Comparative Toxicogenomics Database, 31 genes related to each KE were found to be highly correlated with pulmonary fibrosis, and the top 21 potential stressors were suggested. The AOP for pulmonary fibrosis evaluated in this study will be the basis for the screening of inhaled toxic substances in the environment.

A framework for chemical safety assessment incorporating new approach methodologies within REACH

The long-term investment in new approach methodologies (NAMs) within the EU and other parts of the world is beginning to result in an emerging consensus of how to use information from in silico, in vitro and targeted in vivo sources to assess the safety of chemicals. However, this methodology is being adopted very slowly for regulatory purposes. Here, we have developed a framework incorporating in silico, in vitro and in vivo methods designed to meet the requirements of REACH in which both hazard and exposure can be assessed using a tiered approach. The outputs from each tier are classification categories, safe doses, and risk assessments, and progress through the tiers depends on the output from previous tiers. We have exemplified the use of the framework with three examples. The outputs were the same or more conservative than parallel assessments based on conventional studies. The framework allows a transparent and phased introduction of NAMs in chemical safety assessment and enables science-based safety decisions which provide the same level of public health protection using fewer animals, taking less time, and using less financial and expert resource. Furthermore, it would also allow new methods to be incorporated as they develop through continuous selective evolution rather than periodic revolution.

Food, Beverage, and Feedstock Processing Facility Wastewater: a Unique and Underappreciated Source of Contaminants to U.S. Streams

Process wastewaters from food, beverage, and feedstock facilities, although regulated, are an under-investigated environmental contaminant source. Food process wastewaters (FPWWs) from 23 facilities in 17 U.S. states were sampled and documented for a plethora of chemical and microbial contaminants. Of the 576 analyzed organics, 184 (32%) were detected at least once, with concentrations as large as 143 μg L^-1 (6:2 fluorotelomer sulfonic acid), and as many as 47 were detected in a single FPWW sample. Cumulative per/polyfluoroalkyl substance concentrations up to 185 μg L^-1 and large pesticide transformation product concentrations (e.g., methomyl oxime, 40 μg L^-1; clothianidin TMG, 2.02 μg L^-1) were observed. Despite 48% of FPWW undergoing disinfection treatment prior to discharge, bacteria resistant to third-generation antibiotics were found in each facility type, and multiple bacterial groups were detected in all samples, including total coliforms. The exposure-activity ratios and toxicity quotients exceeded 1.0 in 13 and 22% of samples, respectively, indicating potential biological effects and toxicity to vertebrates and invertebrates associated with the discharge of FPWW. Organic contaminant profiles of FPWW differed from previously reported contaminant profiles of municipal effluents and urban storm water, indicating that FPWW is another important source of chemical and microbial contaminant mixtures discharged into receiving surface waters.

Quantitative non-targeted analysis: Bridging the gap between contaminant discovery and risk characterization

Chemical risk assessments follow a long-standing paradigm that integrates hazard, dose-response, and exposure information to facilitate quantitative risk characterization. Targeted analytical measurement data directly support risk assessment activities, as well as downstream risk management and compliance monitoring efforts. Yet, targeted methods have struggled to keep pace with the demands for data regarding the vast, and growing, number of known chemicals. Many contemporary monitoring studies therefore utilize non-targeted analysis (NTA) methods to screen for known chemicals with limited risk information. Qualitative NTA data has enabled identification of previously unknown compounds and characterization of data-poor compounds in support of hazard identification and exposure assessment efforts. In spite of this, NTA data have seen limited use in risk-based decision making due to uncertainties surrounding their quantitative interpretation. Significant efforts have been made in recent years to bridge this quantitative gap. Based on these advancements, quantitative NTA data, when coupled with other high-throughput data streams and predictive models, are poised to directly support 21st-century risk-based decisions. This article highlights components of the chemical risk assessment process that are influenced by NTA data, surveys the existing literature for approaches to derive quantitative estimates of chemicals from NTA measurements, and presents a conceptual framework for incorporating NTA data into contemporary risk assessment frameworks.

Implementing in vitro bioactivity data to modernize priority setting of chemical inventories

Internationally, there are thousands of existing and newly introduced chemicals in commerce, highlighting the ongoing importance of innovative approaches to identify emerging chemicals of concern. For many chemicals, there is a paucity of hazard and exposure data. Thus, there is a crucial need for efficient and robust approaches to address data gaps and support risk-based prioritization. Several studies have demonstrated the utility of in vitro bioactivity data from the ToxCast program in deriving points of departure (PODs). ToxCast contains data for nearly 1,400 endpoints per chemical, and the bioactivity concentrations, indicative of potential adverse outcomes, can be converted to human-equivalent PODs using high-throughput toxicokinetics (HTTK) modeling. However, data gaps need to be addressed for broader application: the limited chemical space of HTTK and quantitative high-throughput screening data. Here we explore the applicability of in silico models to address these data needs. Specifically, we used ADMET predictor for HTTK predictions and a generalized read-across approach to predict ToxCast bioactivity potency. We applied these models to profile 5,801 chemicals on Canada’s Domestic Substances List (DSL). To evaluate the approach’s performance, bioactivity PODs were compared with in vivo results from the EPA Toxicity Values database for 1,042 DSL chemicals. Comparisons demonstrated that the bioac­tivity PODs, based on ToxCast data or read-across, were conservative for 95% of the chemicals. Comparing bioactivity PODs to human exposure estimates supports the identification of chemicals of potential interest for further work. The bioac­tivity workflow shows promise as a powerful screening tool to support effective triaging of chemical inventories.

Potential Human Health Hazard of Post-Hurricane Harvey Sediments in Galveston Bay and Houston Ship Channel: A Case Study of Using In Vitro Bioactivity Data to Inform Risk Management Decisions

Natural and anthropogenic disasters may be associated with redistribution of chemical contaminants in the environment; however, current methods for assessing hazards and risks of complex mixtures are not suitable for disaster response. This study investigated the suitability of in vitro toxicity testing methods as a rapid means of identifying areas of potential human health concern. We used sediment samples (n = 46) from Galveston Bay and the Houston Ship Channel (GB/HSC) areas after hurricane Harvey, a disaster event that led to broad redistribution of chemically-contaminated sediments, including deposition of the sediment on shore due to flooding. Samples were extracted with cyclohexane and dimethyl sulfoxide and screened in a compendium of human primary or induced pluripotent stem cell (iPSC)-derived cell lines from different tissues (hepatocytes, neuronal, cardiomyocytes, and endothelial) to test for concentration-dependent effects on various functional and cytotoxicity phenotypes (n = 34). Bioactivity data were used to map areas of potential concern and the results compared to the data on concentrations of polycyclic aromatic hydrocarbons (PAHs) in the same samples. We found that setting remediation goals based on reducing bioactivity is protective of both "known" risks associated with PAHs and "unknown" risks associated with bioactivity, but the converse was not true for remediation based on PAH risks alone. Overall, we found that in vitro bioactivity can be used as a comprehensive indicator of potential hazards and is an example of a new approach method (NAM) to inform risk management decisions on site cleanup.

A new approach method for characterizing inter-species toxicodynamic variability

Inter-species differences in toxicodynamics are often a critical source of uncertainty in safety evaluations and typically dealt with using default adjustment factors. In vitro studies that use cells from different species demonstrated some success for estimating the relationships between life span and/or body weight and sensitivity to cytotoxicity; however, no apparent investigation evaluated the utility of these models for risk assessment. It was hypothesized that an in vitro model using dermal fibroblasts derived from diverse species and individuals might be utilized to inform the extent of inter-species and inter-individual variability in toxicodynamics. To test this hypothesis and characterize both inter-species and inter-individual variability in cytotoxicity, concentration-response cytotoxicity screening of 40 chemicals in primary dermal fibroblasts from 68 individuals of 54 diverse species was conducted. Chemicals examined included drugs, environmental pollutants, and food/flavor/fragrance agents; most of these were previously assessed either in vivo or in vitro for inter-species or inter-individual variation. Species included humans, the typical preclinical species and representatives from other orders of mammals and birds. Data demonstrated that both inter-species and inter-individual components of variability contribute to the observed differences in sensitivity to cell death. Further, it was found that the magnitude of the observed inter-species and inter-individual differences was chemical-dependent. This study contributes to the paradigm shift in risk assessment from reliance on in vivo toxicity testing to higher-throughput in vitro or alternative approaches, extending the strategy to replace use of default adjustment factors with experimental characterization of toxicodynamic inter-individual variability and to also address toxicodynamic inter-species variability.

Human Chemical Exposure from Background Emissions in the United States and the Implication for Quantifying Risks from Marginal Emission Increase

The linear dose-response relationship has long been assumed in assessments of health risk from an incremental chemical emission relative to background emissions. In this study, we systematically examine the relevancy of such an assumption with real-world data. We used the reported emission data, as background emissions, from the 2017 U.S. National Emission Inventory for 95 organic chemicals to estimate the central tendencies of exposures of the general U.S. population. Previously published nonlinear dose-response relationships for chemicals were used to estimate health risk from exposure. We also explored and identified four intervals of exposure in which the nonlinear dose-response relationship may be linearly approximated with fixed slopes. Predicted rates of exposure to these 95 chemicals are all within the lowest of the four intervals and associated with low health risk. The health risk may be overestimated if a slope on the dose-response relationship extrapolated from toxicological assays based on high response rates is used for a marginal increase in emission not substantially higher than background emissions. To improve the confidence of human health risk estimates for chemicals, future efforts should focus on deriving a more accurate dose-response relationship at lower response rates and interface it with exposure assessments.

In silico approaches in organ toxicity hazard assessment: current status and future needs in predicting liver toxicity

Hepatotoxicity is one of the most frequently observed adverse effects resulting from exposure to a xenobiotic. For example, in pharmaceutical research and development it is one of the major reasons for drug withdrawals, clinical failures, and discontinuation of drug candidates. The development of faster and cheaper methods to assess hepatotoxicity that are both more sustainable and more informative is critically needed. The biological mechanisms and processes underpinning hepatotoxicity are summarized and experimental approaches to support the prediction of hepatotoxicity are described, including toxicokinetic considerations. The paper describes the increasingly important role of in silico approaches and highlights challenges to the adoption of these methods including the lack of a commonly agreed upon protocol for performing such an assessment and the need for in silico solutions that take dose into consideration. A proposed framework for the integration of in silico and experimental information is provided along with a case study describing how computational methods have been used to successfully respond to a regulatory question concerning non-genotoxic impurities in chemically synthesized pharmaceuticals.

Incorporating human exposure information in a weight of evidence approach to inform design of repeated dose animal studies

Human health risks from chronic exposures to environmental chemicals are typically estimated from potential human exposure estimates and dose-response data obtained from repeated-dose animal toxicity studies. Various criteria are available for selecting the top (highest) dose used in these animal studies. For example, toxicokinetic (TK) and toxicological data provided by shorter-term or dose range finding studies can be evaluated in a weight of evidence approach to provide insight into the dose range that would provide dose-response data that are relevant to human exposures. However, there are concerns that a top dose resulting from the consideration of TK data may be too low compared to other criteria, such as the limit dose or the maximum tolerated dose. In this paper, we address several concerns related to human exposures by discussing 1) the resources and methods available to predict human exposure levels and the associated uncertainty and variability, and 2) the margin between predicted human exposure levels and the dose levels used in repeated-dose animal studies. A series of case studies, ranging from data-rich to data-poor chemicals, are presented to demonstrate that expected human exposures to environmental chemicals are typically orders of magnitude lower than no-observed-adverse-effect levels/lowest-observed-adverse-effect levels (NOAELs/LOAELs) when available (used as conservative surrogates for top doses). The results of these case studies support that a top dose based, in part, on TK data is typically orders of magnitude higher than expected human exposure levels.

Opportunities and challenges related to saturation of toxicokinetic processes: Implications for risk assessment

Top dose selection for repeated dose animal studies has generally focused on identification of apical endpoints, use of the limit dose, or determination of a maximum tolerated dose (MTD). The intent is to optimize the ability of toxicity tests performed in a small number of animals to detect effects for hazard identification. An alternative approach, the kinetically derived maximum dose (KMD), has been proposed as a mechanism to integrate toxicokinetic (TK) data into the dose selection process. The approach refers to the dose above which the systemic exposures depart from being proportional to external doses. This non-linear external-internal dose relationship arises from saturation or limitation of TK process(es), such as absorption or metabolism. The importance of TK information is widely acknowledged when assessing human health risks arising from exposures to environmental chemicals, as TK determines the amount of chemical at potential sites of toxicological responses. However, there have been differing opinions and interpretations within the scientific and regulatory communities related to the validity and application of the KMD concept. A multi-stakeholder working group, led by the Health and Environmental Sciences Institute (HESI), was formed to provide an opportunity for impacted stakeholders to address commonly raised scientific and technical issues related to this topic and, more specifically, a weight of evidence approach is recommended to inform design and dose selection for repeated dose animal studies. Commonly raised challenges related to the use of TK data for dose selection are discussed, recommendations are provided, and illustrative case examples are provided to address these challenges or refute misconceptions.

Potency Ranking of Per- and Polyfluoroalkyl Substances Using High-Throughput Transcriptomic Analysis of Human Liver Spheroids

Per- and polyfluoroalkyl substances (PFAS) are some of the most prominent organic contaminants in human blood. Although the toxicological implications of human exposure to perfluorooctane sulfonate (PFOS) and perfluorooctanoate (PFOA) are well established, data on lesser-understood PFAS are limited. New approach methodologies (NAMs) that apply bioinformatic tools to high-throughput data are being increasingly considered to inform risk assessment for data-poor chemicals. The aim of this study was to compare the potencies (ie, benchmark concentrations: BMCs) of PFAS in primary human liver microtissues (3D spheroids) using high-throughput transcriptional profiling. Gene expression changes were measured using TempO-seq, a templated, multiplexed RNA-sequencing platform. Spheroids were exposed for 1 or 10 days to increasing concentrations of 23 PFAS in 3 subgroups: carboxylates (PFCAs), sulfonates (PFSAs), and fluorotelomers and sulfonamides. PFCAs and PFSAs exhibited trends toward increased transcriptional potency with carbon chain-length. Specifically, longer-chain compounds (7-10 carbons) were more likely to induce changes in gene expression and have lower transcriptional BMCs. The combined high-throughput transcriptomic and bioinformatic analyses support the capability of NAMs to efficiently assess the effects of PFAS in liver microtissues. The data enable potency ranking of PFAS for human liver cell spheroid cytotoxicity and transcriptional changes, and assessment of in vitro transcriptomic points of departure. These data improve our understanding of the possible health effects of PFAS and will be used to inform read-across for human health risk assessment.

Potency Ranking of Per- and Polyfluoroalkyl Substances Using High-Throughput Transcriptomic Analysis of Human Liver Spheroids

Per- and polyfluoroalkyl substances (PFAS) are some of the most prominent organic contaminants in human blood. Although the toxicological implications of human exposure to perfluorooctane sulfonate (PFOS) and perfluorooctanoate (PFOA) are well established, data on lesser-understood PFAS are limited. New approach methodologies (NAMs) that apply bioinformatic tools to high-throughput data are being increasingly considered to inform risk assessment for data-poor chemicals. The aim of this study was to compare the potencies (ie, benchmark concentrations: BMCs) of PFAS in primary human liver microtissues (3D spheroids) using high-throughput transcriptional profiling. Gene expression changes were measured using TempO-seq, a templated, multiplexed RNA-sequencing platform. Spheroids were exposed for 1 or 10 days to increasing concentrations of 23 PFAS in 3 subgroups: carboxylates (PFCAs), sulfonates (PFSAs), and fluorotelomers and sulfonamides. PFCAs and PFSAs exhibited trends toward increased transcriptional potency with carbon chain-length. Specifically, longer-chain compounds (7-10 carbons) were more likely to induce changes in gene expression and have lower transcriptional BMCs. The combined high-throughput transcriptomic and bioinformatic analyses support the capability of NAMs to efficiently assess the effects of PFAS in liver microtissues. The data enable potency ranking of PFAS for human liver cell spheroid cytotoxicity and transcriptional changes, and assessment of in vitro transcriptomic points of departure. These data improve our understanding of the possible health effects of PFAS and will be used to inform read-across for human health risk assessment.

Potency Ranking of Per- and Polyfluoroalkyl Substances Using High-Throughput Transcriptomic Analysis of Human Liver Spheroids

Per- and polyfluoroalkyl substances (PFAS) are some of the most prominent organic contaminants in human blood. Although the toxicological implications of human exposure to perfluorooctane sulfonate (PFOS) and perfluorooctanoate (PFOA) are well established, data on lesser-understood PFAS are limited. New approach methodologies (NAMs) that apply bioinformatic tools to high-throughput data are being increasingly considered to inform risk assessment for data-poor chemicals. The aim of this study was to compare the potencies (ie, benchmark concentrations: BMCs) of PFAS in primary human liver microtissues (3D spheroids) using high-throughput transcriptional profiling. Gene expression changes were measured using TempO-seq, a templated, multiplexed RNA-sequencing platform. Spheroids were exposed for 1 or 10 days to increasing concentrations of 23 PFAS in 3 subgroups: carboxylates (PFCAs), sulfonates (PFSAs), and fluorotelomers and sulfonamides. PFCAs and PFSAs exhibited trends toward increased transcriptional potency with carbon chain-length. Specifically, longer-chain compounds (7-10 carbons) were more likely to induce changes in gene expression and have lower transcriptional BMCs. The combined high-throughput transcriptomic and bioinformatic analyses support the capability of NAMs to efficiently assess the effects of PFAS in liver microtissues. The data enable potency ranking of PFAS for human liver cell spheroid cytotoxicity and transcriptional changes, and assessment of in vitro transcriptomic points of departure. These data improve our understanding of the possible health effects of PFAS and will be used to inform read-across for human health risk assessment.

Extrapolating from acute to chronic toxicity in vitro

Chemical safety assessment requires information on both chronic and acute effects of toxicants. Traditionally, such information has been provided by a set of animal studies conducted over different durations, ranging from a single dose with observation of effects over a few days, to repeat daily dosing and observations made over many months. With the advent of modern mechanistic approaches to toxicology, the role of in vitro studies within alternative approaches has never been more prominent. Typical in vitro experiments are conducted over short durations with measurements of response at a single time point, with a focus on providing effect and concentration-response information as input to hazard and risk assessment. This limits the usefulness of such data since potential chronic effects that cumulate over time are not usually considered. To address this, an experimental design is presented to characterise the toxicodynamics of a response not only in terms of concentration, but also as a function of time. Generation of concentration-time-effect responses allows both the extrapolation of points of departure from an acute to chronic exposure, and the determination of a chronicity index that provides a quantitative measure of a chemical's potential to cause cumulative effects over time. In addition, the approach provides a means to characterise the dynamics of key event relationships for the development of quantitative adverse outcome pathways.

Utilizing High-Throughput Screening Data, Integrative Toxicological Prioritization Index Score, and Exposure-Activity Ratios for Chemical Prioritization: A Case Study of Endocrine-Active Pesticides in Food Crops

Endocrine-active chemicals can directly act on nuclear receptors and trigger the disturbances of metabolism and a homeostatic system, which are important risk factors for complicating chronic diseases in humans. The endocrine-active potentials of pesticides acting on estrogen, androgen, and thyroid hormone receptors have been extensively evaluated for pesticides; however, the effects on other receptors are less understood. This study aims to comprehensively characterize and prioritize the endocrine-active pesticides using an exposure-activity ratio (EAR) method and toxicological prioritization index (ToxPi). The aggregate exposure assessment of pesticides was performed using a computational exposure model [stochastic human exposure and dose simulation high-throughput model (SHEDS-HT)]. Minimum in vitro point of departure values were converted to human oral equivalent doses via in vitro-to-in vivo extrapolation. The overall endocrine-disrupting potentials of pesticides were evaluated via 76 assays, representing 11 nuclear receptors. EARs and ToxPi scores were then derived to prioritize 79 pesticides in food. This case study demonstrates that EAR profiling can inform the regulatory agencies for a relevant chemical prioritization, which would direct in-depth health risk assessments in the future.

Upholding the EU's Commitment to 'Animal Testing as a Last Resort' Under REACH Requires a Paradigm Shift in How We Assess Chemical Safety to Close the Gap Between Regulatory Testing and Modern Safety Science

Animal use for testing chemicals under REACH continues to increase, despite advances in non-animal safety science during the past 15 years. The application of modern science and technology, and the use of 'next generation' weight-of-evidence assessment approaches, are embedded in EU guidance for establishing the safety of cosmetics and foods - and of the ingredients used in these products. However, this is still not the case for the regulation of chemicals. Under the new Chemicals Strategy for Sustainability, thought leaders in human health and environmental protection are calling on the European Commission to quickly embrace the benefits of modern and innovative non-animal safety science, in place of outdated animal testing, if the EU is to be a leader in safe and sustainable innovation under the European Green Deal transformational change ambitions. The European Commission also needs to enable companies to meet their legal obligation to only conduct animal testing as a last resort, by providing a more flexible, science-based and consistent regulatory framework for assuring chemical safety, which supports the integration of data from different sources. We are at a tipping point for closing the gap between regulatory chemicals testing and modern safety science. It is time to join forces, across policy makers, scientists, regulators and lawyers, to lead the paradigm shift needed to deliver what EU citizens want - namely, chemicals and products that are safe and sustainable, without resorting to animal testing.

Sourcing data on chemical properties and hazard data from the US-EPA CompTox Chemicals Dashboard: A practical guide for human risk assessment

For the past six decades, human health risk assessment of chemicals has relied on in vivo data from human epidemiological and experimental animal toxicological studies to inform the derivation of non-cancer toxicity values. The ongoing evolution of this risk assessment paradigm in an environmental landscape of data-poor chemicals has highlighted the need to develop and implement non-testing methods, so-called New Approach Methodologies (NAMs). NAMs include a growing number of in silico and in vitro data streams designed to inform hazard properties of chemicals, including kinetics and dynamics at different levels of biological organization, environmental fate and transport, and exposure. NAMs provide a fit-for-purpose science-basis for human hazard and risk characterization of chemicals ranging from data-gap filling applications to broad evidence-based decision-making. Systematic assembly and delivery of empirical and predicted data for chemicals are paramount to advancing chemical evaluation, and software tools serve an essential role in delivering these data to the scientific community. The CompTox Chemicals Dashboard (from here on referred to as the "Dashboard") is one such tool and is a publicly available web-based application developed by the US Environmental Protection Agency to provide access to chemistry, toxicity and exposure information for ~900,000 chemicals. The Dashboard is increasingly becoming a valuable resource for assessors tasked with the evaluation of potential human health risks associated with chemical exposures. In this context, the significant amount of information present in the Dashboard facilitates: 1) assembly of information on physicochemical properties and environmental fate and transport and exposure parameters and metrics; 2) identification of cancer and non-cancer health effects from extant human and experimental animal studies in the public domain and/or information not available in the public domain (i.e., "grey literature"); 3) systematic literature searching and review for developing cancer and non-cancer hazard evidence bases; and 4) access to mechanistic information that can aid or augment the analysis of traditional toxicology evidence bases, or potentially, serve as the primary basis for informing hazard identification and dose-response when traditional bioassay data are lacking. Finally, in silico predictive tools developed to conduct structure-activity or read-across analyses are also available within the Dashboard. This practical tutorial is intended to address key questions from the human health risk assessment community dealing with chemicals in both food and in the environment. Perspectives for future development or refinement of the Dashboard highlight foreseen activities to further support the research and risk assessment community in cancer and non-cancer chemical evaluations.

High-Throughput Transcriptomic Analysis of Human Primary Hepatocyte Spheroids Exposed to Per- and Polyfluoroalkyl Substances as a Platform for Relative Potency Characterization

Per- and poly-fluoroalkyl substances (PFAS) are widely found in the environment because of their extensive use and persistence. Although several PFAS are well studied, most lack toxicity data to inform human health hazard and risk assessment. This study focused on 4 model PFAS: perfluorooctanoic acid (PFOA; 8 carbon), perfluorobutane sulfonate (PFBS; 4 carbon), perfluorooctane sulfonate (PFOS; 8 carbon), and perfluorodecane sulfonate (PFDS; 10 carbon). Human primary liver cell spheroids (pooled from 10 donors) were exposed to 10 concentrations of each PFAS and analyzed at 4 time points. The approach aimed to: (1) identify gene expression changes mediated by the PFAS, (2) identify similarities in biological responses, (3) compare PFAS potency through benchmark concentration analysis, and (4) derive bioactivity exposure ratios (ratio of the concentration at which biological responses occur, relative to daily human exposure). All PFAS induced transcriptional changes in cholesterol biosynthesis and lipid metabolism pathways, and predicted PPARα activation. PFOS exhibited the most transcriptional activity and had a highly similar gene expression profile to PFDS. PFBS induced the least transcriptional changes and the highest benchmark concentration (ie, was the least potent). The data indicate that these PFAS may have common molecular targets and toxicities, but that PFOS and PFDS are the most similar. The transcriptomic bioactivity exposure ratios derived here for PFOA and PFOS were comparable to those derived using rodent apical endpoints in risk assessments. These data provide a baseline level of toxicity for comparison with other known PFAS using this testing strategy.

Risk Characterization of Environmental Samples Using In Vitro Bioactivity and Polycyclic Aromatic Hydrocarbon Concentrations Data

Methods to assess environmental exposure to hazardous chemicals have primarily focused on quantification of individual chemicals, although chemicals often occur in mixtures, presenting challenges to the traditional risk characterization framework. Sampling sites in a defined geographic region provide an opportunity to characterize chemical contaminants, with spatial interpolation as a tool to provide estimates for non-sampled sites. At the same time, the use of in vitro bioactivity measurements has been shown to be informative for rapid risk-based decisions. In this study, we measured in vitro bioactivity in 39 surface soil samples collected immediately after flooding associated with Hurricane Harvey in Texas in a residential area known to be inundated with polycyclic aromatic hydrocarbon (PAH) contaminants. Bioactivity data were from a number of functional and toxicity assays in 5 human cell types, such as induced pluripotent stem cell-derived hepatocytes, cardiomyocytes, neurons, and endothelial cells, as well as human umbilical vein endothelial cells. Data on concentrations of PAH in these samples were also available and the combination of data sources offered a unique opportunity to assess the joint spatial variation of PAH components and bioactivity. We found significant evidence of spatial correlation of a subset of PAH contaminants and of cell-based phenotypes. In addition, we show that the cell-based bioactivity data can be used to predict environmental concentrations for several PAH contaminants, as well as overall PAH summaries and cancer risk. This study's impact lies in its demonstration that cell-based profiling can be used for rapid hazard screening of environmental samples by anchoring the bioassays to concentrations of PAH. This work sets the stage for identification of the areas of concern and direct quantitative risk characterization based on bioactivity data, thereby providing an important supplement to traditional individual chemical analyses by shedding light on constituents that may be missed from targeted chemical monitoring.

Chemical Characterization of Recycled Consumer Products Using Suspect Screening Analysis

Recycled materials are found in many consumer products as part of a circular economy; however, the chemical content of recycled products is generally uncharacterized. A suspect screening analysis using two-dimensional gas chromatography time-of-flight mass spectrometry (GC × GC-TOFMS) was applied to 210 products (154 recycled, 56 virgin) across seven categories. Chemicals in products were tentatively identified using a standard spectral library or confirmed using chemical standards. A total of 918 probable chemical structures identified (112 of which were confirmed) in recycled materials versus 587 (110 confirmed) in virgin materials. Identified chemicals were characterized in terms of their functional use and structural class. Recycled paper products and construction materials contained greater numbers of chemicals than virgin products; 733 identified chemicals had greater occurrence in recycled compared to virgin materials. Products made from recycled materials contained greater numbers of fragrances, flame retardants, solvents, biocides, and dyes. The results were clustered to identify groups of chemicals potentially associated with unique chemical sources, and identified chemicals were prioritized for further study using high-throughput hazard and exposure information. While occurrence is not necessarily indicative of risk, these results can be used to inform the expansion of existing models or identify exposure pathways currently neglected in exposure assessments.

Paving the way for application of next generation risk assessment to safety decision-making for cosmetic ingredients

Next generation risk assessment (NGRA) is an exposure-led, hypothesis-driven approach that has the potential to support animal-free safety decision-making. However, significant effort is needed to develop and test the in vitro and in silico (computational) approaches that underpin NGRA to enable confident application in a regulatory context. A workshop was held in Montreal in 2019 to discuss where effort needs to be focussed and to agree on the steps needed to ensure safety decisions made on cosmetic ingredients are robust and protective. Workshop participants explored whether NGRA for cosmetic ingredients can be protective of human health, and reviewed examples of NGRA for cosmetic ingredients. From the limited examples available, it is clear that NGRA is still in its infancy, and further case studies are needed to determine whether safety decisions are sufficiently protective and not overly conservative. Seven areas were identified to help progress application of NGRA, including further investments in case studies that elaborate on scenarios frequently encountered by industry and regulators, including those where a ‘high risk’ conclusion would be expected. These will provide confidence that the tools and approaches can reliably discern differing levels of risk. Furthermore, frameworks to guide performance and reporting should be developed.

High-throughput PBTK models for in vitro to in vivo extrapolation

INTRODUCTION

Toxicity data are unavailable for many thousands of chemicals in commerce and the environment. Therefore, risk assessors need to rapidly screen these chemicals for potential risk to public health. High-throughput screening (HTS) for in vitro bioactivity, when used with high-throughput toxicokinetic (HTTK) data and models, allows characterization of these thousands of chemicals.

AREAS COVERED

This review covers generic physiologically based toxicokinetic (PBTK) models and high-throughput PBTK modeling for in vitro-in vivo extrapolation (IVIVE) of HTS data. We focus on 'httk', a public, open-source set of computational modeling tools and in vitro toxicokinetic (TK) data.

EXPERT OPINION

HTTK benefits chemical risk assessors with its ability to support rapid chemical screening/prioritization, perform IVIVE, and provide provisional TK modeling for large numbers of chemicals using only limited chemical-specific data. Although generic TK model design can increase prediction uncertainty, these models provide offsetting benefits by increasing model implementation accuracy. Also, public distribution of the models and data enhances reproducibility. For the httk package, the modular and open-source design can enable the tool to be used and continuously improved by a broad user community in support of the critical need for high-throughput chemical prioritization and rapid dose estimation to facilitate rapid hazard assessments.

Estimating Hepatotoxic Doses Using High-Content Imaging in Primary Hepatocytes

Using in vitro data to estimate point of departure (POD) values is an essential component of new approach methodologies (NAM)-based chemical risk assessments. In this case study, we evaluated a NAM for hepatotoxicity based on rat primary hepatocytes, high-content imaging (HCI), and toxicokinetic modeling. First, we treated rat primary hepatocytes with 10 concentrations (0.2 to 100 µM) of 51 chemicals that produced hepatotoxicity in repeat-dose subchronic and chronic exposures. Second, we used HCI to measure endoplasmic reticulum stress, mitochondrial function, lysosomal mass, steatosis, apoptosis, DNA texture, nuclear size, and cell number at 24, 48, and 72 h and calculated concentrations at 50% maximal activity (AC50). Third, we estimated administered equivalent doses (AEDs) from AC50 values using toxicokinetic modeling. AEDs using physiologically-based toxicokinetic models were 4.1-fold (SD 6.3) and 8.1-fold (SD 15.5) lower than subchronic and chronic lowest observed adverse effect levels (LOAELs), respectively. In contrast, AEDs from ToxCast and Tox21 assays were 89.8-fold (SD 149.5) and 168-fold (SD 323.7) lower than subchronic and chronic LOAELs. Individual HCI end-points also estimated AEDs for specific hepatic lesions that were lower than in vivo PODs. Lastly, AEDs were similar for different in vitro exposure durations, but steady-state toxicokinetic models produced 7.6-fold lower estimates than dynamic physiologically-based ones. Our findings suggest that NAMs from diverse cell types provide conservative estimates of PODs. In contrast, NAMs based on the same species and cell type as the adverse outcome may produce estimates closer to the traditional in vivo PODs.

Quantitative in Vitro-to-in Vivo Extrapolation for Mixtures: A Case Study of Superfund Priority List Pesticides

In vitro cell-based toxicity testing methods generate large amounts of data informative for risk-based evaluations. To allow extrapolation of the quantitative outputs from cell-based tests to the equivalent exposure levels in humans, reverse toxicokinetic (RTK) modeling is used to conduct in vitro-to-in vivo extrapolation (IVIVE) from in vitro effective concentrations to in vivo oral dose equivalents. IVIVE modeling approaches for individual chemicals are well-established; however, the potential implications of chemical-to-chemical interactions in mixture settings on IVIVE remains largely unexplored. We hypothesized that chemical co-exposures could modulate both protein binding efficiency and hepatocyte clearance of the chemicals in a mixture, which would in turn affect the quantitative IVIVE toxicokinetic parameters. To test this hypothesis, we used 20 pesticides from the Agency for Toxic Substances and Disease Registry (ATSDR) Substance Priority List, both individually and as equimolar mixtures, and investigated the concentration-dependent effects of chemical interactions on in vitro toxicokinetic parameters. Plasma protein binding efficiency was determined by using ultracentrifugation, and hepatocyte clearance was estimated in suspensions of cryopreserved primary human hepatocytes. We found that for single chemicals, the protein binding efficiencies were similar at different test concentrations. In a mixture, however, both protein binding efficiency and hepatocyte clearance were affected. When IVIVE was conducted using mixture-derived toxicokinetic data, more conservative estimates of Activity-to-Exposure Ratios (AERs) were produced as compared to using data from single chemical experiments. Because humans are exposed to mixtures of chemicals, this study is significant as it demonstrates the importance of incorporating mixture-derived parameters into IVIVE for in vitro bioactivity data in order to accurately prioritize risks and facilitate science-based decision-making.

Predictive modeling of biological responses in the rat liver using in vitro Tox21 bioactivity: Benefits from high-throughput toxicokinetics

Computational methods are needed to more efficiently leverage data from in vitro cell-based models to predict what occurs within whole body systems after chemical insults. This study set out to test the hypothesis that in vitro high-throughput screening (HTS) data can more effectively predict in vivo biological responses when chemical disposition and toxicokinetic (TK) modeling are employed. In vitro HTS data from the Tox21 consortium were analyzed in concert with chemical disposition modeling to derive nominal, aqueous, and intracellular estimates of concentrations eliciting 50% maximal activity. In vivo biological responses were captured using rat liver transcriptomic data from the DrugMatrix and TG-Gates databases and evaluated for pathway enrichment. In vivo dosing data were translated to equivalent body concentrations using HTTK modeling. Random forest models were then trained and tested to predict in vivo pathway-level activity across 221 chemicals using in vitro bioactivity data and physicochemical properties as predictor variables, incorporating methods to address imbalanced training data resulting from high instances of inactivity. Model performance was quantified using the area under the receiver operator characteristic curve (AUC-ROC) and compared across pathways for different combinations of predictor variables. All models that included toxicokinetics were found to outperform those that excluded toxicokinetics. Biological interpretation of the model features revealed that rather than a direct mapping of in vitro assays to in vivo pathways, unexpected combinations of multiple in vitro assays predicted in vivo pathway-level activities. To demonstrate the utility of these findings, the highest-performing model was leveraged to make new predictions of in vivo biological responses across all biological pathways for remaining chemicals tested in Tox21 with adequate data coverage (n = 6617). These results demonstrate that, when chemical disposition and toxicokinetics are carefully considered, in vitro HT screening data can be used to effectively predict in vivo biological responses to chemicals.

High-Throughput Transcriptomics Platform for Screening Environmental Chemicals

New approach methodologies (NAMs) that efficiently provide information about chemical hazard without using whole animals are needed to accelerate the pace of chemical risk assessments. Technological advancements in gene expression assays have made in vitro high-throughput transcriptomics (HTTr) a feasible option for NAMs-based hazard characterization of environmental chemicals. In this study, we evaluated the Templated Oligo with Sequencing Readout (TempO-Seq) assay for HTTr concentration-response screening of a small set of chemicals in the human-derived MCF7 cell model. Our experimental design included a variety of reference samples and reference chemical treatments in order to objectively evaluate TempO-Seq assay performance. To facilitate analysis of these data, we developed a robust and scalable bioinformatics pipeline using open-source tools. We also developed a novel gene expression signature-based concentration-response modeling approach and compared the results to a previously implemented workflow for concentration-response analysis of transcriptomics data using BMDExpress. Analysis of reference samples and reference chemical treatments demonstrated highly reproducible differential gene expression signatures. In addition, we found that aggregating signals from individual genes into gene signatures prior to concentration-response modeling yielded in vitro transcriptional biological pathway altering concentrations (BPACs) that were closely aligned with previous ToxCast high-throughput screening assays. Often these identified signatures were associated with the known molecular target of the chemicals in our test set as the most sensitive components of the overall transcriptional response. This work has resulted in a novel and scalable in vitro HTTr workflow that is suitable for high-throughput hazard evaluation of environmental chemicals.

Elucidation of the Effects of Bisphenol A and Structural Analogs on Germ and Steroidogenic Cells Using Single Cell High-Content Imaging

Concerns about the potential adverse effects of bisphenol A (BPA) have led to an increase in the use of replacements, yet the toxicity data for several of these chemicals are limited. Using high-content imaging, we compared the effects of BPA, BPAF, BPF, BPS, BPM, and BPTMC in germ (C18-4 spermatogonial) and steroidogenic (MA-10 Leydig and KGN granulosa) cell lines. Effects on cell viability and phenotypic markers were analyzed to determine benchmark concentrations (BMCs) and estimate administered equivalent doses (AEDs). In all 3 cell lines, BPA was one of the least cytotoxic bisphenol compounds tested, whereas BPM and BPTMC were the most cytotoxic. Interestingly, BPF and BPS were cytotoxic only in MA-10 cells. Effects on phenotypic parameters, including mitochondria, lysosomes, lipid droplets, and oxidative stress, were both bisphenol- and cell-line specific. BPA exposure affected mitochondria (BMC: 1.2 μM; AED: 0.09 mg/kg/day) in C18-4 cells. Lysosome numbers were increased in MA-10 cells exposed to BPA or BPAF but decreased in KGN cells exposed to BPAF or BPM. Lipid droplets were decreased in C18-4 cells exposed to BPF and in MA-10 cells exposed to BPTMC but increased in BPF, BPM, and BPTMC-exposed KGN cells. BPA and BPM exposure induced oxidative stress in MA-10 and KGN cells, respectively. In summary, structurally similar bisphenols displayed clear cell-line-specific differences in BMC and AED values for effects on cell viability and phenotypic endpoints. This approach, together with additional data on human exposure, may aid in the selection and prioritization of responsible replacements for BPA. .

FutureTox IV Workshop Summary: Predictive Toxicology for Healthy Children

FutureTox IV, a Society of Toxicology Contemporary Concepts in Toxicology workshop, was held in November 2018. Building upon FutureTox I, II, and III, this conference focused on the latest science and technology for in vitro profiling and in silico modeling as it relates to predictive developmental and reproductive toxicity (DART). Publicly available high-throughput screening data sets are now available for broad in vitro profiling of bioactivities across large inventories of chemicals. Coupling this vast amount of mechanistic data with a deeper understanding of molecular embryology and post-natal development lays the groundwork for using new approach methodologies (NAMs) to evaluate chemical toxicity, drug efficacy, and safety assessment for embryo-fetal development. NAM is a term recently adopted in reference to any technology, methodology, approach, or combination thereof that can be used to provide information on chemical hazard and risk assessment to avoid the use of intact animals (U.S. Environmental Protection Agency [EPA], Strategic plan to promote the development and implementation of alternative test methods within the tsca program, 2018, https://www.epa.gov/sites/production/files/2018-06/documents/epa_alt_strat_plan_6-20-18_clean_final.pdf). There are challenges to implementing NAMs to evaluate chemicals for developmental toxicity compared with adult toxicity. This forum article reviews the 2018 workshop activities, highlighting challenges and opportunities for applying NAMs for adverse pregnancy outcomes (eg, preterm labor, malformations, low birth weight) as well as disorders manifesting postnatally (eg, neurodevelopmental impairment, breast cancer, cardiovascular disease, fertility). DART is an important concern for different regulatory statutes and test guidelines. Leveraging advancements in such approaches and the accompanying efficiencies to detecting potential hazards to human development are the unifying concepts toward implementing NAMs in DART testing. Although use of NAMs for higher level regulatory decision making is still on the horizon, the conference highlighted novel testing platforms and computational models that cover multiple levels of biological organization, with the unique temporal dynamics of embryonic development, and novel approaches for estimating toxicokinetic parameters essential in supporting in vitro to in vivo extrapolation.

Image-based profiling for drug discovery: due for a machine-learning upgrade?

Image-based profiling is a maturing strategy by which the rich information present in biological images is reduced to a multidimensional profile, a collection of extracted image-based features. These profiles can be mined for relevant patterns, revealing unexpected biological activity that is useful for many steps in the drug discovery process. Such applications include identifying disease-associated screenable phenotypes, understanding disease mechanisms and predicting a drug's activity, toxicity or mechanism of action. Several of these applications have been recently validated and have moved into production mode within academia and the pharmaceutical industry. Some of these have yielded disappointing results in practice but are now of renewed interest due to improved machine-learning strategies that better leverage image-based information. Although challenges remain, novel computational technologies such as deep learning and single-cell methods that better capture the biological information in images hold promise for accelerating drug discovery.