Search

Benzo(a)pyrene (BaP) is a ubiquitous carcinogen resulting from incomplete combustion of organic compounds and also present at high levels in cigarette smoke. A wide range of biological effects has been attributed to BaP and its genotoxic metabolite BPDE, but the contribution to BaP toxicity of intermediary metabolites generated along the detoxification path remains unknown. Here, we report for the first time how 3-OH-BaP, 9,10-diol and BPDE, three major BaP metabolites, temporally relate to BaP-induced transcriptomic alterations in HepG2 cells. Since BaP is also known to induce AhR activation, we additionally evaluated TCDD to source the expression of non-genotoxic AhR-mediated patterns. 9,10-Diol was shown to activate several transcription factor networks related to BaP metabolism (AhR), oxidative stress (Nrf2) and cell proliferation (HIF-1α, AP-1) in particular at early time points, while BPDE influenced expression of genes involved in cell energetics, DNA repair and apoptotic pathways. Also, in order to grasp the role of BaP and its metabolites in chemical hepatocarcinogenesis, we compared expression patterns from BaP(-metabolites) and TCDD to a signature set of approximately nine thousand gene expressions derived from hepatocellular carcinoma (HCC) patients. While transcriptome modulation by TCDD appeared not significantly related to HCC, BaP and BPDE were shown to deregulate metastatic markers via non-genotoxic and genotoxic mechanisms and activate inflammatory pathways (NF-κβ signaling, cytokine-cytokine receptor interaction). BaP also showed strong repression of genes involved in cholesterol and fatty acid biosynthesis. Altogether, this study provides new insights into BaP-induced toxicity and sheds new light onto its mechanism of action as a hepatocarcinogen.

Acetaminophen (APAP) is a readily available over-the-counter drug and is one of the most commonly used analgesics/antipyretics worldwide. Large interindividual variation in susceptibility toward APAP-induced liver failure has been reported. However, the exact underlying factors causing this variability in susceptibility are still largely unknown. The aim of this study was to better understand this variability in response to APAP by evaluating interindividual differences in gene expression changes and APAP metabolite formation in primary human hepatocytes (PHH) from several donors (n = 5) exposed in vitro to a non-toxic to toxic APAP dose range. To evaluate interindividual variation, gene expression data/levels of metabolites were plotted against APAP dose/donor. The correlation in APAP dose response between donors was calculated by comparing data points from one donor to the data points of all other donors using a Pearson-based correlation analysis. From that, a correlation score/donor for each gene/metabolite was defined, representing the similarity of the omics response to APAP in PHH of a particular donor to all other donors. The top 1 % highest variable genes were selected for further evaluation using gene set overrepresentation analysis. The biological processes in which the genes with high interindividual variation in expression were involved include liver regeneration, inflammatory responses, mitochondrial stress responses, hepatocarcinogenesis, cell cycle, and drug efficacy. Additionally, the interindividual variation in the expression of these genes could be associated with the variability in expression levels of hydroxyl/methoxy-APAP and C8H13O5N-APAP-glucuronide. The before-mentioned metabolites or their derivatives have also been reported in blood of humans exposed to therapeutic APAP doses. Possibly these findings can contribute to elucidating the causative factors of interindividual susceptibility toward APAP.

Microarray-based transcriptomic analysis has been demonstrated to hold the opportunity to study the effects of human exposure to, e.g., chemical carcinogens at the whole genome level, thus yielding broad-ranging molecular information on possible carcinogenic effects. Since genes do not operate individually but rather through concerted interactions, analyzing and visualizing networks of genes should provide important mechanistic information, especially upon connecting them to functional parameters, such as those derived from measurements of biomarkers for exposure and carcinogenic risk. Conventional methods such as hierarchical clustering and correlation analyses are frequently used to address these complex interactions but are limited as they do not provide directional causal dependence relationships. Therefore, our aim was to apply Bayesian network inference with the purpose of phenotypic anchoring of modified gene expressions. We investigated a use case on transcriptomic responses to cigarette smoking in humans, in association with plasma cotinine levels as biomarkers of exposure and aromatic DNA-adducts in blood cells as biomarkers of carcinogenic risk. Many of the genes that appear in the Bayesian networks surrounding plasma cotinine, and to a lesser extent around aromatic DNA-adducts, hold biologically relevant functions in inducing severe adverse effects of smoking. In conclusion, this study shows that Bayesian network inference enables unbiased phenotypic anchoring of transcriptomics responses. Furthermore, in all inferred Bayesian networks several dependencies are found which point to known but also to new relationships between the expression of specific genes, cigarette smoke exposure, DNA damaging-effects, and smoking-related diseases, in particular associated with apoptosis, DNA repair, and tumor suppression, as well as with autoimmunity.

In order to improve attrition rates of candidate-drugs there is a need for a better understanding of the mechanisms underlying drug-induced hepatotoxicity. We aim to further unravel the toxicological response of hepatocytes to a prototypical cholestatic compound by integrating transcriptomic and metabonomic profiling of HepG2 cells exposed to Cyclosporin A. Cyclosporin A exposure induced intracellular cholesterol accumulation and diminished intracellular bile acid levels. Performing pathway analyses of significant mRNAs and metabolites separately and integrated, resulted in more relevant pathways for the latter. Integrated analyses showed pathways involved in cell cycle and cellular metabolism to be significantly changed. Moreover, pathways involved in protein processing of the endoplasmic reticulum, bile acid biosynthesis and cholesterol metabolism were significantly affected. Our findings indicate that an integrated approach combining metabonomics and transcriptomics data derived from representative in vitro models, with bioinformatics can improve our understanding of the mechanisms of action underlying drug-induced hepatotoxicity. Furthermore, we showed that integrating multiple omics and thereby analyzing genes, microRNAs and metabolites of the opposed model for drug-induced cholestasis can give valuable information about mechanisms of drug-induced cholestasis in vitro and therefore could be used in toxicity screening of new drug candidates at an early stage of drug discovery., (Copyright © 2015 Elsevier Ltd. All rights reserved.)

Large differences in toxicity responses occur within the human population. In this study we evaluate whether interindividual variation in baseline enzyme activity (EA)/gene expression (GE) levels in liver predispose for the variation in toxicity responses by assessing dose-response relationships for several prototypical hepatotoxicants. Baseline levels of cytochrome-P450 (CYP) GE/EA were measured in precision-cut human liver slices. Slices (n=4-5/compound) were exposed to a dose-range of acetaminophen, aflatoxin B1, benzo(α) pyrene or 2-nitrofluorene. Interindividual variation in induced genotoxicity (COMET-assay and CDKN1A/p21 GE) and cytotoxicity (lactate dehydrogenase-leakage), combined with NQO1- and GSTM1-induced GE-responses for oxidative stress and GE-responses of several CYPs was evaluated. The benchmark dose-approach was applied as a tool to model exposure responses on an individual level. Variation in baseline CYP levels, both GE and EA, can explain variation in compound exposure-responses on an individual level. Network analyses enable the definition of key parameters influencing interindividual variation after compound exposure. For 2-nitrofluorene, this analysis suggests involvement of CYP1B1 in the metabolism of this compound, which represents a novel finding. In this study, GSTM1 which is known to be highly polymorphic within the human population, but so far could not be linked to toxicity in acetaminophen-poisoned patients, is suggested to cause interindividual variability in acetaminophen-metabolism, dependent on the individual's gene expression-responses of CYP-enzymes. This study demonstrates that using interindividual variation within network modelling provides a source for the definition of essential and even new parameters involved in compound-related metabolism. This information might enable ways to make more quantitative estimates of human risks., (Copyright © 2014 Elsevier Ireland Ltd. All rights reserved.)

The concordance of RNA-sequencing (RNA-seq) with microarrays for genome-wide analysis of differential gene expression has not been rigorously assessed using a range of chemical treatment conditions. Here we use a comprehensive study design to generate Illumina RNA-seq and Affymetrix microarray data from the same liver samples of rats exposed in triplicate to varying degrees of perturbation by 27 chemicals representing multiple modes of action (MOAs). The cross-platform concordance in terms of differentially expressed genes (DEGs) or enriched pathways is linearly correlated with treatment effect size (R(2)0.8). Furthermore, the concordance is also affected by transcript abundance and biological complexity of the MOA. RNA-seq outperforms microarray (93% versus 75%) in DEG verification as assessed by quantitative PCR, with the gain mainly due to its improved accuracy for low-abundance transcripts. Nonetheless, classifiers to predict MOAs perform similarly when developed using data from either platform. Therefore, the endpoint studied and its biological complexity, transcript abundance and the genomic application are important factors in transcriptomic research and for clinical and regulatory decision making.

The EU FP6 project carcinoGENOMICS explored the combination of toxicogenomics and in vitro cell culture models for identifying organotypical genotoxic- and non-genotoxic carcinogen-specific gene signatures. Here the performance of its gene classifier, derived from exposure of metabolically competent human HepaRG cells to prototypical non-carcinogens (10 compounds) and hepatocarcinogens (20 compounds), is reported. Analysis of the data at the gene and the pathway level by using independent biostatistical approaches showed a distinct separation of genotoxic from non-genotoxic hepatocarcinogens and non-carcinogens (up to 88 % correct prediction). The most characteristic pathway responding to genotoxic exposure was DNA damage. Interlaboratory reproducibility was assessed by blindly testing of three compounds, from the set of 30 compounds, by three independent laboratories. Subsequent classification of these compounds resulted in correct prediction of the genotoxicants. As expected, results on the non-genotoxic carcinogens and the non-carcinogens were less predictive. In conclusion, the combination of transcriptomics with the HepaRG in vitro cell model provides a potential weight of evidence approach for the evaluation of the genotoxic potential of chemical substances.

With the number of new drug candidates increasing every year, there is a need for high-throughput human toxicity screenings. As the liver is the most important organ in drug metabolism and thus capable of generating relatively high levels of toxic metabolites, it is important to find a reliable strategy to screen for drug-induced hepatotoxicity. Microarray-based transcriptomics is a well-established technique in toxicogenomics research and is an ideal approach to screen for drug-induced injury at an early stage. The aim of this study was to prove the principle of classifying known hepatotoxicants and nonhepatotoxicants using their distinctive gene expression profiles in vitro in HepG2 cells. Furthermore, we undertook to subclassify the hepatotoxic compounds by investigating the subclass of cholestatic compounds. Prediction analysis for microarrays was used for classification of hepatotoxicants and nonhepatotoxicants, which resulted in an accuracy of 92% on the training set and 91% on the validation set, using 36 genes. A second model was set up with the goal of finding classifiers for cholestasis, resulting in 12 genes that appeared capable of correctly classifying 8 of the 9 cholestatic compounds, resulting in an accuracy of 93%. We were able to prove the principle that transcriptomic analyses of HepG2 cells can indeed be used to classify chemical entities for hepatotoxicity. Genes selected for classification of hepatotoxicity and cholestasis indicate that endoplasmic reticulum stress and the unfolded protein response may be important cellular effects of drug-induced liver injury. However, the number of compounds in both the training set and the validation set should be increased to improve the reliability of the prediction.

Background: Leukemia incidence has increased in recent decades among European children, suggesting that early-life environmental exposures play an important role in disease development., Objectives: We investigated the hypothesis that childhood susceptibility may increase as a result of in utero exposure to carcinogens and hormonally acting factors. Using cord blood samples from the NewGeneris cohort, we examined associations between a range of biomarkers of carcinogen exposure and hormonally acting factors with micronuclei (MN) frequency as a proxy measure of cancer risk. Associations with gene expression and genotype were also explored., Methods: DNA and protein adducts, gene expression profiles, circulating hormonally acting factors, and GWAS (genome-wide association study) data were investigated in relation to genomic damage measured by MN frequency in lymphocytes from 623 newborns enrolled between 2006 and 2010 across Europe., Results: Malondialdehyde DNA adducts (M1dG) were associated with increased MN frequency in binucleated lymphocytes (MNBN), and exposure to androgenic, estrogenic, and dioxin-like compounds was associated with MN frequency in mononucleated lymphocytes (MNMONO), although no monotonic exposure-outcome relationship was observed. Lower frequencies of MNBN were associated with a 1-unit increase expression of PDCD11, LATS2, TRIM13, CD28, SMC1A, IL7R, and NIPBL genes. Gene expression was significantly higher in association with the highest versus lowest category of bulky and M1dG-DNA adducts for five and six genes, respectively. Gene expression levels were significantly lower for 11 genes in association with the highest versus lowest category of plasma AR CALUX® (chemically activated luciferase expression for androgens) (8 genes), ERα CALUX® (for estrogens) (2 genes), and DR CALUX® (for dioxins). Several SNPs (single-nucleotide polymorphisms) on chromosome 11 near FOLH1 significantly modified associations between androgen activity and MNBN frequency. Polymorphisms in EPHX1/2 and CYP2E1 were associated with MNBN., Conclusion: We measured in utero exposure to selected environmental carcinogens and circulating hormonally acting factors and detected associations with MN frequency in newborns circulating T lymphocytes. The results highlight mechanisms that may contribute to carcinogen-induced leukemia and require further research.

Background: Certain multi-walled carbon nanotubes (MWCNTs) have been shown to elicit asbestos-like toxicological effects. To reduce needs for risk assessment it has been suggested that the physicochemical characteristics or reactivity of nanomaterials could be used to predict their hazard. Fibre-shape and ability to generate reactive oxygen species (ROS) are important indicators of high hazard materials. Asbestos is a known ROS generator, while MWCNTs may either produce or scavenge ROS. However, certain biomolecules, such as albumin - used as dispersants in nanomaterial preparation for toxicological testing in vivo and in vitro - may reduce the surface reactivity of nanomaterials., Methods: Here, we investigated the effect of bovine serum albumin (BSA) and cell culture medium with and without BEAS 2B cells on radical formation/scavenging by five MWCNTs, Printex 90 carbon black, crocidolite asbestos, and glass wool, using electron spin resonance (ESR) spectroscopy and linked this to cytotoxic effects measured by trypan blue exclusion assay. In addition, the materials were characterized in the exposure medium (e.g. for hydrodynamic size-distribution and sedimentation rate)., Results: The test materials induced highly variable cytotoxic effects which could generally be related to the abundance and characteristics of agglomerates/aggregates and to the rate of sedimentation. All carbon nanomaterials were found to scavenge hydroxyl radicals (•OH) in at least one of the solutions tested. The effect of BSA was different among the materials. Two types of long, needle-like MWCNTs (average diameter >74 and 64.2 nm, average length 5.7 and 4.0 μm, respectively) induced, in addition to a scavenging effect, a dose-dependent formation of a unique, yet unidentified radical in both absence and presence of cells, which also coincided with cytotoxicity., Conclusions: Culture medium and BSA affects scavenging/production of •OH by MWCNTs, Printex 90 carbon black, asbestos and glass-wool. An unidentified radical is generated by two long, needle-like MWCNTs and these two CNTs were more cytotoxic than the other CNTs tested, suggesting that this radical could be related to the adverse effects of MWCNTs.

The toxic mechanisms of cisplatin have been frequently studied in many species and in vitro cell models. The Netherlands Toxicogenomics Centre focuses on developing in vitro alternatives using genomics technologies for animal-based assays on, e.g. genotoxic hazards. Models such as human hepatocellular carcinoma cell line (HepG2) cells, mouse primary hepatocytes (PMH) and mouse embryonic stem cells (mESC) are used. Our aim was to identify possibly robust conserved mechanisms between these models using cisplatin as model genotoxic agent. Transcriptomic data newly generated from HepG2 cells and PMH exposed to 7 µM cisplatin for 12, 24 and 48h and 24 and 48h, respectively, were compared with published data from mESC exposed to 5 µM cisplatin for 2-24h. Due to differences in response time between models and marginal changes after shorter exposure periods, we focused on 24 and 48h. At gene level, 44 conserved differentially expressed genes (DEG), involved in processes such as apoptosis, cell cycle, DNA damage response and DNA repair, were found. Functional analysis shows that limited numbers of pathways are conserved. Transcription factor (TF) network analysis indicates 12 common TF networks responding among all models and time points. Four TF, HNF4-α, SP1, c-MYC and p53, capable of regulating ±50% of all DEG, seem of equal importance in all models and exposure periods. Here we showed that transcriptomic responses across several in vitro cell models following exposure to cisplatin are mainly determined by a conserved complex network of 4 TFs. These conserved responses are hypothesised to provide most relevant information for human toxicity prediction and may form the basis for new in vitro alternatives of risk assessment.

Current testing models for predicting drug-induced liver injury are inadequate, as they basically under-report human health risks. We present here an approach towards developing pathways based on hepatotoxicity-associated gene groups derived from two types of publicly accessible hepatotoxicity databases, in order to develop drug-induced liver injury biomarker profiles. One human liver 'omics-based and four text-mining-based databases were explored for hepatotoxicity-associated gene lists. Over-representation analysis of these gene lists with a hepatotoxicant-exposed primary human hepatocytes data set showed that human liver 'omics gene lists performed better than text-mining gene lists and the results of the latter differed strongly between databases. However, both types of databases contained gene lists demonstrating biomarker potential. Visualizing those in pathway format may aid in interpreting the biomolecular background. We conclude that exploiting existing and openly accessible databases in a dedicated manner seems promising in providing venues for translational research in toxicology and biomarker development.

The mycotoxin Ochratoxin A (OTA) is a potent renal carcinogen in male rats. Transcriptomic studies on OTA (4 in vitro, 6 in vivo, 2 in vitro/in vivo) have been reviewed. The aim of 6 of them was mainly mechanistic whereas the rest had mostly predictive (1) or evaluation (5) purposes. An overall tendency towards gene expression downregulation was observed, probably as a result of protein synthesis inhibition. DNA damage response genes were not deregulated in most of the studies. Genes involved in acute renal injury, cell survival and cell proliferation were upregulated in several in vivo studies. Apoptosis genes were deregulated in vitro but less affected in vivo; activation of several MAPKs has been observed. Many genes related to oxidative stress or involved in cell-to-cell interaction pathways (Wnt) or cytoskeleton structure appeared to be deregulated either in vitro or in vivo. Regucalcin was highly downregulated in vivo and other calcium homeostasis genes were significantly deregulated in vitro. Genes related to OTA transport (OATs) and metabolism (CYPs) appeared downregulated in vivo. Overall, the mechanism of action of OTA remains unclear, however transcriptomic data have contributed to new mechanistic hypothesis generation and to in vitro-in vivo comparison., (Copyright © 2013. Published by Elsevier Ltd.)

In vitro models for hepatotoxicity testing are a necessity for advancement of toxicological research. Assessing the in vitro response requires in vivo validated gene sets reflective of the hepatotoxic phenotype. Cholestasis, the impairment of bile flow, is induced in C57BL/6J mice treated with cyclosporine A (CsA) to identify phenotype reflective gene sets. CsA treatment through oral gavage for 25 days induced cholestasis, as confirmed by histopathology and serum chemistry. Over 1, 4, and 11 days of CsA exposure gradual increases in serum markers were correlated to gene expression. This phenotype-directed analysis identified gene sets specific to the onset and progression of cholestasis, such as PPAR related processes and drug metabolism, by circumventing other effects of CsA, such as immunosuppression, found in dose*time group analysis. In vivo gene sets are enriched in publicly available data sets of CsA-treated HepaRG and primary mouse hepatocytes. However, genes identified within these gene sets did not overlap between in vivo and in vitro. In vitro regulated genes represent the initial response to cholestasis, whereas in vivo genes represent the later adaptive response. We conclude that the applicability of in vitro models for hepatotoxicity testing fully depends on a solid in vivo phenotype anchored analysis., (Copyright © 2013 Elsevier Ireland Ltd. All rights reserved.)

Flanders, Belgium, is one of the most densely populated areas in Europe. The Flemish Environment and Health Survey (2002-2006) aimed at determining exposure to pollutants of neonates, adolescents, and older adults and to assess associated biological and health effects. This study investigated genome wide gene expression changes associated with a range of environmental pollutants, including cadmium, lead, PCBs, dioxin, hexachlorobenzene, p,p'-DDE, benzene, and PAHs. Gene expression levels were measured in peripheral blood cells of 20 adults with relatively high and 20 adults with relatively low combined internal exposure levels, all non-smokers aged 50-65. Pearson correlation was used to analyze associations between pollutants and gene expression levels, separately for both genders. Pollutant- and gender-specific correlation analysis results were obtained. For organochlorine pollutants, analysis within genders revealed that genes were predominantly regulated in opposite directions in males and females. Significantly modulated pathways were found to be associated with each of the exposure biomarkers measured. Pathways and/or genes related to estrogen and STAT5 signaling were correlated to organochlorine exposures in both genders. Our work demonstrates that gene expression in peripheral blood is influenced by environmental pollutants. In particular, gender-specific changes are associated with organochlorine pollutants, including gender-specific modulation of endocrine related pathways and genes. These pathways and genes have previously been linked to endocrine disruption related disorders, which in turn have been associated with organochlorine exposure. Based on our results, we recommend that males and females be considered separately when analyzing gene expression changes associated with exposures that may include chemicals with endocrine disrupting properties., (Copyright © 2013 Wiley Periodicals, Inc.)

As the conventional approach to assess the potential of a chemical to cause cancer in humans still includes the 2-year rodent carcinogenicity bioassay, development of alternative methodologies is needed. In the present study, the transcriptomics responses following exposure to genotoxic (GTX) and non-genotoxic (NGTX) hepatocarcinogens and non-carcinogens (NC) in five liver-based in vitro models, namely conventional and epigenetically stabilized cultures of primary rat hepatocytes, the human hepatoma-derived cell lines HepaRG and HepG2 and human embryonic stem cell-derived hepatocyte-like cells, are examined. For full characterization of the systems, several bioinformatics approaches are employed including gene-based, ConsensusPathDB-based and classification analysis. They provide convincingly similar outcomes, namely that upon exposure to carcinogens, the HepaRG generates a gene classifier (a gene classifier is defined as a selected set of characteristic gene signatures capable of distinguishing GTX, NGTX carcinogens and NC) able to discriminate the GTX carcinogens from the NGTX carcinogens and NC. The other in vitro models also yield cancer-relevant characteristic gene groups for the GTX exposure, but some genes are also deregulated by the NGTX carcinogens and NC. Irrespective of the tested in vitro model, the most uniformly expressed pathways following GTX exposure are the p53 and those that are subsequently induced. The NGTX carcinogens triggered no characteristic cancer-relevant gene profiles in all liver-based in vitro systems. In conclusion, liver-based in vitro models coupled with transcriptomics techniques, especially in the case when the HepaRG cell line is used, represent valuable tools for obtaining insight into the mechanism of action and identification of GTX carcinogens.

Drug-induced hepatotoxicity is a leading cause of attrition for candidate pharmaceuticals in development. New preclinical screening methods are crucial to predict drug toxicity prior to human studies. Of all in vitro hepatotoxicity models, primary human hepatocytes are considered as 'the gold standard.' However, their use is hindered by limited availability and inter-individual variation. These barriers may be overcome by using primary mouse hepatocytes. We used differential in gel electrophoresis (DIGE) to study large-scale protein expression of primary mouse hepatocytes. These hepatocytes were exposed to three well-defined hepatotoxicants: acetaminophen, amiodarone, and cyclosporin A. Each hepatotoxicant induces a different hepatotoxic phenotype. Based on the DIGE results, the mRNA expression levels of deregulated proteins from cyclosporin A-treated cells were also analyzed. We were able to distinguish cyclosporin A from controls, as well as acetaminophen and amiodarone-treated samples. Cyclosporin A induced endoplasmic reticulum (ER) stress and altered the ER-Golgi transport. Moreover, liver carboxylesterase and bile salt sulfotransferase were differentially expressed. These proteins were associated with a protective adaptive response against cyclosporin A-induced cholestasis. The results of this study are comparable with effects in HepG2 cells. Therefore, we suggest both models can be used to analyze the cholestatic properties of cyclosporin A. Furthermore, this study showed a conserved response between primary mouse hepatocytes and HepG2 cells. These findings collectively lend support for use of omics strategies in preclinical toxicology, and might inform future efforts to better link preclinical and clinical research in rational drug development.

Whole-genome transcriptome measurements are pivotal for characterizing molecular mechanisms of chemicals and predicting toxic classes, such as genotoxicity and carcinogenicity, from in vitro and in vivo assays. In recent years, deep sequencing technologies have been developed that hold the promise of measuring the transcriptome in a more complete and unbiased manner than DNA microarrays. Here, we applied this RNA-seq technology for the characterization of the transcriptomic responses in HepG2 cells upon exposure to benzo[a]pyrene (BaP), a well-known DNA damaging human carcinogen. Based on EnsEMBL genes, we demonstrate that RNA-seq detects ca 20% more genes than microarray-based technology but almost threefold more significantly differentially expressed genes. Functional enrichment analyses show that RNA-seq yields more insight into the biology and mechanisms related to the toxic effects caused by BaP, i.e., two- to fivefold more affected pathways and biological processes. Additionally, we demonstrate that RNA-seq allows detecting alternative isoform expression in many genes, including regulators of cell death and DNA repair such as TP53, BCL2 and XPA, which are relevant for genotoxic responses. Moreover, potentially novel isoforms were found, such as fragments of known transcripts, transcripts with additional exons, intron retention or exon-skipping events. The biological function(s) of these isoforms remain for the time being unknown. Finally, we demonstrate that RNA-seq enables the investigation of allele-specific gene expression, although no changes could be observed. Our results provide evidence that RNA-seq is a powerful tool for toxicology, which, compared with microarrays, is capable of generating novel and valuable information at the transcriptome level for characterizing deleterious effects caused by chemicals.

Potassium bromate (KBrO(3)) is an oxidising agent that has been widely used in the food and cosmetic industries. It has shown to be both a nephrotoxin and a renal carcinogen in in vivo and in vitro models. Here, we investigated the effects of KBrO(3) in the human and rat proximal tubular cell lines RPTEC/TERT1 and NRK-52E. A genome-wide transcriptomic screen was carried out from cells exposed to a sub-lethal concentration of KBrO(3) for 6, 24 and 72 h. Pathway analysis identified "glutathione metabolism", "Nrf2-mediated oxidative stress" and "tight junction (TJ) signalling" as the most enriched pathways. TJ signalling was less impacted in the rat model, and further studies revealed low transepithelial electrical resistance (TEER) and an absence of several TJ proteins in NRK-52E cells. In RPTEC/TERT1 cells, KBrO(3) exposure caused a decrease in TEER and resulted in altered expression of several TJ proteins. N-Acetylcysteine co-incubation prevented these effects. These results demonstrate that oxidative stress has, in conjunction with the activation of the cytoprotective Nrf2 pathway, a dramatic effect on the expression of tight junction proteins. The further understanding of the cross-talk between these two pathways could have major implications for epithelial repair, carcinogenesis and metastasis.

The γH2AX assay has recently been suggested as a new in vitro assay for detecting genotoxic (GTX) properties of chemicals. This assay is based on the phosphorylation of H2AX histone in response to DNA damage [i.e. induction of double-strand breaks (DSBs)]. Quantification of γH2AX foci using flow cytometry can rapidly detect DNA damage induced by chemicals that cause DNA DSBs. Up to now, only few compounds have been tested with this assay. The main goal of this study was to compare the performance of this automated γH2AX assay with that of standard in vitro genotoxicity assays in predicting in vivo genotoxicity. HepG2 cells were exposed to 64 selected compounds with known GTX properties and subsequently analysed for induction of γH2AX foci. The results of this assay were compared with public data from standard in vitro genotoxicity tests. Accuracy, sensitivity and specificity in predicting in vivo genotoxicity, using the γH2AX assay alone or in combinations with conventional assays, were calculated. Both the γH2AX assay and the bacterial mutagenicity test (Ames) were highly specific for in vivo GTX, whereas chromosomal aberration/micronucleus test (CA/MN) resulted in highest sensitivity. The currently widely used in vitro genotoxicity test battery-Ames test, mouse lymphoma assay (MLA) and CA/MN test-resulted in low accuracy (55-65%) to predict in vivo genotoxicity. Interestingly, the inclusion of γH2AX assay in the standard battery, instead of MLA assay, resulted in higher accuracy (62-70%) compared with other combinations. Advantage of the γH2AX assay in HepG2 cells is its high sensitivity to detect DNA-reactive GTX compounds, although the reduced sensitivity for compounds that require metabolic activation needs to be improved. In conclusion, the automated γH2AX assay can be a useful, fast and cost-effective human cell-based tool for early screening of compounds for in vivo genotoxicity.

The concept of mechanistic toxicogenomics implies that compound-induced changes in gene expression profiles provide valuable information about their mode of action. A growing number of research groups have presented evidence that whole-genome gene expression profiling techniques might be used as tools for in vivo and in vitro generation of gene signatures and elucidation of molecular mechanisms after exposure to toxic compounds. An important issue to be investigated is the in vivo relevance of in vitro-obtained data. In the current study, we compare the gene expression profiles generated in vitro, after exposing conventional and epigenetically stabilized primary rat hepatocytes to well-known genotoxic hepatocarcinogens (aflatoxin B1, 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone and 2-nitrofluorene) with those derived in vivo after oral exposure of rats to these compounds. Similar statistical tools were applied on both sets of data. The major molecular pathways affected in the in vivo setting were DNA damage, detoxification and cell survival response, as previously described. In the conventional hepatocyte cultures, two of the three genotoxicants showed quite similar responses as in vivo with respect to these pathways. The third compound (2-nitrofluorene) revealed in vitro response which was not observed in vivo. In the epigenetically stabilized hepatocytes, in contrast to what was expected, the responses were less relevant for the in vivo situation. This study highlights the importance of in vitro/in vivo comparison of data that are generated using in vitro models and shows that conventional primary rat hepatocyte cultures represent an appropriate in vitro model to retrieve mechanistic information on the exposure to genotoxicants.

Background: It has been suggested that fetal carcinogenic exposure might lead to predisposition to develop cancer during childhood or in later life possibly through modulation of the fetal transcriptome. Because gender effects in the incidence of childhood cancers have been described, we hypothesized differences at the transcriptomic level in cord blood between male and female newborns as a consequence of fetal carcinogenic exposure. The objective was to investigate whether transcriptomic responses to dietary genotoxic and nongenotoxic carcinogens show gender-specific mechanisms-of-action relevant for chemical carcinogenesis., Methods: Global gene expression was applied in umbilical cord blood samples, the CALUX-assay was used for measuring dioxin(-like), androgen(-like), and estrogen(-like) internal exposure, and acrylamide-hemoglobin adduct levels were determined by mass spectrometry adduct-FIRE-procedure(TM). To link gene expression to an established phenotypic biomarker of cancer risk, micronuclei frequencies were investigated., Results: While exposure levels did not differ between sexes at birth, important gender-specific differences were observed in gene expressions associated with these exposures linked with cell cycle, the immune system and more general cellular processes such as posttranslation. Moreover, oppositely correlating leukemia/lymphoma genes between male and female newborns were identified in relation to the different biomarkers of exposure that might be relevant to male-specific predisposition to develop these cancers in childhood. CONCLUSIONS/IMPACT: This study reveals different transcriptomic responses to environmental carcinogens between the sexes. In particular, male-specific TNF-alpha-NF-kB signaling upon dioxin exposure and activation of the Wnt-pathway in boys upon acrylamide exposure might represent possible mechanistic explanations for gender specificity in the incidence of childhood leukemia., (2012 AACR)

At present, substantial efforts are focused on the development of in vitro assays coupled with "omics" technologies for the identification of carcinogenic substances as an alternative to the classical 2-year rodent carcinogenicity bioassay. A prerequisite for the eventual regulatory acceptance of such assays, however, is the in vivo relevance of the observed in vitro findings. In the current study, hepatocarcinogen-induced gene expression profiles generated after the exposure of conventional cultures of primary rat hepatocytes to three non-genotoxic carcinogens (methapyrilene hydrochloride, piperonyl butoxide, and Wy-14643), three genotoxic carcinogens (aflatoxin B1, 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone, and 2-nitrofluorene), and two non-carcinogens (nifedipine and clonidine) are compared with previously obtained in vivo data after oral administration for up to 14 days of the same hepatocarcinogens to rats. In addition to the comparison of deregulated genes and functions per compound between in vivo and in vitro models, the major discriminating cellular pathways found in vivo in livers of exposed rats were examined for deregulation in vitro. Further, in vivo-derived gene signatures for the identification of genotoxic versus non-genotoxic carcinogens are used to classify in vitro-tested hepatocarcinogens and non-carcinogens. In the primary hepatocyte cultures, two out of the three tested genotoxic carcinogens mimicked the in vivo-relevant DNA damage response and were correctly assessed. Exposure to the non-genotoxic hepatocarcinogens, however, triggered a relatively weak response in the in vitro system, with no clear similarities to in vivo. This study contributes to the further optimization of toxicogenomics predictive tools when applied in in vitro settings.

Toxicological studies assessing the safety of compounds for humans frequently use in vitro systems to characterize toxic responses in combination with transcriptomic analyses. Thus far, changes have mostly been investigated at the mRNA level. Recently, microRNAs have attracted attention because they are powerful negative regulators of mRNA levels and, thus, may be responsible for the modulation of important mRNA networks implicated in toxicity. This study aimed to identify possible microRNA-mRNA networks as novel interactions on the gene expression level after a genotoxic insult. We used benzo[a]pyrene (BaP), a polycyclic aromatic hydrocarbon, as a model genotoxic/carcinogenic compound. We analyzed time-dependent effects on mRNA and microRNA profiles in HepG2 cells, a widely used human liver cell line that expresses active p53 and is competent for the biotransformation of BaP. Changes in microRNA expression in response to BaP, in combination with multiple alterations of mRNA levels, were observed. Many of these altered mRNAs are targets of altered microRNAs. Using pathway analysis, we evaluated the relevance of such microRNA deregulations to genotoxicity. This revealed eight microRNAs that appear to participate in specific BaP-responsive pathways relevant to genotoxicity, such as apoptotic signaling, cell cycle arrest, DNA damage response, and DNA damage repair. Our results particularly highlight the potential of microRNA-29b, microRNA-26a-1*, and microRNA-122* as novel players in the BaP response. Therefore, this study demonstrates the added value of an integrated microRNA-mRNA approach for identifying molecular mechanisms induced by BaP in an in vitro human model.

The primary cilium is an immotile sensory and signaling organelle found on the majority of mammalian cell types. Of the multitude of roles that the primary cilium performs, perhaps some of the most important include maintenance of differentiation, quiescence, and cellular polarity. Given that the progression of cancer requires disruption of all of these processes, we have investigated the effects of several carcinogens on the primary cilium of the RPTEC/TERT1 human proximal tubular epithelial cell line. Using both scanning electron microscopy and immunofluorescent labeling of the ciliary markers acetylated tubulin and Arl13b, we confirmed that RPTEC/TERT1 cells express primary cilium upon reaching confluence. Treatment with the carcinogens ochratoxin A (OTA) and potassium bromate (KBrO(3)) caused a significant reduction in the number of ciliated cells, while exposure to nifedipine, a noncarcinogenic renal toxin, had no effect on primary cilium expression. Flow cytometric analysis of the effects of all three compounds on the cell cycle revealed that only KBrO(3) resulted in an increase in the proportion of cells entering the cell cycle. Microarray analysis revealed dysregulation of multiple pathways affecting ciliogenesis and ciliary maintenance following OTA and KBrO(3) exposure, which were unaffected by nifedipine exposure. The primary cilium represents a unique physical checkpoint with relevance to carcinogenesis. We have shown that the renal carcinogens OTA and KBrO(3) cause significant deciliation in a model of the proximal tubule. With KBrO(3), this was followed by reentry into the cell cycle; however, deciliation was not found to be associated with reentry into the cell cycle following OTA exposure. Transcriptomic analysis identified dysregulation of Wnt signaling and ciliary trafficking in response to OTA and KBrO(3) exposure.

The safety assessment for pharmaceuticals includes in vivo repeated dose toxicity tests in laboratory animals. These in vivo studies often generate false negative results and unexpected toxicity. The appearance of this unexpected toxicity is one of the major reasons for the drawback of a drug from the market. The liver is often a target organ in toxicology since it is responsible for the metabolism and elimination of chemical compounds. Therefore, there is need for new screening methods which classify hepatotoxic compounds earlier in development. This will lead to safer drugs and a more efficient drug discovery process. Furthermore, these new screening methods are preferably in vitro test systems, aiming at reducing the use of laboratory animals. In this review the possibilities of proteomics and its promising results for improving current predictive and mechanistic toxicological studies are described. Biomarkers or protein panels for hepatotoxic mechanisms, which reflect the in vivo situation, need to be identified to allow a better toxicity screening. Therefore, in vivo studies and in vitro cell models are discussed and evaluated with regard to the protein expression of their metabolic enzymes, their similarities with liver, their use for analyzing toxicological mechanisms and hepatotoxicity screening. Studies in which proteomics are combined with other omics-technologies are also presented. The results from these integrated data analyses can be used for the development of improved panels of biomarkers for toxicity screening., (Copyright © 2012 Elsevier Ltd. All rights reserved.)

Acetaminophen is the primary cause of acute liver toxicity in Europe/USA, which led the FDA to reconsider recommendations concerning safe acetaminophen dosage/use. Unfortunately, the current tests for liver toxicity are no ideal predictive markers for liver injury, i.e. they only measure acetaminophen exposure after profound liver toxicity has already occurred. Furthermore, these tests do not provide mechanistic information. Here, 'omics techniques (global analysis of metabolomic/gene-expression responses) may provide additional insight. To better understand acetaminophen-induced responses at low doses, we evaluated the effects of (sub-)therapeutic acetaminophen doses on metabolite formation and global gene-expression changes (including, for the first time, full-genome human miRNA expression changes) in blood/urine samples from healthy human volunteers. Many known and several new acetaminophen-metabolites were detected, in particular in relation to hepatotoxicity-linked, oxidative metabolism of acetaminophen. Transcriptomic changes indicated immune-modulating effects (2g dose) and oxidative stress responses (4g dose). For the first time, effects of acetaminophen on full-genome human miRNA expression have been considered and confirmed the findings on mRNA level. 'Omics techniques outperformed clinical chemistry tests and revealed novel response pathways to acetaminophen in humans. Although no definitive conclusion about potential immunotoxic effects of acetaminophen can be drawn from this study, there are clear indications that the immune system is triggered even after intake of low doses of acetaminophen. Also, oxidative stress-related gene responses, similar to those seen after high dose acetaminophen exposure, suggest the occurrence of possible pre-toxic effects of therapeutic acetaminophen doses. Possibly, these effects are related to dose-dependent increases in levels of hepatotoxicity-related metabolites., (Copyright © 2012 Elsevier Inc. All rights reserved.)

Hepatocyte-like cells derived from the differentiation of human embryonic stem cells (hES-Hep) have potential to provide a human relevant in vitro test system in which to evaluate the carcinogenic hazard of chemicals. In this study, we have investigated this potential using a panel of 15 chemicals classified as noncarcinogens, genotoxic carcinogens, and nongenotoxic carcinogens and measured whole-genome transcriptome responses with gene expression microarrays. We applied an ANOVA model that identified 592 genes highly discriminative for the panel of chemicals. Supervised classification with these genes achieved a cross-validation accuracy of > 95%. Moreover, the expression of the response genes in hES-Hep was strongly correlated with that in human primary hepatocytes cultured in vitro. In order to infer mechanistic information on the consequences of chemical exposure in hES-Hep, we developed a computational method that measures the responses of biochemical pathways to the panel of treatments and showed that these responses were discriminative for the three toxicity classes and linked to carcinogenesis through p53, mitogen-activated protein kinases, and apoptosis pathway modules. It could further be shown that the discrimination of toxicity classes was improved when analyzing the microarray data at the pathway level. In summary, our results demonstrate, for the first time, the potential of human embryonic stem cell--derived hepatic cells as an in vitro model for hazard assessment of chemical carcinogenesis, although it should be noted that more compounds are needed to test the robustness of the assay.

Background: We hypothesized that in Flanders (Belgium), the prevalence of at-risk genotypes for genotoxic effects decreases with age due to morbidity and mortality resulting from chronic diseases. Rather than polymorphisms in single genes, the interaction of multiple genetic polymorphisms in low penetrance genes involved in genotoxic effects might be of relevance., Methods: Genotyping was performed on 399 randomly selected adults (aged 50-65) and on 442 randomly selected adolescents. Based on their involvement in processes relevant to genotoxicity, 28 low penetrance polymorphisms affecting the phenotype in 19 genes were selected (xenobiotic metabolism, oxidative stress defense and DNA repair, respectively 13, 6 and 9 polymorphisms). Polymorphisms which, based on available literature, could not clearly be categorized a priori as leading to an 'increased risk' or a 'protective effect' were excluded., Results: The mean number of risk alleles for all investigated polymorphisms was found to be lower in the 'elderly' (17.0 ± 2.9) than the 'adolescent' (17.6 ± 3.1) subpopulation (P = 0.002). These results were not affected by gender nor smoking. The prevalence of a high (> 17 = median) number of risk alleles was less frequent in the 'elderly' (40.6%) than the 'adolescent' (51.4%) subpopulation (P = 0.002). In particular for phase II enzymes, the mean number of risk alleles was lower in the 'elderly' (4.3 ± 1.6 ) than the 'adolescent' age group (4.8 ± 1.9) P < 0.001 and the prevalence of a high (> 4 = median) number of risk alleles was less frequent in the 'elderly' (41.3%) than the adolescent subpopulation (56.3%, P < 0.001). The prevalence of a high (> 8 = median) number of risk alleles for DNA repair enzyme-coding genes was lower in the 'elderly' (37,3%) than the 'adolescent' subpopulation (45.6%, P = 0.017)., Conclusions: These observations are consistent with the hypothesis that, in Flanders, the prevalence of at-risk alleles in genes involved in genotoxic effects decreases with age, suggesting that persons carrying a higher number of at risk alleles (especially in phase II xenobiotic-metabolizing or DNA repair genes) are at a higher risk of morbidity and mortality from chronic diseases. Our findings also suggest that, regarding risk of disease associated with low penetrance polymorphisms, multiple polymorphisms should be taken into account, rather than single ones.

The current method to predict carcinogenicity of chemicals or drugs is the chronic 2-year rodent bioassay, which has disadvantages in duration, animal use, and specificity. An attractive alternative is the DNA repair-deficient Xpa(-/-)p53(+/-) mouse model that is sensitive to both genotoxic and nongenotoxic carcinogens. A next step in alternative carcinogenicity testing is the development of reliable in vitro systems. We investigated the use of primary hepatocytes, isolated from wild-type (WT) and Xpa(-/-)p53(+/-) mice, in combination with transcriptome analyses for their usefulness to predict carcinogenic features of compounds. As a proof of principle, we studied the response of hepatocytes to the genotoxic carcinogen benzo[a]pyrene (B[a]P). Upon treatment, both WT and Xpa(-/-)p53(+/-) hepatocytes appeared to be metabolically active. However, Xpa(-/-)p53(+/-) hepatocytes were more sensitive than WT hepatocytes to B[a]P treatment in terms of cell survival. In B[a]P-treated WT hepatocytes, DNA repair and cell cycle control genes were transcriptionally activated. Xpa(-/-)p53(+/-) hepatocytes were more responsive to B[a]P exposure, resulting in the downregulation of cancer-related pathways. Deregulation of mitogen-activated protein kinase signaling seems to play an essential role in this and might be the underlying reason for the increased susceptibility of Xpa(-/-)p53(+/-) mice toward carcinogens. Our conclusion is that primary hepatocytes combined with transcriptomics are promising to identify the carcinogenic features of chemicals. Furthermore, these cells seem suitable to gain further insight into the molecular mechanisms of the increased sensitivity of Xpa(-/-)p53(+/-) mice toward both genotoxic and nongenotoxic carcinogens.

In the present study, the effect of Trichostatin A (TSA), a histone deacetylase inhibitor, was investigated on the microRNA (miR, miRNA) expression profile in cultured primary rat hepatocytes by means of microarray analysis. Simultaneously, albumin secretory capacity and morphological features of the hepatocytes were evaluated throughout the culture time. In total, 25 out of 348 miRNAs were found to be differentially expressed between freshly isolated hepatocytes and 7-day cultured cells. Nineteen of these miRNAs were connected with 'general metabolism'. miR-21 and miR-126 were shown to be the most up and down regulated miRs upon cultivation and could be linked to the proliferative response triggered in the hepatocytes upon their isolation from the liver. miR-379 and miR-143, on the other hand, were found to be the most up and down regulated miRs upon TSA treatment. Together with the higher expression of miR-122 observed in TSA-treated versus non-treated cultures, we hypothesize that the changes observed for miR-122, miR-143 and miR-379 could be related to the inhibitory effects of TSA on hepatocellular proliferation., (Copyright © 2011 Elsevier Ltd. All rights reserved.)

Diet plays a decisive role in promoting or preventing colon cancer. However, the specific effects of some nutrients remain unclear. The capacity of fruit and vegetables to prevent cancer has been associated with their fiber and antioxidant composition. We investigated whether consumption of a lyophilized red grape pomace containing proanthocyanidin-rich dietary fiber (grape antioxidant dietary fiber, GADF) by female C57BL/6J mice would affect the serum metabolic profile or colon mucosa gene expression using NMR techniques and DNA microarray, respectively. The mice were randomly assigned to 2 groups that for 2 wk consumed a standard rodent diet and were gavaged with 100 mg/kg body weight GADF suspended in water or an equivalent volume of plain tap water (10 mL/kg body weight). The amount of fiber supplemented was calculated to equal the current recommended daily levels of fiber consumption for humans. The inclusion of dietary GADF induced alterations in the expression of tumor suppressor genes and proto-oncogenes as well as the modulation of genes from pathways, including lipid biosynthesis, energy metabolism, cell cycle, and apoptosis. Overexpression of enzymes pertaining to the xenobiotic detoxifying system and endogenous antioxidant cell defenses was also observed. In summary, the genetic and metabolic profiles induced by GADF were consistent with the preventive effects of fiber and polyphenols. On the basis of these observations, we propose that GADF may contribute to reducing the risk of colon cancer.

Background: The integration of different 'omics' technologies has already been shown in several in vivo studies to offer a complementary insight into cellular responses to toxic challenges. Being interested in developing in vitro cellular models as alternative to animal-based toxicity assays, we hypothesize that combining transcriptomics and metabonomics data improves the understanding of molecular mechanisms underlying the effects caused by a toxic compound also in vitro in human cells. To test this hypothesis, and with the focus on non-genotoxic carcinogenesis as an endpoint of toxicity, in the present study, the human hepatocarcinoma cell line HepG2 was exposed to the well-known environmental carcinogen 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD)., Results: Transcriptomics as well as metabonomics analyses demonstrated changes in TCDD-exposed HepG2 in common metabolic processes, e.g. amino acid metabolism, of which some of the changes only being confirmed if both 'omics' were integrated. In particular, this integrated analysis identified unique pathway maps involved in receptor-mediated mechanisms, such as the G-protein coupled receptor protein (GPCR) signaling pathway maps, in which the significantly up-regulated gene son of sevenless 1 (SOS1) seems to play an important role. SOS1 is an activator of several members of the RAS superfamily, a group of small GTPases known for their role in carcinogenesis., Conclusions: The results presented here were not only comparable with other in vitro studies but also with in vivo studies. Moreover, new insights on the molecular responses caused by TCDD exposure were gained by the cross-omics analysis.

The birth cohort BraMat (n = 205; a sub-cohort of the Norwegian Mother and Child Cohort Study (MoBa) conducted by the Norwegian Institute of Public Health) was established to study whether prenatal exposure to toxicants from the maternal diet affects immunological health outcomes in children. We here report on the environmental pollutants polychlorinated biphenyls (PCBs) and dioxins, as well as acrylamide generated in food during heat treatment. The frequency of common infections, eczema or itchiness, and periods of more than 10 days of dry cough, chest tightness or wheeze (called wheeze) in the children during the first year of life was assessed by questionnaire data (n = 195). Prenatal dietary exposure to the toxicants was estimated using a validated food frequency questionnaire from MoBa. Prenatal exposure to PCBs and dioxins was found to be associated with increased risk of wheeze and exanthema subitum, and also with increased frequency of upper respiratory tract infections. We found no associations between prenatal exposure to acrylamide and the health outcomes investigated. Our results suggest that prenatal dietary exposure to dioxins and PCBs may increase the risk of wheeze and infectious diseases during the first year of life., (Copyright © 2011 Elsevier Ltd. All rights reserved.)

Background: In vitro cell systems together with omics methods represent promising alternatives to conventional animal models for toxicity testing. Transcriptomic and proteomic approaches have been widely applied in vitro but relatively few studies have used metabolomics. Therefore, the goal of the present study was to develop an untargeted methodology for performing reproducible metabolomics on in vitro systems. The human liver cell line HepG2, and the well-known hepatotoxic and non-genotoxic carcinogen 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), were used as the in vitro model system and model toxicant, respectively., Results: The study focused on the analysis of intracellular metabolites using NMR, LC-MS and GC-MS, with emphasis on the reproducibility and repeatability of the data. State of the art pre-processing and alignment tools and multivariate statistics were used to detect significantly altered levels of metabolites after exposing HepG2 cells to TCDD. Several metabolites identified using databases, literature and LC-nanomate-Orbitrap analysis were affected by the treatment. The observed changes in metabolite levels are discussed in relation to the reported effects of TCDD., Conclusions: Untargeted profiling of the polar and apolar metabolites of in vitro cultured HepG2 cells is a valid approach to studying the effects of TCDD on the cell metabolome. The approach described in this research demonstrates that highly reproducible experiments and correct normalization of the datasets are essential for obtaining reliable results. The effects of TCDD on HepG2 cells reported herein are in agreement with previous studies and serve to validate the procedures used in the present work.

The 7th amendment to the EU Cosmetics Directive prohibits to put animal-tested cosmetics on the market in Europe after 2013. In that context, the European Commission invited stakeholder bodies (industry, non-governmental organisations, EU Member States, and the Commission's Scientific Committee on Consumer Safety) to identify scientific experts in five toxicological areas, i.e. toxicokinetics, repeated dose toxicity, carcinogenicity, skin sensitisation, and reproductive toxicity for which the Directive foresees that the 2013 deadline could be further extended in case alternative and validated methods would not be available in time. The selected experts were asked to analyse the status and prospects of alternative methods and to provide a scientifically sound estimate of the time necessary to achieve full replacement of animal testing. In summary, the experts confirmed that it will take at least another 7-9 years for the replacement of the current in vivo animal tests used for the safety assessment of cosmetic ingredients for skin sensitisation. However, the experts were also of the opinion that alternative methods may be able to give hazard information, i.e. to differentiate between sensitisers and non-sensitisers, ahead of 2017. This would, however, not provide the complete picture of what is a safe exposure because the relative potency of a sensitiser would not be known. For toxicokinetics, the timeframe was 5-7 years to develop the models still lacking to predict lung absorption and renal/biliary excretion, and even longer to integrate the methods to fully replace the animal toxicokinetic models. For the systemic toxicological endpoints of repeated dose toxicity, carcinogenicity and reproductive toxicity, the time horizon for full replacement could not be estimated.

The murine embryonic stem cell test (EST) is designed to evaluate developmental toxicity based on compound-induced inhibition of embryonic stem cell (ESC) differentiation into cardiomyocytes. The addition of transcriptomic evaluation within the EST may result in enhanced predictability and improved characterization of the applicability domain, therefore improving usage of the EST for regulatory testing strategies. Transcriptomic analyses assessing factors critical for risk assessment (i.e. dose) are needed to determine the value of transcriptomic evaluation in the EST. Here, using the developmentally toxic compound, flusilazole, we investigated the effect of compound concentration on gene expression regulation and toxicity prediction in ESC differentiation cultures. Cultures were exposed for 24 h to multiple concentrations of flusilazole (0.54-54 μM) and RNA was isolated. In addition, we sampled control cultures 0, 24, and 48 h to evaluate the transcriptomic status of the cultures across differentiation. Transcriptomic profiling identified a higher sensitivity of development-related processes as compared to cell division-related processes in flusilazole-exposed differentiation cultures. Furthermore, the sterol synthesis-related mode of action of flusilazole toxicity was detected. Principal component analysis using gene sets related to normal ESC differentiation was used to describe the dynamics of ESC differentiation, defined as the 'differentiation track'. The concentration-dependent effects on development were reflected in the significance of deviation of flusilazole-exposed cultures from this transcriptomic-based differentiation track. Thus, the detection of developmental toxicity in EST using transcriptomics was shown to be compound concentration-dependent. This study provides further insight into the possible application of transcriptomics in the EST as an improved alternative model system for developmental toxicity testing., (Copyright © 2010 Elsevier Inc. All rights reserved.)

Unexpected hepatotoxicity is one of the major reasons of drugs failing in clinical trials. This emphasizes the need for new screening methods that address toxicological hazards early in the drug discovery process. Here, proteomics techniques were used to gain further insight into the mechanistic processes of the hepatotoxic compounds. Drug-induced hepatotoxicity is mainly divided in hepatic steatosis, cholestasis, or necrosis. For each class, a compound was selected, respectively amiodarone, cyclosporin A, and acetaminophen. The changes in protein expressions in HepG2, after exposure to these test compounds, were studied using quantitative two-dimensional differential gel electrophoresis. Identification of differentially expressed proteins was performed by Maldi-TOF/TOF MS and liquid chromatography-tandem mass spectrometry. In this study, 254 differentially expressed protein spots were detected in a two-dimensional proteome map from which 86 were identified, showing that the proteome of HepG2 cells is responsive to hepatotoxic compounds. cyclosporin A treatment was responsible for most differentially expressed proteins and could be discriminated in the hierarchical clustering analysis. The identified differential proteins show that cyclosporin A may induce endoplasmic reticulum (ER) stress and disturbs the ER-Golgi transport, with an altered vesicle-mediated transport and protein secretion as result. Moreover, the differential protein pattern seen after cyclosporin A treatment can be related to cholestatic mechanisms. Therefore, our findings indicate that the HepG2 in vitro cell system has distinctive characteristics enabling the assessment of cholestatic properties of novel compounds at an early stage of drug discovery.

Hepatic systems toxicology is the integrative analysis of toxicogenomic technologies, e.g., transcriptomics, proteomics, and metabolomics, in combination with traditional toxicology measures to improve the understanding of mechanisms of hepatotoxic action. Hepatic toxicology studies that have employed toxicogenomic technologies to date have already provided a proof of principle for the value of hepatic systems toxicology in hazard identification. In the present review, acetaminophen is used as a model compound to discuss the application of toxicogenomics in hepatic systems toxicology for its potential role in the risk assessment process, to progress from hazard identification towards hazard characterization. The toxicogenomics-based parallelogram is used to identify current achievements and limitations of acetaminophen toxicogenomic in vivo and in vitro studies for in vitro-to-in vivo and interspecies comparisons, with the ultimate aim to extrapolate animal studies to humans in vivo. This article provides a model for comparison of more species and more in vitro models enhancing the robustness of common toxicogenomic responses and their relevance to human risk assessment. To progress to quantitative dose-response analysis needed for hazard characterization, in hepatic systems toxicology studies, generation of toxicogenomic data of multiple doses/concentrations and time points is required. Newly developed bioinformatics tools for quantitative analysis of toxicogenomic data can aid in the elucidation of dose-responsive effects. The challenge herein is to assess which toxicogenomic responses are relevant for induction of the apical effect and whether perturbations are sufficient for the induction of downstream events, eventually causing toxicity., (© 2010 Elsevier Inc. All rights reserved.)

The murine embryonic stem cell test (EST) is an alternative testing method designed to assess potential developmental toxicity of compounds. The implementation of transcriptomics in the EST has been shown to reduce the culture duration and improve endpoint evaluation and is expected to result in an enhanced predictability and definition of the applicability domain. We evaluated the identification of developmental toxicity in the EST using two gene sets ("Van&#95;Dartel&#95;heartdiff&#95;24h" and "EST biomarker genes") defined in our earlier studies. Nonexposed embryonic stem cells (ESC) differentiation cultures were sampled 0, 24, and 48 h after initiation of differentiation. Additionally, cultures exposed to 12 diverse well-characterized positive and negative developmental toxicants were isolated 24 h after the onset of exposure. Inhibition of ESC differentiation was evaluated in parallel by morphological scoring on culture day 10. Transcriptomics analysis was conducted using the Affymetrix Gene Chips platform. We applied principal component analysis on the basis of the two predefined gene sets to define the "differentiation track" that represents ESC differentiation. The significance of derivations in the gene expression-based differentiation track because of compound exposures were evaluated to determine developmental toxicity of tested compounds. We successfully predicted developmental toxicity using transcriptomics for 83% (10/12) and 67% (8/12) of the compounds, respectively, using the two predefined gene sets ("Van&#95;Dartel&#95;heartdiff&#95;24h" and "EST biomarker genes"). Our study suggests that the application of transcriptomics may improve the applicability of the EST for the prediction of the developmental toxicity of chemicals.

Investigating the immunotoxic effects of exposure to chemicals usually comprises evaluation of weight and histopathology of lymphoid tissues, various lymphocyte parameters in the circulation, and immune function. Immunotoxicity assessment is time consuming in humans or requires a high number of animals, making it expensive. Furthermore, reducing the use of animals in research is an important ethical and political issue. Immunotoxicogenomics represents a novel approach to investigate immunotoxicity able of overcoming these limitations. The current research, embedded in the European Union project NewGeneris, aimed to retrieve gene expression profiles that are indicative of exposure to immunotoxicants. To this end, whole-genome gene expression was investigated in human peripheral blood mononuclear cells in response to in vitro exposure to a range of immunotoxic chemicals (4-hydroxy-2-nonenal, aflatoxin B1, benzo[a]pyrene, deoxynivalenol, ethanol, malondialdehyde, polychlorinated biphenyl 153, and 2,3,7,8-tetrachlorodibenzo-p-dioxin) and nonimmunotoxic chemicals (acrylamide, dimethylnitrosamine, 2-amino-3-methyl-3H-imidazo[4,5-F]quinoline, and 2-amino-1-methyl-6-phenylimidazo[4,5-b]pyridine). Using Agilent oligonucleotide microarrays, whole-genome gene expression profiles were generated, which were analyzed using Genedata's Expressionist software. Using Recursive Feature Elimination and Support Vector Machine, a set of 48 genes was identified that distinguishes the immunotoxic from the nonimmunotoxic compounds. Analysis for enrichment of biological processes showed the gene set to be highly biologically and immunologically relevant. We conclude that we have identified a promising transcriptomic profile indicative of immunotoxic exposure.

Chemical carcinogens may cause a multitude of effects inside cells, thereby affecting transcript levels of genes by direct activation of transcription factors (TF) or indirectly through the formation of DNA damage. As the temporal profiles of these responses may be profoundly different, examining time-dependent changes may provide new insights in TF networks related to cellular responses to chemical carcinogens. Therefore, we investigated in human hepatoma cells gene expression changes caused by benzo[a]pyrene at 12 time points after exposure, in relation to DNA adduct and cell cycle. Temporal profiles for functional gene sets demonstrate both early and late effects in up- and downregulation of relevant gene sets involved in cell cycle, apoptosis, DNA repair, and metabolism of amino acids and lipids. Many significant transcription regulation networks appeared to be around TF that are proto-oncogenes or tumor suppressor genes. The time series analysis tool Short Time-series Expression Miner (STEM) was used to identify time-dependent correlation of pathways, gene sets, TF networks, and biological parameters. Most correlations are with DNA adduct levels, which is an early response, and less with the later responses on G1 and S phase cells. The majority of the modulated genes in the Reactome pathways can be regulated by several of these TF, e.g., 73% by nuclear factor-kappa B and 34-42% by c-MYC, SRF, AP1, and E2F1. All these TF can also regulate one or more of the others. Our data indicate that a complex network of a few TF is responsible for the majority of the transcriptional changes induced by BaP. This network hardly changes over time, despite that the transcriptional profiles clearly alter, suggesting that also other regulatory mechanisms are involved.

In recent decades, our knowledge of the genetics and functional genomics of drug-metabolizing enzymes has increased and a wealth of data on drug-related 'omics' has become available. Despite the availability of large amounts of biological information on xenobiotic biotransformation, the number of available biotransformation pathway maps that can easily be used for visualization of multiple omics data is limited. Here, we created integrated biotransformation pathway maps suitable for multiple omics analysis using PathVisio. The ease of visualizing data on these maps was demonstrated by using published microarray data from human hepatocyte-like cell models, exemplifying - where a sufficient capacity for metabolizing chemicals is a prerequisite for a suited model - how the biotransformation pathway maps can be used for model selection., (Copyright © 2010 Elsevier Ltd. All rights reserved.)

The embryonic stem cell test (EST) has been designed to predict developmental toxicity based upon compound-induced inhibition of embryonic stem cell (ESC) differentiation. The end point scoring, the test duration, and the definition of the predictivity and the applicability domain require improvements to facilitate implementation of the EST into regulatory testing strategies. The use of transcriptomics to study compound-induced differentiation modulation may improve the EST in each of these aspects. ESC differentiation was induced, and cell samples were collected after 0, 24, and 48 h of differentiation. Additionally, samples were collected that were 24 h exposed to one of five developmentally toxic compounds or a nondevelopmentally toxic compound. All samples were hybridized to Affymetrix GeneChips, and analyses revealed that 26 genes were significantly regulated both during ESC differentiation and by exposure to each of the developmentally toxic compounds tested. Using principal component analysis, we defined a "differentiation track" on the basis of this gene list, which represents ESC differentiation. We showed that significant deviation from the differentiation track was in line with the developmental toxic properties of the compounds. The significance of deviation was analyzed using the leave-one-out cross-validation, which showed a favorable prediction of toxicity in the system. Our findings show that gene expression signatures can be used to identify developmental toxicant-induced differentiation modulation. In addition, studying compound-induced effects at an early stage of differentiation combined with transcriptomics leads to increased objectivity in determining differentiation inhibition and to a reduction of the test duration. Furthermore, this approach may improve the predictivity and applicability domain of the EST.

Direct comparison of the hepatoma cell lines HepG2 and HepaRG has previously been performed by only evaluating a limited set of genes or proteins. In this study, we examined the whole-genome gene expression of both cell lines before and after exposure to the genotoxic (GTX) carcinogens aflatoxin B1 and benzo[a]pyrene and the nongenotoxic (NGTX) carcinogens cyclosporin A, 17beta-estradiol, and 2,3,7,8-tetrachlorodibenzo-para-dioxin for 12 and 48 h. Before exposure, this analysis revealed an extensive network of genes and pathways, which were regulated differentially for each cell line. The comparison of the basal gene expression between HepG2, HepaRG, primary human hepatocytes (PHH), and liver clearly showed that HepaRG resembles PHH and liver the most. After exposure to the GTX and NGTX carcinogens, for both cell lines, common pathways were found that are important in carcinogenesis, for example, cell cycle regulation and apoptosis. However, also clear differences between exposed HepG2 and HepaRG were observed, and these are related to common metabolic processes, immune response, and transcription processes. Furthermore, HepG2 performs better in discriminating between GTX and NGTX carcinogens. In conclusion, these results have shown that HepaRG is a more suited in vitro liver model for biological interpretations of the effects of exposure to chemicals, whereas HepG2 is a more promising in vitro liver model for classification studies using the toxicogenomics approach. Although, it should be noted that only five carcinogens were used in this study.

Reactive oxygen species-induced oxidative stress in the colon is involved in inflammatory bowel diseases and suggested to be associated with colorectal cancer risk. However, our insight in molecular responses to different oxygen radicals is still fragmentary. Therefore, we studied global gene expression by an extensive time series (0.08, 0.25, 0.5, 1, 2, 4, 8, 16, or 24 h) analyses in human colon cancer (caco-2) cells after exposure to H(2)O(2) or the superoxide anion donor menadione. Differences in gene expression were investigated by hybridization on two-color microarrays against nonexposed time-matched control cells. Next to gene expression, correlations with related phenotypic markers (8-oxodG levels and cell cycle arrest) were investigated. Gene expression analysis resulted in 1404 differentially expressed genes upon H(2)O(2) challenge and 979 genes after menadione treatment. Further analysis of gene expression data revealed how these oxidant responses can be discriminated. Time-dependent coregulated genes immediately showed a pulse-like response to H(2)O(2), while the menadione-induced expression is not restored over 24 h. Pathway analyses demonstrated that H(2)O(2) immediately influences pathways involved in the immune function, while menadione constantly regulated cell cycle-related pathways Altogether, this study offers a novel and detailed insight in the similarities and differences of the time-dependent oxidative stress responses induced by the oxidants H(2)O(2) and menadione and show that these can be discriminated regarding their modulation of particular colon carcinogenesis-related mechanisms.

Assessing the potential carcinogenicity of chemicals for humans represents an ongoing challenge. Chronic rodent bioassays predict human cancer risk at only limited reliability and are simultaneously expensive and long lasting. In order to seek for alternatives, the ability of a transcriptomics-based primary mouse hepatocyte model to classify carcinogens by their modes of action was evaluated. As it is obvious that exposure will induce a cascade of gene expression modifications, in particular, the influence of exposure time in vitro on discriminating genotoxic (GTX) carcinogens from nongenotoxic (NGTX) carcinogens class discrimination was investigated. Primary mouse hepatocytes from male C57Bl6 mice were treated for 12, 24, 36, and 48 h with two GTX and two NGTX carcinogens. For validation, two additional GTX compounds were studied at 24 and 48 h. Immunostaining of gammaH2AX foci was applied in order to phenotypically verify DNA damage. It confirmed significant induction of DNA damage after treatment with GTX compounds but not with NGTX compounds. Whole-genome gene expression modifications were analyzed by means of Affymetrix microarrays. When using differentially expressed genes from data sets normalized by Robust Multi-array Average, the two classes and various compounds were better separated from each other by hierarchical clustering when increasing the treatment period. Discrimination of GTX and NGTX carcinogens by Prediction Analysis of Microarray improved with time and resulted in correct classification of the validation compounds. The present study shows that gene expression profiling in primary mouse hepatocytes is promising for discriminating GTX from NGTX compounds and that this discrimination improves with increasing treatment period.

A compound for which marked species differences have been reported in laboratory animals and humans is coumarin. In rats, metabolites of coumarin are highly toxic, whereas in humans, the compound is mainly metabolized to non-toxic metabolites. In the present study, a toxicogenomics-based parallelogram approach was used to compare effects of coumarin on gene expression in human hepatocytes relevant for the situation in vivo. To this purpose, gene expression profiling was performed on human hepatocytes treated with coumarin in a pharmacological relevant and proposed toxic concentration and results were compared to a previously performed coumarin in vivo and in vitro rat toxicogenomics study. No cytotoxicity was observed in human hepatocytes at both concentrations, whereas rats showed clear toxic effects in vitro as well as in vivo. In all three systems, coumarin affected genes involved in the blood coagulation pathway; this indicates relevant responses in cases of human exposure. However, no pathways and processes related to hepatotoxicity in rats were observed in human hepatocytes. Still, repression of energy-consuming biochemical pathways and impairment of mitochondrial function were observed in human hepatocytes treated with the highest concentration of coumarin, possibly indicating toxicity. In conclusion, although species differences in response to coumarin are evident in the present results, the toxicogenomics-based parallelogram approach enables clear discrimination between pharmacological responses at pharmacological doses and proposed toxic responses at high (toxic) doses relevant for humans in vivo.

Metabolites of the human carcinogen 4-aminobiphenyl (4-ABP) form hemoglobin (Hb) adducts, which represent a useful biomarker for exposure. However, not every individual responds to a similar degree to 4-ABP exposure, and variations in 4-ABP-Hb adduct formation might be explained by genetic polymorphisms in genes coding for enzymes involved in 4-ABP metabolism. 4-ABP-Hb adducts were measured in blood samples from 57 smoking and 10 non-smoking volunteers. An association was found between cigarette smoking and 4-ABP-Hb adduct levels in smokers (R(2) = 0.5, P < 0.001). Subsequently, subjects were genotyped for 12 polymorphisms in seven genes involved in biotransformation reactions. From this selection of polymorphisms, a significant impact was found for the CYP1B1 Leu(432)Val polymorphism (P = 0.021), which has been reported to lead to a decrease in enzyme activity. Indeed higher levels of 4-ABP-Hb adducts were observed in homo- and heterozygous carriers of the CYP1B1 (432)Leu as compared to the double CYP1B1 (432)Val genotype. A significant interaction between these CYP1B1 genotypes and the level of exposure was found (P = 0.003). Noteworthy, a saturation effect was observed for 4-ABP-Hb adduct formation at high smoking doses limited to carriers of the CYP1B1 (432)Leu allele. No effect of polymorphisms in other genes were found. This is the first study in humans suggesting a crucial role of the CYP1B1 enzyme in 4-ABP metabolism, indicating a protective effect of the CYP1B1 Leu(432)Val polymorphism against the formation of 4-ABP-Hb adduct levels, depending on the smoking dose.

Paternal exposure to genotoxic compounds is thought to contribute to diseases in their offspring. Therefore, it is of importance to develop biomarkers of male germ cell exposure to genotoxins. Unfortunately, the testis cannot be reached for routine biomonitoring, but mRNA-profiles in spermatozoa may reflect the processes that have occurred in the testis after exposures to genotoxins, since spermatozoa are largely transcriptionally inactive. Therefore, mRNA profiles from sperm in ejaculates of cigarette smokers (N=4) were compared with non-smokers (N= 4). Smoking behaviour was verified by assessing cotinine levels in seminal plasma. High expression of the germ cell specific gene protamine 2 (PRM2) was observed in spermatozoal mRNA isolates by Q-PCR, which was absent in reference mRNA isolates obtained from a pool of other organs. Gene-expression analysis was subsequently performed using microarray technology and a total of 781 genes were found to be differentially expressed in spermatozoa of smokers compared to non-smokers (fold change >40%; p < 0.05). To further limit the number of false positive results, genes were additionally selected on basis of the correlation between their expression levels with cotinine concentrations in seminal plasma (r > 0.80 as arbitrary cut-off value, p < 0.05), and a total of 200 transcripts remained, of which the germ cell specific transcription factor SALF was the highest up-regulated gene (5.4-fold) and the zinc finger encoding gene TRIM26 most down regulated (7.4-fold). Although no altered pathways could be identified for the differentially expressed genes, an enrichment was observed for NF-kappaB regulated genes (46% vs. 27%, p = 0.004) playing a central role in stress response. Accordingly, subsequent analysis of transcription factor networks suggests that apoptosis was inhibited in smokers. These data show the feasibility of using gene-expression profiles in mature sperm to elucidate gene-environment interactions in male testis.