Amar V. Singh

Predictive Models of Prenatal Developmental Toxicity from ToxCast High-Throughput Screening Data

Environmental Protection Agencys ToxCast project is profiling the in vitro bioactivity of chemicals to assess pathway-level and cell-based signatures that correlate with observed in vivo toxicity. We hypothesized that developmental toxicity in guideline animal studies captured in the ToxRefDB database would correlate with cell-based and cell-free in vitro high-throughput screening (HTS) data to reveal meaningful mechanistic relationships and provide models identifying chemicals with the potential to cause developmental toxicity. To test this hypothesis, we built statistical associations based on HTS and in vivo developmental toxicity data from ToxRefDB. Univariate associations were used to filter HTS assays based on statistical correlation with distinct in vivo endpoint. This revealed 423 total associations with distinctly different patterns for rat (301 associations) and rabbit (122 associations) across multiple HTS assay platforms. From these associations, linear discriminant analysis with cross-validation was used to build the models. Species-specific models of predicted developmental toxicity revealed strong balanced accuracy (> 70%) and unique correlations between assay targets such as transforming growth factor beta, retinoic acid receptor, and G-protein-coupled receptor signaling in the rat and inflammatory signals, such as interleukins (IL) (IL1a and IL8) and chemokines (CCL2), in the rabbit. Species-specific toxicity endpoints were associated with one another through common Gene Ontology biological processes, such as cleft palate to urogenital defects through placenta and embryonic development. This work indicates the utility of HTS assays for developing pathway-level models predictive of developmental toxicity.

Modulation of TLR2 Protein Expression by miR-105 in Human Oral Keratinocytes

Mammalian biological processes such as inflammation, involve regulation of hundreds of genes controlling onset and termination. MicroRNAs (miRNAs) can translationally repress target mRNAs and regulate innate immune responses. Our model system comprised primary human keratinocytes, which exhibited robust differences in inflammatory cytokine production (interleukin-6 and tumor necrosis factor-α) following specific Toll-like receptor 2 and 4 (TLR-2/TLR-4) agonist challenge. We challenged these primary cells with Porphyromonas gingivalis (a Gram-negative bacterium that triggers TLR-2 and TLR-4) and performed miRNA expression profiling. We identified miRNA (miR)-105 as a modulator of TLR-2 protein translation in human gingival keratinocytes. There was a strong inverse correlation between cells that had high cytokine responses following TLR-2 agonist challenge and miR-105 levels. Knock-in and knock-down of miR-105 confirmed this inverse relationship. In silico analysis predicted that miR-105 had complementarity for TLR-2 mRNA, and the luciferase reporter assay verified this. Further understanding of the role of miRNA in host responses may elucidate disease susceptibility and suggest new anti-inflammatory therapeutics.

Reprogramming of genetic networks during initiation of the Fetal Alcohol Syndrome

Fetal Alcohol Spectrum Disorders (FASD) are birth defects that result from maternal alcohol use. We used a non a priori approach to prioritize candidate pathways during alcohol‐induced teratogenicity in early mouse embryos. Two C57BL/6 substrains (B6J, B6N) served as the basis for study. Dosing pregnant dams with alcohol (2× 2.9 g/kg ethanol spaced 4 hr on day 8) induced FASD in B6J at a higher incidence than B6N embryos. Counter‐exposure to PK11195 (4 mg/kg) significantly protected B6J embryos but slightly promoted FASD in B6N embryos. Microarray transcript profiling was performed on the embryonic headfold 3 hr after the first maternal alcohol injection (GEO data series accession GSE1074). This analysis revealed metabolic and cellular reprogramming that was substrain‐specific and/or PK11195‐dependent. Mapping ethanol‐responsive KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways revealed down‐regulation of ribosomal proteins and proteasome, and up‐regulation of glycolysis and pentose phosphate pathway in B6N embryos; and significant up‐regulation of tight junction, focal adhesion, adherens junction, and regulation of the actin cytoskeleton (and near‐significant up‐regulation of Wnt signaling and apoptosis) pathways in both substrains. Expression networks constructed computationally from these altered genes identified entry points for EtOH at several hubs (MAPK1, ALDH3A2, CD14, PFKM, TNFRSF1A, RPS6, IGF1, EGFR, PTEN) and for PK11195 at AKT1. Our findings are consistent with the growing view that developmental exposure to alcohol alters common signaling pathways linking receptor activation to cytoskeletal reorganization. The programmatic shift in cell motility and metabolic capacity further implies cell signals and responses that are integrated by the mitochondrial recognition site for PK11195. Developmental Dynamics 236:613–631, 2007.

Activity profiles of 309 ToxCast™ chemicals evaluated across 292 biochemical targets☆

Understanding the potential health risks posed by environmental chemicals is a significant challenge elevated by the large number of diverse chemicals with generally uncharacterized exposures, mechanisms, and toxicities. The present study is a performance evaluation and critical analysis of assay results for an array of 292 high-throughput cell-free assays aimed at preliminary toxicity evaluation of 320 environmental chemicals in EPAs ToxCast™ project (Phase I). The chemicals (309 unique, 11 replicates) were mainly precursors or the active agent of commercial pesticides, for which a wealth of in vivo toxicity data is available. Biochemical HTS (high-throughput screening) profiled cell and tissue extracts using semi-automated biochemical and pharmacological methodologies to evaluate a subset of G-protein coupled receptors (GPCRs), CYP450 enzymes (CYPs), kinases, phosphatases, proteases, HDACs, nuclear receptors, ion channels, and transporters. The primary screen tested all chemicals at a relatively high concentration 25 μM concentration (or 10 μM for CYP assays), and a secondary screen re-tested 9132 chemical-assay pairs in 8-point concentration series from 0.023 to 50 μM (or 0.009-20 μM for CYPs). Mapping relationships across 93,440 chemical-assay pairs based on half-maximal activity concentration (AC50) revealed both known and novel targets in signaling and metabolic pathways. The primary dataset, summary data and details on quality control checks are available for download at http://www.epa.gov/ncct/toxcast/.

Profiling the activity of environmental chemicals in prenatal developmental toxicity studies using the U.S. EPA's ToxRefDB.

As the primary source for regulatory developmental toxicity information, prenatal studies characterize maternal effects and fetal endpoints including malformations, resorptions, and fetal weight reduction. Results from 383 rat and 368 rabbit prenatal studies on 387 chemicals, mostly pesticides, were entered into the U.S. Environmental Protection Agencys (EPA) Toxicity Reference Database (ToxRefDB) using harmonized terminology. An initial assessment of these data was performed with the goal of profiling environmental chemicals based on maternal and fetal endpoints for anchoring in vitro data provided in the EPAs ToxCast research program. Using 30 years worth of standard prenatal studies, maternal and fetal effects were culled from the database and analyzed by target-description fields and lowest effect levels (LELs). Focusing on inter-species comparison, the complexity of fetal target organ response to maternal dosing with environmental chemicals during the period of major organogenesis revealed hierarchical relationships. Of 283 chemicals tested in both species, 53 chemicals (18.7%) had LELs on development (dLEL) that were either specific, with no maternal toxicity (mLEL), or sensitive (dLEL<mLEL) to exposure in one species or another. The primary expressions of developmental toxicity in pregnant rats were fetal weight reduction, skeletal variations and abnormalities, and fetal urogenital defects. General pregnancy/fetal losses were over-represented in the rabbit as were structural malformations to the visceral body wall and CNS. Based upon administered doses, there was a clear hierarchy to the sensitivity and specificity of dLELs in comparing species, with rat development being more sensitive with regards to the number of endpoints affected and the number of active chemicals. Many of these relationships are consistent with previous database studies of developmental toxicology, indicating that they are driven by the biology of the test species. This novel data model provides an important public resource for cross-scale modeling and predictive understanding of developmental processes and toxicities.

Using in Vitro High Throughput Screening Assays to Identify Potential Endocrine-Disrupting Chemicals

Background: Over the past 20 years, an increased focus on detecting environmental chemicals that pose a risk of adverse effects due to endocrine disruption has driven the creation of the U.S. Environmental Protection Agency (EPA) Endocrine Disruptor Screening Program (EDSP). Thousands of chemicals are subject to the EDSP; thus, processing these chemicals using current test batteries could require millions of dollars and decades. A need for increased throughput and efficiency motivated the development of methods using in vitro high throughput screening (HTS) assays to prioritize chemicals for EDSP Tier 1 screening (T1S). Objective: In this study we used U.S. EPA ToxCast HTS assays for estrogen, androgen, steroidogenic, and thyroid-disrupting mechanisms to classify compounds and compare ToxCast results to in vitro and in vivo data from EDSP T1S assays. Method: We implemented an iterative model that optimized the ability of endocrine-related HTS assays to predict components of EDSP T1S and related results. Balanced accuracy was used as a measure of model performance. Results: ToxCast estrogen receptor and androgen receptor assays predicted the results of relevant EDSP T1S assays with balanced accuracies of 0.91 (p < 0.001) and 0.92 (p < 0.001), respectively. Uterotrophic and Hershberger assay results were predicted with balanced accuracies of 0.89 (p < 0.001) and 1 (p < 0.001), respectively. Models for steroidogenic and thyroid-related effects could not be developed with the currently published ToxCast data. Conclusions: Overall, results suggest that current ToxCast assays can accurately identify chemicals with potential to interact with the estrogenic and androgenic pathways, and could help prioritize chemicals for EDSP T1S assays.

Environmental Impact on Vascular Development Predicted by High-Throughput Screening

Background: Understanding health risks to embryonic development from exposure to environmental chemicals is a significant challenge given the diverse chemical landscape and paucity of data for most of these compounds. High-throughput screening (HTS) in the U.S. Environmental Protection Agency (EPA) ToxCast™ project provides vast data on an expanding chemical library currently consisting of > 1,000 unique compounds across > 500 in vitro assays in phase I (complete) and Phase II (under way). This public data set can be used to evaluate concentration-dependent effects on many diverse biological targets and build predictive models of prototypical toxicity pathways that can aid decision making for assessments of human developmental health and disease. Objective: We mined the ToxCast phase I data set to identify signatures for potential chemical disruption of blood vessel formation and remodeling. Methods: ToxCast phase I screened 309 chemicals using 467 HTS assays across nine assay technology platforms. The assays measured direct interactions between chemicals and molecular targets (receptors, enzymes), as well as downstream effects on reporter gene activity or cellular consequences. We ranked the chemicals according to individual vascular bioactivity score and visualized the ranking using ToxPi (Toxicological Priority Index) profiles. Results: Targets in inflammatory chemokine signaling, the vascular endothelial growth factor pathway, and the plasminogen-activating system were strongly perturbed by some chemicals, and we found positive correlations with developmental effects from the U.S. EPA ToxRefDB (Toxicological Reference Database) in vivo database containing prenatal rat and rabbit guideline studies. We observed distinctly different correlative patterns for chemicals with effects in rabbits versus rats, despite derivation of in vitro signatures based on human cells and cell-free biochemical targets, implying conservation but potentially differential contributions of developmental pathways among species. Follow-up analysis with antiangiogenic thalidomide analogs and additional in vitro vascular targets showed in vitro activity consistent with the most active environmental chemicals tested here. Conclusions: We predicted that blood vessel development is a target for environmental chemicals acting as putative vascular disruptor compounds (pVDCs) and identified potential species differences in sensitive vascular developmental pathways.

Gingival epithelial cells heterozygous for Toll-like receptor 4 polymorphisms Asp299Gly and Thr399Ile are hypo-responsive to Porphyromonas gingivalis

The Toll-like receptor (TLR)4 is the major sensor for bacterial lipopolysaccharide and its two common co-segregating polymorphisms, Asp299Gly and Thr399Ile, which occur at a frequency of between 6 and 10%, have been associated with infectious diseases, LPS hypo-responsiveness and cardiovascular disease. Porphyromonas gingivalis is a Gram-negative bacterium implicated in chronic periodontitis and is a known TLR4 and TLR2 agonist. We obtained two gingival epithelial cell primary cultures from subjects heterozygous for the TLR4 polymorphism Asp299Gly and compared response characteristics with similar cells from patients (four) with the wild-type TLR4 genes. Cytokine responses and transcriptome profiles of gingival epithelial cell primary culture cells to TNFα challenge were similar for all primary epithelial cell cultures. P. gingivalis challenge, however, gave markedly different responses for Asp299Gly heterozygous and wild-type epithelial cell cultures. The epithelial cells heterozygous for the TLR4 polymorphism Asp299Gly were functionally hypo-responsive, evidenced by differences in BD-2 mRNA expression, mRNA response profile by microarray analysis and by pro-inflammatory and chemokine cytokines at the protein and mRNA level. These findings emphasize variance in human epithelial cell TLRs, linked with Asp299Gly carriage, which results in a hypo-responsive epithelial cell phenotype less susceptible to Gram-negative diseases and associated systemic conditions.

Evaluation of 309 Environmental Chemicals Using a Mouse Embryonic Stem Cell Adherent Cell Differentiation and Cytotoxicity Assay

The vast landscape of environmental chemicals has motivated the need for alternative methods to traditional whole-animal bioassays in toxicity testing. Embryonic stem (ES) cells provide an in vitro model of embryonic development and an alternative method for assessing developmental toxicity. Here, we evaluated 309 environmental chemicals, mostly food-use pesticides, from the ToxCast™ chemical library using a mouse ES cell platform. ES cells were cultured in the absence of pluripotency factors to promote spontaneous differentiation and in the presence of DMSO-solubilized chemicals at different concentrations to test the effects of exposure on differentiation and cytotoxicity. Cardiomyocyte differentiation (α,β myosin heavy chain; MYH6/MYH7) and cytotoxicity (DRAQ5™/Sapphire700™) were measured by In-Cell Western™ analysis. Half-maximal activity concentration (AC50) values for differentiation and cytotoxicity endpoints were determined, with 18% of the chemical library showing significant activity on either endpoint. Mining these effects against the ToxCast Phase I assays (∼500) revealed significant associations for a subset of chemicals (26) that perturbed transcription-based activities and impaired ES cell differentiation. Increased transcriptional activity of several critical developmental genes including BMPR2, PAX6 and OCT1 were strongly associated with decreased ES cell differentiation. Multiple genes involved in reactive oxygen species signaling pathways (NRF2, ABCG2, GSTA2, HIF1A) were strongly associated with decreased ES cell differentiation as well. A multivariate model built from these data revealed alterations in ABCG2 transporter was a strong predictor of impaired ES cell differentiation. Taken together, these results provide an initial characterization of metabolic and regulatory pathways by which some environmental chemicals may act to disrupt ES cell growth and differentiation.

Identification of Estrogen Target Genes during Zebrafish Embryonic Development through Transcriptomic Analysis

Estrogen signaling is important for vertebrate embryonic development. Here we have used zebrafish (Danio rerio) as a vertebrate model to analyze estrogen signaling during development. Zebrafish embryos were exposed to 1 µM 17β-estradiol (E2) or vehicle from 3 hours to 4 days post fertilization (dpf), harvested at 1, 2, 3 and 4 dpf, and subjected to RNA extraction for transcriptome analysis using microarrays. Differentially expressed genes by E2-treatment were analyzed with hierarchical clustering followed by biological process and tissue enrichment analysis. Markedly distinct sets of genes were up and down-regulated by E2 at the four different time points. Among these genes, only the well-known estrogenic marker vtg1 was co-regulated at all time points. Despite this, the biological functional categories targeted by E2 were relatively similar throughout zebrafish development. According to knowledge-based tissue enrichment, estrogen responsive genes were clustered mainly in the liver, pancreas and brain. This was in line with the developmental dynamics of estrogen-target tissues that were visualized using transgenic zebrafish containing estrogen responsive elements driving the expression of GFP (Tg(5xERE:GFP)). Finally, the identified embryonic estrogen-responsive genes were compared to already published estrogen-responsive genes identified in male adult zebrafish (Gene Expression Omnibus database). The expressions of a few genes were co-regulated by E2 in both embryonic and adult zebrafish. These could potentially be used as estrogenic biomarkers for exposure to estrogens or estrogenic endocrine disruptors in zebrafish. In conclusion, our data suggests that estrogen effects on early embryonic zebrafish development are stage- and tissue- specific.