Targeting sphingolipid metabolism with the sphingosine kinase inhibitor SKI-II overcomes hypoxia-induced chemotherapy resistance in glioblastoma cells: effects on cell death, self-renewal, and invasion
BMC Cancer volume 23, Article number: 762 (2023)
Glioblastoma patients commonly develop resistance to temozolomide chemotherapy. Hypoxia, which supports chemotherapy resistance, favors the expansion of glioblastoma stem cells (GSC), contributing to tumor relapse. Because of a deregulated sphingolipid metabolism, glioblastoma tissues contain high levels of the pro-survival sphingosine-1-phosphate and low levels of the pro-apoptotic ceramide. The latter can be metabolized to sphingosine-1-phosphate by sphingosine kinase (SK) 1 that is overexpressed in glioblastoma. The small molecule SKI-II inhibits SK and dihydroceramide desaturase 1, which converts dihydroceramide to ceramide. We previously reported that SKI-II combined with temozolomide induces caspase-dependent cell death, preceded by dihydrosphingolipids accumulation and autophagy in normoxia. In the present study, we investigated the effects of a low-dose combination of temozolomide and SKI-II under normoxia and hypoxia in glioblastoma cells and patient-derived GCSs.
Drug synergism was analyzed with the Chou-Talalay Combination Index method. Dose–effect curves of each drug were determined with the Sulforhodamine B colorimetric assay. Cell death mechanisms and autophagy were analyzed by immunofluorescence, flow cytometry and western blot; sphingolipid metabolism alterations by mass spectrometry and gene expression analysis. GSCs self-renewal capacity was determined using extreme limiting dilution assays and invasion of glioblastoma cells using a 3D spheroid model.
Temozolomide resistance of glioblastoma cells was increased under hypoxia. However, combination of temozolomide (48 µM) with SKI-II (2.66 µM) synergistically inhibited glioblastoma cell growth and potentiated glioblastoma cell death relative to single treatments under hypoxia. This low-dose combination did not induce dihydrosphingolipids accumulation, but a decrease in ceramide and its metabolites. It induced oxidative and endoplasmic reticulum stress and triggered caspase-independent cell death. It impaired the self-renewal capacity of temozolomide-resistant GSCs, especially under hypoxia. Furthermore, it decreased invasion of glioblastoma cell spheroids.
This in vitro study provides novel insights on the links between sphingolipid metabolism and invasion, a hallmark of cancer, and cancer stem cells, key drivers of cancer. It demonstrates the therapeutic potential of approaches that combine modulation of sphingolipid metabolism with first-line agent temozolomide in overcoming tumor growth and relapse by reducing hypoxia-induced resistance to chemotherapy and by targeting both differentiated and stem glioblastoma cells.
Glioblastoma (GB) is the most common and malignant form of primary brain tumors. The standard treatment for newly diagnosed cases is surgical resection followed by radiotherapy plus concomitant and adjuvant temozolomide (TMZ) chemotherapy . TMZ primarily induces cytotoxicity through methylation of guanine at position O6 (O6-methylguanine) , triggering G2/M cell cycle arrest and, eventually, apoptosis. However, the enzyme O6-methylguanine-DNA methyltransferase (MGMT) is able to repair this DNA adduct . Intrinsic and acquired resistance to TMZ contributes to inevitable tumor recurrence and a short median survival of patients . Currently, there is no established treatment for recurrent GB . Tumor hypoxia, a hallmark of GB and other solid tumors, promotes tumor cell invasion into the brain parenchyma and contributes to therapy resistance . Moreover, hypoxia and TMZ favor the expansion of the glioblastoma stem cell (GSC) population, which are undifferentiated cancer cells with unlimited self-renewal capability and propensity to promote intratumoral heterogeneity [7, 8]. GSCs are considered important drivers of tumor relapse after initial treatment . Thus, prevention of GB recurrence requires the combined targeting of differentiated tumor cells and GSCs. Sphingolipids are involved in a variety of biological functions. Regulation of their metabolism is thus essential for cell homeostasis and fate and, as such, it has evolved as a very attractive therapeutic target . Sphingosine kinases (SK) 1 and 2, are enzymes controlling the balance between the sphingolipids sphingosine-1-phosphate (S1P) and ceramide. Ceramide is catabolized by ceramidases into sphingosine which can be further phosphorylated to generate S1P . S1P can be exported into the extracellular milieu, stimulate five G-protein-coupled receptors (S1PR1-5) and regulate intracellular signaling pathways involved in proliferation, invasion, angiogenesis, and stem cell properties . In contrast, ceramide is known for its pro-apoptotic activities . Ceramide can be generated from sphingosine via ceramide synthases, but also synthetized in the de novo pathway from dihydroceramide via dihydroceramide desaturase 1 (DES1), or generated from the hydrolysis of sphingomyelin via sphingomyelinase. Glioblastoma tissue specimens contain higher S1P and lower ceramide levels relative to healthy tissue, and high SK1 expression correlated with a shorter survival time of GB patients . Various chemotherapeutic agents have been shown to trigger tumor cell death by favoring the production of ceramide . However, in the case of glioblastoma, this increased ceramide production might be counteracted by SK overexpression, leading to conversion of ceramide into S1P. Blocking SK might prevent this conversion and support TMZ efficacy. The dual sphingosine kinase inhibitor SKI-II is orally bioavailable and showed significant in vivo antitumor activity without displaying toxicity . Although first described as a sphingosine kinase inhibitor, it was later demonstrated to have off target effects on DES1 . Reported inhibition constants (Ki) indicate a stronger effect of SKI-II on DES1 (Ki = 0.3 µM) than on SK1 (Ki = 16 µM) . We have previously shown that 10 µM SKI-II combined with 500 µM TMZ induced caspase-dependent cell death in GB cells, preceded by accumulation of dihydrosphingosine and dihydroceramide, oxidative stress, endoplasmic reticulum (ER) stress, and autophagy. In that experimental setup, dihydrosphingolipids, and not ceramide, were the pro-apoptotic sphingolipids, reflecting the inhibitory activities of SKI-II towards both SK1 and DES1 . Other studies have as well demonstrated dihydroceramide to be involved in ER stress, autophagy, and proliferation .
The raising interest in modulating the sphingolipid metabolism and specifically the sphingosine/S1P axis for therapeutic purposes is reflected by the increased number of SK1 and SK2 inhibitors developed in the last years. Whereas some inhibitors are tested in clinical trials (e.g. the SK2 inhibitor ABC294640; see reference in , most SK inhibitors (including SKI-II and more recently described SK1 inhibitors) are still in the preclinical testing phase . This slow translation to the clinic indicates that we need to deepen our knowledge of the effects these inhibitors exert on biological processes and on sphingolipid metabolism for a better evaluation of their therapeutic value and affected pathway(s). In the present study, we further investigated the therapeutic potential of SKI-II combined with TMZ at clinically relevant doses, based on the concentrations detected in the serum and tumor of patients [21, 22], and we compared the effects of the combination under normoxic (21% O2) and hypoxic (3% O2) conditions. We report that a low-dose combination (48 µM TMZ + 2.66 µM SKI-II) has potent anti-tumor activities under hypoxia, such as cell growth inhibition and cell death induction, self-renewal impairment of TMZ-resistant GSCs, and reduction of the invasion capacity of recurrent GB cell-derived spheroids.
Temozolomide, SKI-II, Camptothecin, Z-VAD-FMK, Sulforhodamine B, FCCP and FeTPPS were purchased from Sigma-Aldrich. All drugs were dissolved in DMSO as stock solutions and kept at -20°C. Working dilutions of all drugs were performed in culture media immediately before each experiment. AnnexinV-FITC and 10X Annexin Binding Buffer were purchased from Dianova. MitoStatus TMRE was purchased from BD Pharmingen.
The low-passage-number human glioblastoma cell line NCH82 was generated at the Department of Neurosurgery, Heidelberg University Hospital (Heidelberg, Germany) . NCH82 cells were cultured in complete DMEM (cDMEM) [Dulbecco's Modified Eagle Medium (Sigma),10% (vol/vol) fetal calf serum (Biochrom), 2 mM L-glutamine (Life Technologies), gentamicin (50 μg/ml) (GIBCO)]. NCH82 cells were kept in culture for no longer than six passages after thawing.
The human glioblastoma stem-like cell (GSC) line No. 1080 was generated at the University Medical Center Göttingen (Göttingen, Germany) . GSCs were cultured as spheres in NeuroBasal medium (GIBCO) supplemented with the B27 component (Invitrogen Life technologies), fibroblast growth factor (FGF) (Peprotech, 10 ng/ml) and epidermal growth factor (EGF) (Peprotech, 20 ng/ml). 1080 GSCs were selected with 100 µM TMZ (hereafter referred as TMZ-1080), and the respective DMSO control (hereafter referred as DMSO-1080).
The human glioblastoma cell line U3054MG (mesenchymal subtype) was acquired from the HGCC (Human Glioblastoma Cell Culture Resource) where it was established from a GB that recurred after the patient had undergone surgery and radiotherapy (50.4 Gy in total) . Cells were cultured in laminin-coated Primaria dishes (Corning) in Neurobasal and DMEM/F12 media (1:1 mix) containing N2 and B27 supplements (Invitrogen) and human recombinant FGF and EGF (10 ng/ml, Peprotech).
All cells were cultured at 37°C, 5% CO2 under normoxic (21% O2) conditions. All treatments were performed in culture medium containing 5% FCS.
Cell cultures were regularly tested for mycoplasma contamination (MycoAlert™, Lonza Bioscience). The NCH82 and 1080 cell lines were authenticated via Single Nucleotide Polymorphism (SNP)-profiling.
Sulforhodamine B assay
The sulforhodamine B (SRB) colorimetric assay was used to determine the dose–effect curves of each drug alone and assess the effect of drug combinations. NCH82 cells were seeded in 96-well plates, allowed to attach for 24 h and then treated with each drug and vehicle control (DMSO) in quadruplicates. The cells were further incubated for 5 days at 21% and 3% O2 in dedicated incubators. The SRB assay was performed as described . The absorbance (OD) of each well was read at 510 nm on a microplate reader (TECAN Infinite M200 PRO). Background absorbance was subtracted from all samples. Percentage of control cell growth was calculated using the formula: % of control cell growth = [(mean OD sample at d5 – mean OD at d0) / (mean OD control at d5 – mean OD at d0)] × 100.
Propidium iodide and Annexin V staining of NCH82 cells and GSC cells
Flow cytometry was used to assess drug-induced cell death by performing staining with propidium iodide (PI) and Annexin V-FITC (Dianova). NCH82 cells were seeded in 12-well plates and treated accordingly in duplicates 24 h later. After 3 and 5 days of incubation at 21% and 3% O2, adherent and floating cells were collected and washed twice with 1X Annexin V binding buffer (Dianova). Then, 5 µl of Annexin V-FITC and 1 µg/ml of PI were added to 100 µl of cell suspension and incubated for 15 min in the dark at room temperature. Finally, cells were analyzed with the Attune NxT Flow Cytometer (Thermo Fisher Scientific) at a sample flow rate of 100 µl/min. Flow cytometry data were analyzed with the software FlowJo. The following gating strategy was applied: doublets discrimination (SSC-A vs SSC-H); selection of the population of interest/debris elimination (SSC-H vs FSC-H) and application of the two-color gating/quadrant (BL2-H vs BL1-H). Spheres of the DMSO-1080 and TMZ-1080 GSC cells were dissociated using Trypsin–EDTA 0,05% (Thermo Fischer Scientific) diluted at a 1:1 ratio in NeuroBasal medium, seeded in 12-well plates and treated accordingly in duplicates 24 h later. After 5 days of incubation at 21% O2 and 3% O2, adherent and floating cells were collected and processed for staining and analysis at the flow cytometer as described for the NCH82 cells.
Drug combination analysis
The combination index (CI) theorem of Chou and Talalay was used to quantitatively determine synergism between combinations of TMZ and SKI-II [27,28,29]. A combination is considered to be synergistic when CI is below 1, antagonistic when CI is greater than 1, and to have an additive effect when CI equals 1. The CI was calculated via the third-generation software CompuSyn . First, dose–effect curves of each drug alone were determined by exposing cells to twofold serial dilutions of TMZ and SKI-II for 5 days at 21% and 3% O2. The percentage of cell growth was measured using the SRB assay as previously described. This allowed to determine the parameters Dm and m for each drug, which are required for CI calculation. Dm is the dose required to achieve a 50% effect level (ED50), and m depicts the shape of the dose–effect curve (m = 1 hyperbolic; m > 1 sigmoidal; m < 1 negative sigmoidal). Additionally, the parameter r (linear correlation coefficient of the median effect plot) was obtained and refers to the conformity of the data to the mass-action law (for in vitro studies r > 0.95). Next, combinations of TMZ and SKI-II at non-constant ratios were tested for 5 days at 21% and 3% O2. Two-fold serial dilutions below and above the determined ED50 (21% O2) were performed for SKI-II, and for TMZ only below the ED50, since values above the ED50 (21% and 3% O2) are clinically irrelevant.
In vitro extreme limiting dilution assay
GSCs (DMSO-1080 and TMZ-1080) were seeded in 24-well plates in complete NeuroBasal medium at a density of 10 cells/ml and, subsequently, serial diluted until the theoretical cell density of 0.3125 cells/ml. The number of wells used per cell density for each condition ([DMSO, TMZ, SKI-II and (TMZ + SKI-II)] in every independent experiment was the following: 6 wells for 10 cells/ml; 12 wells for 5 cells/ml; 18 wells for 2.5 cells/ml; 18 cells for 1.25 cells/ml; 12 wells for 0.625 cells/ml.; and 6 wells for 0.3125 cells/ml. Cells were treated 24 h later and further incubated for 4 weeks at 21% and 3% O2. After the incubation period, the number of wells containing at least one sphere (positive cultures) was registered for each dose (number of cells/well). The stem cell frequency was determined using the webtool ELDA provided by the Walter and Eliza Hall Institute Bioinformatics Division (https://bioinf.wehi.edu.au/software/elda/) .
Spheroid invasion assay
U3054 spheroids were generated for 48 h in ultra-low attachment 96-well plates (CellCarrier Spheroid ULA 96-well Microplates, PerkinElmer) at a density of 8000 cells/40 µl/well/spheroid. Spheroids of approximately 400 µm diameter were embedded in a collagen mixture composed of 10X Minimum Essential Medium, 0.1 M NaOH, 3 mg/ml of PureCol® Bovine Collagen solution type I (Advanced Biomatrix) and respective treatment. After 1 h incubation at 37 °C, complete DMEM (5% FCS) containing the respective treatment was added on the surface of the polymerized collagen mixture and light microscopy images were collected from each spheroid (day 0). Embedded spheroids were further incubated at 21% O2 or 3% O2 and images were collected daily. These images were uploaded in the software Spheroid Analyzer from CLADIAC  which automatically determined the invaded area (including the core of the spheroid) for each day and calculated its relative size compared to day 0.
Detection of activated caspase-3 in treated NCH82 cells and GSC cells
NCH82 cells seeded in 12-well plates were treated accordingly in duplicates 24 h later. After drug incubation time at 21% and 3% O2, adherent and floating cells were collected and washed twice with cold 1X PBS, resuspended in BD Cytofix/Cytoperm solution (BD Biosciences) and incubated 20 min on ice. Then, cells were washed twice with 1X BD Perm/Wash buffer (BD Biosciences) at room temperature and incubated for 30 min with 30 µl of FITC Rabbit Anti-Active Caspase-3 antibody (BD Biosciences) in 100 µl of 1X BD Perm/Wash buffer at room temperature. Finally, each sample was washed and resuspended in 1X BD Perm/Wash buffer and analyzed by flow cytometry. Flow cytometry data were analyzed via the software FlowJo. The following gating strategy was applied: doublets discrimination (SSC-A vs SSC-H); selection of the population of interest/debris elimination (SSC-H vs FSC-H) and application of one-color gating (BL1-H vs FSC-H).
Spheres of DMSO-1080 and TMZ-1080 GSC cells were dissociated using Trypsin–EDTA 0,05% diluted at a 1:1 ratio in NeuroBasal medium, plated on ornithin-coated glass coverslips and treated accordingly in triplicates 24 h later. After 5 days of incubation at 21% and at 3% O2, cells were fixed with 4% paraformaldehyde, washed with PBS, incubated in blocking solution (0.1% Triton X-100, 1% bovine serum albumin) and stained overnight at + 4 °C with a mix of antibodies specific for cleaved caspase-3 (Cell Signaling) and human nestin (Abcam). Secondary antibodies were goat α-mouse Alexa Fluor 488 (Invitrogen) and goat α-rabbit Alexa Fluor 555 (Invitrogen). Dapi staining was performed according to standard conditions. Stained cells were analyzed by using the BZ-X fluorescence microscope (Keyence Germany GbmH). For each condition (triplicate), 700 to 1100 cells were counted from a total of 7 to 9 regions. Percentage of cells positive for cleaved caspase-3 was determined by using the Image Tool software.
Measurement of mitochondrial membrane potential
NCH82 cells were seeded in 12-well plates and treated accordingly in triplicates 24 h later. Cells were treated with (TMZ + SKI-II) for 24 and 48 h, or with 50 µM FCCP (carbonyl cyanide 4-(trifluoromethoxy) phenylhydrazone) for 60 min. FCCP was used to dissipate the mitochondrial membrane potential. After incubation time under 21% and 3% O2, the supernatant was removed and replaced with a solution of TMRE (1.2 µM) in Stain Buffer (BD Biosciences). After 30 min incubation at 37°C, the solution was removed and the wells were washed twice with Stain Buffer. Cells were detached with Accutase and collected into polypropylene tubes. Finally, cells were washed twice with Stain Buffer at room temperature and analyzed by flow cytometry. The median fluorescence intensity was calculated with the software FlowJo.
Western blot analysis
After treatment, adherent and floating cells were harvested and washed twice with cold PBS. The cell pellets were lysed in RIPA buffer (150 mM NaCl, 0.1% Triton X-100, 0.5% sodium deoxycholate, 0.1% SDS, 20 mM Tris–HCl pH 7.5) supplemented with cOmplete™ EDTA-free protease inhibitor cocktail (Roche) for 30 min on ice. Lysates were centrifuged at 13,000 RPM for 20 min at 4 °C. Protein concentration was measured with Pierce™ BCA protein assay kit (Thermo Scientific). Samples of lysates were prepared in Laemmli buffer containing 100 mM DTT, boiled for 5 min at 95 °C and resolved in a 12% polyacrylamide gel. The following amount of protein was loaded per lane: 20 µg for Caspase-3, 25 µg for LC3 and p62 and 35 µg for BiP analyses. Transfer was performed onto nitrocellulose membrane via a semi-dry blotting system. Membranes were blocked in 5% nonfat milk powder in TBS-T buffer (Tris-base, NaCl and 0.1% Tween 20) for 1 h at room temperature. Following blocking, the membranes were cut in strips and strips were labeled overnight at 4 °C with the following antibodies against: Cleaved Caspase-3 (Cell Signaling, 1:1000; rabbit), LC3B (Cell Signaling, 1:1000; rabbit), SQSTM1/p62 (Cell Signaling, 1:1000; rabbit), BiP (Cell Signaling, 1:1000; rabbit), GAPDH (BioLegend, 1:1000; mouse), and GAPDH (Thermo Scientific, 1:1000; rabbit). After washing in TBS-T, membranes were incubated with HRP-conjugated secondary antibodies (anti-mouse IgG, Cell Signaling; or anti-rabbit IgG, Cell Signaling) for 1 h at room temperature, washed again and incubated with chemiluminescence reagent (Western Lightning Plus-ECL, Enhanced Chemiluminescence Substrate, PerkinElmer) for 1 min. Detection was performed with the chemiluminescent imaging system FusionCapt Advance FX7.
Sphingolipids were extracted from cell pellets using methyl-tert-butyl ether (MTBE)-based liquid–liquid extraction method (LLE). Cells were pelleted in Precellys tubes (Peqlab) to which steel balls were manually added. 1000 µl of MTBE/methanol (10:3, v/v), which served as an extraction solvent, were added to the tubes and then 250 µl of ice-cold 0.1% formic acid containing 5 µM tetrahydrolipstatin/3′-(aminocarbonyl) [1,1′-biphenyl]3-yl)-cyclohexylcarbamate and 10 µg/ml butyl hydroxytoluene as homogenization solvent. A 10 µl methanolic solution of internal standards: Sphingosine d17:1, Sphingosine-1P d17:1, Ceramide d18:1/17:0, Ceramide-1P d18:1/12:0, and Sphingomyelin d18:1/12:0 was then added to the tubes. The concentration of internal standards in the spike solution was set so as to result in a target concentration in the 100 µl final extracts of: 200 ng/ml for Sphingosine d17:1; 200 ng/ml for Sphingosine-1P d17:1; 500 ng/ml for Ceramide d18:1/17:0; 500 ng/ml for Ceramide-1P d18:1/12:0 and 100 ng/ml for Sphingomyelin d18:1/12:0. Samples were then homogenized for 20 s at 6000 rpm. Homogenates were transferred into 1.5 ml Eppendorf tubes, vortexed for 30 s at 4°C and centrifuged for 10 min at 13,000 rpm. Upper organic phase was recovered in new 1.5 ml tubes, evaporated to dryness under a gentle stream of nitrogen, reconstituted in 90 µl methanol and stored at -20°C till further analysis. Lower aqueous phase was used for protein content determination using the BCA assay. Lipids were analyzed by liquid chromatography/multiple reaction monitoring (LC/MRM) using QTRAP 5500 (AB Sciex) operating in positive/ negative ion mode switching. The ionization and detection parameters as well as LC solvents used for this analysis are as described in  and in  and are depicted in Additional file 1. Calibrants and MRM transitions for sphingomyelins, sphingosine and ceramide 1P species are as described in . For the quantification of Ceramide d18:1/24:1 and Ceramide d18:1/16:0 in cell pellets, Ceramide d18:1/17:0 served as internal standard and Ceramide d18:1/16:0 was used as calibration standard. The MRM transitions for detection/quantification of ceramides were manually inferred by direct infusion analysis. The transitions and MRM parameters are listed in Additional file 2. For analysis an aliquot of 27 µl methanolic solution of lipid extract was mixed with 3 µl water and transferred in vials for liquid-chromatography/multiple reaction monitoring (LC/MRM) analysis. Lipid species were quantified using the Multiquant software version 3.0.3. The determined lipid concentrations were normalized to protein content of the cell pellets.
Real-time quantitative PCR
Isolation of total RNA was performed using RNeasy Mini Kit (Qiagen) according to manufacturers’ instructions. After treatment, cells were collected, washed with 1xPBS and lysed in Buffer RLT (containing 1% ß-mercaptoethanol). All following steps were conducted as described in the manufacturers’ protocol. RNA concentration and quality were determined using a Nanodrop 2200 (ThermoFisher). Only samples showing a 260/280 nm ratio between 1.8 and 2.1 were selected for cDNA transcription, which was performed with the Omniscript RT Kit (Qiagen) and random hexamers (Life Technologies). Quantitative PCR (qPCR) analysis was done using TaqMan® primers and a StepOnePlus System (Applied Biosystems). Briefly, for each well of the 96-well qPCR plate (Sarstedt), 10 µl of TaqMan™ Universal PCR Master Mix (ThermoFisher) were mixed with 10 ng cDNA and 1 µl of the appropriate primer (ThermoFisher). All measurements were performed using three technical replicates. Relative quantification (RQ) of gene expression were determined using the 2−ΔΔCt method. Primer IDs: DDIT3 (Hs00358796_g1), GAPDH (Hs02758991_g1); DEGS1 (Hs00186447_m1); SK1 (Hs00184211_m1); SK2 (Hs01016543_g1). Expression of GAPDH, used as a reference gene, was similar in untreated cells under both 21% and 3% oxygen conditions (Ct values of 20,743 and 20,749 respectively).
Transmission Electron Microscopy (TEM)
NCH82 cells were cryo-fixed within a few milliseconds at a pressure of 2000 bar under liquid nitrogen using a high-pressure freezer Compact 1 (Wohlwend GmbH). Freeze-substitution was conducted using a Leica EM AFS 2 device (Leica Microsystems). Here, the substitution/staining medium (acetone p.a., 0.2% osmium tetroxide, 0.1% uranylacetate and 5% water) was pre-cooled to -90 °C before samples were added. Finally, the samples were embedded in EPON 812 and sectioned at room temperature using a diamond knife. Examination of the thin sections was conducted with a FEI Tecnai F20 transmission electron microscope (FEI) operated at an acceleration voltage of 200 kV. Conventional bright field images were acquired using a Gatan US1000 slow scan CCD camera (Gatan Inc.).
Data were expressed as the mean plus standard deviation (SD) of at least two independent experiments. Statistical analyses were performed using GraphPad Prism software (version 9.2.0). Mean comparison was performed using unpaired t test or ANOVA (one-way or two-way). Tukey and Šídák tests were used for multiple comparison analyses. The used statistical method is referred in the legend of the respective figure. Statistical significance was defined by an alpha of 0.05. * p ≤ 0.05; ** p ≤ 0.01; *** p ≤ 0.001; **** p ≤ 0.0001.
TMZ synergizes with SKI-II to suppress glioblastoma cell growth under normoxia and hypoxia
The Chou-Talalay combination index method [25,26,27] was used to find the combination of TMZ and SKI-II with the strongest degree of synergism at the highest effect level towards the human glioblastoma NCH82 cells. First, dose–effect curves of each drug alone were determined under normoxia (21% O2) and hypoxia (3% O2) by exposing NCH82 cells to several twofold serial dilutions of TMZ and SKI-II (Fig. 1 A, B). Afterwards, the median effect dose (ED50 or Dm), the m value (shape of the dose–effect curve) and the r value (conformity of the data to the mass-action law) were calculated using the software CompuSyn. These parameters are required to quantify synergism of drug combinations. The calculated ED50 for TMZ was 96 µM under normoxia and 428 µM under hypoxia, demonstrating a high resistance of the NCH82 cells to TMZ under low oxygen. On the contrary, SKI-II was as efficient under normoxia as under hypoxia in inhibiting cell growth (ED50 of approximately 1.3 μM in both oxygen conditions) (Additional file 3).
To evaluate the type and degree of interactions between SKI-II and TMZ, we next undertook to calculate the combination index (CI) for 15 combinations, where CI > 1 indicates antagonism, CI = 1 indicates an additive effect, and CI < 1 indicates synergism. SKI-II and TMZ were combined at a non-constant ratio and the concentrations of the 15 combinations were based on twofold serial dilutions above and below the ED50 of each drug (Additional file 4). TMZ concentrations were kept under clinically relevant ranges (3–100 µM; [21, 22]) and, therefore, not tested above 96 µM in combination.
All 15 combinations tested under normoxia had a CI value below 1 and are, accordingly, synergistic (Table 1). Under hypoxia some combinations became nearly additive (combinations A1 and C1-3) or presented a diminished degree of synergism (combinations B1-3) and their effect level was reduced. Interestingly, combinations with higher SKI-II concentrations, such as the combinations D1-3 (2 × ED50SKI-II) and E1-3 (4 × ED50SKI-II), presented a constant degree of synergism and similar effect level between oxygen conditions (Fig. 1 C, D; Table 1).
The synergistic drug combination consisting of 48 μM TMZ and 2.66 μM SKI-II (D2) was able to inhibit cell growth by almost 80% under normoxia and hypoxia (Table 1). The dose-reduction index (DRI), i.e. the measure of how many folds of dose-reduction for each drug at a given effect are allowed in synergistic combination, was obtained from the CompuSyn software (Additional file 5). The DRI of TMZ for this combination was of 84 under normoxia and 252 under hypoxia. Thus, the synergistic (TMZ + SKI-II) combination was able to reduce the needed concentration of TMZ by 252-fold under hypoxia. The DRI of SKI-II was as well favorable (DRI > 1) and comparable under both oxygen conditions (2.2 in normoxia and 1.9 in hypoxia). Combinations with a higher degree of synergism at a higher effect level are the most relevant for cancer therapy . For this reason, and in order to combine a clinically relevant dose of TMZ with the lowest concentration of SKI-II, we decided to use the D2 (TMZ + SKI-II) combination for further studies.
The combination of TMZ and SKI-II potentiates glioblastoma cell death under normoxia and hypoxia
To quantify cell death induced by (TMZ + SKI-II) combination and each drug alone, NCH82 cells were treated for 3 and 5 days under normoxia and hypoxia. After treatment, cells were stained with annexin-V-FITC (AV-FITC) and propidium iodide (PI) and, subsequently, analyzed by flow cytometry (Fig. 2 A). Under normoxia, single treatment with TMZ or SKI-II induced similar levels of cell death (approximately 30% of AV-FITC+/PI− and AV-FITC+/PI+ cells) after 3 and 5 days of treatment, whereas the level of cell death induced by the combination was significantly enhanced by about 30% after 5 days (Fig. 2 B, C). Under hypoxia, the combination and single treatments were similarly cytotoxic after 3 days, and this cytotoxicity increased after 5 days. The combination significantly potentiated cell death in comparison with TMZ alone after 5 days of treatment (Fig. 2 D, E).
Altogether, the present data show that treatment with the (TMZ + SKI-II) combination was, both under normoxia and hypoxia, more efficient in killing the glioblastoma NCH82 cells than single TMZ or SKI-II treatment, and this effect was time-dependent.
The combination induces cell death without caspase-3 activation or mitochondrial membrane potential disruption
Apoptosis can occur via the extrinsic and intrinsic signaling pathways. Both pathways culminate in the activation of effector caspases, such as caspase-3 . Addition of the pan-caspase inhibitor Z-VAD-FMK (zVAD) did not rescue the NCH82 cells from death induced by the combination (Fig. 2 B-E). To confirm this observation, we analyzed the effect of the combination on caspase-3 activation over 72 h of treatment under hypoxia and normoxia. As shown in Fig. 3, cleaved caspase-3 was not detected in cells treated with the combination after 48 h (Fig. 3 A, B, Bi; Additional file 6) and was expressed in approximately 5% or less of the cells analyzed after 72 h of treatment. This indicates that the majority of cells dies via a caspase-independent cell death pathway (Fig. 3 B, Bi). Note that NCH82 cells do activate caspase-3, as shown by the detection of the cleaved caspase-3 after treatment with camptothecin (CPT), a well-known apoptosis inducer (Fig. 3 A, B, Bi; Additional file 6). Altogether, these data suggest that the combination induces a mode of cell death independent of caspase activation.
During apoptosis, mitochondrial outer membrane permeabilization leads to a loss of mitochondrial membrane potential (MMP) – an event occurring before caspase activation. To analyze whether the (TMZ + SKI-II) combination causes alterations in the MMP, we stained treated NCH82 cells with tetramethylrhodamine ethyl ester (TMRE), a fluorescent dye only able to bind mitochondria with an intact MMP. Fluorescence microscopy showed that treatment with (TMZ + SKI-II) did not cause the loss of TMRE fluorescence which was observed in the FCCP-treated cells—a positive control for the dissipation of the MMP (Fig. 3 C). Quantification via flow cytometry confirmed that FCCP treatment caused a significant decrease of approximately 50% of the median fluorescence intensity (MFI) of TMRE. The MFI of cells treated with (TMZ + SKI-II) did not significantly differ from the MFI of control cells (DMSO condition) at 24 h and 48 h, except for an increase observed at 48 h under normoxia (Fig. 3 D). This increased fluorescence, that could be interpreted as an hyperpolarization, is most likely to result from an increase in cell death at late time points leading to an increase in the dye-to-living cells ratio, as described in . These results indicate that the combination does not cause mitochondrial depolarization under both oxygen conditions, and therefore, mitochondria do not seem to have a decisive role in (TMZ + SKI-II)-induced cell death.
The combination of TMZ and SKI-II does not affect autophagy in normoxia and hypoxia
We and others previously reported that TMZ and SKI-II are able to induce autophagy [16, 18, 37]. Therefore, autophagy could have an important role in, or at least, accompany (TMZ + SKI-II)-induced cell death. To monitor autophagic flux, we analyzed LC3-II and p62 levels in the presence and absence of bafilomycin A1 (BA1). BA1 is a specific inhibitor of the vacuolar type H+-ATPase and inhibits the acidification of organelles containing this enzyme (lysosomes and endosomes). Thus, it blocks the process of autophagy at a late stage. The use of BA1 supports the interpretation of LC3 immunoblotting, since an increase in LC3-II levels may represent an increase in autophagy or a block in the autophagy flux .
NCH82 cells were treated with the combination and DMSO control for 24 to 72 h under normoxia and hypoxia. Under normoxia, p62 and LC3-II protein levels remained similar between both conditions over time (Additional file 7 A, Ai, Aii; Additional file 8). The autophagic flux in (TMZ + SKI-II)-treated cells decreased from 24 to 48 h, but there were no differences in comparison with control cells (Additional file 7 Aiii). Under hypoxia, p62 protein levels were not significantly affected by the combination whereas LC3-II levels at 72 h were elevated in (TMZ + SKI-II)-treated cells in comparison with the DMSO condition (Additional file 7 B, Bi, Bii; Additional file 8). Indeed, treatment with the combination tended to decrease autophagic flux over time—however, this lacks statistical significance (Additional file 7 Biii). This trend was not observed in DMSO-treated cells, but as in normoxia, there were no statistically significant differences in the autophagic flux in comparison with the (TMZ + SKI-II)-treated cells. These results suggest that TMZ combined with SKI-II at low-doses does not significantly affect autophagic flux and, therefore, autophagy does not seem to be an essential process in cell death.
SKI-II alone and in combination increases cytoplasmic vacuolization
Microscopy observations indicated that SKI-II alone or in combination with TMZ induces cytoplasmic vacuolization in NCH82 cells (Fig. 4 A). In order to try and quantify this increase, we used flow cytometry to take advantage of the side scatter (SSC) parameter that is a measure of the cell granularity. The SSC of (TMZ + SKI-II)-treated cells was significantly increased in comparison to that of control cells (DMSO vehicle) under normoxia and hypoxia. This increase was evident after 24 h and aggravated over time. It should be noted that the formation of cytoplasmic vacuoles was observed at the microscope as early as 8 h after treatment (data not shown). SSC values observed in combination-treated cells were similar in both oxygen conditions (Fig. 4 B, C). The size of the cells, measured via the forward scatter (FSC, Fig. 4 D, E), was slightly increased after (TMZ + SKII) treatment under normoxia at later time points. These differences were not observed in hypoxia. Altogether these results indicate cytoplasmic vacuolization to be an early effect of sphingosine kinase inhibitor alone and in combination with TMZ.
The combination of TMZ and SKI-II induces ER stress and the Unfolded Protein Response
Based on our previous observations and reported effects of SKI-II on ER dilation and expansion , we further examined the cytoplasmic vacuolization in NCH82 treated cells by transmission electron microscopy (TEM). In control cells (DMSO) and in TMZ-treated cells, the ER presented a normal tubular morphology (Fig. 5 A, B). In cells treated with SKI-II alone and in combination with TMZ, we observed an enlargement but also fragmentation of the rough ER (vesicles with ribosomes) (Fig. 5 C, D). Thus, the cytoplasmic vacuolization observed by light microscopy could represent an expanded and dilated ER. Dilation of the nuclear envelope was also observed in cells treated with SKI-II and with the combination. This is not surprising, since the outer nuclear membrane is directly continuous with the rough ER. These ultrastructural changes were observed under both oxygen conditions (21% O2 not shown). Untreated NCH82 cells commonly presented disarrangement and distortion of mitochondrial cristae. Such abnormalities have been previously described in GB .
These morphological changes suggest that, similar to cells treated with a higher dose combination , cells treated with the (TMZ + SKI-II) combination used in the present study underwent ER stress. To confirm this, we analyzed the protein levels of the ER chaperone BiP (Binding immunoglobulin Protein) also known as GRP78 (glucose regulatory protein 78), and the expression of Unfolded Protein Response (UPR) target gene DDIT3 coding for CHOP (C/EBP Homologous Protein). Under stress conditions, BiP dissociates from ER transmembrane proteins resulting in the activation of the UPR. Subsequently, BiP expression is upregulated to support the re-establishment of ER homeostasis. Prolonged UPR activation leads to the induction of the pro-apoptotic transcription factor CHOP. Treatment of NCH82 cells up to 48 h with the combination resulted in upregulation of BiP expression levels. Co-treatment with Bafilomycin A1 did not affect BiP expression, suggesting that autophagy does not have a role in the mitigation of ER stress (Fig. 6 A-Ai; Additional file 9). In contrast to TMZ treatment, SKI-II alone or in combination elevated DDIT3 expression levels by about sixfold, showing that ER stress levels are increased due to SKI-II treatment solely (Fig. 6 B). Hypoxia does not seem to aggravate ER stress, since BiP and DDIT3 levels were similar in both oxygen conditions. Altogether, these data indicate that SKI-II alone or in combination with TMZ induces ER stress in NCH82 cells.
ER stress and cytotoxicity induced by the combination of TMZ and SKI-II are peroxynitrite-dependent
ER stress and the generation of reactive oxygen and nitrogen species are linked events . Peroxynitrite, a short-lived oxidant species produced by the reaction of nitric oxide with superoxide radicals, can induce ER stress via depletion of ER-Ca2+ . We previously reported that FeTPPS, a decomposition catalyst of peroxynitrite, protects NCH82 cells from ER stress and from death induced by the high-dose (TMZ + SKI-II) combination . To determine the contribution of peroxynitrite to the ER stress and viability loss induced by the lower doses of the (TMZ + SKI-II) combination, we co-treated NCH82 cells with 100 µM FeTPPS. Mitigation of peroxynitrite significantly decreased BiP protein levels in (TMZ + SKI-II)-treated cells (Fig. 6 C-Ci; Additional file 10) as well as the number of dead cells induced by the treatment, under both normoxia and hypoxia (Fig. 6 D). Altogether, these results suggest that oxidative stress has an important role in ER stress and cell death induced by the combination.
SKI-II alone and in combination reduces the levels of ceramide and its metabolites
We next assessed the effects of the treatments on the sphingolipid metabolism by measuring the concentration of relevant metabolites in NCH82 cells. Treatment for 24 h (not shown) and 48 h (Fig. 7) with SKI-II alone and in combination led to a significant decrease in ceramide (d18:1/16:0) in both oxygen conditions; the same trend was observed for ceramide (d18:1/24:1) with a significant decrease under hypoxia. The level of the ceramide precursor dihydrosphingosine (dihydroceramide was not detectable) was not significantly affected by the treatments. The ceramide/dihydrosphingosine ratio was reduced by SKI-II in both oxygen conditions. It was not affected by the combination under normoxia, but was reduced under hypoxia. Sphingomyelin, ceramide-1-phosphate (C1P) and, to a lesser extent, sphingosine, which are all direct metabolites of ceramide, showed a profile similar to that of ceramide under normoxia and hypoxia. Similar to sphingosine, S1P levels were decreased by SKI-II alone and in combination at 21% O2 (although not significantly), and remained unchanged at 3% O2. TMZ did not induce major alterations in the sphingolipid levels. Furthermore, none of the treatments significantly affected the basal level of expression of the genes coding for SK1, SK2, and DES1 in cells treated under either oxygen concentration (Fig. 8). Altogether, these data suggest that changes in sphingolipid levels after treatment with SKI-II alone and in combination were primarily the consequences of changes in ceramide levels.
The combination impairs the self-renewal of resistant glioblastoma stem cells
We performed extreme limiting dilution analysis (ELDA) to determine the effect of the (TMZ + SKI-II) combination on the self-renewal capacity of glioblastoma stem cells (GSCs) under hypoxia and normoxia. Limiting dilution assays allow the quantification of the frequency of biologically active cells in a population. Two GSC lines were used: the control line DMSO-1080 and the line TMZ-1080, selected for its resistance to 100 µM TMZ. Whereas single treatments and the combination decreased the stem cell frequency of the control DMSO-1080 cells under normoxia (Fig. 9), only SKI-II and the combination affected this frequency under hypoxia (Additional file 11). The self-renewal capacity of TMZ-1080 cells was not affected by TMZ, as expected. Treatment with SKI-II and the combination on the contrary impaired this capacity, particularly under hypoxia. The stem cell frequency decreased almost fourfold in comparison with SKI-II alone, yet it remained similar under normoxia (Fig. 9; Additional file 11).
We next examined the cytotoxic potential of the single and combined treatments towards the DMSO-1080 and TMZ-1080 cells. As shown in Fig. 10, a significant cell death induction could not be detected in either GSC line after 5 days of TMZ treatment. A slight though not significant increase in cell death was detectable after treatment with SKI-II alone under normoxia, but not under hypoxia. The combinatorial treatment led to a significant increase in cell death of DMSO-1080 cells at 21% O2, but not in any of the other conditions. We did not detect significant differences in the distribution of early apoptotic (AnnV + /PI-), late apoptotic (AnnV + /PI +) or necrotic (AnnV-/PI +) cells between the conditions. Immunofluorescence analysis of the GSC lines indicated the presence of cells dying via apoptosis as evidenced by morphological features and activation of caspase 3 (Additional file 12). This however represented a low percentage of cells. Altogether these data indicate that SKI-II alone and in combination efficiently reduce the self-renewal capacity of the two GSC lines without exhibiting significant cytotoxic effects.
SKI-II alone or in combination impairs invasion of mesenchymal MGMT-positive GB cells
Next, we analyzed the effect of the combination on the invasion capacity of the human GB U3054 cells under normoxia and hypoxia, using a 3D-spheroid model. Cell arrangement as spheroids allows to better mimic in vivo conditions, such as formation of gradients of nutrients, oxygen and catabolites, that cannot be recapitulated in monolayer cultures. The U3054 line belongs to the transcriptional mesenchymal subtype and is derived from a recurrent GB . Mesenchymal GB is considered to have an increased invasive potential and patients have a poorer prognosis [42, 43]. U3054 cells are MGMT-positive (Additional file 13), suggesting that these cells are resistant to TMZ therapy . Indeed, as shown in Fig. 11, TMZ treatment did not affect the invasion capacity of U3054 spheroids under hypoxia and normoxia. In contrast, treatment with SKI-II alone and in combination reduced invasion by about 30 to 40% after 5 days under hypoxia and normoxia, respectively. Hypoxia per se reduced the invasion capacity of the U3054 cells in the collagen matrix.
Two notorious road-blocks to effective glioblastoma therapy are triggered by low oxygen condition, a hallmark of tumors: chemotherapy resistance and expansion of glioblastoma stem-like cells. In this in vitro study, we demonstrate that combining a clinically relevant dose of TMZ (48 µM) with SKI-II, an inhibitor of three enzymes of the sphingolipid metabolism, can alleviate these blocks.
The hypoxia-induced resistance to TMZ in GB cells was effectively overcome by the combination as demonstrated by the combination index analysis. The latter shows that treatment with 48 µMTMZ combined with 2.66 µM SKI-II (D2 combination, Table 1) under hypoxia allowed to drastically reduce the dose of TMZ that would have been necessary to achieve the same high efficiency if TMZ would have been used as single treatment. TMZ in combination with SKI-II was more effective than TMZ single treatment in inhibiting cell growth, inducing cell death, reducing invasion of differentiated GB cells (NCH82 cells) and impairing the self-renewal of GSCs (DMSO-1080 and TMZ-1080 cells) under hypoxia. SKI-II alone was as effective as in combination in impairing cell invasion of GB cells (U3054 cells) and GSC self-renewal, but not in inhibiting cell growth (as shown by synergism in combination) and inducing cell death (as shown by enhancement in combination). This demonstrates that in combination we leverage the cytotoxic and cytostatic capabilities of TMZ.
Studies have shown that the serum concentration of TMZ is in the range of 20–70 µM , whereas the intratumoral TMZ concentration is in the range of 3–35 µM . In the present in vitro study, we could directly compare the cytotoxic effects of various combinations of doses of TMZ and SKI-II on the NCH82 cell line and we made the following observations. The combination used throughout our study (referred to as D2 in Table 1) was slightly more efficient than the D3 combination (96 µM TMZ and 2.66 µM SKI-II), indicating no advantage of using a higher TMZ dose in NCH82 cells. Combinations with even a very low TMZ dose were able to cause major cytotoxic effects in GB cells: the D1 combination (24 µM TMZ and 2.66 µM SKI-II) induced similar effects as the D2 combination (48 µM TMZ and 2.66 µM SKI-II) and with a comparable degree of synergy. NCH82 cells have low MGMT activity (11.5 ± 0.3 (S.D) fmol/mg) compared to other GB cell lines such as the LN18 (249 ± 10.6 (S.D) fmol/mg; our own unpublished data) and are wild-type for p53 . Thus, lower doses of TMZ are expected to be more effective in these cells . We showed that resistance to TMZ conferred by MGMT expression was also alleviated by the use of the combination, as demonstrated by the reduced invasion capacity of serum-free and MGMT-positive U3054 cells. This reduced capacity could indicate a decrease in motility, or reflect reduced proliferation and/or increased cell death. Altogether the data suggest that using TMZ in combination at half of the dose similar to that used in the clinic (100 µM) might result in an increased therapeutic efficiency in GB cells with (high) MGMT activity.
As already mentioned, a striking feature of the combination is its capacity to impair self-renewal of the two GSC lines we tested and particularly that of TMZ-resistant GSCs under hypoxia. The potential of a cell for self-renewal is taken as a measure of its stemness. Impairment of this biological state can lead to a reduced resistance of cancer stem cells to therapy. The combination however did not exhibit, or only weakly, cytotoxic effects to the DMSO-1080 and TMZ-1080 GSC cell lines. Notwithstanding this observation, we cannot rule out the possibility that the combination, by altering sphingolipid metabolism, increases their sensitivity to chemotherapeutics other than TMZ as a result of the reduction of their stemness status. This certainly would be worth to investigate.
In this study we used the SKI-II inhibitor to alter sphingolipid metabolism. As already mentioned, SKI-II was first discovered as a specific dual inhibitor of SK1 and SK2 with a Ki of 16 µM for SK1 [14, 16]. A few years later, Cingolani et al.  provided evidence that SKI-II inhibits the desaturase 1, a very important enzyme of the de novo sphingolipid pathway, with a Ki (0.3 µM). The use of 10 µM SKI-II in our previous study led to a significant increase in dihydrosphingolipids levels and a strong decrease in S1P production, in line with the capacity of the inhibitor to affect, at that concentration, the enzymatic activities of both SK1 and DES1 . In the present study, treatment with 2.66 µM SKI-II did not cause these alterations. Our data rather suggest that SKI-II at this low concentration mostly affects the desaturase as indicated by the reduced ratio of ceramide/dihydrosphingosine in both oxygen conditions, and leads to ceramide decrease. It only mildly, if at all, affects the sphingosine kinase, leading to modest decrease in S1P. The comparable profile of the levels in ceramide and its metabolites after treatments suggests that the metabolic flux between these sphingolipids is not affected and that the S1P decrease can be interpreted as the result of the lower amount of ceramide available for sphingosine generation. On the basis of these observations, we conclude that SKI-II at 2,66 µM, alone and in combination, mainly affects the desaturase and not the sphingosine kinase. The partial depletion of S1P did not hinder the anti-tumor effects of the low-dose combination and might even be beneficial for cells within the parenchyma, including astrocytes, neurons and tumor-associated cells. Extracellular S1P release by dying cells was reported to serve as a signal for macrophages to phagocytose these dying cells . Suppression of sphingosine kinase in activated microglia was reported to decrease their pro-inflammatory activities , which are desirable to eliminate tumor cells. Thus, strategies that would preserve a low level of SK activity, such as the use of low SKI-II doses, might contribute to an environment that is not only protective to healthy cells, but also supportive to immune cells through the facilitation of their phagocytic and pro-inflammatory response. This suggests that the combination might be applied together with immunotherapies designed to promote the anti-tumor activities of tumor-associated microglia/macrophages.
Oxidative stress seems to have an essential role in cell death induced by the low-dose (present study) and the high-dose (our previous study ) combination of TMZ and SKI-II, as shown by their peroxynitrite-dependent effects (e.g. ER stress and viability). However, whereas the NCH82 cells treated with the high-dose combination increased their autophagic flux and died by caspase-dependent apoptosis , cells treated with the low-dose combination showed no alteration in the autophagic flux and polarization of mitochondria and died through a caspase-independent, potentially non-apoptotic, form of cell death. The lack of accumulation of dihydrosphingolipids in the low-dose combination could correlate with the absence of changes in the autophagic flux. Several studies report autophagy induction by dihydrosphingosine and dihydroceramide . Apoptosis and autophagy were shown to be induced as a linear function of TMZ dose when MGMT is lacking . Therefore, the difference in the dose of TMZ used in our studies (48 µM vs 500 µM) have likely contributed to the dissimilarities observed in autophagy and cell death mechanisms. The lack of PI-positive/Annexin V-negative NCH82 cells suggests that membrane integrity is kept and primary necrosis can be ruled out. The increase in the percentage of Annexin V-positive cells over time in all conditions (except in the vehicle control) demonstrates that phosphatidylserine (PS) exposure is a morphological feature induced by SKI-II and TMZ. Exposure of PS is mediated via scramblases, which are activated by caspase 3/7 cleavage. However, caspase-independent mechanisms of PS exposure were also described in non-apoptotic forms of cell death . Namely, PS exposure was shown to occur due to sustained increase in the concentration of cytoplasmic calcium in a form of paraptosis . The appearance of extensive cytoplasmic vacuolization, ER dilation, and lack of caspase activation and mitochondrial depolarization suggests that the low-dose (TMZ + SKI-II) combination could induce a paraptosis-like cell death in NCH82 cells . Although plausible, this conclusion is a hypothesis that will have to be proven.
The mechanism underlying oxidative stress induction by SKI-II in our experimental system is unclear. Sphingolipids are able to regulate cellular redox homeostasis: ceramide can activate NADPH oxidase, supporting peroxynitrite formation ; whereas, glucosylceramide can inhibit it . We observed a decrease of C16- and C24-ceramide, ceramide-1-phosphate and sphingomyelin in both oxygen conditions. A decreased level of C16-ceramide has been related to ER stress induction as it can promote Ca2 + release from the ER and activate the ATF6/CHOP axis . S1P, which is still produced under treatment with the combination, might contribute to oxidative stress through interaction with its ligand, S1PR1. Indeed, activation of S1PR1 was reported to support peroxynitrite formation by increasing the generation of NADPH oxidase-derived superoxide and nitric oxide synthase-derived nitric oxide . The role of S1P and its receptor in the response of the tumor cells to the combination remains to be elucidated.
Experience of the past decades in delineating cancer therapeutic strategies illustrates that focusing on one singular target is insufficient and ineffective, namely in tumors as complex as glioblastoma. In the present study, we wanted to explore the benefit of combining two drugs affecting two different biological processes in order to demonstrate at a preclinical level the value of the sphingolipid metabolism as a target for cancer therapy, alone and in combination with temozolomide, the first-line agent in glioblastoma therapy. Given their key role in controlling cell fate, sphingosine kinases are prime targets to modulate sphingolipid metabolism  and have been the focus of many in vitro and in vivo studies in the past twenty years [58, 59]. However, SK inhibitors have scarcely been implemented in clinical trials and, to paraphrase Gault et al. , they are “still benched on (their) way to the bedside”. One explanation may lie in the contradictory data that have been reported on the cytotoxic activities of SK inhibitors on various tumor cell models, thus diminishing their potential value. We hope, with the present data, to rekindle interest for SK inhibitors, and, beyond SK inhibitors, for the sphingolipid metabolism as a relevant and attractive therapeutic target. Indeed, we provide evidence that modulation of sphingolipid metabolism, combined with the use of TMZ, has the potential to affect various biological functions of differentiated and stem glioblastoma cells under atmospheric and low oxygen concentration. These observations warrant further experimental exploration to assess the clinical potential of such a combination to reduce the glioblastoma stem cell population in tumors and prevent relapse.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
Dihydroceramide desaturase 1
Median effect dose
Extreme limiting dilution analysis
Glioblastoma stem cell
Reactive oxygen species
S1P receptor 1
Unfolded protein response
Stupp R, Mason WP, van den Bent MJ. Radiotherapy plus Concomitant and Adjuvant Temozolomide for Glioblastoma. Oncol Times. 2005;27:15–6.
Denny BJ, Tsang LLH, Slack JA, Wheelhouse RT, Stevens MFG. NMR and molecular modeling investigation of the mechanism of activation of the antitumor drug temozolomide and its interaction with DNA. Biochemistry. 1994;33:9045–51.
Brennand J, Margison GP. Reduction of the toxicity and mutagenicity of alkylating agents in mammalian cells harboring the Escherichia coli alkyltransferase gene. Proc Natl Acad Sci U S A. 1986;83:6292–6.
Stupp R, Hegi ME, Mason WP, van den Bent MJ, Taphoorn MJ, Janzer RC, et al. Effects of radiotherapy with concomitant and adjuvant temozolomide versus radiotherapy alone on survival in glioblastoma in a randomised phase III study: 5-year analysis of the EORTC-NCIC trial. Lancet Oncol. 2009;10:459–66.
Weller M, Cloughesy T, Perry JR, Wick W. Standards of care for treatment of recurrent glioblastoma-are we there yet? Neuro Oncol. 2013;15:4–27.
Brat DJ, Castellano-Sanchez AA, Hunter SB, Pecot M, Cohen C, Hammond EH, et al. Pseudopalisades in glioblastoma are hypoxic, express extracellular matrix proteases, and are formed by an actively migrating cell population. Cancer Res. 2004;64:920–7.
Auffinger B, Tobias AL, Han Y, Lee G, Guo D, Dey M, et al. Conversion of differentiated cancer cells into cancer stem-like cells in a glioblastoma model after primary chemotherapy. Cell Death Differ. 2014;21:1119–31.
Wang P, Wan WW, Xiong SL, Feng H, Wu N. Cancer stem-like cells can be induced through dedifferentiation under hypoxic conditions in glioma, hepatoma and lung cancer. Cell Death Discov. 2017;3:16105.
Marzagalli M, Fontana F, Raimondi M, Limonta P. Cancer stem cells—key players in tumor relapse. Cancers (Basel). 2021;13:1–23.
Hannun YA, Obeid LM. Sphingolipids and their metabolism in physiology and disease. Nat Rev Mol Cell Biol. 2018;19:175–91.
Zhang H, Desai NN, Olivera A, Seki T, Brooker G, Spiegel S. Sphingosine-1-phosphate, a novel lipid, involved in cellular proliferation. J Cell Biol. 1991;114:155–67.
Obeid LM, Linardic CM, Karolak LA, Hannun YA. Programmed cell death induced by ceramide. Science. 1993;259:1769–71.
Abuhusain HJ, Matin A, Qiao Q, Shen H, Kain N, Day BW, et al. A metabolic shift favoring sphingosine 1-phosphate at the expense of ceramide controls glioblastoma angiogenesis. J Biol Chem. 2013;288:37355–64.
Bose R, Verheij M, Haimovitz-Friedman A, Scotto K, Fuks Z, Kolesnick R. Ceramide synthase mediates daunorubicin-induced apoptosis: an alternative mechanism for generating death signals. Cell. 1995;82:405–14.
French KJ, Upson JJ, Keller SN, Zhuang Y, Yun JK, Smith CD. Antitumor activity of sphingosine kinase inhibitors. J Pharmacol Exp Ther. 2006;318:596–603.
Cingolani F, Casasampere M, Sanllehí P, Casas J, Bujons J, Fabrias G. Inhibition of dihydroceramide desaturase activity by the sphingosine kinase inhibitor SKI II. J Lipid Res. 2014;55:1711–20.
Gao P, Peterson YK, Smith RA, Smith CD. Characterization of isoenzyme-selective inhibitors of human sphingosine kinases. PLoS One. 2012;7(9):e44543. https://doi.org/10.1371/journal.pone.0044543.
Noack J, Choi J, Richter K, Kopp-Schneider A, Régnier-Vigouroux A. A sphingosine kinase inhibitor combined with temozolomide induces glioblastoma cell death through accumulation of dihydrosphingosine and dihydroceramide, endoplasmic reticulum stress and autophagy. Cell Death Dis. 2014;5:e1425.
Casasampere M, Ordóñez YF, Casas J, Fabrias G. Dihydroceramide desaturase inhibitors induce autophagy via dihydroceramide-dependent and independent mechanisms. Biochim Biophys Acta - Gen Subj. 2017;1861:264–75.
Gupta P, Taiyab A, Hussain A, Alajmi MF, Islam A, Hassan MI. Targeting the sphingosine kinase/sphingosine-1-phosphate signaling axis in drug discovery for cancer therapy. Cancers (Basel). 2021;13(8):1898. https://doi.org/10.3390/cancers13081898.
Ostermann S, Csajkag C, Buclin T, Leyvraz S, Lejeune F, Decosterd LA, et al. Plasma and cerebrospinal fluid population pharmacokinetics of temozolomide in malignant glioma patients. Clin Cancer Res. 2004;10:3728–36.
Rosso L, Brock CS, Gallo JM, Saleem A, Price PM, Turkheimer FE, et al. A new model for prediction of drug distribution in tumor and normal tissues: Pharmacokinetics of temozolomide in glioma patients. Cancer Res. 2009;69:120–7.
Karcher S, Steiner H-H, Ahmadi R, Zoubaa S, Vasvari G, Bauer H, et al. Different angiogenic phenotypes in primary and secondary glioblastomas. Int J Cancer. 2006;118:2182–9.
Barrantes-Freer A, Kim E, Bielanska J, Giese A, Mortensen LS, Schulz-Schaeffer WJ, et al. Human glioma-initiating cells show a distinct immature phenotype resembling but not identical to NG2 glia. J Neuropathol Exp Neurol. 2013;72:307–24.
Xie Y, Bergström T, Jiang Y, Johansson P, Marinescu VD, Lindberg N, et al. The Human Glioblastoma Cell Culture Resource: Validated Cell Models Representing All Molecular Subtypes. EBioMedicine. 2015;2:1351–63.
Vichai V, Kirtikara K. Sulforhodamine B colorimetric assay for cytotoxicity screening. Nat Protoc. 2006;1:1112–6.
Chou TC, Talalay P. Quantitative analysis of dose-effect relationships: the combined effects of multiple drugs or enzyme inhibitors. Adv Enzyme Regul. 1984;22 C:27–55.
Chou T-C. Drug Combination Studies and Their Synergy Quantification Using the Chou-Talalay Method. 2010. https://doi.org/10.1158/0008-5472.CAN-09-1947.
Chou TC. Preclinical versus clinical drug combination studies. Leuk Lymphoma. 2008;49:2059–80.
Chou T-C, Martin N. CompuSyn for Drug Combinations and for General Dose-Effect Analysis User’s Guide A Computer Program for Quantitation of Synergism and Antagonism in Drug Combinations, and the Determination of IC 50, ED 50, and LD 50 Values. Paramus, NJ: ComboSyn Inc.; 2005.
Hu Y, Smyth GK. ELDA: Extreme limiting dilution analysis for comparing depleted and enriched populations in stem cell and other assays. J Immunol Methods. 2009;347:70–8.
Cisneros Castillo LR, Oancea AD, Stüllein C, Régnier-Vigouroux A. A novel computer-assisted approach to evaluate multicellular tumor spheroid invasion assay. Sci Rep. 2016;6:35099.
Post JM, Lerner R, Schwitter C, Lutz B, Lomazzo E, Bindila L. Lipidomics and Transcriptomics in Neurological Diseases. J Vis Exp. 2022;(181). https://doi.org/10.3791/59423.
Lerner R, Cuadrado DP, Post JM, Lutz B, Bindila L. Broad lipidomic and transcriptional changes of prophylactic pea administration in adult mice. Front Neurosci. 2019;13:1–16.
McIlwain DR, Berger T, Mak TW. Caspase Functions in cell death and disease. Cold Spring Harb Perspect Biol. 2013;5(4):a008656. https://doi.org/10.1101/cshperspect.a008656.
Christensen ME, Jansen ES, Sanchez W, Waterhouse NJ. Flow cytometry based assays for the measurement of apoptosis-associated mitochondrial membrane depolarisation and cytochrome c release. Methods. 2013;61:138–45.
Würstle S, Schneider F, Ringel F, Gempt J, Lämmer F, Delbridge C, et al. Temozolomide induces autophagy in primary and established glioblastoma cells in an EGFR independent manner. Oncol Lett. 2017;14:322–8.
Klionsky DJ, Abdelmohsen K, Abe A, Abedin MJ, Abeliovich H, Acevedo Arozena A, Adachi H, Adams CM, Adams PD, Adeli K, Adhihetty PJ, Adler SG, Agam G, Agarwal R, Aghi MK, Agnello M, Agostinis P, Aguilar PV, Aguirre-Ghiso J, Airoldi EM, Ait-Si-Ali S, Akemat ZS. Guidelines for use and interpretation of assays for monitoring autophagy (3rd edition). Autophagy. 2016;12:1–222.
Arismendi-Morillo G. Electron microscopy morphology of the mitochondrial network in gliomas and their vascular microenvironment. Biochim Biophys Acta - Bioenerg. 2011;1807:602–8.
Santos CXC, Tanaka LY, Wosniak J, Laurindo FRM. Mechanisms and implications of reactive oxygen species generation during the unfolded protein response: roles of endoplasmic reticulum oxidoreductases, mitochondrial electron transport, and NADPH oxidase. Antioxid Redox Signal. 2009;11:2409–27.
Dickhout JG, Hossain GS, Pozza LM, Zhou J, Lhoták Š, Austin RC. Peroxynitrite causes endoplasmic reticulum stress and apoptosis in human vascular endothelium: implications in atherogenesis. Arterioscler Thromb Vasc Biol. 2005;25:2623–9.
Wang Q, Hu B, Hu X, Kim H, Squatrito M, Scarpace L, et al. Tumor evolution of glioma-intrinsic gene expression subtypes associates with immunological changes in the microenvironment. Cancer Cell. 2017;32:42-56.e6.
Verhaak RGW, Hoadley KA, Purdom E, Wang V, Qi Y, Wilkerson MD, et al. Integrated genomic analysis identifies clinically relevant subtypes of glioblastoma characterized by abnormalities in PDGFRA, IDH1, EGFR, and NF1. Cancer Cell. 2010;17:98–110.
van Nifterik KA, van den Berg J, van der Meide WF, Ameziane N, Wedekind LE, Steenbergen RDM, et al. Absence of the MGMT protein as well as methylation of the MGMT promoter predict the sensitivity for temozolomide. Br J Cancer. 2010;103:29–35.
Di Piazza M, Mader C, Geletneky K, Herrero y Calle M, Weber E, Schlehofer J, et al. Cytosolic Activation of Cathepsins Mediates Parvovirus H-1-Induced Killing of Cisplatin and TRAIL-Resistant Glioma Cells. J Virol. 2007;81:4186–98.
Hermisson M, Klumpp A, Wick W, Wischhusen J, Nagel G, Roos W, et al. O6-methylguanine DNA methyltransferase and p53 status predict temozolomide sensitivity in human malignant glioma cells. J Neurochem. 2006;96:766–76.
Gude DR, Alvarez SE, Paugh SW, Mitra P, Yu J, Griffiths R, et al. Apoptosis induces expression of sphingosine kinase 1 to release sphingosine-1-phosphate as a “come-and-get-me” signal. FASEB J. 2008;22:2629–38.
Wang P, Wan WW, Xiong SL, Feng H, Wu N. Cancer stem-like cells can be induced through dedifferentiation under hypoxic conditions in glioma, hepatoma and lung cancer. Cell Death Discov. 2017;3:132–44.
Siddique MM, Li Y, Chaurasia B, Kaddai VA, Summers SA. Dihydroceramides: From bit players to lead actors. J Biol Chem. 2015;290:15371–9.
He Y, Kaina B. Are there thresholds in glioblastoma cell death responses triggered by temozolomide? Int J Mol Sci. 2019;20(7):1562. https://doi.org/10.3390/ijms20071562.
Ferraro-Peyret C, Quemeneur L, Flacher M, Revillard J-P, Genestier L. Caspase-independent phosphatidylserine exposure during apoptosis of primary t lymphocytes. J Immunol. 2002;169:4805–10.
Jambrina E, Alonso R, Alcalde M, del Rodríguez MC, Serrano A, Martínez AC, et al. Calcium influx through receptor-operated channel induces mitochondria-triggered paraptotic cell death. J Biol Chem. 2003;278:14134–45.
Sperandio S, De Belle I, Bredesen DE. An alternative, nonapoptotic form of programmed cell death. Proc Natl Acad Sci U S A. 2000;97:14376–81.
Zhang DX, Ping Zou A, lan Li P. Ceramide-induced activation of NADPH oxidase and endothelial dysfunction in small coronary arteries. Am J Physiol - Hear Circ Physiol. 2003;284(2):53–2.
Barth BM, Gustafson SJ, Young MM, Fox TE, Shanmugavelandy SS, Kaiser JM, et al. Inhibition of NADPH oxidase by glucosylceramide confers chemoresistance. Cancer Biol Ther. 2010;10:1126–36.
Senkal CE, Ponnusamy S, Bielawski J, Hannun YA, Ogretmen B. Antiapoptotic roles of ceramide-synthase-6-generated C 16 -ceramide via selective regulation of the ATF6/ CHOP arm of ER-stress-response pathways. FASEB J. 2010;24:296–308.
Doyle T, Finley A, Chen Z, Salvemini D. Role for peroxynitrite in sphingosine-1-phosphate-induced hyperalgesia in rats. Pain. 2011;152:643–8.
Gault CR, Obeid LM. Still benched on its way to the bedside: Sphingosine kinase 1 as an emerging target in cancer chemotherapy. Crit Rev Biochem Mol Biol. 2011;46:342–51.
Cao M, Ji C, Zhou Y, Huang W, Ni W, Tong X, et al. Sphingosine kinase inhibitors: a patent review. Int J Mol Med. 2018;41:2450–60.
We acknowledge the help of Bernd Kaina in providing the LN229 and LN229-MGMT cell lines and Lea Beltzig in performing the MGMT analysis.
Open Access funding enabled and organized by Projekt DEAL. No funding to declare.
Ethics approval and consent to participate
In the study, there was NO direct involvement of humans; tissues have not been removed from human patients for the current study purpose.
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.
LC and MS conditions for sphingolipid analysis by liquid chromatography/multiple reaction monitoring (LC/MRM). LC-MS gradient and solvent system are depicted in the table below. The MS ion source parameters were as follows: curtain gas was 40 psi; source heater temperature 550 °C; ion spray voltage was set to−4500 V in negative ion mode and to +5200 V in positive ion mode.
Sphingolipid analysis by liquid chromatography/multiple reaction monitoring (LC/MRM). MRM transitions and parameters used in positive ion mode for sphingolipid profile and quantification are depicted.
Calculation of the median-effect dose. The median-effect dose (ED50) of temozolomide (TMZ) and the sphingosine kinase inhibitor SKI-II was calculated in NCH82 cells at normoxia (21% O2) and hypoxia (3% O2) for 5 days using the median-effect analysis program CompuSyn (Nick Martin, MIT, Cambridge, MA, 2005).
Combination design of temozolomide (TMZ) and the sphingosine kinase inhibitor SKI-II. Fifteen (TMZ + SKI-II) combinations were made based on the ED50 of each drug (SKI-II ED50 = 1.33 µM; TMZ ED50 = 96 µM).
Dose Reduction Index (DRI) of combinations of temozolomide (TMZ) and the sphingosine kinase inhibitor SKI-II at 21% and 3% O2 in NCH82 cells. The dose (µM) of TMZ and SKI-II and respective DRI for each combination was calculated using CompuSyn software. DRI=1, >1, and <1 indicates no dose-reduction, favorable dose-reduction, and not favorable dose-reduction, respectively, for each drug in the combination.
Full-length blots of cleaved caspase-3 and GAPDH detection shown in Fig. 3 A. After the protein transfer, the nitrocellulose membrane was cut below the 30 kDa marker, and each strip was incubated with antibodies against GAPDH or Cleaved Caspase-3. Left panel (21 % O2): exposure times were 2 min for GAPDH, 5 min for Cleaved Caspase-3. Right panel (3% O2): exposure times were 2 min for GAPDH, 5 min for Cleaved Caspase-3.
The combination does not affect autophagic flux in glioblastoma cells. NCH82 cells were treated with 48 µM temozolomide (TMZ) combined with 2.66 µM SKI-II (TMZ+SKI-II) in the presence and absence of 100 nM Bafilomycin A1 (BA1) under 21% O2 and 3% O2. (A, B) Cell extracts were separated on SDS- PAGE and transferred on nitrocellulose membrane. DMSO and (TMZ+SKI-II) gels were run in parallel. Representative blots are shown. The expression levels of p62 (Ai, Bi) and LC3-II (Aii, Bii) were quantified and normalized to GAPDH. Autophagic flux under basal conditions (DMSO) and under treatment (TMZ + SKI-II) was determined by the subtraction of LC3-II levels without BA1 (Ctrl) from LC3-II levels with BA1 (Aiii, Biii). LC3-II (n = 3); p62 (n = 2). The results are shown as mean (±SD). Two-way ANOVA followed by Tukey’s multiple comparisons test was performed for statistical analysis. Full-length blots are presented in Additional file 8.
Full-length blots of LC3, p62 and GAPDH detection shown in Additional file 7 A (panel A) and Additional file 7 B (panel B). After the protein transfer, the nitrocellulose membrane was cut in three pieces, below the 50 and 30 kDa markers and each strip was incubated with antibodies against p62, GAPDH or LC3B. Panel A, left (DMSO) and right (TMZ + SK-II): exposure times were p62, 3 min (left and right); GAPDH, 10 sec (left) and 30 sec (right); LC3B, 4 min (left and right). Panel B, left (DMSO) and right (TMZ + SK-II): exposure times were p62, 4 min (left) and 3 min (right); GAPDH, 1 min (left) and 30 sec (right); LC3B, 4 min (left and right).
Full-length blots of BiP/GRP78 and GAPDH detection shown in Fig. 6 A. Left (21 % O2) and right (3 % O2) panels: After the protein transfer, the nitrocellulose membrane was cut in two pieces, below the 50 kDa marker, and each strip was incubated with antibodies against GAPDH or BiP/GRP78. Exposure times were: GAPDH 1 min (left and right); BiP/GRP78 1 min (left and right).
Full-length blots of BiP/GRP78 and GAPDH detection shown in Fig.6 C. After the protein transfer, the nitrocellulose membrane was cut in two pieces, between the 60 kDa and the 50 kDa markers and each strip was incubated with antibodies against GAPDH or BiP/GRP78. Exposure times were: GAPDH 20 sec; BiP/GRP78 1 min.
ELDA analysis: Full data files.
Immunofluorescence detection of activated caspase-3 in treated GSC. DMSO-1080 and TMZ-1080 GSC cells plated on ornithin-coated glass coverslips were treated with 2.66 µM SKI-II, 48 µM temozolomide (TMZ) and the combination (TMZ+SKI-II) in triplicates. After 5 days of incubation, cells were fixed with 4% paraformaldehyde and stained for cleaved caspase-3, nestin and Dapi. (A) Quantification of cleaved caspase-3 positive cells. Data are presented as percent of total cells counted in 7 to 9 regions of triplicates, and for a total number of cells of 700 to 1100. (B) Representative images of apoptotic cells in 1080-TMZ cells treated with SKI-II under normoxia (left) or hypoxia (right). Magnification 40X, scale bar 50 µM. Left panels in each composite picture: enlarged images of apoptotic cells (marked by arrows in the right panels) manifesting morphological hallmarks of apoptosis such as chromatin condensation and cytoplasmic disintegration. n = 1.
The U3054 cells express the O6-methylguanine-DNA methyltransferase.
About this article
Cite this article
Sousa, N., Geiß, C., Bindila, L. et al. Targeting sphingolipid metabolism with the sphingosine kinase inhibitor SKI-II overcomes hypoxia-induced chemotherapy resistance in glioblastoma cells: effects on cell death, self-renewal, and invasion. BMC Cancer 23, 762 (2023). https://doi.org/10.1186/s12885-023-11271-w
- Glioblastoma therapy
- Drug combination
- Glioblastoma stem cells
- Cell death