Coordinated Repression of Cell Cycle Genes by KDM5A and E2F4 During Differentiation
Total Page:16
File Type:pdf, Size:1020Kb
Coordinated repression of cell cycle genes by KDM5A and E2F4 during differentiation Michael L. Beshiria,1, Katherine B. Holmesa,1, William F. Richtera, Samuel Hessa, Abul B. M. M. K. Islama,b, Qin Yanc, Lydia Planted, Larisa Litovchicke, Nicolas Gévryd, Nuria Lopez-Bigasb,f, William G. Kaelin, Jr.e,2, and Elizaveta V. Benevolenskayaa,2 aDepartment of Biochemistry and Molecular Genetics, University of Illinois, Chicago, IL 60607; bResearch Unit on Biomedical Informatics, Department of Experimental and Health Sciences, Universitat Pompeu Fabra, 08003 Barcelona, Spain; cDepartment of Pathology, Yale University School of Medicine, New Haven, CT 06520; dDépartement de biologie, Faculté des sciences, Université de Sherbrooke, Sherbrooke, QC, Canada J1K 2R1; eDepartment of Medical Oncology, Dana–Farber Cancer Institute and Brigham and Women’s Hospital, Harvard Medical School, Boston, MA 02115; and fInstitució Catalana de Recerca i Estudis Avançats (ICREA), Passeig Lluis Companys, 08010 Barcelona, Spain Contributed by William G. Kaelin, Jr., September 27, 2012 (sent for review August 10, 2012) Epigenetic regulation underlies the robust changes in gene expres- role in human tumors, most of which carry mutations in the RB sion that occur during development. How precisely epigenetic pathway. During cell cycle progression, pRB and the related pro- enzymes contribute to development and differentiation processes teins p107 and p130 (collectively called “pocket” proteins) are pe- is largely unclear. Here we show that one of the enzymes that riodically and reversibly recruited to E2F target genes, accompanied removes the activating epigenetic mark of trimethylated lysine 4 on by recruitment of histone deacetylase Sin3/HDAC and chromatin- histone H3, lysine (K)-specific demethylase 5A (KDM5A), reinforces remodeling activities (2). In the Go/quiescent state, the repression the effects of the retinoblastoma (RB) family of transcriptional of cell cycle genes is mediated by the DP, RB, E2F, and MuvB (DREAM) complex that contains p130 and E2F4 (3). repressors on differentiation. Global location analysis showed that Here we show that the histone demethylase KDM5A is KDM5A cooccupies a substantial portion of target genes with the recruited to genes experiencing strong repression during the E2F4 transcription factor. During ES cell differentiation, knockout of course of differentiation, and that these genes are also occupied KDM5A resulted in derepression of multiple genomic loci that are by components of the DREAM complex. Using high-resolution targets of KDM5A, denoting a direct regulatory function. In termi- ChIP-seq analysis, we found that KDM5A preferentially binds to nally differentiated cells, common KDM5A and E2F4 gene targets the TSS regions. Using cells where KDM5A expression level was were bound by the pRB-related protein p130, a DREAM complex decreased either by RNA inhibition or by using null and condi- component. KDM5A was recruited to the transcription start site tional (floxed) KDM5A alleles, we show that KDM5A is critical for regions independently of E2F4; however, it cooperated with E2F4 H3K4 demethylation and repression of cell cycle genes. Differ- to promote a state of deepened repression at cell cycle genes during entiation time-course experiments showed a prominent recruit- differentiation. These findings reveal a critical role of H3K4 deme- ment of both KDM5A and E2F4 during later stages to specific thylation by KDM5A in the transcriptional silencing of genes that are genomic loci. Importantly, the recruitment of KDM5A and E2F4 suppressed by RB family members in differentiated cells. was independent of one another. Specifically, E2F4 knockdown did not abolish KDM5 recruitment and H3K4 demethylation, and chromatin | histone demethylase | whole-genome sequencing | KDM5A knockdown did not affect E2F4 recruitment and histone histone methylation acetylation. Nevertheless, KDM5A and E2F4 cooperated at cell cycle genes to enforce transcriptional repression. egulation of gene expression is accomplished by transcription Results Rfactors, histone-modifying enzymes, and chromatin remodeling fi machinery. The combined action of these three has been implicated ChIP-Seq Analysis for KDM5A Identi es Common Targets with E2F4. We previously showed that mouse ES cells lacking KDM5A are in a number of biological processes, including cell cycle control, fi development, reprogramming, differentiation, and aging. Deregu- signi cantly impaired with respect to differentiation (4). To de- lation of chromatin-modifying enzymes has been strongly linked to termine KDM5A targets genome-wide, we performed ChIP-seq the development of cancer. For example, two enzymes that regulate experiments in two different clones of mouse ES cells with a con- ditional (floxed) KDM5A allele (Kdm5af/f) and a clone where the methylation at the histone H3 lysine 4 (H3K4) residue, mixed lin- −/− eage leukemia-1 (MLL1) and lysine (K)-specific demethylase 5A KDM5A allele was deleted (Kdm5a ). As expected, we were un- fi able to detect a significant enrichment of bound DNA in ChIP-seq (KDM5A), have been identi ed in translocations associated with −/− fi fi human leukemia. H3K4 histone methylation states exhibit a highly samples from Kdm5a cells, con rming high speci city of the fi KDM5A antibody. We detected 3,093 regions bound by KDM5A distinct distribution pattern in the genome. Speci cally, H3K4 tri- f/f methylation (H3K4me3) is strongly associated with transcriptional in Kdm5a cells (Datasets S1 and S2). Relating KDM5A-bound activation, with the highest levels observed near transcriptional regions to known mouse genes showed that KDM5A is preferen- start sites (TSS). In vitro studies suggest that the four KDM5 tially bound at TSS regions, especially among the top 5% most- enzymes (KDM5A, KDM5B, KDM5C, and KDM5D) are able to highly expressed genes (Fig. 1A). To determine whether direct remove methylation at lysine 4 of histone H3, and in vivo KDM5A and KDM5B may be recruited to common gene regions (1). This finding leads to multiple questions: (i) Do KDM5 proteins play Author contributions: W.G.K. and E.V.B. designed research; M.L.B., K.B.H., W.F.R., S.H., a role in gene regulation by transcription factors? (ii) Are KDM5 and E.V.B. performed research; M.L.B., K.B.H., Q.Y., L.P., N.G., W.G.K., and E.V.B. contrib- GENETICS enzymes nonredundant H3K4 demethylases? uted new reagents/analytic tools; M.L.B., K.B.H., A.B.M.M.K.I., L.L., N.L.-B., and E.V.B. During differentiation, cells exhibit two novel properties: re- analyzed data; and M.L.B., W.G.K., and E.V.B. wrote the paper. pression of cell cycle genes associated with permanent cell cycle The authors declare no conflict of interest. exit and activation of cell type-specific genes. Histone mod- Data deposition: The data reported in this paper have been deposited in the Gene Ex- ifications are thought to be important epigenetic events intimately pression Omnibus (GEO) database, www.ncbi.nlm.nih.gov/geo (accession no. GSE28343). linked to initiation and maintenance of transcriptional changes for 1M.L.B. and K.B.H. contributed equally to this work. both of these processes. Cell cycle exit is associated with repression 2To whom correspondence may be addressed. E-mail: [email protected] or of a large group of genes, many of which are targets of the tran- [email protected]. scription factor complex between the retinoblastoma protein (pRB) This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10. and E2F. Derepression of pRB/E2F target genes plays an important 1073/pnas.1216724109/-/DCSupplemental. www.pnas.org/cgi/doi/10.1073/pnas.1216724109 PNAS | November 6, 2012 | vol. 109 | no. 45 | 18499–18504 Downloaded by guest on October 1, 2021 Fig. 1. KDM5A ChIP-seq analysis in mouse ES cells reveals corecruitment of KDM5A and E2F4. (A) Av- erage KDM5A enrichment near TSS and in genic region. Genic regions are presented as a 3-kb-long metagene surrounded by a 1-kb region upstream TSS and 1-kb downstream region from transcription termination site (TTS). Dotted lines show TSS and TTS. Most KDM5A binding occurs at the TSS. High TSS ChIP signal is observed for KDM5A reads in the top 5% expressed genes compared with the bottom 5% expressed genes (expression data from GEO GSE5914). A drop in the signal occurs near the TTS, probably because of nucleosome depletion. (B) − − KDM5A targets in Kdm5a / ES cells are differen- tially expressed. Gene expression microarray data were generated from ES cells that were induced to differentiate by leukemia inhibitory factor (LIF) withdrawal (4). Change in expression for each gene − − was determined in Kdm5a / (KO) versus normal Kdm5af/f (FF) cells (Dataset S3). Box-and-whisker plots show distribution of Log2 fold-change values for genes associated with the classified KDM5A peaks (in relation to the closest gene as distal, promoter, TSS, gene body, or downstream regions). Solid bars of boxes display the 25–75% of ranked genes, with the median indicated as an intersection. The yellow bars depicting the expression changes that occur in differentiated cells in the group with KDM5A bound at the TSS are shifted toward higher values, indicating activation of this group of genes, which is statistically significant as determined by Wilcoxon rank sum test (Dataset S4). (C) KDM5A and E2F4 bind to many common genes. The Venn diagram shows the overlap between KDM5A tar- gets (present study) and E2F4 targets (5). (D) En- richment analysis for common and distinct target genes shows relations to GO terms and KEGG pathways. Corrected (false-discovery rate) P values are delineated in a colored heatmap, where red signifies overrepresentation of targets in a particu- lar term. All large GO and KEGG groups with over- representation are shown. KDM5A targets are changed in expression in cells lacking KDM5A, “overlap” in Fig. 1D). KDM5A also bound multiple genes that we further analyzed previously acquired gene-expression data (4). were not direct targets of E2F4, yet these genes still participated in We found differences in expression of multiple KDM5A-bound functions linked to genes regulated by both factors.