Search

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.

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.

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.

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.

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.

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.

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_Dartel_heartdiff_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_Dartel_heartdiff_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.

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.

The frequent use of rodent hepatic in vitro systems in pharmacological and toxicological investigations challenges extrapolation of in vitro results to the situation in vivo and interspecies extrapolation from rodents to humans. The toxicogenomics approach may aid in evaluating relevance of these model systems for human risk assessment by direct comparison of toxicant-induced gene expression profiles and infers mechanisms between several systems. In the present study, acetaminophen (APAP) was used as a model compound to compare gene expression responses between rat and human using in vitro cellular models, hepatocytes, and between rat in vitro and in vivo. Comparison at the level of modulated biochemical pathways and biological processes rather than at that of individual genes appears preferable as it increases the overlap between various systems. Pathway analysis by T-profiler revealed similar biochemical pathways and biological processes repressed in rat and human hepatocytes in vitro, as well as in rat liver in vitro and in vivo. Repressed pathways comprised energy-consuming biochemical pathways, mitochondrial function, and oxidoreductase activity. The present study is the first that used a toxicogenomics-based parallelogram approach, extrapolating in vitro to in vivo and interspecies, to reveal relevant mechanisms indicative of APAP-induced liver toxicity in humans in vivo.

Although exposure to polycyclic aromatic hydrocarbons (PAHs) occurs mostly through mixtures, hazard and risk assessment are mostly based on the effects caused by individual compounds. The objective of the current study was to investigate whether interactions between PAHs occur, focusing on gene expression (as measured by cDNA microarrays) and DNA adduct formation. The effects of benzo[a]pyrene or dibenzo[a,h]anthracene (DB[a,h]A) alone and in binary mixtures with another PAH (DB[a,h]A, benzo[b]fluoranthene, fluoranthene or dibenzo[a,l]pyrene) were investigated using precision-cut rat liver slices. All compounds significantly modulated the expression of several genes, but overlap between genes affected by the mixture and by the individual compounds was relatively small. All mixtures showed an antagonistic response on total gene expression profiles. Moreover, at the level of individual genes, mostly antagonism was evident, with additivity and synergism observed for only a few genes. As far as DNA adduct formation is concerned, the binary mixtures generally caused antagonism. The effects in liver slices suggest a lower carcinogenic potency of PAH mixtures than estimated based on additivity of individual compounds.

Polycyclic aromatic hydrocarbons (PAHs) differ markedly in their carcinogenic potencies. Differences in transcriptomic responses upon PAH exposures might improve our current understanding of the differences in carcinogenicity, and therefore gene expression modulation by six PAHs in precision-cut rat liver slices was investigated. Gene expression modulation by benzo[a]pyrene (B[a]P), dibenzo[a,l]pyrene (DB[a,l]P), benzo[b]fluoranthene (B[b]F), fluoranthene (FA), dibenzo[a,h]anthracene (DB[a,h]A) and 1-methylphenanthrene (1-MPA) was assessed after 6- (B[a]P, DB[a,l]P) and 24-h (all compounds) exposure, using oligonucleotide arrays. DNA-adduct formation was determined using (32)P-post-labelling. The effects of PAHs on gene expression and on DNA-adduct formation were much more pronounced after 24-h exposure than after a 6-h exposure. Each compound induced gene expression changes dose-dependently and gene expression profiles were generally compound-specific. B[a]P, B[b]F and DB[a,h]A displayed comparable gene expression profiles, and so did DB[a,l]P, FA and 1-MPA. Only the carcinogenic PAHs (B[a]P, B[b]F, DB[a,l]P and DB[a,h]A) induced the oxidative stress pathway. DNA-adduct levels were: DB[a,l]P >> B[a]P > B[b]F > or = DB[a,h]A > FA > or = 1-MPA. The expression of only a few genes was found to correlate significantly with DNA-adduct formation, carcinogenic potency or Ah-receptor binding capacity (the last two taken from literature). These genes differed between the parameters. Our results indicate that PAHs generally induce a compound-specific response on gene expression and that discrimination of carcinogenic from non-carcinogenic compounds is partly feasible using this approach. Only at a specific pathway level, namely oxidative stress response, PAHs with high and low carcinogenic potency could be discriminated.

In current molecular epidemiology studies, a wide range of methods are used to monitor early biological effects after exposure to xenobiotic agents. Gene expression profiling is considered a promising tool that may provide more sensitive, mechanism-based biomarkers. As a first step toward obtaining information on the applicability of gene expression profiles as a biomarker for early biological effects of carcinogen exposure, we conducted in vitro studies on human peripheral blood mononuclear cells (PBMC). We used cigarette smoke condensate (CSC) and a selection of its genotoxic constituents as model agents, applying cDNA microarray technology to investigate modulated gene expression. In independent experiments using cells from several donors, quiescent PBMC were exposed for 18 h, followed by gene expression analyses on a microarray containing 600 toxicologically relevant genes. The search for candidate biomarker genes was binomial: first we looked for genes responding similarly to all agents; second, for agent-specific genes. Many genes were significantly deregulated by all compounds, but as the direction of deregulation frequently differed per agent, they are not useful as generic biomarkers. Cigarette smoke condensate modulated the expression of many more genes than any of its constituents, with the largest effect in SERPINB2. The affected genes are involved in immune or stress responses, but surprisingly no genes involved in DNA damage response were modulated, and only a few in DNA repair. In conclusion, several genes have been identified as potential biomarkers for population studies on early biological effects caused by cigarette smoke exposure, but no genes were identified that represent a generic biomarker.