- Primary research
- Open Access
Quantitative phosphoproteomic analysis reveals chemoresistance-related proteins and signaling pathways induced by rhIL-6 in human osteosarcoma cells
Cancer Cell International volume 21, Article number: 581 (2021)
IL-6 plays a pivotal role in resistance to chemotherapeutics, including lobaplatin. However, the underlying mechanisms are still unclear. This study was to investigate the changes in phosphoproteins and their related signaling pathways in the process of IL-6-induced chemoresistance to lobaplain in osteosarcoma cells.
We performed a quantitative phosphoproteomic analysis of the response of SaOS-2 osteosarcoma cells to recombinant human IL-6 (rhIL-6) intervention prior to lobaplatin treatment. The cells were divided into the control group (Con), the lobaplatin group (Lob), and the rhIL-6-and-lobaplatin group (IL-6). Three biological replicates of each group were included. The differentially expressed phosphoproteins were subjected to Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. Netphos 3.1 was used for the prediction of kinases, and STRING was used for the visualization of protein–protein interactions. The conserved motifs surrounding the phosphorylated residues were analyzed using the motif-x algorithm. Western blot analysis was performed to verify the differential expression of p-FLNC, its predicted kinase and the related signaling pathway. The results of the bioinformatic analysis were validated by immunohistochemical staining of clinical specimens.
In total, 3373 proteins and 12,183 peptides, including 3232 phosphorylated proteins and 11,358 phosphorylated peptides, were identified and quantified. Twenty-three significantly differentially expressed phosphoproteins were identified in the comparison between the IL-6 and Lob groups, and p-FLNC ranked second among these phosphoproteins. GO and KEGG analyses revealed the pivotal role of mitogen-activated protein kinase signaling in drug resistance induced by rhIL-6. Four motifs, namely, -SPxxK-, -RxxSP-, -SP-, and -SPK-, demonstrated higher expression in the IL-6 group than in the Lob group. The western blot analysis results verified the higher expression of p-FLNC, AKT1, and p-ERK and the lower expression of p-JNK in the IL-6 group than in the Con and Lob groups. The immunohistochemical staining results showed that p-FLNC, AKT1 and p-ERK1/2 were highly expressed in platinum-resistant clinical specimens but weakly expressed in platinum-sensitive specimens, and platinum-resistant osteosarcoma specimens demonstrated weak expression of p-JNK.
This phosphoproteomic study is the first to reveal the signature associated with rhIL-6 intervention before lobaplatin treatment in human osteosarcoma cells. p-FLNC, AKT1, and MAPK signaling contributes to resistance to lobaplatin in osteosarcoma SaOS-2 cells and may represent molecular targets to overcome osteosarcoma chemoresistance.
Osteosarcoma is the most common malignant bone tumor in adolescents and is considered insensitive to radiotherapy and immunotherapy. The clinical treatment for osteosarcoma includes aggressive and complete surgical resection and accompanying preoperative and postoperative chemotherapy [1, 2]. Platinum-based drugs are commonly used chemotherapeutic agents for osteosarcoma and have been developed from the first-generation chemotherapeutic agent cisplatin to the new agent lobaplatin, which exhibits better antitumor activity and low toxicity . However, many patients still develop resistance to lobaplatin after a period of treatment.
IL-6 is a crucial cytokine involved in the pathogenesis of many chronic inflammatory diseases, including cancer. This cytokine can activate STAT3 through the classical pathway and trans-signaling pathway, resulting in either anti-inflammatory or proinflammatory effects . Targeting IL-6 could be an effective strategy for overcoming stroma-induced chemotherapeutic resistance to tocilizumab in gastric cancer . Additionally, according to previous studies, IL-6 is involved in resistance to chemotherapy in osteosarcoma cells in vitro [6,7,8]. However, to date, specific changes in phosphoproteins at the post-translational level have not been discovered.
Proteins should first be modified after translation to perform their specific functions, and site-specific phosphorylation is the most pivotal post-translational modification that can regulate the function of proteins by altering their activities, stability or interactions. Therefore, dysregulation of site-specific phospho-signaling is a hallmark of many cancers . Phosphoproteomic analyses offer great advantages in the investigation of the differences in both phosphoproteins and their related signaling pathways. In this study, we explored the differential expression of phosphoproteins during rhIL-6 intervention before treatment with the chemotherapeutic agent lobaplatin in SaOS-2 osteosarcoma cells. The potential site-specific phosphoproteins and the related signaling pathways that may contribute to lobaplatin resistance in osteosarcoma cells were revealed.
Materials and methods
Osteosarcoma cell culture and treatment
The osteosarcoma cell line SaOS-2 was kindly provided by Procell Life Science & Technology Co., Ltd. and cultured in Dulbecco’s modified Eagle’s medium (Gibco, USA) supplemented with 10% fetal bovine serum (Gibco) and 1% penicillin–streptomycin (HyClone, USA). The cells were grown in a humidified atmosphere with 5% carbon dioxide at 37 °C. Cells grown to 80% confluence were used in the subsequent experiments. Trypsin–EDTA solution (0.25%, HyClone, USA) was used for cell detachment when needed. The rhIL-6 cytokine (PeproTech, USA) was diluted in PBS containing 5% trehalose. Lobaplatin (MCE, USA, Cat# 135558-11-1) was diluted in ddH2O. The cells were divided into the following three groups: the Con group (osteosarcoma cells treated with the solvent of both rhIL-6 and lobaplatin), the Lob group (osteosarcoma cells treated first with a solvent of rhIL-6 and then with lobaplatin at a concentration of 10 μg/mL for 24 h) and the IL-6 group (osteosarcoma cells pretreated with 60 ng/mL rhIL-6 eight hours before the lobaplatin treatment).
The anti-FLNC monoclonal antibody (Cat# ab180941), anti-AKT1 monoclonal antibody (Cat# ab108202), anti-P38 monoclonal antibody (Cat# ab170099), anti-p-P38 monoclonal antibody (Cat# ab178867), and anti-GAPDH monoclonal antibody (Cat# ab8245) were purchased from Abcam; the anti-JNK antibody (Cat# 9252T), anti-p-JNK antibody (Cat# 4668T), anti-ERK1/2 antibody (Cat# 4695T), and anti-p-ERK1/2 antibody (Cat# 4370T) were purchased from Cell Signaling Technology; the anti-p-FLNC antibody (Cat# bs-13183R) was purchased from Bioss Antibodies; and the HRP-conjugated secondary antibodies (Cat# 31430, 31460) were purchased from Thermo Fisher Scientific.
Osteosarcoma cells (9 × 106) were seeded in 15-cm dishes: in total, 18 culture dishes were used to establish three replicates per sample. After the previously described treatment, the cells were washed three times with ice-cold PBS and collected by scraping into several centrifuge tubes. The cells were centrifuged at 4 °C and 1000 rpm three times for 5 min each, and the PBS buffer was removed. Then, the cells were lysed in precooled RIPA lysis and extraction buffer (Thermo Scientific, USA, Cat# 89900) supplemented with 1% protease inhibitor cocktail (Cell Signaling, Cat# 5871) and 1% phosphatase inhibitor (Cell Signaling, Cat# 5872). The lysates were centrifuged at 14,000g for 10 min at 4 °C. The supernatants were transferred to new tubes and quantified using a BCA Protein Assay Kit (Piece, USA). The protein lysates were stored at − 20 °C. Twenty micrograms of protein from each sample was mixed with 6 × loading buffer and boiled for 5 min. The proteins were separated by 10% SDS-PAGE and stained with Coomassie Blue R-250.
In total, 600 μg of protein from each sample was added to 30 μL of SDT buffer containing 100 mM dithiothreitol (DTT), 4% SDS, and 150 mM Tris–HCl. Uric acid (UA) buffer was used to remove DTT and other low-molecular-weight components by repeated ultrafiltration. Then, 100 mM iodoacetamide was added, and the samples were incubated in darkness for 30 min. The filters were washed several times with UA buffer and TEAB buffer. The protein suspensions were digested with trypsin (Promega, USA) overnight at 37 °C, and the filtrates were collected. The peptide content was estimated by measuring the optical density of a 0.1% (g/L) solution under 280 nm UV light irradiation using an extinction coefficient of 1.1 calculated, based on the frequencies of tryptophan and tyrosine in vertebrate proteins.
Enrichment of phosphorylated peptides
A 100 μg peptide mixture of each protein sample was labeled using TMT reagents (Thermo Fisher Scientific, USA). The peptides in the Con group were labeled with TMT-126, TMT-127N, and TMT-127C; the peptides in the Lob group were labeled with TMT-128N, TMT-128C, and TMT-129N; and the peptides in the IL-6 group were labeled with TMT-129C, TMT-130N and TMT-130C. Each group included three biological replicates. The labeled peptides were combined and desalted using a C18 cartridge. The peptide mixtures were processed with a HiSelect TiO2 Phosphopeptide Enrichment Kit (Thermo Fisher Scientific, A32993). The TiO2 flowthrough and wash fractions were pooled, and the phosphopeptides were enriched with a HiSelect Fe-NTA Phosphopeptide Enrichment Kit (Thermo Fisher Scientific, A32992). The TiO2 eluent and Fe-NTA eluent were dried via vacuum centrifugation at 45 °C and then dissolved in 0.1% formic acid buffer.
Each eluent was injected twice for nanoliquid chromatography–tandem mass spectrometry (nanoLC-MS/MS) analysis. The peptide mixtures were loaded onto a C18 reversed-phase analytical column (Thermo Fisher Scientific, USA, P/N164943) in buffer A (0.1% formic acid) and separated with a linear gradient of buffer B (80% acetonitrile and 0.1% formic acid) at a regular flow rate of 300 nL/min. The linear gradient was as follows: 6% buffer B for 5 min, 6–28% buffer B for 215 min, 28–38% buffer B for 10 min, 38–99% buffer B for 5 min, and 99% buffer B for 10 min. A Q Exactive Plus mass spectrometer (Thermo Fisher Scientific, USA) was coupled to an Easy-nLC chromatography system (Thermo Fisher Scientific, USA) and operated in the positive ion mode. The MS data were acquired using a data-dependent method that could dynamically choose the most abundant precursor ions from full scans (ranging from 350 to 1800 m/z) for higher-energy collisional dissociation (HCD) fragmentation. The survey scans were performed at a resolution of 70,000 at 200 m/z with an automatic gain control (AGC) target of 3e6 and a maximum injection time (IT) of 50 ms. The MS2 scans were performed at a resolution of 35,000 for HCD spectra at 200 m/z with an AGC target of 2e5 and a maximum IT of 120 ms with an isolation width of 2 m/z. Only ions with a charge state between 2 and 6 and a minimum intensity of 2e3 were selected for fragmentation. The duration of the dynamic exclusion of the selected ions was 30 s, and the normalized collision energy was set to 30 eV.
The raw MS/MS data files were processed using the MASCOT engine (Matrix Science, London, UK; version 2.6, RRID: SCR_014322) embedded in Proteome Discoverer 2.1 and searched against the UniProt_HomoSapiens_20386_20180905 database. A maximum of 2 missed cleavages was permitted. A precursor mass tolerance of 10 ppm was specified, and a tolerance of 0.05 Da was set for the MS2 fragments. In addition to the TMT labels, carbamidomethyl (C) was set as a fixed modification. The variable modifications included oxidation, acetylation, and phosphorylation. A peptide and protein false discovery rate of 1% was enforced using a reverse database search strategy. Phosphopeptides with a fold change > 1.2, maximum number of missed cleavages ≤ 2, and p-value < 0.05 were considered differentially expressed.
Clustering, GO and KEGG analyses
A bioinformatic analysis was performed to compare the IL-6 group and the Lob group. The quantitative data of the target protein set were normalized. Matplotlib software was used to classify the samples and protein expression in two dimensions simultaneously (distance algorithm: Euclidean, connection method: average linkage), and a hierarchical clustering heatmap was generated for visualization. The GO annotation of the target protein sets was performed with Blast2GO . Specifically, the target protein sets were compared with the appropriate protein sequence database using the sequence comparison tool NCBI BLAST + (ncbi-blast − 2.3.0 +) on the Linux server, and the first 10 aligned sequences with san E-value ≤ 1e−3 were retained for the subsequent analysis. The GO entries associated with the target protein set and the highest bit score alignment sequence were extracted using the Blast2GO Command Line (download link: www.geneontology.org). During the annotation, the similarity between the target protein sequence and the aligned sequence, the reliability of the source of the GO entry, and the structure of the GO-directed acyclic graph were comprehensively considered, and the exact GO entry extracted in the mapping process was assigned to the target protein sequence. After the annotation was completed, to further improve the annotation efficiency, the EBI database was searched for conserved motifs matching the target protein using InterProScan , and the motif-related functional information was annotated to the target protein sequence. ANNEX was used to further supplement the annotation information, and then links between different GO categories were established to improve the accuracy of the annotations.
Orthologous genes and their products with similar functions in the same pathway were grouped and assigned the same KO (or K) label. During the KEGG pathway annotation, the target protein sequences were classified according to their KO labels based on the KEGG gene database, and the information of the pathway in which the sequence participated was obtained according to the KO classification using KEGG Orthology and Links Annotation (KOALA) software .
A motif analysis of the differentially expressed phosphorylated peptides between the IL-6 group and the Lob group was performed using the motif-x algorithm. The prealigned sequences were 13-amino-acid peptides with serine, threonine, or tyrosine arranged in the center as the key residue. The human background was selected, and the score threshold was set to 1.0e-6. Then, the logos of different motifs were visualized using Momo software (version 5.0.1) .
Kinase prediction and PPI network analysis
Kinase prediction was performed using the web resource NetPhos 3.1 (http://www.cbs.dtu.dk/services/NetPhos/), which predicts the serine, threonine, or tyrosine phosphorylation sites in eukaryotic proteins using ensembles of neural networks after inputting the substrate sequences in FASTA format. The residue to be predicted was set to either serine or threonine according to the MS results. For each residue, only the best prediction was chosen. Kinases with a score higher than 0.5 were selected.
The PPIs were retrieved from the STRING database (http://string-db.org/), which shows the experimental data and predicted data by computational algorithms . The protein IDs were input, and Homo sapiens was selected as the species. Protein IDs with a score higher than 0.9 were selected. The degrees of connectivity of each predicted kinase of the differentially expressed phosphoproteins between the IL-6 group and the Lob group were calculated to verify the importance of the kinases in the PPI network. The kinase with connections to the highest number of other predicted kinases was considered the most important kinase.
Western blot analysis
Cell cultures were harvested in RIPA buffer containing protease and phosphatase inhibitors (Cell Signaling, USA; Cat# 5871, 5872). After quantification with a BCA assay kit (Piece, USA), 60 μg of protein was separated on 6–10% SDS-PAGE gels and then transferred to PVDF membranes (Millipore, USA) for the phosphoproteins. The membranes were blocked with 5% nonfat milk powder or BSA and incubated with primary antibodies (anti-FLNC–1:10,000, anti-p-FLNC–1:500, anti-AKT1–1:5000, anti-ERK1/2–1:1000, anti-p-ERK1/2–1:2000, anti-JNK–1:1000, anti-p-JNK–1:1000, and anti-GAPDH–1:7500) at 4 °C overnight. The antibody-bound proteins were incubated with an HRP-conjugated secondary antibody at room temperature for 1 h and then detected using ECL reagents (Pierce, USA) on a chemiluminescence imaging system (BIORAD, USA) as previously described .
The studies involving patient specimens were approved by the Ethics Committee of Tangdu Hospital, and written informed consent was provided by the patients or their relatives. Clinical osteosarcoma specimens were collected from forty osteosarcoma patients who were treated with platinum-based chemotherapy and then underwent resection of visible tumors in the Orthopedic Department of Tangdu Hospital between January 2010 and May 2016. Patients with tumors exhibiting greater than 50% regression were considered chemotherapy-sensitive . In addition, five-year follow-up data were considered; twenty-four specimens were classified into the chemotherapy-resistant group, and the other sixteen specimens were classified into the chemotherapy-sensitive group. The clinicopathological data of the osteosarcoma patients were recorded.
Paraffin-embedded tissue specimens were sliced into 4-μm sections (LEICA, Germany) and heated at 60 °C for two hours. Then, the sliced specimens were dewaxed in xylene and rehydrated in graded ethanol solutions following standard procedures. Citrate buffer (10 mM) was used for antigen retrieval, and 3% H2O2 was used to block endogenous peroxidase activity. The slides were incubated with primary antibodies (anti-FLNC–1:250, anti-p-FLNC–1:400, anti-AKT1–1:500, anti-ERK1/2–1:250, anti-p-ERK1/2–1:200, anti-JNK–1:100, anti-p-JNK–1:50, anti-P38–1:100, and anti-p-P38–1:500) at 4 °C overnight. After incubation with the secondary antibody (GK500710, Gene Tech, China) at room temperature for 30 min, the slides were treated with DAB substrate. The staining images were acquired under an Olympus BX51 microscope.
Each experiment was performed at least three times. The statistical analysis of the differences between two groups was conducted with Student’s t-tests. One-way ANOVA with Tukey’s post hoc test was performed for comparisons among multiple groups. For the analysis of the immunohistochemical staining, the mean optical density (MOD) was calculated using a Mann–Whitney U test. The data were analyzed with GraphPad Prism (version 7, RRID: SCR_002798). P-values < 0.05 were considered indicative of significant differences.
Proteomic screening of phosphorylation in osteosarcoma cells
The total protein expression levels in each group of SaOS-2 osteosarcoma cells were determined by sodium dodecyl sulfate–polyacrylamide gel electrophoresis (SDS-PAGE) and Coomassie Brilliant Blue staining to eliminate the effects of differential protein expression on the comparison of phosphoprotein abundances; the protein expression levels derived from the whole cell lysates of the Con group, the Lob group, and the IL-6 group were similar (Fig. 1a). In addition, base peak chromatograms were generated to assess the correspondence of peptides between different samples, and the findings showed good quality both within and between groups (Fig. 1b). The phosphorylation levels of proteins in each group were then evaluated. The percentages of phosphorylation on the amino acids threonine, serine, and tyrosine are shown in Fig. 1c, indicating the highest modification level on serine (89.7%) and the lowest modification level on tyrosine (0.2%).
A schematic outline of the phosphoproteomic analysis of chemoresistant osteosarcoma cells treated with rhIL-6 and lobaplatin is depicted in a flow chart (Fig. 1d). In total, 3373 proteins and 12,183 peptides were identified, including 3232 phosphorylated proteins and 11,358 phosphorylated peptides. Twenty-three significantly differentially expressed phosphorylated peptides (5 upregulated by more than 1.2-fold and 18 downregulated by less than 0.83-fold, with 39 phosphorylation sites (p < 0.05)) between the IL-6 group and the Lob group were identified. A hierarchical clustering heatmap was generated to show the phosphoproteins and their specific phosphorylation sites and, thus, provide a better visualization of the overall phosphoproteomic changes. p-FLNC ranked as the second-most highly expressed phosphoprotein (Fig. 2, Table 1; Additional file 1).
Functional annotation of the differentially expressed phosphoproteins (DEPs)
To investigate the differences in the differentially expressed phosphoproteins between the Lob group and the IL-6 group, GO and KEGG enrichment analyses and related analyses were performed. Among the significantly upregulated phosphoproteins, the most enriched GO term in the biological process category was protein localization to the plasma membrane (GO: 0072659, enrichment factor = 8.29, p = 0.00116), and in the molecular function category, the most enriched GO term was kinase binding (GO: 0019900, enrichment factor = 5.58, p = 0.0489). Among the significantly downregulated phosphoproteins, the main enriched GO terms were protein localization to the plasma membrane (GO: 0072659, enrichment factor = 8.29, p = 0.00116), toll-like receptor (GO: 0002224, enrichment factor = 27.88, p = 0.00216) and Hippo signaling (GO: 0035329, enrichment factor = 8.76, p = 0.0212) in the biological process category; mitogen-activated protein kinase binding (GO: 0031435, enrichment factor = 17.04, p = 0.00583) and dGTPase activity (GO: 0008832, enrichment factor = 153.32, p = 0.00652) in the molecular function category; and inflammasome complex (GO: 0061702, enrichment factor = 76.66, p = 0.0130), CD95 death-inducing signaling complex (GO: 0,031,265, enrichment factor = 51.11, p = 0.0194), and gap junction (GO: 0005921, enrichment factor = 38.33, p = 0.0258) in the cellular compartment category (Fig. 3a, b; Table 2).
The KEGG pathway analysis revealed that the pathway annotations TNF signaling pathway, NOD-like receptor signaling pathway, apoptosis and MAPK signaling pathway were enriched (including two proteins in each pathway) between the Lob group and IL-6 group, with enrichment factors of 9.29, 7.86, 6.97, and 4.20, respectively (Fig. 3c; Table 2).
Phosphorylation motifs in lobaplatin-resistant osteosarcoma cells
The motifs surrounding phosphorylation sites are generally believed to be specific substrates for kinases. To characterize the motifs identified in the Lob group and IL-6 group, the differentially expressed phosphorylated peptides were analyzed using the motif-x algorithm. Twenty-three motifs (specifically, twenty-two phosphorylated (p)-Ser motifs and one p-Thr motif) were identified among the differentially expressed phosphosites between the Lob group and the IL-6 group (Fig. 4a). Among the twenty-three differentially expressed motifs, nine p-Ser motifs were high-scoring motifs and are listed in Fig. 4b. Among these 9 motifs, the motif -SPxxK- was enriched in the upregulated phosphoproteins, five motifs (-SDDE-, SDxE-, -SPxR-, -RxxSxE-, and -SPE-) were enriched in the downregulated phosphoproteins in the IL-6 group compared with the Lob group, and the motifs -RxxSP-, -SP-, and -SPK- were enriched in both the upregulated and downregulated phosphosites. The motif logos were generated using the sequence logo tool from MEME (http://alternate.meme-suite.org/tools/meme).
Kinase prediction and PPI analysis
The kinases of the differentially expressed phosphoproteins between the IL-6 group and the Lob group were predicted using NetPhos 3.1 while comprehensively considering the phosphorylation sites identified by a tandem mass tag (TMT)-based quantitative phosphoproteomic analysis and the prediction scores. In total, twenty-two kinases were predicted to regulate the differentially expressed phosphosites, and the relationships between the predicted kinases and the significantly differentially expressed phosphoproteins were visualized with Cytoscape (Fig. 5a). The protein–protein interaction (PPI) networks of these predicted kinases were visualized using STRING (http://string-db.org/cgi/input.pl). The most enriched highly connected kinases between the IL-6 group and the Lob group were as follows: rac-alpha serine/threonine protein kinase (UniProtKB-P31749), which interacted with eleven predicted kinases; glycogen synthase kinase-3 beta (UniProtKB-P49841), which interacted with ten predicted kinases; and mitogen-activated protein kinase 14 (UniProtKB-Q16539), which interacted with nine predicted kinases (Fig. 5b; Table 3). Site-specific kinase and substrate predictions and network analysis demonstrated extensive interactions among the predicted kinases.
Verification of the TMT-based phosphoproteomic analysis results by western blot analysis
The phosphoprotein p-filamin-C was differentially upregulated between the IL-6 group and the Lob group, and its predicted kinase rac-alpha serine/threonine protein kinase was analyzed. Western blot analysis showed that the osteosarcoma cells in the IL-6 group exhibited high levels of p-filamin-C, AKT1 and p-extracellular signal-regulated kinase 1/2 (p-ERK1/2), indicating that during rhIL-6 treatment, which results in lobaplatin resistance, elevated AKT1 promotes a high level of p-filamin-C and activates the MAPK signaling pathway (Fig. 6). This finding also verified the KEGG pathway analysis results, showing that filamin-C was enriched in the MAPK signaling pathway (Table 2).
Immunohistochemical (IHC) staining of clinical specimens
In total, forty osteosarcoma specimens (from sixteen patients with resistance to platinum-based chemotherapy and twenty-four chemotherapy-sensitive patients) were obtained to further confirm the quantitative phosphoproteomic analysis results. The patients’ clinical data are recorded in Additional files 2 and 3. Hematoxylin and eosin (HE) staining showed increased necrosis in the chemotherapy-sensitive specimens compared with the chemotherapy-resistant specimens. Immunohistochemical staining of p-filamin-C and its predicted kinase AKT1, p-ERK1/2, p-JNK and p-P38 was performed. Regarding p-filamin-C, twelve chemotherapy-resistant osteosarcoma specimens exhibited strong positive staining, three exhibited moderate positive staining, and one exhibited negative staining. Regarding AKT1 in the same specimens, ten exhibited strong positive staining, four exhibited moderate positive staining and two exhibited negative staining. Regarding p-ERK1/2, thirteen exhibited strong positive staining, and three exhibited moderate staining. Additionally, regarding p-JNK, stronger positive staining was observed in the chemotherapy-sensitive specimens, and no significant difference was detected in p-P38 between the sensitive and resistant specimens (Fig. 7; Additional file 4).
Previous studies revealed the pivotal role of the cytokine IL-6 in the process of chemotherapeutic resistance [17,18,19]. However, no systematic mechanistic investigation has been conducted from a phosphoproteomic perspective. Phosphorylation functions as a switch for protein activity and is the main mechanism of signal transmission. In addition, phosphorylation is a crucial step in the regulation of pathophysiological processes in various cancers, and numerous drugs and drug candidates target protein kinases . Phosphoproteomic analyses can enrich our knowledge of phosphorylated proteins and identify and quantify phosphorylation sites, thereby providing insight into the specific mechanisms of biological processes, such as proliferation, metastasis, and drug resistance . In this study, we explored the differentially expressed site-specific phosphoproteins, their potential kinases, and the related signaling pathways during rhIL-6 intervention before lobaplatin treatment.
This study is the first to use a TMT-based quantitative phosphoproteomics approach to examine the role of post-translational changes in IL-6-induced lobaplatin resistance in osteosarcoma cells. We identified 3,232 phosphorylated proteins and 11,358 phosphorylated peptides. Regarding the IL6-induced changes in phosphoproteins, we compared the IL-6 group with the Lob group and identified twenty-three significantly differentially expressed phosphoproteins that were relatively tightly related to the effects of IL-6 on lobaplatin resistance.
The GO enrichment analysis suggested that IL-6 pretreatment before lobaplatin treatment significantly regulated kinase binding, Toll-like receptors, Hippo signaling, and mitogen-activated protein kinase binding (all p-values < 0.05) in osteosarcoma cells. Toll-like receptors (TLRs) are membrane-bound pattern recognition receptors that can recruit proinflammatory cytokines and costimulatory molecules through NFκB and MAP kinase signaling . Moreover, TLRs mediate immune evasion, which can result in autophagy dysregulation and chemoresistance . Hippo signaling regulates diverse biological processes, and its dysfunction is related to many cancers . Considerable evidence has demonstrated that the deregulation of YAP/TAZ signaling may play a key role in intrinsic and acquired resistance to various chemotherapies . Regarding the enrichment in MAPK binding, emerging evidence suggests that this event can modulate drug resistance and sensitivity in various cancers .
The KEGG pathway analysis provided further insight. The pathway annotations NOD-like receptor (NLR) signaling pathway, TNF signaling pathway, apoptosis and MAPK signaling pathway were enriched. NLRs are considered regulators of tumorigenesis, cancer cell stemness and chemoresistance. The inappropriate activation of NLRs modulates the tissue microenvironment and potentiates the risk for cancers [27, 28]. TNF signaling leads to cell survival, proliferation and differentiation. However, excessive or aberrant activation of TNF-α signaling can lead to chronic inflammation and other pathological complications . The MAPK signaling pathway is a highly conserved signaling pathway involving JNK, ERK1/2, p38, and ERK5 that mediates metabolic reprogramming, cell survival, and cell differentiation. Constitutive activation of this pathway plays a role in cancer development and progression, including autophagy dysregulation and therapeutic resistance [30,31,32]. The MAPK signaling pathway also participates in crosstalk with the NLR and TNF signaling pathways identified above as enriched [33, 34]. Phosphoproteins, such as p-CASP7, p-FLNC, p-TRAF2, and p-TP53BP1, which are involved in the abovementioned enriched pathways, may be targets for lobaplatin resistance.
Among the identified phosphoproteins, p-filamin-C exhibited the second-highest level in the osteosarcoma cells treated with rhIL-6 prior to lobaplatin compared with those treated with lobaplatin alone. Filamin C is a dimeric actin-binding protein that can regulate remodeling of the actin cytoskeleton and act as a scaffold for signaling proteins, including tyrosine kinases, phosphatases, and GTPases. Therefore, filamins are believed to connect adhesive receptors to signal transduction and the cytoskeleton . High filamin-C expression has been reported to enhance the invasiveness of glioblastoma and predict poor outcomes . Moreover, filamin-C can interact with and activate MEK1/2 and ERK1/2 to promote the progression of hepatocellular carcinoma . In contrast, another study showed that the downregulation of filamin-C by acetylated Siah2 increased the invasiveness of gastric cancer cells . Collectively, the above results suggest that filamin-C acts as a double-edged sword with both detrimental and beneficial roles in patients with cancer and that its role may depend on the type of cancer. However, knowledge regarding the role of phosphorylated filamin-C in cancers is limited. In our study, the p-filamin-C level in rhIL-6- and lobaplatin-treated osteosarcoma cells, which are resistant to lobaplatin as described in a previous study, was higher than that in cells treated with lobaplatin alone. In addition, its predicted kinase PKB (also called AKT1) showed the highest connectivity with other predicted kinases of the differentially expressed phosphorylated proteins between the Lob group and the IL-6 group. Moreover, p-filamin-C was enriched in the MAPK signaling pathway (ko04010), as shown by KEGG pathway analysis. MAPKs can phosphorylate substrates with the motif -SP- . Differential expression of this motif between the Lob group and the IL-6 group was observed, indicating that MAPKs might be activated by IL-6-induced chemoresistance to lobaplatin. Therefore, p-filamin-C is expected to be a candidate target for lobaplatin resistance in osteosarcoma. Western blot analysis and immunohistochemical staining were performed to verify the results of the phosphoproteomic and subsequent bioinformatic analyses, and the levels of p-filamin-C, AKT1 and p-ERK1/2 in the rhIL-6-pretreated osteosarcoma cells were found to be higher than those in the cells treated with lobaplatin alone, while lower expression of p-JNK and no significant change in the level of p-P38 were observed. In the clinical specimens from osteosarcoma patients, most chemotherapy-resistant specimens exhibited strong positive staining of p-FLNC, its kinase AKT1, and p-ERK1/2, a key protein in the ERK signaling cascade that regulates the growth and differentiation of cells; in contrast, the chemotherapy-sensitive specimens showed negative or weak positive staining of the above proteins (Additional file 5). Similar to the western blot analysis results, the level of p-JNK was higher in the chemotherapy-sensitive specimens, and no significant difference was found with respect to p-P38. Altogether, the above results suggest that activated MAPK signaling in the process of rhIL-6-induced chemoresistance mainly occurs through p-ERK1/2 and p-JNK, but not p-P38. ERK1/2 inhibitors and the activation of JNK can be considered strategies to overcome chemoresistance to lobaplatin in osteosarcoma.
In summary, in this study, the global phosphorylation patterns of SaOS-2 osteosarcoma cells treated with rhIL-6 prior to lobaplatin were characterized. Based on a previous study that revealed reduced sensitivity in osteosarcoma cells after rhIL-6 treatment, we further found that p-FLNC can function as a lobaplatin resistance-related target in osteosarcoma cells. However, because a phosphoantibody was not available for all upregulated phosphoproteins, we did not detect additional phosphorylated proteins. We aimed to further verify other phosphoproteins after ensuring the reliability of the corresponding antibodies in order to reveal more possibilities related to drug resistance induced by IL-6 in osteosarcoma. Additionally, larger patient cohorts will be used to further evaluate this phosphoprotein and its kinase in lobaplatin-treated osteosarcoma specimens.
In this study, we conducted a comprehensive analysis of the phosphoprotein profiles of chemoresistant and chemosensitive osteosarcoma cells induced by rhIL-6. p-FLNC, its kinase AKT1, and MAPK signaling were identified to play pivotal roles in this process. Further western blot analysis and IHC verification proved that JNK and ERK signaling, which were activated in rhIL-6-induced chemoresistant osteosarcoma cells, were the main contributors, while p38 signaling was not involved. This study provides an experimental basis for chemoresistance to lobaplatin at the post-translational level and reveals potential molecular targets for increasing chemosensitivity in osteosarcoma.
Availability of data and materials
The MS proteomics data were deposited in the ProteomeXchange Consortium via the PRIDE  partner repository with the following ProteomeXchange dataset identifier: PXD025970.
Recombinant human IL-6
Kyoto Encyclopedia of Genes and Genomes
Higher energy collisional dissociation
Automatic gain control
KEGG Orthology and Links Annotation
Sodium dodecyl sulfate–polyacrylamide gel electrophoresis
Differentially expressed phosphoproteins
Tandem mass tag
P-extracellular signal-regulated kinase 1/2
Hematoxylin and eosin
Harrison DJ, Schwartz CL. Osteogenic sarcoma: systemic chemotherapy options for localized disease. Curr Treat Options Oncol. 2017;18(4):24.
Ritter J, Bielack SS. Osteosarcoma. Ann Oncol. 2010;21:vii320–5.
Munnink TO, van der Meer A, de Haan J, Touw D, van Kruchten M. Reversible impaired methotrexate clearance after platinum-based chemotherapy for osteosarcoma. Ther Drug Monit. 2019;41(6):693–5.
Unver N, McAllister F. IL-6 family cytokines: key inflammatory mediators as biomarkers and potential therapeutic targets. Cytokine Growth Factor Rev. 2018;41:10–7.
Ham IH, Oh HJ, Jin H, Bae CA, Jeon SM, Choi KS, Son SY, Han SU, Brekken RA, Lee D, et al. Targeting interleukin-6 as a strategy to overcome stroma-induced resistance to chemotherapy in gastric cancer. Mol Cancer. 2019;18(1):68.
Wang H, Li B, Yan K, Wu Y, Wen Y, Liu Y, Fan P, Ma Q. Protein and signaling pathway responses to rhIL-6 intervention before lobaplatin treatment in osteosarcoma cells. Front Oncol. 2021;11:602712.
Han X, Liu F, Zhang C, Ren Z, Li L, Wang G. SIAH1/ZEB1/IL-6 axis is involved in doxorubicin (Dox) resistance of osteosarcoma cells. Biol Chem. 2019; 400:(9):1241.
Park GB, Kim DJ, Kim YS, Lee HK, Kim CW, Hur DY. Silencing of galectin-3 represses osteosarcoma cell migration and invasion through inhibition of FAK/Src/Lyn activation and β-catenin expression and increases susceptibility to chemotherapeutic agents. Int J Oncol. 2015;46:185–94.
Bekker-Jensen DB, Bernhardt OM, Hogrebe A, Martinez-Val A, Verbeke L, Gandhi T, Kelstrup CD, Reiter L, Olsen JV. Rapid and site-specific deep phosphoproteome profiling by data-independent acquisition without the need for spectral libraries. Nat Commun. 2020;11:787.
Gotz S, Garcia-Gomez JM, Terol J, Williams TD, Nagaraj SH, Nueda MJ, Robles M, Talon M, Dopazo J, Conesa A. High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res. 2008;36(10):3420–35.
Quevillon E, Silventoinen V, Pillai S, Harte N, Mulder N, Apweiler R, Lopez R. InterProScan: protein domains identifier. Nucleic Acids Res. 2005;33:W116–20.
Kanehisa M, Sato Y, Morishima K. BlastKOALA and GhostKOALA: KEGG tools for functional characterization of genome and metagenome sequences. J Mol Biol. 2016;428(4):726–31.
Cheng A, Grant CE, Noble WS, Bailey TL. MoMo: discovery of statistically significant posttranslational modification motifs. Bioinformatics. 2019;35(16):2774–82.
Szklarczyk D, Franceschini A, Wyder S, Forslund K, Heller D, Huerta-Cepas J, Simonovic M, Roth A, Santos A, Tsafou KP, et al. STRING v10: protein–protein interaction networks, integrated over the tree of life. Nucleic Acids Res. 2015;43:D447–52.
Ramakrishnan VG, Miller KC, Macon EP, Kimlinger TK, Haug J, Kumar S, Gonsalves WI, Rajkumar SV, Kumar SK. Histone deacetylase inhibition in combination with MEK or BCL-2 inhibition in multiple myeloma. Hematologica. 2019;104(10):2061–74.
Rodel C, Martus P, Papadoupolos T, Fuzesi L, Klimpfinger M, Fietkau R, Liersch T, Hohenberger W, Raab R, Sauer R, et al. Prognostic significance of tumor regression after preoperative chemoradiotherapy for rectal cancer. J Clin Oncol. 2005;23(34):8688–96.
Pu YS, Hour TC, Chuang SE, Cheng AL, Lai MK, Kuo ML. Interleukin-6 is responsible for drug resistance and anti-apoptotic effects in prostatic cancer cells. Prostate. 2004;60(2):120–9.
Wang Y, Niu XL, Qu Y, Wu J, Zhu YQ, Sun WJ, Li LZ. Autocrine production of interleukin-6 confers cisplatin and paclitaxel resistance in ovarian cancer cells. Cancer Lett. 2010;295(1):110–23.
Xu X, Ye Y, Wang X, Lu B, Guo Z, Wu S. JMJD3-regulated expression of IL-6 is involved in the proliferation and chemosensitivity of acute myeloid leukemia cells. Biol Chem. 2021;402(7):815–24.
Singh V, Ram M, Kumar R, Prasad R, Roy BK, Singh KK. Phosphorylation: implications in cancer. Protein J. 2017;36(1):1–6.
Low TY, Mohtar MA, Lee PY, Omar N, Zhou H, Ye M. Widening the bottleneck of phosphoproteomics: evolving strategies for phosphopeptide enrichment. Mass Spectrom Rev. 2020;00:1–25.
Vidya MK, Kumar VG, Sejian V, Bagath M, Krishnan G, Bhatta R. Toll-like receptors: significance, ligands, signaling pathways, and functions in mammals. Int Rev Immunol. 2018;37(1):20–36.
Wu Q, Wu W, Fu B, Shi L, Wang X, Kuca K. JNK signaling in cancer cell survival. Med Res Rev. 2019;39(6):2082–104.
Zheng Y, Pan D. The hippo signaling pathway in development and disease. Dev Cell. 2019;50(3):264–82.
Nguyen CDK, Yi C. YAP/TAZ signaling and resistance to cancer therapy. Trends Cancer. 2019;5(5):283–96.
Asl ER, Amini M, Najafi S, Mansoori B, Mokhtarzadeh A, Mohammadi A, Lotfinejad P, Bagheri M, Shirjang S, Lotfi Z, et al. Interplay between MAPK/ERK signaling pathway and MicroRNAs: a crucial mechanism regulating cancer cell metabolism and tumor progression. Life Sci. 2021;278:119499.
Liu P, Lu Z, Liu L, Li R, Liang Z, Shen M, Xu H, Ren D, Ji M, Yuan S, et al. NOD-like receptor signaling in inflammation-associated cancers: from functions to targeted therapies. Phytomedicine. 2019;64:152925.
Snaka T, Fasel N. Behind the scenes: nod-like receptor X1 controls inflammation and metabolism. Front Cell Infect Microbiol. 2020;10:609812.
Jang DI, Lee AH, Shin HY, Song HR, Park JH, Kang TB, Lee SR, Yang SH. The role of tumor necrosis factor alpha (TNF-alpha) in autoimmune disease and current TNF-alpha inhibitors in therapeutics. Int J Mol Sci. 2021;22(5):2719.
Lee S, Rauch J, Kolch W. Targeting MAPK signaling in cancer: mechanisms of drug resistance and sensitivity. Int J Mol Sci. 2020;21(3):1102.
Hoang VT, Matossian MD, La J, Hoang K, Ucar DA, Elliott S, Burks HE, Wright TD, Patel S, Bhatt A, et al. Dual inhibition of MEK1/2 and MEK5 suppresses the EMT/migration axis in triple-negative breast cancer through FRA-1 regulation. J Cell Biochem. 2021;122:835.
Xu Y, Huang Y, Chen Y, Cao K, Liu Z, Wan Z, Liao Z, Li B, Cui J, Yang Y, et al. Grape seed proanthocyanidins play the roles of radioprotection on normal lung and radiosensitization on lung cancer via differential regulation of the MAPK signaling pathway. J Cancer. 2021;12(10):2844–54.
Khan S, Zaki H. Crosstalk between NLRP12 and JNK during hepatocellular carcinoma. Int J Mol Sci. 2020;21(2):496.
Ogura K, Yahiro K, Moss J. Cell death signaling pathway induced by cholix toxin, a cytotoxin and eEF2 ADP-ribosyltransferase produced by vibrio cholerae. Toxins (Basel). 2020;13(1):12.
Rosa JP, Raslova H, Bryckaert M. Filamin A: key actor in platelet biology. Blood. 2019;134(16):1279–88.
Kamil M, Shinsato Y, Higa N, Hirano T, Idogawa M, Takajo T, Minami K, Shimokawa M, Yamamoto M, Kawahara K, et al. High filamin-C expression predicts enhanced invasiveness and poor outcome in glioblastoma multiforme. Br J Cancer. 2019;120(8):819–26.
Yang B, Liu Y, Zhao J, Hei K, Zhuang H, Li Q, Wei W, Chen R, Zhang N, Li Y. Ectopic overexpression of filamin C scaffolds MEK1/2 and ERK1/2 to promote the progression of human hepatocellular carcinoma. Cancer Lett. 2017;388:167–76.
Kokate SB, Dixit P, Poirah I, Roy AD, Chakraborty D, Rout N, Singh SP, Ashktorab H, Smoot DT, Bhattacharyya A. Testin and filamin-C downregulation by acetylated Siah2 increases invasiveness of Helicobacter pylori-infected gastric cancer cells. Int J Biochem Cell Biol. 2018;103:14–24.
Wang ZK, Wang J, Liu J, Ying SH, Peng XJ, Feng MG. Proteomic and phosphoproteomic insights into a signaling hub role for cdc14 in asexual development and multiple stress responses in beauveria bassiana. PLoS ONE. 2016;11(4):e0153007.
Perez-Riverol Y, Csordas A, Bai J, Bernal-Llinares M, Hewapathirana S, Kundu DJ, Inuganti A, Griss J, Mayer G, Eisenacher M, et al. The PRIDE database and related tools and resources in 2019: improving support for quantification data. Nucleic Acids Res. 2019;47(D1):D442–50.
We thank Ms. Yingxue Wang for her technical assistance and contribution to amending the language of the manuscript.
Ethics approval and consent to participate
This investigation was performed in accordance with the Declaration of Helsinki. Written informed consent was provided by each included patient. This study was approved by the Institutional Review Board of Tangdu Hospital, Fourth Military Medical University (permission number: TDLL-202105-06).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
. Hierarchical clustering of the differentially expressed phosphoproteins in SaOS-2 osteosarcoma cells between the Lob group and the Con group. Each group contained three biological replicates. In total, 1,815 phosphoproteins with significantly differential expression (specifically, 874 upregulated and 941 downregulated phosphoproteins) were identified (fold change > 1.2, p < 0.05).
. Clinical characteristics of osteosarcoma patients treated with platinum-based chemotherapy.
. Clinical characteristics of osteosarcoma patients and the correlation with chemotherapy sensitivity.
. Quantification of the expression of FLNC, ERK1/2, JNK, and P38 by immunohistochemistry staining.
. Kaplan-Meier overall survival curves comparing osteosarcoma patients with high and low p-FLNC expression levels (n = 40, p < 0.05)
About this article
Cite this article
Zhang, R., Wang, H., Li, E. et al. Quantitative phosphoproteomic analysis reveals chemoresistance-related proteins and signaling pathways induced by rhIL-6 in human osteosarcoma cells. Cancer Cell Int 21, 581 (2021). https://doi.org/10.1186/s12935-021-02286-z
- rhIL-6 intervention
- MAPK signaling