BCR::ABL1-induced malignant transformation and transcriptional reprogramming is associated with genome-wide enhancer activation
To investigate if BCR::ABL1 causes enhancer deregulation to establish its transcriptional program, we first monitored BCR::ABL1-induced malignant transformation and associated transcriptional reprogramming for changes in enhancer activation. In detail, we used an established murine transformation model based on viral transduction of ex vivo cultured primary B-cell precursors (BCPs) with BCR::ABL1-encoding retrovirus21,31,39,42 and monitored transcriptional changes and enhancer activation via RNA-Seq and H3K27ac ChIP-Seq. H3K27ac is a widely used histone marker of active regulatory regions including promoters and enhancers43. To differentiate between active promoters and active enhancers, we specifically analyzed H3K27ac signals at non-promoter regions. We first re-assessed RNA-Seq and H3K27ac ChIP-Seq data that we generated for a previous study39 (Suppl.Fig.1A/B). In this study, BCPs from 53BP1-/- mice were transduced with BCR::ABL1p210 (Fig.1A). Leukemic transformation was accompanied by substantial transcriptional reprogramming visualized by RNA-Seq (Fig.1B, left; Table S1). BCR::ABL1-induced genes included the known targets Xbp1, Bcl2, Etv5, Dusp6 and Ccnd2, while downregulated genes included many with B-cell specific functions as described before44. Changes in gene expression were accompanied by changes in enhancer activation (Fig.1B, middle) including increased signals at the Myc super-enhancer in BCR::ABL1 expressing cells (Fig.1B, right). We repeated our setup using BCPs from C57BL/6 wild-type mice and BCR::ABL1p190-encoding retrovirus (Suppl.Fig.1C/D; Table S2). This again showed that enhancer-deregulation occurs during BCR::ABL1-mediated transformation (Fig.1C). Of note, besides causing specific changes at defined sites, BCR::ABL1 induced a modest increase in H3K27ac, globally (Fig.1D).
Genes regulated by BCR::ABL1 and/or having critical functions in Ph+B-ALL are defined by long-range chromatin interactions with enhancers
To define which genes are enhancer-regulated in Ph+B-ALL, we performed promoter-capture Hi-C (PCHi-C)41,45, which profiles long-range chromatin interactions of promoters with distal DNA regions (i.e., other ends/OE) such as enhancers. This is important as enhancers are often distally located to the promoters they regulate, sometimes being separated by megabases of DNA, and often do not regulate the most proximally located gene3-5. We integrated PCHi-C results with H3K27ac ChIP-Seq data to specifically define interactions between active promoters and active enhancers (Fig.2A). PCHi-C and H3K27ac ChIP-Seq were performed on two Ph+B-ALL cell lines and leukemia cells of one Ph+B-ALL patient that were ex vivo cultured on OP9 feeder cells (Suppl.Fig.2A-E). We identified ~55,000 significant long-range promoter interactions per sample (CHiCAGO score >5) using probes for 18,202 protein-coding and 10,929 non-protein-coding genes. ~19% of these represented interactions of an H3K27ac+ active promoter with another H3K27ac+ distally located site (other end/OE), potentially indicative of an interaction with an enhancer (Fig.2A right, Table S3-8). In line with previous work3,46, most active promoters were connected via long-range chromatin interactions to such distally located active regions (H3K27ac+ OEs), while the opposite was true for inactive promoters (Fig.2B), and active promoters with H3K27ac+ OE interactions displayed a much higher number of total long-range chromatin interactions compared to active promoters that only interacted with H3K27ac-negative OEs (Fig.2C). Notably, a fraction of active genes with H3K27ac+ OE interactions solely interacted with other active promoters (i.e., Promoter-Promoter Interactions/PPIs; 10.7%). However, most of them displayed interactions with H3K27ac+ non-promoter regions (i.e., potential Enhancer-Promoter Interactions/EPIs; 89.3%) either alone or in combination with PPIs (Suppl.Fig.2G). Interestingly, genes with EPIs, EPIs+PPIs, or PPIs appeared to serve different functions, with EPI+ only genes largely relating to signaling, EPI+PPI+ genes relating to gene expression and RNA processing, and PPI+ only genes associating with protein processing and transport (Suppl.Fig.2F). Importantly, while in total only 43% of all active genes displayed EPIs by PCHi-C (Fig.2D), these included many genes known to be regulated by BCR::ABL1 and/or having critical functions in Ph+B-ALL (Fig.2E), with several of them interacting with enhancers with crucial functions previously reported in other cancers47-49.
To test whether the PCHi-C-defined genes with EPIs are indeed enhancer-regulated, we next inhibited the enhancer activator CBP/P30050 in Ph+B-ALL cells with sub-lethal concentrations of A-48551 and monitored resulting changes in gene expression by RNA-Seq. A-485 inhibits the histone acetylase activity of CBP/P300 reported crucial for CBP/P300-mediated enhancer activation50 and causes dissociation of interacting enhancers and promoters52 (Fig.2F) resulting in transcriptional downregulation of enhancer-regulated genes52. As such, we expected A-485 to cause preferential downregulation of PCHi-C-defined EPI+ genes compared to EPI- genes. In agreement with A-485 interfering with transcriptional activation, A-485 predominantly caused gene downregulation (~67% of DEGs; Table S9/10). Downregulated genes included EPI+ Ph+B-ALL key genes (Fig.2G and Suppl.Fig.2H) and lineage/cell type-specific genes as previously described53 (Suppl.Fig.2I). A global analysis further showed that genes with PCHi-C-defined EPIs were indeed generally enriched for A-485-induced downregulation when compared to active genes without EPIs or PPIs (Fig.2H). Genes with PPIs rather showed a trend towards upregulation in response to A-485 (Fig.2H). Genes with both EPIs and PPIs showed a trend towards downregulation (Fig.2H), but this only reached statistical significance in one of the two cell lines tested (Suppl.Fig.J). The latter might relate to the opposite effect of A-485 on PPIs described here, which complicates the response of genes with EPIs+PPIs.
Together, the described results show that the transcriptional program of Ph+B-ALL is largely enhancer-associated and that BCR::ABL1-driven malignant transformation towards Ph+B-ALL is associated with defined changes in enhancer activation.
Enhancer activation and associated enhancer-promoter interactions
define Ph+B-ALL identity
We next assessed if the changes in enhancer activity adapted by Ph+B-ALL cells are specific for this leukemia-subtype. We first assessed enhancer activation itself by comparing non-promoter H3K27ac ChIP-Seq signals present in Ph+B-ALL cell lines and primary leukemic cells to those from other leukemias and leukemia sub-types. H3K27ac signals at non-promoter regions clearly separated Ph+B-ALL from Ph-negative B-ALLs and BCR::ABL1-driven chronic myeloid leukemia (CML) by principal component analysis (PCA; Fig.3A). Likewise, we compared our PCHi-C data from Ph+B-ALL cells to PCHi-C data that we generated for healthy CD19+CD10+ BM BCPs and CML cells (Fig.3B/C). Specifically, we defined Ph+B-ALL ‘core interactions’ (i.e., EPIs and PPIs present in all three Ph+B-ALL samples; Suppl.Fig.3; Table S11) and assessed if these chromatin interactions would discriminate Ph+B-ALL cells from Ph-negative cells. Again, both, PCA and cluster analysis showed efficient separation of Ph+B-ALL cells from healthy BCPs and CML (Fig.3C). Notably, healthy BCPs clustered closer to Ph+B-ALL samples than CML and closer to primary Ph+B-ALL than to Ph+B-ALL cell lines, reflecting their similarities in lineage type and cellular state.
To validate these observations with a different method and samples, we performed HiChIP using H3K27ac antibodies54 on leukemia cells from three patient-derived xenografts (PDX) each for Ph+B-ALL and KMT2A::AFF1+ B-ALL for comparison (Fig.3D-F). Like PCHi-C, H3K27ac HiChIP visualizes EPIs but, additionally, also visualizes enhancer-enhancer interactions (EEIs). In line with the above, H3K27ac HiChIP-defined EPIs and EEIs efficiently separated Ph+B-ALL from the KMT2A::AFF1+ B-ALL cells (Fig.3E). We further compared H3K27ac HiChIP-derived EPI interaction scores per gene with RNA-Seq expression of the respective genes. Specifically, we used log2 fold change (log2FC) values from the comparison of Ph+B-ALL and KMT2A::AFF1+ B-ALL cells to assess if their separation by EPIs and enhancer activation relates to respective differences in their gene expression programs. The comparison indeed showed a trend towards positive correlation of EPIs and differential gene expression (Fig.3F).
Enhancer activity and enhancer-promoter interactions of BCR::ABL1-induced genes depend on BCR::ABL1 kinase activity
Of note, BCR::ABL1-induced transcriptional reprogramming and associated enhancer deregulation required a prolonged period of BCR::ABL1 expression in our murine model (Fig.1), potentially indicating a positive selection for changes in enhancer function during malignant transformation instead of a BCR::ABL1-induced active process. To assess if BCR::ABL1 actively regulates enhancer activation, we thus inhibited BCR::ABL1 function for 24 hours using the tyrosine kinase inhibitor, Ponatinib, (Fig.4A) and monitored its effect on enhancer activation through H3K27ac ChIP-Seq. Ponatinib caused substantial loss of H3K27ac signals at H3K27ac+ non-promoter regions in both BCR::ABL1p190-transformed murine BCPs and human Ph+B-ALL cells (Fig.4B). Ponatinib treatment also increased H3K27ac at other enhancer peaks, possibly reflecting resumed differentiation towards mature B-cells upon BCR::ABL1 inhibition as described previously44,55. Next, we specifically analyzed genes that display PCHi-C-defined EPIs (Fig.4C), focusing on genes that are downregulated upon Ponatinib treatment by RNA-Seq (Table S12-14). Inspection of the BCR::ABL1 target gene CCND2 showed that PCHi-C interaction as well as H3K27ac signals at both promoters and enhancers were reduced upon Ponatinib treatment (Fig.4D). We consecutively analyzed all genes that are downregulated upon Ponatinib for changes in enhancer-promoter interactions, H3K27ac signals at promoters, and H3K27ac signals at interacting enhancers. PCHi-C scores and H3K27ac signals at both promoters and enhancers were indeed significantly reduced upon Ponatinib treatment in all cases (Fig.4E-G).
BCR::ABL1 induces enhancer activation in part through the BCR::ABL1-induced transcription factor STAT5
To define TFs that BCR::ABL1 employs to induce enhancer activation, we next assessed two published mass spectrometry data sets of BCR::ABL1 phosphorylation targets56,57. Within TFs identified by both studies, peptides of STAT family TFs dominated, including known phospho-tyrosine (pY) sites of STAT3, 5 and 6 important for their activation (Fig.5A). For the remainder, respective pYs were either not studied yet or shown functionally redundant (e.g. for NFYAY266 ref.58). For functional analysis, we focused on STAT5B as it was reported the most important STAT in Ph+B-ALL32, most enriched in our analysis (Fig.5A), known to recruit CBP/P300 to target genes59, and linked by previous studies on murine B- and T-cells to enhancer activation60,61. Specifically, we generated a Ph+B-ALL cell line model that allows targetable proteasomal degradation (TPD) of STAT5 using the degron (dTAG) system62,63. In degron cells, a FKBP12F36V tag is added to a protein of interest, which is a ligand for dTAG compounds that recruit E3 ligases to induce proteasomal degradation of FKBP12F36V-tagged proteins. Through gene editing, we rendered both STAT5A alleles in Ph+B-ALL BV-173 cells nonfunctional and added an N-terminal FKBP12F36V tag to STAT5B alleles (Fig.5B). Inclusion of an HA-tag further allow visualization of STAT5 by HA-ChIP-Seq. In parallel, we applied the novel STAT5A/B-specific degrader, AK-229264, on parental BV-173 cells with otherwise germline STAT5A/B. Of note, while dTAGV-1 addition was used in degron cells to induce STAT5 loss, in parental cells, which do not express FKBP12F36V-tagged proteins, we used dTAGV-1 as a ‘non-specific degrader’ control (Fig.5C). With both systems, we achieved an almost complete loss of STAT5A/B (Fig.5C), though STAT5-degron cells exhibited reduced basal expression of STAT5, which was nevertheless further reduced by dTAGV-1 (Fig.5C). We first determined STAT5 binding through HA-ChIP-Seq using our STAT5-degron cells. STAT5-HA-bound sites included known canonical STAT5 targets (Fig.5D) and were highly enriched for STAT5 DNA-binding motifs (Fig.5E). In line with STAT5 being activated by BCR::ABL1 and phosphorylation being a prerequisite for DNA-binding28-30, Ponatinib abolished STAT5-HA DNA-binding almost entirely (Fig.5D/F). STAT5 binding was further enriched at promoters of BCR::ABL1-regulated genes (821/1859 [44%] of all Ponatinib DEGs were STAT5-HA bound), which also displayed increased STAT5-HA signals compared to BCR::ABL1-unregulated genes (Suppl.Fig.4A/B). STAT5 loss through AK-2292 in parental cells or dTAGV-1 in degron cells likewise predominantly affected BCR::ABL1-regulated genes (125/187 [67%] of AK-2292 DEGs overlap with Ponatinib DEGs; Fig.5G/H, Suppl.Fig.4C). In agreement with STAT5 also acting as a repressor65,66, STAT5 loss induced both transcriptional downregulation and upregulation, with AK-2292 largely mirroring Ponatinib-induced transcriptional effects (Fig.5G right). However, notably, only 9% of genes (821/8705) with STAT5 at promoters overlapped with BCR::ABL1-regulated genes (i.e., Ponatinib DEGs, Suppl.Fig.4A), and STAT5 loss through AK-2292 only affected 7% (125/1859) of BCR::ABL1-regulated genes (Fig.5G left, Table S15/16).
Next, we applied these models to investigate STAT5 on BCR::ABL1-induced enhancer activation. Indeed, 69% of STAT5-HA peaks in untreated cells (17,059/24,642) did not overlap with a promoter, supporting a role of STAT5 in enhancer regulation in Ph+B-ALL. As we did not have PCHi-C data on BV-173 cells available, we defined enhancers as H3K27ac+ non-promoter regions within 50 kb proximity of each gene of interest, focusing on STAT5-HA-bound regions at AK-2292/Ponatinib-sensitive genes. AK-2292-induced loss of STAT5 in parental cells and dTAGV-1-induced STAT5 loss in degron cells indeed caused substantial reductions of H3K27ac at enhancers and promoters of the STAT5 target genes CISH and SOCS1 (Fig.5I). Likewise, H3K27ac signals were reduced at STAT5-HA-bound enhancer regions when quantifying all AK-2292/Ponatinib-downregulated genes (Fig.5J). The opposite effect was observed for STAT5-HA-bound enhancers of AK-2292/Ponatinib-upregulated genes (Suppl.Fig.4D/E), and H3K27ac signals at promoters largely followed this trend (Suppl.Fig.4F/G).
Together, these results show that STAT5 function in Ph+B-ALL cells is largely dedicated to the regulation of BCR::ABL1 target genes, including the activation and repression of associated enhancers. However, its contribution to the BCR::ABL1-induced transcriptional program appeared surprisingly small, indicating that BCR::ABL1 employs additional TFs besides STAT5 for transcriptional reprogramming.
Therapeutic targeting of the enhancer-activators CBP/P300 causes cell cycle arrest and consecutive cell death of Ph+B-ALL cells
Given that key genes in Ph+B-ALL cells are enhancer regulated, we next wondered whether Ph+B-ALL cells are especially sensitive to CBP/P300 inhibition, which is essential for enhancer activation50 as described before. Of note, CBP/P300 dependency had previously been described for leukemia and lymphoma cells67. However, B-ALL itself or individual B-ALL subtypes such as Ph+B-ALL had not been assessed in these studies. We therefore first assessed public DepMap68 data, which defines cancer type-specific gene dependencies through negative Chronos/Ceres scores. This analysis showed that B-ALL indeed displays higher dependency for CBP and P300 compared to myeloid leukemia, T-cell leukemia and several B-cell lymphomas (Fig.6A, left), and that CBP/P300 dependencies further vary between individual B-ALL subtypes (Fig.6A, right). However, likely due to the generally low transduction efficiency of Ph+B-ALL cells with large CRISPR/Cas9 vectors that underlies Dependency Mapping68, DepMap did not provide information on Ph+B-ALL cells. Therefore, we investigated CBP/P300 dependency of Ph+B-ALL cells ourselves experimentally using the recently developed CBP/P300-specific degrader/PROTAC, dCBP-167. We chose dCBP-1 over A-485 here as we were not able to achieve complete H3K27ac loss with A-485 in Ph+B-ALL cells (Fig.2F), while dCBP-1 almost completely abolished CBP/P300-dependent H3K27 acetylation (Fig.6B). Notably, dCBP-1 caused a marked G1 cell cycle arrest after 2 days of treatment (Fig.6C) and induced apoptosis in almost all Ph+B-ALL cells by day 5 of treatment (Fig.6D/E). Importantly, Ph+B-ALL cells showed similar sensitivity to dCBP-1 as multiple myeloma (MM) cells (Fig.6E), which so far were described the most CBP/P300-inhibition sensitive cancer cells51,67 and for which CBP/P300-inhibitors are currently in clinical trials.