Abstract
Multi-agent chemotherapeutic regimes remain the cornerstone treatment for Ewing sarcoma, the second most common bone malignancy diagnosed in pediatric and young adolescent populations. We have reached a therapeutic ceiling with conventional cytotoxic agents, highlighting the need to adopt novel approaches that specifically target the drivers of Ewing sarcoma oncogenesis. As KDM1A/lysine-specific demethylase 1 (LSD1) is highly expressed in Ewing sarcoma cell lines and tumors, with elevated expression levels associated with worse overall survival (P = 0.033), this study has examined biomarkers of sensitivity and mechanisms of cytotoxicity to targeted KDM1A inhibition using SP-2509 (reversible KDM1A inhibitor). We report, that innate resistance to SP-2509 was not observed in our Ewing sarcoma cell line cohort (n = 17; IC50 range, 81 –1,593 nmol/L), in contrast resistance to the next-generation KDM1A irreversible inhibitor GSK-LSD1 was observed across multiple cell lines (IC50 > 300 μmol/L). Although TP53/STAG2/CDKN2A status and basal KDM1A mRNA and protein levels did not correlate with SP-2509 response, induction of KDM1B following SP-2509 treatment was strongly associated with SP-2509 hypersensitivity. We show that the transcriptional profile driven by SP-2509 strongly mirrors KDM1A genetic depletion. Mechanistically, RNA-seq analysis revealed that SP-2509 imparts robust apoptosis through engagement of the endoplasmic reticulum stress pathway. In addition, ETS1/HIST1H2BM were specifically induced/repressed, respectively following SP-2509 treatment only in our hypersensitive cell lines. Together, our findings provide key insights into the mechanisms of SP-2509 cytotoxicity as well as biomarkers that can be used to predict KDM1A inhibitor sensitivity in Ewing sarcoma. Mol Cancer Ther; 17(9); 1902–16. ©2018 AACR.
Introduction
Ewing sarcoma (ES) is a highly aggressive pediatric bone tumor that has only seen incremental improvements in survival rates for patients diagnosed with primary metastatic or relapsed disease in recent years (<20% event-free survival). At the genome level, the sole distinguishing feature of ES is the presence of the reciprocal EWS/ETS translocation, with the resulting fusion protein, EWS/FLI in 85% of cases, acting as the molecular driver for tumor development (1, 2). Although this poorly differentiated tumor has a strict dependence on EWS/ETS fusion proteins for oncogenesis, therapeutic targeting of EWS/FLI remains unrealized due to its intrinsically disordered protein nature and lack of enzymatic activity (3). Current chemotherapy backbones for the treatment of primary ES have historically consisted of regimens based on alkylating and anthracycline agents. However, extensive protocol-driven clinical research evaluating dose intensification and schedule optimization of the current five-drug regimen (vincristine, adriamycin, cyclophosphamide, ifosfamide, and etoposide) suggests that further modification will unlikely produce additional benefit. This coupled with the high risk of anthracycline-induced cardiomyopathy and development of secondary neoplasms (4) clearly underscores the urgent requirement to incorporate new therapeutic agents with lower toxicity profiles to improve overall survival for these patients.
For the past several decades, the overwhelming emphasis in cancer biology has been defining and targeting DNA aberrations. However, there is growing appreciation that dysregulation of epigenetic machinery and chromatin modifications are important mechanisms utilized by tumors to favorably modulate DNA repair, cell-cycle control, and apoptosis-promoting genes. This is particularly pertinent for ES, as recent high-throughput screening efforts have shown that this malignancy possesses one of the lowest mutation rates among all cancers (0.15 mutations/Mb; refs. 5–7), yielding a paucity of pharmacologically actionable mutations.
Histone methylation is a major determinant of chromatin structure and function and has shown to play critical roles in transcriptional regulation and genomic stability (8). KDM1A (LSD1/BHC110) was the first flavine adenine dinucleotide (FAD)-dependent lysine-specific demethylase identified to regulate chromatin states through the removal of mono- and dimethyl groups (H3K4 or H3K9; ref. 9). KDM1A is overexpressed in several solid and hematological malignancies including breast, colorectal, prostate, ovarian cancer, and sarcomas with KDM1A inhibition significantly reducing cell proliferation, migration, and invasion in vitro (10). Notably, elevated KDM1A expression is associated with poor clinical outcome (11–14). These observations have led to the active development of several small-molecule KDM1A inhibitors, with tranylcypromine (NCT02273102, NCT02261779), GSK-LSD1 (NCT02177812, NCT02034123), and ORY-1001 (2013-002447-29) currently undergoing phase I/II clinical evaluation in patients with acute myeloid leukemia (AML) and small cell lung carcinoma (SCLC).
In 2012, Bennani-Baiti and colleagues documented KDM1A overexpression in ES patient tumors (15) highlighting the therapeutic utility of targeted KDM1A inhibition for this disease. We previously demonstrated that directed KDM1A inhibition in ES with the use of SP-2509 (previously HCI-2509), a noncompetitive, phenylethylidene-benzohydrazide–reversible KDM1A inhibitor (Ki: 31 nmol/L, KDM1A IC50:13 nmol/L; ref. 16), reverses the EWS/ETS transcriptional signature, impairs several EWS/ETS-associated oncogenic phenotypes, and shows single-agent efficacy in multiple xenograft models of ES (17). However, the precise mechanism by which SP-2509 induces apoptosis remains unknown. In addition, it is imperative that biomarkers of ES sensitivity and resistance to this exciting epigenetic therapeutic are elucidated to guide clinical trial patient stratification. In this report, we have built on our previous work and conducted a series of studies to define mechanisms of ES SP-2509 sensitivity as well as the biological pathways that govern apoptotic cytotoxicity.
Materials and Methods
Cell lines and compounds
Cell lines were sourced and cultured as detailed in Supplementary Table S1. Cell cultures were not used beyond 2 months from initial thawing with culture supernatants tested yearly for mycoplasma infection using a PCR based detection kit (Southern Biotech). Cell lines were authenticated by STR profiling (Genetica LabCorp). SP-2509 was provided by Dr. Sunil Sharma (TGen Clinical Sciences), doxorubicin hydrochloride and GSK-LSD1 (18) were purchased from Cayman Chemical, thapsigargin and AraC (cytosine arabinoside) were sourced from Sigma-Aldrich. [3H]GSK-LSD1 (specific activity, 9 Ci/mmol) was purchased from ViTrax Radiochemicals.
Immunodetection
Whole-cell lysate (35 μg) was run on 4% to 15% Tris-Glycine polyacrylamide gels and transferred onto nitrocellulose membranes using a iBlot2 (Thermo Fisher Scientific) according to the manufacturer's instructions. Membranes were blocked in Odyssey Blocking Buffer (PBS; LI-COR Bioscience) for 1 hour at room temperature followed by overnight 4°C incubation with primary antibodies. Immunodetection was achieved after incubation with infrared (IR)-dye–conjugated 800CW secondary antibodies (LiCor) with bands visualized using the Odyssey Imaging System. Antibodies provided in Supplementary Table S1. Densitometry analysis was performed using ImageJ software (V1.51).
IncuCyte cell proliferation and caspase-3/7 assay
ES cells were seeded (4,000–10,000 cells/well) in 96-clear-microtiter plates (triplicate wells per condition) and left to adhere overnight. If required, compounds of interest were added to each well (SP-2509, GSK-LSD1, doxorubicin) 18 hours after seeding. Real-time apoptosis assays were assessed through the addition of IncuCyte caspase-3/7 apoptosis assay reagent (Essen BioSciences), final concentration of 5 μmol/L. Phase contest and/or green fluorescent images were taken in the IncuCyte ZOOM Kinetic Imaging System (Essen BioScience) at 3 hour intervals for a minimum of 96 hours. Cell confluence (phase contrast) or green fluorescence (green object count per well) was evaluated using IncuCyte ZOOM 2016A software (Essen BioScience). Data represent mean ± SEM from a minimum of three independent experiments.
qRT-PCR
Total RNA was extracted using the RNeasy kit (Qiagen) with on-column DNase digestion. cDNA synthesis and subsequent qRT-PCR were performed with 50 ng of total RNA using iTaq Universal SYBR Green 1 Step Reaction Mix (Bio-Rad), according to the manufacturer's protocol. Primer specificity was determined using melting curve analysis with all amplicons sequenced verified during primer optimization studies. Reactions were processed on a CFX connect Real-Time System (Bio-Rad) with subsequent gene expression quantified using the ΔΔCT method from triplicate reactions. Gene expression was normalized to the internal housekeeping gene RPL19. Primer sequences are listed in Supplementary Table S1.
shRNA generation
shRNAs targeting the 3′UTR of KDM1A, KDM1B, and ERN1 were designed using sirna.wi.mit.edu (Whitehead Institute at MIT). shRNA sequences are listed in Supplementary Table S1. To create shRNA expression retroviral constructs, oligonucleotides were annealed and cloned into Age1 and EcoR1 sites of pMKO.1. To generate retroviruses, 293 EBNA cells were cotransfected (Mirus Bio TransIT-LT1) with retroviral expression plasmids (10 μg), vesicular stomatitus virus-G glycoprotein and gag/pol packaging plasmids, with filtered viral supernatant collected over a 48 hour period. Retroviral-infected ES cell lines were selected in puromycin (0.05–2 mg/mL; Sigma-Aldrich) for a minimum of 48 hours.
Soft-agar assays
Cells were seeded at a density of 7,500 to 12,500 cells per 6 cm plate in 0.8% SeaPlaque GTG agarose (Lonza), in media containing 20% FCS, Iscove's modification of Eagle's media, penicillin/streptomycin/glutamine and puromycin. Duplicate plates per condition were seeded. Colonies were quantified using ImageJ software (V1.51) a minimum of 16 days after seeding.
Viability assays
ES cells were seeded (4,000–10,000 cells/well) in 96-white-microtiter plates (triplicate wells per condition) and left to adhere overnight. Cells were treated with vehicle control, media control, or serial dilutions of SP-2509 18 hours after seeding (0.1% final DMSO concentration). Cell viability was assessed 72 hours after treatment using CellTiter-Glo (Promega) according to the manufacturer's instructions with luminescence read on a GloMax 96 Microplate Luminometer (Promega). Viability was calculated relative to vehicle control cells with IC50 values calculated using GraphPad Prism (Version 7.00).
XBP1 splice variant PCR
To amplify unspliced and spliced XBP1 variants, 50 ng of total RNA was reverse transcribed and amplified using iTaq Universal SYBR Green 1 Step Reaction Mix (Bio-Rad) and the following primers (19): Forward: 5′-TTACGAGAGAAAACTCATGGCC-'3, Reverse 5′-GGGTCCAAGTTGTCCAGAATGC-'3. PCR conditions were as follows; 95°C for an initial 1 minute followed by 50 cycles of denaturation at 95°C for 10 seconds, annealing at 57°C for 30 seconds, and extension at 72°C for 30 seconds. PCR products were electrophoresed on a 2.5% agarose gel. Size difference between unspliced XBP1 (289 bp) and spliced XBP1 (263bp) is 26 nucleotides.
Statistically analysis
P values were calculated using Student t test using Graph Pad Prism (Version 7). Additional methods are provided in Supplementary Table S1.
Results
KDM1A is highly expressed in ES tumors and cell lines
KDM1A represents an ideal therapeutic target for the treatment of ES for two reasons; KDM1A is highly expressed in ES tumors (15) and secondly no mutations in KDM1A have ever been documented in ES tumor samples (10). To confirm elevated KDM1A expression in ES compared with other malignancies, we examined the Broad Institute Cancer Cell Line panel database (1,035 cell lines, 36 distinct cancer subtypes). ES was the second highest KDM1A-expressing malignancy after T-cell AML (Fig. 1A), with expression levels significantly higher than other solid bone sarcomas including osteosarcoma (P = 0.0206) and chondrosarcoma (P = 0.0497). These findings were further supported through immunohistochemical analysis of ES tumor microarrays which demonstrated moderate–strong KDM1A expression in the majority of tumors (62.1%; Fig. 1B).
KDM1A is highly expressed in ES cell lines and tumors. A, Broad Institute Cancer Cell line expression data of KDM1A across 36 distinct cancer entities. B, Representative immunohistochemical analysis of KDM1A in ES patient tumors. C, Event-free and overall survival of ES patients stratified by high (median + 0.5 × MAD), intermediate and low (median −0.5 × MAD) KDM1A expression. Expression data obtained from Postel-Vinay et al. (24). D, Correlation between KDM1A expression and patient gender. Expression data obtained as described in D. Log-rank (Mantel–Cox test) used to determine survival significance. E, Representative images and quantification of soft-agar colonies of A673 and TTC-466 cells stably transduced with iKDM1A or iLuc control shRNA constructs following puromycin selection. F, Real-time live-cell imaging (IncuCyte ZOOM) of A673 and TTC-466 cells transduced with iKDM1A or iLuc shRNA constructs. Cell proliferation measured for 144 and 180 hours, respectively. Dashed line represents 100% confluency (iLuc). Representative phase contrast images are also depicted. Data represent mean ± SEM from three independent experiments for A673 cells and mean ± SD from two independent experiments for TTC-466 cells. Asterisks denote statistical significance (*, P < 0.05; ***, P <0.001).
We next elevated whether KDM1A expression correlates with ES clinical outcome or pathologic features. Analysis of KDM1A expression in six ES microarray/RNA-seq expression platform studies was conducted (6, 20–24) with no correlation between KDM1A expression and pathologic criteria including age, tumor location, EWS/ETS translocation, or presence of metastatic disease observed in any study (Supplementary Table S2). Although females displayed significantly higher KDM1A expression compared with males (P = 0.014) in one large cohort of 117 patients (Fig. 1D), this association was not observed in four smaller cohort studies (n = 14–46 patients; Supplementary Table S2). Notably, high KDM1A expression was significantly correlated with worse overall survival (ref. 24; P = 0.033) and a worse overall survival trend in two additional small studies (ref. 22, P = 0.054; ref. 20, P = 0.052; Fig. 1C; Supplementary Fig. S1A). In contrast, high KDM1A expression was significantly correlated with poor event-free survival only in 1/4 studies (ref. 22, P = 0.016; Supplementary Fig. S1A). Data described in Postel-Vinay and colleagues (24) is the only study that assessed KDM1A expression in chemotherapy/radiotherapy naïve primary tumor samples. The remaining 3 studies all included patient samples from heavily treated and relapsed patients which could account for the discrepancy in survival findings. Together, although the clinical ES cohorts are small, our findings suggest that high KDM1A expression levels are associated worse overall survival, and that targeting KDM1A has broad utility for this aggressive malignancy.
KDM1A is required for ES tumorigenesis
We previously demonstrated that targeted inhibition of KDM1A with SP-2509 significantly impairs the tumorigenic growth properties of ES cell lines (17). As such we sought to investigate whether direct loss of KDM1A recapitulates these findings in vitro. shRNA retroviral knockdown of KDM1A significantly impaired the anchorage-independent growth of A673 and TTC-466 (95.2/18.5-fold reduction in colony number, respectively) ES cell lines in soft agar (Fig. 1E; Supplementary Fig. S1B/C). Correspondingly, silencing of KDM1A significantly reduced the proliferative capacity of both cell lines as shown through IncuCyte live imaging analysis. Following 108 hours of growth, confluency of A673 iLuc shRNA control cells was 100% versus 16.8% for KDM1A knockdown cells (P < 0.0001; Fig. 1F). Analogous results were observed for TTC-466 cells suggesting that both small-molecule inhibition (SP-2509) and direct depletion of KDM1A impair EWS/ETS-driven oncogenic transformation.
Targeted KDM1A inhibition and EWS/FLI knockdown induce similar transcriptional signatures
As we previously established that SP-2509 comprehensively disrupts the transcriptional signature of EWS/ETS (17), we next investigated whether genetic depletion of KDM1A mirrors the gene expression profiles of both SP-2509 treatment and EWS/FLI knockdown. RNA-seq analysis following KDM1A shRNA knockdown in A673 cells revealed 880 genes significantly induced (KDM1A repressed genes; >2Log2 fold) and 952 genes significantly downregulated (KDM1A-activated genes; Supplementary Fig. S2A). Ingenuity pathway analysis (IPA) revealed that the top canonical-enriched pathways associated with KDM1A-activated and -repressed genes were ethanol degradation IV (P = 1.89 × 10−9) and hepatic stellate cell activation (P = 7.51 × 10−9), respectively. We next assessed the similarity between the KDM1A and EWS/FLI-driven transcriptional signature. Using our previous published A673 EWS/FLI knockdown RNA-seq data set (17), significant overlap (χ2; P < 0.0001) between these two transcriptional profiles was observed (Supplementary Fig. S2B/C). In total, 195 repressed genes (11.2%) were commonly induced upon KDM1A-EWS/FLI knockdown, which were highly enriched for hepatic stellate cell activation (P = 3.56 × 10−11) and HMGB1 signaling (P = 7.69 × 10−08). Having established a similarity between EWS/FLI and KDM1A transcriptional profiles, we next compared the global signature of small-molecule KDM1A blockade against KDM1A shNRA knockdown. χ2 analysis revealed a statistically significant overlap (P < 0.0001) between genes activated/repressed following A673 treatment with SP-2509 (2 μmol/L) and KDM1A knockdown. In total, 162 genes were commonly induced upon SP-2509/KDM1A knockdown (KDM1A repressed genes) with 107 genes commonly downregulated (KDM1A activated genes; Supplementary Fig. S2D/E). IPA of these cohorts established that granulocyte/agranulocyte adhesion and diapedesis, HMGB1 signaling, hepatic stellate cell activation, and Il-8 signaling were highly enriched for KDM1A knockdown/SP-2509–induced genes, with genes significantly repressed following SP-2509/KDM1A treatment associated with G2–M DNA damage checkpoint regulation (P = 1.26 × 10−4). We alternatively compared SP-2509 and KDM1A transcriptional profiles using GSEA with a more stringent cutoff of a 3-fold for the KDM1A knockdown gene set (Supplementary Fig. S2F). GSEA showed significant correlation (P < 0.001) for genes both up- and downregulated upon KDM1A knockdown in the SP-2509–regulated gene list. Together, these results suggest that genes regulated by KDM1A in ES correlate with EWS/FLI transcriptional function, with SP-2509 expression profiles strongly mimicking KDM1A genetic depletion.
ES cell lines are hypersensitive to reversible KDM1A inhibition
We next sought to elucidate the molecular basis of SP-2509–induced apoptotic cytotoxicity. Seventeen ES cell lines with varying STAG2 and TP53 mutational status were treated with SP-2509 for 72 hours with viability determined through CellTiter-Glo. Notably, innate resistance to SP-2509 was not observed, with all cell lines sensitive to the cytotoxic effects of KDM1A inhibition, IC50 range 81 to 1,593 nmol/L (Mean IC50 0.621 ± 0.095 μmol/L; Table 1). Importantly, primary hMSC cells, a putative ES cell of origin and TIP5 fibroblasts remained unaffected at these low concentrations, SP-2509 IC50 of >4 μmol/L and 18.2 μmol/L, respectively. To confirm the superior cytotoxicity of SP-2509 in ES, sensitivity of osteosarcoma, rhabdomyosarcoma, prostate (25), and neuroblastoma (26) cell lines to SP-2509 was also assessed. ES cell lines were significantly more sensitive (range, 2.1- to 11.4-fold greater) to SP-2509 compared with all other cancer subtypes (Supplementary Fig. S3A).
Genomic characteristics and SP-2509 sensitivity of the ES cell line cohort
To determine the antitumor activity associated with catalytic inhibition of KDM1A, a panel of ES cell lines was also treated with the irreversible KDM1A inhibitors tranylcypromine (72 hours) and its more specific next-generation analogue derivative GSK-LSD1 (144 hours). Unlike SP-2509, tranylcypromine and GSK-LSD1 both covalently bind the FAD pocket of KDM1A. Strikingly, the low-dose rapid cytotoxicity profile of SP-2509 was not recapitulated with either compound (Supplementary Fig. S3B). GSK-LSD1 concentrations exceeding 300 μmol/L were required to achieve 50% antiproliferative activity in 7/9 cell lines and >400 μmol/L for tranylcypromine in 6/6 cell lines.
To address whether the pronounced disparity in cellular activity between irreversible and reversible KDM1A inhibitors were due in part to intracellular drug accumulation, whole cell and nuclear uptake of tritiated [3H] GSK-LSD1 was assessed. GSK-LSD1–resistant ES cell lines (TC252, A673) and two additional cell lines HL60 and H1417 (AML and SCLC, respectively) previously shown to be hypersensitive to the cytotoxic effects of GSK-LSD1 both in vitro (EC50 <20 nmol/L) and in vivo (18) were treated with [3H] GSK-LSD1 (1 μmol/L). No significant difference in both whole-cell lysate and nuclear uptake between ES and hypersensitive GSK-LSD1 control cell lines was observed following incubation with [3H] GSK-LSD1 (5 minute and 2 hour treatment; Supplementary Fig. S3C). Our data suggests that failure of GSK-LSD1 to impart a cytostatic effect in ES cell lines cannot be attributed to the inability of this KDM1A inhibitor to accumulate within the cell.
Finally, in order to assess whether irreversibly blocking KDM1A enzymatic activity has any effect on SP-2509 sensitivity, A673 and TC252 cells were pretreated with GSK-LSD1 for 48 hours prior to treatment with SP-2509 alone or in combination with GSK-LSD1 for an additional 72 hours. No significant change in SP-2509 sensitivity was observed following pretreatment with GSK-LSD1 in both A673 and TC252 cells, suggesting that blockade of KDM1A enzymatic activity is not required for SP-2509 to induce cytotoxicity in ES (Supplementary Fig. S3D). We propose that ES cells are more susceptible to disruption of the KDM1A-EWS/ETS-chromatin regulatory complex (17) rather than enzymatic irreversible inactivation through direct covalent modification of the KDM1A FAD pocket.
SP-2509 but not GSK-LSD1 significantly impairs the proliferative growth of ES cell lines
To confirm that GSK-LSD1 has no effect on the proliferative growth of ES cell lines, A673, TC252, and ES-2 cells were treated with SP-2509, GSK-LSD1, or the ES chemotherapeutic agent doxorubicin for 96 hours. As expected, doxorubicin significantly impaired the proliferative growth of all three cell lines as early as 18 hours after treatment (ES-2, P = 0.020; Fig. 2) with strong induction of apoptosis observed as quantified through caspase-3/7 activity. Similarly, SP-2509 concentrations exceeding 0.5 μmol/L in hypersensitive cell lines (A673, TC252) and 1 μmol/L in the sensitive cell line ES-2, significantly reduced proliferative capacity as early as 42 hours (A673 and ES-2). In contrast, no significant difference in cellular growth or caspase-3/7 induction between GSK-LSD1 and vehicle control treatment was observed for all three cell lines across all time points, underscoring the superiority of reversible KDM1A inhibitors for the treatment of ES.
SP-2509 but not GSK-LSD1 significantly reduces the proliferative capacity of ES cell lines. A, SP-2509 hypersensitive (A673, TC252) and sensitive (ES-2) Ewing sarcoma cell lines treated with the indicated concentrations of SP-2509, GSK-LSD1, doxorubicin or vehicle control (DMSO) for 96 hours. Proliferative capacity and induction of apoptosis (caspase-3/7 activity) was measured in real-time through IncuCyte ZOOM live-cell imaging. Data represent mean ± SEM confluency or green object count from three independent experiments. B, Statistical analysis of proliferation or caspase-3/7 induction compared with vehicle control cells (ns, not significant). C, Representative phase contrast and green fluorescence (caspase-3/7) images following 72 hours of treatment with the indicated agents.
TP53/STAG2/CDKN2A status and basal KDM1A expression levels do not confer SP-2509 cytotoxicity
Several recent large-scale sequencing studies have confirmed the low mutational burden in ES tumors, with only recurrent deleterious mutations reported in STAG2 (21.5%) and TP53 (6.2%) and focal deletions at the CDKN2A locus (5, 6, 27). As such, we investigated whether STAG2, TP53, and CDKN2A status correlates with SP-2509 sensitivity. Despite that 3/5 hypersensitive SP-2509 cells lines (IC50 < 500 nmol/L; SK-N-MC, A673, and SK-ES-1) harbor null/mutant TP53, no correlation between TP53 mutational status and SP-2509 sensitivity was observed (P = 0.256). This TP53 finding was additionally confirmed in our expanded panel of neuroblastoma, osteosarcoma, prostate, and rhabdomyosarcoma cell lines (Supplementary Fig. S3E). Similarly, neither STAG2, CDKN2A status nor EWS/ETS translocation partner correlated with SP-2509 sensitivity (P = 0.887, P = 0.884, and P = 0.866, respectively; Fig. 3A). We next investigated whether basal KDM1A or KDM1B mRNA expression levels correlate with SP-2509 sensitivity. Although overexpression of KDM1A (>2-fold) was observed in 14/17 cell lines, no correlation between basal KDM1A or KDM1B mRNA levels was demonstrated (R2 = 0.039 and R2 = 0.013, respectively; Fig. 3B and C). In addition, although strong protein expression of KDM1A was evident across the entire cell line cohort, no significant association between SP-2509 sensitivity and KDM1A, KDM1B, FLI, and p53 protein levels as well as cell doubling time (28) or whether the cell line was obtained from chemotherapy-naïve or treated tumors was observed (Fig. 3D; Supplementary Fig. S3F/G).
Basal KDM1A mRNA and protein levels do mediate SP-2509 sensitivity in ES. A, Lack of correlation between TP53/STAG2/CDKN2A status, EWS/ETS translocation partner, and SP-2509 sensitivity (IC50). B, Relative mRNA expression (fold change) of KDM1A, KMD1A exon-2a, and KDM1B in the ES cell line cohort. Expression normalized to IMR90 cells. Data represent mean mRNA expression ± SD from two independent experiments. C, Lack of correlation between KDM1A, KMD1A exon-2a and KDM1B mRNA expression levels and SP-2509 sensitivity (IC50). D, Western blot analysis of KDM1A, KDM1B, FLI, and TP53 protein expression in the cell line cohort. Cell lines are ranked in order of SP-2509 sensitivity (IC50). Data represent mean SP-2509 IC50 ± SEM from three independent experiments. Open symbols denote most (SK-N-MC) and least (TC-71) SP-2509–sensitive cell lines.
Four unique KDM1A isoforms generated through either single or double insertion of two alternatively spliced exons (2a and 8a) in the unstructured amino-terminus and amino oxidase domains of KDM1A, respectively, have been identified (29). As retention of the exon 8a splice variant is restricted to cells of neuronal lineage, mRNA expression levels of exon 2a (60 bp) were evaluated in our cell line panel. Overexpression of exon 2a was observed in 7/17 cell lines and strongly correlated with total KDM1A mRNA levels (R2 = 0.902; Fig. 3B; Supplementary Fig. S3H). Similar to total KDM1A, no association between basal KDM1A exon 2a mRNA expression levels and SP-2509 sensitivity was observed (R2 = 0.004; Fig. 3C). Our data demonstrate that regardless of KDM1A expression levels, ES cells are susceptible to the cytotoxic effects of SP-2509, thereby strengthening the broad utility of targeted KDM1A inhibition for the treatment of ES.
KDM1B expression levels after SP-2509 treatment strongly correlate with SP-2509 sensitivity
To assess whether SP-2509 treatment directly modulates KDM1A, KDM1B, and EWS/FLI mRNA and protein levels, a panel of eight ES cell lines of varying sensitivity was treated with SP-2509 (2 μmol/L) for 48 hours. Although KDM1A mRNA levels were significantly downregulated in 7/8 cell lines, maximum of 52.0% (A673), the degree of KDM1A mRNA reduction did not correlate with SP-2509 sensitivity (Supplementary Fig. S3I; R2 = 0.022). In contrast, KDM1B mRNA levels were significantly induced following SP-2509 treatment in 7/8 cell lines with the greatest induction observed in the hypersensitive SP-2509 cell line A673 (2.6-fold; Supplementary Fig. S3J). KDM1B mRNA levels were strongly predictive of SP-2509 sensitivity (R2 = 0.562), suggesting that hypersensitive ES cell lines positively induce KDM1B following SP-2509 compared with sensitive cell lines. Similarly, KDM1B protein levels were significantly induced (>25%) only in hypersensitive cell lines 44 and 55 hours after SP-2509 treatment (maximum 1.72-fold induction, TC252) (Supplementary Fig. S6A). EWS/FLI mRNA levels were modulated in 4/5 cell lines, with a maximum reduction of 80.5% observed in TC32 cells (Supplementary Fig. S3K). However, this striking effect was only observed at the transcript level, as only a significant reduction in KDM1A but not EWS/FLI protein levels was observed across multiple cell lines, maximum reduction of 26.1% (A673; Supplementary Fig. S3L).
To exclude the possibility that SP-2509′s modulation of the KDM1A, KDM1B, and EWS/FLI signature was simply a nonspecific response to treatment with cytotoxic agents, we sought to determine whether other compounds such as GSK-LSD1 and doxorubicin also induce similar gene expression patterns (Supplementary Fig. S4A/B). No consistent changes (mRNA and protein) in KDM1A, KDM1B, or EWS/FLI were observed across multiple ES cell lines following treatment with doxorubicin or GSK-LSD1. Finally, to elucidate whether small-molecule modulation of EWS/FLI protein levels also affects KDM1A and KDM1B expression, A673, TC252, and ES-2 cells were also treated with the nucleoside analogue ARA-C (cytosine arabinoside), a known modulator of EWS/FLI protein but not transcript levels (30). Although a significant reduction in EWS/FLI protein levels was achieved, no consistent change across all three cell lines was observed for KDM1A and KDM1B protein levels (Supplementary Fig. S4C/D). Together, these findings suggest that the ability of SP-2509 to modulate KDM1A, KDM1B, and EWS/FLI expression levels is a direct effect, and not simply a cellular stress response to cytotoxic agents.
KDM1B mediates SP-2509 sensitivity in hypersensitive ES cell lines
On the premise that induction of KDM1B following SP-2509 treatment is significantly correlated with SP-2509 cytotoxicity, we next addressed the role of KDM1B in mediating SP-2509 sensitivity. Stable retroviral knockdown of KDM1B in the hypersensitive SP-2509 cell line A673 and sensitive cell line EWS-502 was achieved with two KDM1B shRNAs (shRNA #3 and shRNA #7). KDM1B targeting shRNA effectively reduced endogenous KDM1B mRNA expression by a maximum of 76.1% and 86.1%, A673 and EWS-502 cells, respectively (Fig. 4A). Importantly, no significant reduction in KDM1A protein or mRNA levels was observed for both KDM1B hairpins, (Fig. 4B and C). Interestingly, knockdown of KDM1B significantly reduced the proliferative capacity of both A673 and EWS-502 cells compared with non-targeting control (iLuc). Following five-day (120 hours) IncuCyte proliferation analysis, A673 confluency for iLuc, KDM1B shRNA #3, and KDM1B shRNA #7 generated cell lines were as follows, 99.4%, 27.9%, and 41.7%. Analogous results were observed for EWS-502 (Fig. 4D). Finally, to address whether KDM1B has any effect on SP-2509 sensitivity, non-targeting and KDM1B shRNA A673 and EWS-502 cell lines were treated with SP-2509 (0–4 μmol/L). Following 72 hours of treatment, KDM1B knockdown significantly reduced the sensitivity of A673 cells to SP-2509 compared with controls by 2.3- and 4.3-fold, shRNA #3 and shRNA #7, respectively. In contrast, KDM1B knockdown had no effect on EWS-502 sensitivity (IC50: 1.386 μmol/L iLuc, 1.250 μmol/L KDM1B shRNA #3, 1.370 μmol/L KDM1B shRNA #7; Fig. 4E), advocating a potential role for KDM1B driving SP-2509 hypersensitivity.
KDM1B mediates SP-2509 sensitivity in hypersensitive ES cell lines. Relative (A) KDM1B and (B) KDM1A mRNA levels in A673 and EWS-502 cells following transduction with two unique KDM1B knockdown shRNA constructs (#3 and #7) or iLuc shRNA control. C, Representative Western blot analysis of KDM1B and KDM1A protein levels in A673 and EWS-502 cells following KDM1B knockdown. α-Tubulin was used a loading control. KDM1B densitometry normalized to α-tubulin and relative to iLuc shRNA control is depicted. D, Real-time live-cell imaging (IncuCyte ZOOM) of A673 and EWS-502 cells transduced with KDM1B and iLuc shRNA constructs. Cell proliferation measured for 120 hours. E, Cell viability analysis (CellTiter-Glo) of A673 and EWS-502 transduced with the indicated shRNA constructs and treated with SP-2509 (0–4 μmol/L) for 72 hours. Data represent mean ± SEM from three to four independent experiments. Asterisks denote statistical significance (*, P < 0.05; **, P < 0.01; ***, P < 0.001; ****, P < 0.0001).
SP-2509 hypersensitive cell lines share similar basal transcriptomic profiles
To gain additional insight into the molecular mechanisms and pathways that mediate SP-2509 cytotoxicity and sensitivity, three hypersensitive (IC50 < 300 nmol/L: A673, TC252, TC32) and three sensitive (IC50 > 900 nmol/L: EWS-502, ES-2, TC-71) ES cell lines treated with either vehicle control (DMSO) or SP-2509 (2 μmol/L) were submitted for RNA-seq analysis. Unsupervised hierarchical clustering analysis of these transcriptomes showed that those cell lines that were highly sensitive to SP-2509 shared similar basal gene expression profiles and clustered separately from the sensitive cell lines (Supplementary Fig. S5A/B). In all, 424 genes were significantly differentially expressed between the two groups with HIST1H2AJ (Histone H2A type 1-J) being the most significantly differentially expressed gene (P = 1.156 × 10−20). IPA revealed that the top four canonical pathways associated with this core subset of genes included actin cytoskeleton signaling (P = 0.002), CdC42 signaling (P = 0.0161), D-mannose degradation (P = 0.020), and RhoA signaling (P = 0.035).
To complement and extend these analyses, we next sought to identify which core subset of genes were specifically induced and repressed in our hypersensitive but not sensitive ES cell lines following SP-2509 treatment (Supplementary Fig. S5C/D). In total, 44 genes were uniquely induced only in the hypersensitive cohort (2-fold) with ETS1, HIF1A-AS2, RASEF, KCNK9, and CRIM1 being the top five most highly expressed genes. Interestingly, 14/44 (31.8%) of these genes overlapped with known EWS/FLI-regulated targets (17). Similarly, 53 genes were specifically repressed in our hypersensitive cell lines (2-fold) with HIST1H2BM, HIST1H2AJ, CEACAM5, ICAM3, and MORC1 being the top five most highly repressed genes following SP-2509 treatment. Of these genes, 43.4% (23/53) overlapped with known EWS/FLI-regulated targets.
SP-2509 induces endoplasmic reticulum stress in ES cell lines
On the basis that SP-2509 resistance was not observed in our ES cell line cohort, we next examined which genes are commonly induced and repressed across all cell lines to mediate SP-2509 cytotoxicity. In total, 103 genes were commonly induced (>2-fold) with 82 genes universally repressed following SP-2509 treatment (Fig. 5A). Interestingly, 30/103 (29.1%) and 33/82 (40.2%) of these genes overlapped with known EWS/FLI-regulated targets, corroborating our previous findings that SP-2509 reverses the EWS/ETS-driven transcriptional program in ES (Supplementary Table S2). In addition, the majority of genes modulated by SP-2509 were protein coding (80.6% and 69.5%), followed by long coding RNAs (17.5% and 15.9% up/downregulated, respectively; Fig. 5B). Notably, tRNAs (11.0%) were also only specifically repressed following SP-2509 treatment.
SP-2509 cytotoxicity is mediated through the ER-stress response pathway. A, Heat map representation of genes commonly induced/repressed following SP-2509 treatment (2 μmol/L, 48 hours) across 6 ES cell lines compared with vehicle control (DMSO). Three independent replicates per treatment and cell line are depicted. Scale, mean-centered rlog-transformed expression. B, Distribution of SP-2509–regulated genes according to class. C, IPA of the top canonical pathways associated with SP-2509–induced genes (>2-fold increase from vehicle control). D, Normalized RNA-seq mRNA expression of ER-stress response genes following SP-2509 treatment across the ES cell line cohort. (E and F) Relative mRNA expression levels of HSPA5, DDIT3, ERN1, and spliced XBP1 in TC252 cells following treatment with DMSO, thapsigargin (50 nmol/L) or SP-2509 (2 μmol/L) for the indicated time periods. Data represent mean expression ± SEM from triplicate reactions. G, PCR analysis of unspliced (XPB1-US) and spliced XBP1 (XPB1-S) in TC252 cells following treatment as in E. “M” and “W” denote marker and water control, respectively. RPL19 housekeeper was used as a loading control.
IPA of SP-2509–induced genes uncovered the following five highly enriched pathways: unfolded protein response (P = 2.15 × 10−9), endoplasmic reticulum (ER) stress pathway (P = 1.56 × 10−6), EIF2 signaling (P = 1.86 × 10−3), adipogenesis pathway (P = 2.23 × 10−3), and NRF2-mediated oxidative stress response (P = 8.09 × 10−3; Fig. 5C). Indeed, the following genes associated with the ER-stress response pathway; DDIT3, DNAJBP, PPP1R15A, ERN1, HSPA5, XBP1, and SEL1L were significantly induced across all cell lines following SP-2509 treatment (maximum 25.2-fold increase from vehicle control) (Fig. 5D). For this reason, we focused on delineating the role of ER stress in mediating SP-2509 cytotoxicity.
To confirm that SP-2509 activates the ER-stress response pathway, TC252 cells were treated with SP-2509 (2 μmol/L), vehicle control or thapsigargin (50 nmol/L), a chemical agent known to activate the UPR stress response by blocking ER calcium ATPase pumps, leading to depletion of ER calcium stores (31, 32). Significant induction of HSPA5 (BiP) (71.4-fold) and the terminal UPR proapoptotic transcription factor DDIT3 (CHOP; 126.9-fold) occurred rapidly after thapsigargin treatment compared with vehicle control (3 hours; Fig. 5E) with expression weakening (>10-fold induction) after 24 hours of treatment. In contrast, a time-dependent increase in HSPA5 and DDIT3 following SP-2509 treatment was observed, with maximum thapsigargin comparable transactivation of HSPA5 and DDIT3 detected at 48 hours (62.2-fold and 74.3-fold, respectively). Regardless of SP-2509 sensitivity, similar protein induction of ERN1, PERK, and DDIT3 was observed across A673, TC252, ES-2, and EWS-502 cell lines following both thapsigargin and SP-2509 treatment. In contrast, marked protein induction of BiP/HSPA5 (>25% from DMSO control) 48 and 55 hours after SP-2509 treatment was only observed in our hypersensitive cell lines (Supplementary Fig. S6A). As BiP is a major ER-chaperone protein critical for protein quality control of the ER and initiation of the UPR, the lack of protein induction in sensitive SP-2509 cell lines may account for the divergent SP-2509 IC50s documented across our cell line panel. Finally, mRNA levels of 12 known UPR pathway genes were assessed to determine whether hypersensitive cell lines are basally primed for greater robust apoptotic responses following SP-2509 treatment. No significant difference in any UPR mediators was observed between hypersensitive and sensitive SP-2509 cell lines (Supplementary Fig. S6B).
It is well established that during the initiation of the UPR, dissociation of HSPA5 from ERN1 (IRE1α) results in dimerization and activation of ERN1. The endoribonuclease domain of activated ERN1 subsequently cleaves a 26-nucleotide intron from XBP1, causing a frame shift enabling translation and generation of a basic leucine zipper family transcription factor, spliced XBP1. Spliced XBP1 consequently binds numerous UPR and ER-stress response genes leading to transcription of ER-chaperone proteins. Following thapsigargin treatment in TC252 cells, rapid induction of XBP1 spliced mRNA was observed (46.5-fold) as early as 3 hours after treatment, which rapidly declined over a 48-hour period to near pretreatment levels (Fig. 5F). In contrast, similar to HSPA5, DDIT3, and ERN1, significant induction of spliced XBP1 was only observed 24 hours after SP-2509 treatment with peak induction observed 48 hours after drug exposure (22.3-fold increase from vehicle control). PCR analysis of the spliced (263 bp) and unspliced (289 bp) XBP1 variants following thapsigargin and SP-2509 corroborated these findings. The slow kinetic splicing of XBP1 was only visualized 24 and 48 hours after SP-2509 treatment. In contrast, spliced XBP1 was observed as early as 3 hours after thapsigargin treatment, with strong diminution of splicing observed with prolonged drug exposure (Fig. 5G). Taken together, these results reveal that SP-2509 induces the UPR and ER-stress response, although at a slower kinetic rate than known UPR inducing agents.
Knockdown of ERN1 significantly attenuates the cytotoxic activity of SP-2509
The UPR comprises three parallel signaling branches: EIF2AK3/PERK, ATF6α, and ERN1/IRE1α. It is well accepted that the ER-localized transmembrane receptor ERN1 is one of the first prosurvival arms activated upon ER stress, with subsequent attenuation of activity facilitating apoptosis in response to chronic unresolved stress. As ERN1 expression was notably induced across all six ES cell lines following SP-2509 treatment (Fig. 5A/D) and is responsible for the cleavage of XBP1, we next investigated whether stable retroviral knockdown of ERN1 mitigates SP-2509 cytotoxicity. Both ERN1 targeting shRNAs (shRNA #1 and #2) effectively reduced endogenous ERN1 mRNA and protein levels (maximum of 80.7% and 61.5%, respectively; Supplementary Fig. S6C—S6E). Following 72-hour treatment with SP-2509 (0–4 μmol/L), shRNA knockdown of ERN1 (shRNA #2) significantly increased the concentration of SP-2509 required to reduce viability by 50% (2.06-fold increase). These findings corroborate our RNA-seq study and validate the role of the UPR response in instigating SP-2509 cytotoxicity in ES cells.
Both oxidative stress and nutrient deprivation are well-established stimulators of the UPR response. In order to further evaluate how SP-2509 initiates this adaptive pathway, reactive oxygen species (ROS), HIF1α protein levels, and 2-deoxyglucose-6-phosphate (2DG6P) uptake following SP-2509 treatment were assessed (Supplementary Fig. S6F-H). No induction of ROS or HIF1α was observed following SP-2509 treatment. In contrast, a significant dose-dependent reduction in 2-deoxyglucose-6-phosphate (2DG6P) uptake was observed in both A673 and TC252 cells following treatment with SP-2509 (48/72 hours), suggesting that SP-2509 may engage the accumulation of misfolded proteins through nutrient deprivation.
Discussion
Despite global efforts, only minor improvement in overall survival rates for patients afflicted with relapsed and metastatic ES has been achieved over the past 40 years, emphasizing the need to move away from traditional nonspecific agents to drugs that specifically target the underlying drivers of Ewing sarcomagenesis. We clearly show that ES cell lines and tumors express elevated levels KDM1A, with shRNA depletion abolishing the proliferative and tumorigenic properties of ES cells in vitro. Elevated KDM1A expression was associated with poor overall survival, consistent with several other solid cancers including ovarian (14), lung (11), esophageal squamous cell carcinoma (13), and neuroblastoma (12). Thus, in addition to EWS/FLI, KDM1A may represent an additional Achilles' heel of ES, which unlike EWS/FLI can be easily targeted through small-molecule blockade.
In an effort to favorably affect patient outcome, the premise of this study was to elucidate mechanisms of sensitivity to the KDM1A-reversible small-molecule SP-2509, a noncompetitive KDM1A inhibitor we have previously shown to reverse the EWS/ETS transcriptional signature and impair multiple EWS/ETS-associated oncogenic phenotypes both in vitro and in vivo (17). We reveal that proposed biomarkers such as TP53/STAG2/CDKN2A status in addition to basal KDM1A mRNA/protein levels do not correlate with SP-2509 sensitivity. However, our results uncovered a previously unappreciated role for KDM1B/LSD2, the only mammalian homologue of KDM1A, in SP-2509 hypersensitivity. In addition, we demonstrate that induction of the UPR response mediates SP-2509 cytotoxicity in ES, a pathway previously unlinked with KDM1A blockade.
Notably, shRNA targeting of KDM1B significantly reduced the proliferative capacity of ES cell lines and mitigated the sensitivity of hypersensitive SP-2509 cells to small-molecule KDM1A inhibition (Fig. 4). Furthermore, SP-2509–specific induction of KDM1B following treatment was also predictive of SP-2509 hypersensitivity (Supplementary Fig. S3). The reported SP-2509 IC50 for KDM1B is 11.2 μmol/L, well above the SP-2509 IC50 for our entire ES cell line cohort, validating that changes in KDM1B mRNA levels after treatment are not due to direct inhibition of KDM1B. In contrast, no consistent upregulation of KDM1B was observed following 2 μmol/L treatment with GSK-LSD1 (KDM1B IC50 93.7 μmol/L).
Although KDM1A and KDM1B only share 33% amino acid sequence homology in the amine oxidase domain, they are both capable of demethylating the same enzymatic substrate, H3K4me1/2 (33). Emerging evidence has also supported that these two histone demethylases possess distinct functions, primarily due to the fact that KDM1B contains an N-terminal CW-type zinc finger domain, which is absent in KDM1A and lacks a tower domain that is instrumental for binding to the corepressor protein CoREST (33). Genome-wide mapping has also revealed that KDM1B predominantly associates with the gene bodies of actively transcribed genes, but is markedly absent from promoters, which is in direct contrast to KDM1A, which largely associates with promoter regions of genes (34). In addition, functional classification of KDM1B-associated polypeptides suggests that the composition of the KDM1B complex is notably distinct from KDM1A, with KDM1B-associated proteins linked to splicing, DNA replication and damage repair, nucleosome remodeling, and histone modification. Despite this knowledge, the precise role of KDM1B in cancer remains relatively unexplored. Amplification of KDM1B has been observed in breast (35), bladder urothelia carcinoma (36), and neuroendocrine prostate cancer (37), with stable shRNA silencing of KDM1B significantly reducing the colony formation capacity of breast cancer cell lines (38). Although high expression of KDM1B is associated with shorter overall survival in prostate cancer (39), we show that KDM1B expression is not predictive of event-free or overall survival in ES (; Supplementary Fig. S6I/J; refs. 6, 24).
In 2010, Fang and colleagues used MS/MS proteomic analysis to identify KDM1B associated proteins, one of which was PARP-1 (34). It is well established that the EWS/FLI fusion present in ES drives and maintains PARP-1 expression, which in turn further promotes transcriptional activation by EWS/FLI (40, 41). It may be plausible that induction of KD1MB in hypertensive SP-2509 cell lines directly affects the interaction between KDM1B/PARP-1 thereby enhancing the execution of PARP-1–driven apoptosis following SP-2509 treatment. Indeed, our results demonstrated that knockdown of KDM1B only in the hypersensitive SP-2509 cell line A673 but not EWS-502 (sensitive) significantly reduced the cytotoxic effects of SP-2509. Together, our findings lay the foundation for future studies to elucidate the specific role of KDM1B in mediating SP-2509 sensitivity.
To date, the precise mechanism by which SP-2509 drives ES cells toward an apoptotic fate remains unanswered. Significant induction of caspase-3/7 activity, a marker for apoptosis was clearly evident across multiple cell lines as early as 15 hours after SP-2509 treatment (Fig. 2). Previous studies have reported that KDM1A is required for cell-cycle progression as it colocalizes with centrosomes in interphase and spindle poles during mitosis to permit proper chromosome segregation (42). Notably, downregulation of KDM1A leads to abnormal centrosome duplication and significantly impairs nuclear pore complex assembly subsequently leading extended telophase (43). In contrast, our RNA-seq analysis of six Ewing cell lines after SP-2509 treatment revealed a strong enrichment for genes associated with the unfolded protein (UPR) and ER stress response (Fig. 5).
A number of pathophysiological extracellular perturbations (disturbed Ca2+ homeostasis, oxidative stress, metabolic changes, acidosis, and nutrient deprivation) lead to the accumulation of unfolded/misfolded proteins in the ER subsequently triggering ER stress (44). The UPR is a double-edge sword, initially promoting coordinated cytoprotective cascades to reestablish cellular homeostasis and preserve ER protein folding fidelity. However, under conditions of chronic unmitigated stress, UPR signaling shifts away from prosurvival adaptive programs with ER-stress–induced apoptosis ensued. The transcription factor C/EBP homologous protein CHOP (also known as DDIT3/GADD153) is a common point of convergence for all three arms of the UPR as evidenced by the discovery that the DDIT3 promoter contains binding sites for the major UPR induces ATF4, ATF6, and XBP1. In addition, it is well established that DDIT3 promotes UPR apoptotic cell death following unresolved ER stress (45). Indeed, strong induction of DDIT3 was observed across all six ES cell lines; approximately 24 hours after SP-2509 treatment.
Our A673 KDM1A knockdown RNA-seq data revealed that several UPR genes including DDIT3, DNAJB9, and PPP1R15A are upregulated following KDM1A depletion. Similarly, DNAJB9, XBP1, and HSPA5 were identified as direct EWS/ETS repressed targets (17). In addition, to classify regulated proteins associated with EWS/FLI1, recent network i-TRAQ analysis of ES cell lines following EWS/FLI knockdown revealed XBP1 as the highest ranked pathway (46). These findings demonstrate that UPR targets are directly repressed by both KDM1A and EWS/ETS. Interestingly, we show that SP-2509 significantly attenuates glucose uptake in ES cell lines. This observation is supported by multiple lines of evidence advocating the essential role of KDM1A in coordinating glycolytic and mitochondrial metabolism in cancerous cells. Inactivation of KDM1A (chemical and genetic depletion) has been shown to suppress glucose uptake and glycolytic activity, with concurrent activation of mitochondrial respiration (47–50). Collectively, our preliminary findings clearly show for the first time an association between the UPR and the ability of SP-2509 to induce apoptosis in ES cell lines.
In summary, our findings add to the growing body of evidence highlighting the oncogenic role of KDM1A in tumors such as ES and validate the broad therapeutic utility of targeted KDM1A inhibition for the treatment of this aggressive pediatric malignancy. We demonstrate that ES cells are exceptionally dependent upon KDM1A for growth and are exquisitely sensitive to small molecule reversible but not irreversible KDM1A blockade (SP-2509). Epigenetic regulators, such as SP-2509, represent an exciting class of pharmacologic agents for the treatment of ES, with clinical formulations of SP-2509 (SP-2577/Seclidemstat) entering phase I clinical testing in 2018. For this reason, further studies are required to define mechanism(s) of SP-2509 resistance, a key unanswered question, which will be essential to guide combinatorial treatment methods to circumvent patient drug resistance and define more selective inclusion criteria required for KDM1A inhibitor clinical trials.
Disclosure of Potential Conflicts of Interest
S. Sharma reports receiving commercial research grant from and has ownership interest (including stock, patents, etc.) in Salarius Pharmaceuticals. S.L. Lessnick has ownership interest (including stock, patents, etc.) in and is a consultant/advisory board member for Salarius Pharmaceuticals. No potential conflicts of interest were disclosed by the other authors.
Authors' Contributions
Conception and design: K.I. Pishas, E.R. Theisen, M.C. Beckerle, S.L. Lessnick
Development of methodology: K.I. Pishas, E.R. Theisen, I.L Pop
Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): C.D. Drenberg, K.M. Johnson, B.D. Crompton, E.R. Lawlor, F. Tirode, J. Mora, O. Delattre
Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): K.I. Pishas, C.D. Drenberg, C. Taslim, R.S. Saund, I.L Pop, B.D. Crompton, E.R. Lawlor, S. Sharma, S.L. Lessnick
Writing, review, and/or revision of the manuscript: K.I. Pishas, C.D. Drenberg, J. Mora, M.C. Beckerle, S. Sharma, S.L. Lessnick
Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): I.L Pop, O. Delattre
Study supervision: D.F. Callen, S.L. Lessnick
Acknowledgments
We acknowledge experimental support from Iftekhar Showpnil and Julia Selich-Anderson. This study is dedicated to the memory of Tom Wood and Brad Neilsen.
S.L. Lessnick was supported by the NIH/NCI grants R01 CA140394 and R01 CA183776, the NCTN Grant U10 CA180886, Human Specimen Banking Grant U24 CA114766, and the NCTN Statistics and Data Center Grant U10 CA 180899 of the Children's Oncology Group from the NCI/NIH (Bethesda, MD, USA). Additional support for research for S.L. Lessnick was provided by a grant from the WWWW (QuadW) Foundation, Inc. (www.QuadW.org) to the Children's Oncology Group. K.I. Pishas acknowledges financial support from the University of Adelaide Florey Medical Research Foundation Clinical Cancer Research Fellowship, NHMRC CJ Martin Overseas Biomedical Fellowship (APP1111032), and the Alex's Lemonade Stand Young Investigators Grant (APP37138). M.C. Beckerle acknowledges financial support from CureSearch for Children's Cancer and the Huntsman Cancer Institute Cancer Center Support Grant (P30-CA04214).
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
Footnotes
Note: Supplementary data for this article are available at Molecular Cancer Therapeutics Online (http://mct.aacrjournals.org/).
- Received April 10, 2018.
- Revision received May 15, 2018.
- Accepted July 2, 2018.
- Published first July 11, 2018.
- ©2018 American Association for Cancer Research.