Comparative temporal toxicogenomic analysis of TCDD- and TCDF-mediated hepatic effects in immature female C57BL/6 mice

Metatranscriptomic Analysis of the Mouse Gut Microbiome Response to the Persistent Organic Pollutant 2,3,7,8-Tetrachlorodibenzofuran

Persistent organic pollutants (POPs) are important environmental chemicals and continued study of their mechanism of action remains a high priority. POPs, such as 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), 2,3,7,8-tetrachlorodibenzofuran (TCDF), and polychlorinated biphenyls (PCBs), are widespread environmental contaminants that are agonists for the aryl hydrocarbon receptor (AHR). Activation of the AHR modulates the gut microbiome community structure and function, host immunity, and the host metabolome. In the current study, male C57BL6/J mice were exposed, via the diet, to 5 µg/kg body weight (BW) TCDF or 24 µg/kg BW of TCDF every day for 5 days. The functional and structural changes imparted by TCDF exposure to the gut microbiome and host metabolome were explored via 16S rRNA gene amplicon sequencing, metabolomics, and bacterial metatranscriptomics. Significant changes included increases in lipopolysaccharide (LPS) biosynthesis gene expression after exposure to 24 µg/kg BW of TCDF. Increases in LPS biosynthesis were confirmed with metabolomics and LPS assays using serum obtained from TCDF-treated mice. Significant increases in gene expression within aspartate and glutamate metabolism were noted after exposure to 24 µg/kg BW of TCDF. Together, these results suggest that after exposure to 24 µg/kg BW of TCDF, the gut microbiome increases the production of LPS and glutamate to promote localized gut inflammation, potentially using glutamate as a stress response.

2,3,7,8-TCDD enhances the sensitivity of mice to concanavalin A immune-mediated liver injury

Inflammation plays a major role in immune-mediated liver injury, and exposure to environmental pollutants such as 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) has been reported to alter the inflammatory response as well as affect immune cell activity. In this study, we tested the hypothesis that TCDD pretreatment exacerbates hepatotoxicity in a murine model of immune-mediated liver injury induced by concanavalin A (Con A) administration. Mice were pretreated with 30 μg/kg TCDD or vehicle control on day zero and then given either Con A or saline intravenously on day four. Mice treated with TCDD did not develop liver injury; however, TCDD pretreatment increased liver injury resulting from moderate doses of Con A (4-10 mg/kg). TCDD-pretreated mice had altered plasma concentrations of inflammatory cytokines, including interferon gamma (IFNγ), and TCDD/Con A-induced hepatotoxicity was attenuated in IFNγ knockout mice. At various times after treatment, intrahepatic immune cells were isolated, and expression of cell activation markers as well as cytolytic proteins was determined. TCDD pretreatment increased the proportion of activated natural killer T (NKT) cells and the percent of cells expressing Fas ligand (FasL) after Con A administration. In addition FasL knockout mice and mice treated with CD18 antiserum were both protected from TCDD/Con A-induced hepatotoxicity, suggesting a requirement for direct cell-cell interaction between effector immune cells and parenchymal cell targets in the development of liver injury from TCDD/Con A treatment. In summary, exposure to TCDD increased NKT cell activation and exacerbated immune-mediated liver injury induced by Con A through a mechanism involving IFNγ and FasL expression.

Beneficial effect of astaxanthin on 2,3,7,8-tetrachlorodibenzo-p-dioxin-induced liver injury in rats

Exposure to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) represents a potential health risk and hepatotoxicity. Astaxanthin (ASTA) exhibits antioxidant properties and can influence hepatotoxicity. Therefore, the present study was carried out for using ASTA against hepatotoxicity induced by TCDD in the liver of rats. Animals were treated intraperitoneally daily with TCDD (8 µg/kg body weight (b.w.)), ASTA (12.5 mg/kg b.w. and 25 mg/kg b.w.) and TCDD plus ASTA (12.5 and 25 mg/kg b.w.) for 21 days. TCDD significantly decreased the activities of antioxidant enzymes and resulted in serious pathological findings. Moreover, the rate of micronucleus (MN) in hepatocytes increased after treating with TCDD. The activities of enzymes, frequencies of MNs and liver histology in lower dosage group of ASTA remained unchanged compared with the control group. In rats treated with ASTA, at higher dosage alone, the MNs remained unchanged and the activities of antioxidant enzymes significantly increased. The presence of ASTA (except for lower dose) with TCDD alleviated its pathological effects in hepatic tissue. ASTA also prevented the suppression of antioxidant enzymes in the livers of animals exposed to TCDD and displayed a strong protective effect against MNs. Thus, the present findings might provide new insight into the development of therapeutic and preventive approaches of TCDD toxicity.

Characterization of polybrominated diphenyl ether toxicity in Wistar Han rats and use of liver microarray data for predicting disease susceptibilities

The toxicity of polybrominated diphenyl ethers (PBDEs), flame-retardant components, was characterized in offspring from Wistar Han dams exposed by gavage to a PBDE mixture (DE71) starting at gestation day 6 and continuing to weaning on postnatal day (PND) 21. Offspring from the dams underwent PBDE direct dosing by gavage at the same dose as their dams from PND 12 to PND 21, and then after weaning for another thirteen weeks. Liver samples were collected at PND 22 and week 13 for liver gene expression analysis (Affymetrix Rat Genome 230 2.0 Array). Treatment with PBDE induced 1,066 liver gene transcript changes in females and 1,200 transcriptional changes in males at PND 22 (false discovery rate < 0.01), but only 263 liver transcriptional changes at thirteen weeks in male rats (false discovery rate < 0.05). No significant differences in dose response were found between male and female pups. Transcript changes at PND 22 coded for proteins in xenobiotic, sterol, and lipid metabolism, and cell cycle regulation, and overlapped rodent liver transcript patterns after a high-fat diet or phenobarbital exposure. These findings, along with the observed PBDE-induced liver hypertrophy and vacuolization, suggest that long-term PBDE exposure has the potential to modify cell functions that contribute to metabolic disease and/or cancer susceptibilities.

Modulatory effect of l-glutamine on 2,3,7,8 tetrachlorodibenzo-p-dioxin-induced liver injury in rats

The aim of this study was to explore the effectiveness of l-glutamine (Gln) in alleviating the toxicity of 2,3,7,8-tetrachlorodibenzo- p-dioxin (TCDD) in liver of rats. Rats were intraperitoneally administered Gln and TCDD doses daily for 21 days. In the liver of rats, the biochemical tests, pathological examination and micronucleus (MN) test were performed. TCDD significantly decreased the activities of antioxidant enzymes and serious pathological findings. Moreover, the rate of MNs in hepatocytes increased after treatment with dioxin. In rats treated with Gln alone, the MNs remained unchanged, but the ratio of glutathione (GSH) and the activity of glutathione peroxidase (GSH-Px) were significantly increased. Gln also prevented the suppression of GSH-Px (except for superoxide dismutase and catalase) and GSH in the livers of animals exposed to TCDD and displayed a strong protective effect against MNs. Thus, our findings for Gln might provide new insight into the development of therapeutic and preventive approaches in TCDD toxicity.

Comparative metabolomic and genomic analyses of TCDD-elicited metabolic disruption in mouse and rat liver

2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) elicits a broad spectrum of species-specific effects that have not yet been fully characterized. This study compares the temporal effects of TCDD on hepatic aqueous and lipid metabolite extracts from immature ovariectomized C57BL/6 mice and Sprague-Dawley rats using gas chromatography-mass spectrometry and nuclear magnetic resonance-based metabolomic approaches and integrates published gene expression data to identify species-specific pathways affected by treatment. TCDD elicited metabolite and gene expression changes associated with lipid metabolism and transport, choline metabolism, bile acid metabolism, glycolysis, and glycerophospholipid metabolism. Lipid metabolism is altered in mice resulting in increased hepatic triacylglycerol as well as mono- and polyunsaturated fatty acid (FA) levels. Mouse-specific changes included the induction of CD36 and other cell surface receptors as well as lipases- and FA-binding proteins consistent with hepatic triglyceride and FA accumulation. In contrast, there was minimal hepatic fat accumulation in rats and decreased CD36 expression. However, choline metabolism was altered in rats, as indicated by decreases in betaine and increases in phosphocholine with the concomitant induction of betaine-homocysteine methyltransferase and choline kinase gene expression. Results from these studies show that aryl hydrocarbon receptor-mediated differential gene expression could be linked to metabolite changes and species-specific alterations of biochemical pathways.

Acute disruption of bone marrow hematopoiesis by benzo(a)pyrene is selectively reversed by aryl hydrocarbon receptor-mediated processes

Bone marrow (BM) hematopoietic cells are selectively sensitive to polycyclic aromatic hydrocarbons (PAH) in vivo. 7,12-Dimethylbenz(a)anthracene (DMBA), but not benzo(a)pyrene (BP), depletes BM hematopoietic cells in C57BL/6 mice. This difference is due to a BP-selective aryl hydrocarbon receptor (AhR)-mediated recovery. Colony-forming unit assays show suppression of lymphoid progenitors by each PAH within 6 h but a subsequent recovery, exclusively after BP treatment. Suppression of myeloid progenitors (6 h) occurs only for DMBA. Each progenitor responded equally to DMBA and BP in congenic mice expressing the PAH-resistant AhR (AhR(d)). AhR, therefore, mediates this BP recovery in each progenitor type. These PAH suppressions depend on Cyp1b1-mediated metabolism. Paradoxically, few genes responded to DMBA, whereas 12 times more responded to BP. Progenitor suppression by DMBA, therefore, occurs with minimal effects on the general BM population. Standard AhR-mediated stimulations (Cyp1a1, Cyp1b1, Ahrr) were similar for each PAH and for the specific agonist 2,3,7,8-tetrachlorodibenzo-p-dioxin but were absent in AhR(d) mice. A group of 12 such AhR responses was sustained from 6 to 24 h. A second, larger set of BP responses (chemokines, cytokines, cyclooxygenase 2) differed in two respects; DMBA responses were low and BP responses declined extensively from 6 to 24 h. A third cluster exhibited BP-induced increases in protective genes (Nqo1, GST-mu) that appeared only after 12 h. Conversion of BP to quinones contributes oxidative signaling not seen with DMBA. We propose that genes in this second cluster, which share oxidative signaling and AhR activation, provide the AhR-dependent protection of hematopoietic progenitors seen for BP.

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 LD₅₀ 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.

Toxicogenomic analysis of exposure to TCDD, PCB126 and PCB153: identification of genomic biomarkers of exposure to AhR ligands

Background

Two year cancer bioassays conducted by the National Toxicology Program have shown chronic exposure to dioxin-like compounds (DLCs) to lead to the development of both neoplastic and non-neoplastic lesions in the hepatic tissue of female Sprague Dawley rats. Most, if not all, of the hepatotoxic effects induced by DLC's are believed to involve the binding and activation of the transcription factor, the aryl hydrocarbon receptor (AhR). Toxicogenomics was implemented to identify genomic responses that may be contributing to the development of hepatotoxicity in rats.

Results

Through comparative analysis of time-course microarray data, unique hepatic gene expression signatures were identified for the DLCs, 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) (100 ng/kg/day) and 3,3',4,4',5-pentachlorobiphenyl (PCB126) (1000 ng/kg/day) and the non-DLC 2,2',4,4',5,5',-hexachlorobiphenyl (PCB153) (1000 μg/kg/day). A common time independent signature of 41 AhR genomic biomarkers was identified which exhibited at least a 2-fold change in expression following subchronic (13-wk) and chronic (52-wk) p.o. exposure to TCDD and PCB126, but not the non DLC, PCB153. Real time qPCR analysis validated that 30 of these genes also exhibited at least a 2-fold change in hepatic expression at 24 hr following a single exposure to TCDD (5 μg/kg, po). Phenotypic anchoring was conducted which identified forty-six genes that were differently expressed both following chronic p.o. exposure to DLCs and in previously reported studies of cholangiocarcinoma or hepatocellular adenoma.

Conclusions

Together these analyses provide a comprehensive description of the genomic responses which occur in rat hepatic tissue with exposure to AhR ligands and will help to isolate those genomic responses which are contributing to the hepatotoxicity observed with exposure to DLCs. In addition, the time independent gene expression signature of the AhR ligands may assist in identifying other agents with the potential to elicit dioxin-like hepatotoxic responses.

Automated dose-response analysis and comparative toxicogenomic evaluation of the hepatic effects elicited by TCDD, TCDF, and PCB126 in C57BL/6 mice

The toxic equivalency factor (TEF) approach recommended by the World Health Organization is used to quantify dioxin-like exposure concentrations for mixtures of polychlorinated dibenzo-dioxins, -furans, and polychlorinated biphenyls (PCBs), including 2,3,7,8-tetrachlorodibenzofuran (TCDF) and 3,3',4,4',5-pentachlorobiphenyl (PCB126) relative to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). Whole-genome microarrays were used to evaluate the hepatic gene expression potency of TCDF and PCB126 relative to TCDD with complementary histopathology, tissue level analysis, and ethoxyresorufin-O-deethylase (EROD) assay results. Immature ovariectomized C57BL/6 mice were gavaged with 0.001, 0.01, 0.03, 0.1, 0.3, 1, 3, 10, 30, 100, and 300 μg/kg TCDD and TEF-adjusted doses (TEF for TCDF and PCB126 is 0.1) of TCDF or PCB126 (1, 3, 10, 30, 100, 300, 1000, and 3000 μg/kg of TCDF or PCB126) or sesame oil vehicle and sacrificed 24 h post dose. In general, TCDD, TCDF, and PCB126 tissue levels, as well as histopathological effects, were comparable when comparing TEF-adjusted doses. Automated dose-response modeling (ToxResponse Modeler) of the microarray data identified 210 TCDF and 40 PCB126 genes that exhibited sigmoidal dose-response curves with comparable slopes when compared with TCDD. These similar responses were used to calculate a median TCDF gene expression relative potency (REP) of 0.06 and a median PCB126 gene expression REP of 0.02. REPs of 0.02 were also calculated for EROD induction for both compounds. Collectively, these data suggest that differences in the ability of the liganded aryl hydrocarbon receptor:AhR nuclear translocator complex to elicit differential hepatic gene expression, in addition to pharmacokinetic differences between ligands, influence their potency in immature ovariectomized C57BL/6 mice.

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.

PCB153-elicited hepatic responses in the immature, ovariectomized C57BL/6 mice: comparative toxicogenomic effects of dioxin and non-dioxin-like ligands

Polychlorinated biphenyls (PCBs) are ubiquitous contaminants found as complex mixtures of coplanar and non-coplanar congeners. The hepatic temporal and dose-dependent effects of the most abundant non-dioxin-like congener, 2,2',4,4',5,5'-hexachlorobiphenyl (PCB153), were examined in immature, ovariectomized C57BL/6 mice, and compared to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), the prototypical aryl hydrocarbon receptor (AhR) ligand. Animals were gavaged once with 300 mg/kg PCB153 or sesame oil vehicle and sacrificed 4, 12, 24, 72 or 168 h post dose. In the dose-response study, mice were gavaged with 1, 3, 10, 30, 100 or 300 mg/kg PCB153 or sesame oil for 24 h. Significant increases in relative liver weights were induced with 300 mg/kg PCB153 between 24 and 168 h, accompanied by slight vacuolization and hepatocellular hypertrophy. The hepatic differential expression of 186 and 177 genes was detected using Agilent 4 x 44 K microarrays in the time course (|fold change|> or =1.5, P1(t)> or =0.999) and dose-response (|fold change|> or =1.5, P1(t)> or =0.985) studies, respectively. Comparative analysis with TCDD suggests that the differential gene expression elicited by PCB153 was not mediated by the AhR. Furthermore, constitutive androstane and pregnane X receptor (CAR/PXR) regulated genes including Cyp2b10, Cyp3a11, Ces2, Insig2 and Abcc3 were dose-dependently induced by PCB153. Collectively, these results suggest that the hepatocellular effects elicited by PCB153 are qualitatively and quantitatively different from TCDD and suggestive of CAR/PXR regulation.

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.

Automated dose-response analysis of the relative hepatic gene expression potency of TCDF in C57BL/6 mice

Toxic equivalency factors (TEFs) are assigned to dioxin-like chemicals based on relative potency (REP) values of individual adaptive and toxic responses compared to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). Agilent 4x44K oligonucleotide microarrays were used to examine the hepatic gene expression potency of 2,3,7,8-tetrachlorodibenzofuran (TCDF), relative to TCDD with complementary histopathology, TCDD and TCDF tissue level analysis, and ethoxyresorufin-O-deethylase (EROD) assay data. Immature ovariectomized C57BL/6 mice were gavaged with 0.03, 0.1, 0.3, 1, 3, 10, 30, or 100 microg/kg TCDD, the World Health Organization TEF-adjusted doses (10 x TCDD dose) of TCDF (0.3, 1, 3, 10, 30, 100, or 300 microg/kg), or sesame oil vehicle and killed at 72 h. Two thousand two hundred eighty-eight and 1347 genes were differentially expressed (P1(t) > 0.90) at one or more doses by TCDD and TCDF, respectively. Automated dose-response modeling (ToxResponse Modeler) identified a total of 1027 and 837 genes with either a sigmoidal, exponential, linear, Gaussian, or quadratic dose-response relationship 72 h after treatment in TCDD and TCDF, respectively. Two hundred seventy genes exhibited a sigmoidal TCDD-induced dose-response (ED(50s) from 0.08 to 42.2 microg/kg) compared to only 179 sigmoidal responsive genes (ED(50s) from 0.74 to 299.9 microg/kg) elicited by TCDF. Of the 1027 TCDD dose-responsive genes, 654 were not examined further due to the lack of a dose response elicited by TCDF. Of the 373 genes that exhibited a TCDD and TCDF dose response, REPs were calculated for the 83 genes that exhibited comparable sigmoidal curve shapes and slopes. The median REP for these 83 genes was 0.10, with a maximum REP of 0.56 and a minimum of 0.01. REPs of 0.04 were also calculated for EROD and increase in relative liver weight (RLW) at 72 h. Collectively, the lower number of TCDF-induced genes compared to TCDD and the 0.04 REPs for EROD activity and increased RLW are not consistent with the TEF of 0.10 for the hepatotoxicity of TCDF in C57BL/6 mice at 72 h.

Toxicogenomic biomarkers for liver toxicity

Toxicogenomics (TGx) is a widely used technique in the preclinical stage of drug development to investigate the molecular mechanisms of toxicity. A number of candidate TGx biomarkers have now been identified and are utilized for both assessing and predicting toxicities. Further accumulation of novel TGx biomarkers will lead to more efficient, appropriate and cost effective drug risk assessment, reinforcing the paradigm of the conventional toxicology system with a more profound understanding of the molecular mechanisms of drug-induced toxicity. In this paper, we overview some practical strategies as well as obstacles for identifying and utilizing TGx biomarkers based on microarray analysis. Since clinical hepatotoxicity is one of the major causes of drug development attrition, the liver has been the best documented target organ for TGx studies to date, and we therefore focused on information from liver TGx studies. In this review, we summarize the current resources in the literature in regard to TGx studies of the liver, from which toxicologists could extract potential TGx biomarker gene sets for better hepatotoxicity risk assessment.