Edward Dere

In vivo – in vitro toxicogenomic comparison of TCDD-elicited gene expression in Hepa1c1c7 mouse hepatoma cells and C57BL/6 hepatic tissue

BackgroundIn vitro systems have inherent limitations in their ability to model whole organism gene responses, which must be identified and appropriately considered when developing predictive biomarkers of in vivo toxicity. Systematic comparison of in vitro and in vivo temporal gene expression profiles were conducted to assess the ability of Hepa1c1c7 mouse hepatoma cells to model hepatic responses in C57BL/6 mice following treatment with 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD).ResultsGene expression analysis and functional gene annotation indicate that Hepa1c1c7 cells appropriately modeled the induction of xenobiotic metabolism genes in vivo. However, responses associated with cell cycle progression and proliferation were unique to Hepa1c1c7 cells, consistent with the cell cycle arrest effects of TCDD on rapidly dividing cells. In contrast, lipid metabolism and immune responses, representative of whole organism effects in vivo, were not replicated in Hepa1c1c7 cells.ConclusionThese results identified inherent differences in TCDD-mediated gene expression responses between these models and highlighted the limitations of in vitro systems in modeling whole organism responses, and additionally identified potential predictive biomarkers of toxicity.

Comparative Analysis of AhR-Mediated TCDD-Elicited Gene Expression in Human Liver Adult Stem Cells

Time course and dose-response studies were conducted in HL1-1 cells, a human liver cell line with stem cell-like characteristics, to assess the differential gene expression elicited by 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) compared with other established models. Cells were treated with 0.001, 0.01, 0.1, 1, 10, or 100nM TCDD or dimethyl sulfoxide vehicle control for 12 h for the dose-response study, or with 10nM TCDD or vehicle for 1, 2, 4, 8, 12, 24, or 48 h for the time course study. Elicited changes were monitored using a human cDNA microarray with 6995 represented genes. Empirical Bayes analysis identified 144 genes differentially expressed at one or more time points following treatment. Most genes exhibited dose-dependent responses including CYP1A1, CYP1B1, ALDH1A3, and SLC7A5 genes. Comparative analysis of HL1-1 differential gene expression to human HepG2 data identified 74 genes with comparable temporal expression profiles including 12 putative primary responses. HL1-1-specific changes were related to lipid metabolism and immune responses, consistent with effects elicited in vivo. Furthermore, comparative analysis of HL1-1 cells with mouse Hepa1c1c7 hepatoma cell lines and C57BL/6 hepatic tissue identified 18 and 32 commonly regulated orthologous genes, respectively, with functions associated with signal transduction, transcriptional regulation, metabolism and transport. Although some common pathways are affected, the results suggest that TCDD elicits species- and model-specific gene expression profiles.

Differences in TCDD-elicited gene expression profiles in human HepG2, mouse Hepa1c1c7 and rat H4IIE hepatoma cells

Background2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) is an environmental contaminant that elicits a broad spectrum of toxic effects in a species-specific manner. Current risk assessment practices routinely extrapolate results from in vivo and in vitro rodent models to assess human risk. In order to further investigate the species-specific responses elicited by TCDD, temporal gene expression responses in human HepG2, mouse Hepa1c1c7 and rat H4IIE cells were compared.ResultsMicroarray analysis identified a core set of conserved gene expression responses across species consistent with the role of AhR in mediating adaptive metabolic responses. However, significant species-specific as well as species-divergent responses were identified. Computational analysis of the regulatory regions of species-specific and -divergent responses suggests that dioxin response elements (DREs) are involved. These results are consistent with in vivo rat vs. mouse species-specific differential gene expression, and more comprehensive comparative DRE searches.ConclusionsComparative analysis of human HepG2, mouse Hepa1c1c7 and rat H4IIE TCDD-elicited gene expression responses is consistent with in vivo rat-mouse comparative gene expression studies, and more comprehensive comparative DRE searches, suggesting that AhR-mediated gene expression is species-specific.

Comparative Analysis of Temporal and Dose-Dependent TCDD-Elicited Gene Expression in Human, Mouse, and Rat Primary Hepatocytes

2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD)-elicited time- and dose-dependent differential gene expression was compared in human, mouse, and rat primary hepatocytes. Comprehensive time course (10 nM TCDD or dimethyl sulfoxide vehicle control for 1, 2, 4, 8, 12, 24, and 48h) studies identified 495, 2305, and 711 differentially expressed orthologous genes in human, mouse, and rat hepatocytes, respectively. However, only 16 orthologs were differentially expressed across all three species, with the majority of orthologs exhibiting species-specific expression (399 human, 2097 mouse, and 533 rat), consistent with species-specific expression reported in other in vitro and in vivo comparative studies. TCDD also elicited the dose-dependent induction of 397 human, 100 mouse, and 443 rat genes at 12h and 615 human, 426 mouse, and 314 rat genes at 24h. Comparable EC50 values were obtained for AhR battery genes including Cyp1a1 (0.1 nM human, 0.05 nM mouse, 0.08 nM rat at 24h) and Tiparp (0.97 nM human, 0.63 nM mouse, 0.14 nM rat at 12h). Overrepresented functions and pathways included amino acid metabolism in humans, immune response in mice, and energy homeostasis in rats. Differentially expressed genes functionally associated with lipid transport, processing, and metabolism were overrepresented in all three species but exhibited species-specific expression consistent with the induction of hepatic steatosis in mice but not in rats following a single oral gavage of TCDD. Furthermore, human primary hepatocytes showed lipid accumulation following 48h of treatment with TCDD, suggesting that AhR-mediated steatosis in mice more closely resembles human hepatic fat accumulation compared with that in rats. Collectively, these results suggest that species-specific gene expression profiles mediate the species-specific effects of TCDD despite the conservation of the AhR and its signaling mechanism.

Identification of aryl hydrocarbon receptor binding targets in mouse hepatic tissue treated with 2,3,7,8-tetrachlorodibenzo-p-dioxin

Genome-wide, promoter-focused ChIP-chip analysis of hepatic aryl hydrocarbon receptor (AHR) binding sites was conducted in 8-week old female C57BL/6 treated with 30 μg/kg/body weight 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) for 2 h and 24 h. These studies identified 1642 and 508 AHR-bound regions at 2h and 24h, respectively. A total of 430 AHR-bound regions were common between the two time points, corresponding to 403 unique genes. Comparison with previous AHR ChIP-chip studies in mouse hepatoma cells revealed that only 62 of the putative target genes overlapped with the 2 h AHR-bound regions in vivo. Transcription factor binding site analysis revealed an over-representation of aryl hydrocarbon response elements (AHREs) in AHR-bound regions with 53% (2 h) and 68% (24 h) of them containing at least one AHRE. In addition to AHREs, E2f-Myc activator motifs previously implicated in AHR function, as well as a number of other motifs, including Sp1, nuclear receptor subfamily 2 factor, and early growth response factor motifs were also identified. Expression microarray studies identified 133 unique genes differentially regulated after 4 h treatment with TCDD. Of which, 39 were identified as AHR-bound genes at 2 h. Ingenuity Pathway Analysis on the 39 AHR-bound TCDD responsive genes identified potential perturbation in biological processes such as lipid metabolism, drug metabolism, and endocrine system development as a result of TCDD-mediated AHR activation. Our findings identify direct AHR target genes in vivo, highlight in vitro and in vivo differences in AHR signaling and show that AHR recruitment does not necessarily result in changes in target gene expression.

Comparative Temporal Toxicogenomic Analysis of TCDD- and TCDF-Mediated Hepatic Effects in Immature Female C57BL/6 Mice

Temporal analyses were performed on hepatic tissue from immature female C57BL/6 mice in order to compare the gene expression profiles for 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) and 2,3,7,8-tetrachlorodibzofuran (TCDF). Time course studies conducted with a single oral dose of 300 microg/kg TCDF or 30 microg/kg TCDD were used to compare differential gene expression on complementary DNA microarrays containing 13,361 features, representing 8194 genes at 2, 4, 8, 12, 24, 72, 120, and 168 h. One hundred and ninety-five genes were identified as differentially regulated by TCDF, of which 116 genes were in common with TCDD, with 109 exhibiting comparable expression profiles (correlation coefficients > 0.3). In general, TCDF was less effective in eliciting hepatic vacuolization, and differential gene expression compared with TCDD when given at an equipotent dose based on a toxic equivalence factor (TEF) of 0.1 for TCDF, especially 72-h postadministration. For example, the induction of Cyp1a1 messenger RNA by TCDF was less when compared TCDD. Moreover, TCDF induced less severe hepatocyte cytoplasmic vacuolization consistent with lower lipid accumulations which significantly subsided by 120 and 168 h when compared with TCDD. TCDF-elicited responses correlated with their hepatic tissue levels which gradually decreased between 18 and 168 h. Although both compounds elicited comparable gene expression profiles, especially at early time points, the TCDF responses were generally weaker. Collectively, these results suggest that the weaker TCDF responses could be attributed to differences in pharmacokinetics. However, more comprehensive dose-response studies are required at optimal times for each end point of interest in order to investigate the effect of pharmacokinetic differences on relative potencies that are important in establishing TEFs.

Gene Coregulation and Coexpression in the Aryl Hydrocarbon Receptor-mediated Transcriptional Regulatory Network in the Mouse Liver

Tissue-specific network models of chemical-induced gene perturbation can improve our mechanistic understanding of the intracellular events leading to adverse health effects resulting from chemical exposure. The aryl hydrocarbon receptor (AHR) is a ligand-inducible transcription factor (TF) that activates a battery of genes and produces a variety of species-specific adverse effects in response to the potent and persistent environmental contaminant 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). Here we assemble a global map of the AHR gene regulatory network in TCDD-treated mouse liver from a combination of previously published gene expression and genome-wide TF binding data sets. Using Kohonen selforganizing maps and subspace clustering, we show that genes co-regulated by common upstream TFs in the AHR network exhibit a pattern of co-expression. Specifically, directly-bound, indirectly-bound and non-genomic AHR target genes exhibit distinct patterns of gene expression, with the directly bound targets generally associated with highest median expression. Further, among the directly bound AHR target genes, the expression level increases with the number of AHR binding sites in the proximal promoter regions. Finally, we show that co-regulated genes in the AHR network activate distinct groups of downstream biological processes, with AHR-bound target genes enriched for metabolic processes and enrichment of immune responses among AHR-unbound target genes, likely reflecting infiltration of immune cells into the mouse liver upon TCDD treatment. This work describes an approach to the reconstruction and analysis of transcriptional regulatory cascades underlying cellular stress response using bioinformatic and statistical tools.