Transcriptomic analysis of AHR wildtype and Knock-out rat livers supports TCDD's role in AHR/ARNT-mediated circadian disruption and hepatotoxicity.

Single, high doses of TCDD in rats are known to cause wasting, a progressive loss of 30 to 50% body weight and death within several weeks. To identify pathway perturbations at or near doses causing wasting, we examined differentially gene expression (DGE) and pathway enrichment in centrilobular (CL) and periportal (PP) regions of female rat livers following 6 dose levels of TCDD - 0, 3, 22, 100, 300, and 1000 ng/kg/day, 5 days/week for 4 weeks. At the higher doses, rats lost weight, had increased liver/body weight ratios and nearly complete cessation of liver cell proliferation, signs consistent with wasting. DGE curves were left shifted for the CL versus the PP regions. Canonical Phase I and Phase II genes were maximally increased at lower doses and remained elevated at all doses. At lower doses, ≤ 22 ng/kg/day in the CL and ≤ 100 ng/kg/day, upregulated genes showed transcription factor (TF) enrichment for AHR and ARNT. At the mid- and high-dose doses, there was a large number of downregulated genes and pathway enrichment for DEGs which showed downregulation of many cellular metabolism processes including those for steroids, fatty acid metabolism, pyruvate metabolism and citric acid cycle. There was significant TF enrichment of the hi-dose downregulated genes for RXR, ESR1, LXR, PPARalpha. At the highest dose, there was also pathway enrichment with upregulated genes for extracellular matrix organization, collagen formation, hemostasis and innate immune system. TCDD demonstrates most of its effects through binding the aryl hydrocarbon receptor (AHR) while the downregulation of metabolism genes at higher TCDD doses is known to be independent of AHR binding to DREs. Based on our results with DEG, we provide a hypothesis for wasting in which high doses of TCDD shift circadian processes away from the resting state, leading to greatly reduced synthesis of steroids and complex lipids needed for cell growth, and producing gene expression signals consistent with an epithelial-to-mesenchymal transition in hepatocytes.

Application of Cell Painting for chemical hazard evaluation in support of screening-level chemical assessments.

'Cell Painting' is an imaging-based high-throughput phenotypic profiling (HTPP) method in which cultured cells are fluorescently labeled to visualize subcellular structures (i.e., nucleus, nucleoli, endoplasmic reticulum, cytoskeleton, Golgi apparatus / plasma membrane and mitochondria) and to quantify morphological changes in response to chemicals or other perturbagens. HTPP is a high-throughput and cost-effective bioactivity screening method that detects effects associated with many different molecular mechanisms in an untargeted manner, enabling rapid in vitro hazard assessment for thousands of chemicals. Here, 1201 chemicals from the ToxCast library were screened in concentration-response up to ∼100 µM in human U-2 OS cells using HTPP. A phenotype altering concentration (PAC) was estimated for chemicals active in the tested range. PACs tended to be higher than lower bound potency values estimated from a broad collection of targeted high-throughput assays, but lower than the threshold for cytotoxicity. In vitro to in vivo extrapolation (IVIVE) was used to estimate administered equivalent doses (AEDs) based on PACs for comparison to human exposure predictions. AEDs for 18/412 chemicals overlapped with predicted human exposures. Phenotypic profile information was also leveraged to identify putative mechanisms of action and group chemicals. Of 58 known nuclear receptor modulators, only glucocorticoids and retinoids produced characteristic profiles; and both receptor types are expressed in U-2 OS cells. Thirteen chemicals with profile similarity to glucocorticoids were tested in a secondary screen and one chemical, pyrene, was confirmed by an orthogonal gene expression assay as a novel putative GR modulating chemical. Most active chemicals demonstrated profiles not associated with a known mechanism-of-action. However, many structurally related chemicals produced similar profiles, with exceptions such as diniconazole, whose profile differed from other active conazoles. Overall, the present study demonstrates how HTPP can be applied in screening-level chemical assessments through a series of examples and brief case studies.

Combining phenotypic profiling and targeted RNA-Seq reveals linkages between transcriptional perturbations and chemical effects on cell morphology: Retinoic acid as an example.

The United States Environmental Protection Agency has proposed a tiered testing strategy for chemical hazard evaluation based on new approach methods (NAMs). The first tier includes in vitro profiling assays applicable to many (human) cell types, such as high-throughput transcriptomics (HTTr) and high-throughput phenotypic profiling (HTPP). The goals of this study were to: (1) harmonize the seeding density of U-2 OS human osteosarcoma cells for use in both assays; (2) compare HTTr- versus HTPP-derived potency estimates for 11 mechanistically diverse chemicals; (3) identify candidate reference chemicals for monitoring assay performance in future screens; and (4) characterize the transcriptional and phenotypic changes in detail for all-trans retinoic acid (ATRA) as a model compound known for its adverse effects on osteoblast differentiation. The results of this evaluation showed that (1) HTPP conducted at low (400 cells/well) and high (3000 cells/well) seeding densities yielded comparable potency estimates and similar phenotypic profiles for the tested chemicals; (2) HTPP and HTTr resulted in comparable potency estimates for changes in cellular morphology and gene expression, respectively; (3) three test chemicals (etoposide, ATRA, dexamethasone) produced concentration-dependent effects on cellular morphology and gene expression that were consistent with known modes-of-action, demonstrating their suitability for use as reference chemicals for monitoring assay performance; and (4) ATRA produced phenotypic changes that were highly similar to other retinoic acid receptor activators (AM580, arotinoid acid) and some retinoid X receptor activators (bexarotene, methoprene acid). This phenotype was observed concurrently with autoregulation of the RARB gene. Both effects were prevented by pre-treating U-2 OS cells with pharmacological antagonists of their respective receptors. Thus, the observed phenotype could be considered characteristic of retinoic acid pathway activation in U-2 OS cells. These findings lay the groundwork for combinatorial screening of chemicals using HTTr and HTPP to generate complementary information for the first tier of a NAM-based chemical hazard evaluation strategy.

Progress towards an OECD reporting framework for transcriptomics and metabolomics in regulatory toxicology.

Omics methodologies are widely used in toxicological research to understand modes and mechanisms of toxicity. Increasingly, these methodologies are being applied to questions of regulatory interest such as molecular point-of-departure derivation and chemical grouping/read-across. Despite its value, widespread regulatory acceptance of omics data has not yet occurred. Barriers to the routine application of omics data in regulatory decision making have been: 1) lack of transparency for data processing methods used to convert raw data into an interpretable list of observations; and 2) lack of standardization in reporting to ensure that omics data, associated metadata and the methodologies used to generate results are available for review by stakeholders, including regulators. Thus, in 2017, the Organisation for Economic Co-operation and Development (OECD) Extended Advisory Group on Molecular Screening and Toxicogenomics (EAGMST) launched a project to develop guidance for the reporting of omics data aimed at fostering further regulatory use. Here, we report on the ongoing development of the first formal reporting framework describing the processing and analysis of both transcriptomic and metabolomic data for regulatory toxicology. We introduce the modular structure, content, harmonization and strategy for trialling this reporting framework prior to its publication by the OECD.

Bioactivity screening of environmental chemicals using imaging-based high-throughput phenotypic profiling.

The present study adapted an existing high content imaging-based high-throughput phenotypic profiling (HTPP) assay known as "Cell Painting" for bioactivity screening of environmental chemicals. This assay uses a combination of fluorescent probes to label a variety of organelles and measures a large number of phenotypic features at the single cell level in order to detect chemical-induced changes in cell morphology. First, a small set of candidate phenotypic reference chemicals (n = 14) known to produce changes in the cellular morphology of U-2 OS cells were identified and screened at multiple time points in concentration-response format. Many of these chemicals produced distinct cellular phenotypes that were qualitatively similar to those previously described in the literature. A novel workflow for phenotypic feature extraction, concentration-response modeling and determination of in vitro thresholds for chemical bioactivity was developed. Subsequently, a set of 462 chemicals from the ToxCast library were screened in concentration-response mode. Bioactivity thresholds were calculated and converted to administered equivalent doses (AEDs) using reverse dosimetry. AEDs were then compared to effect values from mammalian toxicity studies. In many instances (68%), the HTPP-derived AEDs were either more conservative than or comparable to the in vivo effect values. Overall, we conclude that the HTPP assay can be used as an efficient, cost-effective and reproducible screening method for characterizing the biological activity and potency of environmental chemicals for potential use in in vitro-based safety assessments.

Testing for developmental neurotoxicity using a battery of in vitro assays for key cellular events in neurodevelopment.

Medium- to high-throughput in vitro assays that recapitulate the critical processes of nervous system development have been proposed as a means to facilitate rapid testing and identification of chemicals which may affect brain development. In vivo neurodevelopment is a complex progression of distinct cellular processes. Therefore, batteries of in vitro assays that model and quantify effects on a variety of neurodevelopmental processes have the potential to identify chemicals which may affect brain development at different developmental stages. In the present study, the results of concentration-response screening of 67 reference chemicals in a battery of high content imaging and microplate reader-based assays that evaluate neural progenitor cell proliferation, neural proginitor cell apoptosis, neurite initiation/outgrowth, neurite maturation and synaptogenesis are summarized and compared. The assay battery had a high degree of combined sensitivity (87%) for categorizing chemicals known to affect neurodevelopment as active and a moderate degree of combined specificity (71%) for categorizing chemicals not associated with affects on neurodevelopment as inactive. The combined sensitivity of the assay battery was higher compared to any individual assay while the combined specificity of the assay battery was lower compared to any individual assay. When selectivity of effects for a neurodevelopmental endpoint as compared to general cytotoxicity was taken into account, the combined sensitivity of the assay battery decreased (68%) while the combined specificity increased (93%). The identity and potency of chemicals identified as active varied across the assay battery, underscoring the need for use of a combination of diverse in vitro models to comprehensively screen chemicals and identify those which potentially affect neurodevelopment. Overall, these data indicate that a battery of assays which address many different processes in nervous system development may be used to identify potential developmental neurotoxicants and to distinguish specific from generalized cytotoxic effects with a high degree of success.

Lineage-dependent effects of aryl hydrocarbon receptor agonists contribute to liver tumorigenesis.

UNLABELLED: Rodent cancer bioassays indicate that the aryl hydrocarbon receptor (AHR) agonist, 2,3,7,8-tetracholorodibenzo-p-dioxin (TCDD), causes increases in both hepatocytic and cholangiocytic tumors. Effects of AHR activation have been evaluated on rodent hepatic stem cells (rHpSCs) versus their descendants, hepatoblasts (rHBs), two lineage stages of multipotent, hepatic precursors with overlapping but also distinct phenotypic traits. This was made possible by defining the first successful culture conditions for ex vivo maintenance of rHpScs consisting of a substratum of hyaluronans and Kubota's medium (KM), a serum-free medium designed for endodermal stem/progenitor cells. Supplementation of KM with leukemia inhibitory factor elicited lineage restriction to rHBs. Cultures were treated with various AHR agonists including TCDD, 6-formylindolo-[3,2-b]carbazole (FICZ), and 3-3'-diindolylmethane (DIM) and then analyzed with a combination of immunocytochemistry, gene expression, and high-content image analysis. The AHR agonists increased proliferation of rHpSCs at concentrations producing a persistent AHR activation as indicated by induction of Cyp1a1. By contrast, treatment with TCDD resulted in a rapid loss of viability of rHBs, even though the culture conditions, in the absence of the agonists, were permissive for survival and expansion of rHBs. The effects were not observed with FICZ and at lower concentrations of DIM. CONCLUSION: Our findings are consistent with a lineage-dependent mode of action for AHR agonists in rodent liver tumorigenesis through selective expansion of rHpSCs in combination with a toxicity-induced loss of viability of rHBs. These lineage-dependent effects correlate with increased frequency of liver tumors.

Aryl hydrocarbon receptor knockout rats are insensitive to the pathological effects of repeated oral exposure to 2,3,7,8-tetrachlorodibenzo-p-dioxin.

Sustained activation of the aryl hydrocarbon receptor (AHR) is believed to be the initial key event in AHR receptor-mediated tumorigenesis in the rat liver. The role of AHR in mediating pathological changes in the liver prior to tumor formation was investigated in a 4-week, repeated-dose study using adult female wild-type (WT) and AHR knockout (AHR-KO) rats treated with 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). Beginning at 8 weeks of age, AHR-KO and WT rats were dosed by oral gavage with varying concentrations of TCDD (0, 3, 22, 100, 300 and 1000 ng kg(-1) day(-1) ). Lung, liver and thymus histopathology, hematology, serum chemistry and the distribution of TCDD in liver and adipose tissue were examined. Treatment-related increases in the severity of liver and thymus pathology were observed in WT, but not AHR-KO rats. In the liver, these included hepatocellular hypertrophy, bile duct hyperplasia, multinucleated hepatocytes and inflammatory cell foci. A loss of cellularity in the thymic cortex and thymic atrophy was observed. Treatment-related changes in serum chemistry parameters were also observed in WT, but not AHR-KO rats. Finally, dose-dependent accumulation of TCDD was observed primarily in the liver of WT rats and primarily in the adipose tissue of AHR-KO rats. The results suggest that AHR activation is the initial key event underlying the progression of histological effects leading to liver tumorigenesis following TCDD treatment. Copyright © 2015 John Wiley & Sons, Ltd.

Exploring the effects of experimental parameters and data modeling approaches on in vitro transcriptomic point-of-departure estimates.

Multiple new approach methods (NAMs) are being developed to rapidly screen large numbers of chemicals to aid in hazard evaluation and risk assessments. High-throughput transcriptomics (HTTr) in human cell lines has been proposed as a first-tier screening approach for determining the types of bioactivity a chemical can cause (activation of specific targets vs. generalized cell stress) and for calculating transcriptional points of departure (tPODs) based on changes in gene expression. In the present study, we examine a range of computational methods to calculate tPODs from HTTr data, using six data sets in which MCF7 cells cultured in two different media formulations were treated with a panel of 44 chemicals for 3 different exposure durations (6, 12, 24 hr). The tPOD calculation methods use data at the level of individual genes and gene set signatures, and compare data processed using the ToxCast Pipeline 2 (tcplfit2), BMDExpress and PLIER (Pathway Level Information ExtractoR). Methods were evaluated by comparing to in vitro PODs from a validated set of high-throughput screening (HTS) assays for a set of estrogenic compounds. Key findings include: (1) for a given chemical and set of experimental conditions, tPODs calculated by different methods can vary by several orders of magnitude; (2) tPODs are at least as sensitive to computational methods as to experimental conditions; (3) in comparison to an external reference set of PODs, some methods give generally higher values, principally PLIER and BMDExpress; and (4) the tPODs from HTTr in this one cell type are mostly higher than the overall PODs from a broad battery of targeted in vitro ToxCast assays, reflecting the need to test chemicals in multiple cell types and readout technologies for in vitro hazard screening.

A Comparison of In Vitro Points of Departure with Human Blood Levels for Per- and Polyfluoroalkyl Substances (PFAS).

Per- and polyfluoroalkyl substances (PFAS) are widely used, and their fluorinated state contributes to unique uses and stability but also long half-lives in the environment and humans. PFAS have been shown to be toxic, leading to immunosuppression, cancer, and other adverse health outcomes. Only a small fraction of the PFAS in commerce have been evaluated for toxicity using in vivo tests, which leads to a need to prioritize which compounds to examine further. Here, we demonstrate a prioritization approach that combines human biomonitoring data (blood concentrations) with bioactivity data (concentrations at which bioactivity is observed in vitro) for 31 PFAS. The in vitro data are taken from a battery of cell-based assays, mostly run on human cells. The result is a Bioactive Concentration to Blood Concentration Ratio (BCBCR), similar to a margin of exposure (MoE). Chemicals with low BCBCR values could then be prioritized for further risk assessment. Using this method, two of the PFAS, PFOA (Perfluorooctanoic Acid) and PFOS (Perfluorooctane Sulfonic Acid), have BCBCR values < 1 for some populations. An additional 9 PFAS have BCBCR values < 100 for some populations. This study shows a promising approach to screening level risk assessments of compounds such as PFAS that are long-lived in humans and other species.

Knockout of the aryl hydrocarbon receptor results in distinct hepatic and renal phenotypes in rats and mice.

The aryl hydrocarbon receptor (AHR) is a ligand-activated transcription factor which plays a role in the development of multiple tissues and is activated by a large number of ligands, including 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). In order to examine the roles of the AHR in both normal biological development and response to environmental chemicals, an AHR knockout (AHR-KO) rat model was created and compared with an existing AHR-KO mouse. AHR-KO rats harboring either 2-bp or 29-bp deletion mutation in exon 2 of the AHR were created on the Sprague-Dawley genetic background using zinc-finger nuclease (ZFN) technology. Rats harboring either mutation type lacked expression of AHR protein in the liver. AHR-KO rats were also insensitive to thymic involution, increased hepatic weight and the induction of AHR-responsive genes (Cyp1a1, Cyp1a2, Cyp1b1, Ahrr) following acute exposure to 25 µg/kg TCDD. AHR-KO rats had lower basal expression of transcripts for these genes and also accumulated ~30-45-fold less TCDD in the liver at 7 days post-exposure. In untreated animals, AHR-KO mice, but not AHR-KO rats, had alterations in serum analytes indicative of compromised hepatic function, patent ductus venosus of the liver and persistent hyaloid arteries in the eye. AHR-KO rats, but not AHR-KO mice, displayed pathological alterations to the urinary tract: bilateral renal dilation (hydronephrosis), secondary medullary tubular and uroepithelial degenerative changes and bilateral ureter dilation (hydroureter). The present data indicate that the AHR may play significantly different roles in tissue development and homeostasis and toxicity across rodent species.

Neurotrophic effects of leukemia inhibitory factor on neural cells derived from human embryonic stem cells.

Various growth factor cocktails have been used to proliferate and then differentiate human neural progenitor (NP) cells derived from embryonic stem cells (ESC) for in vitro and in vivo studies. However, the cytokine leukemia inhibitory factor (LIF) has been largely overlooked. Here, we demonstrate that LIF significantly enhanced in vitro survival and promoted differentiation of human ESC-derived NP cells. In NP cells, as well as NP-derived neurons, LIF reduced caspase-mediated apoptosis and reduced both spontaneous and H2O2-induced reactive oxygen species in culture. In vitro, NP cell proliferation and the yield of differentiated neurons were significantly higher in the presence of LIF. In NP cells, LIF enhanced cMyc phosphorylation, commonly associated with self-renewal/proliferation. Also, in differentiating NP cells LIF activated the phosphoinositide 3-kinase and signal transducer and activator of transcription 3 pathways, associated with cell survival and reduced apoptosis. When differentiated in LIF+ media, neurite outgrowth and ERK1/2 phosphorylation were potentiated together with increased expression of gp130, a component of the LIF receptor complex. NP cells, pretreated in vitro with LIF, were effective in reducing infarct volume in a model of focal ischemic stroke but LIF did not lead to significantly improved initial NP cell survival over nontreated NP cells. Our results show that LIF signaling significantly promotes human NP cell proliferation, survival, and differentiation in vitro. Activated LIF signaling should be considered in cell culture expansion systems for future human NP cell-based therapeutic transplant studies.

Integrating Data From In Vitro New Approach Methodologies for Developmental Neurotoxicity.

In vivo developmental neurotoxicity (DNT) testing is resource intensive and lacks information on cellular processes affected by chemicals. To address this, DNT new approach methodologies (NAMs) are being evaluated, including: the microelectrode array neuronal network formation assay; and high-content imaging to evaluate proliferation, apoptosis, neurite outgrowth, and synaptogenesis. This work addresses 3 hypotheses: (1) a broad screening battery provides a sensitive marker of DNT bioactivity; (2) selective bioactivity (occurring at noncytotoxic concentrations) may indicate functional processes disrupted; and, (3) a subset of endpoints may optimally classify chemicals with in vivo evidence for DNT. The dataset was comprised of 92 chemicals screened in all 57 assay endpoints sourced from publicly available data, including a set of DNT NAM evaluation chemicals with putative positives (53) and negatives (13). The DNT NAM battery provides a sensitive marker of DNT bioactivity, particularly in cytotoxicity and network connectivity parameters. Hierarchical clustering suggested potency (including cytotoxicity) was important for classifying positive chemicals with high sensitivity (93%) but failed to distinguish patterns of disrupted functional processes. In contrast, clustering of selective values revealed informative patterns of differential activity but demonstrated lower sensitivity (74%). The false negatives were associated with several limitations, such as the maximal concentration tested or gaps in the biology captured by the current battery. This work demonstrates that this multi-dimensional assay suite provides a sensitive biomarker for DNT bioactivity, with selective activity providing possible insight into specific functional processes affected by chemical exposure and a basis for further research.

An expert-driven literature review of "negative" chemicals for developmental neurotoxicity (DNT) in vitro assay evaluation.

To date, approximately 200 chemicals have been tested in US Environmental Protection Agency (EPA) or Organization for Economic Co-operation and Development (OECD) developmental neurotoxicity (DNT) guideline studies, leaving thousands of chemicals without traditional animal information on DNT hazard potential. To address this data gap, a battery of in vitro DNT new approach methodologies (NAMs) has been proposed. Evaluation of the performance of this battery will increase the confidence in its use to determine DNT chemical hazards. One approach to evaluate DNT NAM performance is to use a set of chemicals to evaluate sensitivity and specificity. Since a list of chemicals with potential evidence of in vivo DNT has been established, this study aims to develop a curated list of "negative" chemicals for inclusion in a "DNT NAM evaluation set". A workflow, including a literature search followed by an expert-driven literature review, was used to systematically screen 39 chemicals for lack of DNT effect. Expert panel members evaluated the scientific robustness of relevant studies to inform chemical categorizations. Following review, the panel discussed each chemical and made categorical determinations of "Favorable", "Not Favorable", or "Indeterminate" reflecting acceptance, lack of suitability, or uncertainty given specific limitations and considerations, respectively. The panel determined that 10, 22, and 7 chemicals met the criteria for "Favorable", "Not Favorable", and "Indeterminate", for use as negatives in a DNT NAM evaluation set. Ultimately, this approach not only supports DNT NAM performance evaluation but also highlights challenges in identifying large numbers of negative DNT chemicals.

Benchmark Dose Modeling Approaches for Volatile Organic Chemicals Using a Novel Air-Liquid Interface In Vitro Exposure System.

Inhalation is the most relevant route of volatile organic chemical (VOC) exposure; however, due to unique challenges posed by their chemical properties and poor solubility in aqueous solutions, in vitro chemical safety testing is predominantly performed using direct application dosing/submerged exposures. To address the difficulties in screening toxic effects of VOCs, our cell culture exposure system permits cells to be exposed to multiple concentrations at air-liquid interface (ALI) in a 24-well format. ALI exposure methods permit direct chemical-to-cell interaction with the test article at physiological conditions. In the present study, BEAS-2B and primary normal human bronchial epithelial cells (pHBEC) are used to assess gene expression, cytotoxicity, and cell viability responses to a variety of volatile chemicals including acrolein, formaldehyde, 1,3-butadiene, acetaldehyde, 1-bromopropane, carbon tetrachloride, dichloromethane, and trichloroethylene. BEAS-2B cells were exposed to all the test agents, whereas pHBECs were only exposed to the latter 4 listed above. The VOC concentrations tested elicited only slight cell viability changes in both cell types. Gene expression changes were analyzed using benchmark dose (BMD) modeling. The BMD for the most sensitive gene set was within one order of magnitude of the threshold-limit value reported by the American Conference of Governmental Industrial Hygienists, and the most sensitive gene sets impacted by exposure correlate to known adverse health effects recorded in epidemiologic and in vivo exposure studies. Overall, our study outlines a novel in vitro approach for evaluating molecular-based points-of-departure in human airway epithelial cell exposure to volatile chemicals.

Comparative sensitivity of human and rat neural cultures to chemical-induced inhibition of neurite outgrowth.

There is a need for rapid, efficient and cost-effective alternatives to traditional in vivo developmental neurotoxicity testing. In vitro cell culture models can recapitulate many of the key cellular processes of nervous system development, including neurite outgrowth, and may be used as screening tools to identify potential developmental neurotoxicants. The present study compared primary rat cortical cultures and human embryonic stem cell-derived neural cultures in terms of: 1) reproducibility of high content image analysis based neurite outgrowth measurements, 2) dynamic range of neurite outgrowth measurements and 3) sensitivity to chemicals which have been shown to inhibit neurite outgrowth. There was a large increase in neurite outgrowth between 2 and 24h in both rat and human cultures. Image analysis data collected across multiple cultures demonstrated that neurite outgrowth measurements in rat cortical cultures were more reproducible and had higher dynamic range as compared to human neural cultures. Human neural cultures were more sensitive than rat cortical cultures to chemicals previously shown to inhibit neurite outgrowth. Parallel analysis of morphological (neurite count, neurite length) and cytotoxicity (neurons per field) measurements were used to detect selective effects on neurite outgrowth. All chemicals which inhibited neurite outgrowth in rat cortical cultures did so at concentrations which did not concurrently affect the number of neurons per field, indicating selective effects on neurite outgrowth. In contrast, more than half the chemicals which inhibited neurite outgrowth in human neural cultures did so at concentrations which concurrently decreased the number of neurons per field, indicating that effects on neurite outgrowth were secondary to cytotoxicity. Overall, these data demonstrate that the culture models performed differently in terms of reproducibility, dynamic range and sensitivity to neurite outgrowth inhibitors. While human neural cultures were more sensitive to neurite outgrowth inhibitors, they also had a lower dynamic range for detecting chemical-induced neurite outgrowth inhibition and greater variability from culture-to-culture as compared to rat primary cortical cultures.

Optimization of Human Neural Progenitor Cells for an Imaging-Based High-Throughput Phenotypic Profiling Assay for Developmental Neurotoxicity Screening.

Studies in in vivo rodent models have been the accepted approach by regulatory agencies to evaluate potential developmental neurotoxicity (DNT) of chemicals for decades. These studies, however, are inefficient and cannot meet the demand for the thousands of chemicals that need to be assessed for DNT hazard. As such, several in vitro new approach methods (NAMs) have been developed to circumvent limitations of these traditional studies. The DNT NAMs, some of which utilize human-derived cell models, are intended to be employed in a testing battery approach, each focused on a specific neurodevelopmental process. The need for multiple assays, however, to evaluate each process can prolong testing and prioritization of chemicals for more in depth assessments. Therefore, a multi-endpoint higher-throughput approach to assess DNT hazard potential would be of value. Accordingly, we have adapted a high-throughput phenotypic profiling (HTPP) approach for use with human-derived neural progenitor (hNP1) cells. HTPP is a fluorescence-based assay that quantitatively measures alterations in cellular morphology. This approach, however, required optimization of several laboratory procedures prior to chemical screening. First, we had to determine an appropriate cell plating density in 384-well plates. We then had to identify the minimum laminin concentration required for optimal cell growth and attachment. And finally, we had to evaluate whether addition of antibiotics to the culture medium would alter cellular morphology. We selected 6,000 cells/well as an appropriate plating density, 20 µg/ml laminin for optimal cell growth and attachment, and antibiotic addition in the culture medium. After optimizing hNP1 cell culture conditions for HTPP, it was then necessary to select appropriate in-plate assay controls from a reference chemical set. These reference chemicals were previously demonstrated to elicit unique phenotypic profiles in various other cell types. Aphidicolin, bafilomycin A1, berberine chloride, and cucurbitacin I induced robust phenotypic profiles as compared to dimethyl sulfoxide vehicle control in the hNP1 cells, and thus can be employed as in-plate assay controls for subsequent chemical screens. We have optimized HTPP for hNP1 cells, and consequently this approach can now be assessed as a potential NAM for DNT hazard evaluation and results compared to previously developed DNT assays.

Comparison of Approaches for Determining Bioactivity Hits from High-Dimensional Profiling Data.

Phenotypic profiling assays are untargeted screening assays that measure a large number (hundreds to thousands) of cellular features in response to a stimulus and often yield diverse and unanticipated profiles of phenotypic effects, leading to challenges in distinguishing active from inactive treatments. Here, we compare a variety of different strategies for hit identification in imaging-based phenotypic profiling assays using a previously published Cell Painting data set. Hit identification strategies based on multiconcentration analysis involve curve fitting at several levels of data aggregation (e.g., individual feature level, aggregation of similarly derived features into categories, and global modeling of all features) and on computed metrics (e.g., Euclidean and Mahalanobis distance metrics and eigenfeatures). Hit identification strategies based on single-concentration analysis included measurement of signal strength (e.g., total effect magnitude) and correlation of profiles among biological replicates. Modeling parameters for each approach were optimized to retain the ability to detect a reference chemical with subtle phenotypic effects while limiting the false-positive rate to 10%. The percentage of test chemicals identified as hits was highest for feature-level and category-based approaches, followed by global fitting, whereas signal strength and profile correlation approaches detected the fewest number of active hits at the fixed false-positive rate. Approaches involving fitting of distance metrics had the lowest likelihood for identifying high-potency false-positive hits that may be associated with assay noise. Most of the methods achieved a 100% hit rate for the reference chemical and high concordance for 82% of test chemicals, indicating that hit calls are robust across different analysis approaches.

High-Throughput Transcriptomics Platform for Screening Environmental Chemicals.

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

Splice variant specific increase in Ca2+/calmodulin-dependent protein kinase 1-gamma mRNA expression in response to acute pyrethroid exposure.

In mammals, pyrethroids are neurotoxicants that interfere with ion channel function in excitable neuronal membranes. Previous work demonstrated increases in the expression of Ca(2+)/calmodulin-dependent protein kinase 1-gamma (Camk1g) mRNA following acute deltamethrin and permethrin exposure. In the rat, this gene is expressed as two distinct splice variants, Camk1g1 and Camk1g2. The present study tests the hypothesis that changes in Camk1g mRNA expression in the rat following acute pyrethroid exposure are due to a specific increase in the Camk1g1 splice variant and not the Camk1g2 splice variant. Long-Evans rats were acutely exposed to permethrin, deltamethrin, or corn oil vehicle. Frontal cortex was collected at 6 h postdosing. In addition, rats were exposed to permethrin (100 mg/kg) or deltamethrin (3 mg/kg), and frontal cortex was collected at 1, 3, 6, 9, 12, or 24 h along with time-matched vehicle controls. Expression of Camk1g1 and Camk1g2 mRNA was measured by quantitative real-time RT-PCR and quantified using the 2(-Delta Delta C)T method. Dose-dependent increases in Camk1g1 mRNA expression were observed for both pyrethroids at 6 h. In addition, a dose-dependent increase in Camk1g2 was observed at 6 h although it was very small in magnitude. The increases in Camk1g1 expression for deltamethrin and permethrin peak between 3 and 6 h postexposure and returns to control levels by 9 h. There was no increase in CAMK1G1 protein as measured with Western blots. The present data demonstrate that pyrethroid-induced changes in Camk1g are driven mainly by increased expression of the Camk1g1 splice variant.