Epigenomic Comparison Reveals Activation of “Seed” Enhancers during Transition from Naive to Primed Pluripotency

SummaryNaive 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

Genome-Wide Analysis of Menin Binding Provides Insights into MEN1 Tumorigenesis

Multiple endocrine neoplasia type I (MEN1) is a familial cancer syndrome characterized primarily by tumors of multiple endocrine glands. The gene for MEN1 encodes a ubiquitously expressed tumor suppressor protein called menin. Menin was recently shown to interact with several components of a trithorax family histone methyltransferase complex including ASH2, Rbbp5, WDR5, and the leukemia proto-oncoprotein MLL. To elucidate menin's role as a tumor suppressor and gain insights into the endocrine-specific tumor phenotype in MEN1, we mapped the genomic binding sites of menin, MLL1, and Rbbp5, to approximately 20,000 promoters in HeLa S3, HepG2, and pancreatic islet cells using the strategy of chromatin-immunoprecipitation coupled with microarray analysis. We found that menin, MLL1, and Rbbp5 localize to the promoters of thousands of human genes but do not always bind together. These data suggest that menin functions as a general regulator of transcription. We also found that factor occupancy generally correlates with high gene expression but that the loss of menin does not result in significant changes in most transcript levels. One exception is the developmentally programmed transcription factor, HLXB9, which is overexpressed in islets in the absence of menin. Our findings expand the realm of menin-targeted genes several hundred-fold beyond that previously described and provide potential insights to the endocrine tumor bias observed in MEN1 patients

Integrative genomic analysis of human ribosomal DNA

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

CHD7 Targets Active Gene Enhancer Elements to Modulate ES Cell-Specific Gene Expression

CHD7 is one of nine members of the chromodomain helicase DNA–binding domain family of ATP–dependent chromatin remodeling enzymes found in mammalian cells. De novo mutation of CHD7 is a major cause of CHARGE syndrome, a genetic condition characterized by multiple congenital anomalies. To gain insights to the function of CHD7, we used the technique of chromatin immunoprecipitation followed by massively parallel DNA sequencing (ChIP–Seq) to map CHD7 sites in mouse ES cells. We identified 10,483 sites on chromatin bound by CHD7 at high confidence. Most of the CHD7 sites show features of gene enhancer elements. Specifically, CHD7 sites are predominantly located distal to transcription start sites, contain high levels of H3K4 mono-methylation, found within open chromatin that is hypersensitive to DNase I digestion, and correlate with ES cell-specific gene expression. Moreover, CHD7 co-localizes with P300, a known enhancer-binding protein and strong predictor of enhancer activity. Correlations with 18 other factors mapped by ChIP–seq in mouse ES cells indicate that CHD7 also co-localizes with ES cell master regulators OCT4, SOX2, and NANOG. Correlations between CHD7 sites and global gene expression profiles obtained from Chd7+/+, Chd7+/−, and Chd7−/− ES cells indicate that CHD7 functions at enhancers as a transcriptional rheostat to modulate, or fine-tune the expression levels of ES–specific genes. CHD7 can modulate genes in either the positive or negative direction, although negative regulation appears to be the more direct effect of CHD7 binding. These data indicate that enhancer-binding proteins can limit gene expression and are not necessarily co-activators. Although ES cells are not likely to be affected in CHARGE syndrome, we propose that enhancer-mediated gene dysregulation contributes to disease pathogenesis and that the critical CHD7 target genes may be subject to positive or negative regulation

Therapeutic targeting of ependymoma as informed by oncogenic enhancer profiling

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)

Myc-regulated microRNAs attenuate embryonic stem cell differentiation

Myc proteins are known to have an important function in stem cell maintenance. As Myc has been shown earlier to regulate microRNAs (miRNAs) involved in proliferation, we sought to determine whether c-Myc also affects embryonic stem (ES) cell maintenance and differentiation through miRNAs. Using a quantitative primer-extension PCR assay we identified miRNAs, including, miR-141, miR-200, and miR-429 whose expression is regulated by c-Myc in ES cells, but not in the differentiated and tumourigenic derivatives of ES cells. Chromatin immunoprecipitation analyses indicate that in ES cells c-Myc binds proximal to genomic regions encoding the induced miRNAs. We used expression profiling and seed homology to identify genes specifically downregulated both by these miRNAs and by c-Myc. We further show that the introduction of c-Myc-induced miRNAs into murine ES cells significantly attenuates the downregulation of pluripotency markers on induction of differentiation after withdrawal of the ES cell maintenance factor LIF. In contrast, knockdown of the endogenous miRNAs accelerate differentiation. Our data show that in ES cells c-Myc acts, in part, through a subset of miRNAs to attenuate differentiation

Otitis Media in a New Mouse Model for CHARGE Syndrome with a Deletion in the Chd7 Gene

Otitis media is a middle ear disease common in children under three years old. Otitis media can occur in normal individuals with no other symptoms or syndromes, but it is often seen in individuals clinically diagnosed with genetic diseases such as CHARGE syndrome, a complex genetic disease caused by mutation in the Chd7 gene and characterized by multiple birth defects. Although otitis media is common in human CHARGE syndrome patients, it has not been reported in mouse models of CHARGE syndrome. In this study, we report a mouse model with a spontaneous deletion mutation in the Chd7 gene and with chronic otitis media of early onset age accompanied by hearing loss. These mice also exhibit morphological alteration in the Eustachian tubes, dysregulation of epithelial proliferation, and decreased density of middle ear cilia. Gene expression profiling revealed up-regulation of Muc5ac, Muc5b and Tgf-β1 transcripts, the products of which are involved in mucin production and TGF pathway regulation. This is the first mouse model of CHARGE syndrome reported to show otitis media with effusion and it will be valuable for studying the etiology of otitis media and other symptoms in CHARGE syndrome

Identification of CIITA Regulated Genetic Module Dedicated for Antigen Presentation

The class II trans-activator CIITA is a transcriptional co-activator required for the expression of Major Histocompatibility Complex (MHC) genes. Although the latter function is well established, the global target-gene specificity of CIITA had not been defined. We therefore generated a comprehensive list of its target genes by performing genome-wide scans employing four different approaches designed to identify promoters that are occupied by CIITA in two key antigen presenting cells, B cells and dendritic cells. Surprisingly, in addition to MHC genes, only nine new targets were identified and validated by extensive functional and expression analysis. Seven of these genes are known or likely to function in processes contributing to MHC-mediated antigen presentation. The remaining two are of unknown function. CIITA is thus uniquely dedicated for genes implicated in antigen presentation. The finding that CIITA regulates such a highly focused gene expression module sets it apart from all other transcription factors, for which large-scale binding-site mapping has indicated that they exert pleiotropic functions and regulate large numbers of genes

A Genetic Locus within the FMN1/GREM1 Gene Region Interacts with Body Mass Index in Colorectal Cancer Risk

Colorectal cancer risk can be impacted by genetic, environmental, and lifestyle factors, including diet and obesity. Geneenvironment interactions (G x E) can provide biological insights into the effects of obesity on colorectal cancer risk. Here, we assessed potential genome-wide G x E interactions between body mass index (BMI) and common SNPs for colorectal cancer risk using data from 36,415 colorectal cancer cases and 48,451 controls from three international colorectal cancer consortia (CCFR, CORECT, and GECCO). The G x E tests included the conventional logistic regression using multiplicative terms (one degree of freedom, 1DF test), the two-step EDGE method, and the joint 3DF test, each of which is powerful for detecting G x E interactions under specific conditions. BMI was associated with higher colorectal cancer risk. The two-step approach revealed a statistically significant GxBMI interaction located within the Formin 1/Gremlin 1 (FMN1/GREM1) gene region (rs58349661). This SNP was also identified by the 3DF test, with a suggestive statistical significance in the 1DF test. Among participants with the CC genotype of rs58349661, overweight and obesity categories were associated with higher colorectal cancer risk, whereas null associations were observed across BMI categories in those with the TT genotype. Using data from three large international consortia, this study discovered a locus in the FMN1/GREM1 gene region that interacts with BMI on the association with colorectal cancer risk. Further studies should examine the potential mechanisms through which this locus modifies the etiologic link between obesity and colorectal cancer

Discovery of common and rare genetic risk variants for colorectal cancer.

To further dissect the genetic architecture of colorectal cancer (CRC), we performed whole-genome sequencing of 1,439 cases and 720 controls, imputed discovered sequence variants and Haplotype Reference Consortium panel variants into genome-wide association study data, and tested for association in 34,869 cases and 29,051 controls. Findings were followed up in an additional 23,262 cases and 38,296 controls. We discovered a strongly protective 0.3% frequency variant signal at CHD1. In a combined meta-analysis of 125,478 individuals, we identified 40 new independent signals at P < 5 × 10-8, bringing the number of known independent signals for CRC to ~100. New signals implicate lower-frequency variants, Krüppel-like factors, Hedgehog signaling, Hippo-YAP signaling, long noncoding RNAs and somatic drivers, and support a role for immune function. Heritability analyses suggest that CRC risk is highly polygenic, and larger, more comprehensive studies enabling rare variant analysis will improve understanding of biology underlying this risk and influence personalized screening strategies and drug development.Goncalo R Abecasis has received compensation from 23andMe and Helix. He is currently an employee of Regeneron Pharmaceuticals. Heather Hampel performs collaborative research with Ambry Genetics, InVitae Genetics, and Myriad Genetic Laboratories, Inc., is on the scientific advisory board for InVitae Genetics and Genome Medical, and has stock in Genome Medical. Rachel Pearlman has participated in collaborative funded research with Myriad Genetics Laboratories and Invitae Genetics but has no financial competitive interest