Analysis of incidence data in developmental toxicity studies: Statistical tests to account for litter effects in fetal defect data.

BACKGROUND: When analyzing fetal defect incidence in laboratory animal studies, correlation in responses within litters (i.e., litter effects) can lead to increased false-positive rates if litter effects are not incorporated into the analysis. Studies of fetal defects require analysis methods that are robust across a broad range of defect types, including those with zero or near-zero incidence rates in control groups. METHODS: A simulation study compared power and false-positive rates for six approaches across a range of background defect rates and litter size distributions. Statistical methods evaluated included ignoring the litter effect as well as parametric and nonparametric approaches based on litter proportions, generalized linear mixed models (GLMMs), the Rao-Scott Cochran-Armitage (RSCA) trend test, and a modification to the RSCA (mRSCA) introduced here to improve estimation at low background rates. These methods were also applied to a common and a rare defect from two prenatal developmental toxicology studies conducted by the National Toxicology Program (NTP). RESULTS: At background defect rates of 1%, the mRSCA and parametric litter proportion methods provided gains in power over the nonparametric litter proportion method, the GLMM method, and the RSCA method. Simulations involving litter loss in high-dose groups showed loss of power for both litter proportion methods. CONCLUSIONS: The mRSCA test developed here compares favorably with other litter-based approaches and is robust across a range of background defect rates and litter size distributions, making it a practical choice for prenatal developmental toxicology studies involving both common and rare fetal defects.

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.

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.

Open source software implementation of an integrated testing strategy for skin sensitization potency based on a Bayesian network.

An open-source implementation of a previously published integrated testing strategy (ITS) for skin sensitization using a Bayesian network has been developed using R, a free and open-source statistical computing language. The ITS model provides probabilistic predictions of skin sensitization potency based on in silico and in vitro information as well as skin penetration characteristics from a published bioavailability model (Kasting et al., 2008). The structure of the Bayesian network was designed to be consistent with the adverse outcome pathway published by the OECD (Jaworska et al., 2011, 2013). In this paper, the previously published data set (Jaworska et al., 2013) is improved by two data corrections and a modified application of the Kasting model. The new data set implemented in the original commercial software package and the new R version produced consistent results. The data and a fully documented version of the code are publicly available (http://ntp.niehs.nih.gov/go/its).

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.

Comparative toxicogenomic responses of mercuric and methyl-mercury.

BACKGROUND: Mercury is a ubiquitous environmental toxicant that exists in multiple chemical forms. A paucity of information exists regarding the differences or similarities by which different mercurials act at the molecular level. RESULTS: Transcriptomes of mixed-stage C. elegans following equitoxic sub-, low- and high-toxicity exposures to inorganic mercuric chloride (HgCl2) and organic methylmercury chloride (MeHgCl) were analyzed. In C. elegans, the mercurials had highly different effects on transcription, with MeHgCl affecting the expression of significantly more genes than HgCl2. Bioinformatics analysis indicated that inorganic and organic mercurials affected different biological processes. RNAi identified 18 genes that were important in C. elegans response to mercurial exposure, although only two of these genes responded to both mercurials. To determine if the responses observed in C. elegans were evolutionarily conserved, the two mercurials were investigated in human neuroblastoma (SK-N-SH), hepatocellular carcinoma (HepG2) and embryonic kidney (HEK293) cells. The human homologs of the affected C. elegans genes were then used to test the effects on gene expression and cell viability after using siRNA during HgCl2 and MeHgCl exposure. As was observed with C. elegans, exposure to the HgCl2 and MeHgCl had different effects on gene expression, and different genes were important in the cellular response to the two mercurials. CONCLUSIONS: These results suggest that, contrary to previous reports, inorganic and organic mercurials have different mechanisms of toxicity. The two mercurials induced disparate effects on gene expression, and different genes were important in protecting the organism from mercurial toxicity.

Caenorhabditis elegans as a model in developmental toxicology.

A number of practical advantages have made the nematode Caenorhabditis elegans a useful model for genetic and developmental biological research. These same advantages, along with conservation of disease and stress response pathways, availability of mutant and transgenic strains, and wealth of biological information, have led to the increased use of C. elegans in toxicological studies. Although the potential to study the mechanisms of developmental toxicology in C. elegans is promising, embryonic and larval growth tests to identify compounds that affect the nematode have remained the primary use of C. elegans in developmental toxicology. Here, we describe a C. elegans larval growth and development assay for medium- and high-throughput screening using the COPAS Biosort flow cytometer and provide descriptions of the data and subsequent analysis.

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.

Medium- and high-throughput screening of neurotoxicants using C. elegans.

The U.S. National Toxicology Program, the U.S. Environmental Protection Agency, and other national and international agencies are committing significant resources towards the development of alternative species to be used as replacements for mammalian models in toxicological studies. Caenorhabditis elegans is a well-characterized soil nematode that is becoming a useful model in the assessment of neurotoxicants. To determine the effects of potential neurotoxicants on C. elegans, four medium-throughput (feeding, growth, reproduction and locomotion) and two high-throughput (growth and reproduction) assays have been developed. Three of these assays use the COPAS Biosort, a flow cytometer capable of rapidly measuring thousands of nematodes in minutes. Medium-throughput feeding, growth, and reproduction assays were used to assess the toxicity of eight suspected neurotoxicants. For several of the neurotoxicants examined, significant effects were observed at similar concentrations between assays. High-throughput reproduction and growth assays were used to estimate the toxicity of thousands of chemicals in two libraries. These assays will prove useful in evaluating the role of alternative toxicological models in tiered toxicity testing of thousands of chemicals.

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.

Characterization of the proneural gene regulatory network during mouse telencephalon development.

BACKGROUND: The proneural proteins Mash1 and Ngn2 are key cell autonomous regulators of neurogenesis in the mammalian central nervous system, yet little is known about the molecular pathways regulated by these transcription factors. RESULTS: Here we identify the downstream effectors of proneural genes in the telencephalon using a genomic approach to analyze the transcriptome of mice that are either lacking or overexpressing proneural genes. Novel targets of Ngn2 and/or Mash1 were identified, such as members of the Notch and Wnt pathways, and proteins involved in adhesion and signal transduction. Next, we searched the non-coding sequence surrounding the predicted proneural downstream effector genes for evolutionarily conserved transcription factor binding sites associated with newly defined consensus binding sites for Ngn2 and Mash1. This allowed us to identify potential novel co-factors and co-regulators for proneural proteins, including Creb, Tcf/Lef, Pou-domain containing transcription factors, Sox9, and Mef2a. Finally, a gene regulatory network was delineated using a novel Bayesian-based algorithm that can incorporate information from diverse datasets. CONCLUSION: Together, these data shed light on the molecular pathways regulated by proneural genes and demonstrate that the integration of experimentation with bioinformatics can guide both hypothesis testing and hypothesis generation.