Following the voluntary phase-out of brominated flame retardants (BFRs) due to their environmental persistence and toxicity, the organophosphorus flame retardants (OPFRs) are emerging replacements. However, there is limited information on the potential human health effects of the OPFRs. Zebrafish embryos are a viable vertebrate model organism with many advantages for high throughput testing toward human hazard assessment. We utilized zebrafish embryos to assess developmental toxicity, neurotoxicity, cardiotoxicity and hepatotoxicity, of eight replacement OPFRs: (triphenyl phosphate [TPHP], isopropylated phenyl phosphate [IPP], 2-ethylhexyl diphenyl phosphate [EHDP], tert-butylated phenyl diphenyl phosphate [BPDP], trimethyl phenyl phosphate [TMPP], isodecyl diphenyl phosphate [IDDP], tris(1,3-dichloroisopropyl) phosphate [TDCIPP], and tris(2-chloroethyl) phosphate [TCEP]) and two BFRs (3,3',5,5'- tetrabromobisphenol A [TBBPA] and 2,2'4,4'-brominated diphenyl ether [BDE-47]). To determine potential effects on teratogenicity, embryos were exposed to flame retardants (FRs) at 4 h post fertilization (hpf) to 4 days post fertilization (dpf) and morphological alterations and corresponding survival were evaluated at 2 and 4 dpf. Internal concentrations were measured in larvae used in this assay by liquid chromatography-mass spectrometry. Locomotor activity was assessed in larvae treated for 48 h (from 3 dpf to 5 dpf), followed by hepatotoxicity evaluation. Finally, alterations in heart rate and rhythmicity were assessed to determine cardiotoxicity in 48 hpf embryos exposed to compounds for 3 h. Results suggest that several OPFRs (BPDP, EHDP; IPP, TMPP; TPHP and TDCIPP) produced adverse effects in multiple target organs at concentrations comparable to the two BFRs. As these OPFRs have the capacity to disrupt an integrated vertebrate model, they potentially have the capacity to affect mammalian biology. Then, we compared the lowest effective levels (LEL) in zebrafish with estimated or measured human plasma concentrations using biomonitoring data (human plasma, breast milk, handwipe samples and house dust) and a high throughput toxicokinetic (HTTK) model. Results indicate that for some compounds, the nominal LELs were within the range of human exposures, while internal LELs in zebrafish are above internal exposures in humans. These findings demonstrate the value of the zebrafish model as a relevant screening tool and support the need for further hazard characterization of the OPFRs.

A main goal of the U.S. Tox21 program is to profile a 10K-compound library for activity against a panel of stress-related and nuclear receptor signaling pathway assays using a quantitative high-throughput screening (qHTS) approach. However, assay artifacts, including nonreproducible signals and assay interference (e.g., autofluorescence), complicate compound activity interpretation. To address these issues, we have developed a data analysis pipeline that includes an updated signal noise–filtering/curation protocol and an assay interference flagging system. To better characterize various types of signals, we adopted a weighted version of the area under the curve (wAUC) to quantify the amount of activity across the tested concentration range in combination with the assay-dependent point-of-departure (POD) concentration. Based on the 32 Tox21 qHTS assays analyzed, we demonstrate that signal profiling using wAUC affords the best reproducibility (Pearson’s r = 0.91) in comparison with the POD (0.82) only or the AC₅₀ (i.e., half-maximal activity concentration, 0.81). Among the activity artifacts characterized, cytotoxicity is the major confounding factor; on average, about 8% of Tox21 compounds are affected, whereas autofluorescence affects less than 0.5%. To facilitate data evaluation, we implemented two graphical user interface applications, allowing users to rapidly evaluate the in vitro activity of Tox21 compounds.

The National Toxicology Program (NTP) receives requests to evaluate chemicals with potential to cause adverse health effects, including developmental neurotoxicity (DNT). Some recent requests have included classes of chemicals such as flame retardants, polycyclic aromatic compounds, perfluoroalkyl substances, and bisphenol A analogs with approximately 20–50 compounds per class, many of which include commercial mixtures. However, all the compounds within a class cannot be tested using traditional DNT animal testing guideline studies due to resource and time limitations. Hence, a rapid and biologically relevant screening approach is required to prioritize compounds for further in vivo testing. Because neurodevelopment is a complex process involving multiple distinct cellular processes, one assay will unlikely address the complexity. Hence, the NTP sought to characterize a battery of in vitro and alternative animal assays to quantify chemical effects on a variety of neurodevelopmental processes. A culmination of this effort resulted in a NTP-hosted collaborative project with approximately 40 participants spanning across domains of academia, industry, government, and regulatory agencies; collaborators presented data on cell-based assays and alternative animal models that was generated using a targeted set of compounds provided by the NTP. The NTP analyzed the assay results using benchmark concentration (BMC) modeling to be able to compare results across the divergent assays. The results were shared with the contributing researchers on a private web application during the workshop, and are now publicly available. This article highlights the overview and goals of the project, and describes the NTP’s approach in creating the chemical library, development of NTPs data analysis strategy, and the structure of the web application. Finally, we discuss key issues with emphasis on the utility of this approach, and knowledge gaps that need to be addressed for its use in regulatory decision making.

Cytotoxicity is a commonly used in vitro endpoint for evaluating chemical toxicity. In support of the U.S. Tox21 screening program, the cytotoxicity of ~10K chemicals was interrogated at 0, 8, 16, 24, 32, & 40 hours of exposure in a concentration dependent fashion in two cell lines (HEK293, HepG2) using two multiplexed, real-time assay technologies. One technology measures the metabolic activity of cells (i.e., cell viability, glo) while the other evaluates cell membrane integrity (i.e., cell death, flor). Using glo technology, more actives and greater temporal variations were seen in HEK293 cells, while results for the flor technology were more similar across the two cell types. Chemicals were grouped into classes based on their cytotoxicity kinetics profiles and these classes were evaluated for their associations with activity in the Tox21 nuclear receptor and stress response pathway assays. Some pathways, such as the activation of H2AX, were associated with the fast-responding cytotoxicity classes, while others, such as activation of TP53, were associated with the slow-responding cytotoxicity classes. By clustering pathways based on their degree of association to the different cytotoxicity kinetics labels, we identified clusters of pathways where active chemicals presented similar kinetics of cytotoxicity. Such linkages could be due to shared underlying biological processes between pathways, for example, activation of H2AX and heat shock factor. Others involving nuclear receptor activity are likely due to shared chemical structures rather than pathway level interactions. Based on the linkage between androgen receptor antagonism and Nrf2 activity, we surmise that a subclass of androgen receptor antagonists cause cytotoxicity via oxidative stress that is associated with Nrf2 activation. In summary, the real-time cytotoxicity screen provides informative chemical cytotoxicity kinetics data related to their cytotoxicity mechanisms, and with our analysis, it is possible to formulate mechanism-based hypotheses on the cytotoxic properties of the tested chemicals.

Over the past decade, the zebrafish is increasingly being used as a model to screen for chemical-mediated toxicities including developmental toxicity (DT) and neurotoxicity (NT). One of the major challenges is lack of harmonization in data analysis approaches, thereby posing difficulty in comparing findings across laboratories. To address this, we sought to establish a unified data analysis strategy for both DT and NT data, by adopting the benchmark concentration (BMC) analysis. There are two critical aspects in the BMC analysis: having a toxicity endpoint amenable for BMC and selecting a proper benchmark response (BMR) for the endpoint. For the former, in addition to the typical endpoints in NT assay (eg, hyper/hypo- response quantified by distance moved), we also used endpoints that assess the differences in movement patterns between chemical-treated embryos and control embryos. For the latter, we standardized the selection of BMR, which is analogous to minimum activity threshold, based on intrinsic response variations in the endpoint. When comparing our BMC results with a traditionally used LOAEL method (lowest-observed-adverse-effect level), we found high active compound concordance (100% for DT vs 74% for NT); generally, the BMC was more sensitive than LOAEL (no. of BMC more sensitive/no. of concordant active compounds, 43/50 for DT vs 16/26 for NT). Using the BMC with standardized toxicity endpoints and an appropriate BMR, we may now have a unified data-analysis approach to comparing results across different zebrafish datasets, for a better understanding of strengths and challenges when using the zebrafish as a screening tool.