Benchmark Concentrations for Untargeted Metabolomics Versus Transcriptomics for Liver Injury Compounds in In Vitro Liver Models.

Interpretation of untargeted metabolomics data from both in vivo and physiologically relevant in vitro model systems continues to be a significant challenge for toxicology research. Potency-based modeling of toxicological responses has served as a pillar of interpretive context and translation of testing data. In this study, we leverage the resolving power of concentration-response modeling through benchmark concentration (BMC) analysis to interpret untargeted metabolomics data from differentiated cultures of HepaRG cells exposed to a panel of reference compounds and integrate data in a potency-aligned framework with matched transcriptomic data. For this work, we characterized biological responses to classical human liver injury compounds and comparator compounds, known to not cause liver injury in humans, at 10 exposure concentrations in spent culture media by untargeted liquid chromatography-mass spectrometry analysis. The analyte features observed (with limited metabolites identified) were analyzed using BMC modeling to derive compound-induced points of departure. The results revealed liver injury compounds produced concentration-related increases in metabolomic response compared to those rarely associated with liver injury (ie, sucrose, potassium chloride). Moreover, the distributions of altered metabolomic features were largely comparable with those observed using high throughput transcriptomics, which were further extended to investigate the potential for in vitro observed biological responses to be observed in humans with exposures at therapeutic doses. These results demonstrate the utility of BMC modeling of untargeted metabolomics data as a sensitive and quantitative indicator of human liver injury potential.

Evaluating Sufficient Similarity of Botanical Dietary Supplements: Combining Chemical and In Vitro Biological Data.

Botanical dietary supplements are complex mixtures with numerous potential sources of variation along the supply chain from raw plant material to the market. Approaches for determining sufficient similarity (ie, complex mixture read-across) may be required to extrapolate efficacy or safety data from a tested sample to other products containing the botanical ingredient(s) of interest. In this work, screening-level approaches for generating both chemical and biological-response profiles were used to evaluate the similarity of black cohosh (Actaea racemosa) and Echinacea purpurea samples to well-characterized National Toxicology Program (NTP) test articles. Data from nontargeted chemical analyses and gene expression of toxicologically important hepatic receptor pathways (aryl hydrocarbon receptor [AhR], constitutive androstane receptor [CAR], pregnane X receptor [PXR], farnesoid X receptor [FXR], and peroxisome proliferator-activated receptor alpha [PPARα]) in primary human hepatocyte cultures were used to determine similarity through hierarchical clustering. Although there were differences in chemical profiles across black cohosh samples, these differences were not reflected in the biological-response profiles. These findings highlight the complexity of biological-response dynamics that may not be reflected in chemical composition profiles. Thus, biological-response data could be used as the primary basis for determining similarity among black cohosh samples. Samples of E. purpurea displayed better correlation in similarity across chemical and biological-response measures. The general approaches described herein can be applied to complex mixtures with unidentified active constituents to determine when data from a tested mixture (eg, NTP test article) can be used for hazard identification of sufficiently similar mixtures, with the knowledge of toxicological targets informing assay selection when possible.

The Power of Resolution: Contextualized Understanding of Biological Responses to Liver Injury Chemicals Using High-throughput Transcriptomics and Benchmark Concentration Modeling.

Prediction of human response to chemical exposures is a major challenge in both pharmaceutical and toxicological research. Transcriptomics has been a powerful tool to explore chemical-biological interactions, however, limited throughput, high-costs, and complexity of transcriptomic interpretations have yielded numerous studies lacking sufficient experimental context for predictive application. To address these challenges, we have utilized a novel high-throughput transcriptomics (HTT) platform, TempO-Seq, to apply the interpretive power of concentration-response modeling with exposures to 24 reference compounds in both differentiated and non-differentiated human HepaRG cell cultures. Our goals were to (1) explore transcriptomic characteristics distinguishing liver injury compounds, (2) assess impacts of differentiation state of HepaRG cells on baseline and compound-induced responses (eg, metabolically-activated), and (3) identify and resolve reference biological-response pathways through benchmark concentration (BMC) modeling. Study data revealed the predictive utility of this approach to identify human liver injury compounds by their respective BMCs in relation to human internal exposure plasma concentrations, and effectively distinguished drug analogs with varied associations of human liver injury (eg, withdrawn therapeutics trovafloxacin and troglitazone). Impacts of cellular differentiation state (proliferated vs differentiated) were revealed on baseline drug metabolizing enzyme expression, hepatic receptor signaling, and responsiveness to metabolically-activated toxicants (eg, cyclophosphamide, benzo(a)pyrene, and aflatoxin B1). Finally, concentration-response modeling enabled efficient identification and resolution of plausibly-relevant biological-response pathways through their respective pathway-level BMCs. Taken together, these findings revealed HTT paired with differentiated in vitro liver models as an effective tool to model, explore, and interpret toxicological and pharmacological interactions.

Comprehensive Analyses and Prioritization of Tox21 10K Chemicals Affecting Mitochondrial Function by in-Depth Mechanistic Studies.

BACKGROUND: A central challenge in toxicity testing is the large number of chemicals in commerce that lack toxicological assessment. In response, the Tox21 program is re-focusing toxicity testing from animal studies to less expensive and higher throughput in vitro methods using target/pathway-specific, mechanism-driven assays. OBJECTIVES: Our objective was to use an in-depth mechanistic study approach to prioritize and characterize the chemicals affecting mitochondrial function. METHODS: We used a tiered testing approach to prioritize for more extensive testing 622 compounds identified from a primary, quantitative high-throughput screen of 8,300 unique small molecules, including drugs and industrial chemicals, as potential mitochondrial toxicants by their ability to significantly decrease the mitochondrial membrane potential (MMP). Based on results from secondary MMP assays in HepG2 cells and rat hepatocytes, 34 compounds were selected for testing in tertiary assays that included formation of reactive oxygen species (ROS), upregulation of p53 and nuclear erythroid 2-related factor 2/antioxidant response element (Nrf2/ARE), mitochondrial oxygen consumption, cellular Parkin translocation, and larval development and ATP status in the nematode Caenorhabditis elegans. RESULTS: A group of known mitochondrial complex inhibitors (e.g., rotenone) and uncouplers (e.g., chlorfenapyr), as well as potential novel complex inhibitors and uncouplers, were detected. From this study, we identified four not well-characterized potential mitochondrial toxicants (lasalocid, picoxystrobin, pinacyanol, and triclocarban) that merit additional in vivo characterization. CONCLUSIONS: The tier-based approach for identifying and mechanistically characterizing mitochondrial toxicants can potentially reduce animal use in toxicological testing. https://doi.org/10.1289/EHP2589.

Editor's Highlight: Comparative Toxicity of Organophosphate Flame Retardants and Polybrominated Diphenyl Ethers to Caenorhabditis elegans.

With the phasing-out of the polybrominated diphenyl ether (PBDE) flame retardants due to concerns regarding their potential developmental toxicity, the use of replacement compounds such as organophosphate flame retardants (OPFRs) has increased. Limited toxicity data are currently available to estimate the potential adverse health effects of the OPFRs. The toxicological effects of 4 brominated flame retardants, including 3 PBDEs and 3,3',5,5'-tetrabromobisphenol A, were compared with 6 aromatic OPFRs and 2 aliphatic OPFRs. The effects of these chemicals were determined using 3 biological endpoints in the nematode Caenorhabditis elegans (feeding, larval development, and reproduction). Because C. elegans development was previously reported to be sensitive to mitochondrial function, results were compared with those from an in vitro mitochondrial membrane permeabilization (MMP) assay. Overall 11 of the 12 flame retardants were active in 1 or more C. elegans biological endpoints, with only tris(2-chloroethyl) phosphate inactive across all endpoints including the in vitro MMP assay. For 2 of the C. elegans endpoints, at least 1 OPFR had similar toxicity to the PBDEs: triphenyl phosphate (TPHP) inhibited larval development at levels comparable to the 3 PBDEs; whereas TPHP and isopropylated phenol phosphate (IPP) affected C. elegans reproduction at levels similar to the PBDE commercial mixture, DE-71. The PBDEs reduced C. elegans feeding at lower concentrations than any OPFR. In addition, 9 of the 11 chemicals that inhibited C. elegans larval development also caused significant mitochondrial toxicity. These results suggest that some of the replacement aromatic OPFRs may have levels of toxicity comparable to PBDEs.

Developmental Effects of the ToxCast™ Phase I and Phase II Chemicals in Caenorhabditis elegans and Corresponding Responses in Zebrafish, Rats, and Rabbits.

BACKGROUND: Modern toxicology is shifting from an observational to a mechanistic science. As part of this shift, high-throughput toxicity assays are being developed using alternative, nonmammalian species to prioritize chemicals and develop prediction models of human toxicity. METHODS: The nematode Caenorhabditis elegans (C. elegans) was used to screen the U.S. Environmental Protection Agency's (EPA's) ToxCast™ Phase I and Phase II libraries, which contain 292 and 676 chemicals, respectively, for chemicals leading to decreased larval development and growth. Chemical toxicity was evaluated using three parameters: a biologically defined effect size threshold, half-maximal activity concentration (AC50), and lowest effective concentration (LEC). RESULTS: Across both the Phase I and Phase II libraries, 62% of the chemicals were classified as active ≤ 200 µM in the C. elegans assay. Chemical activities and potencies in C. elegans were compared with those from two zebrafish embryonic development toxicity studies and developmental toxicity data for rats and rabbits. Concordance of chemical activity was higher between C. elegans and one zebrafish assay across Phase I chemicals (79%) than with a second zebrafish assay (59%). Using C. elegans or zebrafish to predict rat or rabbit developmental toxicity resulted in balanced accuracies (the average value of the sensitivity and specificity for an assay) ranging from 45% to 53%, slightly lower than the concordance between rat and rabbit (58%). CONCLUSIONS: Here, we present an assay that quantitatively and reliably describes the effects of chemical toxicants on C. elegans growth and development. We found significant overlap in the activity of chemicals in the ToxCast™ libraries between C. elegans and zebrafish developmental screens. Incorporating C. elegans toxicological assays as part of a battery of in vitro and in vivo assays provides additional information for the development of models to predict a chemical's potential toxicity to humans. CITATION: Boyd WA, Smith MV, Co CA, Pirone JR, Rice JR, Shockley KR, Freedman JH. 2016. Developmental effects of the ToxCast™ Phase I and II chemicals in Caenorhabditis elegans and corresponding responses in zebrafish, rats, and rabbits. Environ Health Perspect 124:586-593; http://dx.doi.org/10.1289/ehp.1409645.

Use of alternative assays to identify and prioritize organophosphorus flame retardants for potential developmental and neurotoxicity.

Due to their toxicity and persistence in the environment, brominated flame retardants (BFRs) are being phased out of commercial use, leading to the increased use of alternative chemicals such as the organophosphorus flame retardants (OPFRs). There is, however, limited information on the potential health effects of OPFRs. Due to the structural similarity of the OPFRs to organophosphorus insecticides, there is concern regarding developmental toxicity and neurotoxicity. In response, we evaluated a set of 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]) in a battery of cell-based in vitro assays and alternative model organisms and compared the results to those obtained for two classical BFRs (3,3',5,5'-tetrabromobisphenol A [TBBPA] and 2,2'4,4'-brominated diphenyl ether [BDE-47]). The assays used evaluated the effects of chemicals on the differentiation of mouse embryonic stem cells, the proliferation and growth of human neural stem cells, rat neuronal growth and network activity, and development of nematode (Caenorhabditis elegans) and zebrafish (Danio rerio). All assays were performed in a concentration-response format, allowing for the determination of the point of departure (POD: the lowest concentration where a chemically-induced response exceeds background noise). The majority of OPFRs (8/9) were active in multiple assays in the range of 1-10 µM, most of which had comparable activity to the BFRs TBBPA and BDE-47. TCEP was negative in all assays. The results indicate that the replacement OPFRs, with the exception of TCEP, showed comparable activity to the two BFRs in the assays tested. Based on these results, more comprehensive studies are warranted to further characterize the potential hazard of some of these OPFR compounds.

Comparison of the toxicity of fluoridation compounds in the nematode Caenorhabditis elegans.

Fluorides are commonly added to drinking water in the United States to decrease the incidence of dental caries. Silicofluorides, such as sodium hexafluorosilicate (Na2 SiF6 ) and fluorosilicic acid (H2 SiF6 ), are mainly used for fluoridation, although fluoride salts such as sodium fluoride (NaF) are also used. Interestingly, only the toxicity of NaF has been examined and not that of the more often used silicofluorides. In the present study, the toxicities of NaF, Na2 SiF6 , and H2 SiF6 were compared. The toxicity of these fluorides on the growth, feeding, and reproduction in the alternative toxicological testing organism Caenorhabditis elegans was examined. Exposure to these compounds produced classic concentration-response toxicity profiles. Although the effects of the fluoride compounds varied among the 3 biological endpoints, no differences were found between the 3 compounds, relative to the fluoride ion concentration, in any of the assays. This suggests that silicofluorides have similar toxicity to NaF.

Amelioration of metal-induced toxicity in Caenorhabditis elegans: utility of chelating agents in the bioremediation of metals.

The presence of toxic amounts of transition metals in the environment may originate from a range of human activities and natural processes. One method for the removal of toxic levels of metals is through chelation by small molecules. However, chelation is not synonymous with detoxification and may not affect the bioavailability of the metal. To test the bioavailability of chelated metals in vivo, the effects of several metal/chelator combinations were tested in the environmentally relevant organism Caenorhabditis elegans. The effect of metal exposure on nematode growth was used to determine the toxicity of cadmium, copper, nickel, and zinc. The restoration of growth to levels observed in nonexposed nematodes was used to determine the protective effects of the polydentate chelators: acetohydroxamic acid (AHA), cyclam, cysteine, calcium EDTA, desferrioxamine B, 1,2-dimethyl,3-hydroxy,4-pyridinone, and histidine. Cadmium toxicity was removed only by EDTA; copper toxicity was removed by all of the chelators except AHA; nickel toxicity was removed by cyclam, EDTA, and histidine; and zinc toxicity was removed by only EDTA. These results demonstrate the utility of polydentate chelators in the remediation of metal-contaminated systems. They also demonstrate that although the application of a chelator to metal contaminants may be effective, binding alone cannot be used to predict the level of remediation. Remediation depends on a number of factors, including metal complex speciation in the environment.

A high-throughput method for assessing chemical toxicity using a Caenorhabditis elegans reproduction assay.

The National Research Council has outlined the need for non-mammalian toxicological models to test the potential health effects of a large number of chemicals while also reducing the use of traditional animal models. The nematode Caenorhabditis elegans is an attractive alternative model because of its well-characterized and evolutionarily conserved biology, low cost, and ability to be used in high-throughput screening. A high-throughput method is described for quantifying the reproductive capacity of C. elegans exposed to chemicals for 48 h from the last larval stage (L4) to adulthood using a COPAS Biosort. Initially, the effects of exposure conditions that could influence reproduction were defined. Concentrations of DMSO vehicle

Application of a mathematical model to describe the effects of chlorpyrifos on Caenorhabditis elegans development.

BACKGROUND: The nematode Caenorhabditis elegans is being assessed as an alternative model organism as part of an interagency effort to develop better means to test potentially toxic substances. As part of this effort, assays that use the COPAS Biosort flow sorting technology to record optical measurements (time of flight (TOF) and extinction (EXT)) of individual nematodes under various chemical exposure conditions are being developed. A mathematical model has been created that uses Biosort data to quantitatively and qualitatively describe C. elegans growth, and link changes in growth rates to biological events. Chlorpyrifos, an organophosphate pesticide known to cause developmental delays and malformations in mammals, was used as a model toxicant to test the applicability of the growth model for in vivo toxicological testing. METHODOLOGY/PRINCIPAL FINDINGS: L1 larval nematodes were exposed to a range of sub-lethal chlorpyrifos concentrations (0-75 microM) and measured every 12 h. In the absence of toxicant, C. elegans matured from L1s to gravid adults by 60 h. A mathematical model was used to estimate nematode size distributions at various times. Mathematical modeling of the distributions allowed the number of measured nematodes and log(EXT) and log(TOF) growth rates to be estimated. The model revealed three distinct growth phases. The points at which estimated growth rates changed (change points) were constant across the ten chlorpyrifos concentrations. Concentration response curves with respect to several model-estimated quantities (numbers of measured nematodes, mean log(TOF) and log(EXT), growth rates, and time to reach change points) showed a significant decrease in C. elegans growth with increasing chlorpyrifos concentration. CONCLUSIONS: Effects of chlorpyrifos on C. elegans growth and development were mathematically modeled. Statistical tests confirmed a significant concentration effect on several model endpoints. This confirmed that chlorpyrifos affects C. elegans development in a concentration dependent manner. The most noticeable effect on growth occurred during early larval stages: L2 and L3. This study supports the utility of the C. elegans growth assay and mathematical modeling in determining the effects of potentially toxic substances in an alternative model organism using high-throughput technologies.

A discrete time model for the analysis of medium-throughput C. elegans growth data.

BACKGROUND: As part of a program to predict the toxicity of environmental agents on human health using alternative methods, several in vivo high- and medium-throughput assays are being developed that use C. elegans as a model organism. C. elegans-based toxicological assays utilize the COPAS Biosort flow sorting system that can rapidly measure size, extinction (EXT) and time-of-flight (TOF), of individual nematodes. The use of this technology requires the development of mathematical and statistical tools to properly analyze the large volumes of biological data. METHODOLOGY/PRINCIPAL FINDINGS: Findings A Markov model was developed that predicts the growth of populations of C. elegans. The model was developed using observations from a 60 h growth study in which five cohorts of 300 nematodes each were aspirated and measured every 12 h. Frequency distributions of log(EXT) measurements that were made when loading C. elegans L1 larvae into 96 well plates (t = 0 h) were used by the model to predict the frequency distributions of the same set of nematodes when measured at 12 h intervals. The model prediction coincided well with the biological observations confirming the validity of the model. The model was also applied to log(TOF) measurements following an adaptation. The adaptation accounted for variability in TOF measurements associated with potential curling or shortening of the nematodes as they passed through the flow cell of the Biosort. By providing accurate estimates of frequencies of EXT or TOF measurements following varying growth periods, the model was able to estimate growth rates. Best model fits showed that C. elegans did not grow at a constant exponential rate. Growth was best described with three different rates. Microscopic observations indicated that the points where the growth rates changed corresponded to specific developmental events: the L1/L2 molt and the start of oogenesis in young adult C. elegans. CONCLUSIONS: Quantitative analysis of COPAS Biosort measurements of C. elegans growth has been hampered by the lack of a mathematical model. In addition, extraneous matter and the inability to assign specific measurements to specific nematodes made it difficult to estimate growth rates. The present model addresses these problems through a population-based Markov model.