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.

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.

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.

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 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.

ZNF143 mediates CTCF-bound promoter-enhancer loops required for murine hematopoietic stem and progenitor cell function

CCCTC binding factor (CTCF) is an important factor in the maintenance of chromatin-chromatin interactions, yet the mechanism regulating its binding to chromatin is unknown. We demonstrate that zinc finger protein 143 (ZNF143) is a key regulator for CTCF-bound promoter-enhancer loops. In the murine genome, a large percentage of CTCF and ZNF143 DNA binding motifs are distributed 37 bp apart in the convergent orientation. Furthermore, deletion of ZNF143 leads to loss of CTCF binding on promoter and enhancer regions associated with gene expression changes. CTCF-bound promoter-enhancer loops are also disrupted after excision of ZNF143. ZNF143-CTCF-bound promoter-enhancer loops regulate gene expression patterns essential for maintenance of murine hematopoietic stem and progenitor cell integrity. Our data suggest a common feature of gene regulation is that ZNF143 is a critical factor for CTCF-bound promoter-enhancer loops.

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.

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.

HoxC5 and miR-615-3p target newly evolved genomic regions to repress hTERT and inhibit tumorigenesis

The repression of telomerase activity during cellular differentiation promotes replicative aging and functions as a physiological barrier for tumorigenesis in long-lived mammals, including humans. However, the underlying mechanisms remain largely unclear. Here we describe how miR-615-3p represses hTERT expression. mir-615-3p is located in an intron of the HOXC5 gene, a member of the highly conserved homeobox family of transcription factors controlling embryogenesis and development. Unexpectedly, we found that HoxC5 also represses hTERT expression by disrupting the long-range interaction between hTERT promoter and its distal enhancer. The 3'UTR of hTERT and its upstream enhancer region are well conserved in long-lived primates. Both mir-615-3p and HOXC5 are activated upon differentiation, which constitute a feed-forward loop that coordinates transcriptional and post-transcriptional repression of hTERT during cellular differentiation. Deregulation of HOXC5 and mir-615-3p expression may contribute to the activation of hTERT in human cancers.

Long-Range Chromatin Interactions Drive Mutant TERT Promoter Activation

Cancer-specific TERT promoter mutations (-146C>T and -124C>T) have been linked to reactivation of the epigenetically silenced telomerase reverse transcriptase gene (TERT). Understanding how these single-nucleotide alterations drive TERT reactivation is a fundamental unanswered question and is key for making successful therapeutics. We show that unlike wild-type promoters, recruitment of the transcription factor GABPA specifically to mutant TERT promoters mediates long-range chromatin interaction and enrichment of active histone marks, and hence drives TERT transcription. CRISPR-mediated reversal of mutant TERT promoters, or deletion of its long-range interacting chromatin, abrogates GABPA binding and long-range interactions, leading to depletion of active histone marks, loss of POL2 recruitment, and suppression of TERT transcription. In contrast, de novo introduction of a TERT promoter mutation enables GABPA binding and upregulation of TERT via long-range interactions, acquisition of active histone marks, and subsequent POL2 recruitment. This study provides a unifying mechanistic insight into activation of mutant TERT promoters across various human cancers.

SIGNIFICANCE

This study identifies a key mechanism by which cancer-specific mutant TERT promoters cause reactivation of TERT Because the mechanism uncovered here is not utilized by promoters that drive TERT in normal cells, this mechanism could be exploited to make inhibitors which have the potential to block telomerase function and hence the progression of up to 90% of human cancers. Cancer Discov; 6(11); 1276-91. ©2016 AACR.See related commentary by Min and Shay, p. 1212This article is highlighted in the In This Issue feature, p. 1197.

Large-scale functional organization of long-range chromatin interaction networks

Chromatin interactions play important roles in transcription regulation. To better understand the underlying evolutionary and functional constraints of these interactions, we implemented a systems approach to examine RNA polymerase-II-associated chromatin interactions in human cells. We found that 40% of the total genomic elements involved in chromatin interactions converged to a giant, scale-free-like, hierarchical network organized into chromatin communities. The communities were enriched in specific functions and were syntenic through evolution. Disease-associated SNPs from genome-wide association studies were enriched among the nodes with fewer interactions, implying their selection against deleterious interactions by limiting the total number of interactions, a model that we further reconciled using somatic and germline cancer mutation data. The hubs lacked disease-associated SNPs, constituted a nonrandomly interconnected core of key cellular functions, and exhibited lethality in mouse mutants, supporting an evolutionary selection that favored the nonrandom spatial clustering of the least-evolving key genomic domains against random genetic or transcriptional errors in the genome. Altogether, our analyses reveal a systems-level evolutionary framework that shapes functionally compartmentalized and error-tolerant transcriptional regulation of human genome in three dimensions.