Multi-omics profiling of cholangiocytes reveals sex-specific chromatin state dynamics during hepatic cystogenesis in polycystic liver disease

HDI-STARR-seq: Condition-specific enhancer discovery in mouse liver in vivo

Background

STARR-seq and other massively-parallel reporter assays are widely used to discover functional enhancers in transfected cell models, which can be confounded by plasmid vector-induced type-I interferon immune responses and lack the multicellular environment and endogenous chromatin state of complex mammalian tissues.

Results

Here, we describe HDI-STARR-seq, which combines STARR-seq plasmid library delivery to the liver, by hydrodynamic tail vein injection (HDI), with reporter RNA transcriptional initiation driven by a minimal Albumin promoter, which we show is essential for mouse liver STARR-seq enhancer activity assayed 7 days after HDI. Importantly, little or no vector-induced innate type-I interferon responses were observed. Comparisons of HDI-STARR-seq activity between male and female mouse livers and in livers from males treated with an activating ligand of the transcription factor CAR (Nr1i3) identified many condition-dependent enhancers linked to condition-specific gene expression. Further, thousands of active liver enhancers were identified using a high complexity STARR-seq library comprised of ~ 50,000 genomic regions released by DNase-I digestion of mouse liver nuclei. When compared to stringently inactive library sequences, the active enhancer sequences identified were highly enriched for liver open chromatin regions with activating histone marks (H3K27ac, H3K4me1, H3K4me3), were significantly closer to gene transcriptional start sites, and were significantly depleted of repressive (H3K27me3, H3K9me3) and transcribed region histone marks (H3K36me3).

Conclusions

HDI-STARR-seq offers substantial improvements over current methodologies for large scale, functional profiling of enhancers, including condition-dependent enhancers, in liver tissue in vivo, and can be adapted to characterize enhancer activities in a variety of species and tissues by selecting suitable tissue- and species-specific promoter sequences.

HDI-STARR-seq: Condition-specific enhancer discovery in mouse liver in vivo

STARR-seq and other massively-parallel reporter assays are widely used to discover functional enhancers in transfected cell models, which can be confounded by plasmid vector-induced type-I interferon immune responses and lack the multicellular environment and endogenous chromatin state of complex mammalian tissues. Here, we describe HDI-STARR-seq, which combines STARR-seq plasmid library delivery to the liver, by hydrodynamic tail vein injection (HDI), with reporter RNA transcriptional initiation driven by a minimal Albumin promoter, which we show is essential for mouse liver STARR-seq enhancer activity assayed 7 days after HDI. Importantly, little or no vector-induced innate type-I interferon responses were observed. Comparisons of HDI-STARR-seq activity between male and female mouse livers and in livers from males treated with an activating ligand of the transcription factor CAR (Nr1i3) identified many condition-dependent enhancers linked to condition-specific gene expression. Further, thousands of active liver enhancers were identified using a high complexity STARR-seq library comprised of ~50,000 genomic regions released by DNase-I digestion of mouse liver nuclei. When compared to stringently inactive library sequences, the active enhancer sequences identified were highly enriched for liver open chromatin regions with activating histone marks (H3K27ac, H3K4me1, H3K4me3), were significantly closer to gene transcriptional start sites, and were significantly depleted of repressive (H3K27me3, H3K9me3) and transcribed region histone marks (H3K36me3). HDI-STARR-seq offers substantial improvements over current methodologies for large scale, functional profiling of enhancers, including condition-dependent enhancers, in liver tissue in vivo, and can be adapted to characterize enhancer activities in a variety of species and tissues by selecting suitable tissue- and species-specific promoter sequences.

Dissecting microenvironment in cystadenomas and hepatic cysts based on single nucleus RNA-sequencing data

Hepatic cystadenoma is a rare disease, accounting for about 5% of all cystic lesions, with a high tendency of malignant transformation. The preoperative diagnosis of cystadenoma is difficult, and some cystadenomas are easily misdiagnosed as hepatic cysts at first. Hepatic cyst is a relatively common liver disease, most of which are benign, but large hepatic cysts can lead to pressure on the bile duct, resulting in abnormal liver function. To better understand the difference between the microenvironment of cystadenomas and hepatic cysts, we performed single-nuclei RNA-sequencing on cystadenoma and hepatic cysts samples. In addition, we performed spatial transcriptome sequencing of hepatic cysts. Based on nucleus RNA-sequencing data, a total of seven major cell types were identified. Here we described the tumor microenvironment of cystadenomas and hepatic cysts, particularly the transcriptome signatures and regulators of immune cells and stromal cells. By inferring copy number variation, it was found that the malignant degree of hepatic stellate cells in cystadenoma was higher. Pseudotime trajectory analysis demonstrated dynamic transformation of hepatocytes in hepatic cysts and cystadenomas. Cystadenomas had higher immune infiltration than hepatic cysts, and T cells had a more complex regulatory mechanism in cystadenomas than hepatic cysts. Immunohistochemistry confirms a cystadenoma-specific T-cell immunoregulatory mechanism. These results provided a single-cell atlas of cystadenomas and hepatic cyst, revealed a more complex microenvironment in cystadenomas than in hepatic cysts, and provided new perspective for the molecular mechanisms of cystadenomas and hepatic cyst.

Multi-omics approach characterizes the role of Bisphenol F in disrupting hepatic lipid metabolism

Bisphenol F (BPF), a substitute for bisphenol A (BPA), is ubiquitous existed in various environmental media. Exposure to BPF may promote non-alcoholic fatty liver disease (NAFLD), while the potential mechanism is still unknown. In current study, we used in vitro and in vivo model to evaluate its hepatotoxicity and molecular mechanism. Using multi-omics approach, we found that BPF exposure led to changes in hepatic transcriptome, metabolome and chromatin accessible regions that were enriched for binding sites of transcription factors in bZIP family. These alterations were enriched with pathways integral to the endoplasmic reticulum stress and NAFLD. These findings suggested that BPF exposure might reprogram the chromatin accessibility and enhancer landscape in the liver, with downstream effects on genes associated with endoplasmic reticulum stress and lipid metabolism, which relied on bZIP family transcription factors. Overall, our study describes comprehensive molecular alterations in hepatocytes after BPF exposure and provides new insights into the understanding of the hepatoxicity of BPF.

Histone 3 lysine 9 acetylation-specific reprogramming regulates esophageal squamous cell carcinoma progression and metastasis

Dysregulation of histone acetylation is widely implicated in tumorigenesis, yet its specific roles in the progression and metastasis of esophageal squamous cell carcinoma (ESCC) remain unclear. Here, we profiled the genome-wide landscapes of H3K9ac for paired adjacent normal (Nor), primary ESCC (EC) and metastatic lymph node (LNC) esophageal tissues from three ESCC patients. Compared to H3K27ac, we identified a distinct epigenetic reprogramming specific to H3K9ac in EC and LNC samples relative to Nor samples. This H3K9ac-related reprogramming contributed to the transcriptomic aberration of targeting genes, which were functionally associated with tumorigenesis and metastasis. Notably, genes with gained H3K9ac signals in both primary and metastatic lymph node samples (common-gained gene) were significantly enriched in oncogenes. Single-cell RNA-seq analysis further revealed that the corresponding top 15 common-gained genes preferred to be enriched in mesenchymal cells with high metastatic potential. Additionally, in vitro experiment demonstrated that the removal of H3K9ac from the common-gained gene MSI1 significantly downregulated its transcription, resulting in deficiencies in ESCC cell proliferation and migration. Together, our findings revealed the distinct characteristics of H3K9ac in esophageal squamous cell carcinogenesis and metastasis, and highlighted the potential therapeutic avenue for intervening ESCC through epigenetic modulation via H3K9ac.