No abstract

Cell-specific transcriptional responses are lost in the averages of bulk RNA sequencing. We performed single nuclei RNA sequencing (snSeq) on frozen liver samples from male C57BL/6 mice in response to 2,3,7,8tetrachlorodibenzo-p-dioxin (TCDD). Approximately 19,907 hepatic genes were detected across 16,015 sequenced nuclei from control and treated samples. Eleven cell-(sub)types were identified including distinct hepatocyte sub-populations, consistent with the cell diversity of the liver. TCDD increased macrophages from 0.5% to 24.7%, while neutrophils were only present in treated samples. The number of differentially expressed genes correlated with the basal expression level of Ahr. In addition to expected functional enrichments within each cell-(sub)type, RAS signaling was enriched in nonparenchymal cells. snSeq also identified a Kupffer cell subtype highly expressing Gpnmb, consistent with a dietary NASH model. Overall, snSeq distinguished cellspecific transcriptional changes and population shifts consistent with the hepatotoxicity of TCDD.

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 hrs after oral gavage with 30 μg/kg TCDD were integrated with bulk RNA-sequencing (RNAseq) time-course (2 - 72 hrs) and dose-response (0.01 - 30 μg/kg) datasets to assess putative AhR, HNF4α and COUP-TFII interactions associated with differential gene expression. TCDD treatment resulted in the genomic enrichment of 23,701, 11,688, and 9,547 binding regions for AhR, COUP-TFII and HNF4α, respectively, throughout the genome. Functional enrichment analysis of differentially expressed genes (DEGs) identified differential binding enrichment for AhR, COUP-TFII, and HNF4a 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α.