- Research article
- Open Access
Aryl hydrocarbon receptor (AHR)-regulated transcriptomic changes in rats sensitive or resistant to major dioxin toxicities
BMC Genomicsvolume 11, Article number: 263 (2010)
The major toxic effects of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) appear to result from dysregulation of mRNA levels mediated by the aryl hydrocarbon receptor (AHR). Dioxin-like chemicals alter expression of numerous genes in liver, but it remains unknown which lie in pathways leading to major toxicities such as hepatotoxicity, wasting and lethality. To identify genes involved in these responses we exploited a rat genetic model. Rats expressing an AHR splice-variant lacking a portion of the transactivation domain are highly resistant to dioxin-induced toxicities. We examined changes in hepatic mRNA abundances 19 hours after TCDD treatment in two dioxin-resistant rat strains/lines and two dioxin-sensitive rat strains/lines.
Resistant rat strains/lines exhibited fewer transcriptional changes in response to TCDD than did rats with wildtype AHR. However, well-known AHR-regulated and dioxin-inducible genes such as CYP1A1, CYP1A2, and CYP1B1 remained fully responsive to TCDD in all strains/lines. Pathway analysis indicated that the genes which respond differently to TCDD between sensitive and resistant rats are mainly involved in lipid metabolism, cellular membrane function and energy metabolism. These pathways previously have been shown to respond differently to dioxin treatment in dioxin-sensitive versus dioxin-resistant rats at a biochemical level and in the differential phenotype of toxicologic responses.
The transactivation-domain deletion in dioxin-resistant rats does not abolish global AHR transactivational activity but selectively interferes with expression of subsets of genes that are candidates to mediate or protect from major dioxin toxicities such as hepatotoxicity, wasting and death.
Dioxin-like chemicals are exceptionally toxic to a wide variety of birds, fish and mammals including, perhaps, humans. However, susceptibility to dioxin toxicity varies widely among different animal species and between genetic types within a species. Extensive evidence demonstrates that virtually all toxic effects of TCDD and related dioxin-like compounds are mediated by a ligand-dependent transcription factor, the aryl hydrocarbon receptor (AHR) (reviewed in: ). Exposure to dioxins results in major toxicities, including thymic atrophy, teratogenesis, hepatotoxicity, wasting syndrome and death. These toxicities are dependent on both the AHR and its dimerization partner, the aryl hydrocarbon receptor nuclear translocator (ARNT), and require that the AHR have functional structures for nuclear translocation and DNA binding [2–6].
Dioxin binding converts the AHR into an activated ligand:AHR:ARNT complex that regulates transcription either by binding directly to AHRE-I motifs (also known as DREs or XREs) [7–10] or indirectly to AHRE-II motifs [11, 12]. Dioxin toxicities appear to arise from AHR-mediated dysregulation of specific genes [4, 13]. Microarray technologies have accelerated identification of genes that depend on the AHR for constitutive expression or for response to TCDD in vivo[14–23], but the key genes whose dysregulation by dioxin leads to most toxicities remain unknown.
The Han/Wistar(Kuopio) (H/W) rat is an excellent model organism to identify specific AHR-regulated genes whose dysregulation by dioxin may lead to major toxicities. H/W rats are extraordinarily resistant to acute lethality from TCDD, with an LD50 three orders of magnitude higher than for sensitive Long-Evans(Turku/AB) rats (L-E) . Resistance in H/W rats is associated with a point mutation that leads to expression of an aberrant AHR protein missing 38 or 43 amino acids from its transactivation domain (TAD) [25–27]. We recently demonstrated, using transgenic mouse models, that it is the AHR variant with 38 amino acids deleted that is responsible for dioxin resistance . Importantly, this AHR variant is the predominantly-expressed form in the dioxin-resistant H/W rat . Further, dioxin resistance in rats segregates genetically with the AHR locus and is a dominant trait [24, 26]. Multi-generational crosses of L-E and H/W rats and selection for susceptibility or resistance to dioxin lethality  produced two rat lines: Line-A (LnA) and Line-C (LnC). LnA rats harbour the variant AHR and are dioxin-resistant; LnC rats have the wildtype AHR and are dioxin-sensitive. We postulate that the partial deletion of AHR transactivation domain alters toxic responses either by preventing changes in mRNA levels of genes in pro-death pathways or by enhancing responses of genes in pro-survival pathways.
Our strategy for identifying pro-survival or pro-death genes is to contrast changes in mRNA expression profiles following TCDD exposure of dioxin-sensitive and dioxin-resistant rats . To reduce the influence of strain-specific changes not associated with dioxin toxicities, we profiled 4 strains/lines of rats: H/W and LnA which constitute the "resistant collective" and L-E and LnC which constitute the "sensitive collective".
We focused on hepatic mRNA levels because liver displays a broad spectrum of mRNAs that are responsive to dioxins and/or to AHR genotype  and because liver is a prime site of dioxin toxicity, displaying many phenotypic differences between sensitive and resistant rats . We chose a dose of 100 μg/kg TCDD, which produces hepatotoxicity, wasting and death in sensitive rats but no deaths in resistant rats. We previously conducted a smaller-scale transcriptomic study in sensitive versus resistant rats on membrane arrays and using cross-species hybridization to cDNA arrays . Here, we greatly extend those prior studies by assessing transcriptome-wide responses to TCDD which were further validated via real-time RT-PCR. We identify specific biological processes perturbed by TCDD exposure.
Our analysis paints a new picture of dioxin-induced expression changes. Hundreds of genes exhibit responses to TCDD that are specific to individual strains or lines. Our genetic model diminishes this background noise and identifies a small number of genes associated with hepatotoxicity, wasting and death. Genes differentially-expressed between sensitive and resistant strains show functional homogeneity: dioxin-lethality may be associated with broad dysregulation of entire pathways, not just single genes.
To determine the effect of TCDD on mRNA abundances in dioxin-sensitive rats versus dioxin-resistant rats we studied four rat strains/lines at 19 hours after oral administration of TCDD (Figure 1) using Affymetrix RAE230A microarrays. A list of all genes and their responses to TCDD in the four rat strains/lines is given in Additional File 1.
Global differences in mRNA abundance
The total number of genes affected by TCDD varied across rat strains/lines in a manner independent of the statistical threshold (Figure 2A). More genes were affected by TCDD in dioxin-sensitive rats than in dioxin-resistant rats: across all four strains/lines, the order of the number of transcriptional alterations was: L-E > LnC > LnA > H/W. Overall, it appears that the deletion in the AHR transactivation domain reduces the number of genes altered but does not ablate the transcriptional response to TCDD.
Unsupervised clustering analysis (Figure 2B) using all non-constant genes (variance > 0.01) revealed the strongest trend in the dataset was the distinction between TCDD-treated (red annotation bars) and vehicle-treated (white annotation bars) animals. The co-clustering of the two F2 crosses (LnA, dark blue; LnC, light blue) may indicate that their common parentage is a stronger determinant of their transcriptional profiles than is their sensitivity or resistance to TCDD toxicity. Further, this suggests that number of genes causally related to the hepatotoxicity (and possibly acute lethality) of TCDD is small compared with all other changes caused by TCDD in hepatic gene expression.
The extent of the overlap in transcriptional responses to TCDD among sensitive and resistant rat strains/lines was visualized using a two-way table (Table 1) and Venn diagrams (Figure 3). Of the 8605 genes interrogated, 452 (5.3%) responded in at least one strain or line at a 1% false-discovery rate; 8153 genes did not respond to TCDD in any strain/line (Table 1). Of the 452 responsive genes only 144 (31.9%) were altered in more than one strain/line. Only 25 genes responded in all strains/lines, of which 20 were induced while only 5 were genes repressed. This direction of response is concordant with our previous finding that 70% of genes altered in common by TCDD in both mouse and rat are up-regulated .
The sensitive collective exhibited a greater number of responsive genes (412) and a greater overlap of these genes among strains/lines within the collective (Figure 3A) than did the resistant collective (138 responsive genes; Figure 3B). However, the proportion of TCDD-responsive genes that overlapped within each collective to the total genes on the array that responded to TCDD in that collective did not differ between collectives (22.5% in the resistant vs. 26.2% in the sensitive collective, p = 0.44; proportion test).
Classification of Type-I vs. Type-II responses to TCDD
Identifying which of the 452 TCDD-responsive genes are most likely to be involved in major dioxin toxicities is challenging; we therefore exploited our genetic model. Each TCDD-responsive gene was classified according to the Type-I/Type-II TCDD response scheme previously developed for toxic endpoints .
Type-I responses to TCDD are those that are similar between dioxin-sensitive rat strains/lines and dioxin-resistant rat strains/lines. We considered genes that exhibited a statistically significant response to TCDD in all four strains/lines to be Type-I genes. By this definition 25 genes were classified as Type-I (Figure 3C; Score ± 4 in Table 1). These Type-I genes include well-known (Table 2) as well as novel TCDD-responsive genes (Additional File 1). The vast majority of Type-I genes (20/25) were up-regulated by TCDD (score +4 in Table 1). For some Type-I genes the magnitude of induction was very large, including ~90-fold induction of CYP1A1 mRNA and ~75-fold induction of ALDH3A1 mRNA.
Type-II responses to TCDD are those that differ between dioxin-sensitive rats and dioxin-resistant rats. Genes that responded to TCDD exposure in a statistically significant manner in both members of one collective but neither of the strains/lines in the other collective are classified as Type-II genes. By this classification, 46 genes exhibited Type-II responses (Figure 3C; Score ± 2 & 0 in both ways in Table 1). Of these, only three were altered in resistant rats but not in sensitive rats (Il1r1, Phyh, Hacl1). The remaining 43 genes were specifically altered only in dioxin-sensitive rats. In contrast to the 80% upregulation of Type-I genes, 61% (28/46) of Type-II genes were downregulated by TCDD.
Validation by real-time RT-PCR
The validity of our array experiments to identify effects of TCDD on mRNA levels is supported by the facts that well-established dioxin-inducible genes were identified (e.g. CYP1A1, CYP1A2, CYP1B1 and Tiparp; Table 2) and there is overlap between our current list of responsive genes and genes previously reported to be affected by AHR-ligands in other array experiments [15, 17].
To further assess the validity of our array results we employed real-time RT-PCR to evaluate effects of TCDD on mRNA levels. Genes were selected for RT-PCR analysis to span a wide range of magnitudes of response (fold-change in log2: high > 4; medium 2-4; low < 2; or no response). To this end, both array (Table 2) and RT-PCR analyses (Figure 4) showed a high magnitude of TCDD-mediated induction of the prototypical responsive/AHR-activated gene CYP1A1[13, 32, 33]. Further, CYP7A1 and Chka exhibited medium magnitudes of response to TCDD in both collectives as evidenced by RT-PCR assays (Figure 4) as well as, array experiments (Additional File 1). Selenbp1 and Per2 showed low magnitudes of TCDD-mediated induction in both collectives as evidenced by array experiments (Additional File 1) and slightly higher magnitudes of induction by RT-PCR assays (Figure 4). Elov6 exhibited significant induction in the dioxin-sensitive collective but not in the resistant collective as assayed both by array (Table 3) and RT-PCR (Figure 4). TCDD had no significant effect on mRNA levels for Klf10 or Pik3r1 as assayed either by RT-PCR (Figure 4) or by gene array (Additional File 1).
Analysis of the AHR role in regulation
It is well-established that major toxicities of TCDD require the AHR [2, 4, 5, 34, 35]. To determine if the AHR was required for the gene to respond to TCDD we compared hepatic mRNA levels for a few genes in Ahr-null mice (Ahr-/- ) versus mice with wildtype AH receptor (Ahr+/+ ). Mice were treated with a dose of TCDD equitoxic to that in sensitive rats  for a comparable time. As expected, induction of CYP1A1 mRNA was strictly dependent on the AHR (Figure 5). Suppression of Crip2 mRNA levels also was AHR-dependent (Figure 5). Regulation of Chka and Elovl6 by TCDD appears to be species-specific: Chka mRNA was upregulated in rats from both collectives (Additional File 1 and Figure 4) but was significantly downregulated in Ahr-null mice and unaffected in wildtype mice (Figure 5); Elovl6 was upregulated in dioxin-sensitive rats (Table 3 and Figure 4) but Elovl6 was unresponsive to TCDD in both wildtype and Ahr-null mice (Figure 5). The findings with Chka and Elovl6 reinforce recent reports of substantial differences between rat and mouse in transcriptional responses to TCDD [16, 23].
To determine if alterations in mRNA abundances caused by TCDD are functionally coherent, we performed two Gene Ontology (GO) analyses. The first used candidate genes from each rat strain/line to determine the extent of overlap of pathways between rat strains/lines. We found that the pathways dysregulated by TCDD were very similar in the two dioxin-resistant rats, LnA and H/W (Figure 6A) and in the two dioxin-sensitive rats, LnC and L-E (Figure 6B). Further, there was overlap in pathways among all strains/lines: 8 GO terms were enriched in all four strains/lines, while 14 GO terms were specifically enriched in only the sensitive strains (Figure 6C).
The second analysis used both Type-I and Type-II gene lists to determine the specific GO terms enriched in each of the Type-I and Type-II gene lists. This analysis identified 5 GO terms - mostly relating to cytochrome P450 genes - enriched in the Type-I genes (Table 4). Amongst the Type-II genes, the analysis identified 14 GO terms specifically enriched (Table 4). In particular, genes related to the endoplasmic reticulum were present about 4 times as often as expected by chance alone, as were genes related to lipid metabolism.
We performed transcriptional profiling on livers of rats that are sensitive or resistant to major TCDD toxicities. Two key findings arise. First, we show significant inter-strain and inter-species diversity in responses to TCDD. Second, we identify Type-II genes that may be integral to the mechanism(s) of hepatotoxicity, wasting and lethality.
Significant diversity of intra-species and inter-species responses to TCDD
One startling characteristic of the transcriptional response to TCDD across the four rat strains/lines are the dramatic inter-strain differences. Given their substantial genetic relatedness, including at the AHR locus, it might be hypothesized that LnA and H/W animals would have very similar responses, and that LnC would be very similar to L-E. To the contrary, 68.1% (308 of 452) of dioxin-responsive genes were altered in one of the four rat strains/lines. These results clearly demonstrate the importance of genomic context in regulating mRNA responses to dioxin-exposure and mirror an analysis of the basal mRNA levels in these and other rat strains (Boutros et al. submitted). Interestingly, in rats with the AHRH/W genotype, the total number of genes that respond was reduced relative to rats expressing wildtype AHR. However, AHRs from rats with the AHRH/W genotype have similar affinity for TCDD and ability to bind AH response elements as wildtype rats . It is conceivable that AHRH/W rats have a reduced ability to recruit coactivators and interact with transcriptional machinery. Chromatin immuno-precipitation experiments would be valuable in testing this hypothesis directly.
When we attempted to study AHR-dependency of the mRNA changes by comparison with Ahr-/- mice we found that only 2 of 4 mRNA responses to TCDD, measured by RT-PCR, could be compared in both rat and mouse models. This result concords with recent reports of highly divergent transcriptomic responses to TCDD between rat and mouse , and suggests that combining our intra-species rat model with inter-species studies may be a fruitful approach for identifying genes that mediate TCDD-induced toxicities, especially those toxic responses that differ between animal species.
Type-II responsive genes whose regulation by TCDD may be integral to the mechanism(s) of hepatotoxicity, wasting and lethality
In dioxin-sensitive rats, a wasting syndrome commences within the first few days following a single dose of TCDD and is characterized by progressive weight loss (eventually up to 50%) and hypophagia [37–39]. Wasting contributes to lethality starting 2-3 weeks after TCDD exposure. Although wasting accompanies death, wasting per se is not likely to be the sole reason for death since maintenance of body weight by parenteral nutrition does not prevent mortality . Nevertheless, untreated control rats, pair-fed at the same caloric intake as rats treated with lethal TCDD doses, die at much the same time as their TCDD-exposed partners .
The cause of TCDD-induced wasting and ultimate death as well as the key target tissue(s) remain elusive. Since TCDD causes extensive hepatotoxicity in TCDD-sensitive rats (but not in the TCDD-resistant strains), it is reasonable that examination of mechanisms which produce hepatotoxicity may provide clues to mechanisms of wasting and lethality. To this end, we identified 46 Type-II hepatic genes whose TCDD-responsiveness differed between TCDD-resistant and TCDD-sensitive rats. Pathway analysis indicated that these genes are mainly involved in lipid-metabolism, cellular membrane function and energy metabolism (Table 4). These Type-II genes potentially explain why there are greater manifestations of hepatotoxicity in sensitive rats than in resistant rats; for example, a dramatic accumulation of fatty acids (steatosis) and initial liver hypertrophy which switches to atrophy ~1 week later (refer to  for exhaustive list of hepatotoxic responses). Pohjanvirta et al. , in biochemical studies, found that exposure of sensitive rats to TCDD led to steatosis, hypertrophy, liver failure, wasting and eventual death, possibly as the consequence of derailment of energy metabolism due to alterations of (i) lipid homeostasis, (ii) protein metabolism and/or (iii) ATP production/utilization.
(i) Alteration of lipid homeostasis
In sensitive rats only, we previously observed steatosis with the accumulated fatty acids probably originating from redistribution of peripheral fat deposits to liver rather than from increased lipid synthesis within the liver . This is consistent with our current transcriptomic study in which transcripts related to lipid synthesis were not increased. TCDD-induced suppression of hepatic lipid lipogenesis previously has been reported [16, 17, 43].
Pathway analysis indicated that genes involved in the lipid metabolism process were enriched in the Type-II gene list. Of particular interest was the decreased expression of Hsd11b1 and Slc27a5 only in sensitive rats. Hsd11b1 functions in steroid metabolism and colocalizes with the glucocorticoid receptor where it acts as a local amplifier of corticoid responses including the regulation of fuel metabolism during starvation and stress [44, 45]. Hsd11b1 deficiencies in rodents increase energy expenditure, decrease weight gain with chronic high fat feeding, increase weight loss, increase hepatic lipid oxidation while decreasing lipolysis in adipose tissue and display many metabolic deficiencies [45, 46]. Slc27a5 encodes a transporter of long-chain fatty acids into the liver where it is exclusively expressed. Its deletion results in increased de novo biosynthesis of long-chain fatty acids in liver due to inhibited uptake of them. Interestingly, in knockout mice, feed intake is depressed, energy expenditure increased and weight gain suppressed [47, 48]. While further study of these genes is warranted, genes that alter lipid homeostasis may be important in hepatotoxicity and could be involved in pro-death pathways in sensitive rats exposed to TCDD.
(ii) Altered protein metabolism
In the short term, protein catabolism can be a beneficial response to provide amino acids for energy and maintenance of obligatory functions. However, sustained protein catabolism eventually leads to wasting and mortality. Our current study found that TCDD increased expression of genes that facilitate protein breakdown (Derl1, Derl2 and Mug1) but also decreased the expression of Ass1, the key enzyme in the urea cycle. Deficiencies in Ass1 may disrupt the urea cycle resulting in increased accumulation of amino acids, highly toxic ammonia and other toxic byproducts. TCDD previously has been shown to decrease the expression of the Ass1 gene after 24 hour exposure . Increased protein breakdown and deficiencies in the urea cycle are consistent with previous reports of elevated plasma levels of most amino acids and decreased plasma urea in sensitive rats but not in resistant rats 6 days after TCDD exposure . In addition, alteration of the balance between protein synthesis versus protein degradation is a key mechanism in switching hepatocytes from hypertrophy to atrophy, as observed in sensitive rats ~1 week after TCDD exposure  and reported for other wasting diseases (e.g. diabetes and cancer cachexia) . Thus, specific genes involved in protein homeostasis are likely important in TCDD-mediated pro-survival pathways.
(iii) Impaired ATP production/utilization
Derailment of energy metabolism due to impaired ATP production or utilization potentially contributes to manifestations of major TCDD toxicities. After TCDD exposure in sensitive rats, compensatory mechanisms may attempt to increase energy for metabolism by increasing expression of the Atp5c1 gene that resides within the ATP synthesis pathway and is down-regulated in obese subjects . However, TCDD decreased expression levels of Adk which catalyzes the inter-conversion of adenine nucleotides, and plays an important role in cellular energy homeostasis (2 ADP ↔ ATP + AMP); this downregulation may impair use of ATP as an energy source in sensitive rats. Specifically, Adk impairment leads to deficiencies in adenosine nucleotides, including ATP, likely leading to reduced mitochondrial metabolic capacity and impairment of lipid metabolism critical for energy production . Moreover, Adk-deficient mice display hepatic steatosis within 4 days and die within 14 days with fatty liver . Thus, deficiencies in adenosine metabolism are powerful contributors to development of hepatic steatosis and development of lethal fatty liver, processes that are also provoked by TCDD in sensitive rats.
In dioxin-resistant rats, where TCDD does not cause severe hepatotoxicity, wasting or death, there was no alteration of mRNA levels for genes which might derail energy metabolism due to alterations of in the homeostasis of lipids, protein metabolism or ATP production/utilization. Only three genes responded to TCDD in the resistant collective but not in the sensitive collective. Two of these three genes function in lipid metabolism: the expression of Phyh was increased while that of Hacl1 was decreased. The third gene Il1r1, whose levels were increased by TCDD, is a receptor whose responsiveness regulates several biological functions, including adaptive and innate immunity, control of programmed cell death and stress response . The consequences of altered regulation of these genes in mediating potential pro-survival pathways in response to TCDD warrants further investigation.
The mechanisms of dioxin-induced toxicities remain elusive but our transcriptomic approach in an in vivo rat model where there are major phenotypic differences in the toxic response is providing clues to the early events that may trigger toxicity. Compelling evidence shows that the transcriptional activity of the AHR is essential for toxicity. Because altered transcription is central to TCDD toxicity, our group and others have profiled changes of mRNA abundance resulting from exposure to TCDD in several model systems. The results are remarkable: TCDD induces wide-spread alterations in mRNA abundance, but only a very small fraction of these changes are conserved between mouse or rat and, as demonstrated here, within different rat strains. Only 31.9% of dioxin-responsive genes are altered in more than one of the four rat strains/lines.
This diversity of transcriptional responses makes it challenging to identify specific genes responsible for lethality and other major forms of dioxin toxicity. Our results suggest that hepatic toxicity probably is not caused by dysregulation of a single critical gene. Rather, pathways such as lipid metabolism or energy metabolism may be derailed by altered transcription of multiple genes, possibly under coordinate control of the AHR with participation of other regulatory factors. Pathways and individual genes highlighted here are worthy candidates for further mechanistic studies to test their role in mediating or protecting from major dioxin toxicities.
Animals and Treatments for Rat Model
We studied two dioxin-sensitive rat strains/lines expressing wildtype AHR: Long-Evans (Turku/AB) (L-E) and Line-C (LnC). We also studied two dioxin-resistant rat strains/lines expressing the Han/Wistar variant AHR: Han/Wistar (Kuopio) (H/W) and Line-A (LnA) . All animals were males 10-12 weeks of age from breeding colonies of the National Institute for Health and Welfare, Kuopio, Finland. They were housed in groups of 4 (an entire treatment group per cage) in suspended stainless-steel wire-mesh cages with pelleted R36 feed (Lactamin, Stockholm, Sweden) and tap water available ad libitum. The temperature in the animal room was 21 ± 1°C, relative humidity 50 ± 10%, and a 12 hour-light/12 hour-dark cycle. Study plans were approved by the Animal Experiment Committee of the University of Kuopio and the Provincial Government of Eastern Finland. There were four rats per treatment group. Liver was harvested between 8:30 and 11:00 from rats treated by gavage with a single 100 μg/kg dose of TCDD or corn oil vehicle 19 hours previously. The single dose of 100 μg/kg TCDD produces hepatotoxicity, wasting and death in sensitive rats but not in resistant rats.
Animals and Treatment for Mouse Model
Liver tissues were from mice in which we previously mapped AHR-dependent and dioxin-dependent gene batteries by transcriptomic analysis . Briefly, male Ahr-null (Ahr-/- ) mice in a C57BL/6J background (10 weeks old) and C57BL/6 mice carrying wildtype (Ahr+/+ ) (15 weeks old) were given a single dose of 1000 μg/kg TCDD or corn oil vehicle by gavage. Liver was harvested 19 hours after treatment. The single dose of TCDD is lethal to wild-type but not Ahr-null mice and is equitoxic to that given to sensitive rats. We tested 3 TCDD-treated and 3 control mice in the Ahr-/- groups and 4 TCDD-treated and 4 control mice in the Ahr+/+groups.
Total RNA was extracted using Qiagen RNeasy kits according to the manufacturer's instructions (Qiagen, Mississauga, Canada). Total RNA yield was quantified by UV spectrophotometry and RNA integrity was verified using an Agilent 2100 BioAnalyzer (Agilent Technologies, Santa Clara, CA).
Sample labeling and hybridization to Affymetrix RAE230A GeneChips® were performed by The Centre for Applied Genomics (Toronto, Canada) according to the manufacturer's protocols. At each condition four separate animals were profiled, each on an individual RAE230A microarray. Raw array data were examined for spatial and distributional heterogeneity and differential RNA degradation; no arrays were excluded. Array data were loaded into the R statistical environment (v2.9.2) using the affy package (v1.22.1) of the BioConductor open-source library . Array data were pre-processed with the RMA algorithm . Raw and pre-processed array data are available in the Gene Expression Omnibus repository at NCBI (accession GSE10083). An alternative CDF package was used to ensure each ProbeSet was mapped to a single unique Entrez Gene ID (rae230arnentrezgcdf v12.0.0) .
Statistical analysis of array data
The experimental design employed independent pair-wise analyses between treated and control animals for each strain/line (Figure 1). Following quality control and pre-processing of the microarray data, we performed a general linear modeling analysis. For each gene and each strain we determined the magnitude of differential signal intensity between TCDD-exposed and vehicle-treated animals. Gene-lists were derived separately for each strain using the limma package (v2.18.3) in the R statistical environment (v2.9.2) with a condition-specific design matrix and within-strain pair-wise contrasts. An empirical Bayes moderation of the standard error  and false-discovery rate control of multiple-testing were applied . A significance threshold of padjusted < 0.01 was applied to each contrast. We then scored each gene using a scheme described previously . Briefly, each gene was classified as unaltered (Score: 0), statistically significantly repressed by TCDD (Score: -1), or statistically significantly induced by TCDD (Score: +1) in each strain, and these strain-wise scores were summed.
Unsupervised machine-learning was performed using divisive hierarchical clustering with complete linkage in the R statistical environment (v2.9.2) using the cluster package (v1.12.1). Pearson's correlation was used as a similarity metric and within-row scaling was performed. Venn diagrams were produced using custom R code. R visualizations employed the lattice (v0.17-26) and latticeExtra (v0.6-3) packages.
Functional characterization of responsive genes
To determine if genes perturbed by TCDD are enriched for specific Gene Ontology (GO) terms, we first identified groups of genes dysregulated by TCDD within each strain/line (padjusted < 0.01). Then GO terms enriched in these groups were identified using the GoMiner tool . False-discovery rates (FDRs) were estimated using 1000 permutations. Rat-specific annotations were used. GO terms with false-discovery rates below 5% were included in subsequent analyses.
mRNA quantitation by Real-Time RT-PCR
Total RNA (2 μg) was reverse-transcribed into cDNA using oligo-dT primer p(dT)15 (Roche Applied Science, Laval, QC, Canada) and Superscript II RNA polymerase according to the manufacturer's instructions (Invitrogen, Carlsbad, CA). Real-time PCR was performed using in-house designed primers (with 5' fluorogenic probes, as described previously ) or Applied Biosystems gene expression assays, as described by the manufacturer (Applied Biosystems, Forest City, CA). Primer/probe sequences are in Additional File 2.
Normalized expression was calculated as 2−ΔΔC t, where Ct is the threshold cycle for detecting fluorescence. PCR amplification efficiency was determined from a 10-fold serial dilution of a cDNA pool; efficiency ranged from 90-110% for all genes. Data were normalized to either Actb or Gapdh, genes we previously showed to be suitable as normalization standards for dioxin studies . In the rat model, significant differences in mRNA levels were determined using t-tests (two-tailed, unequal variance). Significant differences in mRNA levels in the Ahr-null mouse model were identified using analysis of variance (ANOVA) with Bonferroni post hoc tests (GraphPad version 4.0).
Okey AB: An Aryl Hydrocarbon Receptor Odyssey to the Shores of Toxicology: The Deichmann Lecture, International Congress of Toxicology-XI. Toxicol Sci. 2007, 98 (1): 5-38. 10.1093/toxsci/kfm096.
Mimura J, Yamashita K, Nakamura K, Morita M, Takagi TN, Nakao K, Ema M, Sogawa K, Yasuda M, Katsuki M: Loss of teratogenic response to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) in mice lacking the Ah (dioxin) receptor. Genes Cells. 1997, 2 (10): 645-654. 10.1046/j.1365-2443.1997.1490345.x.
Peters JM, Narotsky MG, Elizondo G, Fernandez-Salguero PM, Gonzalez FJ, Abbott BD: Amelioration of TCDD-induced teratogenesis in aryl hydrocarbon receptor (AhR)-null mice. Toxicol Sci. 1999, 47 (1): 86-92. 10.1093/toxsci/47.1.86.
Bunger MK, Moran SM, Glover E, Thomae TL, Lahvis GP, Lin BC, Bradfield CA: Resistance to 2,3,7,8-Tetrachlorodibenzo-p-dioxin Toxicity and Abnormal Liver Development in Mice Carrying a Mutation in the Nuclear Localization Sequence of the Aryl Hydrocarbon Receptor. J Biol Chem. 2003, 278 (20): 17767-17774. 10.1074/jbc.M209594200.
Walisser JA, Bunger MK, Glover E, Harstad EB, Bradfield CA: Patent ductus venosus and dioxin resistance in mice harboring a hypomorphic Arnt allele. J Biol Chem. 2004, 279 (16): 16326-16331. 10.1074/jbc.M400784200.
Bunger MK, Glover E, Moran SM, Walisser JA, Lahvis GP, Hsu EL, Bradfield CA: Abnormal Liver Development and Resistance to 2,3,7,8-Tetrachlorodibenzo-p-dioxin Toxicity in Mice Carrying a Mutation in the DNA Binding Domain of the Aryl Hydrocarbon Receptor. Toxicol Sci. 2008, 106: 83-92. 10.1093/toxsci/kfn149.
Hankinson O: The aryl hydrocarbon receptor complex. Annu Rev Pharmacol Toxicol. 1995, 35: 307-340. 10.1146/annurev.pa.35.040195.001515.
Whitlock JP: Mechanistic aspects of dioxin action. Chem Res Toxicol. 1993, 6 (6): 754-763. 10.1021/tx00036a003.
Ma Q: Induction of CYP1A1. The AhR/DRE paradigm: transcription, receptor regulation, and expanding biological roles. Curr Drug Metab. 2001, 2: 149-164. 10.2174/1389200013338603.
Denison MS, Fisher JM, Whitlock JP: The DNA recognition site for the dioxin-Ah receptor complex. Nucleotide sequence and functional analysis. J Biol Chem. 1988, 263 (33): 17221-17224.
Sogawa K, Numayama-Tsuruta K, Takahashi T, Matsushita N, Miura C, Nikawa J, Gotoh O, Kikuchi Y, Fujii-Kuriyama Y: A novel induction mechanism of the rat CYP1A2 gene mediated by Ah receptor-Arnt heterodimer. Biochem Biophys Res Commun. 2004, 318 (3): 746-755. 10.1016/j.bbrc.2004.04.090.
Boutros PC, Moffat ID, Franc MA, Tijet N, Tuomisto J, Pohjanvirta R, Okey AB: Dioxin-responsive AHRE-II gene battery: identification by phylogenetic footprinting. Biochem Biophys Res Commun. 2004, 321 (3): 707-715. 10.1016/j.bbrc.2004.06.177.
Okey AB, Boutros PC, Harper PA: Polymorphisms of human nuclear receptors that control expression of drug-metabolizing enzymes. Pharmacogenet Genomics. 2005, 15 (6): 371-379. 10.1097/01213011-200506000-00001.
Vezina CM, Walker NJ, Olson JR: Subchronic exposure to TCDD, PeCDF, PCB126, and PCB153: effect on hepatic gene expression. Environ Health Perspect. 2004, 112 (16): 1636-1644.
Boverhof DR, Burgoon LD, Tashiro C, Chittim B, Harkema JR, Jump DB, Zacharewski TR: Temporal and dose-dependent hepatic gene expression patterns in mice provide new insights into TCDD-mediated hepatotoxicity. Toxicol Sci. 2005, 85: 1048-1063. 10.1093/toxsci/kfi162.
Boverhof DR, Burgoon LD, Tashiro C, Sharratt B, Chittim B, Harkema JR, Mendrick DL, Zacharewski TR: Comparative toxicogenomic analysis of the hepatotoxic effects of TCDD in Sprague Dawley rats and C57BL/6 mice. Toxicol Sci. 2006, 94 (2): 398-416. 10.1093/toxsci/kfl100.
Fletcher N, Wahlstrom D, Lundberg R, Nilsson CB, Nilsson KC, Stockling K, Hellmold H, Hakansson H: 2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) alters the mRNA expression of critical genes associated with cholesterol metabolism, bile acid biosynthesis, and bile transport in rat liver: A microarray study. Toxicol Appl Pharmacol. 2005, 207 (1): 1-24. 10.1016/j.taap.2004.12.003.
Hayes KR, Bradfield CA: Advances in toxicogenomics. Chem Res Toxicol. 2005, 18 (3): 403-414. 10.1021/tx0496690.
Ovando BJ, Vezina CM, McGarrigle BP, Olson JR: Hepatic gene downregulation following acute and subchronic exposure to 2,3,7,8-tetrachlorodibenzo-p-dioxin. Toxicol Sci. 2006, 94 (2): 428-438. 10.1093/toxsci/kfl111.
Tijet N, Boutros PC, Moffat ID, Okey AB, Tuomisto J, Pohjanvirta R: The aryl hydrocarbon receptor regulates distinct dioxin-dependent and dioxin-independent gene batteries. Mol Pharmacol. 2006, 69: 140-153.
Slatter JG, Cheng O, Cornwell PD, de Souza A, Rockett J, Rushmore T, Hartley D, Evers R, He Y, Dai X: Microarray-based compendium of hepatic gene expression profiles for prototypical ADME gene-inducing compounds in rats and mice in vivo. Xenobiotica. 2006, 36 (10-11): 902-937. 10.1080/00498250600861694.
N'Jai A, Boverhof D, Dere E, Burgoon L, Tan Y, Rowlands J, Budinsky R, Stebbins K, Zacharewski T: Comparative Temporal Toxicogenomic Analysis of TCDD- and TCDF-Mediated Hepatic Effects in Immature Female C57BL/6 Mice. Toxicol Sci. 2008, 103 (2): 285-97. 10.1093/toxsci/kfn053.
Boutros PC, Yan R, Pohjanvirta R, Okey AB: Transcriptomic responses to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) in liver: comparison of rat and mouse. BMC Genomics. 2008, 9: 10.1186/1471-2164-9-419.
Pohjanvirta R, Tuomisto J: Short-term toxicity of 2,3,7,8-tetrachlorodibenzo-p-dioxin in laboratory animals: effects, mechanisms, and animal models. Pharmacol Rev. 1994, 46 (4): 483-549.
Pohjanvirta R, Wong JMY, Li W, Harper PA, Tuomisto J, Okey AB: Point mutation in intron sequence causes altered carboxyl-terminal structure in the aryl hydrocarbon receptor of the most 2,3,7,8-tetrachlorodibenzo-p-dioxin-resistant rat strain. Mol Pharmacol. 1998, 54 (1): 86-93.
Tuomisto JT, Viluksela M, Pohjanvirta R, Tuomisto J: The AH receptor and a novel gene determine acute toxic responses to TCDD: segregation of the resistant alleles to different rat lines. Toxicol Appl Pharmacol. 1999, 155 (1): 71-81. 10.1006/taap.1998.8564.
Moffat ID, Roblin S, Harper PA, Okey AB, Pohjanvirta R: Aryl hydrocarbon receptor splice variants in the dioxin-resistant rat: tissue expression and transactivational activity. Mol Pharmacol. 2007, 72 (4): 956-966. 10.1124/mol.107.037218.
Pohjanvirta R: Transgenic mouse lines expressing rat AH receptor variants--a new animal model for research on AH receptor function and dioxin toxicity mechanisms. Toxicol Appl Pharmacol. 2009, 236 (2): 166-182. 10.1016/j.taap.2009.01.005.
Niittynen M, Simanainen U, Syrjala P, Pohjanvirta R, Viluksela M, Tuomisto J, Tuomisto JT: Differences in acute toxicity syndromes of 2,3,7,8-tetrachlorodibenzo-p-dioxin and 1,2,3,4,7,8-hexachlorodibenzo-p-dioxin in rats. Toxicology. 2007, 235 (1-2): 39-51. 10.1016/j.tox.2007.03.012.
Franc MA, Moffat ID, Boutros PC, Tuomisto JT, Tuomisto J, Pohjanvirta R, Okey AB: Patterns of dioxin-altered mRNA expression in livers of dioxin-sensitive versus dioxin-resistant rats. Arch Toxicol. 2008, 82 (11): 809-830. 10.1007/s00204-008-0303-0.
Simanainen U, Tuomisto JT, Tuomisto J, Viluksela M: Structure-activity relationships and dose responses of polychlorinated dibenzo-p-dioxins for short-term effects in 2,3,7,8-tetrachlorodibenzo-p-dioxin-resistant and -sensitive rat strains. Toxicol Appl Pharmacol. 2002, 181 (1): 38-47. 10.1006/taap.2002.9386.
Denison MS, Whitlock JP: Xenobiotic-inducible transcription of cytochrome P450 genes. J Biol Chem. 1995, 270: 18175-18178. 10.1074/jbc.270.43.25733.
Nebert DW, Dalton TP, Okey AB, Gonzalez FJ: Role of aryl hydrocarbon receptor-mediated induction of the CYP1 enzymes in environmental toxicity and cancer. J Biol Chem. 2004, 279 (23): 23847-23850. 10.1074/jbc.R400004200.
Fernandez-Salguero P, Pineau T, Hilbert DM, McPhail T, Lee SS, Kimura S, Nebert DW, Rudikoff S, Ward JM, Gonzalez FJ: Immune system impairment and hepatic fibrosis in mice lacking the dioxin-binding Ah receptor [see comments]. Science. 1995, 268 (5211): 722-726. 10.1126/science.7732381.
Schmidt JV, Su GH, Reddy JK, Simon MC, Bradfield CA: Characterization of a murine Ahr null allele: involvement of the Ah receptor in hepatic growth and development. Proc Natl Acad Sci USA. 1996, 93 (13): 6731-6736. 10.1073/pnas.93.13.6731.
Pohjanvirta R, Viluksela M, Tuomisto JT, Unkila M, Karasinska J, Franc M-A, Holowenko M, Giannone JV, Harper PA, Tuomisto J: Physicochemical differences in the AH receptors of the most TCDD-susceptible and the most TCDD-resistant rat strains. Toxicol Appl Pharmacol. 1999, 155 (1): 82-95. 10.1006/taap.1998.8565.
Pohjanvirta R, Unkila M, Tuomisto J: TCDD-induced hypophagia is not explained by nausea. Pharmacol Biochem Behav. 1994, 47 (2): 273-282. 10.1016/0091-3057(94)90010-8.
Seefeld MD, Corbett SW, Keesey RE, Peterson RE: Characterization of the wasting syndrome in rats treated with 2,3,7,8-tetrachlorodibenzo-p-dioxin. Toxicol Appl Pharmacol. 1984, 73 (2): 311-322. 10.1016/0041-008X(84)90337-5.
Seefeld MD, Keesey RE, Peterson RE: Body weight regulation in rats treated with 2,3,7,8-tetrachlorodibenzo-p-dioxin. Toxicol Appl Pharmacol. 1984, 76 (3): 526-536. 10.1016/0041-008X(84)90357-0.
Gasiewicz TA, Holscher MA, Neal RA: The effect of total parenteral nutrition on the toxicity of 2,3,7,8-tetrachlorodibenzo-p-dioxin in the rat. Toxicol Appl Pharmacol. 1980, 54 (3): 469-488. 10.1016/0041-008X(80)90174-X.
Weber LW, Lebofsky M, Stahl BU, Gorski JR, Muzi G, Rozman K: Reduced activities of key enzymes of gluconeogenesis as possible cause of acute toxicity of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) in rats. Toxicology. 1991, 66 (2): 133-144. 10.1016/0300-483X(91)90214-L.
Pohjanvirta R, Sankari S, Kulju T, Naukkarinen A, Ylinen M, Tuomisto J: Studies on the role of lipid peroxidation in the acute toxicity of TCDD in rats. Pharmacol Toxicol. 1990, 66 (5): 399-408. 10.1111/j.1600-0773.1990.tb00769.x.
Sato S, Shirakawa H, Tomita S, Ohsaki Y, Haketa K, Tooi O, Santo N, Tohkin M, Furukawa Y, FJ G: Low-dose dioxins alter gene expression related to cholesterol biosynthesis, lipogenesis, and glucose metabolism through the aryl hydrocarbon receptor-mediated pathway in mouse liver. Toxicology and Applied Pharmacology. 2008, 229 (1): 10-9. 10.1016/j.taap.2007.12.029.
Walker BR: Extra-adrenal regeneration of glucocorticoids by 11beta-hydroxysteroid dehydrogenase type 1: physiological regulator and pharmacological target for energy partitioning. Proc Nutr Soc. 2007, 66 (1): 1-8. 10.1017/S002966510700523X.
Morton NM, Seckl JR: 11beta-hydroxysteroid dehydrogenase type 1 and obesity. Front Horm Res. 2008, 36: 146-164. full_text.
Berthiaume M, Laplante M, Festuccia W, Gelinas Y, Poulin S, Lalonde J, Joanisse DR, Thieringer R, Deshaies Y: Depot-specific modulation of rat intraabdominal adipose tissue lipid metabolism by pharmacological inhibition of 11beta-hydroxysteroid dehydrogenase type 1. Endocrinology. 2007, 148 (5): 2391-2397. 10.1210/en.2006-1199.
Doege H, Baillie RA, Ortegon AM, Tsang B, Wu Q, Punreddy S, Hirsch D, Watson N, Gimeno RE, Stahl A: Targeted deletion of FATP5 reveals multiple functions in liver metabolism: alterations in hepatic lipid homeostasis. Gastroenterology. 2006, 130 (4): 1245-1258. 10.1053/j.gastro.2006.02.006.
Hubbard B, Doege H, Punreddy S, Wu H, Huang X, Kaushik VK, Mozell RL, Byrnes JJ, Stricker-Krongrad A, Chou CJ: Mice deleted for fatty acid transport protein 5 have defective bile acid conjugation and are protected from obesity. Gastroenterology. 2006, 130 (4): 1259-1269. 10.1053/j.gastro.2006.02.012.
Viluksela M, Unkila M, Pohjanvirta R, Tuomisto JT, Stahl BU, Rozman KK, Tuomisto J: Effects of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) on liver phosphoenolpyruvate carboxykinase (PEPCK) activity, glucose homeostasis and plasma amino acid concentrations in the most TCDD-susceptible and the most TCDD-resistant rat strains. Arch Toxicol. 1999, 73 (6): 323-336. 10.1007/s002040050626.
Tisdale MJ: The ubiquitin-proteasome pathway as a therapeutic target for muscle wasting. J Support Oncol. 2005, 3 (3): 209-217.
Pihlajamaki J, Boes T, Kim EY, Dearie F, Kim BW, Schroeder J, Mun E, Nasser I, Park PJ, Bianco AC: Thyroid hormone-related regulation of gene expression in human fatty liver. J Clin Endocrinol Metab. 2009, 94 (9): 3521-3529. 10.1210/jc.2009-0212.
Boison D, Scheurer L, Zumsteg V, Rulicke T, Litynski P, Fowler B, Brandner S, Mohler H: Neonatal hepatic steatosis by disruption of the adenosine kinase gene. Proc Natl Acad Sci USA. 2002, 99 (10): 6985-6990. 10.1073/pnas.092642899.
Akira S, Takeda K: Toll-like receptor signalling. Nat Rev Immunol. 2004, 4 (7): 499-511. 10.1038/nri1391.
Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004, 5 (10): R80-10.1186/gb-2004-5-10-r80.
Irizarry RA, Bolstad BM, Collin F, Cope LM, Hobbs B, Speed TP: Summaries of Affymetrix GeneChip probe level data. Nucleic Acids Res. 2003, 31 (4): e15-10.1093/nar/gng015.
Dai M, Wang P, Boyd AD, Kostov G, Athey B, Jones EG, Bunney WE, Myers RM, Speed TP, Akil H: Evolving gene/transcript definitions significantly alter the interpretation of GeneChip data. Nucleic Acids Res. 2005, 33 (20): e175-10.1093/nar/gni179.
Smyth GK: Linear Models and Empirical Bayes Methods for Assessing Differential Expression in Microarray Experiments. Statistical Applications in Genetics and Molecular Biology. 2003, 3 (1): 1-26.
Storey JD, Tibshirani R: Statistical significance for genomewide studies. Proc Natl Acad Sci USA. 2003, 100 (16): 9440-9445. 10.1073/pnas.1530509100.
Zeeberg BR, Feng W, Wang G, Wang MD, Fojo AT, Sunshine M, Narasimhan S, Kane DW, Reinhold WC, Lababidi S: GoMiner: a resource for biological interpretation of genomic and proteomic data. Genome Biol. 2003, 4 (4): R28-10.1186/gb-2003-4-4-r28.
Pohjanvirta R, Niittynen M, Linden J, Boutros PC, Moffat ID, Okey AB: Evaluation of various housekeeping genes for their applicability for normalization of mRNA expression in dioxin-treated rats. Chem Biol Interact. 2006, 160 (2): 134-149. 10.1016/j.cbi.2006.01.001.
The authors thank Ms. Ulla Naukkarinen and Mr. Janne Korkalainen for excellent technical assistance. This work was financially supported by grant number MOP57903 from the Canadian Institutes of Health Research (to ABO and PCB) and grant number 123345 from the Academy of Finland to RP. The funding agencies had no role in study design, in collection or interpretation of data, in the writing of the manuscript or the decision to submit the manuscript for publication.
The project was conceived and designed by ABO and RP. Animal treatment was performed by RP. mRNA extraction and RT-PCR analyses were performed by IDM. Microarray data and functional analyses were performed by PCB. Microarray data visualization was performed by PCB and HC. IDM wrote the first draft of the manuscript, which PCB, ABO and RP edited and all authors approved.