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.

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)

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

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.

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.

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.

Primary human and rat hepatocyte cultures are well established in vitro systems used in toxicological studies. However, whereas transgenic mouse models provide an opportunity for studying mechanisms of toxicity, mouse primary hepatocyte cultures are less well described. The potential usefulness of a mouse hepatocyte-based in vitro model was assessed in this study by investigating time-dependent competence for xenobiotic metabolism and gene expression profiles. Primary mouse hepatocytes, isolated using two-step collagenase perfusion, were cultured in a collagen sandwich configuration. Gene expression profiles and the activities of various cytochrome P450 (P450) enzymes were determined after 0, 42, and 90 h in culture. Principal component analysis of gene expression profiles shows that replicates per time point are similar. Gene expression levels of most phase I biotransformation enzymes decrease to approximately 69 and 57% of the original levels at 42 and 90 h, respectively, whereas enzyme activities for most of the studied P450s decrease to 59 and 34%. The decrease for phase II gene expression is only to 96 and 92% of the original levels at 42 and 90 h, respectively. Pathway analysis reveals initial effects at the level of proteins, external signaling pathways, and energy production. Later effects are observed for transcription, translation, membranes, and cell cycle-related gene sets. These results indicate that the sandwich-cultured primary mouse hepatocyte system is robust and seems to maintain its metabolic competence better than that of the rat hepatocyte system.

Background: In gene expression analysis, statistical tests for differential gene expression provide lists of candidate genes having, individually, a sufficiently low p-value. However, the interpretation of each single p-value within complex systems involving several interacting genes is problematic. In parallel, in the last sixty years, game theory has been applied to political and social problems to assess the power of interacting agents in forcing a decision and, more recently, to represent the relevance of genes in response to certain conditions., Results: In this paper we introduce a Bootstrap procedure to test the null hypothesis that each gene has the same relevance between two conditions, where the relevance is represented by the Shapley value of a particular coalitional game defined on a microarray data-set. This method, which is called Comparative Analysis of Shapley value (shortly, CASh), is applied to data concerning the gene expression in children differentially exposed to air pollution. The results provided by CASh are compared with the results from a parametric statistical test for testing differential gene expression. Both lists of genes provided by CASh and t-test are informative enough to discriminate exposed subjects on the basis of their gene expression profiles. While many genes are selected in common by CASh and the parametric test, it turns out that the biological interpretation of the differences between these two selections is more interesting, suggesting a different interpretation of the main biological pathways in gene expression regulation for exposed individuals. A simulation study suggests that CASh offers more power than t-test for the detection of differential gene expression variability., Conclusion: CASh is successfully applied to gene expression analysis of a data-set where the joint expression behavior of genes may be critical to characterize the expression response to air pollution. We demonstrate a synergistic effect between coalitional games and statistics that resulted in a selection of genes with a potential impact in the regulation of complex pathways.

Differences in biological responses to exposure to hazardous airborne substances between children and adults have been reported, suggesting children to be more susceptible. Aim of this study was to improve our understanding of differences in susceptibility in cancer risk associated with air pollution by comparing genome-wide gene expression profiles in peripheral blood of children and their parents. Gene expression analysis was performed in blood from children and parents living in two different regions in the Czech Republic with different levels of air pollution. Data were analyzed by two different approaches: one method first selected significantly differentially expressed genes and analyzed these gene lists for overrepresented biological processes, whereas the other applied the T-profiler tool to directly perform pathway analyses on the total gene set without preselection of significantly modulated gene expressions. In addition, gene expressions in both children and adults were investigated for associations with micronuclei frequencies. Both analysis approaches returned considerably more genes or gene groups and pathways that significantly differed between children from both regions than between parents. Very little overlap was observed between children and adults. The two most important biological processes or molecular functions significantly modulated in children, but not in adults, are nucleosome and immune response related. Our study suggests differences between children and adults in relation to air pollution exposure at the transcriptome level. The findings underline the necessity of implementing environmental health policy measures specifically for protecting children's health.

Background: DNA microarrays are widely used in gene expression analyses. To increase throughput and minimize costs without reducing gene expression data obtained, we investigated whether four mRNA samples can be analyzed simultaneously by applying four different fluorescent dyes., Results: Following tests for cross-talk of fluorescence signals, Alexa 488, Alexa 594, Cyanine 3 and Cyanine 5 were selected for hybridizations. For self-hybridizations, a single RNA sample was labelled with all dyes and hybridized on commercial cDNA arrays or on in-house spotted oligonucleotide arrays. Correlation coefficients for all combinations of dyes were above 0.9 on the cDNA array. On the oligonucleotide array they were above 0.8, except combinations with Alexa 488, which were approximately 0.5. Standard deviation of expression differences for replicate spots were similar on the cDNA array for all dye combinations, but on the oligonucleotide array combinations with Alexa 488 showed a higher variation., Conclusion: In conclusion, the four dyes can be used simultaneously for gene expression experiments on the tested cDNA array, but only three dyes can be used on the tested oligonucleotide array. This was confirmed by hybridizations of control with test samples, as all combinations returned similar numbers of differentially expressed genes with comparable effects on gene expression.

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.

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. [ABSTRACT FROM AUTHOR]

Abstract: 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. [Copyright &y& Elsevier]