- Open Access
Elevated H3K27me3 levels sensitize osteosarcoma to cisplatin
Clinical Epigeneticsvolume 11, Article number: 8 (2019)
In osteosarcoma (OS), chemotherapy resistance has become one of the greatest issues leading to high mortality among patients. However, the mechanisms of drug resistance remain elusive, limiting therapeutic efficacy. Here, we set out to explore the relationship between dynamic histone changes and the efficacy of cisplatin against OS.
First, we found two histone demethylases associated with histone H3 lysine 27 trimethylation (H3K27me3) demethylation, KDM6A, and KDM6B that were upregulated after cisplatin treatment. Consistent with the clinical data, cisplatin-resistant OS specimens showed lower H3K27me3 levels than sensitive specimens. Then, we evaluated the effects of H3K27me3 alteration on OS chemosensitivity. In vitro inhibition of the histone methyltransferase EZH2 in OS cells decreased H3K27me3 levels and led to cisplatin resistance. Conversely, inhibition of the demethylases KDM6A and KDM6B increased H3K27me3 levels in OS and reversed cisplatin resistance in vitro and in vivo. Mechanistically, with the help of RNA sequencing (RNAseq), we found that PRKCA and MCL1 directly participated in the process by altering H3K27me3 on their gene loci, ultimately inactivating RAF/ERK/MAPK cascades and decreasing phosphorylation of BCL2.
Our study reveals a new epigenetic mechanism of OS resistance and indicates that elevated H3K27me3 levels can sensitize OS to cisplatin, suggesting a promising new strategy for the treatment of OS.
Osteosarcoma (OS) is the most common primary malignant bone tumor worldwide, with an incidence of approximately one to three cases per million people per year; the male to female ratio is approximately 1.22:1 , and children/adolescents and those over the age of 75 with Paget’s disease or other bony lesions are the two susceptible populations [2, 3]. The femur (40%), tibia (20%), humerus (10%), and pelvis (8%) together account for almost 80% of the bones involved in OS based on reported cases; the other involved bones include the fibula, jaw, vertebrae, radius, ulna, and ribs . OS patients with metastasis and recurrence have a poor prognosis that is partly attributed to chemoresistance [4, 5]. Although several mechanisms have been reported to mediate chemoresistance and some chemoresistance-reversing agents have even been developed, the mechanisms remain unclear, and almost no substantial survival benefits are conferred by administration of the therapeutic agents . Above all, further exploration of the mechanism of chemoresistance and development of new strategies to reverse chemoresistance in OS are urgently needed.
Aberrant epigenetic reprogramming plays a pivotal role in tumorigenesis . Among epigenetic regulation factors, histone methyltransferases (HMTs) are frequently dysregulated in a spectrum of human tumors, which indicates that HMTs are potential therapeutic targets . Histone H3 lysine 27 trimethylation (H3K27me3) is frequently associated with gene repression and contributes considerably to tissue differentiation, stem cell properties, and cancer occurrence/progression . Increases and decreases in H3K27me3 levels are catalyzed by the histone methyltransferase enzyme enhancer of zeste homolog 2 (EZH2)  and two demethylases belonging to the Jumonji C (JmjC) domain-containing family: ubiquitously transcribed tetratricopeptide repeat on chromosome X (UTX, also called lysine-specific demethylase 6A [KDM6A]) and Jumonji D3 (JMJD3, also called KDM6B) . Thus far, increasing evidence has showed that dynamic methylation of H3K27 frequently occurs during the development of tumors, especially during the modulation of cancer stem cells (CSCs). H3K27me3 has been reported to be a negative modulator of breast CSCs (BCSCs), and upregulation of H3K27me3 content can suppress BCSCs ; the same phenomenon has also observed in ovarian CSCs  and in non-small cell lung cancer . However, it remains unclear whether alterations in histone modifications such as H3K27me3 foster tumorigenesis and chemoresistance in OS.
Here, we illustrate that alterations in H3K27me3 are highly correlated with cisplatin sensitivity in OS and that increased H3K27me3 levels can enhance cisplatin efficacy against OS. This effect may be mediated by direct downregulation of PRKCA and MCL1 expression and subsequent inactivation of RAF/ERK/MAPK cascades as well as decreased phosphorylation of BCL2.
Expression of KDM6A and KDM6B in OS cells and patient specimens
To explore the possible role of epigenetic methyltransferases and demethylases in OS, we first assessed the mRNA expression of 20 methyltransferases in OS cells treated with or without cisplatin (0.5 μM) for 24 h (Fig. 1a and Additional file 1: Figure S1). As shown in Fig. 1a, KDM6A and KDM6B were significantly upregulated after cisplatin treatment, while EZH2 was not induced. Since KDM6A and KDM6B are both demethylases of H3K27me3 [11, 15,16,17], we further investigated whether H3K27me3 levels were correlated with cisplatin resistance. The results indicated that H3K27me3 levels were higher in cisplatin-sensitive patient specimens than in cisplatin-resistant ones (Fig. 1b). We also analyzed the mRNA expression of KDM6A and KDM6B in 20 matched osteosarcoma and peritumoral tissues; the results showed that compared with normal peritumoral tissue, OS tissue expressed higher levels of KDM6A and KDM6B (Fig. 1c). Taken together, the results suggest that the expression of the H3K27 demethylases KDM6A and KDM6B is elevated in OS and that higher H3K27me3 levels indicate better efficacy of cisplatin against OS.
Downregulation of KDM6A and KDM6B sensitizes OS to cisplatin in vitro
To further investigate the roles of KDM6A and KDM6B in OS, we established KDM6A- and KDM6B-depleted 143B and HOS cells and verified the success of knockdown by western blotting (Additional file 2: Figure S2). The results of Cell Counting Kit-8 (CCK-8) and colony formation assays showed that knockdown of KDM6A or KDM6B sensitized OS cells to cisplatin (Fig. 2a, b). The quantitative results revealed that the colony formation efficiency of KDM6A- and KDM6B-deficient cells decreased by 60–90% after cisplatin (0.5 μM) treatment. Additionally, the percentage of cells undergoing apoptosis was higher among KDM6A- and KDM6B-knockdown cells than among control cells, as revealed by Annexin V/7-aminoactinomycin D (7-AAD) staining (Fig. 2c). The quantitative results suggested there were four to sixfold increases in the percentage of apoptotic cells upon KDM6A or KDM6B knockdown. These observation were further confirmed by the western blot results, which showed higher cleaved Caspase 3 in KDM6A- and KDM6B-knockdown cells than in wild-type cells (Additional file 3: Figure S3). Collectively, our results indicated that KDM6A or KDM6B knockdown could sensitize OS to cisplatin by enhancing apoptosis of tumor cells.
Increasing H3K27me3 levels sensitizes OS to cisplatin in vitro
KDM6A and KDM6B are the only two demethylases of H3K27me3 found so far , and KDM6A and KDM6B expression levels are associated with cisplatin sensitivity, as illustrated above. We subsequently verified the correlations between H3K27me3 levels and cisplatin resistance in OS. To achieve this, we used a selective inhibitor of KDM6A and KDM6B, GSK-J4 , and an inhibitor of EZH2 enzymatic activity, EPZ-6438 . The results of CCK-8 and colony formation assays showed that increased H3K27me3 levels (GSK-J4 treatment) sensitized OS cells to cisplatin (Fig. 3a and Additional file 4: Figure S4). In addition, the percentage of cells undergoing apoptosis was significantly higher among high-level H3K27me3 cells treated with cisplatin than among low-level H3K27me3 cells treated with cisplatin (Fig. 3b). As it is known that the histone γH2AX is a marker for persistent DNA damage , we also monitored the efficacy of cisplatin against OS cells alone or in the presence of EPZ-6438 (10 μM) or GSK-J4 (15 μM) for 24 h by analyzing γH2AX foci. We found that compared to that in control cells, the number of γH2AX foci increased greatly after cisplatin treatment in OS cells, especially when they were cotreated with GSK-J4. In contrast, the number of foci decreased when cells were cotreated with EPZ-6438 (Fig. 3c). These results were further verified by western blot analysis, which showed higher γH2AX expression in OS cells cotreated with cisplatin and GSK-J4 and lower γH2AX expression in cells treated with cisplatin and EPZ-6438 than in control cells (Fig. 3d). Altogether, our results indicate that increasing H3K27me3 levels sensitizes OS to cisplatin in vitro.
Upregulating H3K27me3 levels sensitizes OS to cisplatin in vivo
We next determined the relationship between H3K27me3 and cisplatin efficacy in vivo. Twenty-four athymic nude mice received subaxillary injections of 143B cells. On day 5, we randomly divided the mice into four groups: the control (Ctrl) group, the GSK group, the Cis group, and the Cis+GSK group, which were treated with physiological saline, GSK-J4 (50 mg/kg intraperitoneally (ip) once per day (qd) beginning on day 5), cisplatin (8 mg/kg ip twice per week (biw) beginning on day 7), and GSK-J4 (50 mg/kg ip qd beginning on day 5) and cisplatin (8 mg/kg ip biw beginning on day 7), respectively. Most strikingly, there was a dramatic reduction in the size and weight of OS tumors in the group of mice cotreated with cisplatin and GSK-J4 compared with the other groups (Fig. 4a). On average, the tumor volumes and weights in the control group on autopsy were fourfold and threefold greater, respectively, than those in the cotreated group (Fig. 4b, c). Histopathological analysis showed that there was more cell debris (marked by asterisks in Fig. 4d first panel), an indicator of necrosis, in the Cis and Cis+GSK groups, especially in the cotreated group, than in the control group, as visualized by hematoxylin and eosin staining. We also assessed a range of markers in tumors by immunohistochemistry. The results showed that Ki67 expression was lower in the Cis and Cis+GSK groups, especially in the cotreated group, than in the control group, while the number of tumor cells positive for cleaved Caspase 3 or terminal deoxynucleotidyl transferase dUTP nick-end labeling (TUNEL) staining (apoptotic markers) was mostly elevated in the Cis+GSK group (Fig. 4d); there were no significant differences in the expression of CD31 (a marker of endothelial-like cells ) among the different groups. Taken together, these data support and verify the in vitro observations and demonstrate that high levels of H3K27me3 alleviate tumor malignancy and enhance apoptosis in combination with cisplatin in in vivo xenograft models, suggesting that H3K27me3 alterations are important in determining chemosensitivity in OS.
H3K27me3 modification is associated with CSC properties in OS
Cancer is widely believed to develop from mutant premalignant stem and progenitor cells. Accumulated mutations can confer further malignant properties, such as uncontrolled growth and proliferative advantages, on these cells, leading to tumorigenesis. Moreover, the marker H3K27me3 is closely associated with the state of stem cells, and KDM6A and KDM6B are highly expressed in embryonic stem cells . Furthermore, there is a general assumption that CSCs, like their normal stem cell counterparts, are resistant to chemotherapy and/or radiotherapy and are therefore responsible for tumor relapse after therapy . Therefore, we conducted further experiments to determine whether H3K27me3 levels influence CSC properties (stemness) in OS. We performed sphere formation experiments, since it is acknowledged that more spherical clones and larger spherical sizes indicate stronger self-renewal ability, which is a characteristic of CSCs. The results indicated that downregulating H3K27me3 with EPZ-6438 in HOS and 143B cells markedly increased the number and size of spheres, while upregulating H3K27me3 with GSK-J4 decreased the number and size of spheres (Fig. 5a). Furthermore, quantitative polymerase chain reaction (qPCR) and western blot results showed that the expression of both SOX9 and CD117 (which are associated with stem cell maintenance [24, 25]) increased when cells were exposed to EPZ-6438 but decreased when cells were exposed to GSK-J4 (Fig. 5b, c). Collectively, our data indicate that H3K27me3 modification is closely associated with CSC properties in OS and that higher H3K27me3 levels are correlated with reduced CSC properties and lower expression of CSC-related genes.
Alterations in H3K27me3 regulate cell apoptosis- and stemness-related genes involved in OS chemosensitivity to cisplatin
To gain insight into the mechanism by which higher H3K27me3 levels sensitize OS cells to cisplatin, we performed RNAseq on control (Ctrl), cisplatin-treated (Cis), cisplatin- and EPZ-6438-cotreated (Cis+EPZ), and cisplatin- and GSK-J4-cotreated (Cis+GSK) 143B cells. The results indicated that a number of gene expression profiles were changed by exposure to cisplatin with or without EPZ-6438 or GSK-J4, especially genes closely related to apoptosis and stem cell characteristics (Fig. 6a); these findings are similar to the results of a functional protein association network analysis  (Fig. 6b). To find the key molecules involved in this process, we mined the RNAseq data more deeply and discovered that 882 genes were differentially expressed in the Cis group compared with the Ctrl group (circle a), 159 genes were upregulated in the Cis+EPZ group compared with the Cis group (circle b), and 331 genes were downregulated in the Cis+GSK group compared with the Cis group (circle c). Interestingly, we found that the expression of PRKCA, one of 23 genes shared by circles a, b, and c, was changed significantly (Fig. 6c). Real-time (RT)-qPCR analysis further verified that higher H3K27me3 levels downregulated the expression of PRKCA, MCL1, SOX9, CD117, and KLF4 but upregulated the expression of BMF (Fig. 6d). Altogether, our data revealed that alterations in H3K27me3 regulate apoptosis- and stemness-related genes involved in OS chemosensitivity to cisplatin.
H3K27me3 alterations affect OS chemosensitivity through PRKCA via phosphorylation of BCL2 and modulation of RAF/ERK/MAPK cascades
It is widely acknowledged that H3K27me3 methylation is associated with gene suppression, while demethylation of H3K27me3 is associated with gene activation. Since PRKCA, MCL1, SOX9, CD117, and KLF4 were upregulated in response to reduced H3K27me3 levels, we sought to determine whether the expression profile changes at these gene loci were directly associated with H3K27me3 alterations. Therefore, we conducted chromatin immunoprecipitation (ChIP)-qPCR, and the results showed that there were significant increases in H3K27me3 signals at the PRKCA and MCL1 loci when cells were treated with GSK-J4; these results were reversed when cells were treated with EPZ-6438 (Fig. 7a). PRKCA, also known as protein kinase C alpha, has been reported to mediate a variety of cellular functions including proliferation, apoptosis, and differentiation [27,28,29,30,31]. In general, PRKCA exerts antiapoptotic effects and may serve as a survival factor in some types of cells . In addition, in our study, overexpression of PRKCA in OS cells (PRKCA-OE cells), which was verified by western blot analysis (Additional file 5: Figure S5A), reversed the chemosensitization effects of GSK-J4 with regard to cisplatin, as demonstrated by flow cytometry (Additional file 5: Figure S5B). These findings provide further evidence for the role of PRKCA in mediating this process. Notably, the antiapoptotic effects of PRKCA have been shown to be mediated via activation of the RAF/ERK/MAPK cascade  and/or phosphorylation of the antiapoptotic protein BCL2 at serine 70 , and MCL1 is a member of the antiapoptotic BCL2 family that suppresses caspase activation . In this study, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis indicated that 51 pathway terms were associated with the genes differentially expressed between the Cis group and the Ctrl group (circle a), 32 pathway terms were associated with those between the Cis+EPZ group and the Cis group (circle b), and 46 pathway terms were associated with those between the Cis group and the Cis+GSK group (circle c) (Additional file 6: Table S5). In addition, seven terms were shared by circles a, b, and c (Fig. 7b), including the MAPK signaling pathway and the Rap1 signaling pathway (Fig. 7c). Based on the results mentioned above, we postulate that PRKCA regulates OS chemosensitivity by modulating BCL2 phosphorylation and the RAF/ERK/MAPK pathway. As expected, western blot analysis revealed that reductions in H3K27me3 levels were linked to upregulation of PRKCA, phosphorylated BCL2, phosphorylated RAF, and phosphorylated ERK, while total BCL2, RAF, and ERK remained stable (Fig. 7d). Taken together, our results illustrate that reductions in H3K27me3 levels induce OS chemoresistance via direct upregulation of PRKCA and MCL1 expression and that PRKCA subsequently phosphorylates BCL2 and activates RAF/ERK/MAPK cascades (Fig. 7e, f).
OS is a rare bone malignancy that primarily affects children and adolescents. At present, patients are typically treated with surgery and intensive adjuvant chemotherapy. However, unfortunately, the 5-year survival outcome for recurrent or metastatic OS is poor. Since the survival time of OS has not changed over the past two decades , new strategies and agents urgently need to be introduced to provide greater benefit to patients.
Epigenetic alterations are known to contribute greatly to carcinogenesis in humans ; therefore, epigenetic modulation could be a novel and promising strategy for fighting cancers, especially OS. Based on this possibility, we screened a panel of DNA methyltransferases and demethylases that could be induced by cisplatin in OS cell lines. The results indicated that KDM6A and KDM6B, which are both demethylases of H3K27, were upregulated by cisplatin treatment. We then sought to explore the correlations between H3K27me3 levels and chemosensitivity in OS and the underlying mechanism. The results indicated that the H3K27 demethylases KDM6A and KDM6B were upregulated in OS and that higher H3K27me3 levels indicated better efficacy of cisplatin against OS. In addition, KDM6A/KDM6B knockdown or upregulation of H3K27me3 levels sensitized OS to cisplatin by enhancing apoptosis in tumor cells.
It is natural to associate the link between epigenetic alterations and cancer chemosensitivity with stem cell properties, since epigenetic modification events are pivotal for stem cell renewal and differentiation and since CSCs are recognized to be more resistant than normal stem cells to chemotherapeutic agents . Hence, in our study, we further tested CSC properties upon up- or downregulation of H3K27me3.
Our findings revealed that higher H3K27me3 levels were correlated with reduced CSC properties and reduced expression of CSC-related genes such as SOX9 and CD117. Although we did not observe a significant change in the expression of SOX2 (a master regulator of stemness phenotype), we found that SOX9, which was also a SOX family member, showed a significant downregulation upon GSK-J4 treatment. SOX9 (SRY (sex-determining region Y)-box 9) is a member of the SOX family of transcription factors that plays critical roles in embryonic development, lineage commitment, and stem cell maintenance . Likewise, CD117/c-kit is a tyrosine kinase receptor associated with cancer progression and normal stem cell maintenance , which also significantly downregulated upon GSK-J4 treatment.
In our study, RNAseq data indicated that PRKCA, MCL1, and other genes are also involved in this process. More importantly, through ChIP-qPCR analysis, we found that the expression of PRKCA and MCL1 was directly modulated by changes in H3K27me3 levels at the loci of these genes. The PKCs, which include PRKCA, are pivotal components of many signaling pathways that regulate diverse cellular functions including proliferation, angiogenesis, migration, invasion, metastasis, and resistance to apoptosis . PRKCA (the α-isozyme of PKC) is widely expressed in various tissues, and abnormal levels of PRKCA have been found in many transformed cell lines and in several human tumors . PRKCA has also been reported to be a dominant tumor progression factor that is associated with tumor proliferation and survival . Pharmacological inhibition of PRKCA can reduce the growth and survival of tumors, promote apoptosis, and sensitize tumor cells to chemotherapy . Mechanistically, the antiapoptotic action of PRKCA has been shown to be mediated by activation of RAF/ERK/MAPK cascades  and/or by phosphorylation of the antiapoptotic protein BCL2 . Moreover, MCL1, is a member of the antiapoptotic BCL2 family, which can suppress apoptosis . As the RAF/ERK/MAPK pathway is an important signaling pathway involved in cell survival, apoptosis inhibition, cell cycle progression, and proliferation, we reasonably associated PRKCA with the RAF/ERK/MAPK pathway and the phosphorylation of BCL2. Fortunately, in line with our hypothesis, RAF/ERK/MAPK pathway inactivation was found to participate in the process of OS sensitization to chemotherapy by increasing H3K27me3 levels on PRKCA and MCL1 gene loci.
It remains unclear whether histone modification dynamics are the driving force of tumorigenesis and chemoresistance. PRKCA may also perform proapoptotic functions ; therefore, further investigation into the mechanism that influences the PRKCA functional mode is needed in the future. In addition, increasing numbers of studies have also demonstrated that epigenetic therapies enhance the efficacy of adoptive immunotherapies  and immune checkpoint inhibitors . Thus, whether the immune system is involved in the process of chemosensitization also needs further exploration.
In summary, we have shown, for the first time, that increasing H3K27me3 levels can promote the efficacy of cisplatin against OS by upregulating BMF and downregulating PRKCA, MCL1, SOX9, CD117, and KLF4, which subsequently dephosphorylates BCL2 and deactivates RAF/ERK/MAPK cascades. These findings may provide a promising strategy for enhancing chemotherapeutic efficacy against OS and improving prognoses of patients.
Patients with histologically confirmed OS who were involved in this study provided their informed consent, and all studies were performed following the guidelines of the ethics committee of Shanghai Ninth People’s Hospital. Twenty patients (six females and 14 males; there was no significant difference in composition regarding age or sex) were enrolled in this study (the clinical information is summarized in Additional file 6: Tables S1 and S2) and were divided into two groups based on their cisplatin efficacy outcomes: a cisplatin-sensitive group (cisplatin sensitivity 50–100%) and a cisplatin-resistant group (cisplatin sensitivity 0–50%).
After deparaffinization, rehydration, antigen retrieval, and blockade of endogenous peroxidase, tissues were incubated with anti-TUNEL (Roche), anti-cleaved Caspase 3 (Servicebio), anti-Ki67 (Servicebio), and anti-CD31 (Abways) primary antibodies at 4 °C overnight. Next, secondary antibodies were applied, and diaminobenzidine (DAB) (Dako) solution was used as a chromogen. Finally, the sections were counterstained with hematoxylin (Sigma-Aldrich) to identify nuclei.
Cell lines and cell culture
The human OS cell lines 143B and HOS were purchased from the Shanghai Institute of Biochemistry and Cell Biology (Shanghai, China). Cells were grown in Dulbecco’s modified Eagle’s medium (DMEM) (HyClone) supplemented with 10% fetal bovine serum (FBS; Gibco) and antibiotics (100 U/ml penicillin and 100 μg/ml streptomycin) at 37 °C in a humidified atmosphere with 5% CO2.
Cell viability assay
A cell viability assay was performed as previously described . Briefly, cells were added into 96-well plates at a density of 5 × 103 cells/well. The cells were incubated with or without the indicated agents for 72 h, and then CCK-8 reagent was added into each well at a volume of 1:10 of that of the culture medium. The cells were then cultured for 2 h while protected from light. Data were collected by reading the optical density (OD) at 450 nm.
RNA extraction and RT-qPCR
RNA extraction and RT-qPCR were performed as previously reported . In brief, RNA was isolated from cell lines using TRIzol Reagent (Invitrogen), and RNA samples (1 μg) were subjected to RT-qPCR using an RT-PCR kit (Takara) according to the manufacturer’s protocols. Relative quantification (RQ) was calculated as RQ = 2−ΔΔCt. The sequences of the primers are listed in Additional file 6: Table S3.
Protein extraction and western blot analysis
Cells were dissolved in lysis buffer containing protease inhibitors. Then, the proteins were subjected to 10% sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) and transferred to 0.22 μm polyvinylidene fluoride (PVDF) membranes. The membranes were blocked with 5% nonfat milk at room temperature for 2 h and then incubated with primary antibodies, including antibodies against GAPDH (Cell Signaling Technology), H3K27me3 (Active Motif), H3 (Active Motif), PRKCA (Abways), BCL2 (Cell Signaling Technology), pBCL2 (Cell Signaling Technology), ERK (Cell Signaling Technology), pERK (Cell Signaling Technology), RAF (Abways), pRAF (Abways), γH2AX (Abways), SOX9 (Abways), CD117 (Abways), and cleaved Caspase3 (Abways), overnight at 4 °C. After three washes in Tris-buffered saline with Tween 20 (TBST), the membranes were incubated with corresponding secondary antibodies for 2 h at room temperature.
Colony formation assay
Cells were seeded into six-well plates at a density of 1000 cells/well. After cell attachment, different agents were added. On day 7, the colonies were fixed and stained with 0.1% crystal violet. The number of colonies was counted manually.
Flow cytometry analysis
A flow cytometry assay was performed as described previously . In brief, the percentage of cells undergoing apoptosis was evaluated by flow cytometry (Beckman Gallios). Cells were incubated with different treatments for 48 h, harvested and resuspended in 500 μl of staining buffer with Annexin V-allophycocyanin (APC) and 7-AAD. The data were analyzed using FlowJo software.
Tumorsphere formation assay
Single-cell suspensions were added to ultra-low-attachment multiwell plates (Costar) in serum-free chemically defined medium  with different treatments. Half of the volume was replaced with fresh medium every second day. Images were obtained on day 10.
shRNA transfection and lentivirus packaging
KDM6A and KDM6B shRNA and the full protein-coding PRKCA open reading frame were synthesized by Sangon Biotech (Shanghai, China) and cloned into the pLVX plasmid (Open Biosystems) (sequences listed in Additional file 6: Table S4). Empty vectors were used as a negative control. Then, lentivirus was packaged using psPAX2 and pMD2G vectors. To obtain stable cell lines, lentivirus supernatant was added to the cells, and the cells were cultured with 2 μg/ml puromycin for 2 weeks.
Xenograft experiments and histologic analysis
All animals were housed and maintained in specific pathogen-free (SPF) conditions. All animal experimental protocols were approved by the ethics committee of Shanghai Ninth People’s Hospital. Twenty-four athymic nude mice (6–8 weeks old) each were injected subcutaneously with 1 × 106 143B cells. On day 5 after injection, the mice were randomly divided into four groups (n = 6): the Ctrl group, the GSK group, the Cis group, and the Cis+GSK group. Beginning on day 5, the GSK group and the Cis+GSK group received ip injections of GSK-J4 (50 mg/kg) every day. Beginning on day 7, the Cis group and the Cis+GSK group received ip injections of cisplatin (8 mg/kg biw). The volume of each tumor was measured as (length × width2/2) weekly after implantation until the mice were sacrificed.
Cells were fixed and stained with antibodies against γH2AX (Abways) or H3K27me3 (Active Motif) after being treated with different agents for 24 h. After being washed in phosphate-buffered saline (PBS) three times, the cells were stained with the antibodies Alexa Fluor 555 (Invitrogen) and Alexa Fluor 488 (Invitrogen) for 1 h at room temperature. Images were obtained with a Cell Observer (ZEISS).
RNAseq and KEGG analysis
RNA samples were obtained from control (Ctrl), cisplatin-treated (Cis) (cisplatin 1 μM, 24 h), cisplatin- and EPZ-6438-cotreated (Cis+EPZ) (cisplatin 1 μM, 24 h and EPZ-6438 10 μM, 24 h), and cisplatin- and GSK-J4-cotreated (Cis+GSK) (cisplatin 1 μM, 24 h and GSK-J4 15 μM, 24 h) 143B cells. After the proper pretreatments, the samples were subjected to RNAseq on a HiSeq 2500 system by Oebiotech Co. Ltd. Gene expression levels were analyzed with the programs htseq-count  and Cufflinks . All differentially expressed gene lists were generated with DESeq software .
Control and EPZ-6438 or GSK-J4 treated 143B cells (2 × 107) were crosslinked, lysed, and sheared using UCD-300 (Bioruptor) to ~ 200–700 base pairs in length. Then, ChIP was performed using an EZ ChIP kit (Millipore; 17–371) according to the manufacturer’s instructions, followed by quantification of ChIP-enriched DNA by RT-qPCR. The primer sequences are listed in Additional file 6: Table S3.
Each sample was analyzed in triplicate, and the experiments were repeated three times. The mean, standard error of the mean (SEM), and P values based on two-tailed t tests were calculated with Excel (Microsoft). Differences were considered significant at P < 0.05 (*P < 0.05, **P < 0.01).
Histone H3 lysine 27 trimethylation
Mirabello L, Troisi RJ, Savage S. Osteosarcoma incidence and survival improvement. Cancer. 2009;115:1531–43.
Kansara M, Teng MW, Smyth MJ, Thomas DM. Translational biology of osteosarcoma. Nat Rev Cancer. 2014;14:722–35. https://doi.org/10.1038/nrc3838 Nature Publishing Group.
Moore DD, Luu HH. Osteosarcoma. In: Peabody TD, Attar S, editors. Orthop Oncol Prim Metastatic Tumors Skelet Syst [Internet]. Cham: Springer International Publishing; 2014. p. 65–92. https://doi.org/10.1007/978-3-319-07323-1_4.
Yang J, Zhang W. New molecular insights into osteosarcoma targeted therapy. Curr Opin Oncol. 2013;25:398–406 Available from: http://content.wkhealth.com/linkback/openurl?sid=WKPTLP:landingpage&an=00001622-201307000-00011. Cited 30 Apr 2018.
Dieudonné F-X, Marion A, Haÿ E, Marie PJ, Modrowski D. High Wnt signaling represses the proapoptotic proteoglycan syndecan-2 in osteosarcoma cells. Cancer Res. 2010;70:5399–408 American Association for Cancer Research. Available from: http://www.ncbi.nlm.nih.gov/pubmed/20530678. Cited 30 Apr 2018.
Wang J, Seebacher N, Shi H, Kan Q, Duan Z. Novel strategies to prevent the development of multidrug resistance (MDR) in cancer. Oncotarget. 2015; Available from: http://www.oncotarget.com/fulltext/19187 http://www.oncotarget.com/fulltext/19187.
Xiao ZG, Shen J, Zhang L, Li LF, Li MX, Hu W, et al. The roles of histone demethylase UTX and JMJD3 (KDM6B) in cancers: current progress and future perspectives. Curr Med Chem. 2016;23:3687–96.
van Haaften G, Dalgliesh GL, Davies H, Chen L, Bignell G, Greenman C, et al. Somatic mutations of the histone H3K27 demethylase gene UTX in human cancer. Nat Genet. 2009;41:521–3 Available from: http://www.ncbi.nlm.nih.gov/pubmed/19330029. Cited 29 Apr 2018.
Mikkelsen TS, Ku M, Jaffe DB, Issac B, Lieberman E, Giannoukos G, et al. Genome-wide maps of chromatin state in pluripotent and lineage-committed cells. Nature. 2007;448:553–60 Available from: http://www.ncbi.nlm.nih.gov/pubmed/17603471. Cited 29 Apr 2018.
Cao R, Wang L, Wang H, Xia L, Erdjument-Bromage H, Tempst P, et al. Role of histone H3 lysine 27 methylation in Polycomb-group silencing. Science (80- ). 2002;298:1039–43 Available from: http://www.ncbi.nlm.nih.gov/pubmed/12351676. Cited 29 Apr 2018.
Agger K, Cloos PAC, Christensen J, Pasini D, Rose S, Rappsilber J, et al. UTX and JMJD3 are histone H3K27 demethylases involved in HOX gene regulation and development. Nature. 2007;449:731–4.
Yan N, Xu L, Wu X, Zhang L, Fei X, Cao Y, et al. GSKJ4, an H3K27me3 demethylase inhibitor, effectively suppresses the breast cancer stem cells. Exp Cell Res. 2017;359:405–14. https://doi.org/10.1016/j.yexcr.2017.08.024 Elsevier Inc.
Sakaki H, Okada M, Kuramoto K, Takeda H, Watarai H, Suzuki S, et al. GSKJ4, a selective jumonji H3K27 demethylase inhibitor, effectively targets ovarian cancer stem cells. Anticancer Res. 2015;35:6607–14 International Institute of Anticancer Research. Available from: http://www.ncbi.nlm.nih.gov/pubmed/26637876. Cited 19 Nov 2018.
Watarai H, Okada M, Kuramoto K, Takeda H, Sakaki H, Suzuki S, et al. Impact of H3K27 demethylase inhibitor GSKJ4 on NSCLC cells alone and in combination with metformin. Anticancer Res. 2016;36:6083–92 International Institute of Anticancer Research. Available from: http://www.ncbi.nlm.nih.gov/pubmed/27793936. Cited 19 Nov 2018.
Lan F, Bayliss PE, Rinn JL, Whetstine JR, Wang JK, Chen S, et al. A histone H3 lysine 27 demethylase regulates animal posterior development. Nature. 2007;449:689–94 Available from: http://www.ncbi.nlm.nih.gov/pubmed/17851529. Cited 29 Apr 2018.
Lee MG, Villa R, Trojer P, Norman J, Yan K-P, Reinberg D, et al. Demethylation of H3K27 regulates polycomb recruitment and H2A ubiquitination. Science (80- ). 2007;318:447–50 Available from: http://www.ncbi.nlm.nih.gov/pubmed/17761849. Cited 29 Apr 2018.
De Santa F, Totaro MG, Prosperini E, Notarbartolo S, Testa G, Natoli G. The histone H3 Lysine-27 demethylase Jmjd3 links inflammation to inhibition of polycomb-mediated gene silencing. Cell. 2007;130:1083–94 Available from: http://www.ncbi.nlm.nih.gov/pubmed/17825402. Cited 29 Apr 2018.
Kruidenier L, Chung C, Cheng Z, Liddle J, Che K, Joberty G, et al. A selective jumonji H3K27 demethylase inhibitor modulates the proinflammatory macrophage response. Nature. 2012;488:404–8 Available from: http://www.ncbi.nlm.nih.gov/pubmed/22842901. Cted 29 Apr 2018.
Knutson SK, Warholic NM, Wigle TJ, Klaus CR, Allain CJ, Raimondi A, et al. Durable tumor regression in genetically altered malignant rhabdoid tumors by inhibition of methyltransferase EZH2. Proc Natl Acad Sci. 2013;110:7922–7 Available from: http://www.ncbi.nlm.nih.gov/pubmed/23620515. Cited 29 Apr 2018.
Reinhardt HC, Aslanian AS, Lees JA, Yaffe MB. p53-deficient cells rely on ATM- and ATR-mediated checkpoint signaling through the p38MAPK/MK2 pathway for survival after DNA damage. Cancer Cell. 2007;11:175–89 NIH Public Access. Available from: http://www.ncbi.nlm.nih.gov/pubmed/17292828. Cited 30 Apr 2018.
Sobczyńska-Rak A, Silmanowicz P, Polkowska I. Tumor angiogenesis – factors influencing the development of a tumor vascular network and assessment of neoangiogenesis in histopathological samples. Med Weter. 2016;72:542–8 Available from: http://www.medycynawet.edu.pl/index.php/archives/382/5563-summary-med-weter-72-9-542-548-2016.
Boyer LA, Plath K, Zeitlinger J, Brambrink T, Medeiros LA, Lee TI, et al. Polycomb complexes repress developmental regulators in murine embryonic stem cells. Nature. 2006;441:349–53 Available from: http://www.ncbi.nlm.nih.gov/pubmed/16625203. Cited 29 Apr 2018.
Beck B, Blanpain C. Unravelling cancer stem cell potential. Nat Rev Cancer. 2013:727–38 Available from: http://www.ncbi.nlm.nih.gov/pubmed/24060864. Cited 29 Apr 2018.
Bongiovanni L, Caposano F, Romanucci M, Grieco V, Malatesta D, Brachelente C, et al. Survivin and Sox9: potential Stem cell markers in canine Normal, hyperplastic, and neoplastic canine prostate. Vet Pathol. 2018. https://doi.org/10.1177/0300985818794161.
Foster B, Zaidi D, Young T, Mobley M, Kerr B. CD117/c-kit in cancer stem cell-mediated progression and therapeutic resistance. Biomedicines. 2018;6:31 Available from: http://www.mdpi.com/2227-9059/6/1/31.
Szklarczyk D, Morris JH, Cook H, Kuhn M, Wyder S, Simonovic M, et al. The STRING database in 2017: quality-controlled protein–protein association networks, made broadly accessible. Nucleic Acids Res. 2017;45:D362–8 Oxford University Press. Available from: https://academic.oup.com/nar/article-lookup/doi/10.1093/nar/gkw937. Cited 10 Aug 2018.
Nakashima S. JB minireview-protein kinase C isotypes and their specific function protein kinase ca (PKCa): regulation and biological function. J Biockem. 2002; Available from: https://academic.oup.com/jb/article-abstract/132/5/669/919783.
Murakami M, Horowitz A, Tang S, Anthony Ware J, Simons M. Protein kinase C (PKC) δ regulates PKCα activity in a syndecan-4-dependent manner. J Biol Chem. 2002;277:20367–71.
Dempsey EC, Newton AC, Mochly-Rosen D, Fields AP, Reyland ME, Insel PA, et al. Protein kinase C isozymes and the regulation of diverse cell responses. Am J Physiol Lung Cell Mol Physiol. 2000;279(3):L429–38.
Romanova LY, Alexandrov IA, Nordan RP, Blagosklonny MV, Mushinski JF. Cross-talk between protein kinase C-α (PKC-α) and -δ (PKC-δ): PKC- α elevates the PKC-δ protein level, altering its mRNA transcription and degradation. Biochemistry. 1998;37:5558–65.
Deucher A, Efimova T, Eckert RL. Calcium-dependent involucrin expression is inversely regulated by protein kinase C (PKC) α and PKCδ. J Biol Chem. 2002;277:17032–40.
Kolch W, Heidecker G, Kochs G, Hummel R, Vahidi H, Mischak H, et al. Protein kinase Cα activates RAF-1 by direct phosphorylation. Nature. 1993;364:249–52.
Ruvolo PP, Deng X, Carr BK, May WS. A functional role for mitochondrial protein kinase Cα in Bcl2 phosphorylation and suppression of apoptosis. J Biol Chem. 1998;273:25436–42.
Budihardjo I, Oliver H, Lutter M, Luo X, Wang X. Biochemical pathways of caspase activation during apoptosis. Annu Rev Cell Dev Biol. 1999;15:269–90 Annual Reviews 4139 El Camino Way, P.O. Box 10139, Palo Alto, CA 94303–0139, USA. Available from: http://www.annualreviews.org/doi/10.1146/annurev.cellbio.15.1.269. Cited 18 Aug 2018.
Gill J, Ahluwalia MK, Geller D, Gorlick R. New targets and approaches in osteosarcoma. Pharmacol Ther. 2013;137:89–99 Pergamon. Available from: https://www.sciencedirect.com/science/article/pii/S0163725812002057?via%3Dihub. Cited 12 Aug 2018.
Dawson MA. The cancer epigenome: concepts, challenges, and therapeutic opportunities. Science. 2017;355:1147–52 American Association for the Advancement of Science. Available from: http://www.ncbi.nlm.nih.gov/pubmed/28302822. Cited 12 Aug 2018.
Liu B, Ma W, Jha RK, Gurung K. Cancer stem cells in osteosarcoma: recent progress and perspective. Acta Oncol (Madr). 2011;50:1142–50.
Jo A, Denduluri S, Zhang B, Wang Z, Yin L, Yan Z, et al. The versatile functions of Sox9 in development, stem cells, and human diseases. Genes Dis. 2014;1:149–61. Elsevier Ltd. https://doi.org/10.1016/j.gendis.2014.09.004.
Rosse C, Linch M, Kermorgant S, Cameron AJM, Boeckeler K, Parker PJ. PKC and the control of localized signal dynamics. Nat Rev Mol Cell Biol. 2010;11:103–12.
Denham DW, Franz MG, Denham W, Zervos EE, Gower WRJ, Rosemurgy AS, et al. Directed antisense therapy confirms the role of protein kinase C-α in the tumorigenicity of pancreatic cancer. Surgery. 1998;124:218–24.
Carter C. Protein kinase C as a drug target implications for drug or diet prevention and treatment of cancer. Curr Drug Targets. 2000;1:163–83 Available from: http://www.eurekaselect.com/openurl/content.php?genre=article&issn=1389-4501&volume=1&issue=2&spage=163. Cited 9 Aug 2018.
Wang X-Y, Repasky E, Liu H-T. Antisense inhibition of protein kinase Cα reverses the transformed phenotype in human lung carcinoma cells. Exp Cell Res. 1999;250:253–63.
Schönwasser DC, Marais RM, Marshall CJ, Parker PJ. Activation of the mitogen-activated protein kinase/extracellular signal- regulated kinase pathway by conventional, novel and atypical protein kinase C isotypes. Mol Cell Biol. 1998;18:790–8.
Kagoya Y, Nakatsugawa M, Yamashita Y, Ochi T, Guo T, Anczurowski M, et al. BET bromodomain inhibition enhances T cell persistence and function in adoptive immunotherapy models. J Clin Invest. 2016;126:3479–94.
Zhu H, Bengsch F, Svoronos N, Rutkowski MR, Bitler BG, Allegrezza MJ, et al. BET bromodomain inhibition promotes anti-tumor immunity by suppressing PD-L1 expression. Cell Rep. 2016;16:2829–37.
He C, Wu T, Hao Y. Anlotinib induces hepatocellular carcinoma apoptosis and inhibits proliferation via Erk and Akt pathway. Biochem Biophys Res Commun. 2018;503:3093–9 Available from: https://linkinghub.elsevier.com/retrieve/pii/S0006291X18317820. Elsevier Ltd.
Jiang C, He C, Wu Z, Li F, Xiao J. Histone methyltransferase SETD2 regulates osteosarcoma cell growth and chemosensitivity by suppressing Wnt/β-catenin signaling. Biochem Biophys Res Commun. 2018;502:382–8 Academic Press. Available from: https://www.sciencedirect.com/science/article/pii/S0006291X18312531?via%3Dihub. Cited 29 July 2018.
Gibbs CP, Kukekov VG, Reith JD, Tchigrinova O, Suslov ON, Scott EW, et al. Stem-like cells in bone sarcomas: implications for tumorigenesis. Neoplasia. 2005;7:967–76 Available from: http://www.ncbi.nlm.nih.gov/pubmed/16331882. Cited 8 Apr 2018.
Anders S, Pyl PT, Huber W. HTSeq--a Python framework to work with high-throughput sequencing data. Bioinformatics. 2015;31:166–9 Oxford University Press. Available from: https://academic.oup.com/bioinformatics/article-lookup/doi/10.1093/bioinformatics/btu638. Cited 25 May 2018.
Roberts A, Trapnell C, Donaghey J, Rinn JL, Pachter L. Improving RNA-Seq expression estimates by correcting for fragment bias. Genome Biol. 2011;12:R22 BioMed Central. Available from: http://genomebiology.biomedcentral.com/articles/10.1186/gb-2011-12-3-r22. Cited 25 May 2018.
Anders S, Huber W. Differential expression of RNA-Seq data at the gene level – the DESeq package. [cited 25 May 2018]; Available from: http://bioconductor.org/packages/release/bioc/vignettes/DESeq/inst/doc/DESeq.pdf
We thank Oebiotech Co. Ltd. for their help.
This study was supported by grants from the National Key R&D Program of China (2016YFC1100600); the Shanghai Ninth People’s Hospital Lin+ Program (JYLJ012); and the National High Technology R&D Program (863 Program) (2015AA033603).
Availability of data and materials
All original data are available upon request.
Ethics approval and consent to participate
This study was approved by the ethics committee of Shanghai Ninth People’s Hospital and was performed in accordance with the relevant guidelines.
Consent for publication
No parts of this manuscript are being considered for publication elsewhere.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. mRNA expression of methyltransferases and demethylases in OS cells with or without cisplatin treatment. (TIF 1114 kb)
Figure S2. Knockdown efficacy of KDM6A and KDM6B as measured by western blot analysis. (TIF 217 kb)
Figure S3. Knockdown of KDM6A or KDM6B enhances cisplatin-induced apoptosis in OS cells. Cleaved Caspase 3 expression in KDM6A- or KDM6B-knockdown and control OS cells with cisplatin treatment. (TIF 147 kb)
Figure S4. Upregulation of H3K27me3 levels sensitizes OS to cisplatin. Colony formation ability of OS cells subjected to the indicated treatments. *P < 0.05. (TIF 1853 kb)
Figure S5. Overexpression of PRKCA reverses the chemosensitization effects of GSK-J4 with regards to cisplatin. (A) Overexpression efficacy of PRKCA as measured by western blot analysis. (B) Apoptosis in PRKCA-overexpressing 143B and HOS cells and their parental cells as determined by flow cytometry. *P < 0.05, **P < 0.01. (TIF 881 kb)
Supplementary Tables. (DOC 138 kb)