A Quantitative Model of Systemic Toxicity Using ToxCast and ToxRefDB

Presented at the 9th Worl

A New Statistical Approach to Characterize Chemical-Elicited Behavioral Effects in High-Throughput Studies Using Zebrafish

<div><p>Zebrafish have become an important alternative model for characterizing chemical bioactivity, partly due to the efficiency at which systematic, high-dimensional data can be generated. However, these new data present analytical challenges associated with scale and diversity. We developed a novel, robust statistical approach to characterize chemical-elicited effects in behavioral data from high-throughput screening (HTS) of all 1,060 Toxicity Forecaster (ToxCast™) chemicals across 5 concentrations at 120 hours post-fertilization (hpf). Taking advantage of the immense scale of data for a global view, we show that this new approach reduces bias introduced by extreme values yet allows for diverse response patterns that confound the application of traditional statistics. We have also shown that, as a summary measure of response for local tests of chemical-associated behavioral effects, it achieves a significant reduction in coefficient of variation compared to many traditional statistical modeling methods. This effective increase in signal-to-noise ratio augments statistical power and is observed across experimental periods (light/dark conditions) that display varied distributional response patterns. Finally, we integrated results with data from concomitant developmental endpoint measurements to show that appropriate statistical handling of HTS behavioral data can add important biological context that informs mechanistic hypotheses.</p></div

Application of ToxCast and ToxRefDB to Develop a Quantitative Model of Systemic Toxicity

Application of ToxCast and ToxRefDB to Develop a Quantitative Model of Systemic Toxicit

Quantitative Model of Systemic Toxicity Using ToxCast and ToxRef DB

Presented at the Annual Society of Toxicolog

Behavioral response patterns.

<p>Differential entropy of each concentration was plotted across experimental time. Color key is shown in <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0169408#pone.0169408.g001" target="_blank">Fig 1</a>. Y axis: Differential entropy (Nats); X axis: Time (m). Red segments represent light condition from 3m to 9m. <b>A: Inactive:</b> TX000888 (Terbacil) was inactive at all concentrations. <b>B: Hypoactivity (L) and Hypoactivity (D):</b> TX001406 (Cyclanilide) shows significant hypoactivities at 64uM for both light and dark intervals. <b>C: Inactive (L) and Hypoactivity (D):</b> TX001412 (Fipronil) is inactive at light interval and shows significant hypoactivity at dark interval at 0.064uM, 0.64uM, 6.4uM, and 64uM. <b>D: Hyperactivity (L) and Inactive (D):</b> TX007214 (Dieldrin) shows significant hyperactivity at light interval but it is inactive at dark at 64uM. <b>E: Hypoactivity (L) and Inactive (D):</b> TX003357 (44’-Oxydianiline) shows significant hypoactivity at light interval and inactive pattern at dark interval at 0.064uM, 6.4uM, and 64uM. <b>F: Hyperactivity (L) and Hypoactivity (D):</b> TX006644 (Haloperidol) shows significant hyperactivity at light interval and significant hypoactivity at dark interval at both 6.4uM and 64uM. In addition, at 0.64uM, it shows significant hyperactivity at light interval. <b>G: Hyperactivity (L) and Hyperactivity (D):</b> TX005098 (4-Pentylaniline) shows significant hyperactivity at 64uM for both light and dark conditions. <b>H: Inactive (L) and Hyperactivity (D):</b> TX005080 (44’4”-Ethane-111-triyltriphenol) is inactive at light and shows significant hyperactivity at dark at 6.4uM and 64uM.</p

Framework for a Quantitative Systemic Toxicity Model

Framework for a Quantitative Systemic Toxicity Mode

Performance of Our Novel Statistical Modeling Method.

<p><b>A:</b> Example controls (having similar survival rates) illustrate the transformation. These two separate controls were plotted by different colors. Blue: TX000769 (Propoxycarbazone-sodium); Black: TX000900 (Methamidophos). Top: Movement index of each time point, and the line was drawn by connecting the mean movement indexes. Y axis: Movement index; X axis: Time. Bottom: Lines were drawn using our method, which connects the differential entropy of each time point. Y axis: Differential entropy (Nats); X axis: Time. <b>B:</b> All 1,060 control groups were plotted. Top: Each line represents a chemical. Line was drawn by connecting mean movement indexes. Y axis: Movement index; X axis: Time. Bottom: Each line represents a chemical. Line was drawn by connecting differential entropy across time. Y axis: Differential entropy (Nats); X axis: Time. <b>C:</b> Coefficient variation of each time point using all control groups and various statistical modeling methods. Y axis: Coefficient of Variation; X axis: Time.</p

Summary of results by the whole experimental system.

<p>This Venn diagram provides the summary of the total number of significant chemicals detected by each assay. It also provides statistics regarding the benefits of including all assays. Missing rate: the number of chemicals that would have been missed using a subset of these assays.</p

Statistical Workflow.

<p>Step 1: Visualize movement index; Step 2: Remove annotated dead fish for every concentration of a chemical; Step 3: Propose a statistical modeling method; Step 4: Check for artifacts, such as technical issues, global plate and position effect, and remove any bad plates; P1: Plate 1; P2: Plate 2; C1: Column 1; C12: Column 12; Step 5: Apply statistical modeling method to provide dose-response patterns for analysis; Step 6: Statistical analysis pipeline; Step 7: Assess reproducibility of our computational framework.</p

Relationship between morphological profiles and behavioral profiles.

<p>Disulfiram significantly affected 13 endpoints starting at 0.64 uM. Disulfiram also caused significant hypoactivity in both intervals with a lowest effect level of 0.64 uM. For morphological profiles, the panels represent (from top left) Aggregate Entropy, mortality, summation of any endpoint, then each of the specific endpoints (see ‘<a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0169408#sec002" target="_blank">Methods</a>‘). The X axes show concentration (0uM, 0.0064uM, 0.064uM, 0.64uM, 6.4uM, 64uM from left to right). The Y axes show Aggregate Entropy for the first panel, then incidence counts for all other panels. Red indicates statistical significance for each measure (p < 0.05).</p