Extensive promoter-centered chromatin interactions provide a topological basis for transcription regulation.

Higher-order chromosomal organization for transcription regulation is poorly understood in eukaryotes. Using genome-wide Chromatin Interaction Analysis with Paired-End-Tag sequencing (ChIA-PET), we mapped long-range chromatin interactions associated with RNA polymerase II in human cells and uncovered widespread promoter-centered intragenic, extragenic, and intergenic interactions. These interactions further aggregated into higher-order clusters, wherein proximal and distal genes were engaged through promoter-promoter interactions. Most genes with promoter-promoter interactions were active and transcribed cooperatively, and some interacting promoters could influence each other implying combinatorial complexity of transcriptional controls. Comparative analyses of different cell lines showed that cell-specific chromatin interactions could provide structural frameworks for cell-specific transcription, and suggested significant enrichment of enhancer-promoter interactions for cell-specific functions. Furthermore, genetically-identified disease-associated noncoding elements were found to be spatially engaged with corresponding genes through long-range interactions. Overall, our study provides insights into transcription regulation by three-dimensional chromatin interactions for both housekeeping and cell-specific genes in human cells.

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.

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.

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.

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.

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.

DeepCPP: a deep neural network based on nucleotide bias information and minimum distribution similarity feature selection for RNA coding potential prediction.

The development of deep sequencing technologies has led to the discovery of novel transcripts. Many in silico methods have been developed to assess the coding potential of these transcripts to further investigate their functions. Existing methods perform well on distinguishing majority long noncoding RNAs (lncRNAs) and coding RNAs (mRNAs) but poorly on RNAs with small open reading frames (sORFs). Here, we present DeepCPP (deep neural network for coding potential prediction), a deep learning method for RNA coding potential prediction. Extensive evaluations on four previous datasets and six new datasets constructed in different species show that DeepCPP outperforms other state-of-the-art methods, especially on sORF type data, which overcomes the bottleneck of sORF mRNA identification by improving more than 4.31, 37.24 and 5.89% on its accuracy for newly discovered human, vertebrate and insect data, respectively. Additionally, we also revealed that discontinuous k-mer, and our newly proposed nucleotide bias and minimal distribution similarity feature selection method play crucial roles in this classification problem. Taken together, DeepCPP is an effective method for RNA coding potential prediction.

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.

Chromatin Interactions and Regulatory Elements in Cancer: From Bench to Bedside.

Chromatin interactions regulate gene expression by bringing distal regulatory elements, such as super-enhancers, to promoters in close spatial proximity. It has been recognized that in cancer, chromatin interactions can be dysregulated, leading to aberrant oncogene expression. Chromatin interactions may potentially serve as biomarkers, or be modulated via CRISPR therapy and small molecule inhibitors against transcription. However, these methods face challenges that must be resolved and raise questions for further research. Understanding chromatin interactions is essential for safety aspects of anticancer therapies, such as the mechanism of action of epigenetic regulators and transcription factors in cancer, and potential off-target effects arising from targeting super-enhancers and promoters. In this review article, we discuss how chromatin interactions and regulatory elements may become dysregulated in cancer, potential methods to target them for clinical therapy, and outline outstanding questions that require addressing before epigenetic therapies can translate to the clinic safely and effectively.

TP63, SOX2, and KLF5 Establish a Core Regulatory Circuitry That Controls Epigenetic and Transcription Patterns in Esophageal Squamous Cell Carcinoma Cell Lines.

BACKGROUND & AIMS: We investigated the transcriptome of esophageal squamous cell carcinoma (ESCC) cells, activity of gene regulatory (enhancer and promoter regions), and the effects of blocking epigenetic regulatory proteins. METHODS: We performed chromatin immunoprecipitation sequencing with antibodies against H3K4me1, H3K4me3, and H3K27ac and an assay for transposase-accessible chromatin to map the enhancer regions and accessible chromatin in 8 ESCC cell lines. We used the CRC_Mapper algorithm to identify core regulatory circuitry transcription factors in ESCC cell lines, and determined genome occupancy profiles for 3 of these factors. In ESCC cell lines, expression of transcription factors was knocked down with small hairpin RNAs, promoter and enhancer regions were disrupted by CRISPR/Cas9 genome editing, or bromodomains and extraterminal (BET) family proteins and histone deacetylases (HDACs) were inhibited with ARV-771 and romidepsin, respectively. ESCC cell lines were then analyzed by whole-transcriptome sequencing, immunoprecipitation, immunoblots, immunohistochemistry, and viability assays. Interactions between distal enhancers and promoters were identified and verified with circular chromosome conformation capture sequencing. NOD-SCID mice were given injections of modified ESCC cells, some mice where given injections of HDAC or BET inhibitors, and growth of xenograft tumors was measured. RESULTS: We identified super-enhancer-regulated circuits and transcription factors TP63, SOX2, and KLF5 as core regulatory factors in ESCC cells. Super-enhancer regulation of ALDH3A1 mediated by core regulatory factors was required for ESCC viability. We observed direct interactions between the promoter region of TP63 and functional enhancers, mediated by the core regulatory circuitry transcription factors. Deletion of enhancer regions from ESCC cells decreased expression of the core regulatory circuitry transcription factors and reduced cell viability; these same results were observed with knockdown of each core regulatory circuitry transcription factor. Incubation of ESCC cells with BET and HDAC disrupted the core regulatory circuitry program and the epigenetic modifications observed in these cells; mice given injections of HDAC or BET inhibitors developed smaller xenograft tumors from the ESCC cell lines. Xenograft tumors grew more slowly in mice given the combination of ARV-771 and romidepsin than mice given either agent alone. CONCLUSIONS: In epigenetic and transcriptional analyses of ESCC cell lines, we found the transcription factors TP63, SOX2, and KLF5 to be part of a core regulatory network that determines chromatin accessibility, epigenetic modifications, and gene expression patterns in these cells. A combination of epigenetic inhibitors slowed growth of xenograft tumors derived from ESCC cells in mice.

Identification of a novel enhancer of CEBPE essential for granulocytic differentiation.

CCAAT/enhancer binding protein ε (CEBPE) is an essential transcription factor for granulocytic differentiation. Mutations of CEBPE occur in individuals with neutrophil-specific granule deficiency (SGD), which is characterized by defects in neutrophil maturation. Cebpe-knockout mice also exhibit defects in terminal differentiation of granulocytes, a phenotype reminiscent of SGD. Analysis of DNase I hypersensitive sites sequencing data revealed an open chromatin region 6 kb downstream of the transcriptional start site of Cebpe in murine myeloid cells. We identified an interaction between this +6-kb region and the core promoter of Cebpe using circular chromosome conformation capture sequencing (4C-seq). To understand the role of this putative enhancer in transcriptional regulation of Cebpe, we targeted it using catalytically inactive Cas9 fused to Krüppel-associated box (KRAB) domain and observed a significant downregulation of transcript and protein levels of CEBPE in cells expressing guide RNA targeting the +6-kb region. To further investigate the role of this novel enhancer further in myelopoiesis, we generated mice with deletion of this region using CRISPR/Cas9 technology. Germline deletion of the +6-kb enhancer resulted in reduced levels of CEBPE and its target genes and caused a severe block in granulocytic differentiation. We also identified binding of CEBPA and CEBPE to the +6-kb enhancer, which suggests their role in regulating the expression of Cebpe In summary, we have identified a novel enhancer crucial for regulating expression of Cebpe and required for normal granulocytic differentiation.

The enhancer RNA ARIEL activates the oncogenic transcriptional program in T-cell acute lymphoblastic leukemia.

The oncogenic transcription factor TAL1 regulates the transcriptional program in T-ALL. ARID5B is one of the critical downstream targets of TAL1, which further activates the oncogenic regulatory circuit in T-ALL cells. Here, we elucidated the molecular functions of the noncoding RNA, ARID5B-inducing enhancer associated long noncoding RNA (ARIEL), in T-ALL pathogenesis. We demonstrated that ARIEL is specifically activated in TAL1 + T-ALL cases, and its expression is associated with ARID5B enhancer activity. ARIEL recruits mediator proteins to the ARID5B enhancer, promotes enhancer-promoter interactions, and activates the expression of ARID5B, thereby positively regulating the TAL1-induced transcriptional program and the MYC oncogene. The TAL1 complex coordinately regulates the expression of ARIEL Knockdown of ARIEL inhibits cell growth and survival of T-ALL cells in culture and blocks disease progression in a murine xenograft model. Our results indicate that ARIEL plays an oncogenic role as an enhancer RNA in T-ALL.

Robust CTCF-Based Chromatin Architecture Underpins Epigenetic Changes in the Heart Failure Stress-Gene Response.

BACKGROUND: The human genome folds in 3 dimensions to form thousands of chromatin loops inside the nucleus, encasing genes and cis-regulatory elements for accurate gene expression control. Physical tethers of loops are anchored by the DNA-binding protein CTCF and the cohesin ring complex. Because heart failure is characterized by hallmark gene expression changes, it was recently reported that substantial CTCF-related chromatin reorganization underpins the myocardial stress-gene response, paralleled by chromatin domain boundary changes observed in CTCF knockout. METHODS: We undertook an independent and orthogonal analysis of chromatin organization with mouse pressure-overload model of myocardial stress (transverse aortic constriction) and cardiomyocyte-specific knockout of Ctcf. We also downloaded published data sets of similar cardiac mouse models and subjected them to independent reanalysis. RESULTS: We found that the cardiomyocyte chromatin architecture remains broadly stable in transverse aortic constriction hearts, whereas Ctcf knockout resulted in ≈99% abolition of global chromatin loops. Disease gene expression changes correlated instead with differential histone H3K27-acetylation enrichment at their respective proximal and distal interacting genomic enhancers confined within these static chromatin structures. Moreover, coregulated genes were mapped out as interconnected gene sets on the basis of their multigene 3D interactions. CONCLUSIONS: This work reveals a more stable genome-wide chromatin framework than previously described. Myocardial stress-gene transcription responds instead through H3K27-acetylation enhancer enrichment dynamics and gene networks of coregulation. Robust and intact CTCF looping is required for the induction of a rapid and accurate stress response.

Aberrant hyperediting of the myeloma transcriptome by ADAR1 confers oncogenicity and is a marker of poor prognosis.

DNA alterations have been extensively reported in multiple myeloma (MM); however, they cannot yet fully explain all the biological and molecular abnormalities in MM, which remains to this day an incurable disease with eventual emergence of refractory disease. Recent years have seen abnormalities at the RNA levels being reported to possess potential biological relevance in cancers. ADAR1-mediated A-to-I editing is an important posttranscriptional mechanism in human physiology, and the biological implication of its abnormality, especially at the global level, is underexplored in MM. In this study, we define the biological implications of A-to-I editing and how it contributes to MM pathogenesis. Here, we identified that the MM transcriptome is aberrantly hyperedited because of the overexpression of ADAR1. These events were associated with patients' survival independent of 1q21 amplifications and could affect patients' responsiveness to different treatment regimes. Our functional assays established ADAR1 to be oncogenic, driving cellular growth and proliferation in an editing-dependent manner. In addition, we identified NEIL1 (base-excision repair gene) as an essential and a ubiquitously edited ADAR1 target in MM. The recoded NEIL1 protein showed defective oxidative damage repair capacity and loss-of-function properties. Collectively, our data demonstrated that ADAR1-mediated A-to-I editing is both clinically and biologically relevant in MM. These data unraveled novel insights into MM molecular pathogenesis at the global RNA level.

Super-Enhancer-Driven Long Non-Coding RNA LINC01503, Regulated by TP63, Is Over-Expressed and Oncogenic in Squamous Cell Carcinoma.

BACKGROUND & AIMS: Long non-coding RNAs (lncRNAs) are expressed in tissue-specific pattern, but it is not clear how these are regulated. We aimed to identify squamous cell carcinoma (SCC)-specific lncRNAs and investigate mechanisms that control their expression and function. METHODS: We studied expression patterns and functions of 4 SCC-specific lncRNAs. We obtained 113 esophageal SCC (ESCC) and matched non-tumor esophageal tissues from a hospital in Shantou City, China, and performed quantitative reverse transcription polymerase chain reaction assays to measure expression levels of LINC01503. We collected clinical data from patients and compared expression levels with survival times. LINC01503 was knocked down using small interfering RNAs and oligonucleotides in TE7, TE5, and KYSE510 cell lines and overexpressed in KYSE30 cells. Cells were analyzed by chromatin immunoprecipitation sequencing, luciferase reporter assays, colony formation, migration and invasion, and mass spectrometry analyses. Cells were injected into nude mice and growth of xenograft tumors was measured. LINC01503 interaction with proteins was studied using fluorescence in situ hybridization, RNA pulldown, and RNA immunoprecipitation analyses. RESULTS: We identified a lncRNA, LINC01503, which is regulated by a super enhancer and is expressed at significantly higher levels in esophageal and head and neck SCCs than in non-tumor tissues. High levels in SCCs correlated with shorter survival times of patients. The transcription factor TP63 bound to the super enhancer at the LINC01503 locus and activated its transcription. Expression of LINC01503 in ESCC cell lines increased their proliferation, colony formation, migration, and invasion. Knockdown of LINC01503 in SCC cells reduced their proliferation, colony formation, migration, and invasion, and the growth of xenograft tumors in nude mice. Expression of LINC01503 in ESCC cell lines reduced ERK2 dephosphorylation by DUSP6, leading to activation of ERK signaling via MAPK. LINC01503 disrupted the interaction between EBP1 and the p85 subunit of PI3K, increasing AKT signaling. CONCLUSIONS: We identified an lncRNA, LINC01503, which is increased in SCC cells compared with non-tumor cells. Increased expression of LINC01503 promotes ESCC cell proliferation, migration, invasion, and growth of xenograft tumors. It might be developed as a biomarker of aggressive SCCs in patients.

Landscape of transcription in human cells.

Eukaryotic cells make many types of primary and processed RNAs that are found either in specific subcellular compartments or throughout the cells. A complete catalogue of these RNAs is not yet available and their characteristic subcellular localizations are also poorly understood. Because RNA represents the direct output of the genetic information encoded by genomes and a significant proportion of a cell's regulatory capabilities are focused on its synthesis, processing, transport, modification and translation, the generation of such a catalogue is crucial for understanding genome function. Here we report evidence that three-quarters of the human genome is capable of being transcribed, as well as observations about the range and levels of expression, localization, processing fates, regulatory regions and modifications of almost all currently annotated and thousands of previously unannotated RNAs. These observations, taken together, prompt a redefinition of the concept of a gene.

An oestrogen-receptor-alpha-bound human chromatin interactome.

Genomes are organized into high-level three-dimensional structures, and DNA elements separated by long genomic distances can in principle interact functionally. Many transcription factors bind to regulatory DNA elements distant from gene promoters. Although distal binding sites have been shown to regulate transcription by long-range chromatin interactions at a few loci, chromatin interactions and their impact on transcription regulation have not been investigated in a genome-wide manner. Here we describe the development of a new strategy, chromatin interaction analysis by paired-end tag sequencing (ChIA-PET) for the de novo detection of global chromatin interactions, with which we have comprehensively mapped the chromatin interaction network bound by oestrogen receptor alpha (ER-alpha) in the human genome. We found that most high-confidence remote ER-alpha-binding sites are anchored at gene promoters through long-range chromatin interactions, suggesting that ER-alpha functions by extensive chromatin looping to bring genes together for coordinated transcriptional regulation. We propose that chromatin interactions constitute a primary mechanism for regulating transcription in mammalian genomes.

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.