G9a/GLP-dependent H3K9me2 patterning alters chromatin structure at CpG islands in hematopoietic progenitors
© Schones et al.; licensee BioMed Central Ltd. 2014
Received: 21 June 2014
Accepted: 1 September 2014
Published: 10 September 2014
The formation of chromatin domains is an important step in lineage commitment. In human hematopoietic stem and progenitor cells (HSPCs), G9a/GLP-dependent H3K9me2 chromatin territories form de novo during lineage specification and are nucleated at punctate sites during lineage commitment. Here, we examined the patterning of G9a/GLP-dependent H3K9me2 in HSPCs and the consequences for chromatin structure.
We profiled chromatin accessibility across the genome of HSPCs treated with either a small molecule inhibitor of G9a/GLP or DMSO. We observed that chromatin accessibility is dramatically altered at the regions of H3K9me2 nucleation. We have characterized the regions of H3K9me2 nucleation, with our analysis revealing that H3K9me2 is nucleated in HSPCs at CpG islands (CGIs) and CGI-like sequences across the genome. Our analysis furthermore revealed a bias of H3K9me2 nucleation towards regions with low rates of C- > T deamination, which typically lack DNA methylation. Lastly, we examined the interaction of H3K9me2 and DNA methylation and determined that chromatin accessibility changes upon loss of H3K9me2 are dependent on the presence of DNA methylation.
These results indicate that H3K9me2 nucleation is established at specific sequences that have base composition similar to CGIs. Our results furthermore indicate that H3K9me2 nucleation leads to local changes in chromatin accessibility and that H3K9me2 and DNA methylation work synergistically to regulate chromatin accessibility.
KeywordsChromatin accessibility CpG island G9a GLP H3K9me2 Hematopoietic progenitor
G9a/EHMT2 and GLP/EHMT1 are conserved protein lysine methyltransferases that play key roles in regulating gene expression and chromosome structure during mammalian development through de novo mono- and di-methylation of histone H3 lysine 9 (H3K9me1/2), histone marks associated with transcriptional repression [1–5]. During embryogenesis, large G9a/GLP-dependent H3K9me2 chromatin territories arise that have been proposed to reinforce lineage choice by determining higher order chromatin structure .
We recently observed that in adult human hematopoietic stem and progenitor cells (HSPCs), H3K9me2 chromatin territories are absent in primitive cells and are formed de novo during lineage commitment . In committed HSPCs, G9a/GLP activity nucleates H3K9me2 marks at CpG islands (CGIs) and other genomic sites, and this mark then spreads to form larger domains during differentiation . A recently developed small molecule inhibitor of G9a and GLP, UNC0638, inhibits the methyltransferase activity of both proteins by blocking substrate access to the SET domains . We have shown that treatment of HSPCs with UNC0638 results in a genome-wide loss of H3K9me2, a less dramatic reduction in H3K9me1 and no effect on H3K9me3 or the expression of G9a . These results are consistent with previous studies showing that loss of G9a leads to loss of H3K9me1/me2 [8, 9]. We furthermore observed that HSPCs treated with UNC0638, a G9a/GLP small molecular inhibitor, better retain stem cell-like phenotypes and function during in vitro expansion and increased expression of lineage-affiliated genes and certain gene clusters, suggestive of changes in regulation of chromatin structure .
Primitive hematopoietic stem cells (HSCs) have been hypothesized to have a more “open” chromatin structure that might help maintain a multipotent state by, for example, allowing transcriptional priming of lineage-affiliated genes [10–12]. One possible interpretation from our previous data is that G9a/GLP-H3K9me2 patterning helps restrict chromatin accessibility to reinforce lineage commitment. To investigate this, we examined the consequences of G9a/GLP-dependent H3K9me2 patterning on chromatin structure in HSPCs using FAIRE-seq (Formaldehyde Assisted Isolation of Regulatory Elements Sequencing)  to map accessible chromatin in CD34+ HSPCs treated with UNC0638 or dimethyl sulfoxide (DMSO) control (see Methods and Additional file 1: Table S1). We furthermore investigated the sequence features of sites of H3K9me2 nucleation. Our results indicate that H3K9me2 is nucleated at CGI-like sites across the genome, with a bias towards regions with low rates of C- > T deamination. Our results further demonstrate that H3K9me2 nucleation is associated with loss of chromatin accessibility and that changes in chromatin accessibility corresponding to loss of H3K9me2 are dependent on the presence of DNA methylation.
Results and discussion
To examine this on a genome-wide scale, we identified all sites of H3K9me2 nucleation (“peaks”) in CD34+ cells and counted the FAIRE reads in DMSO- and UNC0638-treated cells. This analysis revealed increased chromatin accessibility as the predominant behavior at H3K9me2 sites upon treatment with UNC0638 (Figure 1D). To evaluate these results in terms of the background level of chromatin changes across the genome, we randomly sampled sites from the genome and calculated the fold change in FAIRE signal at these randomly chosen regions. Compared to regions of H3K9me2 nucleation, randomly sampled regions had significantly smaller changes (P = 2.03577 × 10-250; Wilcoxon rank sum test) in chromatin accessibility (Figure 1E). We further evaluated the background changes in chromatin accessibility by sliding 1 kb windows in 50 bp increments across the genome and calculating the fold change of FAIRE read density in UNC0638-treated cells versus control cells. This analysis revealed that changes in FAIRE-seq read densities upon UNC0638 treatment were largely specific to H3K9me2 nucleation sites, indicating that changes in chromatin structure are specific to sites of H3K9me2 nucleation (see Additional file 3: Figure S2).
We next stratified CGIs based on whether they are promoter associated or “orphan” sites, as defined by Illingworth et al.  and analyzed UNC0638-driven changes in chromatin. Of the CGIs in the human genome, approximately half are associated with promoters and may play roles in facilitating transcriptional regulation while the other half are found in inter- and intragenic regions and have unknown functions (so-called orphan CGIs) . Analysis of FAIRE read counts across both promoter and orphan CGIs indicates that both sets of CGIs display similar behavior to all sites of H3K9me2 nucleation, with dramatic increases in chromatin accessibility in response to UNC0638 (Figure 2C).
In summary, we find that H3K9me2 patterning regulates chromatin structure at promoter and orphan CGIs and other sites of H3K9me2 nucleation, specifically promoting “closed” chromatin states. These results support the notion that G9a/GLP-H3K9me2 participates in global changes in chromatin structure in addition to histone patterning during HSC lineage formation. However, the biological significance of this patterning remains a question. One possibility is that H3K9me2 patterning helps reinforce chromatin states at sites of transcription during lineage specification, which may need to be reset in certain lineages and re-formed de novo. To our knowledge, this is the first observation of coordination between H3K9me2 patterning, promoter and orphan CGIs, DNA methylation, and chromatin structure.
Cell culture and treatment
Human CD34+ cells from healthy adults were purchased from the Fred Hutchinson Cancer Research Center Cell Processing Shared Resource, as described previously . Unfractionated CD34+ cells were treated with 2 μM of UNC0638 or 0.02% DMSO for 48 h, as described previously .
FAIRE was performed as previously described . Paired-end sequencing (100 × 100) was performed in replicate on a HiSeq 2500 to obtain ~50 M reads per replicate. Sequenced reads were aligned to the hg19 build of the human genome (hg19; GRCh37) using bowtie2  with local read alignment. Aligned reads were further filtered to exclude improperly paired reads and duplicate reads. Wiggle tracks were prepared for visualization on the UCSC Genome Browser  by sliding 10 bp windows across each chromosome and counting the sequenced fragments overlapping each window; reproducibility of FAIRE tracks was assessed visually and replicate libraries were combined to make final bed files for each condition. Peaks of FAIRE-seq were called with F-seq  using default parameters and a 200 bp feature length. Irreproducible Discover Rate analysis  was performed to identify reproducible peaks.
Aligned bam files for H3K9me2 ChIP-seq data were obtained from . Wiggle tracks were generated for visualization on the UCSC Genome Browser . Visual examination of CD34 HSPC H3K9me2 peaks indicated punctate peaks and regions of enrichment were identified using MACS with a P value threshold of 1 × 10-10. This analysis resulted in 43,159 peaks.
DNA methylation results were obtained from . Visualization of DNA methylation levels in CD90+ HSCs and CD34+ HSPCs (see Additional file 5: Figure S4) revealed a bimodal distribution. CGIs were considered methylated with a methylation score >0.75 and unmethylated with a methylation score <0.25. All regions in between were considered indeterminate.
Simulation of random sites for Figure 2E was performed by randomly choosing 1,000 regions of 1 kb 10,000 times and calculating the fold change of FAIRE signal for UNC0638 over DMSO at each region.
The data discussed in this publication have been deposited in NCBI’s Gene Expression Omnibus  and are accessible through GEO Series accession number GSE59749.
Formaldehyde Assisted Isolation of Regulatory Elements
Hematopoietic stem cell
Hematopoietic stem and progenitor cells
Histone H3 Lysine 9 di-methylation.
This work was supported by K22HL101950 (DS), U01HL099993 (PP), U01HL099997 (PP), American Cancer Society Research Scholar Grant (ACS RSG-14-056-01-LIB) (PP) and the HHMI/UW Molecular Medicine Scholar award (XC). Research reported in this publication included work performed in the Integrative Genomics Core of the City of Hope supported by the National Cancer Institute of the National Institutes of Health under award number P30CA33572 and work support by the Core Center of Excellence in Hematology (FHCRC) under award number P30 DK56465-14. The authors would like to thank Arthur Riggs, Amy Leung, Beverly Torok-Storb, Matthew Fero, David Emery, Mark Groudine, MA Bender, and the Schones and Paddison labs for helpful discussions.
- Litt MD, Simpson M, Gaszner M, Allis CD, Felsenfeld G: Correlation between histone lysine methylation and developmental changes at the chicken beta-globin locus. Science. 2001, 293: 2453-2455.View ArticlePubMedGoogle Scholar
- Noma K, Allis CD, Grewal SI: Transitions in distinct histone H3 methylation patterns at the heterochromatin domain boundaries. Science. 2001, 293: 1150-1155.View ArticlePubMedGoogle Scholar
- Su RC, Brown KE, Saaber S, Fisher AG, Merkenschlager M, Smale ST: Dynamic assembly of silent chromatin during thymocyte maturation. Nat Genet. 2004, 36: 502-506.View ArticlePubMedGoogle Scholar
- Wen B, Wu H, Shinkai Y, Irizarry RA, Feinberg AP: Large histone H3 lysine 9 dimethylated chromatin blocks distinguish differentiated from embryonic stem cells. Nat Genet. 2009, 41: 246-250.PubMed CentralView ArticlePubMedGoogle Scholar
- Shankar SR, Bahirvani AG, Rao VK, Bharathy N, Ow JR, Taneja R: G9a, a multipotent regulator of gene expression. Epigenetics. 2013, 8: 16-22.PubMed CentralView ArticlePubMedGoogle Scholar
- Chen X, Skutt-Kakaria K, Davison J, Ou YL, Choi E, Malik P, Loeb K, Wood B, Georges G, Torok-Storb B, Paddison PJ: G9a/GLP-dependent histone H3K9me2 patterning during human hematopoietic stem cell lineage commitment. Genes Dev. 2012, 26: 2499-2511.PubMed CentralView ArticlePubMedGoogle Scholar
- Vedadi M, Barsyte-Lovejoy D, Liu F, Rival-Gervier S, Allali-Hassani A, Labrie V, Wigle TJ, Dimaggio PA, Wasney GA, Siarheyeva A, Dong A, Tempel W, Wang SC, Chen X, Chau I, Mangano TJ, Huang XP, Simpson CD, Pattenden SG, Norris JL, Kireev DB, Tripathy A, Edwards A, Roth BL, Janzen WP, Garcia BA, Petronis A, Ellis J, Brown PJ, Frye SV, Arrowsmith CH, Jin J: A chemical probe selectively inhibits G9a and GLP methyltransferase activity in cells. Nat Chem Biol. 2011, 7: 566-574.PubMed CentralView ArticlePubMedGoogle Scholar
- Tachibana M, Sugimoto K, Nozaki M, Ueda J, Ohta T, Ohki M, Fukuda M, Takeda N, Niida H, Kato H, Shinkai Y: G9a histone methyltransferase plays a dominant role in euchromatic histone H3 lysine 9 methylation and is essential for early embryogenesis. Genes Dev. 2002, 16: 1779-1791.PubMed CentralView ArticlePubMedGoogle Scholar
- Rice JC, Briggs SD, Ueberheide B, Barber CM, Shabanowitz J, Hunt DF, Shinkai Y, Allis CD: Histone methyltransferases direct different degrees of methylation to define distinct chromatin domains. Mol Cell. 2003, 12: 1591-1598.View ArticlePubMedGoogle Scholar
- Hu M, Krause D, Greaves M, Sharkis S, Dexter M, Heyworth C, Enver T: Multilineage gene expression precedes commitment in the hemopoietic system. Genes Dev. 1997, 11: 774-785.View ArticlePubMedGoogle Scholar
- Mansson R, Hultquist A, Luc S, Yang L, Anderson K, Kharazi S, Al-Hashmi S, Liuba K, Thoren L, Adolfsson J, Buza-Vidas N, Qian H, Soneji S, Enver T, Sigvardsson M, Jacobsen SE: Molecular evidence for hierarchical transcriptional lineage priming in fetal and adult stem cells and multipotent progenitors. Immunity. 2007, 26: 407-419.View ArticlePubMedGoogle Scholar
- Miyamoto T, Iwasaki H, Reizis B, Ye M, Graf T, Weissman IL, Akashi K: Myeloid or lymphoid promiscuity as a critical step in hematopoietic lineage commitment. Dev Cell. 2002, 3: 137-147.View ArticlePubMedGoogle Scholar
- Giresi PG, Kim J, McDaniell RM, Iyer VR, Lieb JD: FAIRE (Formaldehyde-Assisted Isolation of Regulatory Elements) isolates active regulatory elements from human chromatin. Genome Res. 2007, 17: 877-885.PubMed CentralView ArticlePubMedGoogle Scholar
- Gardiner-Garden M, Frommer M: CpG islands in vertebrate genomes. J Mol Biol. 1987, 196: 261-282.View ArticlePubMedGoogle Scholar
- Kent WJ, Sugnet CW, Furey TS, Roskin KM, Pringle TH, Zahler AM, Haussler D: The human genome browser at UCSC. Genome Res. 2002, 12: 996-1006.PubMed CentralView ArticlePubMedGoogle Scholar
- Illingworth RS, Gruenewald-Schneider U, Webb S, Kerr AR, James KD, Turner DJ, Smith C, Harrison DJ, Andrews R, Bird AP: Orphan CpG islands identify numerous conserved promoters in the mammalian genome. PLoS Genet. 2010, 6: e1001134.PubMed CentralView ArticlePubMedGoogle Scholar
- Deaton AM, Bird A: CpG islands and the regulation of transcription. Genes Dev. 2011, 25: 1010-1022.PubMed CentralView ArticlePubMedGoogle Scholar
- Bird A, Taggart M, Frommer M, Miller OJ, Macleod D: A fraction of the mouse genome that is derived from islands of nonmethylated, CpG-rich DNA. Cell. 1985, 40: 91-99.View ArticlePubMedGoogle Scholar
- Cohen NM, Kenigsberg E, Tanay A: Primate CpG islands are maintained by heterogeneous evolutionary regimes involving minimal selection. Cell. 2011, 145: 773-786.View ArticlePubMedGoogle Scholar
- Langmead B, Salzberg SL: Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012, 9: 357-359.PubMed CentralView ArticlePubMedGoogle Scholar
- Boyle AP, Guinney J, Crawford GE, Furey TS: F-Seq: a feature density estimator for high-throughput sequence tags. Bioinformatics. 2008, 24: 2537-2538.PubMed CentralView ArticlePubMedGoogle Scholar
- Li QH, Brown JB, Huang HY, Bickel PJ: Measuring reproducibility of high-throughput experiments. Ann Appl Stat. 2011, 5: 1752-1779.View ArticleGoogle Scholar
- Edgar R, Domrachev M, Lash AE: Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002, 30: 207-210.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.