SET is over-expressed in KMT2A-R cell lines and primary samples
To characterise the role of SET across various hematopoietic contexts, we used public data repositories to analyse its mRNA expression in human HSCs and progenitors (n=34) as well as in multiple independent human AML primary samples (n=384) covering the main cytogenetic subsets (10, 26). SET was expressed at high levels in both HSC and myeloid progenitors (MP) compared with mature monocytes and myelocytes (Fig. 1A and supplementary Figure 1A), thereby indicating that it is a gene expressed during early hematopoietic development. In silico evaluation of SET mRNA expression across all AML was uniformly high and it did not show segregation with any distinct molecular group as previously suggested (Fig. 1B and supplementary Fig 1A-C) (19, 20). A PrognoScan database-based Kaplan–Meier analysis of the overall survival of 117 AML patients by high (n=11) and low (n=106) SET levels (28), revealed that high SET expression positively correlated with poor overall survival in human AML (Fig. 1C), consistent with previous reports (20). To further substantiate the role of SET in AML we evaluated its protein levels by western blot in a panel of KMT2A-R-AML (THP1, MV411, ML2, MOLM13, NOMO-1) and KMT2A-R-ALL (SEM, Hb1119, KOPN8, RS411) cell lines and primary samples (PS), KMT2A wild-type (wt) cell lines (K562, BCR::ABL+ erythroleukemia, Kasumi1 AML1::ETO+ AML, REH TEL::AML+ ALL, U937 CALM::AF10+ AML), mononuclear cells isolated from the bone marrow (BM) and peripheral blood (PB) of healthy adult volunteers. SET was significantly up-regulated in all the leukemic cells lines and in all KMT2A-R primary samples, irrespective of tumor lineage, compared to BM controls (Fig. 1D-E). Given the poor outcomes of the KMT2A-R-leukemias and the relatively well characterised cellular context driving leukemogenesis, we further investigated the mechanistic role of SET in this group. As SET oncoprotein has several distinct roles depending on its subcellular localisation (22), we analysed its subcellular localisation in two KMT2A-R-cell lines (THP1 and MV411) and one KMT2A-wt cell line (K562) by nuclear/cytoplasm fractionation followed by western blot and showed that SET was relatively more abundant in the cytoplasm than in the nucleus of these cells (Fig. 1F). As phosphorylation of SER9 and SER25 of SET inhibits its nuclear import (29, 30), we also investigated the phosphorylation status of SET. As no specific antibodies against phospho-SET are commercially available, SET was immunoprecipitated and protein samples were analysed by western blot with anti-pSER antibodies. Consistent with the prevalent cytoplasmatic localization, our results showed that, in KMT2A-R-AML cell lines, SET is phosphorylated on Serine residues (Fig. 1G). Overall, these data indicate that SET is over-expressed, phosphorylated and abundantly localised in the cytosol in KMT2A-R and -wt cells.
SET is essential for KMT2A-R leukemia self-renewal
To understand the role of SET in leukemic stem cells (LSC), we analysed the expression of SET mRNA in 12 high LSC frequency mouse KMT2A-R AMLs (KMT2A::MLLT3 and KMT2A:MLLT1) and 22 low LSC frequency KMT2A-R-AMLs (KMT2A::AFF1p, KMT2A::AF10 and KMT2A::GAS7), using public data repositories (32). Interestingly, we found that SET mRNA expression was significantly higher in high LSC frequency KMT2A-R AML than in low LSC frequency KMT2A-R AML, (Fig. 2A), suggesting a potential role of SET in KMT2A-R leukemia self-renewal. Therefore, to determine whether SET has a functional role in KMT2A-R- leukemia, we knocked down SET gene in human cell lines by using RNA interference. SET knockdown (KD) completely abolished the clonogenic ability of KMT2A-R leukemic cell lines (for THP1 and MV411 almost 100% reduction in colony number, for SEM 90.5% reduction) (Fig 2B, E, F and H and supplementary 2), whereas it had little or no effect on the colony formation of three independent KMT2A-wt leukemic cell lines (K562 shScramble 554± 126 vs shSET 604± 83 p>0.05; Kasumi1 shScramble 227± 54 vs shSET 167±140 p>0.05; REH shScramble 479±54 vs shSET 271±115 p>0.05) (Fig. 2B, C, D, G). SET KD attenuated the proliferation of K562 (eGFP-K562 fold ratio shSET vs shScramble 0.89 p<0.05) (Fig 2I, K and supplementary Figure 3), in line with a previous report (18) and similarly to KMT2A-wt AML cell line Kasumi1 (eGFP-Kasumi fold ratio shSET vs shScramble 0.55 p<0.0001) (Fig. 2J-K and supplementary Figure 3). Overall, these data indicate that SET KD specifically abolishes the self-renewal of KMT2A-R-leukemic cells and, consistently, it attenuates the proliferation of KMT2A-wt leukemic cells.
The SET inhibitor FTY720 induces cell cycle arrest and drives cell death in KMT2A-R leukemic cells
We next tested pharmacological modulation of SET by FTY720 (Fingolimod) (31), a FDA-approved immunosuppressive drug that has gained further attention as anti-cancer and PP2A activating drug, due to its ability to disrupt the binding between SET and PP2A (32-34). We first carried out a dose-response titration assay and determined that the half maximal inhibitory concentration of FTY720 in vitro ranged between 1 and 5 mM (Supplementary Fig. 4A), which were reported as non-toxic to healthy bone marrow mononuclear cells (23, 35). We then assessed the effect of FTY720 on proliferation, cell cycle and cell death of KMT2A-wt and KMT2A-R-leukemic cells. As in our model the GFP fluorescence is proportional to the number of alive cells, we used GFP as a reporter of cell viability (Fig. 3A). We observed that 5 µM FTY720 had a variable but significant effect on the proliferation of all the analysed cell lines, ranging from modest for the KMT2A-wt-cell lines (vehicle vs FTY720 fold ratio eGFP-K562= 1.19; eGFP-Kasumi1= 1.85) to severe for the KMT2A-R-cell lines (vehicle vs FTY720 fold ratio eGFP-THP1= 3.17; eGFP-MV411= 4.18; eGFP-SEM= 2.32). In addition, FTY720 severely halted the proliferation of KMT2A-wt eGFP-REH (vehicle vs FTY720 fold ratio= 6.17), as previously reported (36). Cell cycle analyses revealed that treatment with FTY720 for 48 hours induced a significant increase of cells in G1 and a reduction in S and G2-M phase, in KMT2A-R-cells (Fig. 3B and supplementary Fig. 4B), indicating cell cycle failure. Moreover, we investigated the faction of leukemic cells undergoing apoptosis upon FTY720 treatment by FACS, by gating the GFP negative (GFP-) cells. FTY720 induced a statistically significant increase in apoptosis in two KMT2A-wt-cell lines (fraction of GFP- cells in FTYT720 treated vs vehicle for eGFP-REH 0.9± 0.028 vs 0.016± 0.002; for eGFP-K562: 0.393±0.059 vs 0.147±0.027), whereas this was not significant for the KMT2A-wt cell line eGFP-Kasumi1 (0.124± 0.036 vs 0.147±0.062). In contrast, FTY720 induced a consistent and statistically significant increase in apoptosis in KMT2A-R-cells (eGFP-THP1:0.43±0.08 vs 0.10±0.061; eGFP-MV411: 0.285±0.136 vs 0.12±0.072; eGFP-SEM: 0.12±0.01 vs 0.04±0.004; eGFP-Hb1119: 0.391± 0.028 vs 0.038± 0.004) (Fig. 3C and supplementary Fig. 4C). These results indicate that the SET inhibitor FTY720 induces heterogenous effects in leukemic cells; specifically, in KMT2A-R-leukemic cells, FTY720 causes cell cycle arrest in G1 and significantly increases the rate of apoptosis.
The effect of FTY720 is dependent on PP2A activation
FTY720 is a SET inhibitor able to rescue the activity of PP2A towards its target pathways (31, 37). By immunoprecipitation, we showed that FTY720 treatment for 24 hours disrupted the binding between SET and PP2A in KMT2A-R cells, confirming the molecular mechanism reported in other leukemic models (Fig.4A) (23, 32-34). To investigate whether the observed effects of FTY720 on KMT2A-R cells were due to the activation of PP2A, we analyzed the phosphorylation of some of PP2A targeted pathways (38) by western blot. 24 hours after treatment we observed a reduction in the abundance of phospho-AKT1(SER473) in KMT2A-wt cells (Fig. 4B); in Kasumi1 cells, the reduction in phospho-AKT1 was sustained over 48 hours, whereas in K562, phospho-AKT1 expression returned to the levels of vehicle-treated cells, indicating that the effect of FTY720 was temporary. In K562 cells, we also observed reduction in phospho-ERK1 (THR202/TYR04), which was instead unchanged in Kasumi1; in these cell lines, also phospho- GSK3β (SER9) expression decreased. Notably, in KMT2A-R-leukemic cells, we observed a dramatic and stable decrement in the expression of the PP2A targets phospho- AKT1 and phospho- ERK1 (Fig. 4B). To confirm that the effect of FTY720 was specifically dependent on PP2A activation, we performed the same experiments by pre-treating cells with the phosphatase inhibitor okadaic acid (OA) used at a concentration that inhibits PP2A but no other phosphatases (19). We first tested OA in K562 cells and confirmed that the range of concentrations used did not impact the cell proliferation (Supplementary Fig. 5A). As expected, OA treatment caused a significant increase in the PP2A targets, phospho-GSK3β, phospho-ERK1 and phospho-AKT1, with maximal expression after treatment with 2.5nM OA for 4 hours and 5nM for 2 hours (Supplementary Fig. 5B). We therefore used the lower concentration of OA for our combination experiments with FTY720. Western blot analysis of PP2A targets in the cells pre-treated with OA for 4 hours and then treated with FTY720 for 20 hours, revealed that pre-treatment with OA restored phospho-AKT1 and phospho-ERK1 levels in all the leukemic cells analysed (Fig. 4C). Whereas pre-treatment with OA did not have any effect on the percentage of apoptotic cells in eGFP-K562 and eGFP-Kasumi1, it significantly decreased the percentage of apoptosis in KMT2A-R-cells (fraction GFP- cells FTY720+ OA vs FTY720 eGFP-THP1: 0.11 ± 0.047 vs 0.32 ± 0.090; eGFP-MV411 0.10 ± 0.039 vs 0.26 ± 0.14) (Fig. 4D and supplementary Fig. 6), suggesting that OA rescues the effect of FTY720 in KMT2A-R-leukemic cells. A similar effect was obtained by knocking down PPP2CA, the gene encoding for the catalytic a subunit of PP2A (Fig 4E-F). These data indicate that the effects of FTY720 on KMT2A-R cells are dependent on PP2A activation.
The phospho-proteome reveals FTY720-mediated effects on cell division, apoptosis and gene transcription
To reveal the global impact of FTY720 on leukemic cells, we performed a phospho-proteomic analysis on two KMT2A-R cell lines (eGFP-THP1 and eGFP-MV411) treated with FTY720, using in gel digestion and liquid chromatography–tandem mass spectrometry (LC–MS/MS) (25). LC–MS/MS analysis showed a differential pattern of phospho-protein abundance in eGFP-THP1 and eGFP-MV411 cells, with 2276 phosphosites up-regulated and 1862 phosphosites down-regulated (>2 fold and p<0.01) in eGFP-THP1, and with 1428 phosphosites up-regulated and 743 phosphosites down-regulated ( >2 fold and <0.01) in eGFP-MV411 (Fig 5A and B and Supplementary Fig. 7), in comparison to vehicle-treated cells. The complete list of identified proteins and phospho-peptides are provided in Supplementary Table 1. To gain deeper biological insights, we categorized the phospho-proteins using Gene Ontology (GO) (Fig 5C), revealing that treatment with FTY720 led to a robust decrease in cell division and an increase in apoptosis-related kinase signalling in eGFP-THP1 cells (Fig 5C and Supplementary Fig. 7E). In contrast, the data indicate a strong decrease in cell division-related kinase signalling, but a subtler increase in apoptosis eGFP-MV411 (Fig 5C and Supplementary Fig 7F). In addition, FTY720-modulated phosphosites implicated in transcription regulation, chromatin organisation, DNA damage repair, mRNA processing, microtubule and actin cytoskeletal organization (Fig 5C). We then used Kinase-Substrate Enrichment Analysis (KSEA) for the characterisation of kinase activity from the phospho-proteomic dataset. The mitosis-regulating kinase Aurora Kinase B was the kinase most significantly impacted by FTY720 in both cell lines (Fig. 5D- E). Indeed, FTY720 decreased the phosphorylation of several Aurora Kinase B targets, including PLK1 (THR210) (Fig 5 F-G-H). As expected, FTY720 inhibited MAPK3 (ERK1) and phosphorylation of ERK1 on THR202 and TYR204 was significantly reduced in both cell lines (Fig Supplementary Fig 7G); furthermore, FTY720 significantly inhibited Abl1 and CDK2 in eGFP-MV411 (Fig 5I and supplementary Figure 7I). Although the overall effect of FTY720 resulted in a significant increase in GSK3b activity, the analysis of the phosphosites regulated by this kinase indicated that this effect was mostly due to the increased phosphorylation of a single target RCAN1 (SER163), whereas GSK3b-dependent phosphorylation of MYC on THR58 was significantly inhibited by FTY720 (Fig 5J and supplementary Fig 7J). FTY720 significantly increased the activity of DNA damage kinase ATM, with an overall activation of targets implicated in DNA repair by Non homologous end joining (NHEJ), such as PRKDC (DNA-PK), and inhibition of sensors involved in transcription and DNA repair by homologous recombination, such as BRCA1 (Fig 5K-L and supplementary Fig7K). Taken together, these data indicate that FTY720 reduced the activity of phospho-signalling associated to cell division and increased apoptosis-related and DNA damage kinase signalling in KMT2A-R cells.
FTY720 affects the core transcriptome of KMT2A-R leukemia
As the phospho-proteomics indicated phosphorylation changes in targets involved in transcription (Fig 5C), we performed RNA-seq analysis to determine the gene expression profile of KMT2A-R cell line THP1 treated with FTY720 for 24 hours. Volcano plot filtering was used to identify differentially expressed genes between vehicle control group and FTY720 treated group (Fig. 6A). According to the results, 980 genes were significantly up-regulated and 898 genes were significantly downregulated by FTY720 (fold change >1.3, padj<0.05). Functional annotation clustering of these differentially expressed genes using the Gene ontology (GO), Kyoto Encyclopedia of gene and genomes (KEGG) and Reactome annotation databases indicated that sets of down-regulated genes were highly enriched in functional groups that related to ribosome biogenesis and rRNA processing (Supplementary Fig. 8A), whereas up-regulated genes were highly enriched in functional groups that related to myeloid cell activation (Supplementary Fig. 8B). In addition, upon treatment with FTY720, genes involved in autophagy, intrinsic apoptotic signalling pathway, extrinsic apoptotic signalling pathway and cell cycle arrest resulted upregulated (Supplementary Fig. 8B). Thus, these data corroborate our previous findings on cell cycle arrest and apoptosis and also indicate that FTY720 can induce cell death by multiple mechanisms, as reported in other models (36, 39-41). More importantly, several genes over-expressed in cancer, including SET and MYC were down-regulated (Supplementary table 2), as well as genes associated with histone methyltransferase activity, among which several members of the KMT2A- fusion epigenetic complex, and HOXA9/MEIS1 target genes (3, 8-10) (Supplementary table 3). We performed RT-qPCR and western blot to validate the decreased expression of SET and MYC upon FTY720 treatment (Fig. 6B and supplementary Fig. 9A-B). Interestingly, western blot analyses indicated that SET protein was significantly reduced only in KMT2A-R-cells (Fig. 6C). RT-qPCR also confirmed a specific decrease in the expression of the KMT2A target genes HOXA9 and HOXA10, in KMT2A-R-cells (Fig. 6B); these genes were also down-regulated in KMT2A-R-cells, when SET was knocked down (Fig. 6D), suggesting that SET might regulate the expression of these genes in KMT2A-R-cells. Previous reports indicated that SET and the oncofusion protein SET::NUP214 modulate the expression of HOXA gene cluster in HeLa cell line and in T-ALL primary samples (42-44). To identify whether SET was enriched on HOXA9 and HOXA10 promoters, we performed chromatin immunoprecipitation (ChIP) experiments. As expected, KMT2A localized on the promoters of both HOXA9 and HOXA10; in contrast, SET was enriched only on HOXA10 promoter (Fig. 6E and supplementary Fig 9C and 10). Immunoprecipitation experiments revealed that SET interacted with both KMT2A and KMT2A-fusion protein (Fig. 6F and supplementary Fig. 9D). Collectively, our results indicate that genetic and pharmacological modulation of SET rewires the KMT2A-R- core gene expression signature and reduces the expression of key genes critical for sustaining this disease.
FTY720 leads to increased sensitivity of KMT2A-R leukemia to chemotherapy
Daunorubicin is an anthracycline included in the intensive multiagent chemotherapy to induce remission in KMT2A-R-AML patients (1). We investigated whether SET targeting via FTY720 could enhance daunorubicin-induced cytotoxicity in KMT2A-R-cells. To this aim, we tested 5µM FTY720 with 10nM daunorubicin, a concentration that had shown a very modest effect on the proliferation and survival of AML cells, including those carrying KMT2A-translocations. Whereas the combination treatment FTY720 + Daunorubicin did not have any effect on eGFP-K562, the percentage of apoptotic cells was significantly high for all the AML cell lines (Figure 7A). The apoptosis increase was modest only for eGFP-Kasumi1 cell line (fraction GFP- cells in FTY720 vs FTY720 +DNR 0.12±0.06 vs 0.22± 0.082), but the effect was very strong in KMT2A-R-cell lines, namely eGFP-THP1 (0.25±0.13 vs 0.55± 0.23, p<0.0001) and eGFP-MV411 (0.24 ±0.18 vs 0.57±0.19, p<0.001) (Fig. 7A and supplementary Fig. 11). To assess the biological importance and therapeutic relevance of SET targeting via FTY720 in KMT2A-R-leukemia, we tested FTY720 in combination with daunorubicin in two KMT2A-R-patient-derived-xenograft (PDX) models in vitro. Combination treatment resulted in a significant reduction of colony number compared with vehicle and with single drug treatment of all samples (74% reduction in comparison to vehicle) (Fig. 7B-C), suggesting that FTY720 treatment enhances the response to daunorubicin.