Yan Zi Au1,2,3, Muxin Gu2, Etienne De Braekeleer2, Malgorzata Gozdecka2,4, Demetrios Aspris2, Yusuke Tarumoto5, Jonathan Cooper2, Jason Yu1,6, Swee Hoe Ong1, Xi Chen7, Konstantinos Tzelepis2,8, Brian J P Huntly4,9,10, George Vassiliou11,12,13, Kosuke Yusa14,15. 1. Stem Cell Genetics, Wellcome Sanger Institute, Hinxton, Cambridge, UK. 2. Haematological Cancer Genetics, Wellcome Sanger Institute, Hinxton, Cambridge, UK. 3. Dana-Farber Cancer Institute, Boston, MA, USA. 4. Wellcome Trust-MRC Stem Cell Institute, Cambridge Biomedical Campus, University of Cambridge, Cambridge, UK. 5. Stem Cell Genetics, Institute for Frontier Life and Medical Sciences, Kyoto University, Kyoto, Japan. 6. Department of Cell Biology, The Francis Crick Institute, London, UK. 7. Gene Expression Genomics, Wellcome Sanger Institute, Hinxton, Cambridge, UK. 8. Gurdon Institute, University of Cambridge, Tennis Court Road, Cambridge, UK. 9. Department of Haematology, Cambridge University Hospitals NHS Trust, Cambridge, UK. 10. Cambridge Institute for Medical Research, University of Cambridge, Cambridge, UK. 11. Haematological Cancer Genetics, Wellcome Sanger Institute, Hinxton, Cambridge, UK. gsv20@sanger.ac.uk. 12. Wellcome Trust-MRC Stem Cell Institute, Cambridge Biomedical Campus, University of Cambridge, Cambridge, UK. gsv20@sanger.ac.uk. 13. Department of Haematology, Cambridge University Hospitals NHS Trust, Cambridge, UK. gsv20@sanger.ac.uk. 14. Stem Cell Genetics, Wellcome Sanger Institute, Hinxton, Cambridge, UK. k.yusa@infront.kyoto-u.ac.jp. 15. Stem Cell Genetics, Institute for Frontier Life and Medical Sciences, Kyoto University, Kyoto, Japan. k.yusa@infront.kyoto-u.ac.jp.
Abstract
Histone acetyltransferases (HATs) catalyze the transfer of an acetyl group from acetyl-CoA to lysine residues of histones and play a central role in transcriptional regulation in diverse biological processes. Dysregulation of HAT activity can lead to human diseases including developmental disorders and cancer. Through genome-wide CRISPR-Cas9 screens, we identified several HATs of the MYST family as fitness genes for acute myeloid leukemia (AML). Here we investigate the essentiality of lysine acetyltransferase KAT7 in AMLs driven by the MLL-X gene fusions. We found that KAT7 loss leads to a rapid and complete loss of both H3K14ac and H4K12ac marks, in association with reduced proliferation, increased apoptosis, and differentiation of AML cells. Acetyltransferase activity of KAT7 is essential for the proliferation of these cells. Mechanistically, our data propose that acetylated histones provide a platform for the recruitment of MLL-fusion-associated adaptor proteins such as BRD4 and AF4 to gene promoters. Upon KAT7 loss, these factors together with RNA polymerase II rapidly dissociate from several MLL-fusion target genes that are essential for AML cell proliferation, including MEIS1, PBX3, and SENP6. Our findings reveal that KAT7 is a plausible therapeutic target for this poor prognosis AML subtype.
Histone acetyltransferases (HATs) catalyze the transfer of an acetyl group from acetyl-CoA to lysine residues of histones and play a central role in transcriptional regulation in diverse biological processes. Dysregulation of HAT activity can lead to human diseases including developmental disorders and cancer. Through genome-wide CRISPR-Cas9 screens, we identified several HATs of the MYST family as fitness genes for acute myeloid leukemia (AML). Here we investigate the essentiality of lysine acetyltransferase KAT7 in AMLs driven by the MLL-X gene fusions. We found that KAT7 loss leads to a rapid and complete loss of both H3K14ac and H4K12ac marks, in association with reduced proliferation, increased apoptosis, and differentiation of AML cells. Acetyltransferase activity of KAT7 is essential for the proliferation of these cells. Mechanistically, our data propose that acetylated histones provide a platform for the recruitment of MLL-fusion-associated adaptor proteins such as BRD4 and AF4 to gene promoters. Upon KAT7 loss, these factors together with RNA polymerase II rapidly dissociate from several MLL-fusion target genes that are essential for AML cell proliferation, including MEIS1, PBX3, and SENP6. Our findings reveal that KAT7 is a plausible therapeutic target for this poor prognosis AML subtype.
Acute myeloid leukemia (AML) is an aggressive malignancy of haematopoietic stem and progenitor cells. Translocations involving the mixed lineage leukemia gene (MLL or KMT2A) characterize an aggressive subtype of the disease and are associated with an intermediate or poor prognosis [1-3]. Currently, chemotherapy based on a combination of anthracyclines and purine analogues is the standard of care for AML, yet patients with leukemias driven by MLL-X gene fusions commonly become refractory to such treatments [4]. Advances in our understanding of the molecular pathogenesis of acute leukemias driven by MLL-X fusion oncogenes have led to the development of new drugs, including inhibitors of bromodomain proteins [5-7], DOT1L [8,9] and MENIN [10]; however, none of these have demonstrated a significant clinical benefit for patients as yet. In order to identify new therapeutic approaches in AML, we and others have performed genome-wide CRISPR-Cas9 fitness screens in different subtypes of the disease including the MLL-X subtype [11,12]. Through our own screen data, we recently showed that GCN5 (KAT2A) and SRPK1 are novel vulnerability in AML and that sensitivity to SRPK1 inhibition is particularly associated with AML driven by MLL-X fusion oncogenes [11,13]. Systematic detailed analysis of our CRISPR screen dataset revealed that several members of the MYST family of histone acetyltransferases (HATs) are required for the survival of AML cell lines, with KAT7 displaying strong essentiality for AMLs driven by MLL-X oncogenes.HATs are classified based on structural and sequence homology into distinct groups including the p300/CBP, MYST and GCN5 families [14]. They play a central role in transcriptional regulation through their function to catalyze the transfer of acetyl groups from acetyl-CoA to the ε-amino group of lysine residues of histones. Dysregulation of HATs is known to be associated with human diseases including developmental disorders and cancer [14-16]. The MYST family of HATs, characterized by their conserved MYST catalytic domain, includes KAT5 (TIP60), KAT6A (MOZ and MYST3), KAT6B (MORF and MYST4), KAT7 (HBO1 and MYST2) and KAT8 (MOF) [15,16]. While KAT6A and KAT6B are known targets of chromosomal translocations that drive AML [17-21] and KAT8 can play a role in some cases of AML [22], it is currently not known whether KAT7 plays a role in leukemogenesis. As KAT7 lacks chromatin binding domains, it relies on forming complexes with a scaffold protein (JADE or BRPF), along with EAF6 and ING4/5, for efficient histone acetylation [23]. The choice of the scaffold protein primarily determines target histones, namely H3 or H4. The KAT7-BRPF acetylates histone H3 tails (K14 and K23) [23-26], whereas KAT7-JADE complexes acetylates histone H4 tails (K5, K8, and K12) [27], respectively. In development, abolishing KAT7 results in major reduction of H3K14ac in both mouse fetal liver erythroblasts [24] and mouse embryo [25]; the latter associates with lethality.Here we investigate the molecular mechanism underlying the requirement of KAT7 in MLL-X-driven AMLs. We report that loss of KAT7 leads to reduced proliferation and enhanced apoptosis and differentiation. We also show that KAT7 is essential for maintaining the transcriptional program driven by MLL-AF9, through the recruitment of BRD4 and other MLL-fusion associated proteins such as AFF1 to the promoters of a subset of MLL-AF9 target genes. Together, our findings propose KAT7 as a plausible therapeutic target in MLL-X AML and also provide novel mechanistic insights into the molecular basis of MLL-X-driven transcriptional dysregulation.
Materials and Methods
Cell culture
Human cell lines (MOLM-13, MV4-11, THP-1, HL-6, K562, Nomo-1, OCI-AML2, OCI-AML3, MA9-ITD and MA9-RAS) and a mouse cell line (RN2) were used in this study. All these cell lines stably express Cas9 (refs. [11,28,29]) and are confirmed to be mycoplasma negative. Details in Supplementary Methods.
Proliferation assay
Cells were transduced with gRNA-expressing lentivirus and the percentage of BFP-positive (knock-out) cells was determined every 2 days between day 4 and day 14 using flow cytometry. gRNA sequences used were listed in Table S1. Details in Supplementary Methods.
Differentiation and apoptosis assays
Cells were transduced with gRNA-expressing lentivirus and BFP-positive (knock-out) cells were collected by cell sorting. For differentiation analysis, cells were stained with APC-conjugated CD11b (eBioscience) 7 days post transduction. For apoptosis assays, cells were analysed 9 days post transduction using Annexin V Apoptosis Detection Kit APC (eBioscience). Details in Supplementary Methods.
Generation of Auxin Inducible Degron (AID) KAT7 protein degradation system and treatment with indole-3-acetic acid (IAA)
MOLM-13 cells were used to generate cells expressing KAT7-AID. The resulting cells were treated with 500 μM IAA for the duration specified for each assay. Details in Supplementary Methods.
Western blot analysis
Western blot analysis was performed according to the manufacturer’s instructions. Antibodies used were listed in Table S2. Details in Supplementary Methods.
In vivo mouse work
MOLM-13-Cas9-luc cells [11] were transduced with gRNA-expressing lentivirus and BFP-positive cells were sorted 2 days post transduction. On the following day, 5 x 105 cells were injected into NSGW41 male mice. All animal studies were carried out at the Wellcome Sanger Institute under UK Home Office License PBF095404. Details in Supplementary Methods.
ChIP-seq and ChIP-qPCR
MOLM-13, OCI-AML3, THP-1 and OCI-AML2 were used for ChIP-seq and/or ChIP-qPCR analysis using antibodies listed in Table S2 and qPCR primers listed in Table S3. Details in Supplementary Methods.
RNA extraction and RNA-seq processing
RNA was extracted from AML cells with RNeasy Plus Mini Kit (Qiagen) according to manufacturer’s instructions and sequenced on Illumina HiSeq v4 platform with 75-bp paired-end sequencing. Data processing and analysis were detailed in Supplementary Methods.
Statistical analysis
Normal distribution was first tested using F-test. Student’s t-test was used for statistical testing unless stated otherwise. Mean was calculated from at least three replicates, as indicated in each figure, and error bars represent the standard deviation. P values ≤0.05 were considered statistically significant.
Results
Loss of KAT7 induces myeloid differentiation and apoptosis in AML cell lines with MLL fusion oncoprotein
We previously catalogued fitness genes in 5 AML cell lines (MOLM-13, MV4-11, OCI-AML2, OCI-AML3 and HL-60) by genome-scale CRISPR-KO screening and characterized promising therapeutic targets including the histone lysine acetyltransferase KAT2A and SRPK1 [11,13]. We further looked into genes encoding histone modifying enzymes and found that several histone acetyl transferases (HATs) were required for AML cell proliferation (Fig. S1A). Notably, two of the MYST family HATs, KAT6A and KAT7, were essential for AML cell lines driven by MLL-X fusion genes: MOLM-13 (driven by MLL-AF9), MV4-11 (MLL-AF4) and OCI-AML2 (MLL-AF6). In order to validate these fitness defects, we first investigated the effects on proliferation in MOLM-13 and MV4-11 by gRNA-mediated KO and found that KAT7 loss exhibited a consistent strong anti-proliferative effect in both cell lines (Fig. S1B). We then tested the impact of genetic disruption of KAT7 on proliferation of all cell lines used in our CRISPR-KO screens and 3 additional human leukaemia cell lines, and found that, consistent with MOLM-13 and MV4-11, KAT7 loss decreased the proliferation of other MLL-X cell lines, namely OCI-AML2, THP-1 (MLL-AF9) and Nomo-1 (MLL-AF9), but not of MLL-WT cell lines such as OCI-AML3, HL-60 and K562 (Figs. 1A, S1C,D). KAT7 dependency of MLL-X leukemic cells was also confirmed in an independently generated THP-1 Cas9 line [28], human CD34+ cells transformed by MLL-AF9 (refs. [28,30]) and a mouse AML cell line bearing MLL-AF9 (refs. [6,29]) (Fig. S1E-G). The proliferation phenotype was associated with increased expression of the myeloid differentiation marker CD11b in MOLM-13 and OCI-AML2 (Figs. 1B, S1H). KAT7 loss also led to an increased apoptosis in MLL-X cell lines (Fig. 1C). Collectively, these results indicate that KAT7 is required for in vitro survival and proliferation of leukemic cells driven by MLL-X fusion oncogenes. To assess if KAT7 inactivation has also anti-proliferative effects in vivo, we injected MOLM-13 cells previously transduced with a KAT7 gRNA (KAT7-KO) into NSGW41 male mice. Mice transplanted with KAT7-KO MOLM-13 cells showed significantly slower AML progression and increased survival, compared to those injected with MOLM-13 expressing a control gRNA (Figs. 1D, S1I), indicating the requirement of KAT7 for in vivo proliferation.
Figure 1
Loss of KAT7 exhibits anti-leukemic effects in vitro and in vivo
A) Proliferation of KAT7 KO using two sgRNAs or empty control in MLL-X AML cell lines. Percentages of BFP-positive (KAT7-KO) cells were assayed at the indicated time point and were normalized to day 4. Data are shown as mean ± SD (n= 3). B) CD11b staining in WT (empty gRNA) and KAT7 KO (gKAT7-A10) cells 7 days post transduction. C) Annexin V/PI staining of WT (Empty) and KAT7 KO (gKAT7-A10) cells 9 days post transduction. Data are shown as mean ± SD (n= 3). Two-tailed t-test was performed (N.S., not significant; **, P < 0.01; ***, P < 0.001). D) Xenograft analysis of KAT7-KO MOLM-13. Kaplan-Meier survival analysis were performed (n= 5 in each arm).
KAT7 regulates global H3K14 and H4K12 acetylation
KAT7 is known to form two distinct HAT complexes that acetylate either H3 or H4 (ref. [23]). We therefore investigated the acetylation status of H3 and H4 tails in KAT7 KO MOLM-13 cells. Among the potential KAT7 acetylation sites, acetylation of Lysine-14 residue of H3 (K3K14ac) and Lysine-12 of H4 (H4K12ac) were completely abolished upon KAT7 KO, whereas other acetylation sites were not affected (Fig. 2A). We additionally investigated major histone marks for active transcription, namely H3K4me3, H3K9ac and H3K27ac, and found that KAT7 loss did not affect the global level of these marks (Fig. S2A). Interestingly, global loss of H3K14ac and H4K12ac upon KAT7 loss was observed in all AML cell lines tested, irrespective of the MLL translocation status or proliferation phenotype (Fig. S2B).
Figure 2
Catalytic activity of KAT7 is required for leukemic maintenance.
A) Western blot analysis of potential KAT7 acetylation sites on H3 and H4 in MOLM-13 cells. B) Western blot analysis of the KAT7 acetylation sites in MOLM-13 expressing exogenous WT-KAT7 or HAT-dead KAT7 (E508Q) with or without endogenous KAT7 KO by lentiviral gRNA expression. C-D) Proliferation (C) and CD11b staining (D) of MOLM-13 expressing exogenous KAT7 WT or E508Q mutant with endogenous KAT7 disrupted by gKAT7. Data are shown as mean ± SD (n= 3).
Acetyltransferase activity of KAT7 is essential for leukemic maintenance
Since both MLL-X and MLL-WT AML cell lines showed globally reduced H3K14 and H4K12 acetylation upon KAT7 KO but only MLL-X AML cell lines showed reduced proliferation, we considered that non-catalytic roles of KAT7 complexes may have played a role in the observed proliferative defect. To rule out such possibility, we performed cDNA reconstitution experiments using gRNA-resistant WT and HAT-dead KAT7 cDNAs. For the HAT-dead KAT7 mutant, we utilized a previously characterized E508Q MYST domain mutant [31]. To replace endogenous KAT7 with a cDNA-derived protein, we first expressed a gRNA-resistant KAT7 cDNA and then knocked out the endogenous KAT7 by lentiviral transduction of KAT7 gRNA. Western blot analyses showed that control (GFP only) MOLM-13 cells lost KAT7 as well as acetylation marks on both H3K14 and H4K12 upon the endogenous KAT7 KO, whereas cells carrying the WT KAT7 cDNA retained KAT7 protein and both acetylation marks, even after gRNA-targeting of the endogenous KAT7 gene (Fig. 2B), indicating the appropriate reconstitution of KAT7 activity. In contrast, cells expressing HAT-dead E508Q mutant KAT7 protein completely lost both acetylation marks after endogenous KAT7 KO (Fig. 2B). We noted that cells expressing both endogenous KAT7 and exogenous HAT-dead KAT7 showed considerably reduced H4K12ac, although H3K14ac was maintained at the level similar to that of the control cells. This suggest a dominant negative effect of HAT-dead KAT7 protein; however, proliferation of MOLM-13 cells expressing either exogenous KAT7 was not significantly different (Fig. S2C). We then analyzed proliferation of KAT7-reconstituted cells and found that KAT7-WT-expressing cells did not show a proliferation defect, ruling out the possibility that off-target effects of KAT7 gRNAs might have been responsible for the observed phenotype (Fig. 2C, centre). In sharp contrast, the cells expressing HAT-dead KAT7 showed proliferation defect(Fig. 2C, right), similar to that in endogenous KAT7 KO cells (Fig. 2C, left). This proliferation defect was associated with increased CD11b expression (Fig. 2D). Thus, the cells expressing HAT-dead KAT7 phenocopied the endogenous KAT7 KO cells in histone acetylation, proliferation and differentiation. These results strongly indicate that KAT7-mediated acetylation plays crucial roles in the maintenance of the leukemic program driven by the MLL-X fusion protein.
KAT7 loss leads to down-regulation of a specific set of MLL-AF9 target genes
In order to analyze transcriptomic changes upon KAT7 KO, we performed RNA-seq analyses of cells transduced with lentivirus expressing KAT7 or control gRNA (Supplementary Dataset 1). We chose MOLM-13 and OCI-AML3 as the MLL-X positive and negative cell lines, respectively, in this analysis. On day 3 after transduction, we observed a total of 544 differentially expressed (DE) genes (317 up- and 227 down-regulation) in MOLM-13 (Fig. 3A). The number of DE genes increased further on day 5 to 4093 (1954 up- and 2139 down-regulation). In contrast, the KAT7-independent OCI-AML3 showed only 11 genes on day 3 and 501 on day 5. Consistent with myeloid differentiation observed by flow cytometry (CD11b upregulation, Fig. 1B), gene set enrichment analysis (GSEA) revealed that a gene set consisting of genes upregulated upon myeloid differentiation was significantly enriched on days 3 and 5 (Fig. 3B, left). Concomitantly, the MLL-AF9 target gene set was significantly depleted on days 3 and 5 (Fig. 3B, centre).
Figure 3
Transcriptomic profiling of KAT7 KO in MOLM-13 (MLL-AF9), OCI-AML3 (MLL WT).
A) Numbers of differentially expressed (DE) genes in MOLM-13 and OCI-AML3 on day 3 and day 5 after gRNA-mediated KAT7 KO. B) GSEA on day 3 (top) and day 5 (bottom) transcriptomes in MOLM-13. Gene sets consisting of genes upregulated upon myeloid differentiation (left), MLL-AF9 target genes (centre), or MLL-AF9 spreading genes (right; Table S4) were used. C) MA plot of day 3 MOLM-13 comparing KAT7 WT and KO. Black, DE genes (adj.P<0.05); yellow, MLL-AF9 spreading genes. Key MLL-AF9 spreading genes, namely MEIS1, PBX3, JMJD1C, SENP6 and MEF2C were highlighted.
Through inspection of differentially expressed genes in MOLM-13 on day 3, we identified that a small fraction of the MLL-AF9 target genes, most notably MEIS1, PBX3, JMJD1C, SENP6 and MEF2C, were significantly down-regulated (Fig. 3C). Together with HOXA9, MEIS1 and PBX3 are known to form a complex and promote transcription of HOXA target genes [32]. JMJD1C and MEF2C are also known to play an important role in MLL-AF9-induced leukemogenesis [28,33]. In addition, these genes have recently been shown to be among genes, of which MLL-AF9 binding is not limited to their promoter region but spreads into their gene body [34]. These MLL-X oncoprotein “spreading” genes are collectively more susceptible to DOT1L inhibition than MLL-X “non-spreading” genes [34], indicating that the spreading genes comprise the core downstream network of the MLL-X-driven leukemic program. We further performed GSEA on the expression profiles with a gene set consisting of protein-coding KAT7-bound MLL-AF9 spreading genes (see below, Table S4) and found significant enrichment (Fig. 3B, right). These enrichment was not observed in the MLL-WT, KAT7-independent OCI-AML3 cell line.In order to more accurately capture transcriptomic changes immediately after KAT7 loss, we employed the auxin-inducible degron (AID) system [35]. For this, we first stably expressed KAT7-AID and then knocked out the endogenous KAT7 with gRNA. We confirmed that proliferation of MOLM-13 KAT7-AID was similar to that of parental MOLM-13 cells (Fig. S3A), indicating that AID-tagged KAT7 is functional. We then introduced OsTIR1 and treated cells with auxin indole-3-acetic acid (IAA) to deplete KAT7 protein. The KAT7 protein level rapidly decreased and plateaued 2 h after treatment (Fig. 4A). This was associated with complete loss of H3K14ac over the same time course (Fig. 4A). Cells treated with IAA showed decreased proliferation, increased levels of apoptosis and myeloid differentiation (CD11b) at 48 h (Figs. 4B, S3B,C). Therefore, the phenotype after IAA treatment faithfully mirrored those seen with CRISPR-mediated KAT7 KO (Figs. 1A-C and 2A), indicating that the KAT7-AID system can be used to study the direct effects of KAT7 loss.
Figure 4
Transcriptomic profiling of the acute phase of KAT7 depletion in MOLM-13 using the degron system.
A) Time-course western blot analysis of KAT7 and H3K14ac in IAA-treated MOLM-13 KAT7-AID cells. B) Cell proliferation of IAA-treated KAT7-AID and parental MOLM-13. Data are shown as mean ± S.D. (n= 3). Two-tailed t-test was performed (***, P ≤ 0.001). C) GSEA on transcriptome of MOLM-13 KAT7-AID cells treated with IAA for 2- (top) or 4 h (bottom). The same genesets as Fig. 3B were used.
We performed RNA-seq analysis using this MOLM-13 KAT7-AID line (Supplementary Dataset 1) and subsequently GSEA. After both 2- and 4-h of treatment with IAA, there was no significant myeloid gene upregulation (Fig. 4C, left). In contrast, enrichment of MLL-AF9 spreading genes was detectable as early as 2 h, which preceded enrichment of the entire MLL-AF9 target gene set (Fig. 4C, centre and right). These results indicate that KAT7 loss disrupts a leukemia maintenance program through down-regulation of the MLL-AF9 spreading genes and subsequently leads to myeloid differentiation.
KAT7 binds to promoters of active genes, especially a subset of MLL-AF9 spreading genes
We next investigated the genomic location of KAT7 binding sites using chromatin-immunoprecipitation followed by sequencing (ChIP-seq) analysis (Supplementary Dataset 2). We again chose MOLM-13 and OCI-AML3 for this analysis and identified 9214 and 60153 KAT7 peaks, respectively. KAT7 was most enriched at promoters in both cell lines; 63.3 % of the KAT7 signals in MOLM-13 and 48.5 % in OCI-AML3 were found in promoter regions (±2 kb from TSSs) (Figs. 5A,B, S4A). Consistent with previously reported ChIP-seq studies [16,23,24], KAT7 occupancy in the AML cell lines is significantly correlated with the expression levels of target genes (Figs. 5A,B, S4B). We then focused our analysis on the MLL-AF9 spreading gene set. These genes are expressed not only in the MLL-AF9-bearing MOLM-13 but also in MLL-WT OCI-AML3 with a significant correlation of their expression levels between the two cell lines (Fig. S4C, r2=0.80). As KAT7 marks expressed genes, the majority of the MLL-AF9 spreading genes were bound by KAT7 in both cell lines (Fig. S4D). We then compared expression levels and KAT7 promoter occupancy of the spreading genes (Fig. S4E) and found that, consistent with the genome-wide patterns, OCI-AML3 showed good correlation between the two parameters (Fig. S4E, right); however, MOLM-13 did not show clear correlation, but KAT7 rather seemed to be highly bound to a small fraction of the MLL-AF9 spreading genes (Fig. S4E, left). We further compared KAT7 promoter occupancy with expression fold changes on day 3 after KAT7 KO and found that in MOLM-13 genes with higher KAT7 occupancy were more susceptible to and down-regulated upon KAT7 loss (Fig. 5C). In particular, MEIS1, PBX3, JMJD1C and SENP6, which were significantly down-regulated on day 3 (Fig. 3C), showed particularly high KAT7 promoter occupancy (Fig. 5C). There were significant differences in KAT7 promoter occupancy between the down-regulated and the up-regulated (n=4) or not-differentially expressed (n=126) genes (Fig. S4F), and KAT7 was markedly enriched in promoters of the KAT7-responsive MLL-AF9 spreading genes (Fig. S4G). In sharp contrast and consistent with its KAT7 independence, in OCI-AML3 KAT7 loss did not show a significant impact on expression on day 3, even for genes with KAT7 promoter occupancy comparable to those down-regulated in MOLM-13 (Figs. 5D, S4F). In both cell lines, for non-MLL-AF9 target genes, there was no significant difference in KAT7 promoter occupancy between up-regulated, down-regulated and not-differentially expressed genes (Fig. S5A-D). Taken together, these results suggest that although KAT7 is not required for general transcription in AML, KAT7 and its acetylase activity are required for MLL-AF9 oncoprotein to sufficiently up-regulate and maintain expression of a specific set of the MLL-AF9 spreading genes.
Figure 5
KAT7 binds to and is required for expression of a subset of MLL-AF9 targets.
A-B) Average occupancy of KAT7 per gene per base at the TSS of highly active genes (>10 FPKM), active genes (>0, ≤10 FPKM) and inactive genes (FPKM=0) in MOLM-13 (A) and OCI-AML3 (B). C-D) Comparison between KAT7 promoter occupancy and gene expression changes 3 days after KAT7 KO in MOLM-13 (C) and OCI-AML3 (D) for MLL-AF9 spreading genes. DE genes are highlighted in blue (up) or red (down).
MLL-fusion associated machineries are evicted from their target loci upon KAT7 loss
To study the consequence of KAT7 loss at the chromatin level, we utilized the KAT7-degron system (Figs. 4, S3) and performed ChIP followed by quantitative PCR (ChIP-qPCR) analysis. We first investigated whether KAT7 is required for MLL-AF9 fusion protein to maintain its binding to the promoters of key MLL-AF9 spreading genes, namely JMJD1C, SENP6, PBX3 and MEIS1, which were down-regulated significantly on day 3 after KAT7 KO (Fig. 3C) and bound by KAT7 with markedly high occupancy (Fig. 5C). ChIP-qPCR analysis using antibodies against MLL N-terminal and AF9 C-terminal regions 24 h after IAA treatment showed no significant changes in their occupancy at the promoter region (Fig. S6A,B). This observation is consistent with previous findings that MLL-AF9 recruitment is dependent on interactions with MENIN [36], LEDGF [37] and the PAF complex [38].We next investigated the binding of the histone acetylation “reader” BRD4 at gene promoters. BRD4 is a key interactor of MLL-fusion proteins and a target of anti-leukemic drugs [5,6]. Globally, BRD4 occupancy at promoter regions is highly correlated with KAT7 occupancy (Fig. S6C). At the selected key MLL-AF9 spreading genes, KAT7, MLL-AF9 and BRD4 highly co-localized (Fig. S6D). Upon KAT7 depletion, BRD4 markedly dissociated from the promoters of the key spreading genes (Fig. 6A). This was associated with dissociation of AFF1, a scaffold protein of the super elongation complex (SEC), from these promoters (Fig. 6B). In addition, RNA polymerase II Serine 5 phosphorylation (PolII-pS5) occupancy was also dramatically reduced at these loci (Fig. 6C). BRD4 dissociation from the promoters was also confirmed in THP-1 and OCI-AML2 (Fig. S6E,F). These results suggest that KAT7-mediated histone acetylation serves as a scaffold for BRD4 binding to chromatin, which is required for sustained recruitment of transcriptional activators (Fig. S7).
Figure 6
KAT7-dependent recruitment of BRD4 and SEC complex to a subset of MLL-AF9 spreading genes.
A-C) ChIP-qPCR analysis of BRD4 (A), AFF1 (B) and pS5 Pol II (C) after 24 hours of IAA treatment at the promoter region of the key MLL-AF9 spreading genes. Data are shown as mean ± S.D. (n= 3). Two-tailed t-test was performed (*, P ≤ 0.05; **, P ≤ 0.01; ***, P ≤ 0.001).
Discussion
Despite advances in understanding its genomics and molecular pathogenesis, AML remains lethal to the majority of sufferers [39] and anti-AML mainstream therapies have not changed significantly for several decades. Among different AML subtypes, cases driven by MLL-X fusion genes continue to represent an intermediate and poor prognostic category [1-3] and despite recent developments in the field [5,6,8,10,40] clinical progress is still lacking, emphasizing the need for new therapies. Here, we demonstrate that KAT7 represents a promising novel therapeutic target for MLL-X AML and provide insights into its function in the maintenance of these leukemias.KAT7-containing complexes have been shown to acetylates lysine residues of histone H3 and H4 tails. Specificity for the target histone is determined by the scaffold subunits, BRPF and JADE. In in vitro HAT assays using nucleosomes, BRPF-containing complex acetylates H3K14 and K23, while JADE-containing complex acetylates H4K5, K8 and K12 (ref. [23]). These scaffold subunits can also form a HAT complex with KAT6A and KAT6B and show a similar specificity against H3 and H4 (ref. [23]). In principle, this complex formation is interchangeable, but seems to be differentially regulated in different cell types. For example, in HeLa cells, purified KAT7-containing HAT complex contained JADE but not BRPF [27], and BRPF-containing HAT complexes predominantly contain KAT6A/B [23]. In keeping with this, knock-down of KAT7 in HeLa cells results in reduction of all of the three acetylation marks on H4 (ref. [27]). In contrast, KAT7-KO in mouse embryonic fibroblasts and KAT7 knock-down in erythroblasts led to reduction specifically of H3K14ac among the 5 potential acetylation sites and had no effect on histone H4 (ref. [24,25]). In the present study, we found that KAT7 KO in AML cells resulted in complete loss of both H3K14ac and H4K12ac, a pattern that has not been described before. This strongly suggests that KAT7 is solely responsible for acetylation of these 2 sites and that KAT6A does not compensate for KAT7 loss. It has been shown that in mouse embryos KAT6A plays a specific role in regulating local H3K9 acetylation at the Hoxa and Hoxb loci and is not involved in global H3K9 and H3K14 acetylation [41]. This might be the case in AML cells, indicating that KAT6A and KAT7 play non-overlapping roles in histone acetylation. KAT6A also showed fitness defects in MLL-X AML cell lines in our analysis. It would be interesting to investigate the molecular basis of KAT6A dependency and the mechanisms by which cell-type-specific complex formation is regulated.Besides the general correlation between expression level and KAT7 promoter occupancy, we found that KAT7 strongly bound to a small subset of MLL-AF9 spreading genes in MOLM-13. The expression of this subset is especially susceptible to KAT7 loss and it includes well-known leukemia maintenance genes in MLL-X-driven AML such as MEIS1, PBX3, JMJD1C and MEF2C. MLL-X target genes are known to be associated with high level of H3K79me2, which is deposited by histone methyltransferase DOT1L [8,40]. In MLL-X spreading genes, H3K79me2 marks also spread into gene bodies, showing broader peaks [34]. These spreading genes were shown to be more sensitive to DOT1L inhibition than the other genes (non-spreading and no-bound genes); the entire set of the spreading genes are downregulated upon DOT1L inhibition [34]. This pattern of down-regulation is in sharp contrast to the case of KAT7 loss, in which only a subset of the spreading genes that are highly bound by KAT7 are down-regulated. Therefore, recruitment of KAT7 complex to promoters of this subset might be regulated differently from the rest of the promoters.As a potential reader of H3K14ac and/or H4K12ac by KAT7, we investigated BRD4 and showed that BRD4 dissociated from the promoter of the key spreading genes upon KAT7 loss. BRD4 is known to play crucial roles in maintaining MLL-X-driven leukemic programs [5,6]. A prominent role of BRD4 in AML is to activate super-enhancers at the MYC locus and maintain high MYC expression [42]. Bromodomain inhibitors such as JQ1 and i-BET151 disrupt the interaction between BRD4 and acetylated histones, and evict BRD4 from chromatin, thereby leading to rapid downregulation of MYC. In addition to the role at enhancers, BRD4 is known to play a role in promoting Pol II elongation at promoter-proximal regions by recruiting SEC including P-TEFb [42]. Accordingly, we observed the reduced binding of AFF1 scaffold protein of SEC upon KAT7 loss. This suggests that Pol II elongation is severely affected, causing downregulation of target genes. However, since the role of BRD4 at promoters is generic, it is still unclear why only a small subset of MLL-AF9 spreading genes are affected by KAT7 loss in MOLM-13 and why gene expression is relatively stable, even in the absence of KAT7, in the KAT7-independent OCI-AML3 cells. It is plausible that additional acetylated histone readers are involved in KAT7-dependent transcriptional activation for the MLL-AF9 spreading genes. From this perspective, it is interesting to observe that pS5-Pol II reduced its occupancy at the selected promoters. Serine 5 of Pol II is phosphorylated by TFIIH after Pol II loading onto promoters [43]. Subsequent Serine 2 phosphorylation on the C-terminal domain of Pol II by P-TEFb promotes Pol II elongation [43]. Therefore, KAT7-mediated histone acetylation may also be required for the initiation phase of Pol II transcription. Further studies are required to elucidate the molecular basis of KAT7-dependent transcription of MLL-AF9 target genes.Interestingly, a study by Sauer et al. [44] reported that KAT7 was downregulated in patient AML samples, yet when knocked down by shRNA, increased cell proliferation of leukemic cell lines. Analysis using a recent large-scale RNA-seq dataset [45] revealed no significant difference between CD34+ cells and AML cells (Fig. S8). Moreover, amongst 23 TCGA cancer types with matched normal samples from the same tissue as the cancer, there were only 2 cancer types, which showed significant differences in KAT7 expression between them, namely cholangiocarcinoma (KAT7 high in tumour) and chromophobe renal cell carcinoma (KAT7 low in tumour). These observations suggest no significant correlation between KAT7 expression levels and neoplasia or cell proliferation. In addition, during the revision of this paper, Dawson and colleagues reported that KAT7 is required for the maintenance of leukemia stem cells [46]. Their conclusion is consistent with ours and further highlights the importance of KAT7 function in leukemia maintenance.Collectively, our findings reveal that KAT7 acts upstream of BRD4 to recruit MLL-fusion associated machineries to the promoter of key MLL-AF9 target genes, which represent an alternative therapeutic target for MLL-X leukemia. We anticipate that our work will motivate the development of small-molecule inhibitors that are specific to KAT7 and highlight KAT7 as a potential therapeutic target for MLL-fusion AML.
Authors: Jolanta Grembecka; Shihan He; Aibin Shi; Trupta Purohit; Andrew G Muntean; Roderick J Sorenson; Hollis D Showalter; Marcelo J Murai; Amalia M Belcher; Thomas Hartley; Jay L Hess; Tomasz Cierpicki Journal: Nat Chem Biol Date: 2012-01-29 Impact factor: 15.040
Authors: Elli Papaemmanuil; Moritz Gerstung; Hartmut Döhner; Peter J Campbell; Lars Bullinger; Verena I Gaidzik; Peter Paschka; Nicola D Roberts; Nicola E Potter; Michael Heuser; Felicitas Thol; Niccolo Bolli; Gunes Gundem; Peter Van Loo; Inigo Martincorena; Peter Ganly; Laura Mudie; Stuart McLaren; Sarah O'Meara; Keiran Raine; David R Jones; Jon W Teague; Adam P Butler; Mel F Greaves; Arnold Ganser; Konstanze Döhner; Richard F Schlenk Journal: N Engl J Med Date: 2016-06-09 Impact factor: 91.245
Authors: Scott R Daigle; Edward J Olhava; Carly A Therkelsen; Aravind Basavapathruni; Lei Jin; P Ann Boriack-Sjodin; Christina J Allain; Christine R Klaus; Alejandra Raimondi; Margaret Porter Scott; Nigel J Waters; Richard Chesworth; Mikel P Moyer; Robert A Copeland; Victoria M Richon; Roy M Pollock Journal: Blood Date: 2013-06-25 Impact factor: 22.113
Authors: Scott R Daigle; Edward J Olhava; Carly A Therkelsen; Christina R Majer; Christopher J Sneeringer; Jeffrey Song; L Danielle Johnston; Margaret Porter Scott; Jesse J Smith; Yonghong Xiao; Lei Jin; Kevin W Kuntz; Richard Chesworth; Mikel P Moyer; Kathrin M Bernt; Jen-Chieh Tseng; Andrew L Kung; Scott A Armstrong; Robert A Copeland; Victoria M Richon; Roy M Pollock Journal: Cancer Cell Date: 2011-07-12 Impact factor: 31.743
Authors: Johannes Zuber; Junwei Shi; Eric Wang; Amy R Rappaport; Harald Herrmann; Edward A Sison; Daniel Magoon; Jun Qi; Katharina Blatt; Mark Wunderlich; Meredith J Taylor; Christopher Johns; Agustin Chicas; James C Mulloy; Scott C Kogan; Patrick Brown; Peter Valent; James E Bradner; Scott W Lowe; Christopher R Vakoc Journal: Nature Date: 2011-08-03 Impact factor: 49.962
Authors: Konstantinos Tzelepis; Hiroko Koike-Yusa; Etienne De Braekeleer; Yilong Li; Emmanouil Metzakopian; Oliver M Dovey; Annalisa Mupo; Vera Grinkevich; Meng Li; Milena Mazan; Malgorzata Gozdecka; Shuhei Ohnishi; Jonathan Cooper; Miten Patel; Thomas McKerrell; Bin Chen; Ana Filipa Domingues; Paolo Gallipoli; Sarah Teichmann; Hannes Ponstingl; Ultan McDermott; Julio Saez-Rodriguez; Brian J P Huntly; Francesco Iorio; Cristina Pina; George S Vassiliou; Kosuke Yusa Journal: Cell Rep Date: 2016-10-18 Impact factor: 9.423
Authors: Maheswaran Solayappan; Adam Azlan; Kang Zi Khor; Mot Yee Yik; Matiullah Khan; Narazah Mohd Yusoff; Emmanuel Jairaj Moses Journal: Front Genet Date: 2022-01-28 Impact factor: 4.599