From multi-omics integration towards novel genomic interaction networks to identify key cancer cell line characteristics.

Cancer is a complex disease where cancer cells express epigenetic and transcriptomic mechanisms to promote tumor initiation, progression, and survival. To extract relevant features from the 2019 Cancer Cell Line Encyclopedia (CCLE), a multi-layer nonnegative matrix factorization approach is used. We used relevant feature genes and DNA promoter regions to construct genomic interaction network to study gene-gene and gene-DNA promoter methylation relationships. Here, we identified a set of gene transcripts and methylated DNA promoter regions for different clusters, including one homogeneous lymphoid neoplasms cluster. In this cluster, we found different methylated transcription factors that affect transcriptional activation of EGFR and downstream interactions. Furthermore, the hippo-signaling pathway might not function properly because of DNA hypermethylation and low gene expression of both LATS2 and YAP1. Finally, we could identify a potential dysregulation of the CD28-CD86-CTLA4 axis. Characterizing the interaction of the epigenome and the transcriptome is vital for our understanding of cancer cell line behavior, not only for deepening insights into cancer-related processes but also for future disease treatment and drug development. Here we have identified potential candidates that characterize cancer cell lines, which give insight into the development and progression of cancers.

Valproic acid promotes mitochondrial dysfunction in primary human hepatocytes in vitro; impact of C/EBPα-controlled gene expression.

Valproic acid (VPA) is a frequently prescribed anti-epileptic drug which is known to cause liver toxicity and steatosis through mitochondrial dysfunction. Nevertheless the mechanisms underlying these adverse effects are incompletely understood. In this study, we determined the effect of relatively short (3 h) or prolonged (72 h) exposure to VPA on mitochondrial function in primary human hepatocytes (PHHs). While 3 h VPA exposure did not affect oxygen consumption rates (OCRs) in PHHs, prolonged exposure (24-72 h) significantly reduced basal and maximal OCRs. Given that in particular prolonged VPA exposure is required to cause mitochondrial dysfunction, we investigated gene expression data after VPA exposure for 24, 48, 72 h and 72 h VPA followed by a 72 h washout period. We were able to reduce the comprehensive gene expression changes into a more comprehensible set of 18 TFs that were predicted to be persistently activated after 72 h of VPA exposure. Lentiviral knock-down of one of the candidate TFs, C/EBPα, partly rescued VPA-induced mitochondrial dysfunction. Furthermore, RNA-Seq analysis of shC/EBPα and shGFP control PHHs identified 24 genuine C/EBPα target genes that are regulated in response to prolonged VPA exposure in PHHs. Altogether this provides new insights on the involvement of C/EBPα in driving VPA-induced mitochondrial dysfunction in human liver cells. This hub gene, with its downstream regulators involved in this deregulation, thus represent potential new biomarkers for VPA-induced mitochondrial dysfunction.

Circulating microRNAs as potential biomarkers for psychiatric and neurodegenerative disorders.

Circulating microRNAs (cimiRNAs) are a class of non-encoding RNAs found in bodily fluids such as blood, cerebrospinal fluid (CSF) and tears. CimiRNAs have been implicated as promising biomarkers for central nervous system (CNS) disorders because they are actively secreted as messengers and are profoundly involved in fine-tuning of developmental and differentiation processes. Furthermore, they are attractive biomarkers because they are extremely stable, tissue enriched and can be determined in a quantitative manner. This review aims to provide a comprehensive assessment on the current progress regarding the potential value of cimiRNAs as CNS biomarkers. Within this framework five CNS disorders are explored which share a common pathological hallmark namely cognitive impairment. The CNS disorders include Major depression disorder (MDD), Bipolar disorder (BD), Schizophrenia (SZ), Alzheimer's disease (AD) and Parkinson disease (PD). The similarities and differences between altered cimiRNAs in the different disorders are described. The miR-29 family, miR-34a-5p and miR-132-3p are discussed as common dysregulated cimiRNAs found in the CNS disorders. Furthermore, it is shown that the type of bodily fluid used for measuring cimiRNAs is important as inconsistencies in cimiRNAs expression directions are found when comparing CSF, blood cell-free and blood cell-bound samples.

DynOVis: a web tool to study dynamic perturbations for capturing dose-over-time effects in biological networks.

BACKGROUND: The development of high throughput sequencing techniques provides us with the possibilities to obtain large data sets, which capture the effect of dynamic perturbations on cellular processes. However, because of the dynamic nature of these processes, the analysis of the results is challenging. Therefore, there is a great need for bioinformatics tools that address this problem. RESULTS: Here we present DynOVis, a network visualization tool that can capture dynamic dose-over-time effects in biological networks. DynOVis is an integrated work frame of R packages and JavaScript libraries and offers a force-directed graph network style, involving multiple network analysis methods such as degree threshold, but more importantly, it allows for node expression animations as well as a frame-by-frame view of the dynamic exposure. Valuable biological information can be highlighted on the nodes in the network, by the integration of various databases within DynOVis. This information includes pathway-to-gene associations from ConsensusPathDB, disease-to-gene associations from the Comparative Toxicogenomics databases, as well as Entrez gene ID, gene symbol, gene synonyms and gene type from the NCBI database. CONCLUSIONS: DynOVis could be a useful tool to analyse biological networks which have a dynamic nature. It can visualize the dynamic perturbations in biological networks and allows the user to investigate the changes over time. The integrated data from various online databases makes it easy to identify the biological relevance of nodes in the network. With DynOVis we offer a service that is easy to use and does not require any bioinformatics skills to visualize a network.

Human 3D multicellular microtissues: An upgraded model for the in vitro mechanistic investigation of inflammation-associated drug toxicity.

Inflammation is one of the factors that may increase the sensitivity of hepatic cells to acetaminophen (APAP) induced toxicity. To investigate the mechanisms, we exposed 3-dimensional (3D) Human Liver Microtissues, a co-culture of primary human hepatocytes (PHH) and Kupffer cells (KCs), to 0, 0.5 (low), 5 (median) and 10 mM (high dose) APAP for 24 h, with/without lipopolysaccharide (LPS). Microarray-technology was used to evaluate the transcriptome changes. In the presence of LPS, the median-dose of APAP is sufficient to inhibit the expression of respiratory chain- and antioxidant-related genes, suggesting the involvement of reactive oxygen species (ROS) and oxidative stress. Furthermore, the median- and high-dose of APAP inhibited the expression of Fc fragment receptor (FcγR)-coding genes, regardless of the presence of LPS. The toll-like receptor 4 (TLR4) expression, however, was continuously elevated after the LPS/APAP co-exposures, which may result in reduced KC-phagocytosis and unbalanced cytokine patterns. Compared to the treatment with LPS only, LPS/APAP co-exposures induced the production of interleukin (IL)-8, a pro-inflammatory cytokine, but suppressed the secretion of IL-6, a cytokine regulating hepatic regeneration, along with the increase in APAP dosages. In addition to the disrupted mitochondrial functions, the presence of LPS exacerbated APAP toxicity. These findings suggest that 3D Microtissues are a suitable model for the mechanistic exploration of inflammation-associated drug toxicity.

DMSO induces drastic changes in human cellular processes and epigenetic landscape in vitro.

Though clinical trials for medical applications of dimethyl sulfoxide (DMSO) reported toxicity in the 1960s, later, the FDA classified DMSO in the safest solvent category. DMSO became widely used in many biomedical fields and biological effects were overlooked. Meanwhile, biomedical science has evolved towards sensitive high-throughput techniques and new research areas, including epigenomics and microRNAs. Considering its wide use, especially for cryopreservation and in vitro assays, we evaluated biological effect of DMSO using these technological innovations. We exposed 3D cardiac and hepatic microtissues to medium with or without 0.1% DMSO and analyzed the transcriptome, proteome and DNA methylation profiles. In both tissue types, transcriptome analysis detected >2000 differentially expressed genes affecting similar biological processes, thereby indicating consistent cross-organ actions of DMSO. Furthermore, microRNA analysis revealed large-scale deregulations of cardiac microRNAs and smaller, though still massive, effects in hepatic microtissues. Genome-wide methylation patterns also revealed tissue-specificity. While hepatic microtissues demonstrated non-significant changes, findings from cardiac microtissues suggested disruption of DNA methylation mechanisms leading to genome-wide changes. The extreme changes in microRNAs and alterations in the epigenetic landscape indicate that DMSO is not inert. Its use should be reconsidered, especially for cryopreservation of embryos and oocytes, since it may impact embryonic development.

The idiosyncratic drug-induced gene expression changes in HepG2 cells.

The inflammatory stress has been associated with an increase in susceptibility to idiosyncratic drug-induced liver injury (DILI). However, the molecular mechanisms of this inflammation-associated idiosyncratic drug hepatotoxicity remain unknown. We exposed HepG2 cells with high and low doses of three idiosyncratic (I) and three non-idiosyncratic (N) compounds, in the presence (I+ and N+) or absence (I- and N-) of a cytokine mix for 6, 12 and 24 h. To investigate the genome-wide expression patterns, microarray was performed using the Agilent 4×44K Whole Human Genome chips. The data presented in this DIB include the expression of genes participating in the ceramide metabolism, ER stress, apoptosis and cell survival pathways. The functions of these genes were illustrated in our associated article (Jiang et al., 2017) [1]. Raw and normalized gene expression data are available through NCBI GEO (accession number GSE102006).

Inferring transcription factor activity from microarray data reveals novel targets for toxicological investigations.

Transcription factors (TFs) are important modulators of the inducible portion of the transcriptome, and therefore relevant in the context of exposure to exogenous compounds. Current approaches to predict the activity of TFs in biological systems are usually restricted to a few entities at a time due to low-throughput techniques targeting a limited fraction of annotated human TFs. Therefore, high-throughput alternatives may help to identify new targets of mechanistic and predictive value in toxicological investigations. In this study, we inferred the activity multiple TFs using publicly available microarray data from primary human hepatocytes exposed to hundreds of chemicals and evaluated these molecular profiles using multiple correspondence analysis. Our results demonstrate that the lowest dose and latest exposure time (24h) in a subset of chemicals generates a signature indicative of carcinogenicity possibly due to DNA-damaging properties. Furthermore, profiles from the earliest exposure time (2h) and highest dose creates clusters of chemicals implicated in the development of diverse forms of drug-induced liver injury (DILI). Both approaches yielded a number of TFs with similar activity across groups of chemicals, including TFs known in toxicological responses such as AhR, NFE2L2 (Nrf2), NF-κB and PPARG. FOXM1, IRF1 and E2F4 were some of the TFs identified that may be relevant in genotoxic carcinogenesis. SMADs (SMAD1, SMAD2, SMAD5) and KLF5 were identified as some of potentially new TFs whose inferred activities were linked to acute and progressive outcomes in DILI. In conclusion this study offers a novel mechanistic approach targeting TF activity during chemical exposure.

Omics-based identification of the combined effects of idiosyncratic drugs and inflammatory cytokines on the development of drug-induced liver injury.

The mechanisms of idiosyncratic drug-induced hepatotoxicity remain largely unclear. It has demonstrated that the drug idiosyncrasy is potentiated in the context of inflammation and intracellular ceramides may play a role in this process. To study the mechanisms, HepG2 cells were co-treated with high and low doses of three idiosyncratic (I) and three non-idiosyncratic (N) compounds, with (I+ and N+) or without (I- and N-) a cytokine mix. Microarray, lipidomics and flow cytometry were performed to investigate the genome-wide expression patterns, the intracellular ceramide levels and the induction of apoptosis. We found that all I+ treatments significantly influenced the immune response- and response to stimulus-associated gene ontology (GO) terms, but the induction of apoptotic pathways, which was confirmed by flow cytometry, only appeared to be induced after the high-dose treatment. The ceramide signaling-, ER stress-, NF-kB activation- and mitochondrial activity-related pathways were biologically involved in apoptosis induced by the high-dose I+. Additionally, genes participating in ceramide metabolism were significantly altered resulting in a measurable increase in ceramide levels. The increases in ceramide concentrations may induce ER stress and activate the JNK pathway by affecting the expression of the related genes, and eventually trigger the mitochondria-independent apoptosis in hepatocytes. Overall, our study provides a potential mechanism to explain the role of inflammation in idiosyncratic drug reactions.

The exposome in practice: Design of the EXPOsOMICS project.

EXPOsOMICS is a European Union funded project that aims to develop a novel approach to the assessment of exposure to high priority environmental pollutants, by characterizing the external and the internal components of the exposome. It focuses on air and water contaminants during critical periods of life. To this end, the project centres on 1) exposure assessment at the personal and population levels within existing European short and long-term population studies, exploiting available tools and methods which have been developed for personal exposure monitoring (PEM); and 2) multiple "omic" technologies for the analysis of biological samples (internal markers of external exposures). The search for the relationships between external exposures and global profiles of molecular features in the same individuals constitutes a novel advancement towards the development of "next generation exposure assessment" for environmental chemicals and their mixtures. The linkage with disease risks opens the way to what are defined here as 'exposome-wide association studies' (EWAS).

ToxDB: pathway-level interpretation of drug-treatment data.

MOTIVATION: Extensive drug treatment gene expression data have been generated in order to identify biomarkers that are predictive for toxicity or to classify compounds. However, such patterns are often highly variable across compounds and lack robustness. We and others have previously shown that supervised expression patterns based on pathway concepts rather than unsupervised patterns are more robust and can be used to assess toxicity for entire classes of drugs more reliably. RESULTS: We have developed a database, ToxDB, for the analysis of the functional consequences of drug treatment at the pathway level. We have collected 2694 pathway concepts and computed numerical response scores of these pathways for 437 drugs and chemicals and 7464 different experimental conditions. ToxDB provides functionalities for exploring these pathway responses by offering tools for visualization and differential analysis allowing for comparisons of different treatment parameters and for linking this data with toxicity annotation and chemical information.Database URL:http://toxdb.molgen.mpg.de.

Inter-laboratory study of human in vitro toxicogenomics-based tests as alternative methods for evaluating chemical carcinogenicity: a bioinformatics perspective.

The assessment of the carcinogenic potential of chemicals with alternative, human-based in vitro systems has become a major goal of toxicogenomics. The central read-out of these assays is the transcriptome, and while many studies exist that explored the gene expression responses of such systems, reports on robustness and reproducibility, when testing them independently in different laboratories, are still uncommon. Furthermore, there is limited knowledge about variability induced by the data analysis protocols. We have conducted an inter-laboratory study for testing chemical carcinogenicity evaluating two human in vitro assays: hepatoma-derived cells and hTERT-immortalized renal proximal tubule epithelial cells, representing liver and kidney as major target organs. Cellular systems were initially challenged with thirty compounds, genome-wide gene expression was measured with microarrays, and hazard classifiers were built from this training set. Subsequently, each system was independently established in three different laboratories, and gene expression measurements were conducted using anonymized compounds. Data analysis was performed independently by two separate groups applying different protocols for the assessment of inter-laboratory reproducibility and for the prediction of carcinogenic hazard. As a result, both workflows came to very similar conclusions with respect to (1) identification of experimental outliers, (2) overall assessment of robustness and inter-laboratory reproducibility and (3) re-classification of the unknown compounds to the respective toxicity classes. In summary, the developed bioinformatics workflows deliver accurate measures for inter-laboratory comparison studies, and the study can be used as guidance for validation of future carcinogenicity assays in order to implement testing of human in vitro alternatives to animal testing.

Cell line-specific oxidative stress in cellular toxicity: A toxicogenomics-based comparison between liver and colon cell models.

Imbalance between high reactive oxygen species formation and antioxidant capacity in the colon and liver has been linked to increased cancer risk. However, knowledge about possible cell line-specific oxidative stress-mechanisms is limited. To explore this further, gene expression data from a human liver and colon cell line (HepG2/Caco-2), both exposed to menadione and H2O2 at six time points (0.5-1-2-4-8 and 24h) were compared in association with cell cycle distribution. In total, 3164 unique- and 1827 common genes were identified between HepG2 and Caco-2 cells. Despite the higher number of unique genes, most oxidative stress-related genes such as CAT, OGG1, NRF2, NF-κB, GCLC, HMOX1 and GSR were differentially expressed in both cell lines. However, cell-specific regulation of genes such as KEAP1 and GCLM, or of the EMT pathway, which are of pathophysiological importance, indicates that oxidative stress induces different transcriptional effects and outcomes in the two selected cell lines. In addition, expression levels and/or -direction of common genes were often different in HepG2 and Caco-2 cells, and this led to very diverse downstream effects as confirmed by correlating pathways to cell cycle changes. Altogether, this work contributes to obtaining a better molecular understanding of cell line-specific toxicity upon exposure to oxidative stress-inducing compounds.

Prediagnostic transcriptomic markers of Chronic lymphocytic leukemia reveal perturbations 10 years before diagnosis.

BACKGROUND: B-cell lymphomas are a diverse group of hematological neoplasms with differential etiology and clinical trajectories. Increased insights in the etiology and the discovery of prediagnostic markers have the potential to improve the clinical course of these neoplasms. METHODS: We investigated in a prospective study global gene expression in peripheral blood mononuclear cells of 263 incident B-cell lymphoma cases, diagnosed between 1 and 17 years after blood sample collection, and 439 controls, nested within two European cohorts. RESULTS: Our analyses identified only transcriptomic markers for specific lymphoma subtypes; few markers of multiple myeloma (N = 3), and 745 differentially expressed genes in relation to future risk of chronic lymphocytic leukemia (CLL). The strongest of these associations were consistently found in both cohorts and were related to (B-) cell signaling networks and immune system regulation pathways. CLL markers exhibited very high predictive abilities of disease onset even in cases diagnosed more than 10 years after blood collection. CONCLUSIONS: This is the first investigation on blood cell global gene expression and future risk of B-cell lymphomas. We mainly identified genes in relation to future risk of CLL that are involved in biological pathways, which appear to be mechanistically involved in CLL pathogenesis. Many but not all of the top hits we identified have been reported previously in studies based on tumor tissues, therefore suggesting that a mixture of preclinical and early disease markers can be detected several years before CLL clinical diagnosis.

Increased frequency of micronuclei in mononucleated lymphocytes and cytome analysis in healthy newborns as an early warning biomarkers of possible future health risks.

Impact of intrauterine development on health risks during adolescence and adulthood still needs to be investigated. The aim of study was to compare genome damage in newborns and mothers using the cytokinesis blocked micronucleus assay, nuclear division index (NDI), and centromere fluorescence in situ hybridization. The study was performed on 92 mothers and their newborns. Results showed that micronucleus frequency in binuclear T-lymphocytes (MNBN) in newborns was significantly lower than in mothers but higher in mononuclear T-lymphocytes (MNMONO). The NDI in the mothers was significantly higher than in the newborns. In newborns with <2500g birth weight, NDI was similar to the mothers'. Mothers have significantly more centromere negative micronuclei than newborns. A significantly higher NDI and MNBN was found in newborns with ≥2 MNMONO/1000 than in newborns with <2 MNMONO/1000. It is suggested that MOMN and NDI might be good candidates for biomarkers of health risks in newborns.

Baseline and genotoxic compound induced gene expression profiles in HepG2 and HepaRG compared to primary human hepatocytes.

Efforts are put into developing toxicogenomics-based toxicity testing methods using in vitro human cell models for improving human risk assessment/replacing animal models. Human in vitro liver models include HepG2, HepaRG and primary human hepatocytes (PHH). Studies on comparability/applicability of these cell types mainly focus on assessing baseline biotransformation capacities/cytochrome P450-inducibility, but compound-induced gene expression profiles are at least as important. Therefore, we compared baseline and aflatoxin B1- and benzo(α)pyrene-induced gene expression profiles in HepG2, HepaRG and PHH (11-13 donors). At baseline, all liver models differ from each other with respect to whole genome gene expression levels. PHH show profound inter-individual differences, and are most similar to HepaRG. After compound exposure, induced gene expression profiles are more similar between cell models, especially for benzo(α)pyrene. Pathways involved in compound metabolism are induced in all 3 models, while others are more pronounced in a specific cell model. Examples are transcriptomic modifications of carbohydrate-related genes (HepaRG) and of receptor-related genes (PHH) after benzo(α)pyrene exposure, and of cell cycle-related genes (HepG2) after aflatoxin B1 exposure. PHH gene expression responses are the most heterogeneous. In conclusion, at base line level PHH are more similar to HepaRG than to HepG2, but for toxicogenomics applications both cell lines perform equally well in comparison to PHH.

Time series analysis of oxidative stress response patterns in HepG2: a toxicogenomics approach.

Oxidative stress plays an important role in chemically induced liver injury, however, our insight into molecular responses to different oxygen radicals is fragmentary. Since these cellular responses will differ over time, examining time-dependent changes in gene expression, and correlating these with markers for oxidative stress, may provide new insights into responses to oxidants. We used the human hepatoma cell line HepG2 to investigate the effects of oxidative stress on the transcriptome level by micro-arrays at seven time points (0.5, 1, 2, 4, 6, 8 and 24h) following exposure to the oxidants menadione, hydrogen peroxide and tert-butyl hydroperoxide including the effects on cell cycle and apoptosis by flow cytometry, protein carbonyl formation by spectrophotometry and oxidative DNA damage by FPG-comet. In total, 3429 genes were differentially expressed, including 136 genes that were significantly modified by all oxidants. Time-dependent biological pathway analysis showed that these genes were particularly involved in inflammatory responses, cell cycle processes and glutathione signaling. These responses were confirmed and supported by phenotypic anchoring to the different cellular endpoints. In addition, using an innovative temporal analysis we established an oxidative stress-related gene expression time cluster. Altogether, this study provides new insights in temporal oxidative stress mechanisms and demonstrates sequential cellular responses that may contribute to a better hazard identification and the mechanisms of toxicological responses in the liver induced by oxidative stress.

Benzo[a]pyrene-induced transcriptomic responses in primary hepatocytes and in vivo liver: toxicokinetics is essential for in vivo-in vitro comparisons.

The traditional 2-year cancer bioassay needs replacement by more cost-effective and predictive tests. The use of toxicogenomics in an in vitro system may provide a more high-throughput method to investigate early alterations induced by carcinogens. Recently, the differential gene expression response in wild-type and cancer-prone Xpa (-/-) p53 (+/-) primary mouse hepatocytes after exposure to benzo[a]pyrene (B[a]P) revealed downregulation of cancer-related pathways in Xpa (-/-) p53 (+/-) hepatocytes only. Here, we investigated pathway regulation upon in vivo B[a]P exposure of wild-type and Xpa (-/-) p53 (+/-) mice. In vivo transcriptomics analysis revealed a limited gene expression response in mouse livers, but with a significant induction of DNA replication and apoptotic/anti-apoptotic cellular responses in Xpa (-/-) p53 (+/-) livers only. In order to be able to make a meaningful in vivo-in vitro comparison we estimated internal in vivo B[a]P concentrations using DNA adduct levels and physiologically based kinetic modeling. Based on these results, the in vitro concentration that corresponded best with the internal in vivo dose was chosen. Comparison of in vivo and in vitro data demonstrated similarities in transcriptomics response: xenobiotic metabolism, lipid metabolism and oxidative stress. However, we were unable to detect cancer-related pathways in either wild-type or Xpa (-/-) p53 (+/-) exposed livers, which were previously found to be induced by B[a]P in Xpa (-/-) p53 (+/-) primary hepatocytes. In conclusion, we showed parallels in gene expression responses between livers and primary hepatocytes upon exposure to equivalent concentrations of B[a]P. Furthermore, we recommend considering toxicokinetics when modeling a complex in vivo endpoint with in vitro models.

Toxicogenomic profiles in relation to maternal immunotoxic exposure and immune functionality in newborns.

A crucial period for the development of the immune system occurs in utero. This results in a high fetal vulnerability to immunotoxic exposure, and indeed, immunotoxic effects have been reported, demonstrating negative effects on immune-related health outcomes and immune functionality. Within the NewGeneris cohort BraMat, a subcohort of the Norwegian Mother and Child Cohort Study (MoBa), immunotoxicity was demonstrated for polychlorinated biphenyls and dioxins, showing associations between estimated maternal intake levels and reduced measles vaccination responses in the offspring at the age of 3. The present study aimed to investigate this link at the transcriptomic level within the same BraMat cohort. To this end, whole-genome gene expression in cord blood was investigated and found to be associated with maternal Food Frequency Questionnaires-derived exposure estimates and with vaccination responses in children at 3 years of age. Because the literature reports gender specificity in the innate, humoral, and cell-mediated responses to viral vaccines, separate analysis for males and females was conducted. Separate gene sets for male and female neonates were identified, comprising genes significantly correlating with both 2,3,7,8-tetrachlorodibenzodioxin (TCDD) and polychlorinated biphenyls (PCB) exposure and with measles vaccination response. Noteworthy, genes correlating negatively with exposure in general show positive correlations with antibody levels and vice versa. For both sexes, these included immune-related genes, suggesting immunosuppressive effects of maternal exposure to TCDD and PCB at the transcriptomic level in neonates in relation to measles vaccination response 3 years later.

Compound-specific effects of diverse neurodevelopmental toxicants on global gene expression in the neural embryonic stem cell test (ESTn).

Alternative assays for developmental toxicity testing are needed to reduce animal use in regulatory toxicology. The in vitro murine neural embryonic stem cell test (ESTn) was designed as an alternative for neurodevelopmental toxicity testing. The integration of toxicogenomic-based approaches may further increase predictivity as well as provide insight into underlying mechanisms of developmental toxicity. In the present study, we investigated concentration-dependent effects of six mechanistically diverse compounds, acetaldehyde (ACE), carbamazepine (CBZ), flusilazole (FLU), monoethylhexyl phthalate (MEHP), penicillin G (PENG) and phenytoin (PHE), on the transcriptome and neural differentiation in the ESTn. All compounds with the exception of PENG altered ESTn morphology (cytotoxicity and neural differentiation) in a concentration-dependent manner. Compound induced gene expression changes and corresponding enriched gene ontology biological processes (GO-BP) were identified after 24h exposure at equipotent differentiation-inhibiting concentrations of the compounds. Both compound-specific and common gene expression changes were observed between subsets of tested compounds, in terms of significance, magnitude of regulation and functionality. For example, ACE, CBZ and FLU induced robust changes in number of significantly altered genes (≥ 687 genes) as well as a variety of GO-BP, as compared to MEHP, PHE and PENG (≤ 55 genes with no significant changes in GO-BP observed). Genes associated with developmentally related processes (embryonic morphogenesis, neuron differentiation, and Wnt signaling) showed diverse regulation after exposure to ACE, CBZ and FLU. In addition, gene expression and GO-BP enrichment showed concentration dependence, allowing discrimination of non-toxic versus toxic concentrations on the basis of transcriptomics. This information may be used to define adaptive versus toxic responses at the transcriptome level.