Loss of p27 accelerates CCA formation and predisposes tumor cells to DNA damage
To study the impact of p27 in cholangiocarcinoma we used a previously established mouse model in which Notch and AKT expressing vectors (Notch-ICD; myr-AKT) are stably delivered into the liver by hydrodynamic tail-vein injection (HDTVi) [10]. We had previously shown that aberrant expression of NOTCH receptors is often found in human CCA and that mice which overexpress NOTCH-ICD develop CCA [3]. Importantly NOTCH-ICD will induce transcription of the Cyclin E gene which in turn results in the formation of genetically unstable cells and cellular transformation [3].
To ablate p27 expression in combination with Notch/Akt expression we generated an additional vector system co-expressing a p27 specific sgRNA together with the CAS9 endonuclease (Supplemental Fig. S1A) (further denoted as Notch/AKT/sgp27 mice). Transfection efficacy, evaluated seven days post HDTVi, was unaltered independent of the p27 genotype (Supplemental Fig. S1 B and C). After injection of these plasmids into C57BL/6J mice, we performed in vivo imaging by MRI and observed formation of CCA with a latency of 4 to 5 weeks. As expected, loss of p27 resulted in a higher incidence of hepatic lesions 4 weeks after HDTVi with 68% of the Notch/AKT/sgp27 showing MRI-detectable tumors at 4 weeks, while only 29% of Notch/AKT mice showed tumors (Fig. 1A). No difference in liver weight, body weight, liver/body weight ratio or serum liver enzymes was found (Supplemental Fig. S2). Chemotherapy involving gemcitabine constitutes the current standard of care for advanced CCA. To test the role of p27 in the response to treatment, we initiated treatment by intraperitoneal injection of gemcitabine (0.1 mg/ gram mouse body weight) twice per week. To ensure that the tumor incidence in both groups was comparable, we only started the treatment when the mice were tumor-positive by MRI. No difference in survival of vehicle treated mice caring wtp27 or p27 deficient tumors were indicated (Fig. 1B). Surprisingly, mice with p27 deficient tumors survived significantly longer compared to p27 wildtype tumor mice (52 vs. 24 days CI: 0.15–1.43, p = 0.002) (Fig. 1C) after gemcitabine treatment. Histological analyses (H&E) and immune histochemistry for Ki67 revealed no significant differences in proliferation rate between tumors in both groups after Gemcitabine treatment (Fig. 1D and E). Also, no differences in liver weight, body weight, liver/body weight ratio or serum liver enzymes were found (Supplemental Fig. S3).
However, expression of the apoptosis marker cleaved caspase 3 and the DNA damage marker γH2AX were significantly higher in the p27 deficient compared to the wildtype mouse tumors. After treatment of both groups with gemcitabine this difference became even more apparent. To ensure that these differences were due to changes in the tumor cells we corroborated the results using immunohistochemical staining (Fig. 1F and G). Together these results pointed towards a defect in DNA repair in tumor cells which had lost p27.
Loss of p27 promotes an aberrant homologous recombination-dependent DNA damage response during G1
To understand the molecular mechanism by which loss of p27 results in increased levels of DNA damage we turned to a mouse embryonic fibroblast cell line with targeted deletions of the p27 gene. The use of MEF cells also allowed us to assess the function of p27 as a cyclin kinase inhibitor by using a cell line derived from a knock-in mouse strain which expresses a form of p27 which cannot bind to cyclin-cdk complexes
(p27c − k−)[11].
These three cell lines were synchronized and irradiated in different stages of the cell cycle (G1, S-phase). We found a significant increase in the number of γH2AX foci 24h after irradiation in p27 knockout cell which were irradiated during the G1-phase as compared to WT cells. No difference in γH2AX foci formation was found in cell lines which were irradiated during S-phase (Fig. 2A). Importantly, cells expressing the non-cyclin/cdk binding form of p27 (p27c − k−) behaved like wildtype cells after irradiation pointing towards a cdk independent mechanism by which p27 is involved in the DNA damage response. To ensure that this phenotype was not restricted to fibroblasts, cholangiocarcinoma cells (TFK1) in which p27 expression was ablated by siRNA were irradiated after synchronization. As in fibroblasts we detected a significant difference in the levels of γH2AX comparing p27siRNA to p27 wildtype cells when cells were irradiated in G1 but not after irradiation in S phase (Fig. 2B).
As p27 deficient cells accumulated γH2AX foci specifically when the damage occurred during G1 we first analyzed whether non-homologous end joining (NHEJ) dependent DNA repair was altered after loss of p27. The protein p53BP1 restricts end resection of DNA double strand breaks and promotes NHEJ pathway choice [12]. When cholangiocarcinoma cells (TFK1) cells were irradiated in G1 or S-phase no differences in the number of p53BP1 foci were detected in p27 deficient or control cells (Fig. 2C) indicating, that NHEJ was not affected.
Next, we tested if homologous recombination (HR) dependent DNA repair was affected by loss of p27. BRCA1 has been shown to promote HR and should therefore be largely restricted to cells which have undergone damage during S or G2-phases [13]. Surprisingly, we found that p27 depleted CCA cells in which DNA damage was induced during G1 displayed significantly more pBRCA1 positive foci than p27 expressing cells. However, no difference was detected in cells which were damaged in S phase (Fig. 2D). In agreement with an aberrant activation of HR, we also detected a significant increase in phosphorylation of the BRCA1 protein in p27 deficient TFK1 cells which were irradiated in G1 as compared to wildtype cells (Fig. 2E) [14].
To further corroborate the finding that HR dependent DNA double strand repair is activated during the G1 phase in p27 deficient cell we used pRAD51 immunofluorescence staining. RAD51 is often used as an indicator of HR as its loading on DBS has been shown to be essential for this process [15]. As shown in Fig. 2F p27 depleted TFK1cells showed significantly more pRAD51 foci compared to WT cell lines when damaged during G1. No difference in pRAD51 foci formation was observed, when the damage occurred during S phase (Fig. 2F and Supplemental Fig. S4A). This finding was also not dependent on the ability of p27 to act a as cyclin kinase inhibitor as p27C − K− MEF behaved like p27 wildtype cells while p27 deficient cells showed a significant increase in RAD51 foci formation (Supplemental Fig. S4B)
To directly test whether the HR-dependent DNA repair mechanism is indeed active in p27−/− cells and independent of cyclin kinase inhibitory activity, we created p27+/+, p27−/−, and p27c − k− cell lines with a stably integrated plasmid expressing a pDR-GFP recombination reporter which allows the detection of cells undergoing HR in cells [16] (Supplemental Fig. S5A). pDR-GFP expressing cells were transfected with the pCβA-SceI plasmid to allow the transient expression of the SceI endonuclease and synchronized in G1 to measure HR during this stage in the cell cycle (Supplemental Fig.S5B). p27−/− cells showed significantly increased GFP expression as compared to p27+/+ and p27c − k− cells indicating that these cells were able to activate a HR response after damage (Fig. 3A). HR is initiated with the resection of the 5´ strands to generate 3´ single-stranded DNA (ssDNA), which is required for RAD51 binding and strand invasion. We therefore compared 5´strand resection in p27−/− cells to p27+/+ cells using a previously described method [17] which uses the ER-AsiSI restriction enzyme to induce DSB at sequence specific sites. A higher percentage of ssDNA was found in p27−/− than WT HeLa cells synchronized in G1 upon expression of ER-AsiSI (Fig. 3B)
These data suggests that cells which have lost p27 are undergoing an aberrant attempt to repair DNA double strand breaks which occur during the G1 phase by homologous recombination. Given that during G1 no sister chromatids are available for recombination, we wondered if the increased propensity to initiate HR dependent repair might also increase the rate of interhomolog recombination. To test this, we used the previously described S/P cell line which allows for the detection of interhomolog recombination events [18]. To test the effect of p27 on homologous recombination we deleted p27 using CRISPR/Cas and derived two independent p27−/− S/P cell lines (Supplemental Fig.S6). Upon transient expression of the I-SceI endonuclease and selection of the surviving cells we found that loss of p27 led to a significant increase in interhomolog recombination (Fig. 3C).
Taken together, these findings indicate that p27 expression is necessary to block the initiation of HR during the G1 stage of the cell cycle. Aberrant activation of HR during G1 results in persistent DNA damage and cell death. However, loss of p27 also promotes recombination between homologous chromosomes. Surprisingly, the HR suppressive activity is independent of the ability of p27 to bind and inhibit cyclin-CDK complexes. We therefore tried to identify the underlying molecular mechanism which requires p27 to suppress HR before S-phase.
The interaction of p27 and RAD17 is necessary to prevent HR-dependent DNA double strand repair during G1
In an independent shRNA screen, which we conducted to identify resistance mechanisms to p27 stabilizing compounds [19], we found an interaction between p27 and the RAD17 protein [20, 21]. RAD17 is an essential component of the DDR machinery required to allow binding of the 9-1-1 (RAD9-HUS1-RAD1) clamp loading complex to damaged DNA [22–24]. Recently, RAD17 was shown to be also involved in the recruitment of the MRN (MRE11-RAD50-NBS1) complex to sites of DNA damage [25, 26]. Upon irradiation, RAD17 is phosphorylated by the ATM kinase at two conserved serines (S645 and S635) [27, 28] which is believed to be essential for the interaction of RAD17 with the MRN complex.
To analyze p27, RAD17 and ATM binding we synchronized CCA cells and measured complex formation of the endogenous proteins at different stages of the cell cycle
in irradiated and control cells. Immunoprecipitation experiments showed that p27 and RAD17 form a complex during the G1 phase which is not disrupted when cells are irradiated. However, RAD17 was no longer in complex with p27 when irradiated cells reached the G1/S transition and p27 levels dropped (Fig. 4A). Mutation of RAD17 phosphorylation sites (S635, S645) prevented dissociation of p27 and RAD17 after irradiation (Supplemental Fig. S7) indicating that the cell cycle regulated decrease in p27 levels at the G1/S transition and the phosphorylation of RAD17 regulate the release of RAD17 from p27. During S phase p27 is not expressed and hence no complex formation was detectable. As expected, ATM became phosphorylated at S1981 in response to irradiation irrespective of the cell cycle stage (Fig. 4A). A complex of RAD17 and ATM or phosphorylated ATM in irradiated p27 expressing cells was not detectable before the cells had reached the G1/S transition (Fig. 4A). In line with this observation, RAD17 phosphorylation (S635/S645) was only weakly detectable in irradiated G0 or G1 p27 WT cells while RAD17 was strongly phosphorylated when it was bound to ATM at the G1/S transition (Supplemental Fig. S8). Importantly knockdown of p27 expression promoted binding of RAD17 to ATM and allowed phosphorylation of RAD17 at S635 and S645 before cells had reached the G1/S transition (Supplemental Fig. S8). This data suggested that p27 expression during G1 is necessary to prevent RAD17 binding to and phosphorylation by ATM.
The initiation of double strand repair depends on the recruitment of the MRN (MRE11-RAD50, NBS1) complex onto damaged DNA. RAD17 was shown to possess an additional function in allowing MRN recruitment [25].
We expressed RAD17 wildtype or a non-phosphorylatable mutant (RAD17AA) in Hela cells in which the endogenous RAD17 was knocked down using siRNAs. At the same time, we depleted p27 by siRNA treatment to obtain cells which differ in their p27 status (wildtype or knockdown) and in their RAD17 (wildtype or phospho-mutant) status. Upon synchronization of these cells in G1 and irradiation with 4Gy we again detected complexes between ATM and RAD17 in cells in which p27 was depleted but not in wildtype expressing cells during the G1 phase (Supplemental Fig. S8).
To analyze complex formation between RAD17 and MRN components we immunoprecipitated RAD17 from damaged chromatin and measured the amount of bound MRE11, Rad50 and NBS1. In agreement with the unscheduled formation of ATM and RAD17 complexes we also detected significantly increased binding of MRN components to RAD17 on damaged chromatin (Fig. 4B). Importantly loss of RAD17 phosphorylation did not prevent ATM binding but prevented binding to MRN components almost completely. Moreover, loss of RAD17 phosphorylation prevented yH2AX appearance completely (Supplemental Fig. S9). This role of p27 in preventing premature MRN interaction with RAD17 on chromatin during G0 and G1 phase was independent of p27 ability to inhibit cyclin cdk complexes (Supplemental Fig. S10).
This data suggests that phosphorylation of RAD17 by the ATM kinase is critically important for the observed phenotype of unscheduled HR in p27 deficient cells. However, binding of ATM to RAD17 was not affected by RAD17 phosphorylation. We therefore tested if ATM kinase activity itself is important for complex formation. For this we synchronized irradiated CCA cells in G1 and measured the effect of caffeine (nonspecific kinase inhibitor), ATM inhibitor, CHK2, ATR and CHK1 inhibitor in p27 depleted or wildtype expressing cells. While loss of p27 resulted in increased binding of ATM and RAD17, caffeine and ATM inhibition significantly reduced complex formation. All other inhibitors had no effect (Fig. 4C).
Finally, we measured the levels of RAD17 phosphorylation in primary mouse cholangiocarcinoma tissues derived from our treatment study in p27 WT and deficient mice. In agreement with the results from cell lines loss of p27 resulted in a striking increase in RAD17 and ATM activation in tumors from p27 deficient mice which had been treated with gemicitabine as compared to wildtype p27 mice (Fig. 4D).
Using embryonic fibroblast and CCA cell lines as well as primary mouse tumor tissues this data shows that loss of p27 promotes premature binding of ATM to RAD17 during the G1 phase which results in phosphorylation of RAD17 which in turn leads to premature binding to the MRN complex to Rad17 onto damaged chromatin.
Loss of p27 sensitizes tumor cells against drugs that target the DDR response
As our data suggested a critical role for ATM inhibition, we tested a series of pharmacological inhibitors of the DDR for their ability to synergize with loss of p27 to affect the viability of cholangiocarcinoma derived mouse organoids. We tested ATM, ATR, DNA-PK, Chk1, Chk2 and PARP inhibitors in combination with either gemcitabine or irradiation on CCA organoids derived from p27 deficient or wildtype mice. We found that already 48hr after treatment with ATM or PARP inhibitors organoids derived from p27 deficient mice started to undergo apoptosis as measured by caspase 3 cleavage when combined with irradiation or gemcitabine treatment (Supplemental Fig. S11). Organoids derived from p27 wildtype mice were much more resistant to induction of cell death under these conditions (Fig. 5A, 5B). We then quantitated cell viability of the different organoid cell lines in response to different amounts of ATM or PARP inhibitor. While both inhibitors alone lead only to small changes in viability (Supplemental Fig. S12), combination with irradiation or gemcitabine resulted in significantly increased levels of cell death in the p27 deficient cholangiocarcinoma organoids (Fig. 4C,D) as compared to wildtype organoids. Moreover, organoids derived from p27 wild type mice were much more resistant to the combination of DNA damage and inhibition of PARP activity as determined via morphological examination (Fig. 4E). To ensure that the selective sensitivity to combination treatment was not caused by a defect in cdk binding we treated MEFs derived from p27+/+, p27−/− and p27c − k− mice with the PARP inhibitor Olaparib alone or in combination with irradiation. Only p27 knockout cells respond to combination treatment with apoptotic cell death while the other cell lines are largely resistant, reinforcing the notion that the DDR phenotype observed in p27 knockout cells is independent of cyclin kinase complex inhibition (Supplemental Fig. S13).
These results show that the aberrant induction of HR-dependent DNA damage repair in cells which have lost p27, sensitizes both cell lines and CCA organoids to drugs that target the DDR in combination with irradiation or gemcitabine. Moreover, our results also point towards a treatment option for patients suffering from cholangiocarcinoma. Therefore, we finally asked if a iCCA patient subgroup exists that corresponds to our murine iCCA model being characterized by a highly aggressive/poorly differentiated phenotype with low levels of p27 and high appearance of DNA damage markers such as γH2AX. To this end, we examined the presence of both markers by immunohistochemistry (IHC) in 178 iCCA patients using TissueMicroArrays (TMAs) containing well (G1), moderately (G2), and poorly differentiated (G3) tumors. We identified three distinct subgroups by unsupervised model-based clustering using p27 and γH2AX IHC scores as numerical variables and the histologic grade as categorical variable. Among these, cluster 3 corresponds very well to the murine CCA model exclusively consisting of poorly differentiated tumors (G3) showing low p27 and higher γH2AX IHC scores (Fig. 6A,B Supplemental Table S1) for detailed information regarding cluster characteristics). Notably, ~ 30% (53/178) of all analyzed patients fall into cluster 3 and could, based on our data, benefit from a PARP inhibitor-containing combinatorial treatment. Further in line with our findings is the fact that none of the analyzed iCCAs showed both, high p27 and high γH2AX levels (Fig. 6C red box).