Search

CHARGE syndrome, a rare multiple congenital anomaly condition, is caused by haploinsufficiency of the chromatin remodeling protein gene CHD7 (Chromodomain helicase DNA binding protein 7). Brain abnormalities and intellectual disability are commonly observed in individuals with CHARGE, and neuronal differentiation is reduced in CHARGE patient-derived iPSCs and conditional knockout mouse brains. However, the mechanisms of CHD7 function in nervous system development are not well understood. In this study, we asked whether CHD7 promotes gene transcription in neural progenitor cells via changes in chromatin accessibility. We used Chd7 null embryonic stem cells (ESCs) derived from Chd7 mutant mouse blastocysts as a tool to investigate roles of CHD7 in neuronal and glial differentiation. Loss of Chd7 significantly reduced neuronal and glial differentiation. Sholl analysis showed that loss of Chd7 impaired neuronal complexity and neurite length in differentiated neurons. Genome-wide studies demonstrated that loss of Chd7 leads to modified chromatin accessibility (ATAC-seq) and differential nascent expression (Bru-Seq) of neural-specific genes. These results suggest that CHD7 acts preferentially to alter chromatin accessibility of key genes during the transition of NPCs to neurons to promote differentiation. Our results form a basis for understanding the cell stage-specific roles for CHD7-mediated chromatin remodeling during cell lineage acquisition.

The degree of intrinsic and interpatient phenotypic heterogeneity and its role in tumor evolution is poorly understood. Phenotypic drifts can be transmitted via inheritable transcriptional programs. Cell-type specific transcription is maintained through the activation of epigenetically defined regulatory regions including promoters and enhancers. Here we have annotated the epigenome of 47 primary and metastatic estrogen-receptor (ERα)-positive breast cancer clinical specimens and inferred phenotypic heterogeneity from the regulatory landscape, identifying key regulatory elements commonly shared across patients. Shared regions contain a unique set of regulatory information including the motif for transcription factor YY1. We identify YY1 as a critical determinant of ERα transcriptional activity promoting tumor growth in most luminal patients. YY1 also contributes to the expression of genes mediating resistance to endocrine treatment. Finally, we used H3K27ac levels at active enhancer elements as a surrogate of intra-tumor phenotypic heterogeneity to track the expansion and contraction of phenotypic subpopulations throughout breast cancer progression. By tracking the clonality of SLC9A3R1-positive cells, a bona fide YY1-ERα-regulated gene, we show that endocrine therapies select for phenotypic clones under-represented at diagnosis. Collectively, our data show that epigenetic mechanisms significantly contribute to phenotypic heterogeneity and evolution in systemically treated breast cancer patients.

A hallmark of craniofacial development is the differentiation of multiple cell lineages in close proximity to one another. The mouse skull bones and overlying dermis are derived from the cranial mesenchyme (CM). Cell fate selection of the embryonic cranial bone and dermis in the CM requires Wnt/β-catenin signaling, and loss of β-catenin leads to an ectopic chondrogenic cell fate switch. The mechanism by which Wnt/β-catenin activity suppresses the cartilage fate is unclear. Upon conditional deletion of β-catenin in the CM, several key determinants of the cartilage differentiation program, including Sox9, become differentially expressed. Many of these differentially expressed genes are known targets of the Polycomb Repressive Complex 2 (PRC2). Thus, we hypothesized that PRC2 is required for Wnt/β-catenin-mediated repression of chondrogenesis in the embryonic CM. We find that β-catenin can physically interact with PRC2 components in the CM in vivo. However, upon genetic deletion of Enhancer of Zeste homolog 2 (EZH2), the catalytic component of PRC2, chondrogenesis remains repressed and the bone and dermis cell fate is preserved in the CM. Furthermore, loss of β-catenin does not alter either the H3K27me3 enrichment levels genome-wide or on cartilage differentiation determinants, including Sox9. Our results indicate that EZH2 is not required to repress chondrogenesis in the CM downstream of Wnt/β-catenin signaling.

Metastasis results from a complex set of traits acquired by tumor cells, distinct from those necessary for tumorigenesis. Here, we investigate the contribution of enhancer elements to the metastatic phenotype of osteosarcoma. Through epigenomic profiling, we identify substantial differences in enhancer activity between primary and metastatic human tumors and between near isogenic pairs of highly lung metastatic and nonmetastatic osteosarcoma cell lines. We term these regions metastatic variant enhancer loci (Met-VELs). Met-VELs drive coordinated waves of gene expression during metastatic colonization of the lung. Met-VELs cluster nonrandomly in the genome, indicating that activity of these enhancers and expression of their associated gene targets are positively selected. As evidence of this causal association, osteosarcoma lung metastasis is inhibited by global interruptions of Met-VEL-associated gene expression via pharmacologic BET inhibition, by knockdown of AP-1 transcription factors that occupy Met-VELs, and by knockdown or functional inhibition of individual genes activated by Met-VELs, such as that encoding coagulation factor III/tissue factor (F3). We further show that genetic deletion of a single Met-VEL at the F3 locus blocks metastatic cell outgrowth in the lung. These findings indicate that Met-VELs and the genes they regulate play a functional role in metastasis and may be suitable targets for antimetastatic therapies.

Genomic sequencing has driven precision-based oncology therapy; however, the genetic drivers of many malignancies remain unknown or non-targetable, so alternative approaches to the identification of therapeutic leads are necessary. Ependymomas are chemotherapy-resistant brain tumours, which, despite genomic sequencing, lack effective molecular targets. Intracranial ependymomas are segregated on the basis of anatomical location (supratentorial region or posterior fossa) and further divided into distinct molecular subgroups that reflect differences in the age of onset, gender predominance and response to therapy1,2,3. The most common and aggressive subgroup, posterior fossa ependymoma group A (PF-EPN-A), occurs in young children and appears to lack recurrent somatic mutations2. Conversely, posterior fossa ependymoma group B (PF-EPN-B) tumours display frequent large-scale copy number gains and losses but have favourable clinical outcomes1,3. More than 70% of supratentorial ependymomas are defined by highly recurrent gene fusions in the NF-κB subunit gene RELA (ST-EPN-RELA), and a smaller number involve fusion of the gene encoding the transcriptional activator YAP1 (ST-EPN-YAP1)1,3,4. Subependymomas, a distinct histologic variant, can also be found within the supratetorial and posterior fossa compartments, and account for the majority of tumours in the molecular subgroups ST-EPN-SE and PF-EPN-SE. Here we describe mapping of active chromatin landscapes in 42 primary ependymomas in two non-overlapping primary ependymoma cohorts, with the goal of identifying essential super-enhancer-associated genes on which tumour cells depend. Enhancer regions revealed putative oncogenes, molecular targets and pathways; inhibition of these targets with small molecule inhibitors or short hairpin RNA diminished the proliferation of patient-derived neurospheres and increased survival in mouse models of ependymomas. Through profiling of transcriptional enhancers, our study provides a framework for target and drug discovery in other cancers that lack known genetic drivers and are therefore difficult to treat., This work was supported by an Alex's Lemonade Stand Young Investigator Award (S.C.M.), The CIHR Banting Fellowship (S.C.M.), The Cancer Prevention Research Institute of Texas (S.C.M., RR170023), Sibylle Assmus Award for Neurooncology (K.W.P.), the DKFZ-MOST (Ministry of Science, Technology & Space, Israel) program in cancer research (H.W.), James S. McDonnell Foundation (J.N.R.) and NIH grants: CA154130 (J.N.R.), R01 CA169117 (J.N.R.), R01 CA171652 (J.N.R.), R01 NS087913 (J.N.R.) and R01 NS089272 (J.N.R.). R.C.G. is supported by NIH grants T32GM00725 and F30CA217065. M.D.T. is supported by The Garron Family Chair in Childhood Cancer Research, and grants from the Pediatric Brain Tumour Foundation, Grand Challenge Award from CureSearch for Children’s Cancer, the National Institutes of Health (R01CA148699, R01CA159859), The Terry Fox Research Institute and Brainchild. M.D.T. is also supported by a Stand Up To Cancer St. Baldrick’s Pediatric Dream Team Translational Research Grant (SU2C-AACR-DT1113).

Commonly-mutated genes have been found for many cancers, but less is known about mutations in cis-regulatory elements. We leverage gains in tumor-specific enhancer activity, coupled with allele-biased mutation detection from H3K27ac ChIP-seq data, to pinpoint potential enhancer-activating mutations in colorectal cancer (CRC). Analysis of a genetically-diverse cohort of CRC specimens revealed that microsatellite instable (MSI) samples have a high indel rate within active enhancers. Enhancers with indels show evidence of positive selection, increased target gene expression, and a subset is highly recurrent. The indels affect short homopolymer tracts of A/T and increase affinity for FOX transcription factors. We further demonstrate that signature mismatch-repair (MMR) mutations activate enhancers using a xenograft tumor metastasis model, where mutations are induced naturally via CRISPR/Cas9 inactivation of MLH1 prior to tumor cell injection. Our results suggest that MMR signature mutations activate enhancers in CRC tumor epigenomes to provide a selective advantage.

// Callie R. Merry 1, 2 , Sarah McMahon 1 , Megan E. Forrest 1 , Cynthia F. Bartels 1 , Alina Saiakhova 1 , Courtney A. Bartel 3 , Peter C. Scacheri 1, 3 , Cheryl L. Thompson 3, 5 , Mark W. Jackson 3 , Lyndsay N. Harris 3, 4 , Ahmad M. Khalil 1, 2, 3 1 Department of Genetics and Genome Sciences, Case Western Reserve University, Cleveland, OH 44106, USA 2 Department of Biochemistry, Case Western Reserve University, Cleveland, OH 44106, USA 3 Case Comprehensive Cancer Center, Case Western Reserve University, Cleveland, OH 44106, USA 4 Department of Medicine and Case Western Reserve University, Cleveland, OH 44106, USA 5 Department of Nutrition, Case Western Reserve University, Cleveland, OH 44106, USA Correspondence to: Ahmad M. Khalil, email: Dr.Ahmad.Khalil@gmail.com Keywords: breast cancer, HER2, drug resistance, trastuzumab-resistance, cancer therapy Received: March 11, 2016 Accepted: July 09, 2016 Published: July 16, 2016 ABSTRACT Approximately, 25–30% of early-stage breast tumors are classified at the molecular level as HER2-positive, which is an aggressive subtype of breast cancer. Amplification of the HER2 gene in these tumors results in a substantial increase in HER2 mRNA levels, and consequently, HER2 protein levels. HER2, a transmembrane receptor tyrosine kinase (RTK), is targeted therapeutically by a monoclonal antibody, trastuzumab (Tz), which has dramatically improved the prognosis of HER2-driven breast cancers. However, ~30% of patients develop resistance to trastuzumab and recur; and nearly all patients with advanced disease develop resistance over time and succumb to the disease. Mechanisms of trastuzumab resistance (TzR) are not well understood, although some studies suggest that growth factor signaling through other receptors may be responsible. However, these studies were based on cell culture models of the disease, and thus, it is not known which pathways are driving the resistance in vivo . Using an integrative transcriptomic approach of RNA isolated from trastuzumab-sensitive and trastuzumab-resistant HER2+ tumors, and isogenic cell culture models, we identified a small set of mRNAs and lincRNAs that are associated with trastuzumab-resistance (TzR). Functional analysis of a top candidate gene, S100P, demonstrated that inhibition of S100P results in reversing TzR. Mechanistically, S100P activates the RAS/MEK/MAPK pathway to compensate for HER2 inhibition by trastuzumab. Finally, we demonstrated that the upregulation of S100P appears to be driven by epigenomic changes at the enhancer level. Our current findings should pave the path toward new therapies for breast cancer patients.

Commonly-mutated genes have been found for many cancers, but less is known about mutations incis-regulatory elements. We leverage gains in tumor-specific enhancer activity, coupled with allele-biased mutation detection from H3K27ac ChIP-seq data, to pinpoint potential enhancer-activating mutations in colorectal cancer (CRC). Analysis of a genetically-diverse cohort of CRC specimens revealed that microsatellite instable (MSI) samples have a high indel rate within active enhancers. Enhancers with indels show evidence of positive selection, increased target gene expression, and a subset is highly recurrent. The indels affect short homopolymer tracts of A/T and increase affinity for FOX transcription factors. We further demonstrate that signature mismatch-repair (MMR) mutations activate enhancers using a xenograft tumor metastasis model, where mutations are induced naturally via CRISPR/Cas9 inactivation ofMLH1prior to tumor cell injection. Our results suggest that MMR signature mutations activate or augment enhancers in CRC tumor epigenomes to provide a selective advantage.

Long non-coding RNAs (lncRNAs) are frequently dysregulated in many human cancers. We sought to identify candidate oncogenic lncRNAs in human colon tumors by utilizing RNA sequencing data from 22 colon tumors and 22 adjacent normal colon samples from The Cancer Genome Atlas (TCGA). The analysis led to the identification of ~200 differentially expressed lncRNAs. Validation in an independent cohort of normal colon and patient-derived colon cancer cell lines identified a novel lncRNA, lincDUSP, as a potential candidate oncogene. Knockdown of lincDUSP in patient-derived colon tumor cell lines resulted in significantly decreased cell proliferation and clonogenic potential, and increased susceptibility to apoptosis. The knockdown of lincDUSP affects the expression of ~800 genes, and NCI pathway analysis showed enrichment of DNA damage response and cell cycle control pathways. Further, identification of lincDUSP chromatin occupancy sites by ChIRP-Seq demonstrated association with genes involved in the replication-associated DNA damage response and cell cycle control. Consistent with these findings, lincDUSP knockdown in colon tumor cell lines increased both the accumulation of cells in early S-phase and γH2AX foci formation, indicating increased DNA damage response induction. Taken together, these results demonstrate a key role of lincDUSP in the regulation of important pathways in colon cancer.

The degree of intrinsic and interpatient phenotypic heterogeneity and its role in tumour evolution is poorly understood. Phenotypic divergence can be achieved via the inheritance of alternative transcriptional programs1,2. Cell-type specific transcription is maintained through the activation of epigenetically-defined regulatory regions including promoters and enhancers1,3,4. In this work, we annotated the epigenome of 47 primary and metastatic oestrogen-receptor (ERα)-positive breast cancer specimens from clinical samples, and developed strategies to deduce phenotypic heterogeneity from the regulatory landscape, identifying key regulatory elements commonly shared across patients. Highly shared regions contain a unique set of regulatory information including the motif for the transcription factor YY1. In vitro work shows that YY1 is essential for ERα transcriptional activity and defines the critical subset of functional ERα binding sites driving tumor growth in most luminal patients. YY1 also control the expression of genes that mediate resistance to endocrine treatment. Finally, we show that H3K27ac levels at active enhancer elements can be used as a surrogate of intra-tumor phenotypic heterogeneity, and to track expansion and contraction of phenotypic subpopulations throughout breast cancer progression. Tracking YY1 and SLC9A3R1 positive clones in primary and metastatic lesions, we show that endocrine therapies drive the expansion of phenotypic clones originally underrepresented at diagnosis. Collectively, our data show that epigenetic mechanisms significantly contribute to phenotypic heterogeneity and evolution in systemically treated breast cancer patients.

Metastasis results from a complex set of traits acquired by tumor cells, distinct from those necessary for tumorigenesis. Here, we investigate the contribution of enhancer elements to the metastatic phenotype of osteosarcoma. Through epigenomic profiling, we identify substantial differences in enhancer activity between primary and metastatic tumors in human patients as well as nearisogenic pairs of high and low lung-metastatic osteosarcoma cells. We term these regions Metastatic Variant Enhancer Loci (Met-VELs). We demonstrate that these Met-VELs drive coordinated waves of gene expression during metastatic colonization of the lung. Met-VELs cluster non-randomly, indicating that activity of these enhancers and their associated gene targets are positively selected. As evidence of this causal association, osteosarcoma lung metastasis is inhibited by global interruptions of Met-VEL-associated gene expression via pharmacologic BET inhibition, by knockdown of AP-1 transcription factors that occupy Met-VELs, and by knockdown or functional inhibition of individual genes activated by Met-VELs, such as F3. We further show that genetic deletion of a single Met-VEL at the F3 locus blocks metastatic cell outgrowth in the lung. These findings indicate that Met-VELs and the genes they regulate play a functional role in metastasis and may be suitable targets for anti-metastatic therapies.

In addition to mutations in genes, aberrant enhancer element activity at non-coding regions of the genome is a key driver of tumorigenesis. Here, we perform epigenomic enhancer profiling of a cohort of more than forty genetically diverse human colorectal cancer (CRC) specimens. Using normal colonic crypt epithelium as a comparator, we identify enhancers with recurrently gained or lost activity across CRC specimens. Of the enhancers highly recurrently activated in CRC, most are constituents of super enhancers, are occupied by AP-1 and cohesin complex members, and originate from primed chromatin. Many activate known oncogenes, and CRC growth can be mitigated through pharmacologic inhibition or genome editing of these loci. Nearly half of all GWAS CRC risk loci co-localize to recurrently activated enhancers. These findings indicate that the CRC epigenome is defined by highly recurrent epigenetic alterations at enhancers which activate a common, aberrant transcriptional programme critical for CRC growth and survival., Active enhancers are defined by the presence of post-translational modifications of histones. Here, the authors use these marks to identify enhancers recurrently activated in colorectal cancer and find that these enhancers turn on oncogenes and are associated with known risk loci for developing the disease.

Summary Naive mouse embryonic stem cells (mESCs) and primed epiblast stem cells (mEpiSCs) represent successive snapshots of pluripotency during embryogenesis. Using transcriptomic and epigenomic mapping we show that a small fraction of transcripts are differentially expressed between mESCs and mEpiSCs and that these genes show expected changes in chromatin at their promoters and enhancers. Unexpectedly, the cis -regulatory circuitry of genes that are expressed at identical levels between these cell states also differs dramatically. In mESCs, these genes are associated with dominant proximal enhancers and dormant distal enhancers, which we term seed enhancers. In mEpiSCs, the naive-dominant enhancers are lost, and the seed enhancers take up primary transcriptional control. Seed enhancers have increased sequence conservation and show preferential usage in downstream somatic tissues, often expanding into super enhancers. We propose that seed enhancers ensure proper enhancer utilization and transcriptional fidelity as mammalian cells transition from naive pluripotency to a somatic regulatory program.

Trithorax (TRX) antagonizes epigenetic silencing by Polycomb group (PcG) proteins, stimulates enhancer-dependent transcription, and establishes a ‘cellular memory’ of active transcription of PcG-regulated genes. The mechanisms underlying these TRX functions remain largely unknown, but are presumed to involve its histone H3K4 methyltransferase activity. We report that the SET domains of TRX and TRX-related (TRR) have robust histone H3K4 monomethyltransferase activity in vitro and that Tyr3701 of TRX and Tyr2404 of TRR prevent them from being trimethyltransferases. The trxZ11 missense mutation (G3601S), which abolishes H3K4 methyltransferase activity in vitro, reduces the H3K4me1 but not the H3K4me3 level in vivo. trxZ11 also suppresses the impaired silencing phenotypes of the Pc3 mutant, suggesting that H3K4me1 is involved in antagonizing Polycomb silencing. Polycomb silencing is also antagonized by TRX-dependent H3K27 acetylation by CREB-binding protein (CBP). We show that perturbation of Polycomb silencing by TRX overexpression requires CBP. We also show that TRX and TRR are each physically associated with CBP in vivo, that TRX binds directly to the CBP KIX domain, and that the chromatin binding patterns of TRX and TRR are highly correlated with CBP and H3K4me1 genome-wide. In vitro acetylation of H3K27 by CBP is enhanced on K4me1-containing H3 substrates, and independently altering the H3K4me1 level in vivo, via the H3K4 demethylase LSD1, produces concordant changes in H3K27ac. These data indicate that the catalytic activities of TRX and CBP are physically coupled and suggest that both activities play roles in antagonizing Polycomb silencing, stimulating enhancer activity and cellular memory.

DNA variants (SNPs) that predispose to common traits often localize within noncoding regulatory elements such as enhancers. Moreover, loci identified by genome-wide association studies (GWAS) often contain multiple SNPs in linkage disequilibrium (LD), any of which may be causal. Thus, determining the effect of these multiple variant SNPs on target transcript levels has been a major challenge. Here, we provide evidence that for six common autoimmune disorders (rheumatoid arthritis, Crohn's disease, celiac disease, multiple sclerosis, lupus, and ulcerative colitis), the GWAS association arises from multiple polymorphisms in LD that map to clusters of enhancer elements active in the same cell type. This finding suggests a “multiple enhancer variant” hypothesis for common traits, where several variants in LD impact multiple enhancers and cooperatively affect gene expression. Using a novel method to delineate enhancer–gene interactions, we show that multiple enhancer variants within a given locus typically target the same gene. Using available data from HapMap and B lymphoblasts as a model system, we provide evidence at numerous loci that multiple enhancer variants cooperatively contribute to altered expression of their gene targets. The effects on target transcript levels tend to be modest and can be either gain- or loss-of-function. Additionally, the genes associated with multiple enhancer variants encode proteins that are often functionally related and enriched in common pathways. Overall, the multiple enhancer variant hypothesis offers a new paradigm by which noncoding variants can confer susceptibility to common traits.

The transcription of ribosomal RNA (rRNA) is critical to life. Despite its importance, ribosomal DNA (rDNA) is not included in current genome assemblies and, consequently, genomic analyses to date have excluded rDNA. Here, we show that short sequence reads can be aligned to a genome assembly containing a single rDNA repeat. Integrated analysis of ChIP-seq, DNase-seq, MNase-seq and RNA-seq data reveals several novel findings. First, the coding region of active rDNA is contained within nucleosome-depleted open chromatin that is highly transcriptionally active. Second, histone modifications are located not only at the rDNA promoter but also at novel sites within the intergenic spacer. Third, the distributions of active modifications are more similar within and between different cell types than repressive modifications. Fourth, UBF, a positive regulator of rRNA transcription, binds to sites throughout the genome. Lastly, the insulator binding protein CTCF associates with the spacer promoter of rDNA, suggesting that transcriptional insulation plays a role in regulating the transcription of rRNA. Taken together, these analyses confirm and expand the results of previous ChIP studies of rDNA and provide novel avenues for exploration of chromatin-mediated regulation of rDNA.

Metastasis results from a complex set of traits acquired by tumor cells, distinct from those necessary for tumorigenesis. Here, we investigate the contribution of enhancer elements to the metastatic phenotype of osteosarcoma. Through epigenomic profiling, we identify substantial differences in signature enhancer-histone marks between near-isogenic pairs of high and low lung-metastatic osteosarcoma cells. We term these regions Metastatic Variant Enhancer Loci (Met-VELs). Met-VELs drive coordinated waves of gene expression during metastatic colonization of the lung. Met-VELs cluster non-randomly, indicating that activity of these enhancers and their associated gene targets is positively selected. Osteosarcoma lung metastasis is inhibited by global interruptions of Met-VEL associated gene expression via pharmacologic BET inhibition, by knockdown of AP-1 transcription factors whose motifs are enriched in Met-VELs, and by knockdown of individual genes activated by Met-VELs. These observations have implications for the discovery and development of targeted anti-metastatic therapies. Citation Format: James J. Morrow, Tyler E. Miller, Alina Saiakhova, Michael M. Lizardo, Cynthia F. Bartels, Ian Bayles, Stevephen Hung, Arnulfo Mendoza, Jay T. Myers, Frederick Allen, Analisa DiFeo, Brian P. Rubin, Alex Y. Huang, Paul S. Meltzer, Lee J. Helman, Chand Khanna, Peter C. Scacheri. Positively selected enhancer elements endow tumor cells with metastatic competence. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr LB-151.

Colorectal Cancer Signature The mutations and genome aberrations that characterize cancer result in often dramatically altered gene and protein expression patterns. It is these altered expression patterns that directly and indirectly drive progression of the disease. In human primary colorectal cancer cells, Akhtar-Zaidi et al. (p. 736 , published online 12 April) analyzed the pattern of epigenetically modified chromatin at “enhancer” sequences that are known to be critical in the control of gene expression. An epigenetic enhancer signature was defined that was specifically associated with colorectal cancer cells.