The increased use of flammable plastics and electronic devices along with stricter fire safety standards has led to the heavy use of flame retardant chemicals in many consumer, commercial, and industrial products. Although flame retardant use has increased, a great deal of uncertainty surrounds their safety with some evidence showing toxicity and risk to human and environmental health. Recent efforts have focused on designing high-throughput biological platforms with nonmammalian models to evaluate and prioritize chemicals with limited hazard information. To complement these efforts, this study used a new morphological and behavioral testing platform with embryonic zebrafish to characterize the developmental toxicity of 44 halogenated and organophosphate flame retardants, including several of their known metabolites. Zebrafish were exposed to flame retardants from 6 to 120 h post fertilization (hpf) across concentrations spanning 4 orders of magnitude (eg, 6.4 nM to 64 µM). Flame retardant effects on survival and development were evaluated at 24 and 120 hpf, and neurobehavioral changes were measured using 2 photomotor response (PMR) assays. Compared to controls, 93% (41/44) of flame retardants studied elicited adverse effects among one or more of the bioassays and concentrations tested with the aryl phosphate ester (APE)-based mono-isopropylated triaryl phosphate and the brominated-bisphenol-A analog tetrabromobisphenol-A producing the greatest array of malformations. Hierarchical clustering showed that APE flame retardants with isopropyl, butyl, and cresyl substituents on phenyl rings clustered tightly and were particularly potent. Both PMR assays were highly predictive of morphological defects supporting their use as nonlethal means of evaluating teratogenicity that could allow for additional evaluations of long-term or delayed effects in older animals. Taken together, evidence presented here indicates that zebrafish neurodevelopment is highly sensitive to many flame retardants currently in use and can be used to understand potential vulnerabilities to human health.

New approach methodologies (NAMs) that efficiently provide information about chemical hazard without using whole animals are needed to accelerate the pace of chemical safety 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 the present 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 (HTS) 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.

The U.S. Environmental Protection Agency Endocrine Disruptor Screening Program and the Organization for Economic Co-operation and Development (OECD) have used the human adrenocarcinoma (H295R) cell-based assay to predict chemical perturbation of androgen and estrogen production. Recently, a high-throughput H295R (HT-H295R) assay was developed as part of the ToxCast program that includes measurement of 11 hormones, including progestagens, corticosteroids, androgens, and estrogens. To date, 2012 chemicals have been screened at 1 concentration; of these, 656 chemicals have been screened in concentration-response. The objectives of this work were to: (1) develop an integrated analysis of chemical-mediated effects on steroidogenesis in the HT-H295R assay and (2) evaluate whether the HT-H295R assay predicts estrogen and androgen production specifically via comparison with the OECD-validated H295R assay. To support application of HT-H295R assay data to weight-of-evidence and prioritization tasks, a single numeric value based on Mahalanobis distances was computed for 654 chemicals to indicate the magnitude of effects on the synthesis of 11 hormones. The maximum mean Mahalanobis distance (maxmMd) values were high for strong modulators (prochloraz, mifepristone) and lower for moderate modulators (atrazine, molinate). Twenty-five of 28 reference chemicals used for OECD validation were screened in the HT-H295R assay, and produced qualitatively similar results, with accuracies of 0.90/0.75 and 0.81/0.91 for increased/decreased testosterone and estradiol production, respectively. The HT-H295R assay provides robust information regarding estrogen and androgen production, as well as additional hormones. The maxmMd from this integrated analysis may provide a data-driven approach to prioritizing lists of chemicals for putative effects on steroidogenesis.