Epigenetic regulators play a critical role in normal and malignant hematopoiesis. Deregulation, including epigenetic deregulation, of the HOXA gene cluster drives transformation of about 50% of acute myeloid leukemia. We recently showed that the Histone 3 Lysine 9 methyltransferase SETDB1 negatively regulates the expression of the pro-leukemic genes Hoxa9 and its cofactor Meis1 through deposition of promoter H3K9 trimethylation in MLL-AF9 leukemia cells. Here, we investigated the biological impact of altered SETDB1 expression and changes in H3K9 methylation on acute myeloid leukemia. We demonstrate that SETDB1 expression is correlated to disease status and overall survival in acute myeloid leukemia patients. We recapitulated these findings in mice, where high expression of SETDB1 delayed MLL-AF9 mediated disease progression by promoting differentiation of leukemia cells. We also explored the biological impact of treating normal and malignant hematopoietic cells with an H3K9 methyltransferase inhibitor, UNC0638. While myeloid leukemia cells demonstrate cytotoxicity to UNC0638 treatment, normal bone marrow cells exhibit an expansion of cKit+ hematopoietic stem and progenitor cells. Consistent with these data, we show that bone marrow treated with UNC0638 is more amenable to transformation by MLL-AF9. Next generation sequencing of leukemia cells shows that high expression of SETDB1 induces repressive changes to the promoter epigenome and downregulation of genes linked with acute myeloid leukemia, including Dock1 and the MLL-AF9 target genes Hoxa9, Six1, and others. These data reveal novel targets of SETDB1 in leukemia that point to a role for SETDB1 in negatively regulating pro-leukemic target genes and suppressing acute myeloid leukemia.
Epigenetic deregulation has emerged as an important contributor to oncogenesis and disease progression in a variety of malignancies, including leukemia.1 Deep sequencing has revealed that genes encoding epigenetic modifying proteins are mutated in over 70% of acute myeloid leukemia (AML) patients.2 H3K9me2/3 marks large regions of condensed transcriptionally inactive chromatin, such as pericentric heterochromatin.3 H3K9me2/3 also plays a functional role in the dynamic repression of genes in euchromatic regions of the genome.4 Two families of proteins are associated with deposition of H3K9 methylation: the SUV39 family and PRDM family. The SUV39 family of H3K9 methyltransferases consists of SUV39H1/2, EHMT1/2, and SETDB1/2.4 Our lab and others have previously demonstrated that members of the SUV39 family of H3K9 methyltransferases bind to the polymerase associated factor complex (PAF1c).5,6 SETDB1, G9a (EHMT2), and GLP (EHMT1), were identified in a proteomics study exploring the interactome of the PAF1c in AML.5 The PAF1c is an epigenetic regulator complex that physically associates with RNA polymerase II (RNAPII) and both positively and negatively regulates gene transcription.7–11 In AML, the PAF1c is critical for the regulation of a proleukemic HOXA gene program in AML cells through the recruitment of MLL and MLL-fusion proteins to the HOXA locus via direct physical interactions.12–15 Hoxa9 and its co-factor MEIS1 are upregulated in about 50% of AML and associated with a poor patient prognosis.16 Given our recent data linking H3K9 methyltransferases with Hoxa9 and Meis1 repression along with altered H3K9me3 in AML patients compared to CD34+ cells,17 it is important to understand the epigenetic and biological impact of H3K9 methyltransferases on AML.
SETDB1 is a H3K9 mono/di/tri-methyltransferase involved in heterochromatin regulation and euchromatic gene silencing.18 SETDB1 was shown to bind gene loci associated with development in mouse embryonic stem (ES) cells, such as the Hoxd cluster of genes.19 SETDB1 has been implicated as an oncogene in melanoma, breast cancer, liver cancer, and lung cancer.4 Importantly, Ceol et al. reported amplification of SETDB1 in melanoma results in aberrant binding and regulation of the HOXA locus.20 In contrast to these oncogenic roles for SETDB1, Avgustinova and colleagues report that depletion of the H3K9 methyltransferase G9a in squamous tumors leads to a delayed, but more aggressive phenotype due to expanded cancer progenitor pools with increased genomic instability. 21 In the hematopoietic system, the methyltransferase activity of G9a is required for leukemogenesis due to a physical interaction with HOXA9.22 Importantly, loss of G9a has no effect on hematopoietic stem cells.22,23 SETDB1, however, is required for both hematopoietic stem and progenitor cell (HSPC) maintenance and leukemic stem cells.23 Further, Cuellar and colleagues show that SETDB1 mediated silencing of endogenous retroviral elements is required for the growth of AML cell lines.24 Together, these studies suggest that therapeutic targeting of SETDB1 may benefit AML patients. However, we recently demonstrated that SETDB1 negatively regulates the expression of the pro-leukemic Hoxa9 and Meis1 genes in MLL-AF9 transformed AML cells through association with the PAF1c, which localizes to Hoxa and Meis1 loci. The PAF1c-SETDB1 interaction mediates promoter H3K9me3 and repression of Hoxa9 and Meis1 expression.5 Further, SETDB1 expression is inversely correlated with HOXA9 and MEIS1 expression in AML patient samples.5 These data imply a more complex role for H3K9 methylation in AML similar to skin tumors whereby H3K9 methyltransferases display both oncogenic and suppressive roles.20,21 Thus, further investigation into the role of H3K9 methyltransferases in AML is required.
Here we show that AML patients with higher expression of SETDB1 display a better prognosis, consistent with repression of HOXA9 and MEIS1 expression. SETDB1 overexpression induces cellular differentiation and delays disease onset in a mouse model of AML, recapitulating AML patient survival. We also investigated the utility of inhibiting H3K9 methyltransferases in AML cells and HSPC, demonstrating that inhibition of H3K9 methylation in HSPC leads to retention of self-renewal capacity in HSPC and more efficient transformation by the MLL-AF9 fusion protein. Finally, we show that SETDB1 regulates gene expression by inducing changes in the epigenetic landscape and chromatin accessibility at gene targets critical to leukemogenesis.
Patient sample data
Data for patient gene expression relative to normal hematopoietic cells were mined from BloodPool on BloodSpot database.25 Bloodspot assigns AML patient samples to a closest normal hematopoietic counterpart using transcriptomic profiles.26 AML patient RNA sequencing (RNA-seq) and survival data were mined from The Cancer Genome Atlas.2
Cell line generation
Cell lines were generated from C57Bl/6 (Taconic Farms) mouse bone marrow or from SETDB1floxed27 mice. Platinum-E viral packaging cells were transfected with the indicated constructs: MSCVneo-FLAG-MLL-AF9 (MA9), MSCVneo-FLAG-E2A-HLF (EHF), MSCVhygro-FLAG-EHMT2 (G9a) (Ge lab; Addgene #41721), MSCVpuro-HA-SETDB1 or empty vector (EV) controls. Cells were spinfected with viral supernatants and 5 ug/mL polybrene (Millipore), selected with 1 mg/mL G418 (Invitrogen) and 1 ug/mL puromycin (Invitrogen) or 200 ug/mL hygromycin (Invitrogen) and cultured in IMDM with 15% stem cell fetal bovine serum (Millipore), 1% penicillin/ streptomycin (Invitrogen), 10 ng/mL interleukin-3 (IL-3) and 100 ng/mL stem cell factor (R&D).
Primary MLL-AF9 mouse leukemia cells were retrovirally transduced with MSCVpuro-HA-SETDB1 or EV and selected in 2 ug/mL puromycin for 4 days. 100,000 cells were injected intravenously into lethally irradiated (950 rads) C57Bl/6 mice. Mice were monitored for survival, moribund mice were euthanized, and bone marrow, spleen, and liver were harvested. Animal studies were approved by the University of Michigan’s Committee on Use and Care of Animals and Unit for Laboratory Medicine.
Quantitative PCR (qPCR)
RNA was harvested using the Qiagen RNeasy mini plus kit. cDNA synthesis was performed using the SuperScript III kit (Invitrogen). qPCR was performed using fast SYBR-green mastermix (Thermo Fisher). Primers sequences are listed in the Online Supplementary Table S2.
Sequencing libraries preparation
MA9+ empty vector (EV) and MA9+SETDB1 cells were harvested for next generation sequencing library preparation for RNAseq, ChIP-seq, and ATAC-seq. ChIP-seq antibodies were validated using Epicypher histone peptide arrays (Online Supplementary Figure S6). Library preparation details are in the Online Supplementary Materials and Methods.
Sequencing data is available via the Gene Expression Omnibus, accession GSE136850.
SETDB1 expression is correlated with AML patient prognosis
Given our data linking H3K9 methyltransferases with HOXA9 and MEIS1 expression, we investigated their expression in AML patient sample data. Interestingly, SETDB1, SUV39H1, and SUV39H2 exhibit lower expression in AML patient samples when compared with normal hematopoietic cells, with median expression levels that are 66%, 56%, and 41% relative to their nearest normal hematopoietic counterparts, respectively (Figure 1A).25,26 The downregulation of these genes in AML samples was consistent regardless of AML karyotype (Figure 1B and Online Supplementary Figure S1A-B).25,26 We next tested whether H3K9 methyltransferase gene expression significantly correlated with patient survival using publicly available RNA-seq data2 and found that only SETDB1 expression significantly correlated with patient survival (Figure 1C). Median survival for patients with SETDB1 expression above the median was 26.3 months and 9.5 months in patients with SETDB1 expression below the median (Figure 1D).2 Univariate and multivariate analyses reveal that higher expression of SETDB1 is associated with a higher overall survival rate with a P<0.003 and a lower expected hazard ratio of 0.29 per two-fold change in expression (Figure 1D-E).2
SETDB1 or G9A expression reduces AML growth and colony formation
We next explored the biological effects of SETDB1 expression on the transformation and growth of AML cells. First, we explored whether there is a difference in Setdb1 expression in mouse AML relative to normal HSPC by isolating lineage negative (Lin–) cKit+ cells from mouse bone marrow. qPCR demonstrates that Setdb1 expression is reduced in murine MLL-AF9 and CALM-AF10 leukemias compared to normal HSPC, consistent with the patient sample data (Figure 2A). We next performed colony replating assays where Lin– mouse bone marrow cells were retrovirally co-transduced with the MLL-AF9 fusion oncogene with and without SETDB1 and plated in semi-solid methylcellulose. Overexpression of human SETDB1 significantly reduced MLL-AF9 mediated colony formation (Figure 2B-C) while not affecting the expression of the exogenous MLL-AF9 (Online Supplementary Figure S2A).24 Colony formation driven by a separate leukemic fusion protein, E2A-HLF, was also reduced in the presence of SETDB1, suggesting a general effect on AML transformation (Online Supplementary Figure S2B-C). Ex vivo proliferation assays demonstrate that overexpression of SETDB1 in MLL-AF9 or E2A-HLF transformed AML cells leads to a significant reduction in cellular proliferation (Figure 2D and Online Supplementary Figure S2D). Interestingly, we observed a strong selective pressure to reduce exogenous SETDB1 expression in cultured MLLAF9+ SETDB1 cells (Online Supplementary Figure S2E) resulting in rescue of MLL-AF9 cellular proliferation (Online Supplementary Figure S2F). MLL-AF9 cells that overexpress SETDB1 undergo morphological changes consis-tent with differentiation (Figure 2E). Genes associated with hematopoietic differentiation, including Id2, Cd80, Nab2, and Itgam have significantly increased expression upon overexpression of SETDB1 (Figure 2F). Additionally, overexpression of SETDB1 in MLL-AF9 cells leads to increased apoptosis (Figure 2G and Online Supplementary Figure S2G). We next examined whether SETDB1 expression is associated with reduced chemoresistance. MA9+SETDB1 do not exhibit reduced apoptosis after treatment with Daunorubicin at their IC50 (Figure 2H and Online Supplementary Figure S2H-I), suggesting SETDB1 expression is not associated with increased Daunorubicin sensitivity. Overexpression of another H3K9 methyltransferase, human G9A, also reduces colony formation and proliferation of MLL-AF9 cells (Online Supplementary Figure S3A-C) and induces morphological changes consistent with differentiation (Online Supplementary Figure S3D). These data demonstrate that expression of multiple H3K9 methyltransferases reduces AML cell proliferation and colony forming potential and induces AML differentiation.
SETDB1 expression delays MLL-AF9 mediated AML
To examine the effects of SETDB1 in vivo, we transplanted primary mouse MLL-AF9 AML cells retrovirally transduced with or without SETDB1 into sublethally irradiated syngeneic recipient mice and monitored survival. Consistent with AML patient data, overexpression of SETDB1 significantly delays MLL-AF9 mediated leukemogenesis in vivo (Figure 2I). All moribund mice from both the control MLL-AF9 group and MLL-AF9+SETDB1 group exhibited splenomegaly, leukemic infiltration in the liver (Figure 2J and Online Supplementary Figure S3E) and similar MLL-AF9 expression levels (Online Supplementary Figure S3F). We measured expression of exogenous SETDB1 in MLL-AF9 leukemic cells before injection and in bone marrow of moribund mice and observed a decrease in exogenous SETDB1 expression in 6 of 8 mice (Figure 2K). These data demonstrate that SETDB1 expression suppresses MLL-AF9 mediated leukemic progression in vivo.
H3K9 methylation impairs HSPC colony formation and suppresses leukemic transformation
We examined how loss of Setdb1 affects normal HSPC and AML growth and transformation. First, using inducible CreER-mediated knock-out of Setdb1, we confirmed that Setdb1 is required for MLL-AF9 cell growth23 and demonstrated that Setdb1 is also required for the growth of E2A-HLF leukemic cells (Online Supplementary Figure S4A-B). We next sought to determine the effect of reduced expression but not deletion of Setdb1 on AML cell growth. However, heterozygous deletion of Setdb1 in MLL-AF9-Setdb1fl/+-CreER cells did not reduce Setdb1 protein (Online Supplementary Figure S4C). Therefore, we utilized short hairpin RNA (shRNA) mediated knockdown of Setdb1 using a doxycycline inducible MLL-AF9 Tet-on cell line in competitive growth assays. Interestingly, our results were inconclusive where two shRNA constructs with efficient knockdown displayed negative selection, while two less efficient shRNA showed a positive selection compared to MLL-AF9-Tet-on cells relative to the untreated control (Online Supplementary Figure S4D-E). This may suggest Setdb1 is maintained within a narrow expression window in AML cells.
To test the effect of Setdb1 knock-down on HSPC colony formation we transduced Lin– bone marrow with shRNA targeting Setdb1. Our results show a trend for increased HSPC colony formation without changes in Hoxa9 or Meis1 expression following Setdb1 knockdown (Online Supplementary Figure S4F, data not shown). Because inhibition of H3K9 methyltransferases has been proposed as a therapeutic option to treat AML,22,23 we explored the effects of chemical inhibition of H3K9 methylation in normal and malignant hematopoietic cells. Without a selective SETDB1 small molecule inhibitor, we utilized the G9a inhibitor UNC0638, which addresses the function of H3K9 methyltransferases more broadly. Lin– mouse bone marrow cells treated with UNC0638 exhibit a reduction in both H3K9me2 and H3K9me3 (Figure 3A).28,29 Given our data demonstrating Hoxa9 and Meis1 are repressed by H3K9 methylation,5 we tested the effect of UNC0638 on HSPC self-renewal. Previous reports have demonstrated that mouse Lin–Ska1+cKit+ cells (LSK) are preserved in culture following treatment with UNC0638.30 To explore this further, Lin– cells were treated with increasing doses of UNC0638 for five days prior to plating in semi-solid methylcellulose in the presence of SCF and IL-3. Interestingly, treatment with UNC0638 significantly increased colony formation capacity of Lin– cells in a dosedependent manner (Figure 3B-C). Further, human CD34+ cells isolated from mobilized peripheral blood and treated with increasing doses of UNC0638 demonstrate increased colony formation capacity (Figure 3D-E), consistent with reports that UNC0638 preserves CD34+ cells in culture.31 Thus, chemical or genetic inhibition of H3K9 methyltransferases preserves self-renewal capacity of bone marrow cells.
We also examined the effects of UNC0638 chemical inhibition of H3K9 methylation on AML cells. Consistent with previous studies, treatment with UNC0638 results in reduced cellular proliferation of MLL-AF9 cells (Online Supplementary Figure S5A).22 Since UNC0638 preserves self-renewal of HSPC (Figure 3), we asked whether inhibition of H3K9me alters MLL-AF9 mediated transformation of bone marrow cells. First, we found that UNC0638 treatment of Lin– bone marrow cells preserve more primitive cKit+ and Cd11b– populations (Figure 4A-D and Online Supplementary Figure S5B-C). Lin– cells display significantly increased Hoxa9 expression in response to UNC0638 treatment (Figure 4E). Lin–cKit+ (LK) and related primitive hematopoietic cell populations are more amenable to transformation than more differentiated subtypes.32 To explore whether this expansion of LK cells results in greater AML transformation capacity, we pretreated Lin– bone marrow cells with UNC0638 for four days then retrovirally transduced them with MigR1-MLL-AF9 and monitored green fluorescent protein (GFP) by flow cytometry. We observed a more rapid expansion of GFP+ MLLAF9 cells following treatment with UNC0638 compared to vehicle, with a 1.4-1.7 fold increase in GFP+ cells before both populations reached 100% GFP positivity (Figure 4FG and Online Supplementary Figure S5D). One of four replicates demonstrated a more rapid expansion of vehicle treated cells; however, GFP+ cells were increased two-fold in vehicle treated cells at day 1 suggesting this is attributable to differences in transduction rates. Our combined data suggests that H3K9 methyltransferases can suppress leukemic transformation and may point to a narrow window of H3K9 methylation that is optimal for leukemic transformation and cell growth.
SETDB1 regulates oncogenic gene programs in AML
To explore the gene programs regulated by SETDB1 in AML, we performed RNA-sequencing experiments on MLL-AF9 cells overexpressing SETDB1. 2,285 genes are upregulated and 1,771 genes are downregulated by SETDB1 overexpression with a fold-change of 1.5 or more and an false discovery rate (FDR) of <0.05 (Figure 5A and Online Supplementary Table S4). Notably, many genes found in the Hoxa cluster were significantly downregulated, including Hoxa3, Hoxa5, Hoxa6, Hoxa9, and the Hoxa9 cofactor Meis1; while only one Hoxa gene was upregulated, the long non-coding RNA Hoxa11os (Figure 5A, highlighted genes). In fact, gene programs that are upregulated by forced expression of HOXA9 and MEIS1 in mouse cells are significantly downregulated by SETDB1 overexpression using gene set analysis33 (Figure 5B). We show significant upregulation of genes that exhibit increased expression in mature blood cells compared to HSPC, consistent with the differentiation observed upon SETDB1 overexpression in AML34,35 (Figure 5C-D and Figure 2F-G). Because SETDB1 binds the PAF1c, which is required for localization of MLL fusion proteins,5 we asked how SETDB1 expression affects direct targets of MLL-AF9. Interestingly, genes bound by regulation of MLL-AF9 were significantly downregulated upon SETDB1 overexpression, 36 suggesting H3K9me3 regulation MLL-AF9 gene programs in leukemic cells (Figure 5E). 193 genes downregulated by SETDB1 overexpression are reported to have reduced promoter H3K9 methylation in AML relative to normal CD34+ cells,17 suggesting SETDB1 may be responsible for regulating a subset of these genes, including Kit, Cbl, Ptpn11, Six1, and other genes that are important in AML (Figure 5F and Online Supplementary Table S4). There is significant overlap between genes downregulated by SETDB1 and genes upregulated by Crispr/Cas9 mediated knockdown of SETDB1 in human THP-1 AML cells harboring an MLL-AF9 fusion (Figure 5E and Online Supplementary Table S4).24 This suggests SETDB1 regulates conserved pro-leukemic gene programs in leukemic cells including direct MLL fusion targets.
SETDB1 regulates the epigenome to affect changes in chromatin accessibility and gene expression
We have shown that overexpression of SETDB1 in MLL-AF9 cells leads to global increases in H3K9me3.5 To understand the specific epigenomic changes induced by SETDB1 in AML, we performed ChIP-seq for H3K9me3 in MLL-AF9 and MLL-AF9+SETDB1 leukemic cells. We also performed ATAC-seq to assess changes in chromatin accessibility. We first explored differences at the HOXA locus due to its importance in a large subset of AML, including MLL leukemias.5,37 Overexpression of SETDB1 reduced chromatin accessibility and increased H3K9me3 at posterior Hoxa genes, which results in reduced transcription of Hoxa9 (Figure 6A). We sought to define whole genome epigenetic regulation mediated by SETDB1 but observed only 552 consensus H3K9me3 peaks (Online Supplementary Figure S7A and Online Supplementary Table S5), which included regions enriched for H3K9me3, such as the zinc finger protein cluster on chromosome 7 (Online Supplementary Figure S7B). Repetitive elements may consume H3K9me3 sequencing depth making it difficult to map with stringency. We circumvented this by performing ChIP-Seq for H3K9ac, which is mutually exclusive of H3K9me3 and associated with gene activation, in contrast to H3K9me3.38 We performed H3K9ac ChIP-seq on MLLAF9 and MLL-AF9+SETDB1 cells and saw changes at over 6,000 promoter regions, including both increased and decreased H3K9ac signal (Online Supplementary Figure S7C and Online Supplementary Table S5). Gene ontology analysis reveals that genes with decreased promoter H3K9ac in MLL-AF9+SETDB1 cells were associated with cell cycle and RNA binding, whereas genes with increased promoter H3K9ac were associated with signaling pathways and negative regulation of proliferation (Figure 6B). Overlapping downregulated genes with genes that exhibit reduced promoter ATAC-seq and H3K9ac ChIP-seq signal in MLL-AF9+SETDB1 cells and scoring by their combined fold changes reveals several interesting targets including Six1 and Mefv, which are implicated as biomarkers in AML patients39,40 (Figure 6C-D and Online Supplementary Table S5). We also observed significant loss of promoter H3K9ac and chromatin accessibility and reduced gene expression of Dock1 (Figure 6E), which is a prognostic marker of AML that displays changes in DNA methylation in AML patient samples relative to normal HSPC.41,42 We identified several direct binding targets of MLL-AF9 that underwent epigenetic remodeling and expression changes with overexpression of SETDB1. Thus, we asked whether SETDB1 impacts H3K79me2, which is deposited by DOT1L and associated with MLL-fusion proteins,36 using ChIP-seq for H3K79me2. We found a marked decrease of H3K9ac and ATAC-seq signal at H3K79me2 peaks in MLL-AF9+SETDB1 cells suggesting a role for H3K9 modifications in regulating genes marked with H3K79me2 (Figure 7A-B). To further explore the role of SETDB1 in regulating MLL-AF9 targets, we performed an overlap analysis for genes that, upon SETDB1 overexpression, are downregulated, lose promoter H3K9ac, compact chromatin, lose gene body H3K79me2 signal, and are bound by MLL-AF9.36 This defined a target list of genes that may be coregulated by SETDB1 and MLL-AF9. Included in this group is Gfi1, which has been shown to affect AML cell growth;43 Rap1gds1, a nucleotide exchange factor; Arid1b, a member of the SWI/SNF complex; and Six1, which promotes formation of leukemic stem cells40 (Figure 7C). Six1 is of particular interest given its role in promoting leukemogenesis and the striking reductions observed in H3K9ac, H3K79me2, ATAC-seq signal, and gene expression at this locus (Figure 7D). We investigated H3K9me3 levels at specific loci by performing ChIP-qPCR at the promoter of Six1 and Dock1. We detected increased H3K9me3 in MLL-AF9 cells overexpressing SETDB1 (Figure 7E). Additionally, ENCODE data from human K562 cells shows SETDB1 binding at SIX1 (Figure 7F), GFI1 (Figure S7D), RAP1GDS1, and ARID1B (data not shown), but not DOCK1 (data not shown), suggesting we identified both direct and indirect targets. To determine whether SETDB1 expression inhibits cell growth primarily through Hoxa9 repression, we asked whether HOXA9/MEIS1 driven colony formation is affected by SETDB1 overexpression. A modest but insignificant decrease in colony formation is observed in HOXA9/MEIS1 transformed cells following SETDB1 overexpression (Figure 8A). However, SETDB1 overexpression is accompanied by significant upregulation of exogenous Meis1 (Figure 8B and Online Supplementary Figure S8A-B) that may account for the modest effects. Thus, Hoxa9, Meis1, and their downstream targets are likely affected by SETDB1 to influence AML cell growth. Taken together, our data demonstrate profound effects of SETDB1 on the epigenome and transcriptome that affects genes critical for AML.
Here we demonstrate that SETDB1 and H3K9 methylation suppresses AML disease progression in vivo through the repression of pro-leukemic genes including direct MLL-fusion protein targets. We found that AML patient samples exhibit lower SETDB1 expression compared to normal hematopoietic cells and that higher SETDB1 expression correlates strongly with better overall AML patient survival (Figure 1). We recapitulated these findings in mice where forced expression of SETDB1 in MLL-AF9 driven AML induces differentiation of AML cells and increases disease latency (Figure 2). These data suggest SETDB1 suppresses AML cell growth and self-renewal by relieving the block in differentiation.
We attribute the phenotypes in AML cells to altered H3K9 methylation. We altered H3K9 methylation levels genetically (SETDB1 and G9a) and through small molecule inhibition (UNC0638). Similar to our results with SETDB1, manipulation of G9a suggests H3K9 methylation can suppress AML progression by promoting differentiation (Figures 3-4, and Online Supplementary Figure S3). Thus, H3K9 methylation may have a more general effect on AML initiation and progression. Interestingly, Lehnertz and colleagues reported G9a overexpression accelerated Hoxa9/Meis1 mediated leukemia in vivo. We found Hoxa9/Meis1 mediated transformation in vitro was inhibited by SETDB1, but to a lesser extent than MLL-AF9 (Figure 8). Unique experimental strategies or functions for SETDB1 and G9a may account for these differences.22 Consistent with a role for H3K9 methylation in suppressing hematopoietic transformation, deletion of the H3K9 methyltransferase, SUV39H1 (and to a lesser degree SUV39H2), leads to the development of B-cell lymphomas in mice.44 Additionally, SETDB2 resides in a region of chromosome 13 that is commonly deleted in chronic lymphocytic leukemia (CLL).45 Thus, H3K9 methylation is likely exquisitely regulated in hematopoietic cells and performs context dependent functions that require further investigation to fully understand its role in AML.
Mechanistically, we found that SETDB1 is linked with altered H3K9 methylation and acetylation, decreased chromatin accessibility and transcriptional repression of critical AML oncogenes (Figures 5-7). These genes included several that have been implicated in myelodysplastic syndromes (MDS) and AML.39 We show SETDB1 regulates Dock1 expression, which is correlated with leukemic stem cell gene signatures and a poor prognosis in AML patients.41,46 We also observed that SETDB1 represses genes associated with AML, such as Hoxa9 and Six1, which are direct targets of MLL-AF9.16,37,40 Interestingly, Six1 was recently shown to be important for leukemic stem cell maintenance where depletion of Six1 leads to increased disease latency.47 This suggests repression of Six1 may contribute to SETDB1 mediated extension of leukemic disease latency. These data point to SETDB1 negatively regulating a pro-leukemic gene program, many of which are potential therapeutic targets. Thus, understanding the mechanisms regulating SETDB1 at the transcriptional and post-translational level may be a valuable therapeutic approach for AML. For example, miRNA29 is a critical mediator of SETDB1 expression.48 Another potential mediator of H3K9 methylation is the PAF1c. We showed SETDB1 binds to the PAF1c and mediates promoter H3K9me3 of the Hoxa9 and Meis1 loci.5 Further, we and others identified G9a and SUV39H1 as interacting partners of the PAF1c.5,6 Interestingly, SETDB1, G9a, GLP and SUV39H1 form a complex that directs H3K9 methylation to euchromatic gene promoters.49 Thus, the PAF1c may recruit H3K9 methyltransferases to specific targets to mediate gene repression. The PAF1c is a critical regulator of transcription of several pro-leukemic genes in AML cells through direct physical interaction with wild-type MLL and MLL-fusion proteins.13,15 It will be interesting to consider the biochemical interplay between H3K9 methyltransferases and MLL-fusion proteins with the PAF1c.
Previous studies have demonstrated that SETDB1 and G9a are required for AML initiation and progression.22–24 Our current data demonstrating that SETDB1 suppresses AML growth may suggest AML cells maintain a narrow SETDB1 expression level. We show increased SETDB1 expression induces differentiation of AML cells through H3K9me3 and repression of self-renewal genes. Conversely, loss of SETDB1 is detrimental to leukemic cells due to derepression of endogenous retroviral elements (ERV) and inhibition of HOXA9 transcriptional activity22,24 (Figure 8C). Given the essential role for SETDB1 in leukemia, small molecule inhibition of H3K9 methyltransferases has been proposed as a therapeutic option.22,23 However, a recent study shows depletion of G9a increased cancer progenitor cell populations that initiate a delayed but more aggressive disease state.21 Thus, it is critical to fully understand the effects of chemically inhibiting of H3K9 methylation as a treatment for AML. Further investigation into the roles of SETDB1, G9a and more generally H3K9 methylation levels will likely shed light on the precise role of these methyltransferases in normal and malignant hematopoiesis and determine the value of these epigenetic modifiers as therapeutic targets.
- Received April 4, 2019
- Accepted September 25, 2019
This work was supported by NIH grants R01-HL-136420 (AGM), T32CA140044 (JR) and P30CA046592 (University of Michigan Flow Cytometry Core).
The authors thank Dr. Jianyong Shou for a SETDB1 expression construct and Dr. Schahram Akbarian for Setdb1fl/fl mice. We also thank Dr. Russell Ryan, Dr. Emmalee Adelman, Dr. Maria Figueroa, and Dr. Sami Malek for helpful discussion.
- Flavahan WA, Gaskell E, Bernstein BE. Epigenetic plasticity and the hallmarks of cancer. Science. 2017; 357(6348)https://doi.org/10.1126/science.aal2380PubMedPubMed CentralGoogle Scholar
- Cancer Genome Atlas Research Network, Ley TJ, Miller C. Genomic and epigenomic landscapes of adult de novo acute myeloid leukemia. N Engl J Med. 2013; 368(22):2059-2074. https://doi.org/10.1056/NEJMoa1301689PubMedPubMed CentralGoogle Scholar
- Wang J, Jia ST, Jia S. New insights into the regulation of heterochromatin. Trends Genet. 2016; 32(5):284-294. https://doi.org/10.1016/j.tig.2016.02.005PubMedPubMed CentralGoogle Scholar
- Rao VK, Pal A, Taneja R. A drive in SUVs: from development to disease. Epigenetics. 2017; 12(3):177-186. https://doi.org/10.1080/15592294.2017.1281502PubMedPubMed CentralGoogle Scholar
- Ropa J, Saha N, Chen Z. PAF1 complex interactions with SETDB1 mediate promoter H3K9 methylation and transcriptional repression of Hoxa9 and Meis1 in acute myeloid leukemia. Oncotarget. 2018; 9(31):22123-22136. https://doi.org/10.18632/oncotarget.25204PubMedPubMed CentralGoogle Scholar
- Yang Y-J, Han J-W, Youn H-D, Cho E-J. The tumor suppressor, parafibromin, mediates histone H3 K9 methylation for cyclin D1 repression. Nucleic Acids Res. 2010; 38(2):382-390. https://doi.org/10.1093/nar/gkp991PubMedPubMed CentralGoogle Scholar
- Jaehning JA. The Paf1 complex: platform or player in RNA polymerase II transcription?. Biochim Biophys Acta. 2010; 1799(5-6):379-388. https://doi.org/10.1016/j.bbagrm.2010.01.001PubMedPubMed CentralGoogle Scholar
- Van Oss SB, Cucinotta CE, Arndt KM. Emerging insights into the roles of the Paf1 complex in gene regulation. Trends Biochem Sci. 2017; 42(10):788-798. https://doi.org/10.1016/j.tibs.2017.08.003PubMedPubMed CentralGoogle Scholar
- Chen FX, Xie P, Collings CK. PAF1 regulation of promoter-proximal pause release via enhancer activation. Science. 2017; 357(6357):1294-1298. https://doi.org/10.1126/science.aan3269PubMedPubMed CentralGoogle Scholar
- Krogan NJ, Dover J, Wood A. The Paf1 complex is required for histone H3 methylation by COMPASS and Dot1p: linking transcriptional elongation to histone methylation. Mol Cell. 2003; 11(3):721-729. https://doi.org/10.1016/S1097-2765(03)00091-1Google Scholar
- Yu M, Yang W, Ni T. RNA polymerase II-associated factor 1 regulates the release and phosphorylation of paused RNA polymerase II. Science. 2015; 350(6266):1383-1386. https://doi.org/10.1126/science.aad2338PubMedGoogle Scholar
- Serio J, Ropa J, Chen W. The PAF complex regulation of Prmt5 facilitates the progression and maintenance of MLL fusion leukemia. Oncogene. 2018; 37(4):450-460. https://doi.org/10.1038/onc.2017.337PubMedPubMed CentralGoogle Scholar
- Milne TA, Kim J, Wang GG. Multiple interactions recruit MLL1 and MLL1 fusion proteins to the HOXA9 locus in leukemogenesis. Mol Cell. 2010; 38(6):853-863. https://doi.org/10.1016/j.molcel.2010.05.011PubMedPubMed CentralGoogle Scholar
- Muntean AG, Chen W, Jones M, Granowicz EM, Maillard I, Hess JL. MLL fusion protein-driven AML is selectively inhibited by targeted disruption of the MLL-PAFc interaction. Blood. 2013; 122(11):1914-1922. https://doi.org/10.1182/blood-2013-02-486977PubMedPubMed CentralGoogle Scholar
- Muntean AG, Tan J, Sitwala K. The PAF complex synergizes with MLL fusion proteins at HOX loci to promote leukemogenesis. Cancer Cell. 2010; 17(6):609-621. https://doi.org/10.1016/j.ccr.2010.04.012PubMedPubMed CentralGoogle Scholar
- Collins CT, Hess JL. Role of HOXA9 in leukemia: dysregulation, cofactors and essential targets. Oncogene. 2016; 35(9):1090-1098. https://doi.org/10.1038/onc.2015.174PubMedPubMed CentralGoogle Scholar
- Müller-Tidow C, Klein H-U, Hascher A. Profiling of histone H3 lysine 9 trimethylation levels predicts transcription factor activity and survival in acute myeloid leukemia. Blood. 2010; 116(18):3564-3571. https://doi.org/10.1182/blood-2009-09-240978PubMedPubMed CentralGoogle Scholar
- Schultz DC, Ayyanathan K, Negorev D, Maul GG, Rauscher FJ. SETDB1: a novel KAP-1-associated histone H3, lysine 9-specific methyltransferase that contributes to HP1-mediated silencing of euchromatic genes by KRAB zinc-finger proteins. Genes Dev. 2002; 16(8):919-932. https://doi.org/10.1101/gad.973302PubMedPubMed CentralGoogle Scholar
- Bilodeau S, Kagey MH, Frampton GM, Rahl PB, Young RA. SetDB1 contributes to repression of genes encoding developmental regulators and maintenance of ES cell state. Genes Dev. 2009; 23(21):2484-2489. https://doi.org/10.1101/gad.1837309PubMedPubMed CentralGoogle Scholar
- Ceol CJ, Houvras Y, Jane-Valbuena J. The histone methyltransferase SETDB1 is recurrently amplified in melanoma and accelerates its onset. Nature. 2011; 471(7339):513-517. https://doi.org/10.1038/nature09806PubMedPubMed CentralGoogle Scholar
- Avgustinova A, Symeonidi A, Castellanos A. Loss of G9a preserves mutation patterns but increases chromatin accessibility, genomic instability and aggressiveness in skin tumours. Nat Cell Biol. 2018; 20(12):1400-1409. https://doi.org/10.1038/s41556-018-0233-xPubMedGoogle Scholar
- Lehnertz B, Pabst C, Su L. The methyltransferase G9a regulates HoxA9- dependent transcription in AML. Genes Dev. 2014; 28(4):317-327. https://doi.org/10.1101/gad.236794.113PubMedPubMed CentralGoogle Scholar
- Koide S, Oshima M, Takubo K. Setdb1 maintains hematopoietic stem and progenitor cells by restricting the ectopic activation of nonhematopoietic genes. Blood. 2016; 128(5):638-649. https://doi.org/10.1182/blood-2016-01-694810PubMedGoogle Scholar
- Cuellar TL, Herzner A-M, Zhang X. Silencing of retrotransposons by SETDB1 inhibits the interferon response in acute myeloid leukemia. J Cell Biol. 2017; 216(11):3535-3549. https://doi.org/10.1083/jcb.201612160PubMedPubMed CentralGoogle Scholar
- Bagger FO, Sasivarevic D, Sohi SH. BloodSpot: a database of gene expression profiles and transcriptional programs for healthy and malignant haematopoiesis. Nucleic Acids Res. 2016; 44(D1):D917-24. https://doi.org/10.1093/nar/gkv1101PubMedPubMed CentralGoogle Scholar
- Rapin N, Bagger FO, Jendholm J. Comparing cancer vs normal gene expression profiles identifies new disease entities and common transcriptional programs in AML patients. Blood. 2014; 123(6):894-904. https://doi.org/10.1182/blood-2013-02-485771PubMedGoogle Scholar
- Jiang Y, Loh Y-HE, Rajarajan P. The methyltransferase SETDB1 regulates a large neuron-specific topological chromatin domain. Nat Genet. 2017; 49(8):1239-1250. Google Scholar
- Kim Y, Lee H-M, Xiong Y. Targeting the histone methyltransferase G9a activates imprinted genes and improves survival of a mouse model of Prader-Willi syndrome. Nat Med. 2017; 23(2):213-222. https://doi.org/10.1038/nm.4257PubMedPubMed CentralGoogle Scholar
- Vedadi M, Barsyte-Lovejoy D, Liu F. A chemical probe selectively inhibits G9a and GLP methyltransferase activity in cells. Nat Chem Biol. 2011; 7(8):566-574. https://doi.org/10.1038/nchembio.599PubMedPubMed CentralGoogle Scholar
- Ugarte F, Sousae R, Cinquin B. Progressive chromatin condensation and H3K9 methylation regulate the differentiation of embryonic and hematopoietic stem cells. Stem Cell Rep. 2015; 5(5):728-740. https://doi.org/10.1016/j.stemcr.2015.09.009PubMedPubMed CentralGoogle Scholar
- Chen X, Skutt-Kakaria K, Davison J. G9a/GLP-dependent histone H3K9me2 patterning during human hematopoietic stem cell lineage commitment. Genes Dev. 2012; 26(22):2499-2511. https://doi.org/10.1101/gad.200329.112PubMedPubMed CentralGoogle Scholar
- Krivtsov AV, Figueroa ME, Sinha AU. Cell of origin determines clinically relevant subtypes of MLL-rearranged AML. Leukemia. 2013; 27(4):852-860. https://doi.org/10.1038/leu.2012.363PubMedPubMed CentralGoogle Scholar
- Hess JL, Bittner CB, Zeisig DT. c-Myb is an essential downstream target for homeobox-mediated transformation of hematopoietic cells. Blood. 2006; 108(1):297-304. https://doi.org/10.1182/blood-2005-12-5014PubMedPubMed CentralGoogle Scholar
- Ivanova NB, Dimos JT, Schaniel C, Hackney JA, Moore KA, Lemischka IR. A stem cell molecular signature. Science. 2002; 298(5593):601-604. https://doi.org/10.1126/science.1073823PubMedGoogle Scholar
- Laslo P, Spooner CJ, Warmflash A. Multilineage transcriptional priming and determination of alternate hematopoietic cell fates. Cell. 2006; 126(4):755-766. https://doi.org/10.1016/j.cell.2006.06.052PubMedGoogle Scholar
- Bernt KM, Zhu N, Sinha AU. MLLrearranged leukemia is dependent on aberrant H3K79 methylation by DOT1L. Cancer Cell. 2011; 20(1):66-78. https://doi.org/10.1016/j.ccr.2011.06.010PubMedPubMed CentralGoogle Scholar
- Faber J, Krivtsov AV, Stubbs MC. HOXA9 is required for survival in human MLL-rearranged acute leukemias. Blood. 2009; 113(11):2375-2385. https://doi.org/10.1182/blood-2007-09-113597PubMedPubMed CentralGoogle Scholar
- Ghare SS, Joshi-Barve S, Moghe A. Coordinated histone H3 methylation and acetylation regulate physiologic and pathologic fas ligand gene expression in human CD4+ T cells. J Immunol. 2014; 193(1):412-421. https://doi.org/10.4049/jimmunol.1400055PubMedPubMed CentralGoogle Scholar
- Celik S, Tangi F, Oktenli C. Increased frequency of Mediterranean fever gene variants in multiple myeloma. Oncol Lett. 2014; 8(4):1735-1738. https://doi.org/10.3892/ol.2014.2407PubMedPubMed CentralGoogle Scholar
- Wang Q-F, Wu G, Mi S. MLL fusion proteins preferentially regulate a subset of wild-type MLL target genes in the leukemic genome. Blood. 2011; 117(25):6895-6905. https://doi.org/10.1182/blood-2010-12-324699PubMedPubMed CentralGoogle Scholar
- Lee S-H, Chiu Y-C, Li Y-H. High expression of dedicator of cytokinesis 1 (DOCK1) confers poor prognosis in acute myeloid leukemia. Oncotarget. 2017; 8(42):72250-72259. https://doi.org/10.18632/oncotarget.19706PubMedPubMed CentralGoogle Scholar
- Jung N, Dai B, Gentles AJ, Majeti R, Feinberg AP. An LSC epigenetic signature is largely mutation independent and implicates the HOXA cluster in AML pathogenesis. Nat Commun. 2015; 6:8489. https://doi.org/10.1038/ncomms9489PubMedPubMed CentralGoogle Scholar
- Volpe G, Walton DS, Grainger DE. Prognostic significance of high GFI1 expression in AML of normal karyotype and its association with a FLT3-ITD signature. Sci Rep. 2017; 7(1):11148. https://doi.org/10.1038/s41598-017-11718-8PubMedPubMed CentralGoogle Scholar
- Peters AH, O’Carroll D, Scherthan H. Loss of the Suv39h histone methyltransferases impairs mammalian heterochromatin and genome stability. Cell. 2001; 107(3):323-337. https://doi.org/10.1016/S0092-8674(01)00542-6Google Scholar
- Falandry C, Fourel G, Galy V. CLLD8/KMT1F is a lysine methyltransferase that is important for chromosome segregation. J Biol Chem. 2010; 285(26):20234-20241. https://doi.org/10.1074/jbc.M109.052399PubMedPubMed CentralGoogle Scholar
- Zhang Y, Li Y, Zhang L, Fang X, Zhen C, Wang X. Prognostic importance of DOCK1 transcript levels, and biologic insights from DOCK1 -associated gene and microrna expression signatures in de novo acute myeloid leukemia. Blood. 2017; 130(Suppl 1):3952. Google Scholar
- Chu Y, Chen Y, Li M. Six1 regulates leukemia stem cell maintenance in acute myeloid leukemia. Cancer Sci. 2019; 110(7):2200-2210. https://doi.org/10.1111/cas.14033PubMedPubMed CentralGoogle Scholar
- Wong C-M, Wei L, Law C-T. Up-regulation of histone methyltransferase SETDB1 by multiple mechanisms in hepatocellular carcinoma promotes cancer metastasis. Hepatology. 2016; 63(2):474-487. https://doi.org/10.1002/hep.28304PubMedGoogle Scholar
- Fritsch L, Robin P, Mathieu JRR. A subset of the histone H3 lysine 9 methyltransferases Suv39h1, G9a, GLP, and SETDB1 participate in a multimeric complex. Mol Cell. 2010; 37(1):46-56. https://doi.org/10.1016/j.molcel.2009.12.017PubMedGoogle Scholar
Figures & Tables
This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License.