Nuclear RNA catabolism controls endogenous retroviruses, gene expression asymmetry, and dedifferentiation.

Endogenous retroviruses (ERVs) are remnants of ancient parasitic infections and comprise sizable portions of most genomes. Although epigenetic mechanisms silence most ERVs by generating a repressive environment that prevents their expression (heterochromatin), little is known about mechanisms silencing ERVs residing in open regions of the genome (euchromatin). This is particularly important during embryonic development, where induction and repression of distinct classes of ERVs occur in short temporal windows. Here, we demonstrate that transcription-associated RNA degradation by the nuclear RNA exosome and Integrator is a regulatory mechanism that controls the productive transcription of most genes and many ERVs involved in preimplantation development. Disrupting nuclear RNA catabolism promotes dedifferentiation to a totipotent-like state characterized by defects in RNAPII elongation and decreased expression of long genes (gene-length asymmetry). Our results indicate that RNA catabolism is a core regulatory module of gene networks that safeguards RNAPII activity, ERV expression, cell identity, and developmental potency.

The transcription factor NFIL3/E4BP4 regulates the developmental stage-specific acquisition of basophil function.

BACKGROUND: Basophils are rare but important effector cells in many allergic disorders. Contrary to their early progenitors, the terminal developmental processes of basophils in which they gain their unique functional properties are unknown. OBJECTIVE: We sought to identify a novel late-stage basophil precursor and a transcription factor regulating the terminal maturation of basophils. METHODS: Using flow cytometry, transcriptome analysis, and functional assays, we investigated the identification and functionality of the basophil precursors as well as basophil development. We generated mice with basophil-specific deletion of nuclear factor IL-3 (NFIL3)/E4BP4 and analyzed the functional impairment of NFIL3/E4BP4-deficient basophils in vitro and in vivo using an oxazolone-induced murine model of allergic dermatitis. RESULTS: We report a new mitotic transitional basophil precursor population (referred to as transitional basophils) that expresses the FcεRIα chain at higher levels than mature basophils. Transitional basophils are less responsive to IgE-linked degranulation but produce more cytokines in response to IL-3, IL-33, or IgE cross-linking than mature basophils. In particular, we found that the expression of NFIL3/E4BP4 gradually rises as cells mature from the basophil progenitor stage. Basophil-specific deletion of NFIL3/E4BP4 reduces the expression of genes necessary for basophil function and impairs IgE receptor signaling, cytokine secretion, and degranulation in the context of murine atopic dermatitis. CONCLUSIONS: We discovered transitional basophils, a novel late-stage mitotic basophil precursor cell population that exists between basophil progenitors and postmitotic mature basophils. We demonstrated that NFIL3/E4BP4 augments the IgE-mediated functions of basophils, pointing to a potential therapeutic regulator for allergic diseases.

Impact of media compositions and culture systems on the immunophenotypes of patient-derived breast cancer cells.

BACKGROUND: Heterogeneous tumor cells are thought to be a significant factor in the failure of endocrine therapy in estrogen receptor-positive (ER+) cancers. Culturing patient-derived breast cancer cells (PDBCCs) provides an invaluable tool in pre-clinical and translational research for the heterogeneity of cancer cells. This study aimed to investigate the effects of different media components and culture methods on the BCSC-associated immunophenotypes and gene expression in ER + PDBCCs. METHODS: Ten patients with ER + breast cancer were employed in this study, six of whom had neoadjuvant chemotherapy and four of whom did not. PDBCCs were isolated by enzymatic methods using collagen I and hyaluronidase. PDBCCs were grown as monolayers in mediums with different compositions and as multicellular spheroid in a suspended condition. Collagen I-coated plate and ultralow attachment plate coated with polymer-X were used for monolayer and spheroid culture. Flow cytometry, immunofluorescent staining, RT-PCR, and RNA-sequencing were employed to examine the immunophenotype and genetic profile of PDBCCs. RESULTS: More than 95% of PDBCCs sustain EpCAM high/+/fibroblast marker- phenotypes in monolayer conditions by subculturing 3-4 times. A83-01 removal induced senescent cells with high ß-galactosidase activity. PDBCCs grown as monolayers were characterized by the majority of cells having an EpCAM+/CD49f + phenotype. Compared to full media in monolayer culture, EGF removal increased EpCAM+/CD49f - phenotype (13.8-fold, p = 0.028), whereas R-spondin removal reduced it (0.8-fold, p = 0.02). A83-01 removal increased EpCAM+/CD24 + phenotype (1.82-fold, p = 0.023) and decreased EpCAM low/-/CD44+/CD24- phenotype (0.45-fold, p = 0.026). Compared to monolayer, spheroid resulted in a significant increase in the population with EpCAM-/CD49+ (14.6-fold, p = 0.006) and EpCAM low/-/CD44+/CD24- phenotypes (4.16-fold, p = 0.022) and ALDH high activity (9.66-fold, p = 0.037). ALDH1A and EMT-related genes were upregulated. In RNA-sequencing analysis between spheroids and monolayers, a total of 561 differentially expressed genes (2-fold change, p < 0.05) were enriched in 27 KEGG pathways including signaling pathways regulating pluripotency of stem cells. In a recurrence-free survival analysis based on the Kaplan-Meier Plotter database of the up-and down-regulated genes identified in spheroids, 15 up-, and 14 down-regulated genes were associated with poor prognosis of breast cancer patients. CONCLUSION: The media composition and spheroid culture method change in the BCSCs and EMT markers of PDBCCs, implying the importance of defining the media composition and culture method for studying PDBCCs in vitro.

Octopus-toolkit: a workflow to automate mining of public epigenomic and transcriptomic next-generation sequencing data.

Octopus-toolkit is a stand-alone application for retrieving and processing large sets of next-generation sequencing (NGS) data with a single step. Octopus-toolkit is an automated set-up-and-analysis pipeline utilizing the Aspera, SRA Toolkit, FastQC, Trimmomatic, HISAT2, STAR, Samtools, and HOMER applications. All the applications are installed on the user's computer when the program starts. Upon the installation, it can automatically retrieve original files of various epigenomic and transcriptomic data sets, including ChIP-seq, ATAC-seq, DNase-seq, MeDIP-seq, MNase-seq and RNA-seq, from the gene expression omnibus data repository. The downloaded files can then be sequentially processed to generate BAM and BigWig files, which are used for advanced analyses and visualization. Currently, it can process NGS data from popular model genomes such as, human (Homo sapiens), mouse (Mus musculus), dog (Canis lupus familiaris), plant (Arabidopsis thaliana), zebrafish (Danio rerio), fruit fly (Drosophila melanogaster), worm (Caenorhabditis elegans), and budding yeast (Saccharomyces cerevisiae) genomes. With the processed files from Octopus-toolkit, the meta-analysis of various data sets, motif searches for DNA-binding proteins, and the identification of differentially expressed genes and/or protein-binding sites can be easily conducted with few commands by users. Overall, Octopus-toolkit facilitates the systematic and integrative analysis of available epigenomic and transcriptomic NGS big data.

Nuclear aconitase antagonizes heterochromatic silencing by interfering with Chp1 binding to DNA.

In Schizosaccharomyces pombe, there are two aconitases, Aco1 and Aco2, involved in the Krebs cycle in mitochondria. Interestingly, Aco2 is localized to nucleus as well. Here, we investigated the nuclear role of Aco2 by deleting its nuclear localization signal. The aco2ΔNLS mutation suppressed the gene-silencing defects of RNAi mutants at the centromere, where heterochromatin formation depends on RNAi pathway. In Δago1, the aco2ΔNLS mutation restored heterochromatin through elevating Chp1 binding. Aco2 physically interacted with Chp1 via the N-terminal chromodomain that binds to methylated histone H3K9. In the sub-telomeric region, where heterochromatin forms independent of RNAi pathway, the single aco2ΔNLS mutation caused extra gene silencing via elevating Chp1 binding, without increasing histone methylation. The anti-silencing effect did not require the catalytic function of aconitase. Taken together, Aco2 functions as an epigenetic regulator of gene expression, through associating with chromodomain of Chp1 to maintain heterochromatin.

The 19S proteasome is directly involved in the regulation of heterochromatin spreading in fission yeast.

Cumulative evidence suggests that non-proteolytic functions of the proteasome are involved in transcriptional regulation, mRNA export, and ubiquitin-dependent histone modification and thereby modulate the intracellular levels of regulatory proteins implicated in controlling key cellular functions. To date, the non-proteolytic roles of the proteasome have been mainly investigated in euchromatin; their effects on heterochromatin are largely unknown. Here, using fission yeast as a model, we randomly mutagenized the subunits of the 19S proteasome subcomplex and sought to uncover a direct role of the proteasome in heterochromatin regulation. We identified a mutant allele, rpt4-1, that disrupts a non-proteolytic function of the proteasome, also known as a non-proteolytic allele. Experiments performed using rpt4-1 cells revealed that the proteasome is involved in the regulation of heterochromatin spreading to prevent its uncontrolled invasion into neighboring euchromatin regions. Intriguingly, the phenotype of the non-proteolytic rpt4-1 mutant resembled that of epe1Δ cells, which lack the Epe1 protein that counteracts heterochromatin spreading. Both mutants exhibited variegated gene-silencing phenotypes across yeast colonies, spreading of heterochromatin, bypassing of the requirement for RNAi in heterochromatin formation at the outer repeat region (otr), and up-regulation of RNA polymerase II. Further analysis revealed Mst2, another factor that antagonizes heterochromatin spreading, may function redundantly with Rpt4. These observations suggest that the 19S proteasome may be involved in modulating the activities of Epe1 and Mst2. In conclusion, our findings indicate that the proteasome appears to have a heterochromatin-regulating function that is independent of its canonical function in proteolysis.

Chd1p recognizes H3K36Ac to maintain nucleosome positioning near the transcription start site.

In Saccharomyces cerevisiae, the ATP-dependent chromatin remodeler, Chd1p, globally affects nucleosome positioning at coding regions, where nucleosomes are specifically and directionally aligned with respect to the transcription start site (TSS). Various auxiliary domains of remodelers play critical roles by performing specialized functions that are unique to the type of remodeler. Here, we report that yeast Chd1p directly binds to acetylated histone H3K36 (H3K36Ac) via its chromodomain, and that H3K36Ac stimulates the nucleosome sliding activity of Chd1p in vitro. Furthermore, we use genome-wide analysis to demonstrate that H3K36Ac promotes the remodeling activity of Chd1p to maintain chromatin stability at the 5' ends of genes in vivo. Our work linking Chd1p with H3K36Ac provides novel insights into how the nucleosome remodeling activity of Chd1p is controlled near the TSS.

The 19S proteasome regulates subtelomere silencing and facultative heterochromatin formation in fission yeast.

Accumulating evidence shows that non-proteolytic functions of the proteasome are as crucial as its well-known proteolytic function in regulating cellular activities. In our recent work, we showed that the 19S proteasome mediates the heterochromatin spreading of centromeric heterochromatin in non-proteolytic manner. However, the involvement of the proteasome in other heterochromatin regions remained largely unknown. In the present study, we investigated the non-proteolytic role of the 19S proteasome in subtelomere and facultative heterochromatin regions. Using the non-proteolytic mutant, rpt4-1, we show that the 19S proteasome is involved in regulating subtelomere silencing and facultative heterochromatin formation in fission yeast. In addition to this proteasome-related regulation, we also observed a distinct pathway that regulates subtelomere silencing and facultative heterochromatin formation through the Paf1 complex subunit, Leo1. Our comparison of the two pathways revealed a new group of heterochromatin domains that are regulated exclusively by the proteasome pathway. Taken together, our findings reveal that the proteasome is involved in the global regulation of facultative and constitutive heterochromatin.

Dynamic changes in DNA methylation and hydroxymethylation when hES cells undergo differentiation toward a neuronal lineage.

DNA methylation and hydroxymethylation have been implicated in normal development and differentiation, but our knowledge is limited about the genome-wide distribution of 5-methylcytosine (5 mC) and 5-hydroxymethylcytosine (5 hmC) during cellular differentiation. Using an in vitro model system of gradual differentiation of human embryonic stem (hES) cells into ventral midbrain-type neural precursor cells and terminally into dopamine neurons, we observed dramatic genome-wide changes in 5 mC and 5 hmC patterns during lineage commitment. The 5 hmC pattern was dynamic in promoters, exons and enhancers. DNA hydroxymethylation within the gene body was associated with gene activation. The neurogenesis-related genes NOTCH1, RGMA and AKT1 acquired 5 hmC in the gene body and were up-regulated during differentiation. DNA methylation in the promoter was associated with gene repression. The pluripotency-related genes POU5F1, ZFP42 and HMGA1 acquired 5 mC in their promoters and were down-regulated during differentiation. Promoter methylation also acted as a locking mechanism to maintain gene silencing. The mesoderm development-related genes NKX2-8, TNFSF11 and NFATC1 acquired promoter methylation during neural differentiation even though they were already silenced in hES cells. Our findings will help elucidate the molecular mechanisms underlying lineage-specific differentiation of pluripotent stem cells during human embryonic development.

Hrp3 controls nucleosome positioning to suppress non-coding transcription in eu- and heterochromatin.

The positioning of the nucleosome by ATP-dependent remodellers provides the fundamental chromatin environment for the regulation of diverse cellular processes acting on the underlying DNA. Recently, genome-wide nucleosome mapping has revealed more detailed information on the chromatin-remodelling factors. Here, we report that the Schizosaccharomyces pombe CHD remodeller, Hrp3, is a global regulator that drives proper nucleosome positioning and nucleosome stability. The loss of Hrp3 resulted in nucleosome perturbation across the chromosome, and the production of antisense transcripts in the hrp3Δ cells emphasized the importance of nucleosome architecture for proper transcription. Notably, perturbation of the nucleosome in hrp3 deletion mutant was also associated with destabilization of the DNA-histone interaction and cell cycle-dependent alleviation of heterochromatin silencing. Furthermore, the effect of Hrp3 in the pericentric region was found to be accomplished via a physical interaction with Swi6, and appeared to cooperate with other heterochromatin factors for gene silencing. Taken together, our data indicate that a well-positioned nucleosome by Hrp3 is important for the spatial-temporal control of transcription-associated processes.

Both H4K20 mono-methylation and H3K56 acetylation mark transcription-dependent histone turnover in fission yeast.

Nucleosome dynamics facilitated by histone turnover is required for transcription as well as DNA replication and repair. Histone turnover is often associated with various histone modifications such as H3K56 acetylation (H3K56Ac), H3K36 methylation (H3K36me), and H4K20 methylation (H4K20me). In order to correlate histone modifications and transcription-dependent histone turnover, we performed genome wide analyses for euchromatic regions in G2/M-arrested fission yeast. The results show that transcription-dependent histone turnover at 5' promoter and 3' termination regions is directly correlated with the occurrence of H3K56Ac and H4K20 mono-methylation (H4K20me1) in actively transcribed genes. Furthermore, the increase of H3K56Ac and H4K20me1 and antisense RNA production was observed in the absence of the histone H3K36 methyltransferase Set2 and histone deacetylase complex (HDAC) that are involved in the suppression of histone turnover within the coding regions. These results together indicate that H4K20me1 as well as H3K56Ac are bona fide marks for transcription-dependent histone turnover in fission yeast.

A histone methyltransferase inhibits seed germination by increasing PIF1 mRNA expression in imbibed seeds.

Phytochrome-interacting factor 1 (PIF1) inhibits light-dependent seed germination. The specific function of PIF1 in seed germination is partly due to its high level of expression in imbibed seeds, but the associated regulatory factors have not been identified. Here we show that mutation of the early flowering in short days (EFS) gene, encoding an H3K4 and H3K36 methyltransferase, decreases the level of H3K36me2 and H3K36me3 but not H3K4me3 at the PIF1 locus, reduces the targeting of RNA polymerase II to the PIF1 locus, and reduces mRNA expression of PIF1 in imbibed seeds. Consistently, the efs mutant geminated even under the phyBoff condition, and had an expression profile of PIF1 target genes similar to that of the pif1 mutant. Introduction of an EFS transgene into the efs mutant restored the level of H3K36me2 and H3K36me3 at the PIF1 locus, the high-level expression of PIF1 mRNA, the expression pattern of PIF1 target genes, and the light-dependent germination of these seeds. Introduction of a PIF1 transgene into the efs mutant also restored the expression pattern of PIF1 target genes and light-dependent germination in imbibed seeds, but did not restore the flowering phenotype. Taken together, our results indicate that EFS is necessary for high-level expression of PIF1 mRNA in imbibed seeds.

H2B monoubiquitylation is a 5'-enriched active transcription mark and correlates with exon-intron structure in human cells.

H2B monoubiquitylation (H2Bub1), which is required for multiple methylations of both H3K4 and H3K79, has been implicated in gene expression in numerous organisms ranging from yeast to human. However, the molecular crosstalk between H2Bub1 and other modifications, especially the methylations of H3K4 and H3K79, remains unclear in vertebrates. To better understand the functional role of H2Bub1, we measured genome-wide histone modification patterns in human cells. Our results suggest that H2Bub1 has dual roles, one that is H3 methylation dependent, and another that is H3 methylation independent. First, H2Bub1 is a 5'-enriched active transcription mark and co-occupies with H3K79 methylations in actively transcribed regions. Second, this study shows for the first time that H2Bub1 plays a histone H3 methylations-independent role in chromatin architecture. Furthermore, the results of this work indicate that H2Bub1 is largely positioned at the exon-intron boundaries of highly expressed exons, and it demonstrates increased occupancy in skipped exons compared with flanking exons in the human and mouse genomes. Our findings collectively suggest that a potentiating mechanism links H2Bub1 to both H3K79 methylations in actively transcribed regions and the exon-intron structure of highly expressed exons via the regulation of nucleosome dynamics during transcription elongation.

Genome-wide analyses reveal the extent of opportunistic STAT5 binding that does not yield transcriptional activation of neighboring genes.

Signal Transducers and Activators of Transcription (STAT) 5A/B regulate cytokine-inducible genes upon binding to GAS motifs. It is not known what percentage of genes with GAS motifs bind to and are regulated by STAT5. Moreover, it is not clear whether genome-wide STAT5 binding is modulated by its concentration. To clarify these issues we established genome-wide STAT5 binding upon growth hormone (GH) stimulation of wild-type (WT) mouse embryonic fibroblasts (MEFs) and MEFs overexpressing STAT5A more than 20-fold. Upon GH stimulation, 23 827 and 111 939 STAT5A binding sites were detected in WT and STAT5A overexpressing MEFs, respectively. 13 278 and 71 561 peaks contained at least one GAS motif. 1586 and 8613 binding sites were located within 2.5 kb of promoter sequences, respectively. Stringent filtering revealed 78 genes in which the promoter/upstream region (-10 kb to +0.5 kb) was recognized by STAT5 both in WT and STAT5 overexpressing MEFs and 347 genes that bound STAT5 only in overexpressing cells. Genome-wide expression analyses identified that the majority of STAT5-bound genes was not under GH control. Up to 40% of STAT5-bound genes were not expressed. For the first time we demonstrate the magnitude of opportunistic genomic STAT5 binding that does not translate into transcriptional activation of neighboring genes.

Human histone H3K79 methyltransferase DOT1L protein [corrected] binds actively transcribing RNA polymerase II to regulate gene expression.

Histone-modifying enzymes play a pivotal role in gene expression and repression. In human, DOT1L (Dot1-like) is the only known histone H3 lysine 79 methyltransferase. hDOT1L is associated with transcriptional activation, but the general mechanism connecting hDOT1L to active transcription remains largely unknown. Here, we report that hDOT1L interacts with the phosphorylated C-terminal domain of actively transcribing RNA polymerase II (RNAPII) through a region conserved uniquely in multicellular DOT1 proteins. Genome-wide profiling analyses indicate that the occupancy of hDOT1L largely overlaps with that of RNAPII at actively transcribed genes, especially surrounding transcriptional start sites, in embryonic carcinoma NCCIT cells. We also find that C-terminal domain binding or H3K79 methylations by hDOT1L is important for the expression of target genes such as NANOG and OCT4 and a marker for pluripotency in NCCIT cells. Our results indicate that a functional interaction between hDOT1L and RNAPII targets hDOT1L and subsequent H3K79 methylations to actively transcribed genes.

High-resolution crystal structure of the catalytic domain of human dual-specificity phosphatase 26.

Dual-specificity phosphatases (DUSPs) play an important role in regulating cellular signalling pathways governing cell growth, differentiation and apoptosis. Human DUSP26 inhibits the apoptosis of cancer cells by dephosphorylating substrates such as p38 and p53. High-resolution crystal structures of the DUSP26 catalytic domain (DUSP26-C) and its C152S mutant [DUSP26-C (C152S)] have been determined at 1.67 and 2.20 Å resolution, respectively. The structure of DUSP26-C showed a novel type of domain-swapped dimer formed by extensive crossover of the C-terminal α7 helix. Taken together with the results of a phosphatase-activity assay, structural comparison with other DUSPs revealed that DUSP26-C adopts a catalytically inactive conformation of the protein tyrosine phosphate-binding loop which significantly deviates from that of canonical DUSP structures. In particular, a noticeable difference exists between DUSP26-C and the active forms of other DUSPs at the hinge region of a swapped C-terminal domain. Additionally, two significant gaps were identified between the catalytic core and its surrounding loops in DUSP26-C, which can be exploited as additional binding sites for allosteric enzyme regulation. The high-resolution structure of DUSP26-C may thus provide structural insights into the rational design of DUSP26-targeted anticancer drugs.

Identification of DNA methylation markers for lineage commitment of in vitro hepatogenesis.

Hepatocytes that have differentiated from human embryonic stem cells (hESCs) have great potential for the treatment of liver disease as well as for drug testing. Moreover, in vitro hepatogenesis is a powerful model system for studying the molecular mechanisms underlying liver development. DNA methylation is an important epigenetic mechanism that influences differential gene expression during embryonic development. We profiled gene expression and DNA methylation of three cell states of in vitro hepatogenesis-hESC, definitive endoderm and hepatocyte-using microarray analysis. Among 525 state-specific expressed genes, 67 showed significant negative correlation between gene expression and DNA methylation. State-specific expression and methylation of target genes were validated by quantitative reverse transcription-polymerase chain reaction and pyrosequencing, respectively. To elucidate genome-scale methylation changes beyond the promoter, we also performed high-throughput sequencing of methylated DNA captured by the MBD2 protein. We found dynamic methylation changes in intergenic regions of the human genome during differentiation. This study provides valuable methylation markers for the lineage commitment of in vitro hepatogenesis and should help elucidate the molecular mechanisms underlying stem cell differentiation and liver development.

Decoding the genome with an integrative analysis tool: combinatorial CRM Decoder.

The identification of genome-wide cis-regulatory modules (CRMs) and characterization of their associated epigenetic features are fundamental steps toward the understanding of gene regulatory networks. Although integrative analysis of available genome-wide information can provide new biological insights, the lack of novel methodologies has become a major bottleneck. Here, we present a comprehensive analysis tool called combinatorial CRM decoder (CCD), which utilizes the publicly available information to identify and characterize genome-wide CRMs in a species of interest. CCD first defines a set of the epigenetic features which is significantly associated with a set of known CRMs as a code called 'trace code', and subsequently uses the trace code to pinpoint putative CRMs throughout the genome. Using 61 genome-wide data sets obtained from 17 independent mouse studies, CCD successfully catalogued ∼12 600 CRMs (five distinct classes) including polycomb repressive complex 2 target sites as well as imprinting control regions. Interestingly, we discovered that ∼4% of the identified CRMs belong to at least two different classes named 'multi-functional CRM', suggesting their functional importance for regulating spatiotemporal gene expression. From these examples, we show that CCD can be applied to any potential genome-wide datasets and therefore will shed light on unveiling genome-wide CRMs in various species.

The versatility of the proteasome in gene expression and silencing: Unraveling proteolytic and non-proteolytic functions.

The 26S proteasome consists of a 20S core particle and a 19S regulatory particle and critically regulates gene expression and silencing through both proteolytic and non-proteolytic functions. The 20S core particle mediates proteolysis, while the 19S regulatory particle performs non-proteolytic functions. The proteasome plays a role in regulating gene expression in euchromatin by modifying histones, activating transcription, initiating and terminating transcription, mRNA export, and maintaining transcriptome integrity. In gene silencing, the proteasome modulates the heterochromatin formation, spreading, and subtelomere silencing by degrading specific proteins and interacting with anti-silencing factors such as Epe1, Mst2, and Leo1. This review discusses the proteolytic and non-proteolytic functions of the proteasome in regulating gene expression and gene silencing-related heterochromatin formation. This article is part of a special issue on the regulation of gene expression and genome integrity by the ubiquitin-proteasome system.

Cell line-specific features of 3D chromatin organization in hepatocellular carcinoma.

Liver cancer, particularly hepatocellular carcinoma (HCC), poses a significant global threat to human lives. To advance the development of innovative diagnostic and treatment approaches, it is essential to examine the hidden features of HCC, particularly its 3D genome architecture, which is not well understood. In this study, we investigated the 3D genome organization of four HCC cell lines-Hep3B, Huh1, Huh7, and SNU449-using in situ Hi-C and assay for transposase-accessible chromatin sequencing. Our findings revealed that HCC cell lines had more long-range interactions, both intra-and interchromosomal, compared to human mammary epithelial cells (HMECs). Unexpectedly, HCC cell lines displayed cell line-specific compartmental modifications at the megabase (Mb) scale, which could potentially be leveraged in determining HCC subtypes. At the sub-Mb scale, we observed decreases in intra-TAD (topologically associated domain) interactions and chromatin loops in HCC cell lines compared to HMECs. Lastly, we discovered a correlation between gene expression and the 3D chromatin architecture of SLC8A1, which encodes a sodium-calcium antiporter whose modulation is known to induce apoptosis by comparison between HCC cell lines and HMECs. Our findings suggest that HCC cell lines have a distinct 3D genome organization that is different from those of normal and other cancer cells based on the analysis of compartments, TADs, and chromatin loops. Overall, we take this as evidence that genome organization plays a crucial role in cancer phenotype determination. Further exploration of epigenetics in HCC will help us to better understand specific gene regulation mechanisms and uncover novel targets for cancer treatment.