Abstract: EPA's 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 endpoints. 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 (BA > 70%), and unique correlations between assay targets such as transforming growth factor beta (TGFβ), retinoic acid receptor (RAR), and G-protein-coupled receptor (GPCR) signaling in the rat and inflammatory signals, such as interleukins (IL1a and IL8) and chemokines (CCL2) in the rabbit. Species specific toxicity endpoints were associated with one another through common Gene Ontology (GO) 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.
Abstract: 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 (ACâ â) 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.
Abstract: 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 EPA's 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/.
Abstract: 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 EPA's ToxCastTM project provides vast data on a chemical library consisting of 309 unique compounds across over 500 in vitro assays. Objective: This public dataset can be used to evaluate concentration-dependent effects on many diverse biological targets and build predictive models of prototypical toxicity pathways that can help in decision-making for assessments to human developmental health and disease. Methods: The present study mined the ToxCast dataset and identified signatures for potential chemical disruption of blood vessel formation and remodeling. Results: Targets in inflammatory chemokine signaling, the vascular endothelial growth factor pathway, and the plasminogen activating system were strongly perturbed, and positive correlations were found with developmental effects from the ToxRefDB in vivo database containing prenatal rat and rabbit guideline studies. Distinctly different correlative patterns were observed 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 anti-angiogenic thalidomide analogs and additional in vitro vascular targets showed in vitro activity consistent with the most active environmental chemicals tested here. Conclusions: These predictions implicate the embryonic vasculature as a target for environmental chemicals acting as putative Vascular Disruptor Compounds (pVDCs) and illuminates potential species differences in sensitive vascular developmental pathways.
Abstract: The present study was performed to determine experimental conditions for thalidomide induction of fetal malformations and to understand the molecular mechanisms underlying thalidomide teratogenicity in cynomolgus monkeys. Cynomolgus monkeys were orally administered thalidomide at 15 or 20mg/kg-d on days 26-28 of gestation, and fetuses were examined on day 100-102 of gestation. Limb defects such as micromelia/amelia, paw/foot hyperflexion, polydactyly, syndactyly, and brachydactyly were observed in seven of eight fetuses. Cynomolgus monkeys were orally administered thalidomide at 20mg/kg on day 26 of gestation, and whole embryos were removed from the dams 6h after administration. Three embryos each were obtained from the thalidomide-treated and control groups. Total RNA was isolated from individual embryos, amplified to biotinylated cRNA and hybridized to a custom Non-Human Primate (NHP) GeneChip((R)) Array. Altered genes were clustered into genes that were up-regulated (1281 genes) and down-regulated (1081 genes) in thalidomide-exposed embryos. Functional annotation by Gene Ontology (GO) categories revealed up-regulation of actin cytoskeletal remodeling and insulin signaling, and down-regulation of pathways for vasculature development and the inflammatory response. These findings show that thalidomide exposure perturbs a general program of morphoregulatory processes in the monkey embryo. Bioinformatics analysis of the embryonic transcriptome following maternal thalidomide exposure has now identified many key pathways implicated in thalidomide embryopathy, and has also revealed some novel processes that can help unravel the mechanism of this important developmental phenotype.
Abstract: 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 Agency's (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 EPA's 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.
Abstract: 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-alpha) 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.
Abstract: Development of the fetal gubernaculum is a prerequisite for testicular descent and dependent on insulin-like 3 and androgen, but knowledge of downstream effectors is limited. We analyzed transcript profiles in gubernaculum and testis to address changes occurring during normal and abnormal testicular descent in Long Evans wild-type (wt) and cryptorchid (orl) fetuses. Total RNA from male wt and orl gubernacula (gestational days [GD]18-20), wt female gubernacula (GD18), and testis (GD17 and 19) was hybridized to Affymetrix GeneChips. Statistical analysis of temporal, gender, and strain-specific differences in gene expression was performed with the use of linear models analysis with empirical Bayes statistics and analysis of variance (gubernaculum) and linear analysis (testis). Overrepresented common gene ontology functional categories and pathways were identified in groups of differentially expressed genes with the Database for Annotation, Visualization, and Integrated Discovery. Transcript profiles were dynamic in wt males between GD18-19 and GD20, comparatively static in orl GD18-20 gubernaculum, and similar in wt and orl testis. Functional analysis of differentially expressed genes in wt and orl gubernaculum identified categories related to metabolism, cellular biogenesis, small GTPase-mediated signal transduction, cytoskeleton, muscle development, and insulin signaling. Genes involved in androgen receptor signaling, regulated by androgens, or both were overrepresented in differentially expressed gubernaculum and testis gene groups. Quantitative reverse transcription polymerase chain reaction (RT-PCR) confirmed differential expression of genes related to muscle development, including Myog, Tnnt2, Fst, Igf1, Igfbp5, Id2, and Msx1. These data suggest that the orl mutation results in a primary gubernacular defect that affects muscle development and cytoskeletal function and might alter androgen-regulated pathways.
Abstract: Microarrays provide a means of studying expression level of tens of thousands of genes by providing one or more oligonucleotide probe(s) for each transcript studied. Affymetrix(R) GeneChiptrade mark platforms historically pair each 25-base perfect match (PM) probe with a mismatch probe (MM) differing by a complementary base located in the 13(th) position to quantify and deflate effects of cross-hybridization. Analytical routines for analyzing these arrays take into account difference in expression levels of MM and PM probes to determine which ones are useful for further study. If a single nucleotide polymorphism (SNP) occurs at the 13(th) base, a probe with a higher MM expression level may be incorrectly omitted. In order to examine SNP affects on PM and MM expression levels, known human SNPs from dbSNP were mapped to probe sets within the Affymetrix(R) HG-U133A platform. Probe sets containing one or more probe pairs with a single SNP at the 13(th) position were extracted. A set of twelve microarray experiments were analyzed for the PM and MM expression levels for these probe sets. Over 6,000,000 human SNPs and their flanking regions were extracted from dbSNP. These sequences were aligned against each of the 247,965 probe pair sequences from the Affymetrix(R) HG-U133A platform. A total of 915 probe sets containing a single probe sequence with a SNP mapped to the 13(th) base were extracted. A subset containing 166 probe sets result in complementary base SNPs. Comparison of gene expression levels for the SNP to non-SNP PM and MM probes does not yield a significant difference using chi2 analysis. Thus, omission of probes with MM expression levels higher than PM expression levels does not appear to result in a loss of information concerning SNPs for these regions.
Abstract: To perform a large-scale gene profiling of the liver in a mouse model of fatty liver induced by high carbohydrate (sucrose) diet (HCD) to gain a deeper insight into potential mechanisms of diet-induced hepatic steatosis.
Abstract: Cells in the developing embryo must integrate complex signals from the genome and environment to make decisions about their behavior or fate. The ability to understand the fundamental biology of the decision-making process, and how these decisions may go awry during abnormal development, requires a systems biology paradigm. Presently, the ability to build models with predictive capability in birth defects research is constrained by an incomplete understanding of the fundamental parameters underlying embryonic susceptibility, sensitivity, and vulnerability. Key developmental milestones must be parameterized in terms of system structure and dynamics, the relevant control methods, and the overall design logic of metabolic and regulatory networks. High-content data from genome-based studies provide some comprehensive coverage of these operational processes but a key research challenge is data integration. Analysis can be facilitated by data management resources and software to reveal the structure and function of bionetwork motifs potentially associated with an altered developmental phenotype. Borrowing from applied mathematics and artificial intelligence, we conceptualize a system that can help address the new challenges posed by the transformation of birth defects research into a data-driven science.
Abstract: 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 (2x 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.
Abstract: Monitoring global gene expression provides insight into how genes and regulatory signals work together to guide embryo development. The fields of developmental biology and teratology are now confronted with the need for automated access to a reference library of gene-expression signatures that benchmark programmed (genetic) and adaptive (environmental) regulation of the embryonic transcriptome. Such a library must be constructed from highly-distributed microarray data. Birth Defects Systems Manager (BDSM), an open access knowledge management system, provides custom software to mine public microarray data focused on developmental health and disease. The present study describes tools for seamless data integration in the BDSM library (MetaSample, MetaChip, CIAeasy) using the QueryBDSM module. A field test of the prototype was run using published microarray data series derived from a variety of laboratories, experiments, microarray platforms, organ systems, and developmental stages. The datasets focused on several developing systems in the mouse embryo, including preimplantation stages, heart and nerve development, testis and ovary development, and craniofacial development. Using BDSM data integration tools, a gene-expression signature for 346 genes was resolved that accurately classified samples by organ system and developmental sequence. The module builds a potential for the BDSM approach to decipher a large number developmental processes through comparative bioinformatics analysis of embryological systems at-risk for specific defects, using multiple scenarios to define the range of probabilities leading from molecular phenotype to clinical phenotype. We conclude that an integrative analysis of global gene-expression of the developing embryo can form the foundation for constructing a reference library of signaling pathways and networks for normal and abnormal regulation of the embryonic transcriptome. These tools are available free of charge from the web-site http://systemsanalysis.louisville.edu requiring only a short registration process.
Abstract: Birth defects and developmental disabilities remain an important public health issue worldwide. With the availability of genomic sequences from a growing number of human and model organisms and the rapid expansion of the public repositories holding large-scale gene expression datasets, a computational systems analysis of developmental toxicology can incorporate this vast digital information toward the realization of predictive models for complex disease. Here we describe the initial design, development and implementation of a Birth Defects Systems Manager (BDSM). The project was motivated by the need for a computational-bioinformatics infrastructure to manage vast digital information from functional genomics and for a new knowledge environment specifically engineered for the analysis of developmental processes and toxicities. Proof-of-concept tested BDSM using meta-analysis of gene expression data collected from different laboratories, technology platforms, and study models. The composite dataset incorporated 232 microarray comparisons of RNA samples by single or dual microarray platforms, cDNA or oligonucleotide based probes, and human or mouse sequence information. Preliminary results identified system-level features in the embryonic transcriptome as it reacted to various developmental-teratological stimuli. BDSM is open access through the worldwide web (http://systemsanalysis.louisville.edu/) and can be integrated with other bioinformatics tools and resources to advance the pace of discovery in birth defects research.
Abstract: Gene expression arrays reveal the potential linkage of altered gene expression with specific adverse effects leading to disease phenotypes. But how closely do microarray data reflect early physiological or pharmacological measures that predict toxic event(s)? To explore this issue, we have undertaken experiments in early mouse embryos exposed to various teratogens during neurulation stages with the aim of correlating large-scale changes in gene expression across the critical period during exposure. This study reports some of the large-scale changes in gene expression that can be detected in the optic rudiment of the developing mouse and rat embryo across the window of development during which the eye is exceedingly sensitive to teratogen-induced micro-/anophthalmia. Microarray analysis was performed on RNA from the headfold or ocular region at the optic vesicle and optic cup stages when the ocular primordium is enriched for Pax-6, a master control gene for eye morphogenesis. Statistical selection of differentially regulated genes and various clustering techniques identified groups of genes in upward or downward trajectories in the normal optic primordium during early eye development in mouse and rat species. We identified 165 genes with significant differential expression during eye development, and a smaller subset of 58 genes that showed a tight correlation between mouse-rat development. Significantly over-represented functional categories included fatty acid metabolism (up-regulated) and glycolysis (down-regulated). From studies such as these that benchmark large-scale gene expression during normal embryonic development, we may be able to identify the panel of biomarkers that best correlate with species differences and the risks for developmental toxicity.
Abstract: The need for a computational bioinformatics infrastructure to manage the vast digital information from functional genomics and proteomics motivated us to develop Birth Defects Systems Manager (BDSM) as an open resource to facilitate analysis and discovery in developmental biology and developmental toxicity. This report describes the design, development and implementation of the data loading module of BDSM, referred to as LoadBDSM. It includes a shared data directory resource that can be granted various levels of security for different research groups or investigators to manage experimental datasets individually or in groups. LoadBDSM allows the upload of data and experiment details using controlled semantics for developmental exposure (toxicant, dosing scenario, intervention), biological sample (species, tissue, stage) and disease outcome (time, risk, phenotype). It adheres to existing controlled vocabulary plus rules of inference (ontologies) for experiment, data and metadata annotations. LoadBDSM extends the capabilities of BDSM to support the emergence of "embryo-formatics" defined here as the data, information and knowledge from genomic sciences applied to, or derived from, an embryological context. This includes, but is not limited to, delineating pathways and biological regulatory networks for specific chemicals or classes of developmental toxicants, developing novel biomarkers indicative of exposure and/or predictive of adverse effects, and integrating modern computing and information technology with data from molecular biology.
