The expression of CSNK2A1 is associated with poor prognosis of osteosarcoma patients
To evaluate the clinicopathological significance of CSNK2A1 expression in human osteosarcomas, we performed immunohistochemical staining for CSNK2A1. Representative images of the immunohistochemical expression pattern of CSNK2A1 are presented in Fig. 1a. The positivity of the immunohistochemical expression of CSNK2A1 was determined with receiver operating characteristic curve analysis (Fig. 1b). The cut-off point was determined at the point with the highest area under the curve to predict the death of osteosarcoma patients (Fig. 1b). The cut-off point was twelve, and the cases with immunohistochemical staining scores equal or greater than twelve were considered positive for CSNK2A1 immunostaining (Fig. 1b). With this cut-off point, CSNK2A1-positivity was significantly associated with sex (P = 0.046), higher tumor stage (P = 0.006), higher T category (P = 0.028), higher M category (P = 0.047), and latent distant metastasis (P = 0.025) (Table 1).
In univariate survival analysis, age (OS; P = 0.039, RFS; P = 0.018), tumor size (OS; P = 0.015, RFS; P = 0.019), tumor stage (OS; P = 0.015, RFS; P = 0.005), T category (OS; overall P = 0.059, RFS; overall P = 0.020), M category (OS; P = 0.007, RFS; P = 0.018), and CSNK2A1 expression (OS; P = 0.002, RFS; P < 0.001) were significantly associated with OS or RFS (Table 2). CSNK2A1-positivity predicted a 10.081-fold (95% confidence interval [95% CI]; 2.307–44.054) greater risk of death and a 12.179-fold (95% CI; 2.777–53.407) greater risk of relapse or death of osteosarcoma patients (Table 2). The Kaplan-Meier survival curves for OS and RFS of CSNK2A1 expression are presented in Fig. 1c.
Table 2
Univariate analysis with Cox proportional hazards regression analysis for the survival of 37 osteosarcoma patients
Characteristics | No. | OS | | | RFS | |
| | HR (95% CI) | P | | HR (95% CI) | P |
Age, years, ≥ 30 (vs. <30) | 13/37 | 2.599 (1.050–6.433) | 0.039 | | 2.902 (1.201–7.014) | 0.018 |
Sex, male (vs. female) | 25/37 | 0.825 (0.294–2.312) | 0.714 | | 0.657 (0.235–1.836) | 0.423 |
Tumor size, ≥ 8 cm (vs. < 8cm) | 18/37 | 3.359 (1.263–8.933) | 0.015 | | 3.076 (1.207–7.838) | 0.019 |
Stage, III & IV (vs. I & II) | 11/37 | 3.161 (1.255–7.961) | 0.015 | | 3.647 (1.472–9.034) | 0.005 |
T category, 1 | 17/37 | 1 | 0.059 | | 1 | 0.020 |
2 | 16/37 | 3.335 (1.150–9.668) | 0.027 | | 3.455 (1.188–10.043) | 0.023 |
3 and 4 | 4/37 | 3.946 (0.938–16.606) | 0.061 | | 5.895 (1.559–22.287) | 0.009 |
N category, N1 (vs. N0) | 3/37 | 4.841 (0.957–24.486) | 0.057 | | 2.800 (0.599–13.095) | 0.191 |
M category, M1 (vs. M0) | 8/37 | 3.973 (1.464–10.784) | 0.007 | | 3.349 (1.229–9.125) | 0.018 |
CSNK2A1, positive (vs. negative) | 21/37 | 10.081 (2.307–44.054) | 0.002 | | 12.179 (2.777–53.407) | < 0.001 |
OS, overall survival; RFS, relapse-free survival; HR, hazard ratio; 95% CI, 95% confidence interval. |
Multivariate analysis for OS and RFS was performed with the inclusion of age, tumor size, stage, T category, N category, M category, and CSNK2A1 expression. In multivariate analysis, CSNK2A1 expression was an independent indicator of OS and RFS (Table 3). CSNK2A1-positivity predicted a 10.147-fold (95% CI; 2.320-44.385, P = 0.002) greater risk of death and a 12.179-fold (95% CI; 2.777–53.407, P < 0.001) greater risk of relapse or death of osteosarcoma patients (Table 2).
Table 3
Multivariate Cox proportional hazards regression analysis for the survival of 37 osteosarcoma patients.
Characteristics | OS | | | RFS | |
| HR (95% CI) | P | | HR (95% CI) | P |
N category, N1 (vs. N0) | 5.099 (0.945–27.506) | 0.058 | | | |
CSNK2A1, positive (vs. negative) | 10.147 (2.320-44.385) | 0.002 | | 12.179 (2.777–53.407) | < 0.001 |
OS, overall survival; RFS, relapse-free survival; HR, hazard ratio; 95% CI, 95% confidence interval. The variables included in multivariate analysis were age, tumor size, stage, T category, N category, M category, and CSNK2A1 expression. |
In the univariate analysis of 26 osteosarcoma patients who received adjuvant chemotherapy, CSNK2A1 expression was significantly associated with OS (P = 0.008) and RFS (P = 0.003) (Table 4). In the multivariate analysis performed with the inclusion of age, tumor size, stage, T category, N category, M category, and CSNK2A1 expression, N category was an independent indicator of OS (P = 0.035), and CSNK2A1 expression was an independent indicator of OS (P = 0.006) and RFS (P = 0.003) (Table 5). CSNK2A1 expression predicted a 9.619-fold (95% CI; 1.932–47.877) greater risk in OS analysis and a 10.374-fold (95% CI; 2.244–47.968) greater risk in RFS analysis in osteosarcoma patients who received adjuvant chemotherapy (Table 4).
Table 4
Univariate and multivariate Cox proportional hazards regression analysis of the survival of 26 osteosarcoma patients who received adjuvant chemotherapy
Characteristics | No. | OS | | | RFS | |
| | | HR (95% CI) | P | | HR (95% CI) | P |
Univariate analysis | | | | | | |
| CSNK2A1, positive (vs. negative) | 13/26 | 7.741 (1.693–35.408) | 0.008 | | 10.374 (2.224–47.968) | 0.003 |
Multivariate analysis | | | | | | |
| N category, N1 (vs. N0) | | 18.317 (1.223-274.238) | 0.035 | | | |
| CSNK2A1, positive (vs. negative) | | 9.619 (1.932–47.877) | 0.006 | | 10.374 (2.244–47.968) | 0.003 |
OS, overall survival; RFS, relapse-free survival; HR, hazard ratio; 95% CI, 95% confidence interval. The variables included in multivariate analysis were age, tumor size, stage, T category, N category, M category, and CSNK2A1 expression. |
The expression of CSNK2A1 is involved in the resistance to the anti-proliferative effect of doxorubicin
In human osteosarcomas, especially in the patients who received adjuvant chemotherapy, the expression of CSNK2A1 was significantly associated with shorter survival. Therefore, we evaluated the effect of CSNK2A1 in doxorubicin cytotoxicity on osteosarcoma cells. In U2OS and KHOS/NP osteosarcoma cells, overexpression of CSNK2A1 has not affected the proliferation of cells (Fig. 3a and 3b). However, under doxorubicin treatment, overexpression of WT-CSNK2A1 attenuated the anti-proliferative effect of doxorubicin (Fig. 3a and 3b). Under doxorubicin treatment, the cellular proliferation of CSNK2A1-overexpressing cells was significantly higher than the cells transfected with empty vector (Fig. 3a and 3b). In contrast, knock-down of CSNK2A1 sensitized U2OS and KHOS/NP osteosarcoma cells to doxorubicin in the CCK8 assay and colony-forming assay (Fig. 3c and 3d). The proliferation of osteosarcoma cells that had a knock-down of CSNK2A1 was significantly lower compared with cells transfected with empty vector under treatment with doxorubicin (Fig. 3c and 3d).
CSNK2A1 induces resistance to doxorubicin-mediated apoptosis of osteosarcoma cells
When the U2OS and KHOS/NP osteosarcoma cells were treated with doxorubicin, the expression of cleaved PARP1, cleaved caspase 3, and BAX increased, and the expression of BCL2 decreased (Fig. 4a). Under of doxorubicin treatment of U2OS and KHOS/NP osteosarcoma cells, the expression levels of cleaved PARP1, cleaved caspase 3, and BAX increased, and the expression of BCL2 decreased with knock-down of CSNK2A1 (Fig. 4a). Overexpression of CSNK2A1 decreased the expression levels of cleaved PARP1, cleaved caspase 3, and BAX and increased the expression of BCL2 under doxorubicin treatment (Fig. 4a). In flow-cytometry apoptotic analysis, apoptosis of U2OS and KHOS/NP cells significantly increased with knock-down of CSNK2A1 compared with controls under treatment with doxorubicin (Fig. 4b).
CSNK2A1-mediated resistance to doxorubicin is associated with SIRT6 phosphorylation-mediated activation of the DNA damage repair pathway
When U2OS and KHOS/NP osteosarcoma cells were treated with doxorubicin, the expression of pSIRT6, pATM, pChk2, and γH2AX were increased (Fig. 5a). In addition, under treatment with doxorubicin in U2OS and KHOS/NP cells, knock-down of CSNK2A1 decreased the expression of pSIRT6, pATM, and pChk2, and overexpression of CSNK2A1 increased expression of pSIRT6, pATM, and pChk2 (Fig. 5a). The expression of γH2AX increased with the knock-down of CSNK2A1 and decreased with overexpression of CSNK2A1 (Fig. 5a). Based on the role of CSNK2A1 on the phosphorylation of SIRT6 on Ser338 [10] and the role of SIRT6 on the induction of the DNA damage repair pathway [21], we evaluated the effect of a mutation of the Ser338 phosphorylation site of SIRT6 on CSNK2A1-mediated activation of the DNA damage repair pathway triggered by treatment with doxorubicin. Overexpression of WT-CSNK2A1 or WT-CSNK2A1/WT-SIRT6 increased expression of pSIRT6, pATM, and pChk2, and decreased expression of γH2AX with doxorubicin treatment in KHOS/NP cells (Fig. 5b). However, despite overexpression of CSNK2A1, transfection of SIRT6-S338A mutant decreased expression of pSIRT6, pATM, and pChk2 that were induced by doxorubicin and overexpression of CSNK2A1 (Fig. 5b). In CCK8 and colony-forming assays, cells overexpressing WT-CSNK2A1 or WT-CSNK2A1/WT-SIRT6 were resistant to doxorubicin, and the resistance to doxorubicin was attenuated with transfection of SIRT6-S338A mutant (Fig. 5c and 5d). The number of cells and colonies were significantly decreased in the cells transfected with SIRT6-S338A mutant compared with the cells overexpressing WT-CSNK2A1 or WT-CSNK2A1/WT-SIRT6 under treatment with doxorubicin (Fig. 5c and 5d).
In an orthotopic xenograft model, the growth of tumor in vivo was significantly inhibited with doxorubicin treatment (4mg/kg, once a week, intraperitoneal injection) (Fig. 6a and 6b). However, the anti-tumor effect of doxorubicin was attenuated via overexpression of WT-CSNK2A1 or WT-CSNK2A1/WT-SIRT6 (Fig. 6a and 6b). The growth of tumors in vivo was significantly higher with overexpression of WT-CSNK2A1 or WT-CSNK2A1/WT-SIRT6 compared with transfection of empty vector under treatment with doxorubicin (Fig. 6a and 6b). The resistance to doxorubicin induced by overexpression of WT-CSNK2A1 and WT-SIRT6 was attenuated by mutation of SIRT6 on Ser338. In vivo tumor growth of KHOS/NP cells was also significantly decreased with transfection of SIRT6-S338A mutant compared with overexpression of WT-CSNK2A1 or WT-CSNK2A1/WT-SIRT6 (Fig. 6a and b, and 6c). Furthermore, pulmonary metastasis was significantly decreased in the KHOS/NP cells transfected SIRT6-S338A mutant (number of metastatic nodule per mouse: mean ± standard deviation, 0 ± 0) compared with cells with induced overexpression of WT-CSNK2A1 (number of metastatic nodule per mouse: mean ± standard deviation, 3.0 ± 1.2) or cells induced to overexpress WT-CSNK2A1/WT-SIRT6 (number of metastatic nodule per mouse: mean ± standard deviation, 3.5 ± 2.4) (Fig. 6d). There were no metastasis on the liver or kidney in all experimental groups.
Inhibition of CSNK2A1 with emodin potentiates the cytotoxic effects of doxorubicin
Based on the role of CSNK2A1 in resistance to doxorubicin, we evaluated the effects of emodin, a CSNK2A1 inhibitor, on osteosarcoma cells. Treatment with emodin inhibited proliferation of U2OS and KHOS/NP cells in a dose- and time-dependent manner (Fig. 7a). As shown in Fig. 3, knock-down or overexpression of CSNK2A1 did not affect the proliferation of osteosarcoma cells without doxorubicin treatment. However, in contrast, emodin showed significant anti-proliferative effects with 24 hours of treatment at a 0.4 mM concentration (Fig. 7a). Moreover, similar to when cells underwent a knock-down of CSNK2A1, 0.5 mM emodin treatment potentiated the cytotoxic effects of doxorubicin (Fig. 7b and 7c). Co-treatment with emodin and doxorubicin significantly inhibited proliferation of U2OS and KHOS/NP cells compared with treatment with either doxorubicin or emodin alone (Fig. 7b and 7c).