Genome-wide CRISPR screen identifies PIK3R1 associated with resistance to pharmacological Notch inhibition in T-ALL
The efficacy of targeted therapies for treatment of cancer patients is often limited by development of drug resistance(11). Potential resistance mechanisms to pharmacological Notch1 inhibition mediated by GSI or CB-103 in T-ALL are currently unclear. Thus, we performed a genome-wide loss-of-function (LoF) CRISPR/Cas9 screen(12) to identify genes responsible for resistance to Notch inhibition and novel combination therapies for efficient treatment of human T-ALL. We used a Notch-dependent human T-ALL cell line, DND-41, which responds moderately to both GSI and CB-103 treatment in vitro(9). DND-41 cells stably expressing Cas9 were infected with human GeCKO v2 CRIPSR libraries, containing 123,411 sgRNAs targeting 19,050 genes and treated with either vehicle, GSI or CB-103 for 21 days enabling both positive and negative selection of sgRNAs (Fig. 1A).
sgRNAs targeting 293 (GSI-treated) and 131 (CB-103-treated) genes were identified as significantly depleted (P < 0.05, log2FC<-1) in GSI- and CB-103-treated T-ALL cells compared to vehicle control (Fig. 1B, C). Negatively selected sgRNAs indicate genes that, when inhibited, might function synergistically with Notch inhibition to effectively eradicate T-ALL cells. Pathway analysis revealed that significantly depleted genes were regulating MYC- and E2F signaling, as well as G2M checkpoint and mTOR signaling pathways (Fig. S1).
Conversely, sgRNAs targeting 178 (GSI-treated) and 76 (CB-103) genes were identified as significantly enriched (P < 0.05, log2FC > 1) in GSI- and CB-103-treated cells compared to vehicle control, indicating that the loss of these genes could confer resistance to Notch inhibition. Robust rank aggregation (RRA) method was used to identify genes preferentially lost in response to Notch inhibition (Fig. 1D, E). Among these genes, Phosphoinositide-3-Kinase regulatory subunit 1 (PIK3R1) was identified at the top of the list in the GSI versus DMSO screen and was also identified in the screen of CB-103 versus DMSO-treated T-ALL cells. The PIK3R1 gene encodes for the p85α regulatory subunit of PI3K, which contains an SH2 domain that binds to and inhibits the catalytic subunit (p110) of PI3Ks. Interestingly, PIK3R1 mutations were identified as drivers of tumorigenesis in ovarian cancer(13), endometrial cancer(14) and breast cancer(15). PIK3R1 hotspot mutations in the SH2 domain were also recently reported in pediatric T-ALL patients(5, 6). In addition, we noticed that the positive regulatory subunit of PI3K (PIK3CD, leukocyte-restricted catalytic p110δ subunit) was depleted in Notch inhibitor treated cells. Taken together, these observations suggested a key role of PI3K signaling in acquired resistance to Notch1 inhibition.
Loss of PIK3R1 renders T-ALL cells resistant to pharmacological Notch inhibition
To validate the screening results, we generated multiple PIK3R1 knockout clones in two different NOTCH1-driven T-ALL cell lines (DND-41, RPMI-8402) and stable knockdown clones in the NOTCH3-driven cell line TALL-1 (Fig. S2A-C). Loss of PIK3R1 in several T-ALL cell lines led to a mild growth advantage compared to non-targeting control sgRNA clones (NT) or scrambled shRNA controls (scr). In contrast, cell growth of all GSI- and CB-103-treated PIK3R1 knock-out (KO) or knock-down (KD) clones was significantly enhanced compared to NT or scr (Fig. 2A and Fig. S2D). We observed a significant decrease in the percentage of cells in S phase in NT or scr T-ALL clones when treated with GSI, confirming that GSI induces cell cycle arrest in T-ALL cells(7). However, this effect was alleviated in all PIK3R1 KO and KD cell lines under the same treatment conditions (Fig. 2B and Fig. S2E). Interestingly, we also observed cell cycle arrest in RPMI-8402 and TALL-1 control lines treated with CB-103 and the effect was significantly decreased when PIK3R1 was lost (Fig. 2C and Fig. S2F). Previously, we showed that CB-103 induces apoptosis in T-ALL cells(9). Correspondingly, CB-103 treatment for three days induced substantial apoptosis significantly in all three control cell lines. However, loss of PIK3R1 significantly ablated this effect (Fig. 2D and Fig. S2G). Altogether, these results suggest that loss of PIK3R1 confers resistance of T-ALL cells to Notch inhibition by protecting them from both drug-induced apoptosis and cell cycle arrest.
PIK3R1 deficiency leads to elevated gene expression of proliferation and pro-survival pathways in response to Notch inhibition
To gain insights how loss of PIK3R1 confers resistance to pharmacological Notch inhibition in T-ALL cells, we performed gene expression analysis (Fig. S3A). Treatment of RPMI-8402 cells for 24hrs with CB-103 resulted in significantly down regulation of genes associated to hallmark pathways including NOTCH signaling, MYC targets, and E2F targets (Fig. S3B) as previously reported(9), whereas GSI treatment resulted in significantly down-regulated MYC targets and MTOR signaling (Fig. S3B). We did not observe significant enrichment of hallmark pathways analyzing the gene expression differences of RPMI-8402 PIK3R1 KO versus NT cells, albeit a moderate trend of increased expression of PI3K-AKT and KRAS hallmark pathway genes (Fig. S3C). This might explain the mild growth advantage observed under normal culture conditions due to loss of PIK3R1 in T-ALL cells (Fig. 2A). Interestingly, Gene Set Enrichment Analysis (GSEA) from CB-103-treated KO versus NT cells revealed an enrichment in multiple hallmarks including E2F targets, MYC targets, PI3K-AKT-MTOR signaling, G2M checkpoint and Apoptosis pathways (Fig. 3A). Increased expression of MYC target genes was also observed in GSI-treated KO versus NT cells (Fig. S3D). Specifically, up-regulation of key E2F family transcriptional activators including E2F1, E2F2, E2F3, cell cycle regulators CCND2, CCND3, and down-regulation of the transcriptional repressor E2F5 were observed. In addition, we detected significant up-regulation of anti-apoptotic genes such as BCL2 and BCL-xL (Fig. 3B, C). In contrast, typical Notch target genes including MYC, HES1 or DTX1 were equally down-regulated in CB-103-treated PIK3R1 KO and NT cells (Fig. S3E). These results are consistent with the increased proliferation and survival observed in drug-treated PIK3R1 KO vs NT cells (Fig. 2B-D), and provide some mechanistic insight for Notch inhibitor resistance.
Notch-inhibited PIK3R1-mutant T-ALL cells reveal major phosphorylation changes in the cell cycle and spliceosome machinery
The p85 protein encoded by PIK3R1 is part of an important kinase signaling complex. Its loss may lead to immediate altered signaling events. Therefore, we performed total- and phospho-proteome analysis of NT and PIK3R1 KO cells treated with DMSO or CB-103 (Fig. S4A, B). Across samples, we quantified 29904 peptides corresponding to 7886 protein groups and 25221 phosphopeptides, of which 21601 were categorized as class I phosphosites(16) originating from 5531 phosphoproteins (Fig. 4C). At the total protein level, we observed 54 (NT, CB-103 versus vehicle), 215 (PIK3R1 KO versus NT) and 206 (PIK3R1 KO CB-103 versus NT CB-103) significant changes (Fig. 4D). The comparisons at the phosphorylation level revealed 2983 (NT, CB-103 versus vehicle), 2636 (PIK3R1 KO versus NT) and 3731 (PIK3R1 KO CB-103 versus NT CB-103) significant changes (Fig. S4E). Thus, changes occurring at the level of phosphorylation profiles are much more pronounced compared to changes of the total proteome. KEGG analysis of total protein changes of CB-103-treated PIK3R1 KO versus NT cells, identified cell cycle regulation as the most significantly affected pathway (Fig. S4F), which corroborated observations from the RNA-seq data. Similar analysis at the phosphoproteome level pointed to cell cycle and spliceosome as the most significant alterations (Fig. 4A).
To dissect kinase regulation in more detail, we performed Kinase-substrate Enrichment Analysis (KSEA) on the differential phosphorylation profiles of our comparison groups (Fig. 4B). The analysis of CB-103 versus vehicle revealed that CB-103 treatment led to decreased AKT1, MTOR, and S6K signaling, whereas the PIK3R1 versus NT comparison showed the expected reciprocal outcome, with increased PKC family, AKT signaling, due to loss of p85, (Fig. S5A). Importantly, comparison of PIKR1 KO CB-103-treatment versus NT CB-103-treatment showed increased activating phosphorylation events for AKT 1/2/3, PKC family, and S6K, which were maintained and no longer downregulated by CB-103 treatment (Fig. 4B).
Subsequently, we examined interactions among key proteins (Fig. 4A) using experimentally validated knowledge from the STRING database (Fig. 4C) and highlighted phosphorylation changes on these proteins (Fig. 4D). This detailed phopho-mapping provides insights regarding functionally established phosphorylation events such as S780 for RB as well as less examined events including T451 on AKT2, which has previously been associated with oncogenic signaling (Fig. 4D). Immunoblotting validated key phosphorylation events for AKTs, S6K, RB1 and BAD, which are important regulators of proliferation and cell survival (Fig. 4E). CB-103 treatment resulted in marked downregulation of NICD and c-MYC(9), along with reduced total AKT levels and more pronounced reduced phosphorylation at residues T308, S473 and T450. Yet, these effects were largely ablated in p85-deficient cells (Fig. 4E and Fig. S5B). Similarly, the phosphorylation of ribosome protein S6 kinase (p-S6K, T389, T421/S424) and its downstream substrate S6 (p-S6, S235/236) were down-regulated by CB-103 treatment but not in p85-deficient cells. Thus, loss of PIK3R1 indeed helps to maintain proteins involved in protein translation under CB-103 treatment. In addition, all phosphorylation sites of RB tested (S780, S795 and S807/811) were down regulated in CB-103 sensitive compared to the resistant cells (RB). The same holds true for BCL2 and BCL-xL, whereas BAD and p-BAD levels (pro-survival) remained comparable. These results confirm that p85-deficient T-ALL cells are able to cope with Notch inhibition through increased AKT signaling and maintain protein translation, cell proliferation and pro-survival pathways.
Interestingly, LoF PIK3R1 led to prominent phosphorylation changes in proteins involved in the spliceosome and RNA processing in cells treated with pharmacological Notch inhibitors (Fig. 4A and Fig. S6A, B). This analysis allowed to establish changes in phosphorylation profiles of splicing factors upon altered PI3K signaling and highlighted a wide spectrum of so far uncharacterized phosphorylation sites. A recent report linked oncogenic PI3K signaling with splicing alterations in breast cancer on the transcriptional level(17). Thus, we reanalyzed our RNA-seq data for differentially expressed transcripts, which were indeed associated with genes involved in cell cycle and regulation of apoptosis signaling pathways (Fig. S6C, D). Furthermore, we assessed the differential exon usage using DEXSeq and identified a spectrum of genes with alternative exon usage events in PIK3R1 deficient cells in response to Notch inhibition compared to NT cells, including transcripts of ELL Associated Factor 1 (EAF1) (Fig. S6E).
Our results show that loss of PIK3R1 in T-ALL cells led to increased PI3K-AKT signaling, causing major phosphorylation changes in the cell cycle and spliceosome machinery changes that resulted in downstream activation of cell cycle progression, increased cell proliferation, E2F gene activation, increased protein synthesis and cell survival. Changes in the spliceosome at phosphorylation levels correlated also with differential splicing at the transcriptional level. Consequently, these mechanisms contribute to resistance to Notch inhibition in T-ALL (Fig. 4F).
Pharmacological Notch inhibitors synergize with targeted therapies in human T-ALL cells
The advantage using a CRSIPR/Cas9 screen in T-ALL cells under drug selection is that it not only allows for identification of candidate genes mediating drug resistance such as PIK3R1, but also genes and pathways crucial for survival under drug selection. This opens avenues to identify novel combination therapies. Preferentially depleted sgRNAs in GSI- and CB-103-treated T-ALL cells pointed to well established signaling components within T-ALL, including components of the IL7/JAK pathway (IL7R, JAK1), regulators of the cell cycle machinery (CDK6:CCND3), and the key gene encoding the PI3K catalytic subunit (PIK3CD) (Fig. 1B, C).
We validated these candidates using available FDA-approved inhibitors against CDK4/6 (PD-0332991), JAK1/2 (Ruxolitinib), and PIK3δ (CAL-101). We first established in vitro sensitivity profiles, and observed that the single agent IC50 of CB-103 for DND-41 cells was 4.3mM and 0.1mM for PD-0332991. We then tested a combination treatment administrating CB-103 and PD-0332991 at three fixed ratios of their corresponding IC50 (1:1, 1:2.5 and 1:0.5) and established dose response curves (Fig. 5A). The combination treatment increased sensitivity of the cells to CB-103 by lowering its IC50 to approximately 0.1µM, which is 43-fold lower than single agent treatment (Fig. 5A). Similarly, combination of PD-0332991 and GSI lowered the IC50 of GSI approximately 100-fold (Fig. 5A). The Combination Index(18) (CI) was 0.06 for CB-103 plus PD-0332991 and 0.0183 for GSI plus PD-0332991, both of which are below 0.1 indicating very strong synergism (Fig. 5B). In addition, combination treatment induces the downregulation of c-MYC which is downstream of Notch and p-RB as key cell cycle regulator in two independent T-ALL cell lines (Fig. 5C). Similarly, we observed very strong synergism combining Notch inhibitors with a JAK1/2 inhibitor or a PI3Kδ inhibitor (Fig. 5B). These findings suggest that Notch inhibition in combination with FDA-approved compounds targeting CDK4/6, IL7R signaling, or PI3K/AKT pathway should be more efficacious compared to single agent treatment.
These promising in vitro results prompted us to assess their efficacy in xenotransplantation assays. RPMI-8402 T-ALL cells expressing a luciferase reporter were transplanted into NSG mice to monitor tumor growth and progression of disease over time. Animals with established tumors were treated with single agent compounds (vehicle, CB-103, GSI, PD-0332991) or with combination therapy (CB-103 or GSI plus PD-0332991) for two weeks (Fig. 6A). The kinetics of tumor progression showed a moderate and statistically significant reduction in tumor burden for both single agent treatments of CB-103 or GSI compared to vehicle (Fig. 6B, C). Single agent treatment of PD-0332991 revealed a robust reduction in tumor burden. However, the strongest reduction in tumor burden was observed, when mice were treated with combination of PD-0332991 and either CB-103 or GSI (Fig. 6C). To test whether combination treatment led to an increase in overall survival of experimental animals, treatment was ceased after 2 weeks and tumor relapse and survival rates were monitored. Despite the short treatment window, the dual agent treatment of GSI plus PD-033291 translated into significant prolonged overall survival compared to other treatment regiments (Fig. 6D).
The PI3K-AKT axis was identified as a main switch of downstream signaling events responsible for resistance to Notch inhibition in the CRSIPR/Cas screen, RNA-seq and proteomics data. Therefore, we also tested dual treatment of the AKT inhibitor (MK-2206) combined with CB-103 and observed significant prolongation of overall survival with combination compared to single agent therapy (Fig. 6E).
In light of increased BCL2 expression in our RNA-seq data and a recent report on complete clinical response of a relapse refractory T-ALL patient, treated with the BCL2 inhibitor Venetoclax and CB-10310, we proceeded to assess the efficacy of combining CB-103 and Venetoclax in our model. Indeed, this combination treatment significantly extended overall survival compared to single agent treatment in a comparable range as with CB-103 plus MK-2206 (Fig. 6F). Overall, the CRSIPR/Cas9 screen in T-ALL cells unveiled potentially novel avenues of combination therapies.