Search

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.

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.

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.