Comparative toxicoproteogenomics of mouse and rat liver identifies TCDD-resistance genes.

The aryl hydrocarbon receptor (AHR) mediates many toxic effects of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). However, the AHR alone does not explain the widely different outcomes among organisms. To identify the other factors involved, we evaluated three transgenic mouse lines, each expressing a different rat AHR isoform (rWT, DEL, and INS) providing widely different resistance to TCDD toxicity, as well as C57BL/6 and DBA/2 mice which exhibit a ~ tenfold divergence in TCDD sensitivity (exposures of 5-1000 µg/kg TCDD). We supplement these with whole-genome sequencing, together with transcriptomic and proteomic analyses of the corresponding rat models, Long-Evans (L-E) and Han/Wistar (H/W) rats (having a ~ 1000-fold difference in their TCDD sensitivities; 100 µg/kg TCDD), to identify genes associated with TCDD-response phenotypes. Overall, we identified up to 50% of genes with altered mRNA abundance following TCDD exposure are associated with a single AHR isoform (33.8%, 11.7%, 5.2% and 0.3% of 3076 genes altered unique to rWT, DEL, C57BL/6 and INS respectively following 1000 µg/kg TCDD). Hepatic Pxdc1 was significantly repressed in all three TCDD-sensitive animal models (C57BL/6 and rWT mice, and L-E rat) after TCDD exposure. Three genes, including Cxxc5, Sugp1 and Hgfac, demonstrated different AHRE-1 (full) motif occurrences within their promoter regions between rat strains, as well as different patterns of mRNA abundance. Several hepatic proteins showed parallel up- or downward alterations with their RNAs, with three genes (SNRK, IGTP and IMPA2) showing consistent, strain-dependent changes. These data show the value of integrating genomic, transcriptomic and proteomic evidence across multi-species models in toxicologic studies.

Compendium of TCDD-mediated transcriptomic response datasets in mammalian model systems.

BACKGROUND: 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) is the most potent congener of the dioxin class of environmental contaminants. Exposure to TCDD causes a wide range of toxic outcomes, ranging from chloracne to acute lethality. The severity of toxicity is highly dependent on the aryl hydrocarbon receptor (AHR). Binding of TCDD to the AHR leads to changes in transcription of numerous genes. Studies evaluating the transcriptional changes brought on by TCDD may provide valuable insight into the role of the AHR in human health and disease. We therefore compiled a collection of transcriptomic datasets that can be used to aid the scientific community in better understanding the transcriptional effects of ligand-activated AHR. RESULTS: Specifically, we have created a datasets package - TCDD.Transcriptomics - for the R statistical environment, consisting of 63 unique experiments comprising 377 samples, including various combinations of 3 species (human derived cell lines, mouse and rat), 4 tissue types (liver, kidney, white adipose tissue and hypothalamus) and a wide range of TCDD exposure times and doses. These datasets have been fully standardized using consistent preprocessing and annotation packages (available as of September 14, 2015). To demonstrate the utility of this R package, a subset of "AHR-core" genes were evaluated across the included datasets. Ahrr, Nqo1 and members of the Cyp family were significantly induced following exposure to TCDD across the studies as expected while Aldh3a1 was induced specifically in rat liver. Inmt was altered only in liver tissue and primarily by rat-AHR. CONCLUSIONS: Analysis of the "AHR-core" genes demonstrates a continued need for studies surrounding the impact of AHR-activity on the transcriptome; genes believed to be consistently regulated by ligand-activated AHR show surprisingly little overlap across species and tissues. Until now, a comprehensive assessment of the transcriptome across these studies was challenging due to differences in array platforms, processing methods and annotation versions. We believe that this package, which is freely available for download ( http://labs.oicr.on.ca/boutros-lab/tcdd-transcriptomics ) will prove to be a highly beneficial resource to the scientific community evaluating the effects of TCDD exposure as well as the variety of functions of the AHR.

2,3,7,8 Tetrachlorodibenzo-p-dioxin-induced RNA abundance changes identify Ackr3, Col18a1, Cyb5a and Glud1 as candidate mediators of toxicity.

2,3,7,8 Tetrachlorodibenzo-p-dioxin (TCDD) is an aromatic, long-lived environmental contaminant. While the pathogenesis of TCDD-induced toxicity is poorly understood, it has been shown that the aryl hydrocarbon receptor (AHR) is required. However, the specific transcriptomic changes that lead to toxic outcomes have not yet been identified. We previously identified a panel of 33 genes that respond to TCDD treatment in two TCDD-sensitive rodent species. To identify genes involved in the onset of hepatic toxicity, we explored 25 of these in-depth using liver from two rat strains: the TCDD-resistant Han/Wistar (H/W) and the TCDD-sensitive Long-Evans (L-E). Time course and dose-response analyses of mRNA abundance following TCDD insult indicate that eight genes are similarly regulated in livers of both strains of rat, suggesting that they are not central to the severe L-E-specific TCDD-induced toxicities. The remaining 17 genes exhibited various divergent mRNA abundances between L-E and H/W strains after TCDD treatment. Several genes displayed a biphasic response where the initial response to TCDD treatment was followed by a secondary response, usually of larger magnitude in L-E liver. This secondary response was most often an exaggeration of the original TCDD-induced response. Only cytochrome b5 type A (microsomal) (Cyb5a) had equivalent TCDD sensitivity to the prototypic AHR-responsive cytochrome P450, family 1, subfamily a, polypeptide 1 (Cyp1a1), while six genes were less sensitive. Four genes showed an early inter-strain difference that was sustained throughout most of the time course (atypical chemokine receptor 3 (Ackr3), collagen, type XVIII, alpha 1 (Col18a1), Cyb5a and glutamate dehydrogenase 1 (Glud1)), and of those genes examined in this study, are most likely to represent genes involved in the pathogenesis of TCDD-induced hepatotoxicity in L-E rats.

Transcriptional profiling of rat white adipose tissue response to 2,3,7,8-tetrachlorodibenzo-ρ-dioxin.

Polychlorinated dibenzodioxins are environmental contaminants commonly produced as a by-product of industrial processes. The most potent of these, 2,3,7,8-tetrachlorodibenzo-ρ-dioxin (TCDD), is highly lipophilic, leading to bioaccumulation. White adipose tissue (WAT) is a major site for energy storage, and is one of the organs in which TCDD accumulates. In laboratory animals, exposure to TCDD causes numerous metabolic abnormalities, including a wasting syndrome. We therefore investigated the molecular effects of TCDD exposure on WAT by profiling the transcriptomic response of WAT to 100µg/kg of TCDD at 1 or 4days in TCDD-sensitive Long-Evans (Turku/AB; L-E) rats. A comparative analysis was conducted simultaneously in identically treated TCDD-resistant Han/Wistar (Kuopio; H/W) rats one day after exposure to the same dose. We sought to identify transcriptomic changes coinciding with the onset of toxicity, while gaining additional insight into later responses. More transcriptional responses to TCDD were observed at 4days than at 1day post-exposure, suggesting WAT shows mostly secondary responses. Two classic AHR-regulated genes, Cyp1a1 and Nqo1, were significantly induced by TCDD in both strains, while several genes involved in the immune response, including Ms4a7 and F13a1 were altered in L-E rats alone. We compared genes affected by TCDD in rat WAT and human adipose cells, and observed little overlap. Interestingly, very few genes involved in lipid metabolism exhibited altered expression levels despite the pronounced lipid mobilization from peripheral fat pads by TCDD in L-E rats. Of these genes, the lipolysis-associated Lpin1 was induced slightly over 2-fold in L-E rat WAT on day 4.

Transcriptional profiling of rat hypothalamus response to 2,3,7,8-tetrachlorodibenzo-ρ-dioxin.

In some mammals, halogenated aromatic hydrocarbon (HAH) exposure causes wasting syndrome, defined as significant weight loss associated with lethal outcomes. The most potent HAH in causing wasting is 2,3,7,8-tetrachlorodibenzo-ρ-dioxin (TCDD), which exerts its toxic effects through the aryl hydrocarbon receptor (AHR). Since TCDD toxicity is thought to predominantly arise from dysregulation of AHR-transcribed genes, it was hypothesized that wasting syndrome is a result of to TCDD-induced dysregulation of genes involved in regulation of food-intake. As the hypothalamus is the central nervous systems' regulatory center for food-intake and energy balance. Therefore, mRNA abundances in hypothalamic tissue from two rat strains with widely differing sensitivities to TCDD-induced wasting syndrome: TCDD-sensitive Long-Evans rats and TCDD-resistant Han/Wistar rats, 23h after exposure to TCDD (100µg/kg) or corn oil vehicle. TCDD exposure caused minimal transcriptional dysregulation in the hypothalamus, with only 6 genes significantly altered in Long-Evans rats and 15 genes in Han/Wistar rats. Two of the most dysregulated genes were Cyp1a1 and Nqo1, which are induced by TCDD across a wide range of tissues and are considered sensitive markers of TCDD exposure. The minimal response of the hypothalamic transcriptome to a lethal dose of TCDD at an early time-point suggests that the hypothalamus is not the predominant site of initial events leading to hypophagia and associated wasting. TCDD may affect feeding behaviour via events upstream or downstream of the hypothalamus, and further work is required to evaluate this at the level of individual hypothalamic nuclei and subregions.

Cross-species transcriptomic analysis elucidates constitutive aryl hydrocarbon receptor activity.

BACKGROUND: Research on the aryl hydrocarbon receptor (AHR) has largely focused on variations in toxic outcomes resulting from its activation by halogenated aromatic hydrocarbons. But the AHR also plays key roles in regulating pathways critical for development, and after decades of research the mechanisms underlying physiological regulation by the AHR remain poorly characterized. Previous studies identified several core genes that respond to xenobiotic AHR ligands across a broad range of species and tissues. However, only limited inferences have been made regarding its role in regulating constitutive gene activity, i.e. in the absence of exogenous ligands. To address this, we profiled transcriptomic variations between AHR-active and AHR-less-active animals in the absence of an exogenous agonist across five tissues, three of which came from rats (hypothalamus, white adipose and liver) and two of which came from mice (kidney and liver). Because AHR status alone has been shown sufficient to alter transcriptomic responses, we reason that by contrasting profiles amongst AHR-variant animals, we may elucidate effects of the AHR on constitutive mRNA abundances. RESULTS: We found significantly more overlap in constitutive mRNA abundances amongst tissues within the same species than from tissues between species and identified 13 genes (Agt, Car3, Creg1, Ctsc, E2f6, Enpp1, Gatm, Gstm4, Kcnj8, Me1, Pdk1, Slc35a3, and Sqrdl) that are affected by AHR-status in four of five tissues. One gene, Creg1, was significantly up-regulated in all AHR-less-active animals. We also find greater overlap between tissues at the pathway level than at the gene level, suggesting coherency to the AHR signalling response within these processes. Analysis of regulatory motifs suggests that the AHR mostly mediates transcriptional regulation via direct binding to response elements. CONCLUSIONS: These findings, though preliminary, present a platform for further evaluating the role of the AHR in regulation of constitutive mRNA levels and physiologic function.

TCDD dysregulation of 13 AHR-target genes in rat liver.

Despite several decades of research, the complete mechanism by which 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) and other xenobiotic agonists of the aryl hydrocarbon receptor (AHR) cause toxicity remains unclear. While it has been shown that the AHR is required for all major manifestations of toxicity, the specific downstream changes involved in the development of toxic phenotypes remain unknown. Here we examine a panel of 13 genes that are AHR-regulated in many species and tissues. We profiled their hepatic mRNA abundances in two rat strains with very different sensitivities to TCDD: the TCDD-sensitive Long-Evans (Turku/AB; L-E) and the TCDD-resistant Han/Wistar (Kuopio; H/W). We evaluated doses ranging from 0 to 3000µg/kg at 19h after TCDD exposure and time points ranging from 1.5 to 384h after exposure to 100µg/kg TCDD. Twelve of 13 genes responded to TCDD in at least one strain, and seven of these showed statistically significant inter-strain differences in the time course analysis (Aldh3a1, Cyp1a2, Cyp1b1, Cyp2a1, Fmo1, Nfe2l2 and Nqo1). Cyp2s1 did not respond to TCDD in either rat strain. Five genes exhibited biphasic responses to TCDD insult (Ahrr, Aldh3a1, Cyp1b1, Nfe2l2 and Nqo1), suggesting a secondary event, such as association with additional transcriptional modulators. Of the 12 genes that responded to TCDD during the dose-response analysis, none had an ED50 equivalent to that of Cyp1a1, the most sensitive gene in this study, while nine genes responded to doses at least 10-100 fold higher, in at least one strain (Ahrr (L-E), Aldh3a1 (both), Cyp1a2 (both), Cyp1b1 (both), Cyp2a1 (L-E), Inmt (both), Nfe2l2 (L-E), Nqo1 (L-E) and Tiparp (both)). These data shed new light on the association of the AHR target genes with TCDD toxicity, and in particular the seven genes exhibiting strain-specific differences represent strong candidate mediators of Type-II toxicities.

Systematic evaluation of medium-throughput mRNA abundance platforms.

Profiling of mRNA abundances with high-throughput platforms such as microarrays and RNA-seq has become an important tool in both basic and biomedical research. However, these platforms remain prone to systematic errors and have challenges in clinical and industrial applications. As a result, it is standard practice to validate a subset of key results using alternate technologies. Similarly, clinical and industrial applications typically involve transitions from a high-throughput discovery platform to medium-throughput validation ones. These medium-throughput validation platforms have high technical reproducibility and reduced sample input needs, and low sensitivity to sample quality (e.g., for processing FFPE specimens). Unfortunately, while medium-throughput platforms have proliferated, there are no comprehensive comparisons of them. Here we fill that gap by comparing two key medium-throughput platforms--NanoString's nCounter Analysis System and ABI's OpenArray System--to gold-standard quantitative real-time RT-PCR. We quantified 38 genes and positive and negative controls in 165 samples. Signal:noise ratios, correlations, dynamic range, and detection accuracy were compared across platforms. All three measurement technologies showed good concordance, but with divergent price/time/sensitivity trade-offs. This study provides the first detailed comparison of medium-throughput RNA quantification platforms and provides a template and a standard data set for the evaluation of additional technologies.

Inter-strain heterogeneity in rat hepatic transcriptomic responses to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD).

The biochemical and toxic effects of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) have been the subject of intense study for decades. It is now clear that essentially all TCDD-induced toxicities are mediated by DNA-protein interactions involving the Aryl Hydrocarbon Receptor (AHR). Nevertheless, it remains unknown which AHR target genes cause TCDD toxicities. Several groups, including our own, have developed rodent model systems to probe these questions. mRNA expression profiling of these model systems has revealed significant inter-species heterogeneity in rodent hepatic responses to TCDD. It has remained unclear if this variability also exists within a species, amongst rodent strains. To resolve this question, we profiled the hepatic transcriptomic response to TCDD of diverse rat strains (L-E, H/W, F344 and Wistar rats) and two lines derived from L-E×H/W crosses, at consistent age, sex, and dosing (100 µg/kg TCDD for 19 h). Using this uniquely consistent dataset, we show that the majority of TCDD-induced alterations in mRNA abundance are strain/line-specific: only 11 genes were affected by TCDD across all strains, including well-known dioxin-responsive genes such as Cyp1a1 and Nqo1. Our analysis identified two novel universally dioxin-responsive genes as well as 4 genes induced by TCDD in dioxin-sensitive rats only. These 6 genes are strong candidates to explain TCDD-related toxicities, so we validated them using 152 animals in time-course (0 to 384 h) and dose-response (0 to 3000 µg/kg) experiments. This study reveals that different rat strains exhibit dramatic transcriptional heterogeneity in their hepatic responses to TCDD and that inter-strain comparisons can help identify candidate toxicity-related genes.

mRNA levels in control rat liver display strain-specific, hereditary, and AHR-dependent components.

Rat is a major model organism in toxicogenomics and pharmacogenomics. Hepatic mRNA profiles after treatment with xenobiotic chemicals are used to predict and understand drug toxicity and mechanisms. Surprisingly, neither inter- and intra-strain variability of mRNA abundances in control rats nor the heritability of rat mRNA abundances yet been established. We address these issues by studying five populations: the popular Sprague-Dawley strain, sub-strains of Long-Evans and Wistar rats, and two lines derived from crosses between the Long-Evans and Wistar sub-strains. Using three independent techniques--variance analysis, linear modelling, and unsupervised pattern recognition--we characterize extensive intra- and inter-strain variability in mRNA levels. We find that both sources of variability are non-random and are enriched for specific functional groups. Specific transcription-factor binding-sites are enriched in their promoter regions and these genes occur in "islands" scattered throughout the rat genome. Using the two lines generated by crossbreeding we tested heritability of hepatic mRNA levels: the majority of rat genes appear to exhibit directional genetics, with only a few interacting loci. Finally, a comparison of inter-strain heterogeneity between mouse and rat orthologs shows more heterogeneity in rats than mice; thus rat and mouse heterogeneity are uncorrelated. Our results establish that control hepatic mRNA levels are relatively homogeneous within rat strains but highly variable between strains. This variability may be related to increased activity of specific transcription-factors and has clear functional consequences. Future studies may take advantage of this phenomenon by surveying panels of rat strains.

Hepatic transcriptomic responses to TCDD in dioxin-sensitive and dioxin-resistant rats during the onset of toxicity.

The dioxin congener 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) causes a wide range of toxic effects in rodent species, all of which are mediated by a ligand-dependent transcription-factor, the aryl hydrocarbon receptor (AHR). The Han/Wistar (Kuopio) (H/W) strain shows exceptional resistance to many TCDD-induced toxicities; the LD50 of > 9600 µg/kg for H/W rats is higher than for any other wild-type mammal known. We previously showed that this resistance primarily results from H/W rats expressing a variant AHR isoform that has a substantial portion of the AHR transactivation domain deleted. Despite this large deletion, H/W rats are not entirely refractory to the effects of TCDD; the variant AHR in these animals remains fully competent to up-regulate well-known dioxin-inducible genes. TCDD-sensitive (Long-Evans, L-E) and resistant (H/W) rats were treated with either corn-oil (with or without feed-restriction) or 100 µg/kg TCDD for either four or ten days. Hepatic transcriptional profiling was done using microarrays, and was validated by RT-PCR analysis of 41 genes. A core set of genes was altered in both strains at all time points tested, including CYP1A1, CYP1A2, CYP1B1, Nqo1, Aldh3a1, Tiparp, Exoc3, and Inmt. Outside this core, the strains differed significantly in the breadth of response: three-fold more genes were altered in L-E than H/W rats. At ten days almost all expressed genes were dysregulated in L-E rats, likely reflecting emerging toxic responses. Far fewer genes were affected by feed-restriction, suggesting that only a minority of the TCDD-induced changes are secondary to the wasting syndrome.

Flavin-containing monooxygenase-3: induction by 3-methylcholanthrene and complex regulation by xenobiotic chemicals in hepatoma cells and mouse liver.

Flavin-containing monooxygenases often are thought not to be inducible but we recently demonstrated aryl hydrocarbon receptor (AHR)-dependent induction of FMO mRNAs in mouse liver by 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) (Celius et al., Drug Metab Dispos 36:2499, 2008). We now evaluated FMO induction by other AHR ligands and xenobiotic chemicals in vivo and in mouse Hepa1c1c7 hepatoma cells (Hepa-1). In mouse liver, 3-methylcholanthrene (3MC) induced FMO3 mRNA 8-fold. In Hepa-1 cells, 3MC and benzo[a]pyrene (BaP) induced FMO3 mRNA >30-fold. Induction by 3MC and BaP was AHR dependent but, surprisingly, the potent AHR agonist, TCDD, did not induce FMO3 mRNA in Hepa-1 cells nor did chromatin immunoprecipitation assays detect recruitment of AHR or ARNT to Fmo3 regulatory elements after exposure to 3MC in liver or in Hepa-1 cells. However, in Hepa-1, 3MC and BaP (but not TCDD) caused recruitment of p53 protein to a p53 response element in the 5'-flanking region of the Fmo3 gene. We tested the possibility that FMO3 induction in Hepa-1 cells might be mediated by Nrf2/anti-oxidant response pathways, but agents known to activate Nrf2 or to induce oxidative stress did not affect FMO3 mRNA levels. The protein synthesis inhibitor, cycloheximide (which causes "superinduction" of CYP1A1 mRNA in TCDD-treated cells), by itself caused dramatic upregulation (>300-fold) of FMO3 mRNA in Hepa-1 suggesting that cycloheximide prevents synthesis of a labile protein that suppresses FMO3 expression. Although FMO3 mRNA is highly induced by 3MC or TCDD in mouse liver and in Hepa-1 cells, FMO protein levels and FMO catalytic function showed only modest elevation.

Aryl hydrocarbon receptor (AHR)-regulated transcriptomic changes in rats sensitive or resistant to major dioxin toxicities.

BACKGROUND: 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. RESULTS: 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. CONCLUSION: 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.

Red meat intake, doneness, polymorphisms in genes that encode carcinogen-metabolizing enzymes, and colorectal cancer risk.

Colorectal cancer literature regarding the interaction between polymorphisms in carcinogen-metabolizing enzymes and red meat intake/doneness is inconsistent. A case-control study was conducted to evaluate the interaction between red meat consumption, doneness, and polymorphisms in carcinogen-metabolizing enzymes. Colorectal cancer cases diagnosed 1997 to 2000, ages 20 to 74 years, were identified through the population-based Ontario Cancer Registry and recruited by the Ontario Family Colorectal Cancer Registry. Controls were sex-matched and age group-matched random sample of Ontario population. Epidemiologic and food questionnaires were completed by 1,095 cases and 1,890 controls; blood was provided by 842 and 1,251, respectively. Multivariate logistic regression was used to obtain adjusted odds ratio (OR) estimates. Increased red meat intake was associated with increased colorectal cancer risk [OR (> 5 versus < or = 2 servings/wk), 1.67 (1.36-2.05)]. Colorectal cancer risk also increased significantly with well-done meat intake [OR (> 2 servings/wk well-done versus < or = 2 servings/wk rare-regular), 1.57 (1.27-1.93)]. We evaluated interactions between genetic variants in 15 enzymes involved in the metabolism of carcinogens in overcooked meat (cytochrome P450, glutathione S-transferase, UDP-glucuronosyltransferases, SULT, NAT, mEH, and AHR). CYP2C9 and NAT2 variants were associated with colorectal cancer risk. Red meat intake was associated with increased colorectal cancer risk regardless of genotypes; however, CYP1B1 combined variant and SULT1A1-638G>A variant significantly modified the association between red meat doneness intake and colorectal cancer risk. In conclusion, well-done red meat intake was associated with an increased risk of colorectal cancer regardless of carcinogen-metabolizing genotype, although our data suggest that persons with CYP1B1 and SULT1A1 variants had the highest colorectal cancer risk.

Transcriptomic responses to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) in liver: comparison of rat and mouse.

BACKGROUND: Mouse and rat models are mainstays in pharmacology, toxicology and drug development -- but differences between strains and between species complicate data interpretation and application to human health. Dioxin-like polyhalogenated aromatic hydrocarbons represent a major class of environmentally and economically relevant toxicants. In mammals dioxin exposure leads to a broad spectrum of adverse affects, including hepatotoxicity of varying severity. Several studies have shown that dioxins extensively alter hepatic mRNA levels. Surprisingly, though, analysis of a limited portion of the transcriptome revealed that rat and mouse responses diverge greatly (Boverhof et al. Toxicol Sci 94:398-416, 2006). RESULTS: We employed oligonucleotide arrays to compare the response of 8,125 rat and mouse orthologs. We confirmed that there is limited inter-species overlap in dioxin-responsive genes. Rat-specific and mouse-specific genes are enriched for specific functional groups which differ between species, conceivably accounting for species-specificities in liver histopathology. While no evidence for the involvement of copy-number variation was found, extensive inter-species variation in the transcriptional-regulatory network was identified; Nr2f1 and Fos emerged as candidates to explain species-specific and species-independent responses, respectively. CONCLUSION: Our results suggest that a small core of genes is responsible for mediating the similar features of dioxin hepatotoxicity in rats and mice but non-overlapping pathways are simultaneously at play to result in distinctive histopathological outcomes. The extreme divergence between mouse and rat transcriptomic responses appears to reflect divergent transcriptional-regulatory networks. Taken together, these data suggest that both rat and mouse models should be used to screen the acute hepatotoxic effects of drugs and toxic compounds.

Aryl hydrocarbon receptor-dependent induction of flavin-containing monooxygenase mRNAs in mouse liver.

Flavin-containing monooxygenases (FMOs) are important in detoxication but generally are considered not to be inducible by xenobiotics. Our recent microarray studies revealed induction of FMO2 and FMO3 mRNAs by 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) in liver of mice with wild-type aryl hydrocarbon receptor (AHR) but not in Ahr-null mice. The aim of the present study was to delineate mechanisms of FMO regulation. In adult male mice, basal FMO3 mRNA is low but was induced 6-fold at 4 h and 6000-fold at 24 h. The ED50 was approximately 1 microg/kg for FMO2 and FMO3, similar to that for the classic AHR-regulated gene, Cyp1a1. In adult female mice basal FMO3 mRNA is high and was not induced at 4 h but was elevated 8-fold at 24 h. FMO5 mRNA was significantly down-regulated by TCDD in both male and female adult mice. Juvenile mice show no sex difference in response to TCDD; FMO3 was induced 4 to 6-fold by TCDD in both sexes. Chromatin immunoprecipitation demonstrated recruitment of AHR and aryl hydrocarbon nuclear translocator proteins to Fmo3 regulatory regions, suggesting that induction by TCDD is a primary AHR-mediated event. Although FMO2 and FMO3 mRNAs were highly induced by TCDD in adult males, overall FMO catalytic activity increased only modestly. In contrast to the striking up-regulation of FMO2 and FMO3 in mouse liver, TCDD has little effect on FMO mRNA in rat liver. However, FMO2 and FMO3 mRNAs were highly induced in transgenic mice that express wild-type rat AHR, indicating that lack of induction in rat is not due to an incompetent AHR in this species.

Patterns of dioxin-altered mRNA expression in livers of dioxin-sensitive versus dioxin-resistant rats.

Dioxins exert their major toxicologic effects by binding to the aryl hydrocarbon receptor (AHR) and altering gene transcription. Numerous dioxin-responsive genes previously were identified both by conventional biochemical and molecular techniques and by recent mRNA expression microarray studies. However, of the large set of dioxin-responsive genes the specific genes whose dysregulation leads to death remain unknown. To identify specific genes that may be involved in dioxin lethality we compared changes in liver mRNA levels following exposure to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) in three strains/lines of dioxin-sensitive rats with changes in three dioxin-resistant rat strains/lines. The three dioxin-resistant strains/lines all harbor a large deletion in the transactivation domain of the aryl hydrocarbon receptor (AHR). Despite this deletion, many genes exhibited a "Type-I" response-that is, their responses were similar in dioxin-sensitive and dioxin-resistant rats. Several genes that previously were well established as being dioxin-responsive or under AHR regulation emerged as Type-I responses (e.g. CYP1A1, CYP1A2, CYP1B1 and Gsta3). In contrast, a relatively small number of genes exhibited a Type-II response-defined as a difference in responsiveness between dioxin-sensitive and dioxin-resistant rat strains. Type-II genes include: malic enzyme 1, ubiquitin C, cathepsin L, S-adenosylhomocysteine hydrolase and ferritin light chain 1. In silico searches revealed that AH response elements are conserved in the 5'-flanking regions of several genes that respond to TCDD in both the Type-I and Type-II categories. The vast majority of changes in mRNA levels in response to 100 microg/kg TCDD were strain-specific; over 75% of the dioxin-responsive clones were affected in only one of the six strains/lines. Selected genes were assessed by quantitative RT-PCR in dose-response and time-course experiments and responses of some genes were assessed in Ahr-null mice to determine if their response was AHR-dependent. Type-II genes may lie in pathways that are central to the difference in susceptibility to TCDD lethality in this animal model.

Genome-wide effects of acute progressive feed restriction in liver and white adipose tissue.

Acute progressive feed restriction (APFR) represents a specific form of caloric restriction in which feed availability is increasingly curtailed over a period of a few days to a few weeks. It is often used for control animals in toxicological and pharmacological studies on compounds causing body weight loss to equalize weight changes between experimental and control groups and thereby, intuitively, to also set their metabolic states to the same phase. However, scientific justification for this procedure is lacking. In the present study, we analyzed by microarrays the impact on hepatic gene expression in rats of two APFR regimens that caused identical diminution of body weight (19%) but differed slightly in duration (4 vs. 10 days). In addition, white adipose tissue (WAT) was also subjected to the transcriptomic analysis on day-4. The data revealed that the two regimens led to distinct patterns of differentially expressed genes in liver, albeit some major pathways of energy metabolism were similarly affected (particularly fatty acid and amino acid catabolism). The reason for the divergence appeared to be entrainment by the longer APFR protocol of peripheral oscillator genes, which resulted in derailment of circadian rhythms and consequent interaction of altered diurnal fluctuations with metabolic adjustments in gene expression activities. WAT proved to be highly unresponsive to the 4-day APFR as only 17 mRNA levels were influenced by the treatment. This study demonstrates that body weight is a poor proxy of metabolic state and that the customary protocols of feed restriction can lead to rhythm entrainment.

microRNAs in adult rodent liver are refractory to dioxin treatment.

Dioxin-like chemicals are well known for their ability to upregulate expression of numerous genes via the AH receptor (AHR). However, recent transcriptomic analyses in several laboratories indicate that dioxin-like chemicals or AHR genotype itself also can downregulate levels of mRNAs encoded by numerous genes. The mechanism responsible for such downregulation is unknown. We hypothesized that microRNAs (miRNAs), which have emerged as powerful negative regulators of mRNA levels in several systems, might be responsible for mRNA downregulation in dioxin/AHR pathways. We used two miRNA array platforms as well as quantitative reverse transcriptase-polymerase chain reaction to measure miRNA levels in wild-type (WT) versus Ahr-null mice, in dioxin-sensitive Long-Evans (L-E; Turku/AB) rats versus dioxin-resistant Han/Wistar (H/W; Kuopio) rats and in rat 5L and mouse Hepa-1 hepatoma cells in culture. Treatment with 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) in vivo caused few changes in miRNA levels in mouse or rat livers, and those changes that were statistically significant were of modest magnitude. Hepatoma cells in culture also exhibited few changes in miRNA levels in response to TCDD. AHR genotype had little effect on hepatic miRNA levels, either in constitutive expression or in response to TCDD-only a few miRNAs differed in expression between Ahr-null mice compared to mice with WT AHR or between L-E rats (that have WT AHR) compared to H/W rats (whose AHR has a large deletion in the transactivation domain). It is unlikely that mRNA downregulation by dioxins is mediated by miRNAs, nor are miRNAs likely to play a significant role in dioxin toxicity in adult rodent liver.

Aryl hydrocarbon receptor splice variants in the dioxin-resistant rat: tissue expression and transactivational activity.

The AHR locus encodes the aryl hydrocarbon receptor (AHR), a transcriptional regulator of multiple drug-metabolizing enzymes and mediator of toxicity of dioxin-like chemicals. The Han/Wistar (Kuopio) rat strain (H/W) is remarkably resistant to lethal effects of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) because of a point mutation in the exon/intron 10 boundary in AHR genomic structure that leads to use of 3 alternative cryptic splice sites, potentially creating 3 alternative transcripts and 2 protein products. The deletion variant (DV), which lacks 43 amino acids in the transactivation domain, has the highest intrinsic transactivation activity in vitro; amino acids 766 to 783 suppress transactivation function. However, DV expression levels in H/W rats in vivo are low in liver, lung, thymus, kidney, and testis; insertion variant mRNAs (IVs) are the dominant mRNA forms in H/W rats in which wild-type AHR mRNA is undetectable. In dioxin-sensitive rat strains and lines that are homozygous for wild-type AHR alleles, wild-type AHR mRNA is the most abundant transcript but some IV transcripts are detectable. TCDD treatment in vivo increases transcript levels for both the DV and IVs in H/W rats and increases wild-type transcript levels in dioxin-sensitive rats but does not alter which transcript forms are expressed. In silico modeling indicates that the DV mRNA has lost considerable secondary structure, whereas at the protein level, the transactivation domain of the IV in the dioxin-resistant H/W rat has greater alpha-helical content and a more hydrophobic terminus than wild-type AHR, which may produce a protein conformation that is less amenable to interaction with other regulatory proteins.