Developmental neurotoxicity (DNT) and many forms of reproductive toxicity (RT) often manifest themselves in functional deficits that are not necessarily based on cell death, but rather on minor changes relating to cell differentiation or communication. The fields of DNT/RT would greatly benefit from in vitro tests that allow the identification of toxicant-induced changes of the cellular proteostasis, or of its underlying transcriptome network. Therefore, the ‘human embryonic stem cell (hESC)-derived novel alternative test systems (ESNATS)’ European commission research project established RT tests based on defined differentiation protocols of hESC and their progeny. Valproic acid (VPA) and methylmercury (MeHg) were used as positive control compounds to address the following fundamental questions: (1) Does transcriptome analysis allow discrimination of the two compounds? (2) How does analysis of enriched transcription factor binding sites (TFBS) and of individual probe sets (PS) distinguish between test systems? (3) Can batch effects be controlled? (4) How many DNA microarrays are needed? (5) Is the highest non-cytotoxic concentration optimal and relevant for the study of transcriptome changes? VPA triggered vast transcriptional changes, whereas MeHg altered fewer transcripts. To attenuate batch effects, analysis has been focused on the 500 PS with highest variability. The test systems differed significantly in their responses (<20 % overlap). Moreover, within one test system, little overlap between the PS changed by the two compounds has been observed. However, using TFBS enrichment, a relatively large ‘common response’ to VPA and MeHg could be distinguished from ‘compound-specific’ responses. In conclusion, the ESNATS assay battery allows classification of human DNT/RT toxicants on the basis of their transcriptome profiles.Electronic supplementary materialThe online version of this article (doi:10.1007/s00204-012-0967-3) contains supplementary material, which is available to authorized users.

Information on design principles governing transcriptome changes upon transition from safe to hazardous drug concentrations or from tolerated to cytotoxic drug levels are important for the application of toxicogenomics data in developmental toxicology. Here, we tested the effect of eight concentrations of valproic acid (VPA; 25–1000 μM) in an assay that recapitulates the development of human embryonic stem cells to neuroectoderm. Cells were exposed to the drug during the entire differentiation process, and the number of differentially regulated genes increased continuously over the concentration range from zero to about 3000. We identified overrepresented transcription factor binding sites (TFBS) as well as superordinate cell biological processes, and we developed a gene ontology (GO) activation profiler, as well as a two-dimensional teratogenicity index. Analysis of the transcriptome data set by the above biostatistical and systems biology approaches yielded the following insights: (i) tolerated (≤25 μM), deregulated/teratogenic (150–550 μM), and cytotoxic (≥800 μM) concentrations could be differentiated. (ii) Biological signatures related to the mode of action of VPA, such as protein acetylation, developmental changes, and cell migration, emerged from the teratogenic concentrations range. (iii) Cytotoxicity was not accompanied by signatures of newly emerging canonical cell death/stress indicators, but by catabolism and decreased expression of cell cycle associated genes. (iv) Most, but not all of the GO groups and TFBS seen at the highest concentrations were already overrepresented at 350–450 μM. (v) The teratogenicity index reflected this behavior, and thus differed strongly from cytotoxicity. Our findings suggest the use of the highest noncytotoxic drug concentration for gene array toxicogenomics studies, as higher concentrations possibly yield wrong information on the mode of action, and lower drug levels result in decreased gene expression changes and thus a reduced power of the study.

The currently available techniques for the safety evaluation of candidate drugs are usually cost-intensive and time-consuming and are often insufficient to predict human relevant cardiotoxicity. The purpose of this study was to develop an in vitro repeated exposure toxicity methodology allowing the identification of predictive genomics biomarkers of functional relevance for drug-induced cardiotoxicity in human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs). The hiPSC-CMs were incubated with 156 nM doxorubicin, which is a well-characterized cardiotoxicant, for 2 or 6 days followed by washout of the test compound and further incubation in compound-free culture medium until day 14 after the onset of exposure. An xCELLigence Real-Time Cell Analyser was used to monitor doxorubicin-induced cytotoxicity while also monitoring functional alterations of cardiomyocytes by counting of the beating frequency of cardiomyocytes. Unlike single exposure, repeated doxorubicin exposure resulted in long-term arrhythmic beating in hiPSC-CMs accompanied by significant cytotoxicity. Global gene expression changes were studied using microarrays and bioinformatics tools. Analysis of the transcriptomic data revealed early expression signatures of genes involved in formation of sarcomeric structures, regulation of ion homeostasis and induction of apoptosis. Eighty-four significantly deregulated genes related to cardiac functions, stress and apoptosis were validated using real-time PCR. The expression of the 84 genes was further studied by real-time PCR in hiPSC-CMs incubated with daunorubicin and mitoxantrone, further anthracycline family members that are also known to induce cardiotoxicity. A panel of 35 genes was deregulated by all three anthracycline family members and can therefore be expected to predict the cardiotoxicity of compounds acting by similar mechanisms as doxorubicin, daunorubicin or mitoxantrone. The identified gene panel can be applied in the safety assessment of novel drug candidates as well as available therapeutics to identify compounds that may cause cardiotoxicity.Electronic supplementary materialThe online version of this article (doi:10.1007/s00204-015-1623-5) contains supplementary material, which is available to authorized users.

Embryonic development can be partially recapitulated in vitro by differentiating human embryonic stem cells (hESCs). Thalidomide is a developmental toxicant in vivo and acts in a species-dependent manner. Besides its therapeutic value, thalidomide also serves as a prototypical model to study teratogenecity. Although many in vivo and in vitro platforms have demonstrated its toxicity, only a few test systems accurately reflect human physiology. We used global gene expression and proteomics profiling (two dimensional electrophoresis (2DE) coupled with Tandem Mass spectrometry) to demonstrate hESC differentiation and thalidomide embryotoxicity/teratogenecity with clinically relevant dose(s). Proteome analysis showed loss of POU5F1 regulatory proteins PKM2 and RBM14 and an over expression of proteins involved in neuronal development (such as PAK2, PAFAH1B2 and PAFAH1B3) after 14 days of differentiation. The genomic and proteomic expression pattern demonstrated differential expression of limb, heart and embryonic development related transcription factors and biological processes. Moreover, this study uncovered novel possible mechanisms, such as the inhibition of RANBP1, that participate in the nucleocytoplasmic trafficking of proteins and inhibition of glutathione transferases (GSTA1, GSTA2), that protect the cell from secondary oxidative stress. As a proof of principle, we demonstrated that a combination of transcriptomics and proteomics, along with consistent differentiation of hESCs, enabled the detection of canonical and novel teratogenic intracellular mechanisms of thalidomide.

Test systems to identify developmental toxicants are urgently needed. A combination of human stem cell technology and transcriptome analysis was to provide a proof of concept that toxicants with a related mode of action can be identified and grouped for read-across. We chose a test system of developmental toxicity, related to the generation of neuroectoderm from pluripotent stem cells (UKN1), and exposed cells for 6 days to the histone deacetylase inhibitors (HDACi) valproic acid, trichostatin A, vorinostat, belinostat, panobinostat and entinostat. To provide insight into their toxic action, we identified HDACi consensus genes, assigned them to superordinate biological processes and mapped them to a human transcription factor network constructed from hundreds of transcriptome data sets. We also tested a heterogeneous group of ‘mercurials’ (methylmercury, thimerosal, mercury(II)chloride, mercury(II)bromide, 4-chloromercuribenzoic acid, phenylmercuric acid). Microarray data were compared at the highest non-cytotoxic concentration for all 12 toxicants. A support vector machine (SVM)-based classifier predicted all HDACi correctly. For validation, the classifier was applied to legacy data sets of HDACi, and for each exposure situation, the SVM predictions correlated with the developmental toxicity. Finally, optimization of the classifier based on 100 probe sets showed that eight genes (F2RL2, TFAP2B, EDNRA, FOXD3, SIX3, MT1E, ETS1 and LHX2) are sufficient to separate HDACi from mercurials. Our data demonstrate how human stem cells and transcriptome analysis can be combined for mechanistic grouping and prediction of toxicants. Extension of this concept to mechanisms beyond HDACi would allow prediction of human developmental toxicity hazard of unknown compounds with the UKN1 test system.Electronic supplementary materialThe online version of this article (doi:10.1007/s00204-015-1573-y) contains supplementary material, which is available to authorized users.