Enhancer infestation drives tumorigenic activation of inactive B compartment in Epstein-Barr virus-positive nasopharyngeal carcinoma.

BACKGROUND: Nasopharyngeal carcinoma (NPC) is an Epstein-Barr virus (EBV)-associated malignant epithelial tumor endemic to Southern China and Southeast Asia. While previous studies have revealed a low frequency of gene mutations in NPC, its epigenomic aberrations are not fully elucidated apart from DNA hypermethylation. Epigenomic rewiring and enhancer dysregulation, such as enhancer hijacking due to genomic structural changes or extrachromosomal DNA, drive cancer progression. METHODS: We conducted Hi-C, 4C-seq, ChIP-seq, and RNA-seq analyses to comprehensively elucidate the epigenome and interactome of NPC using C666-1 EBV(+)-NPC cell lines, NP69T immortalized nasopharyngeal epithelial cells, clinical NPC biopsy samples, and in vitro EBV infection in HK1 and NPC-TW01 EBV(-) cell lines. FINDINGS: In C666-1, the EBV genome significantly interacted with inactive B compartments of host cells; the significant association of EBV-interacting regions (EBVIRs) with B compartment was confirmed using clinical NPC and in vitro EBV infection model. EBVIRs in C666-1 showed significantly higher levels of active histone modifications compared with NP69T. Aberrant activation of EBVIRs after EBV infection was validated using in vitro EBV infection models. Within the EBVIR-overlapping topologically associating domains, 14 H3K4me3(+) genes were significantly upregulated in C666-1. Target genes of EBVIRs including PLA2G4A, PTGS2 and CITED2, interacted with the enhancers activated in EBVIRs and were highly expressed in NPC, and their knockdown significantly reduced cell proliferation. INTERPRETATION: The EBV genome contributes to NPC tumorigenesis through "enhancer infestation" by interacting with the inactive B compartments of the host genome and aberrantly activating enhancers. FUNDING: The funds are listed in the Acknowledgements section.

Core-binding factor fusion downregulation of ADAR2 RNA editing contributes to AML leukemogenesis.

Adenosine-to-inosine RNA editing, which is catalyzed by adenosine deaminases acting on RNA (ADAR) family of enzymes, ADAR1 and ADAR2, has been shown to contribute to multiple cancers. However, other than the chronic myeloid leukemia blast crisis, relatively little is known about its role in other types of hematological malignancies. Here, we found that ADAR2, but not ADAR1 and ADAR3, was specifically downregulated in the core-binding factor (CBF) acute myeloid leukemia (AML) with t(8;21) or inv(16) translocations. In t(8;21) AML, RUNX1-driven transcription of ADAR2 was repressed by the RUNX1-ETO additional exon 9a fusion protein in a dominant-negative manner. Further functional studies confirmed that ADAR2 could suppress leukemogenesis specifically in t(8;21) and inv16 AML cells dependent on its RNA editing capability. Expression of 2 exemplary ADAR2-regulated RNA editing targets coatomer subunit α and component of oligomeric Golgi complex 3 inhibits the clonogenic growth of human t(8;21) AML cells. Our findings support a hitherto, unappreciated mechanism leading to ADAR2 dysregulation in CBF AML and highlight the functional relevance of loss of ADAR2-mediated RNA editing to CBF AML.

Emerging regulatory mechanisms of noncoding RNAs in topologically associating domains.

Topologically associating domains (TADs) are integral to spatial genome organization, instructing gene expression, and cell fate. Recently, several advances have uncovered roles for noncoding RNAs (ncRNAs) in the regulation of the form and function of mammalian TADs. Phase separation has also emerged as a potential arbiter of ncRNAs in the regulation of TADs. In this review we discuss the implications of these novel findings in relation to how ncRNAs might structurally and functionally regulate TADs from two perspectives: moderating loop extrusion through interactions with architectural proteins, and facilitating TAD phase separation. Additionally, we propose future studies and directions to investigate these phenomena.

Identification of mechanism of cancer-cell-specific reactivation of hTERT offers therapeutic opportunities for blocking telomerase specifically in human colorectal cancer.

Transcriptional reactivation of hTERT is the limiting step in tumorigenesis. While mutations in hTERT promoter present in 19% of cancers are recognized as key drivers of hTERT reactivation, mechanisms by which wildtype hTERT (WT-hTERT) promoter is reactivated, in majority of human cancers, remain unknown. Using primary colorectal cancers (CRC) we identified Tert INTeracting region 2 (T-INT2), the critical chromatin region essential for reactivating WT-hTERT promoter in CRCs. Elevated ß-catenin and JunD level in CRC facilitates chromatin interaction between hTERT promoter and T-INT2 that is necessary to turn on hTERTexpression. Pharmacological screens uncovered salinomycin, which inhibits JunD mediated hTERT-T-INT2 interaction that is required for the formation of a stable transcription complex on the hTERT promoter. Our results showed for the first time how known CRC alterations, such as APC, lead to WT-hTERT promoter reactivation during stepwise-tumorigenesis and provide a new perspective for developing cancer-specific drugs.

Healthy and cancer cells harbor the same DNA sequence, but reactivation of the Human Telomerase Reverse Transcriptase (hTERT) gene is observed only in cancer cells. How does that happen was not known for over three decades of research? This study identifies a specific DNA structure that forms only in cancer cells and brings the necessary molecular machinery into the correct position to activate the hTERT gene. The detailed mechanism of hTERT activation provided in this study will be instrumental in designing cancer cell-specific hTERT inhibitors, especially since all the other ways of inhibiting telomerase failed in the clinic.

MCI-frcnn: A deep learning method for topological micro-domain boundary detection.

Chromatin structural domains, or topologically associated domains (TADs), are a general organizing principle in chromatin biology. RNA polymerase II (RNAPII) mediates multiple chromatin interactive loops, tethering together as RNAPII-associated chromatin interaction domains (RAIDs) to offer a framework for gene regulation. RAID and TAD alterations have been found to be associated with diseases. They can be further dissected as micro-domains (micro-TADs and micro-RAIDs) by clustering single-molecule chromatin-interactive complexes from next-generation three-dimensional (3D) genome techniques, such as ChIA-Drop. Currently, there are few tools available for micro-domain boundary identification. In this work, we developed the MCI-frcnn deep learning method to train a Faster Region-based Convolutional Neural Network (Faster R-CNN) for micro-domain boundary detection. At the training phase in MCI-frcnn, 50 images of RAIDs from Drosophila RNAPII ChIA-Drop data, containing 261 micro-RAIDs with ground truth boundaries, were trained for 7 days. Using this well-trained MCI-frcnn, we detected micro-RAID boundaries for the input new images, with a fast speed (5.26 fps), high recognition accuracy (AUROC = 0.85, mAP = 0.69), and high boundary region quantification (genomic IoU = 76%). We further applied MCI-frcnn to detect human micro-TADs boundaries using human GM12878 SPRITE data and obtained a high region quantification score (mean gIoU = 85%). In all, the MCI-frcnn deep learning method which we developed in this work is a general tool for micro-domain boundary detection.

Accurate Identification of DNA Replication Origin by Fusing Epigenomics and Chromatin Interaction Information.

DNA replication initiation is a complex process involving various genetic and epigenomic signatures. The correct identification of replication origins (ORIs) could provide important clues for the study of a variety of diseases caused by replication. Here, we design a computational approach named iORI-Epi to recognize ORIs by incorporating epigenome-based features, sequence-based features, and 3D genome-based features. The iORI-Epi displays excellent robustness and generalization ability on both training datasets and independent datasets of K562 cell line. Further experiments confirm that iORI-Epi is highly scalable in other cell lines (MCF7 and HCT116). We also analyze and clarify the regulatory role of epigenomic marks, DNA motifs, and chromatin interaction in DNA replication initiation of eukaryotic genomes. Finally, we discuss gene enrichment pathways from the perspective of ORIs in different replication timing states and heuristically dissect the effect of promoters on replication initiation. Our computational methodology is worth extending to ORI identification in other eukaryotic species.

MCIBox: a toolkit for single-molecule multi-way chromatin interaction visualization and micro-domains identification.

The emerging ligation-free three-dimensional (3D) genome mapping technologies can identify multiplex chromatin interactions with single-molecule precision. These technologies not only offer new insight into high-dimensional chromatin organization and gene regulation, but also introduce new challenges in data visualization and analysis. To overcome these challenges, we developed MCIBox, a toolkit for multi-way chromatin interaction (MCI) analysis, including a visualization tool and a platform for identifying micro-domains with clustered single-molecule chromatin complexes. MCIBox is based on various clustering algorithms integrated with dimensionality reduction methods that can display multiplex chromatin interactions at single-molecule level, allowing users to explore chromatin extrusion patterns and super-enhancers regulation modes in transcription, and to identify single-molecule chromatin complexes that are clustered into micro-domains. Furthermore, MCIBox incorporates a two-dimensional kernel density estimation algorithm to identify micro-domains boundaries automatically. These micro-domains were stratified with distinctive signatures of transcription activity and contained different cell-cycle-associated genes. Taken together, MCIBox represents an invaluable tool for the study of multiple chromatin interactions and inaugurates a previously unappreciated view of 3D genome structure.

SETDB1 acts as a topological accessory to Cohesin via an H3K9me3-independent, genomic shunt for regulating cell fates.

SETDB1 is a key regulator of lineage-specific genes and endogenous retroviral elements (ERVs) through its deposition of repressive H3K9me3 mark. Apart from its H3K9me3 regulatory role, SETDB1 has seldom been studied in terms of its other potential regulatory roles. To investigate this, a genomic survey of SETDB1 binding in mouse embryonic stem cells across multiple libraries was conducted, leading to the unexpected discovery of regions bereft of common repressive histone marks (H3K9me3, H3K27me3). These regions were enriched with the CTCF motif that is often associated with the topological regulator Cohesin. Further profiling of these non-H3K9me3 regions led to the discovery of a cluster of non-repeat loci that were co-bound by SETDB1 and Cohesin. These regions, which we named DiSCs (domains involving SETDB1 and Cohesin) were seen to be proximal to the gene promoters involved in embryonic stem cell pluripotency and lineage development. Importantly, it was found that SETDB1-Cohesin co-regulate target gene expression and genome topology at these DiSCs. Depletion of SETDB1 led to localized dysregulation of Cohesin binding thereby locally disrupting topological structures. Dysregulated gene expression trends revealed the importance of this cluster in ES cell maintenance as well as at gene 'islands' that drive differentiation to other lineages. The 'unearthing' of the DiSCs thus unravels a unique topological and transcriptional axis of control regulated chiefly by SETDB1.

Super-Enhancers, Phase-Separated Condensates, and 3D Genome Organization in Cancer.

3D chromatin organization plays an important role in transcription regulation and gene expression. The 3D genome is highly maintained by several architectural proteins, such as CTCF, Yin Yang 1, and cohesin complex. This structural organization brings regulatory DNA elements in close proximity to their target promoters. In this review, we discuss the 3D chromatin organization of super-enhancers and their relationship to phase-separated condensates. Super-enhancers are large clusters of DNA elements. They can physically contact with their target promoters by chromatin looping during transcription. Multiple transcription factors can bind to enhancer and promoter sequences and recruit a complex array of transcriptional co-activators and RNA polymerase II to effect transcriptional activation. Phase-separated condensates of transcription factors and transcriptional co-activators have been implicated in assembling the transcription machinery at particular enhancers. Cancer cells can hijack super-enhancers to drive oncogenic transcription to promote cell survival and proliferation. These dysregulated transcriptional programs can cause cancer cells to become highly dependent on transcriptional regulators, such as Mediator and BRD4. Moreover, the expression of oncogenes that are driven by super-enhancers is sensitive to transcriptional perturbation and often occurs in phase-separated condensates, supporting therapeutic rationales of targeting SE components, 3D genome organization, or dysregulated condensates in cancer.

Pan-cancer pervasive upregulation of 3' UTR splicing drives tumourigenesis.

Most mammalian genes generate messenger RNAs with variable untranslated regions (UTRs) that are important post-transcriptional regulators. In cancer, shortening at 3' UTR ends via alternative polyadenylation can activate oncogenes. However, internal 3' UTR splicing remains poorly understood as splicing studies have traditionally focused on protein-coding alterations. Here we systematically map the pan-cancer landscape of 3' UTR splicing and present this in SpUR ( http://www.cbrc.kaust.edu.sa/spur/home/ ). 3' UTR splicing is widespread, upregulated in cancers, correlated with poor prognosis and more prevalent in oncogenes. We show that antisense oligonucleotide-mediated inhibition of 3' UTR splicing efficiently reduces oncogene expression and impedes tumour progression. Notably, CTNNB1 3' UTR splicing is the most consistently dysregulated event across cancers. We validate its upregulation in hepatocellular carcinoma and colon adenocarcinoma, and show that the spliced 3' UTR variant is the predominant contributor to its oncogenic functions. Overall, our study highlights the importance of 3' UTR splicing in cancer and may launch new avenues for RNA-based anti-cancer therapeutics.

Chromatin Rewiring by Mismatch Repair Protein MSH2 Alters Cell Adhesion Pathways and Sensitivity to BET Inhibition in Gastric Cancer.

Mutations in the DNA mismatch repair gene MSH2 are causative of microsatellite instability (MSI) in multiple cancers. Here, we discovered that besides its well-established role in DNA repair, MSH2 exerts a novel epigenomic function in gastric cancer. Unbiased CRISPR-based mass spectrometry combined with genome-wide CRISPR functional screening revealed that in early-stage gastric cancer MSH2 genomic binding is not randomly distributed but rather is associated specifically with tumor-associated super-enhancers controlling the expression of cell adhesion genes. At these loci, MSH2 genomic binding was required for chromatin rewiring, de novo enhancer-promoter interactions, maintenance of histone acetylation levels, and regulation of cell adhesion pathway expression. The chromatin function of MSH2 was independent of its DNA repair catalytic activity but required MSH6, another DNA repair gene, and recruitment to gene loci by the SWI/SNF chromatin remodeler SMARCA4/BRG1. Loss of MSH2 in advanced gastric cancers was accompanied by deficient cell adhesion pathway expression, epithelial-mesenchymal transition, and enhanced tumorigenesis in vitro and in vivo. However, MSH2-deficient gastric cancers also displayed addiction to BAZ1B, a bromodomain-containing family member, and consequent synthetic lethality to bromodomain and extraterminal motif (BET) inhibition. Our results reveal a role for MSH2 in gastric cancer epigenomic regulation and identify BET inhibition as a potential therapy in MSH2-deficient gastric malignancies. SIGNIFICANCE: DNA repair protein MSH2 binds and regulates cell adhesion genes by enabling enhancer-promoter interactions, and loss of MSH2 causes deficient cell adhesion and bromodomain and extraterminal motif inhibitor synthetic lethality in gastric cancer.

Long-Distance Repression by Human Silencers: Chromatin Interactions and Phase Separation in Silencers.

Three-dimensional genome organization represents an additional layer in the epigenetic regulation of gene expression. Active transcription controlled by enhancers or super-enhancers has been extensively studied. Enhancers or super-enhancers can recruit activators or co-activators to activate target gene expression through long-range chromatin interactions. Chromatin interactions and phase separation play important roles in terms of enhancer or super-enhancer functioning. Silencers are another major type of cis-regulatory element that can mediate gene regulation by turning off or reducing gene expression. However, compared to active transcription, silencer studies are still in their infancy. This review covers the current knowledge of human silencers, especially the roles of chromatin interactions and phase separation in silencers. This review also proposes future directions for human silencer studies.

p53-NEIL1 co-abnormalities induce genomic instability and promote synthetic lethality with Chk1 inhibition in multiple myeloma having concomitant 17p13(del) and 1q21(gain).

Recurrent cytogenetic abnormalities are the main hallmark of multiple myeloma (MM) and patients having 2 or more high-risk prognostic events are associated with extremely poor outcome. 17p13(del) and 1q21(gain) are critical and independent high-risk cytogenetic markers, however, the biological significance underlying the poor outcome in MM patients having co-occurrence of both these chromosomal aberrations has never been interrogated. Herein, we identified that patients harbouring concomitant 17p13(del) with 1q21(gain) demonstrated the worst prognosis as compared to patients with single- (either 17p13(del) or 1q21(gain)) and with no chromosomal events (WT for both chromosomal loci); and they are highly enriched for genomic instability (GI) signature. We discovered that the GI feature in the patients with concomitant 17p13(del)-1q21(gain) was recapitulating the biological properties of myeloma cells with co-existing p53-deficiency and NEIL1 mRNA-hyper-editing (associated with chromosome 17p and 1q, respectively) that have inherent DNA damage response (DDR) and persistent activation of Chk1 pathway. Importantly, this became a vulnerable point for therapeutic targeting whereby the cells with this co-abnormalities demonstrated hyper-sensitivity to siRNA- and pharmacological-mediated-Chk1 inhibition, as observed at both the in vitro and in vivo levels. Mechanistically, this was attributable to the synthetic lethal relationship between p53-NEIL1-Chk1 abnormalities. The Chk1 inhibitor (AZD7762) tested showed good synergism with standard-of-care myeloma drugs, velcade and melphalan, thus further reinforcing the translational potential of this therapeutic approach. In summary, combination of NEIL1-p53 abnormalities with an ensuing Chk1 activation could serve as an Achilles heel and predispose MM cells with co-existing 1q21(gain) and 17p13(del) to therapeutic vulnerability for Chk1 inhibition.

MYC overexpression leads to increased chromatin interactions at super-enhancers and MYC binding sites.

The MYC oncogene encodes for the MYC protein and is frequently dysregulated across multiple cancer cell types, making it an attractive target for cancer therapy. MYC overexpression leads to MYC binding at active enhancers, resulting in a global transcriptional amplification of active genes. Because super-enhancers are frequently dysregulated in cancer, we hypothesized that MYC preferentially invades into super-enhancers and alters the cancer genome organization. To that end, we performed ChIP-seq, RNA-seq, circular chromosome conformation capture (4C-seq), and Spike-in Quantitative Hi-C (SIQHiC) on the U2OS osteosarcoma cell line with tetracycline-inducible MYC MYC overexpression in U2OS cells modulated histone acetylation and increased MYC binding at super-enhancers. SIQHiC analysis revealed increased global chromatin contact frequency, particularly at chromatin interactions connecting MYC binding sites at promoters and enhancers. Immunofluorescence staining showed that MYC molecules formed punctate foci at these transcriptionally active domains after MYC overexpression. These results demonstrate the accumulation of overexpressed MYC at promoter-enhancer hubs and suggest that MYC invades into enhancers through spatial proximity. At the same time, the increased protein-protein interactions may strengthen these chromatin interactions to increase chromatin contact frequency. CTCF siRNA knockdown in MYC-overexpressed U2OS cells demonstrated that removal of architectural proteins can disperse MYC and abrogate the increase in chromatin contacts. By elucidating the chromatin landscape of MYC-driven cancers, we can potentially target MYC-associated chromatin interactions for cancer therapy.

Super Enhancer-Mediated Upregulation of HJURP Promotes Growth and Survival of t(4;14)-Positive Multiple Myeloma.

Multiple myeloma is an incurable malignancy with marked clinical and genetic heterogeneity. The cytogenetic abnormality t(4;14) (p16.3;q32.3) confers aggressive behavior in multiple myeloma. Recently, essential oncogenic drivers in a wide range of cancers have been shown to be controlled by super-enhancers (SE). We used chromatin immunoprecipitation sequencing of the active enhancer marker histone H3 lysine 27 acetylation (H3K27ac) to profile unique SEs in t(4;14)-translocated multiple myeloma. The histone chaperone HJURP was aberrantly overexpressed in t(4;14)-positive multiple myeloma due to transcriptional activation by a distal SE induced by the histone lysine methyltransferase NSD2. Silencing of HJURP with short hairpin RNA or CRISPR interference of SE function impaired cell viability and led to apoptosis. Conversely, HJURP overexpression promoted cell proliferation and abrogated apoptosis. Mechanistically, the NSD2/BRD4 complex positively coregulated HJURP transcription by binding the promoter and active elements of its SE. In summary, this study introduces SE profiling as an efficient approach to identify new targets and understand molecular pathogenesis in specific subtypes of cancer. Moreover, HJURP could be a valuable therapeutic target in patients with t(4;14)-positive myeloma. SIGNIFICANCE: A super-enhancer screen in t(4;14) multiple myeloma serves to identify genes that promote growth and survival of myeloma cells, which may be evaluated in future studies as therapeutic targets.

Profiling of 3D Genome Organization in Nasopharyngeal Cancer Needle Biopsy Patient Samples by a Modified Hi-C Approach.

Nasopharyngeal cancer (NPC), a cancer derived from epithelial cells in the nasopharynx, is a cancer common in China, Southeast Asia, and Africa. The three-dimensional (3D) genome organization of nasopharyngeal cancer is poorly understood. A major challenge in understanding the 3D genome organization of cancer samples is the lack of a method for the characterization of chromatin interactions in solid cancer needle biopsy samples. Here, we developed Biop-C, a modified in situ Hi-C method using solid cancer needle biopsy samples. We applied Biop-C to characterize three nasopharyngeal cancer solid cancer needle biopsy patient samples. We identified topologically associated domains (TADs), chromatin interaction loops, and frequently interacting regions (FIREs) at key oncogenes in nasopharyngeal cancer from the Biop-C heatmaps. We observed that the genomic features are shared at some important oncogenes, but the patients also display extensive heterogeneity at certain genomic loci. On analyzing the super enhancer landscape in nasopharyngeal cancer cell lines, we found that the super enhancers are associated with FIREs and can be linked to distal genes via chromatin loops in NPC. Taken together, our results demonstrate the utility of our Biop-C method in investigating 3D genome organization in solid cancers.

Chromatin interaction neural network (ChINN): a machine learning-based method for predicting chromatin interactions from DNA sequences.

Chromatin interactions play important roles in regulating gene expression. However, the availability of genome-wide chromatin interaction data is limited. We develop a computational method, chromatin interaction neural network (ChINN), to predict chromatin interactions between open chromatin regions using only DNA sequences. ChINN predicts CTCF- and RNA polymerase II-associated and Hi-C chromatin interactions. ChINN shows good across-sample performances and captures various sequence features for chromatin interaction prediction. We apply ChINN to 6 chronic lymphocytic leukemia (CLL) patient samples and a published cohort of 84 CLL open chromatin samples. Our results demonstrate extensive heterogeneity in chromatin interactions among CLL patient samples.

Chromatin loop anchors predict transcript and exon usage.

Epigenomics and transcriptomics data from high-throughput sequencing techniques such as RNA-seq and ChIP-seq have been successfully applied in predicting gene transcript expression. However, the locations of chromatin loops in the genome identified by techniques such as Chromatin Interaction Analysis with Paired End Tag sequencing (ChIA-PET) have never been used for prediction tasks. Here, we developed machine learning models to investigate if ChIA-PET could contribute to transcript and exon usage prediction. In doing so, we used a large set of transcription factors as well as ChIA-PET data. We developed different Gradient Boosting Trees models according to the different tasks with the integrated datasets from three cell lines, including GM12878, HeLaS3 and K562. We validated the models via 10-fold cross validation, chromosome-split validation and cross-cell validation. Our results show that both transcript and splicing-derived exon usage can be effectively predicted with at least 0.7512 and 0.7459 of accuracy, respectively, on all cell lines from all kinds of validations. Examining the predictive features, we found that RNA Polymerase II ChIA-PET was one of the most important features in both transcript and exon usage prediction, suggesting that chromatin loop anchors are predictive of both transcript and exon usage.

Interplay and cooperation between SREBF1 and master transcription factors regulate lipid metabolism and tumor-promoting pathways in squamous cancer.

Squamous cell carcinomas (SCCs) comprise one of the most common histologic types of human cancer. Transcriptional dysregulation of SCC cells is orchestrated by tumor protein p63 (TP63), a master transcription factor (TF) and a well-researched SCC-specific oncogene. In the present study, both Gene Set Enrichment Analysis (GSEA) of SCC patient samples and in vitro loss-of-function assays establish fatty-acid metabolism as a key pathway downstream of TP63. Further studies identify sterol regulatory element binding transcription factor 1 (SREBF1) as a central mediator linking TP63 with fatty-acid metabolism, which regulates the biosynthesis of fatty-acids, sphingolipids (SL), and glycerophospholipids (GPL), as revealed by liquid chromatography tandem mass spectrometry (LC-MS/MS)-based lipidomics. Moreover, a feedback co-regulatory loop consisting of SREBF1/TP63/Kruppel like factor 5 (KLF5) is identified, which promotes overexpression of all three TFs in SCCs. Downstream of SREBF1, a non-canonical, SCC-specific function is elucidated: SREBF1 cooperates with TP63/KLF5 to regulate hundreds of cis-regulatory elements across the SCC epigenome, which converge on activating cancer-promoting pathways. Indeed, SREBF1 is essential for SCC viability and migration, and its overexpression is associated with poor survival in SCC patients. Taken together, these data shed light on mechanisms of transcriptional dysregulation in cancer, identify specific epigenetic regulators of lipid metabolism, and uncover SREBF1 as a potential therapeutic target and prognostic marker in SCC.