Marcel van Herwijnen

Comparison of HepG2 and HepaRG by Whole-Genome Gene Expression Analysis for the Purpose of Chemical Hazard Identification

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.

In vitro effects of coal fly ashes: hydroxyl radical generation, iron release, and DNA damage and toxicity in rat lung epithelial cells

Oxygen radical generation due to surface radicals, inflammation, and iron release has been suggested as the mechanism of adverse effects of quartz, such as emphysema, fibrosis, and carcinogenic effects. Therefore, we measured iron release, acellular generation of hydroxyl radicals, and oxidative DNA damage and cytotoxicity in rat lung epithelial (RLE) cells by different coal fly ashes (CFA) that contain both quartz and iron. Seven samples of CFA with different particle size and quartz content (up to 14.1%) were tested along with silica (alpha-quartz), ground coal, and coal mine dust (respirable) as positive control particles, and fine TiO(2) (anatase) as a negative control. Five test samples were pulverized fuel ashes (PFA), two samples were coal gasification (SCG) ashes (quartz content <0.1%), and one sample was a ground coal. No marked differences between SCG and PFA fly ashes were observed, and toxicity did not correlate with physicochemical characteristics or effect parameters. Stable surface radicals were only detected in the reference particles silica and coal mine dust, but not in CFA. On the other hand, hydroxyl radical generation by all fly ashes was observed in the presence of hydrogen peroxide, which was positively correlated with iron mobilization and inhibited by deferoxamine, but not correlated with iron or quartz content. Also a relationship between acellular hydroxyl radical generation and oxidative DNA damage in RLE cells by CFA was observed. Differences in hydroxyl radical generation and oxidative damage by the CFA were not related to iron and quartz content, but the respirable ashes (MAT023, 38, and 41) showed a very extensive level of hydroxyl radical generation in comparison to nonrespirable fly ashes and respirable references. This radical generation was clearly related to the iron mobilization from these particles. In conclusion, the mechanisms by which CFA and the positive references (silica, coal mine dust) affect rat lung epithelial cells seem to be different, and the data suggest that quartz in CFA does not act the same as quartz in silica or coal mine dust. On the other hand, the results indicate an important role for size and iron release in generation and subsequent effects of reactive oxygen species caused by CFA.

Discrimination for genotoxic and nongenotoxic carcinogens by gene expression profiling in primary mouse hepatocytes improves with exposure time.

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.

Genomic analysis suggests higher susceptibility of children to air pollution

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

Parallelogram approach using rat-human in vitro and rat in vivo toxicogenomics predicts acetaminophen-induced hepatotoxicity in humans.

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.

Embryotoxicant-specific transcriptomic responses in rat postimplantation whole-embryo culture.

Rat postimplantation whole-embryo culture (WEC) is a promising alternative test for the assessment of developmental toxicity. Toxicogenomic-based approaches may improve the predictive ability of the WEC model by providing a means to identify compound-specific mechanistic responses associated with embryotoxicity in vivo. Furthermore, alterations in gene expression may serve as a sensitive, objective, and robust marker, which precedes the observation of classical developmental toxicity endpoints in time. In this study, in combination with morphological developmental assessments, we studied transcriptomic responses associated with four distinct teratogens (caffeine [CAF], methylmercury [MM], monobutyl phthalate, and methoxyacetic acid) after 4 h of exposure, well before apparent embryotoxicity in WEC. We evaluated gene expression changes associated with similar levels of induced morphological embryotoxicity for each teratogen (as determined by total morphological score), evaluating for functional enrichment and quantitative changes in response. Concentrations selected for each of the four teratogens used induced a number of common effects on embryonic development (neural tube closure and optic/otic system). Despite inducing common morphological effects, our analysis suggests limited overlap in terms of toxicogenomic response at the gene expression level and at the level of biological processes across all four test chemicals. Many unique responses associated with each chemical correlated with previously hypothesized modes of developmental toxicity. For example, alterations in developmental signaling and cholesterol metabolism were observed with MM and CAF, respectively. This initial study suggests that distinct chemically induced toxicogenomic responses precede morphological effects in WEC and that these responses are relevant with mechanisms of toxicity previously observed in vivo.

Time Series Analysis of Benzo[A]Pyrene-Induced Transcriptome Changes Suggests That a Network of Transcription Factors Regulates the Effects on Functional Gene Sets

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.

Classification of hepatotoxicants using HepG2 cells: A proof of principle study.

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.

A toxicogenomics-based parallelogram approach to evaluate the relevance of coumarin-induced responses in primary human hepatocytes in vitro for humans in vivo.

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.

Global Gene Expression Analysis Reveals Differences in Cellular Responses to Hydroxyl- and Superoxide Anion Radical-Induced Oxidative Stress in Caco-2 Cells

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.