Inhibition of the urea cycle by the environmental contaminant 2,3,7,8-tetrachlorodibenzo-p-dioxin increases serum ammonia levels in mice.

The aryl hydrocarbon receptor is a ligand-activated transcription factor known for mediating the effects of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) and related compounds. TCDD induces nonalcoholic fatty liver disease (NAFLD)-like pathologies including simple steatosis that can progress to steatohepatitis with fibrosis and bile duct proliferation in male mice. Dose-dependent progression of steatosis to steatohepatitis with fibrosis by TCDD has been associated with metabolic reprogramming, including the disruption of amino acid metabolism. Here, we used targeted metabolomic analysis to reveal dose-dependent changes in the level of ten serum and eleven hepatic amino acids in mice upon treatment with TCDD. Bulk RNA-seq and protein analysis showed TCDD repressed CPS1, OTS, ASS1, ASL, and GLUL, all of which are associated with the urea cycle and glutamine biosynthesis. Urea and glutamine are end products of the detoxification and excretion of ammonia, a toxic byproduct of amino acid catabolism. Furthermore, we found that the catalytic activity of OTC, a rate-limiting step in the urea cycle was also dose dependently repressed. These results are consistent with an increase in circulating ammonia. Collectively, the repression of the urea and glutamate-glutamine cycles increased circulating ammonia levels and the toxicity of TCDD.

Generative modeling of single-cell gene expression for dose-dependent chemical perturbations.

Single-cell sequencing reveals the heterogeneity of cellular response to chemical perturbations. However, testing all relevant combinations of cell types, chemicals, and doses is a daunting task. A deep generative learning formalism called variational autoencoders (VAEs) has been effective in predicting single-cell gene expression perturbations for single doses. Here, we introduce single-cell variational inference of dose-response (scVIDR), a VAE-based model that predicts both single-dose and multiple-dose cellular responses better than existing models. We show that scVIDR can predict dose-dependent gene expression across mouse hepatocytes, human blood cells, and cancer cell lines. We biologically interpret the latent space of scVIDR using a regression model and use scVIDR to order individual cells based on their sensitivity to chemical perturbation by assigning each cell a "pseudo-dose" value. We envision that scVIDR can help reduce the need for repeated animal testing across tissues, chemicals, and doses.

A Case for Accelerating Standards to Achieve the FAIR Principles of Environmental Health Research Experimental Data.

BACKGROUND: Funding agencies, publishers, and other stakeholders are pushing environmental health science investigators to improve data sharing; to promote the findable, accessible, interoperable, and reusable (FAIR) principles; and to increase the rigor and reproducibility of the data collected. Accomplishing these goals will require significant cultural shifts surrounding data management and strategies to develop robust and reliable resources that bridge the technical challenges and gaps in expertise. OBJECTIVE: In this commentary, we examine the current state of managing data and metadata-referred to collectively as (meta)data-in the experimental environmental health sciences. We introduce new tools and resources based on in vivo experiments to serve as examples for the broader field. METHODS: We discuss previous and ongoing efforts to improve (meta)data collection and curation. These include global efforts by the Functional Genomics Data Society to develop metadata collection tools such as the Investigation, Study, Assay (ISA) framework, and the Center for Expanded Data Annotation and Retrieval. We also conduct a case study of in vivo data deposited in the Gene Expression Omnibus that demonstrates the current state of in vivo environmental health data and highlights the value of using the tools we propose to support data deposition. DISCUSSION: The environmental health science community has played a key role in efforts to achieve the goals of the FAIR guiding principles and is well positioned to advance them further. We present a proposed framework to further promote these objectives and minimize the obstacles between data producers and data scientists to maximize the return on research investments. https://doi.org/10.1289/EHP11484.

Consequences of reprogramming acetyl-CoA metabolism by 2,3,7,8-tetrachlorodibenzo-p-dioxin in the mouse liver.

2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) is a persistent environmental contaminant that induces the progression of steatosis to steatohepatitis with fibrosis in mice. Furthermore, TCDD reprograms hepatic metabolism by redirecting glycolytic intermediates while inhibiting lipid metabolism. Here, we examined the effect of TCDD on hepatic acetyl-coenzyme A (acetyl-CoA) and ß-hydroxybutyrate levels as well as protein acetylation and ß-hydroxybutyrylation. Acetyl-CoA is not only a central metabolite in multiple anabolic and catabolic pathways, but also a substrate used for posttranslational modification of proteins and a surrogate indicator of cellular energy status. Targeted metabolomic analysis revealed a dose-dependent decrease in hepatic acetyl-CoA levels coincident with the phosphorylation of pyruvate dehydrogenase (E1), and the induction of pyruvate dehydrogenase kinase 4 and pyruvate dehydrogenase phosphatase, while repressing ATP citrate lyase and short-chain acyl-CoA synthetase gene expression. In addition, TCDD dose-dependently reduced the levels of hepatic ß-hydroxybutyrate and repressed ketone body biosynthesis gene expression. Moreover, levels of total hepatic protein acetylation and ß-hydroxybutyrylation were reduced. AMPK phosphorylation was induced consistent with acetyl-CoA serving as a cellular energy status surrogate, yet subsequent targets associated with re-establishing energy homeostasis were not activated. Collectively, TCDD reduced hepatic acetyl-CoA and ß-hydroxybutyrate levels eliciting starvation-like conditions despite normal levels of food intake.

Single-cell transcriptomics shows dose-dependent disruption of hepatic zonation by TCDD in mice.

2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) dose-dependently induces the development of hepatic fat accumulation and inflammation with fibrosis in mice initially in the portal region. Conversely, differential gene and protein expression is first detected in the central region. To further investigate cell-specific and spatially resolved dose-dependent changes in gene expression elicited by TCDD, single-nuclei RNA sequencing and spatial transcriptomics were used for livers of male mice gavaged with TCDD every 4 days for 28 days. The proportion of 11 cell (sub)types across 131 613 nuclei dose-dependently changed with 68% of all portal and central hepatocyte nuclei in control mice being overtaken by macrophages following TCDD treatment. We identified 368 (portal fibroblasts) to 1339 (macrophages) differentially expressed genes. Spatial analyses revealed initial loss of portal identity that eventually spanned the entire liver lobule with increasing dose. Induction of R-spondin 3 (Rspo3) and pericentral Apc, suggested dysregulation of the Wnt/ß-catenin signaling cascade in zonally resolved steatosis. Collectively, the integrated results suggest disruption of zonation contributes to the pattern of TCDD-elicited NAFLD pathologies.

A Physiologically Based Pharmacokinetic (PBPK) Modeling Framework for Mixtures of Dioxin-like Compounds.

Humans are exposed to persistent organic pollutants, such as dioxin-like compounds (DLCs), as mixtures. Understanding and predicting the toxicokinetics and thus internal burden of major constituents of a DLC mixture is important for assessing their contributions to health risks. PBPK models, including dioxin models, traditionally focus on one or a small number of compounds; developing new or extending existing models for mixtures often requires tedious, error-prone coding work. This lack of efficiency to scale up for multi-compound exposures is a major technical barrier toward large-scale mixture PBPK simulations. Congeners in the DLC family, including 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), share similar albeit quantitatively different toxicokinetic and toxicodynamic properties. Taking advantage of these similarities, here we reported the development of a human PBPK modeling framework for DLC mixtures that can flexibly accommodate an arbitrary number of congeners. Adapted from existing TCDD models, our mixture model contains the blood and three diffusion-limited compartments-liver, fat, and rest of the body. Depending on the number of congeners in a mixture, varying-length vectors of ordinary differential equations (ODEs) are automatically generated to track the tissue concentrations of the congeners. Shared ODEs are used to account for common variables, including the aryl hydrocarbon receptor (AHR) and CYP1A2, to which the congeners compete for binding. Binary and multi-congener mixture simulations showed that the AHR-mediated cross-induction of CYP1A2 accelerates the sequestration and metabolism of DLC congeners, resulting in consistently lower tissue burdens than in single exposure, except for the liver. Using dietary intake data to simulate lifetime exposures to DLC mixtures, the model demonstrated that the relative contributions of individual congeners to blood or tissue toxic equivalency (TEQ) values are markedly different than those to intake TEQ. In summary, we developed a mixture PBPK modeling framework for DLCs that may be utilized upon further improvement as a quantitative tool to estimate tissue dosimetry and health risks of DLC mixtures.

Dioxin-elicited decrease in cobalamin redirects propionyl-CoA metabolism to the ß-oxidation-like pathway resulting in acrylyl-CoA conjugate buildup.

2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) is a persistent environmental contaminant that induces diverse biological and toxic effects, including reprogramming intermediate metabolism, mediated by the aryl hydrocarbon receptor. However, the specific reprogramming effects of TCDD are unclear. Here, we performed targeted LC-MS analysis of hepatic extracts from mice gavaged with TCDD. We detected an increase in S-(2-carboxyethyl)-L-cysteine, a conjugate from the spontaneous reaction between the cysteine sulfhydryl group and highly reactive acrylyl-CoA, an intermediate in the cobalamin (Cbl)-independent ß-oxidation-like metabolism of propionyl-CoA. TCDD repressed genes in both the canonical Cbl-dependent carboxylase and the alternate Cbl-independent ß-oxidation-like pathways as well as inhibited methylmalonyl-CoA mutase (MUT) at lower doses. Moreover, TCDD decreased serum Cbl levels and hepatic cobalt levels while eliciting negligible effects on gene expression associated with Cbl absorption, transport, trafficking, or derivatization to 5'-deoxy-adenosylcobalamin (AdoCbl), the required MUT cofactor. Additionally, TCDD induced the gene encoding aconitate decarboxylase 1 (Acod1), the enzyme responsible for decarboxylation of cis-aconitate to itaconate, and dose-dependently increased itaconate levels in hepatic extracts. Our results indicate MUT inhibition is consistent with itaconate activation to itaconyl-CoA, a MUT suicide inactivator that forms an adduct with adenosylcobalamin. This adduct in turn inhibits MUT activity and reduces Cbl levels. Collectively, these results suggest the decrease in MUT activity is due to Cbl depletion following TCDD treatment, which redirects propionyl-CoA metabolism to the alternate Cbl-independent ß-oxidation-like pathway. The resulting hepatic accumulation of acrylyl-CoA likely contributes to TCDD-elicited hepatotoxicity and the multihit progression of steatosis to steatohepatitis with fibrosis.

Genome-Wide ChIPseq Analysis of AhR, COUP-TF, and HNF4 Enrichment in TCDD-Treated Mouse Liver.

The aryl hydrocarbon receptor (AhR) is a ligand-activated transcription factor known for mediating the toxicity of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) and related compounds. Although the canonical mechanism of AhR activation involves heterodimerization with the aryl hydrocarbon receptor nuclear translocator, other transcriptional regulators that interact with AhR have been identified. Enrichment analysis of motifs in AhR-bound genomic regions implicated co-operation with COUP transcription factor (COUP-TF) and hepatocyte nuclear factor 4 (HNF4). The present study investigated AhR, HNF4α and COUP-TFII genomic binding and effects on gene expression associated with liver-specific function and cell differentiation in response to TCDD. Hepatic ChIPseq data from male C57BL/6 mice at 2 h after oral gavage with 30 µg/kg TCDD were integrated with bulk RNA-sequencing (RNAseq) time-course (2-72 h) and dose-response (0.01-30 µg/kg) datasets to assess putative AhR, HNF4α and COUP-TFII interactions associated with differential gene expression. Functional enrichment analysis of differentially expressed genes (DEGs) identified differential binding enrichment for AhR, COUP-TFII, and HNF4α to regions within liver-specific genes, suggesting intersections associated with the loss of liver-specific functions and hepatocyte differentiation. Analysis found that the repression of liver-specific, HNF4α target and hepatocyte differentiation genes, involved increased AhR and HNF4α binding with decreased COUP-TFII binding. Collectively, these results suggested TCDD-elicited loss of liver-specific functions and markers of hepatocyte differentiation involved interactions between AhR, COUP-TFII and HNF4α.

Benchmarking of a Bayesian single cell RNAseq differential gene expression test for dose-response study designs.

The application of single-cell RNA sequencing (scRNAseq) for the evaluation of chemicals, drugs, and food contaminants presents the opportunity to consider cellular heterogeneity in pharmacological and toxicological responses. Current differential gene expression analysis (DGEA) methods focus primarily on two group comparisons, not multi-group dose-response study designs used in safety assessments. To benchmark DGEA methods for dose-response scRNAseq experiments, we proposed a multiplicity corrected Bayesian testing approach and compare it against 8 other methods including two frequentist fit-for-purpose tests using simulated and experimental data. Our Bayesian test method outperformed all other tests for a broad range of accuracy metrics including control of false positive error rates. Most notable, the fit-for-purpose and standard multiple group DGEA methods were superior to the two group scRNAseq methods for dose-response study designs. Collectively, our benchmarking of DGEA methods demonstrates the importance in considering study design when determining the most appropriate test methods.

Genetics-Based Approach to Identify Novel Genes Regulated by the Aryl Hydrocarbon Receptor in Mouse Liver.

The aryl hydrocarbon receptor (AHR) is a ligand-activated transcription factor in the Per-Arnt-Sim superfamily of environmental sensors that is linked to several metabolic diseases, including nonalcoholic fatty liver disease. Much remains unknown regarding the impact of genetic variation in AHR-driven disease, as past studies have focused on a small number of inbred strains. Recently, the presence of a wide range of interindividual variability amongst humans was reported in response to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), the prototypical ligand of the AHR. In this study, a panel of 14 diverse mouse strains was exposed to TCDD for 10 days to characterize the AHR-mediated response across genetic backgrounds. Responses to TCDD are heavily dependent on genetic background. Although mice carry 1 of 4 Ahr alleles known to impact the affinity to AHR-ligands, we observed significant intra-allelic variability suggesting the presence of novel genetic modifiers of AHR signaling. A regression-based approach was used to scan for genes regulated by the AHR and/or associated with TCDD-induced phenotypes. The approach identified 7 genes, 2 of which are novel, that are likely regulated by the AHR based on association with hepatic TCDD burden (p ≤ .05). Finally, we identified 1 gene, Dio1, which was associated with change in percent body fat across the diverse set of strains (p ≤ .05). Overall, the results in this study exemplify the power of genetics-based approaches in identifying novel genes that are putatively regulated by the AHR.

Single-Nuclei RNA Sequencing Assessment of the Hepatic Effects of 2,3,7,8-Tetrachlorodibenzo-p-dioxin.

BACKGROUND AND AIMS: Characterization of cell specific transcriptional responses to hepatotoxicants is lost in the averages of bulk RNA-sequencing (RNA-seq). Single-cell/nuclei RNA-seq technologies enable the transcriptomes of individual cell (sub)types to be assessed within the context of in vivo models. METHODS: Single-nuclei RNA-sequencing (snSeq) of frozen liver samples from male C57BL/6 mice gavaged with sesame oil vehicle or 30 µg/kg 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) every 4 days for 28 days was used to demonstrate the application of snSeq for the evaluation of xenobiotics. RESULTS: A total of 19,907 genes were detected across 16,015 nuclei from control and TCDD-treated livers. Eleven cell (sub)types reflected the expected cell diversity of the liver including distinct pericentral, midzonal, and periportal hepatocyte subpopulations. TCDD altered relative proportions of cell types and elicited cell-specific gene expression profiles. For example, macrophages increased from 0.5% to 24.7%, while neutrophils were only present in treated samples, consistent with histological evaluation. The number of differentially expressed genes (DEGs) in each cell type ranged from 122 (cholangiocytes) to 7625 (midcentral hepatocytes), and loosely correlated with the basal expression level of Ahr, the canonical mediator of TCDD and related compounds. In addition to the expected functions within each cell (sub)types, RAS signaling and related pathways were specifically enriched in nonparenchymal cells while metabolic process enrichment occurred primarily in hepatocytes. snSeq also identified the expansion of a Kupffer cell subtype highly expressing Gpnmb, as reported in a dietary NASH model. CONCLUSIONS: We show that snSeq of frozen liver samples can be used to assess cell-specific transcriptional changes and population shifts in models of hepatotoxicity when examining freshly isolated cells is not feasible.

2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) dysregulates hepatic one carbon metabolism during the progression of steatosis to steatohepatitis with fibrosis in mice.

2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD), a persistent environmental contaminant, induces steatosis that can progress to steatohepatitis with fibrosis, pathologies that parallel stages in the development of non-alcoholic fatty liver disease (NAFLD). Coincidently, one carbon metabolism (OCM) gene expression and metabolites are often altered during NAFLD progression. In this study, the time- and dose-dependent effects of TCDD were examined on hepatic OCM in mice. Despite AhR ChIP-seq enrichment at 2 h, OCM gene expression was not changed within 72 h following a bolus dose of TCDD. Dose-dependent repression of methionine adenosyltransferase 1A (Mat1a), adenosylhomocysteinase (Achy) and betaine-homocysteine S-methyltransferase (Bhmt) mRNA and protein levels following repeated treatments were greater at 28 days compared to 8 days. Accordingly, levels of methionine, betaine, and homocysteic acid were dose-dependently increased, while S-adenosylmethionine, S-adenosylhomocysteine, and cystathionine exhibited non-monotonic dose-dependent responses consistent with regulation by OCM intermediates and repression of glycine N-methyltransferase (Gnmt). However, the dose-dependent effects on SAM-dependent metabolism of polyamines and creatine could not be directly attributed to alterations in SAM levels. Collectively, these results demonstrate persistent AhR activation disrupts hepatic OCM metabolism at the transcript, protein and metabolite levels within context of TCDD-elicited progression of steatosis to steatohepatitis with fibrosis.

2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) alters hepatic polyunsaturated fatty acid metabolism and eicosanoid biosynthesis in female Sprague-Dawley rats.

2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) is a potent aryl hydrocarbon receptor (AhR) agonist that elicits a broad spectrum of dose-dependent hepatic effects including lipid accumulation, inflammation, and fibrosis. To determine the role of inflammatory lipid mediators in TCDD-mediated hepatotoxicity, eicosanoid metabolism was investigated. Female Sprague-Dawley (SD) rats were orally gavaged with sesame oil vehicle or 0.01-10 µg/kg TCDD every 4 days for 28 days. Hepatic RNA-Seq data was integrated with untargeted metabolomics of liver, serum, and urine, revealing dose-dependent changes in linoleic acid (LA) and arachidonic acid (AA) metabolism. TCDD also elicited dose-dependent differential gene expression associated with the cyclooxygenase, lipoxygenase, and cytochrome P450 epoxidation/hydroxylation pathways with corresponding changes in ω-6 (e.g. AA and LA) and ω-3 polyunsaturated fatty acids (PUFAs), as well as associated eicosanoid metabolites. Overall, TCDD increased the ratio of ω-6 to ω-3 PUFAs. Phospholipase A2 (Pla2g12a) was induced consistent with increased AA metabolism, while AA utilization by induced lipoxygenases Alox5 and Alox15 increased leukotrienes (LTs). More specifically, TCDD increased pro-inflammatory eicosanoids including leukotriene LTB4, and LTB3, known to recruit neutrophils to damaged tissue. Dose-response modeling suggests the cytochrome P450 hydroxylase/epoxygenase and lipoxygenase pathways are more sensitive to TCDD than the cyclooxygenase pathway. Hepatic AhR ChIP-Seq analysis found little enrichment within the regulatory regions of differentially expressed genes (DEGs) involved in eicosanoid biosynthesis, suggesting TCDD-elicited dysregulation of eicosanoid metabolism is a downstream effect of AhR activation. Overall, these results suggest alterations in eicosanoid metabolism may play a key role in TCDD-elicited hepatotoxicity associated with the progression of steatosis to steatohepatitis.

A toxicogenomic approach for the risk assessment of the food contaminant acetamide.

Acetamide (CAS 60-35-5) is detected in common foods. Chronic rodent bioassays led to its classification as a group 2B possible human carcinogen due to the induction of liver tumors in rats. We used a toxicogenomics approach in Wistar rats gavaged daily for 7 or 28 days at doses of 300 to 1500 mg/kg/day (mkd) to determine a point of departure (POD) and investigate its mode of action (MoA). Ki67 labeling was increased at doses ≥750 mkd up to 3.3-fold representing the most sensitive apical endpoint. Differential gene expression analysis by RNA-Seq identified 1110 and 1814 differentially expressed genes in male and female rats, respectively, following 28 days of treatment. Down-regulated genes were associated with lipid metabolism while up-regulated genes included cell signaling, immune response, and cell cycle functions. Benchmark dose (BMD) modeling of the Ki67 labeling index determined the BMD10 lower confidence limit (BMDL10) as 190 mkd. Transcriptional BMD modeling revealed excellent concordance between transcriptional POD and apical endpoints. Collectively, these results indicate that acetamide is most likely acting through a mitogenic MoA, though specific key initiating molecular events could not be elucidated. A POD value of 190 mkd determined for cell proliferation is suggested for risk assessment purposes.

Comparison of Hepatic NRF2 and Aryl Hydrocarbon Receptor Binding in 2,3,7,8-Tetrachlorodibenzo-p-dioxin-Treated Mice Demonstrates NRF2-Independent PKM2 Induction.

2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) induces hepatic oxidative stress following activation of the aryl hydrocarbon receptor (AhR). Our recent studies showed TCDD induced pyruvate kinase muscle isoform 2 (Pkm2) as a novel antioxidant response in normal differentiated hepatocytes. To investigate cooperative regulation between nuclear factor, erythroid derived 2, like 2 (Nrf2) and the AhR in the induction of Pkm2, hepatic chromatin immunoprecipitation sequencing (ChIP-seq) analyses were integrated with RNA sequencing (RNA-seq) time-course data from mice treated with TCDD for 2-168 hours. ChIP-seq analysis 2 hours after TCDD treatment identified genome-wide NRF2 enrichment. Approximately 842 NRF2-enriched regions were located in the regulatory region of differentially expressed genes (DEGs), whereas 579 DEGs showed both NRF2 and AhR enrichment. Sequence analysis of regions with overlapping NRF2 and AhR enrichment showed over-representation of either antioxidant or dioxin response elements, although 18 possessed both motifs. NRF2 exhibited negligible enrichment within a closed Pkm chromatin region, whereas the AhR was enriched 29-fold. Furthermore, TCDD induced Pkm2 in primary hepatocytes from wild-type and Nrf2-null mice, indicating NRF2 is not required. Although NRF2 and AhR cooperate to regulate numerous antioxidant gene expression responses, the induction of Pkm2 by TCDD is independent of reactive oxygen species-mediated NRF2 activation.

Loss of liver-specific and sexually dimorphic gene expression by aryl hydrocarbon receptor activation in C57BL/6 mice.

The aryl hydrocarbon receptor (AhR) is a highly conserved transcription factor that mediates a broad spectrum of species-, strain-, sex-, age-, tissue-, and cell-specific responses elicited by structurally diverse ligands including 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). Dose-dependent effects on liver-specific and sexually dimorphic gene expression were examined in male and female mice gavaged with TCDD every 4 days for 28 or 92 days. RNA-seq data revealed the coordinated repression of 181 genes predominately expressed in the liver including albumin (3.7-fold), α-fibrinogen (14.5-fold), and ß-fibrinogen (17.4-fold) in males with corresponding AhR enrichment at 2 hr. Liver-specific genes exhibiting sexually dimorphic expression also demonstrated diminished divergence between sexes. For example, male-biased Gstp1 was repressed 3.0-fold in males and induced 4.5-fold in females, which were confirmed at the protein level. Disrupted regulation is consistent with impaired GHR-JAK2-STAT5 signaling and inhibition of female specific CUX2-mediated transcription as well as the repression of other key transcriptional regulators including Ghr, Stat5b, Bcl6, Hnf4a, Hnf6, Foxa1/2/3, and Zhx2. Attenuated liver-specific and sexually dimorphic gene expression was concurrent with the induction of fetal genes such as alpha-fetoprotein. The results suggest AhR activation causes the loss of liver-specific and sexually dimorphic gene expression producing a functionally "de-differentiated" hepatic phenotype.

Dose-Dependent Metabolic Reprogramming and Differential Gene Expression in TCDD-Elicited Hepatic Fibrosis.

We have previously shown that in response to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD)-elicited NAFLD progression, central carbon, glutaminolysis, and serine/folate metabolism are reprogrammed to support NADPH production and ROS defenses. To further investigate underlying dose-dependent responses associated with TCDD-induced fibrosis, female C57BL/6 mice were gavaged with TCDD every 4 days (d) for 28 d or 92 d. RNA-Seq, ChIP-Seq (2 h), and 28 d metabolomic (urine, serum, and hepatic extract) analyses were conducted with complementary serum marker assessments at 92 d. Additional vehicle and 30 µg/kg treatment groups were allowed to recover for 36 d following the 92-d treatment regimen to examine recovery from TCDD-elicited fibrosis. Histopathology revealed dose-dependent increases in hepatic fat accumulation, inflammation, and periportal collagen deposition at 92 days, with increased fibrotic severity in the recovery group. Serum proinflammatory and profibrotic interleukins-1ß, -2, -4, -6, and -10, as well as TNF-α and IFN-γ, exhibited dose-dependent induction. An increase in glucose tolerance was observed with a concomitant 3.0-fold decrease in hepatic glycogen linked to increased ascorbic acid biosynthesis and proline metabolism, consistent with increased fibrosis. RNA-Seq identified differential expression of numerous matrisome genes including an 8.8-fold increase in Tgfb2 indicating myofibroblast activation. Further analysis suggests reprogramming of glycogen, ascorbic acid, and amino acid metabolism in support of collagen deposition and the use of proline as a substrate for ATP production via the proline cycle. In summary, we demonstrate that glycogen, ascorbic acid, and amino acid metabolism are also reorganized to support remodeling of the extracellular matrix, progressing to hepatic fibrosis in response to chronic injury from TCDD.

RNA-Seq versus oligonucleotide array assessment of dose-dependent TCDD-elicited hepatic gene expression in mice.

BACKGROUND: Dose-dependent differential gene expression provides critical information required for regulatory decision-making. The lower costs associated with RNA-Seq have made it the preferred technology for transcriptomic analysis. However, concordance between RNA-Seq and microarray analyses in dose response studies has not been adequately vetted. RESULTS: We compared the hepatic transcriptome of C57BL/6 mice following gavage with sesame oil vehicle, 0.01, 0.03, 0.1, 0.3, 1, 3, 10, or 30 µg/kg TCDD every 4 days for 28 days using Illumina HiSeq RNA-Sequencing (RNA-Seq) and Agilent 4 × 44 K microarrays using the same normalization and analysis approach. RNA-Seq and microarray analysis identified a total of 18,063 and 16,403 genes, respectively, that were expressed in the liver. RNA-Seq analysis for differentially expressed genes (DEGs) varied dramatically depending on the P1(t) cut-off while microarray results varied more based on the fold change criteria, although responses strongly correlated. Verification by WaferGen SmartChip QRTPCR revealed that RNA-Seq had a false discovery rate of 24% compared to 54% for microarray analysis. Dose-response modeling of RNA-Seq and microarray data demonstrated similar point of departure (POD) and ED50 estimates for common DEGs. CONCLUSIONS: There was a strong correspondence between RNA-Seq and Agilent array transcriptome profiling when using the same samples and analysis strategy. However, RNA-Seq provided superior quantitative data, identifying more genes and DEGs, as well as qualitative information regarding identity and annotation for dose response modeling in support of regulatory decision-making.

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.

Use of Xenopus laevis as a model for investigating in vitro and in vivo endocrine disruption in amphibians.

The estrogenic activity of 17beta-estradiol (E2), alpha-zearalenol (alpha-ZEA), genistein (GEN), and 4-t-octylphenol (4-t-OP) was investigated using Xenopus laevis-based assays. All test compounds competed with [3H]E2 for binding to a recombinant Xenopus estrogen receptor (xER) with the following relative affinities: E2 > alpha-ZEA > 4-t-OP > GEN. The ability of these compounds to induce xER-mediated reporter gene expression was then assessed in MCF-7 human breast cancer cells cotransfected with a Gal4-xERdef chimeric estrogen receptor and a Gal4-regulated luciferase reporter gene. Luciferase activity was increased 30- to 50-fold by 10 nM E2 relative to that in solvent control. Maximal reporter gene activity induced by 10 nM alpha-ZEA was 54% of that induced by E2; however, the activity did not increase following doses of up to 10 microM. A dose of 1 microM 4-t-OP induced 23% of the maximal reporter gene activity induced by E2, whereas 10 microM GEN induced activity to the same level as E2. A dose-dependent increase in vitellogenin (VTG) mRNA expression was observed in Xenopus treated intraperitoneally with E2 at 0.05 to 5 mg/kg/d for three consecutive days, with the maximal induction observed in the group receiving 1 mg/kg/d. The alpha-ZEA, GEN, and 4-t-OP also significantly induced VTG mRNA expression, although at higher doses. These results demonstrate the utility of X. laevis as an amphibian model to assess the estrogenic activity of endocrine disruptors.