Human Pluripotent Stem Cell Based Developmental Toxicity Assays for Chemical Safety Screening and Systems Biology Data Generation

Efficient protocols to differentiate human pluripotent stem cells to various tissues in combination with -omics technologies opened up new horizons for in ...

Technical university of Dortmund

Technical University of Dortmund

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.

Human embryonic stem cell-derived test systems for developmental neurotoxicity: a transcriptomics approach.

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.

Design principles of concentration-dependent transcriptome deviations in drug-exposed differentiating stem cells.

The superordinate principles governing the transcriptome response of differentiating cells exposed to drugs are still unclear. Often, it is assumed that toxicogenomics data reflect the immediate mode of action (MoA) of drugs. Alternatively, transcriptome changes could describe altered differentiation states as indirect consequence of drug exposure. We used here the developmental toxicants valproate and trichostatin A to address this question. Neurally differentiating human embryonic stem cells were treated for 6 days. Histone acetylation (primary MoA) increased quickly and returned to baseline after 48 h. Histone H3 lysine methylation at the promoter of the neurodevelopmental regulators PAX6 or OTX2 was increasingly altered over time. Methylation changes remained persistent and correlated with neurodevelopmental defects and with effects on PAX6 gene expression, also when the drug was washed out after 3-4 days. We hypothesized that drug exposures altering only acetylation would lead to reversible transcriptome changes (indicating MoA), and challenges that altered methylation would lead to irreversible developmental disturbances. Data from pulse-chase experiments corroborated this assumption. Short drug treatment triggered reversible transcriptome changes; longer exposure disrupted neurodevelopment. The disturbed differentiation was reflected by an altered transcriptome pattern, and the observed changes were similar when the drug was washed out during the last 48 h. We conclude that transcriptome data after prolonged chemical stress of differentiating cells mainly reflect the altered developmental stage of the model system and not the drug MoA. We suggest that brief exposures, followed by immediate analysis, are more suitable for information on immediate drug responses and the toxicity MoA.

From transient transcriptome responses to disturbed neurodevelopment: role of histone acetylation and methylation as epigenetic switch between reversible and irreversible drug effects.

A long-term goal of numerous research projects is to identify biomarkers for in vitro systems predicting toxicity in vivo. Often, transcriptomics data are used to identify candidates for further evaluation. However, a systematic directory summarizing key features of chemically influenced genes in human hepatocytes is not yet available. To bridge this gap, we used the Open TG-GATES database with Affymetrix files of cultivated human hepatocytes incubated with chemicals, further sets of gene array data with hepatocytes from human donors generated in this study, and publicly available genome-wide datasets of human liver tissue from patients with non-alcoholic steatohepatitis (NASH), cirrhosis, and hepatocellular cancer (HCC). After a curation procedure, expression data of 143 chemicals were included into a comprehensive biostatistical analysis. The results are summarized in the publicly available toxicotranscriptomics directory ( http://wiki.toxbank.net/toxicogenomics-map/ ) which provides information for all genes whether they are up- or downregulated by chemicals and, if yes, by which compounds. The directory also informs about the following key features of chemically influenced genes: (1) Stereotypical stress response. When chemicals induce strong expression alterations, this usually includes a complex but highly reproducible pattern named 'stereotypical response.' On the other hand, more specific expression responses exist that are induced only by individual compounds or small numbers of compounds. The directory differentiates if the gene is part of the stereotypical stress response or if it represents a more specific reaction. (2) Liver disease-associated genes. Approximately 20 % of the genes influenced by chemicals are up- or downregulated, also in liver disease. Liver disease genes deregulated in cirrhosis, HCC, and NASH that overlap with genes of the aforementioned stereotypical chemical stress response include CYP3A7, normally expressed in fetal liver; the phase II metabolizing enzyme SULT1C2; ALDH8A1, known to generate the ligand of RXR, one of the master regulators of gene expression in the liver; and several genes involved in normal liver functions: CPS1, PCK1, SLC2A2, CYP8B1, CYP4A11, ABCA8, and ADH4. (3) Unstable baseline genes. The process of isolating and the cultivation of hepatocytes was sufficient to induce some stress leading to alterations in the expression of genes, the so-called unstable baseline genes. (4) Biological function. Although more than 2,000 genes are transcriptionally influenced by chemicals, they can be assigned to a relatively small group of biological functions, including energy and lipid metabolism, inflammation and immune response, protein modification, endogenous and xenobiotic metabolism, cytoskeletal organization, stress response, and DNA repair. In conclusion, the introduced toxicotranscriptomics directory offers a basis for a rationale choice of candidate genes for biomarker evaluation studies and represents an easy to use source of background information on chemically influenced genes.

Toxicogenomics directory of chemically exposed human hepatocytes.

Endometrial carcinoma differential 3 (EDI3) was the first member of the glycerophosphodiesterase (GDE) protein family shown to be associated with cancer. Our initial work demonstrated that endometrial and ovarian cancer patients with primary tumors overexpressing EDI3 had a higher risk of developing metastasis and decreased survival. Further analysis indicated that EDI3 cleaves glycerophosphocholine to choline and glycerol-3-phosphate, increases the levels of active PKC, and enhances the migratory activity of tumor cells. Despite these initial findings, EDI3 remained mainly uncharacterized. Therefore, to obtain an overview of processes in which EDI3 may be involved, gene array analysis was performed using MCF-7 breast cancer cells after EDI3 knockdown compared with a non-targeting control siRNA. Several biological motifs were altered, including an enrichment of genes involved in integrin-mediated signaling. More specifically, silencing of EDI3 in MCF-7 and OVCAR-3 cells was associated with reduced expression of the key receptor subunit integrin β1, leading to decreased cell attachment and spreading accompanied by delayed formation of cell protrusions. To confirm these results, we stably overexpressed EDI3 in MCF-7 cells which led to elevated integrin β1 expression associated with enhanced cell attachment and spreading - two processes critical for metastasis. In conclusion, our data provide further insight into the role of EDI3 during cancer progression.

Eugen Rempel

Department of Statistics

Technical university of Dortmund

Human Pluripotent Stem Cell Based Developmental Toxicity Assays for Chemical Safety Screening and Systems Biology Data Generation

Eugen Rempel

.css-o250i3{font-size:18px;font-family:Open Sans,sans-serif;line-height:21.6px;}Department of Statistics

Technical university of Dortmund

Human Pluripotent Stem Cell Based Developmental Toxicity Assays for Chemical Safety Screening and Systems Biology Data Generation

Department of Statistics