Epigenetic modifiers frequently harbor loss-of-function mutations in lung cancer, but their tumor-suppressive roles are poorly characterized. Histone methyltransferase KMT2D (a COMPASS-like enzyme, ...also called MLL4) is among the most highly inactivated epigenetic modifiers in lung cancer. Here, we show that lung-specific loss of Kmt2d promotes lung tumorigenesis in mice and upregulates pro-tumorigenic programs, including glycolysis. Pharmacological inhibition of glycolysis preferentially impedes tumorigenicity of human lung cancer cells bearing KMT2D-inactivating mutations. Mechanistically, Kmt2d loss widely impairs epigenomic signals for super-enhancers/enhancers, including the super-enhancer for the circadian rhythm repressor Per2. Loss of Kmt2d decreases expression of PER2, which regulates multiple glycolytic genes. These findings indicate that KMT2D is a lung tumor suppressor and that KMT2D deficiency confers a therapeutic vulnerability to glycolytic inhibitors.
Display omitted
•Lung-specific Kmt2d loss in mice promotes lung tumorigenesis•Kmt2d loss impairs enhancers, including a super-enhancer for the tumor suppressor Per2•KMT2D activates Per2 expression and thereby represses glycolytic genes•Glycolysis inhibition impedes the growth of KMT2D-mutant lung cancer
Histone methyltransferase KMT2D is frequently mutated in lung tumors, and Alam et al. identify KMT2D as a lung tumor suppressor. KMT2D deficiency induces aberrant metabolic reprogramming via super-enhancer impairment, conferring sensitivity to glycolytic inhibitors in lung cancer with KMT2D-inactivating mutations.
The roles of Plant Homeodomain (PHD) fingers in catalysis of histone modifications are unknown. We demonstrated that the PHD finger of Ubiquitin Protein Ligase E3 Component N-Recognin7 (UBR7) harbors ...E3 ubiquitin ligase activity toward monoubiquitination of histone H2B at lysine120 (H2BK120Ub). Purified PHD finger or full-length UBR7 monoubiquitinated H2BK120 in vitro, and loss of UBR7 drastically reduced H2BK120Ub genome-wide binding sites in MCF10A cells. Low UBR7 expression was correlated with occurrence of triple-negative breast cancer and metastatic tumors. Consistently, UBR7 knockdown enhanced the invasiveness, induced epithelial-to-mesenchymal transition and promoted metastasis. Conversely, ectopic expression of UBR7 restored these cellular phenotypes and reduced tumor growth. Mechanistically, UBR7 loss reduced H2BK120Ub levels on cell adhesion genes, including CDH4, and upregulated the Wnt/β-Catenin signaling pathway. CDH4 overexpression could partially revert UBR7-dependent cellular phenotypes. Collectively, our results established UBR7 as a histone H2B monoubiquitin ligase that suppresses tumorigenesis and metastasis of triple-negative breast cancer.
A major goal of biomedical research is to identify molecular features associated with a biological or clinical class of interest. Differential expression analysis has long been used for this purpose; ...however, conventional methods perform poorly when applied to data with high within class heterogeneity.
To address this challenge, we developed EMDomics, a new method that uses the Earth mover's distance to measure the overall difference between the distributions of a gene's expression in two classes of samples and uses permutations to obtain q-values for each gene. We applied EMDomics to the challenging problem of identifying genes associated with drug resistance in ovarian cancer. We also used simulated data to evaluate the performance of EMDomics, in terms of sensitivity and specificity for identifying differentially expressed gene in classes with high within class heterogeneity. In both the simulated and real biological data, EMDomics outperformed competing approaches for the identification of differentially expressed genes, and EMDomics was significantly more powerful than conventional methods for the identification of drug resistance-associated gene sets. EMDomics represents a new approach for the identification of genes differentially expressed between heterogeneous classes and has utility in a wide range of complex biomedical conditions in which sample classes show within class heterogeneity.
The R package is available at http://www.bioconductor.org/packages/release/bioc/html/EMDomics.html.
We reported that histone H3 lysine (K) 4 methyltransferase, KMT2D, serves as a potent tumor-suppressor in melanoma, which was identified via in vivo epigenome-focused RNA interference (RNAi) screen. ...KMT2D-deficient tumors show substantial reprogramming of key metabolic pathways including glycolysis via reduction of H3K4me1 (Histone H3K4 mono-methylation)-marked active enhancers, conferring sensitivity to inhibitors of glycolysis and IGFR (Insulin Growth Factor Receptor) pathway.
Desmoplastic small round cell tumor (DSRCT) is an aggressive, usually incurable sarcoma subtype that predominantly occurs in post-pubertal young males. Recent evidence suggests that the androgen ...receptor (AR) can promote tumor progression in DSRCTs. However, the mechanism of AR-induced oncogenic stimulation remains undetermined. Herein, we demonstrate that enzalutamide and AR-directed antisense oligonucleotides (AR-ASO) block 5α-dihydrotestosterone (DHT)-induced DSRCT cell proliferation and reduce xenograft tumor burden. Gene expression analysis and chromatin immunoprecipitation sequencing (ChIP-seq) were performed to elucidate how AR signaling regulates cellular epigenetic programs. Remarkably, ChIP-seq revealed novel DSRCT-specific AR DNA binding sites adjacent to key oncogenic regulators, including WT1 (the C-terminal partner of the pathognomonic fusion protein) and FOXF1. Additionally, AR occupied enhancer sites that regulate the Wnt pathway, neural differentiation, and embryonic organ development, implicating AR in dysfunctional cell lineage commitment. Our findings have direct clinical implications given the widespread availability of FDA-approved androgen-targeted agents used for prostate cancer.
Histone methyltransferase KMT2D harbors frequent loss-of-function somatic point mutations in several tumor types, including melanoma. Here, we identify KMT2D as a potent tumor suppressor in melanoma ...through an in vivo epigenome-focused pooled RNAi screen and confirm the finding by using a genetically engineered mouse model (GEMM) based on conditional and melanocyte-specific deletion of KMT2D. KMT2D-deficient tumors show substantial reprogramming of key metabolic pathways, including glycolysis. KMT2D deficiency aberrantly upregulates glycolysis enzymes, intermediate metabolites, and glucose consumption rates. Mechanistically, KMT2D loss causes genome-wide reduction of H3K4me1-marked active enhancer chromatin states. Enhancer loss and subsequent repression of IGFBP5 activates IGF1R-AKT to increase glycolysis in KMT2D-deficient cells. Pharmacological inhibition of glycolysis and insulin growth factor (IGF) signaling reduce proliferation and tumorigenesis preferentially in KMT2D-deficient cells. We conclude that KMT2D loss promotes tumorigenesis by facilitating an increased use of the glycolysis pathway for enhanced biomass needs via enhancer reprogramming, thus presenting an opportunity for therapeutic intervention through glycolysis or IGF pathway inhibitors.
Display omitted
•KMT2D is a tumor suppressor in melanoma•KMT2D rewires metabolic pathways through enhancer reprogramming•KMT2D loss impairs IGFBP5 enhancers and thereby deprives repression to glycolytic genes•KMT2D mutant melanomas are preferentially sensitive to glycolysis and IGFR inhibition
Through an in vivo epigenome-focused pooled RNAi screen, Maitituoheti et al. identify KMT2D as a tumor suppressor in melanoma. KMT2D-deficient tumors show substantial reprogramming of key metabolic pathways by reduction of H3K4me1-marked active enhancers, conferring sensitivity to glycolysis and IGFR inhibitors in melanoma with KMT2D-inactivating mutations.
The extent and nature of epigenomic changes associated with melanoma progression is poorly understood. Through systematic epigenomic profiling of 35 epigenetic modifications and transcriptomic ...analysis, we define chromatin state changes associated with melanomagenesis by using a cell phenotypic model of non-tumorigenic and tumorigenic states. Computation of specific chromatin state transitions showed loss of histone acetylations and H3K4me2/3 on regulatory regions proximal to specific cancer-regulatory genes in important melanoma-driving cell signaling pathways. Importantly, such acetylation changes were also observed between benign nevi and malignant melanoma human tissues. Intriguingly, only a small fraction of chromatin state transitions correlated with expected changes in gene expression patterns. Restoration of acetylation levels on deacetylated loci by histone deacetylase (HDAC) inhibitors selectively blocked excessive proliferation in tumorigenic cells and human melanoma cells, suggesting functional roles of observed chromatin state transitions in driving hyperproliferative phenotype. Through these results, we define functionally relevant chromatin states associated with melanoma progression.
Display omitted
•Epigenomic alterations in non-tumorigenic to tumorigenic transition in melanoma•Chromatin states harboring acetylation and H3K4me2/3 are lost in tumorigenic cells•Chromatin state transitions preferentially occur on melanoma-regulatory pathways•HDAC inhibition preferentially impacts proliferative ability of tumorigenic cells
Using comprehensive profiling of 35 epigenetic marks and determination of chromatin state transitions between non-tumorigenic and tumorigenic systems, Fiziev et al. find that, in tumorigenic cells, loss of histone acetylation and H3K4 methylation occur on regulatory regions proximal to specific cancer-regulatory genes.
Histone modifications constitute a major component of the epigenome and play important regulatory roles in determining the transcriptional status of associated loci. In addition, the presence of ...specific modifications has been used to determine the position and identity non-coding functional elements such as enhancers. In recent years, chromatin immunoprecipitation followed by next generation sequencing (ChIP-seq) has become a powerful tool in determining the genome-wide profiles of individual histone modifications. However, it has become increasingly clear that the combinatorial patterns of chromatin modifications, referred to as Chromatin States, determine the identity and nature of the associated genomic locus. Therefore, workflows consisting of robust high-throughput (HT) methodologies for profiling a number of histone modification marks, as well as computational analyses pipelines capable of handling myriads of ChIP-Seq profiling datasets, are needed for comprehensive determination of epigenomic states in large number of samples. The HT-ChIP-Seq workflow presented here consists of two modules: 1) an experimental protocol for profiling several histone modifications from small amounts of tumor samples and cell lines in a 96-well format; and 2) a computational data analysis pipeline that combines existing tools to compute both individual mark occupancy and combinatorial chromatin state patterns. Together, these two modules facilitate easy processing of hundreds of ChIP-Seq samples in a fast and efficient manner. The workflow presented here is used to derive chromatin state patterns from 6 histone mark profiles in melanoma tumors and cell lines. Overall, we present a comprehensive ChIP-seq workflow that can be applied to dozens of human tumor samples and cancer cell lines to determine epigenomic aberrations in various malignancies.
Enhancer aberrations are beginning to emerge as a key epigenetic feature of colorectal cancers (CRC), however, a comprehensive knowledge of chromatin state patterns in tumour progression, ...heterogeneity of these patterns and imparted therapeutic opportunities remain poorly described.
We performed comprehensive epigenomic characterisation by mapping 222 chromatin profiles from 69 samples (33 colorectal adenocarcinomas, 4 adenomas, 21 matched normal tissues and 11 colon cancer cell lines) for six histone modification marks: H3K4me3 for Pol II-bound and CpG-rich promoters, H3K4me1 for poised enhancers, H3K27ac for enhancers and transcriptionally active promoters, H3K79me2 for transcribed regions, H3K27me3 for polycomb repressed regions and H3K9me3 for heterochromatin.
We demonstrate that H3K27ac-marked active enhancer state could distinguish between different stages of CRC progression. By epigenomic editing, we present evidence that gains of tumour-specific enhancers for crucial oncogenes, such as
and
was required for excessive proliferation. Consistently, combination of MEK plus bromodomain inhibition was found to have synergistic effects in CRC patient-derived xenograft models. Probing intertumour heterogeneity, we identified four distinct enhancer subtypes (EPIgenome-based Classification, EpiC), three of which correlate well with previously defined transcriptomic subtypes (consensus molecular subtypes, CMSs). Importantly, CMS2 can be divided into two EpiC subgroups with significant survival differences. Leveraging such correlation, we devised a combinatorial therapeutic strategy of enhancer-blocking bromodomain inhibitors with pathway-specific inhibitors (PARPi, EGFRi, TGFβi, mTORi and SRCi) for EpiC groups.
Our data suggest that the dynamics of active enhancer underlies CRC progression and the patient-specific enhancer patterns can be leveraged for precision combination therapy.