Cell-specific AHR-driven differential gene expression in the mouse liver cell following acute TCDD exposure.

2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) is a persistent environmental contaminant that disrupts hepatic function leading to steatotic liver disease (SLD)-like pathologies, such as steatosis, steatohepatitis, and fibrosis. These effects are mediated by the aryl hydrocarbon receptor following changes in gene expression. Although diverse cell types are involved, initial cell-specific changes in gene expression have not been reported. In this study, differential gene expression in hepatic cell types was examined in male C57BL/6 mice gavaged with 30 µg/kg of TCDD using single-nuclei RNA-sequencing. Ten liver cell types were identified with the proportions of most cell types remaining unchanged, except for neutrophils which increased at 72 h. Gene expression suggests TCDD induced genes related to oxidative stress in hepatocytes as early as 2 h. Lipid homeostasis was disrupted in hepatocytes, macrophages, B cells, and T cells, characterized by the induction of genes associated with lipid transport, steroid hormone biosynthesis, and the suppression of ß-oxidation, while linoleic acid metabolism was altered in hepatic stellate cells (HSCs), B cells, portal fibroblasts, and plasmacytoid dendritic cells. Pro-fibrogenic processes were also enriched, including the induction retinol metabolism genes in HSCs and the early induction of anti-fibrolysis genes in hepatocytes, endothelial cells, HSCs, and macrophages. Hepatocytes also had gene expression changes consistent with hepatocellular carcinoma. Collectively, these findings underscore the effects of TCDD in initiating SLD-like phenotypes and identified cell-specific gene expression changes related to oxidative stress, steatosis, fibrosis, cell proliferation and the development of HCC.

Integrins play a role in stress relaxation of perivascular adipose tissue.

Perivascular adipose tissue (PVAT) is known for being anti-contractile in healthy tissues. We discovered a new function of PVAT, the ability to stress relax and maintain a tone in response to a stretch. This is of note because stress relaxation has been attributed to smooth muscle, of which PVAT has none that is organized in a functional layer. We test the hypothesis the interactions of integrins with collagen play a role in stress relaxation. Our model is the thoracic aorta of the male Dahl SS rat. The PVAT and aorta were physically separated for most assays. Results from single nuclei RNA sequencing (snRNAseq) experiments, histochemistry and isometric contractility were also used. Masson Trichrome staining made evident the expression of collagen in PVAT. From snRNA seq experiments of the PVAT, mRNA for multiple collagen and integrin isoforms were detected: the α1 and ß1 integrin were most highly expressed. Pharmacological inhibition of integrin/collagen interaction was effected by the specific α1ß1 distintegrin obtustatin or general integrin inhibitor RGD peptide. RGD peptide but not obtustatin increased the stress relaxation. Cell-cell communication inference identified integrins αv and α5, two major RGD motif containing isoforms, as potential signaling partners of collagens. Collectively, these findings validate that stress relaxation can occur in a non-smooth muscle tissue, doing so in part through integrin-collagen interactions that may not include α1ß1 heterodimers. The importance of this lies in considering PVAT as a vascular layer that possesses mechanical functions.

Early life adversity promotes gastrointestinal dysfunction through a sex-dependent phenotypic switch in enteric glia.

Irritable bowel syndrome and related disorders of gut-brain interaction (DGBI) are common and exhibit a complex, poorly understood etiology that manifests as abnormal gut motility and pain. Risk factors such as biological sex, stressors during critical periods, and inflammation are thought to influence DGBI vulnerability by reprogramming gut-brain circuits, but the specific cells affected are unclear. Here, we used a model of early life stress to understand cellular mechanisms in the gut that produce DGBIs. Our findings identify enteric glia as a key cellular substrate in which stress and biological sex converge to dictate DGBI susceptibility. Enteric glia exhibit sexual dimorphism in genes and functions related to cellular communication, inflammation, and disease susceptibility. Experiencing early life stress has sex-specific effects on enteric glia that cause a phenotypic switch in male glia toward a phenotype normally observed in females. This phenotypic transformation is followed by physiological changes in the gut, mirroring those observed in DGBI in humans. These effects are mediated, in part, by alterations to glial prostaglandin and endocannabinoid signaling. Together, these data identify enteric glia as a cellular integration site through which DGBI risk factors produce changes in gut physiology and suggest that manipulating glial signaling may represent an attractive target for sex-specific therapeutic strategies in DGBIs.

A cell atlas of thoracic aortic perivascular adipose tissue: a focus on mechanotransducers.

Perivascular adipose tissue (PVAT) is increasingly recognized for its function in mechanotransduction. However, major gaps remain in our understanding of the cells present in PVAT, as well as how different cells contribute to mechanotransduction. We hypothesized that snRNA-seq would reveal the expression of mechanotransducers, and test one (PIEZO1) to illustrate the expression and functional agreement between single-nuclei RNA sequencing (snRNA-seq) and physiological measurements. To contrast two brown tissues, subscapular brown adipose tissue (BAT) was also examined. We used snRNA-seq of the thoracic aorta PVAT (taPVAT) and BAT from male Dahl salt-sensitive (Dahl SS) rats to investigate cell-specific expression mechanotransducers. Localization and function of the mechanostransducer PIEZO1 were further examined using immunohistochemistry (IHC) and RNAscope, as well as pharmacological antagonism. Approximately 30,000 nuclei from taPVAT and BAT each were characterized by snRNA-seq, identifying eight major cell types expected and one unexpected (nuclei with oligodendrocyte marker genes). Cell-specific differential gene expression analysis between taPVAT and BAT identified up to 511 genes (adipocytes) with many (≥20%) being unique to individual cell types. Piezo1 was the most highly, widely expressed mechanotransducer. The presence of PIEZO1 in the PVAT but not the adventitia was confirmed by RNAscope and IHC in male and female rats. Importantly, antagonism of PIEZO1 by GsMTX4 impaired the PVAT's ability to hold tension. Collectively, the cell compositions of taPVAT and BAT are highly similar, and PIEZO1 is likely a mechanotransducer in taPVAT.NEW & NOTEWORTHY This study describes the atlas of cells in the thoracic aorta perivascular adipose tissue (taPVAT) of the Dahl-SS rat, an important hypertension model. We show that mechanotransducers are widely expressed in these cells. Moreover, PIEZO1 expression is shown to be restricted to the taPVAT and is functionally implicated in stress relaxation. These data will serve as the foundation for future studies investigating the role of taPVAT in this model of hypertensive disease.

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.

Characterizing the impact of simvastatin co-treatment of cell specific TCDD-induced gene expression and systemic toxicity.

2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) is associated with metabolic syndrome (MetS) in humans and elicits pathologies in rodents that resemble non-alcoholic fatty liver disease (NAFLD) in humans through activation of the aryl hydrocarbon receptor (AHR) pathway. Dysregulation of cholesterol homeostasis, an aspect of MetS, is linked to NAFLD pathogenesis. TCDD exposure is also linked to the suppression of genes that encode key cholesterol biosynthesis steps and changes in serum cholesterol levels. In a previous experiment, treating mice with TCDD in the presence of simvastatin, a 3-Hydroxy-3-Methylglutaryl-CoA Reductase competitive inhibitor, altered lipid and glycogen levels, AHR-battery gene expression, and liver injury in male mice compared to TCDD alone. The aim of this study was to deduce a possible mechanism(s) for the metabolic changes and increased injury using single-nuclei RNA sequencing in mouse liver. We demonstrated that co-treated mice experienced wasting and increased AHR activation compared to TCDD alone. Furthermore, relative proportions of cell (sub)types were different between TCDD alone and co-treated mice including important mediators of NAFLD progression like hepatocytes and immune cell populations. Analysis of non-overlapping differentially expressed genes identified several pathways where simvastatin co-treatment significantly impacted TCDD-induced changes, which may explain the differences between treatments. Overall, these results demonstrate a connection between dysregulation of cholesterol homeostasis and toxicant-induced metabolic changes.

A Cell Atlas of Thoracic Aortic Perivascular Adipose Tissue: a focus on mechanotransducers.

Perivascular adipose tissue (PVAT) is increasingly recognized for its function in mechanotransduction. To examine the cell-specificity of recognized mechanotransducers we used single nuclei RNA sequencing (snRNAseq) of the thoracic aorta PVAT (taPVAT) from male Dahl SS rats compared to subscapular brown adipose tissue (BAT). Approximately 30,000 nuclei from taPVAT and BAT each were characterized by snRNAseq, identifying 8 major cell types expected and one unexpected (nuclei with oligodendrocyte marker genes). Cell-specific differential gene expression analysis between taPVAT and BAT identified up to 511 genes (adipocytes) with many (≥20%) being unique to individual cell types. Piezo1 was the most highly, widely expressed mechanotransducer. Presence of PIEZO1 in the PVAT was confirmed by RNAscope® and IHC; antagonism of PIEZO1 impaired the PVAT's ability to hold tension. Collectively, the cell compositions of taPVAT and BAT are highly similar, and PIEZO1 is likely a mechanotransducer in taPVAT.

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.

Cystine/Glutamate Xc- Antiporter Induction Compensates for Transsulfuration Pathway Repression by 2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) to Ensure Cysteine for Hepatic Glutathione Biosynthesis.

Exposure to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) has been associated with the induction of oxidative stress and the progression of steatosis to steatohepatitis with fibrosis. It also disrupts metabolic pathways including one-carbon metabolism (OCM) and the transsulfuration pathway with possible consequences on glutathione (GSH) levels. In this study, complementary RNAseq and metabolomics data were integrated to examine the hepatic transsulfuration pathway and glutathione biosynthesis in mice following treatment with TCDD every 4 days for 28 days. TCDD dose-dependently repressed hepatic cystathionine ß-synthase (CBS) and cystathionine γ-lyase (CTH) mRNA and protein levels. Reduced CBS and CTH levels are also correlated with dose-dependent decreases in hepatic extract hydrogen sulfide (H2S). In contrast, cysteine levels increased consistent with the induction of Slc7a11, which encodes for the cystine/glutamate Xc- antiporter. Cotreatment of primary hepatocytes with sulfasalazine, a cystine/glutamate Xc- antiporter inhibitor, decreased labeled cysteine incorporation into GSH with a corresponding increase in TCDD cytotoxicity. Although reduced and oxidized GSH levels were unchanged following treatment due to the induction of GSH/GSSG efflux transporter by TCDD, the GSH:GSSG ratio decreased and global protein S-glutathionylation levels in liver extracts increased in response to oxidative stress along with the induction of glutamate-cysteine ligase catalytic subunit (Gclc), glutathione synthetase (Gss), glutathione disulfide reductase (Gsr), and glutathione transferase π (Gstp). Furthermore, levels of ophthalmic acid, a biomarker of oxidative stress indicating GSH consumption, were also increased. Collectively, the data suggest that increased cystine transport due to cystine/glutamate Xc- antiporter induction compensated for decreased cysteine production following repression of the transsulfuration pathway to support GSH synthesis in response to TCDD-induced oxidative stress.

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.

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.

Thioesterase induction by 2,3,7,8-tetrachlorodibenzo-p-dioxin results in a futile cycle that inhibits hepatic ß-oxidation.

2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD), a persistent environmental contaminant, induces steatosis by increasing hepatic uptake of dietary and mobilized peripheral fats, inhibiting lipoprotein export, and repressing ß-oxidation. In this study, the mechanism of ß-oxidation inhibition was investigated by testing the hypothesis that TCDD dose-dependently repressed straight-chain fatty acid oxidation gene expression in mice following oral gavage every 4 days for 28 days. Untargeted metabolomic analysis revealed a dose-dependent decrease in hepatic acyl-CoA levels, while octenoyl-CoA and dicarboxylic acid levels increased. TCDD also dose-dependently repressed the hepatic gene expression associated with triacylglycerol and cholesterol ester hydrolysis, fatty acid binding proteins, fatty acid activation, and 3-ketoacyl-CoA thiolysis while inducing acyl-CoA hydrolysis. Moreover, octenoyl-CoA blocked the hydration of crotonyl-CoA suggesting short chain enoyl-CoA hydratase (ECHS1) activity was inhibited. Collectively, the integration of metabolomics and RNA-seq data suggested TCDD induced a futile cycle of fatty acid activation and acyl-CoA hydrolysis resulting in incomplete ß-oxidation, and the accumulation octenoyl-CoA levels that inhibited the activity of short chain enoyl-CoA hydratase (ECHS1).

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.