Next Article in Journal
Iridoid Derivatives as Anticancer Agents: An Updated Review from 1970–2022
Next Article in Special Issue
Routine EWS Fusion Analysis in the Oncology Clinic to Identify Cancer-Specific Peptide Sequence Patterns That Span Breakpoints in Ewing Sarcoma and DSRCT
Previous Article in Journal
Combination, Modulation and Interplay of Modern Radiotherapy with the Tumor Microenvironment and Targeted Therapies in Pancreatic Cancer: Which Candidates to Boost Radiotherapy?
Previous Article in Special Issue
Current State of Immunotherapy and Mechanisms of Immune Evasion in Ewing Sarcoma and Osteosarcoma
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Exploiting the Stemness and Chemoresistance Transcriptome of Ewing Sarcoma to Identify Candidate Therapeutic Targets and Drug-Repurposing Candidates

by
Elizabeth Ann Roundhill
1,*,
Pan Pantziarka
2,
Danielle E. Liddle
1,
Lucy A. Shaw
1,
Ghadeer Albadrani
1 and
Susan Ann Burchill
1,*
1
Children’s Cancer Research Group, Leeds Institute of Medical Research, St James’s University Hospital, Beckett Street, Leeds LS9 7TF, UK
2
Anticancer Fund, Brusselsesteenweg 11, 1860 Meise, Belgium
*
Authors to whom correspondence should be addressed.
Cancers 2023, 15(3), 769; https://doi.org/10.3390/cancers15030769
Submission received: 7 November 2022 / Revised: 8 January 2023 / Accepted: 17 January 2023 / Published: 26 January 2023

Abstract

:

Simple Summary

Ewing sarcoma is a cancer arising most frequently in teenagers and young adults. For many patients, outcomes are the same today as they were 30 years ago, emphasising the need for more effective treatments to eradicate the cells responsible for progression and relapse. These cells responsible for progression and relapse have been identified using assays that evaluate functional characteristics and expression of cell surface markers. For the first time, we reveal ABCG1 as an additional potential cell surface marker of progression. In rare cancers like Ewing sarcoma, commercial development of new drugs is seldom a priority, reflecting the small number of patients and lack of well-characterised molecular subtypes. Therefore, rather than creating new drugs, which can take 20 to 30 years, repurposing of existing drugs may be an efficient cost-effective strategy to accelerate novel molecularly targeted therapy into clinical trials for patients with Ewing sarcoma. To identify candidate molecular targets, we have used a combination of functional assays and transcriptomic analyses to characterize the cells responsible for progression and relapse. We have then applied a bespoke in silico pipeline to find drugs with known safety profiles that bind to these targets. In the future, after preclinical validation of efficacy and specificity in Ewing sarcoma, some of these drugs may be assessed as combination treatments in clinical trials, with the goal of improving outcomes.

Abstract

Outcomes for most patients with Ewing sarcoma (ES) have remained unchanged for the last 30 years, emphasising the need for more effective and tolerable treatments. We have hypothesised that using small-molecule inhibitors to kill the self-renewing chemotherapy-resistant cells (Ewing sarcoma cancer stem-like cells; ES-CSCs) responsible for progression and relapse could improve outcomes and minimise treatment-induced morbidities. For the first time, we demonstrate that ABCG1, a potential oncogene in some cancers, is highly expressed in ES-CSCs independently of CD133. Using functional models, transcriptomics and a bespoke in silico drug-repurposing pipeline, we have prioritised a group of tractable small-molecule inhibitors for further preclinical studies. Consistent with the cellular origin of ES, 21 candidate molecular targets of pluripotency, stemness and chemoresistance were identified. Small-molecule inhibitors to 13 of the 21 molecular targets (62%) were identified. POU5F1/OCT4 was the most promising new therapeutic target in Ewing sarcoma, interacting with 10 of the 21 prioritised molecular targets and meriting further study. The majority of small-molecule inhibitors (72%) target one of two drug efflux proteins, p-glycoprotein (n = 168) or MRP1 (n = 13). In summary, we have identified a novel cell surface marker of ES-CSCs and cancer/non-cancer drugs to targets expressed by these cells that are worthy of further preclinical evaluation. If effective in preclinical models, these drugs and drug combinations might be repurposed for clinical evaluation in patients with ES.

1. Introduction

Ewing sarcoma (ES) arises in the bone or soft tissue [1], most frequently presenting in young people aged 10–25 years [2]. Standard of care treatment, including a combination of multi-agent chemotherapy, surgery and radiotherapy, has improved outcomes for some patients [1]. However, around 40% of these patients develop multidrug-resistant (MDR) metastatic disease [2,3,4,5], leading to relapse and survival rates normally associated with metastatic disease (5-year survival 10% [5,6,7]). Late relapses and chemotherapy-induced morbidities are additional enduring burdens for patients, their families and carers. A major unmet clinical need is therefore the introduction of molecularly targeted therapy to minimise treatment-related toxicity and improve outcomes.
Progression and relapse are driven by subpopulations of cells capable of self-renewal and migration within tumours that are resistant to current treatments. These cancer stem-like cells (CSCs) have been identified in a range of adult and paediatric solid tumours based on expression of cell surface markers, most frequently CD133 (also known as prominim-1) [8,9,10,11,12]. Ewing sarcoma cancer stem-like cells (ES-CSCs) have been isolated based on expression of CD133 [13,14,15]. However, distinct CD133-negative CSCs are present in some cancers [8,9,10,11,12], including Ewing sarcoma [15]. Therefore, additional approaches, including formation of clones from single cells and generation of three-dimensional (3D) spheroids, have been used to improve the identification of ES-CSC [16,17] and CSC in other cancer types [18,19,20].
In this study, we have investigated the completeness of CD133 as a cell surface marker of ES-CSCs and identified ABCG1 as an additional marker of these cells using 3D spheroids and self-renewal from single cells. The expression and prognostic potential of this ABC transporter protein has been evaluated in patient-derived cells and tumours in the online dataset GSE17618. Genes that regulate pluripotency, stemness and the MDR ABC transporter proteins have been identified by comparing the transcriptomes of substrate adherent two-dimensional (2D) non-CSC and 3D spheroid-derived ES-CSCs. We combined genes that were differentially expressed in ES cells grown in 3D spheroids and cells grown in 2D with genes previously reported in patient-derived ES-CSCs [17] to pinpoint candidate molecular targets expressed by ES-CSCs. We then developed a bespoke pipeline to identify small-molecule inhibitors of these targets for further preclinical studies. If these drugs are effective in preclinical models of ES, they may in the future be accelerated into clinical trials for evaluation in patients. This repurposing strategy seeks to reuse existing licensed drugs to treat new indications and is a complementary strategy to de novo drug development [21,22].

2. Materials and Methods

2.1. Cell Culture

ES cell lines (A673, RD-ES, SKES-1, SK-N-MC, TC-32 and TTC-466) were cultured as previously described [23] and purchased from the American Type Culture Collection, Manassas, VA, except for the following cells that were kind gifts: TC-32 cells from Dr. J. Toretsky (Division of Pediatrics, University of Maryland, Baltimore, MD, USA), the TTC-466 cells from Dr. P. Sorenson (British Columbia Children’s Hospital, Vancouver, BC, Canada). Primary ES cell cultures and daughter ES-CSCs were cultured as previously described [17]. The embryonic stem cell (ESC) culture SHEF-4 (RRID:CVCL_9791) was a gift from Professor P. Andrews, Centre for Stem Cell Biology, University of Sheffield, Sheffield, UK [24] and used as a positive control for CD133 and ATP binding cassette subfamily G member 1 (ABCG1). The glioblastoma cell line (T98G) was a gift from Professor M. Knowles, University of Leeds, and was a positive control for multidrug resistance-associated protein 1 (MRP1, gene name = ABCC1) [23]. The human embryonic fibroblast line KMST-6 (cultured in MEM containing 10% FBS and 2 mM glutamine) and Jurkat cells (cultured in RPMI 1640 containing 10% FBS and 2 mM glutamine) were gifts from Dr E. Morrison, University of Leeds, and used as the positive control and to generate the calibration curve, respectively, for the measurement of telomere length.

2.2. Western Blotting

Western blotting (WB) was performed as previously described [23]. Equal loading of proteins was confirmed using α-tubulin [23] or β-actin (0.4 µg/mL, A5441, Sigma-Aldrich, Paisley, UK). For the detection of CD133, protein extracts were heated at 95 °C for 5 min before cooling on ice and WB (CD133, 1 µg/mL, W6B3C1, Miltenyi Biotech, Surrey, UK). After incubation with primary antibodies (MRP1 [23] or ABCG1 (1:1000, ab36969; Abcam Plc., Cambridge, UK)) and secondary antibodies [23], proteins were visualised by image capture (at different exposure times depending on the intensity of signal) and quantified using the Li-cor Odyssey infrared imaging system (Li-cor Biosciences, Lincoln, NE, USA).

2.3. Flow Cytometry

2.3.1. Cell Surface Expression of CD133

Cells (5 × 105) were incubated in FcR blocking buffer (Miltenyi Biotech) and either anti-CD133/2 (4.5 µg/mL, clone 293C3, Miltenyi Biotech) or the isotype control IgG2b-PE (4.5 µg/mL, Miltenyi Biotech) antibodies in the dark at 4 °C for 10 min. Cells (10,000 per sample) were then analysed by flow cytometry using the FACSCalibur (BD Biosciences, Berkshire, UK).

2.3.2. Co-Expression of ABCG1 and CD133

SK-N-MC cells were incubated in normal goat serum (1:10, Dako, Agilent Technologies, Santa Clara, CA, USA) in BD Perm/Wash™ Buffer for 30 min at 4 °C. Cells were then incubated with CD133 (anti-CD133/2, 4.5 µg/mL, clone 293C3) and/or ABCG1 rabbit polyclonal antibody (100 µg/mL, PA5-13462, Thermo-Fisher Scientific, Paisley, UK) in BD Perm/Wash™ Buffer for 1h at 4 °C in the dark. Control cells were incubated with IgG2b-PE (4.5 µg/mL, Miltenyi Biotech) or rabbit IgG isotype control antibody (100 µg/mL, Dako). Cells were then incubated with the secondary antibody (2 μg/mL goat anti-rabbit IgG FITC, A31556, Thermo-Fisher Scientific) in BD Perm/Wash™ Buffer for 30 min at 4 °C in the dark. Cells (10,000 per sample) were analysed by flow cytometry using the CytoFLEX (Beckman Coulter, High Wycombe, UK).

2.4. Self-Renewing Ability

Growth of progeny from a single cell as an adherent culture was determined as previously described [17,18]. A single cell (Poisson distribution probability of λ < 1 = 0.9) was seeded into each well of 10 PrimariaTM 96-well plates (Corning) and the number of wells containing ≥5 cells was recorded after 21 days by light microscopy (Olympus CKX41). Where possible, single cell self-renewing cell populations were propagated to establish daughter cell cultures; these are subsequently referred to as ES-CSCs. To examine colony forming efficiency in soft agar, a single cell suspension (1.8 × 105 cells) in cell-specific media containing 0.3% agar (Aldrich, Poole, UK), was overlaid on a solid agar bed (0.5% agar in media). After 21 days, colonies were stained with 8 mM iodonitrotetrazolium chloride (made up in ddH2O (w/v); Sigma-Aldrich) for 16 h and colony number counted by light microscopy. Colony forming efficiency = [the number of colonies formed in the field of view/number of cells seeded] × 100.
For spheroid formation, a single cell was seeded into each well of an ultra-low attachment plate (Corning, UK) in cell line-specific media. The number of spheroids at 21 days was recorded to calculate the spheroid forming efficiency (SFE); SFE = [number of wells containing a spheroid/the number of wells seeded with a single cell] × 100. Spheroids were imaged by light microscopy (Olympus CKX41). For reverse transcriptase quantitative polymerase chain reaction (RTqPCR), Western blot, immunocytochemistry (ICC) and flow cytometry, spheroids formed from single cells were collected at 21 days.

2.5. Magnetically Activated Cell Sorting for the Enrichment of CD133-Positive Cells

ES cells (1 × 108) were incubated with 300 µL of CD133 microbeads and 100 µL of FcR blocking buffer (Miltenyi Biotech). CD133-positive cells were isolated using LS columns (Miltenyi Biotech). CD133-negative cells were depleted of labelled (CD133-positive) cells by passing the cells through two LD columns (Miltenyi Biotech). Cell surface expression of CD133 was confirmed by flow cytometry immediately after separation: >90% of the CD133-positive selection expressed CD133, and CD133 was expressed by <5% of the CD133-negative selected cells. CD133-positive A673 and TC-32 cells remained positive over 15 passages.

2.6. RNA Expression of Markers of Pluripotency and Differentiation, the Wnt Signalling Pathway and ABC Transporter Proteins

RNA was extracted using the RNeasy Micro Kit (Qiagen, Düsseldorf, Germany) and RNA quantity and quality measured using the Nanodrop and Agilent 2100 bioanalyzer. RNA (with a RIN > 8) was converted to cDNA by reverse transcription [23]. The mRNA expression was evaluated using the TaqMan® Human Stem Cell Pluripotency Array, TaqMan® Human Wnt Pathway and TaqMan® Human ABC Transporter Array (Applied Biosystems, Waltham, MA, USA) [23]. To allow direct comparison across the 3 array platforms target Ct values were normalised using the global mean [25,26,27]. ABC transporter and pluripotency mRNAs (n = 140, excluding endogenous control mRNAs) more highly expressed than other mRNAs (Ct values < 25), expressed (Ct values 25–35) and not expressed (Ct values >35, [18,28,29]) in 3D spheroids from TTC 466 and SK-N-MC cells were identified. Unique and shared mRNAs in each group were analysed using the Search Tool for Retrieval of Interacting Genes/Proteins (STRING) database (http://string-db.org, [17,30]) to identify GO terms. Significant differences in mRNA expression were determined using Linear Models for Microarray Data (LIMMA) [18]. Target mRNAs were validated using individual RTqPCR assays if the Q value was <0.1, there was at least a mean Ct difference of >2 between the compared populations and Ct values were <35 [18]. For mRNA validation, total RNA (20 ng) was reverse-transcribed and cDNA added to the PCR mix containing sequence specific reverse and forward primers and probe for PPIA (the housekeeping gene) or ABCG1 (Thermo Fisher Scientific; ABCG1 Hs00245154_m1) and 1 × TaqMan® Universal PCR Master Mix (Thermo Fisher Scientific). RNA expression was calculated using the comparative Ct method [31].
ABCG1 transcripts were characterised in total RNA extracted from SK-N-MC cells grown as 2D cultures or 3D spheroids, sequenced and aligned as previously described [17]. Normalised read counts were identified in total RNA sequencing data using DESeq2 [32], adjusted p values of <0.01 were considered significant.

2.7. Immunohistochemistry (IHC)

Formalin-fixed paraffin-embedded (FFPE) sections (4 µm) of spheroids were deparaffinised in xylene (2 min) and rehydrated before antigen retrieval in citric acid buffer (10 mM in ddH2O, pH6, [33]). Endogenous peroxidases were blocked using 3% hydrogen peroxide for 5 min and endogenous biotin, biotin receptors and avidin-binding sites blocked using an avidin/biotin blocking kit (Invitrogen, Waltham, MA, USA). Sections were incubated for 1 h with ABCG1 primary antibody (10 µg/mL, PA5-13462, Thermo-Fisher Scientific) or rabbit IgG control (10 µg/mL, Dako) at room temperature, followed by incubation with the secondary antibody. Sections were then incubated with streptavidin–peroxidase (Abcam Plc.), followed by DAB substrate (Dako) for 10 min and nuclei counterstained with haematoxylin.

2.8. siRNA Knockdown of ABCG1

Cells (5 × 104) were seeded and allowed to adhere overnight. Media were replaced with Accell siRNA Delivery Media (B-005000-500, Dharmacon, Lafayette, CO, USA) alone or containing ABCG1 siRNA (1.5 µM, E-008615-00-0005, SMARTpool:Accell ABCG1 siRNA, containing 4 siRNAs targeting exon 6, present in all canonical and novel transcripts (Dharmacon) or non-targeting control siRNAs (1.5 µM, D-001910-10-20, Dharmacon) for 72 h. ABCG1 knockdown was confirmed by RTqPCR.

2.9. Apoptosis

Cells (5 × 104) were harvested following treatment with ABCG1 and control siRNA for 24–72 h and apoptosis measured by annexin V and propidium iodide (PI) labelling of cells (annexin V–FITC apoptosis detection kit, BD Biosciences, [34]).

2.10. Statistical Analyses

Differences in proliferation and viable cell numbers were log transformed and analysed by linear regression. Differences in the gradients of each plot were compared using the extra sum of squares F test. For all other experiments, data were analysed by analysis of variance (ANOVA followed by a Bonferroni’s or Tukey’s post hoc test), or a non-parametric Mann–Whitney or unpaired two-tailed t-test. Correlations were determined using a Pearson’s correlation coefficient (r). p values of <0.05 were considered significantly different. Statistical analyses were performed using GraphPad PRISM 7.03 (GraphPad Software, San Diego, CA, USA).

2.11. Identification of Candidate Drugs for Prioritised Molecular Targets

Gene lists were analysed using the STRING database. To identify drugs reported to target these molecular candidates, we interrogated DrugBank (version 5.1.9 [35]) and DGI database (DGIdb, version v4.2.0 [36]) to access information on FDA-approved drugs and their molecular targets. Additional drug-target information was derived from the literature-based repurposing drugs in oncology (ReDO) database [37]. A polypharmacology approach identified targets for each drug candidate. Those targets previously associated with ES using the Open Targets Platform [38] and DisGeNET [39] were employed as the identification source. Data from DGIdb and Open Targets characterise the strength of target–disease and drug–target associations. Using DGIdb, the interaction score is computed from: (publication count + data source count) × (average known gene partners for all drugs/known gene partners for candidate drug) × (average known drug partners for all genes/known drug partners for target gene). The Dir DGI score therefore incorporates data on the number of molecular targets of the drug, the number of drugs that target a gene and the number of publications and data sources supporting the association. The Dir Assoc score is generated from the Open Targets database and is a measure of the relationship between the molecular target and “cancer” as the disease term.
Licensed cancer drugs were identified from the list of drug candidates addressing >1 ES molecular target using the Cancer Drugs database [40]. ES trials were identified using clinicaltrials.gov, EU Clinical Trials Register and WHO International Clinical Trials Registry Platform. For non-cancer drugs previously identified as oncological repurposing candidates (ReDO database), clinical trial activity in any cancer was extracted using the ReDO_Trials database of active repurposing trials in oncology [41]. For these repurposing candidates, a support score was calculated from data in the ReDO database to characterise the range of data available illustrating the anticancer effects of the drugs (e.g., in vitro, in vivo, case reports, observational data and clinical trial data).

3. Results

3.1. ES Cell Lines Produce Spheroids and Clones from a Single Cell

All ES cell lines (6/6) formed clones or spheroids from a single cell when cultured in soft agar or on ultra-low attachment plates (Table 1, Figure 1), consistent with previous reports that ES cell cultures contain a self-renewing cell population [8,9,10,11,12,13,14,17].
The formation of single cell-derived clones and spheroids in soft agar (29 ± 2%, p < 0.05) and ultra-low attachment conditions (75 ± 7%, p < 0.05; Table 1, Figure 1) was most efficient in the SKES-1 cell line. Spheroid formation was 100% in all cell lines when combining 100 cells or more (Table 1). Larger spheroids (>400 µm), produced cell line-specific morphology (Figure 1). SKES-1 spheroids produced a relatively uniform sphere of disseminated or dissociated cells, consistent with spheroid formation in a range of solid tumour cell types [42,43,44,45,46]. For the first time, we identified 3D projections developing from the central core in five of six ES cell lines (Figure 1). The biological significance of these projections requires further investigation.

3.2. CD133 Identifies Some Self-Renewing Drug-Resistant ES-CSCs

CD133 protein was detected in all ES cell lines, with the exception of SK-N-MC cells (Figure 2A and Supplementary Data S1). A673 and TC-32 CD133-positive cells formed significantly more colonies in soft agar than the CD133-negative SK-N-MC cells (Supplementary Data S1–S3). However, there was no difference in proliferation, cell cycle status or telomere length across the two cell lines (Supplementary Data S1–S3), phenotypes frequently associated with CSCs and self-renewing ability [47]. Moreover, expression and activity of the multidrug-resistant protein MRP1 was increased in CD133-positive TC-32 cells, but not the A673 CD133-positive population (Supplementary Data S1–S3). The shared CSC phenotype of both CD133-positive and CD133-negative ES populations demonstrates that CD133 expression alone is not sufficient to enrich for all the ES-CSC population.
Consistently with the cellular origin of ES and the high level of stemness markers expressed by these tumour cells [48,49,50], no significant differentially expressed genes associated with pluripotency, ABC transporter and Wnt signalling pathways were identified in CD133-positive and CD133-negative cells (Figure 2B (TC-32) and Figure 2C (A673), Supplementary Data S4). These observations are also consistent with the premise that CD133 can be used to identify some cells with characteristics of ES-CSCs, although pathways classically associated with the CSC phenotype in other cancer types are also expressed in CD133-negative ES cells. To investigate the CD133-independent ES-CSC phenotype, we went on to investigate ES-CSCs enriched through spheroid formation in ES cells with no or low CD133 expression.

3.3. Gene Expression Profile of CD133 Low or Negative ES 3D Spheroids and 2D Cultures

As TTC 466 and SK-N-MC cells had low or no detectable CD133 protein (Figure 2A), we compared the transcriptome of TTC 466 and SK-N-MC 3D spheroids and 2D cultures (Figure 3A). Consistently with the premise that spheroid formation is a feature of CSCs, expression of the stemness markers LEFTB and LIN28 were significantly increased in SK-N-MC spheroids compared to cells in 2D culture, whereas expression of LAMA1, COL2A1, ACTC, GCG, SEMA3A and PODXL were significantly decreased (Q value < 0.1, Figure 3B,D, Supplementary Data S5A). In TTC 466 spheroids, no markers were significantly increased; however, expression of FLT1 and GATA6 was significantly decreased compared to cells in 2D culture (Figure 3C,E, Supplementary Data S5A).
Expression of 20 ABC transporter genes was significantly different in SK-N-MC spheroids compared to the cells in 2D cultures, although only ABCG1 and CFTR were differentially expressed greater than twofold (Figure 4A–C). However, ABC transporter genes, including ABCG1, were not significantly differentially expressed in TTC 466 cells grown in 2D or as 3D spheroids (Figure 4D–F). Comparison of the mRNAs in SK-N-MC and TTC 466 spheroids revealed that 35% of genes were undetected in both populations and an additional seven (5%) and six (4%) unique mRNAs were undetected in TTC 466 and SK-N-MC spheroids respectively (Figure 4G). Fifty-six percent of mRNAs were detected in both spheroid populations and 8% of these were highly expressed (Figure 4G). The six mRNAs uniquely expressed in TTC 466 spheroids are involved in cell differentiation (GO:0030154), which may explain why increased expression of pluripotency mRNAs was not observed in TTC 466 spheroids. Previous studies have shown knockdown of EWSR1-FLI1 decreases the expression of stemness markers and stem-like properties of ES cells [51]. Whether the lack of shared pluripotency mRNAs increased in the TTC 466 and SK-N-MC spheroids reflects the different transcriptional activators of TTC 466 (EWSR1-ERG) and SK-N-MC cells (EWSR1-FLI1) remains to be seen.
The decrease in expression of the CFTR gene was not validated using RTqPCR or Western blot (results not shown). However, the increase in ABCG1 expression was confirmed in SK-N-MC spheroids at the mRNA (p = 0.0015, Figure 4H) and protein (Figure 4I) level. ABCG1 is heterogeneously expressed in all ES cell lines (6/6; Supplementary Data S5B). By IHC, ABCG1 expression was confirmed in the outer proliferating region of 3D SK-N-MC spheroids but not the hypoxic region or inner necrotic core (Figure 4J; [34]), suggesting expression of ABCG1 may be regulated by hypoxia and have a functional role in ES cells under these conditions.

3.4. Functional Role and Characterisation of ABCG1 in ES-CSCs

Knockdown of ABCG1 using siRNA had no effect on the viability or apoptosis of SK-N-MC cells in culture (10 ± 0.8% in non-targeting and 10 ± 1.1% in ABCG1 siRNA-treated cells, p > 0.05; Figure 5A). There was also no effect on proliferation (0.93 ± 0.1 and 0.94 ± 0.05 in non-targeting and ABCG1 siRNA treated cells respectively, p > 0.05) or self-renewing ability from a single cell in soft agar (60% in non-targeting and 58% in ABCG1 siRNA-treated cells). However, after knockdown of ABCG1 spheroid production was significantly reduced; 65 ± 3% in non-targeting and 33 ± 5% in ABCG1 siRNA-treated cells (p < 0.0001, n = 10). These data suggest that although ABCG1 may not have a role in homeostasis of ES cells in 2D, it may affect cell–cell interactions and components of the tumour microenvironment important in the development of 3D spheroids and possibly tumours. This hypothesis requires further investigation.
Analysis of ABCG1 RNA in SK-N-MC spheroids and 2D cultures using total RNA sequencing, revealed expression of 11 previously reported transcripts (www.ensembl.org, accessed on 27 April 2016; Figure 5B,C) and two novel ABCG1 RNA species (Figure 5C). Canonical ABCG1 generates the multiple transcripts by alternative splicing (Figure 5C). Novel transcript 1 is most similar to transcript 4 (ENST00000450121.5), both sequences missing exon 5 (Figure 5C). However, novel transcript 1 included exons 8–15 and had an extended 3′UTR region predicted to produce a larger molecular weight protein than transcript 4. The RNA sequence of novel transcript 2 is most similar to transcript 5 (ENST00000361802.6), the additional 3′UTR sequence unlikely to produce a unique protein product. Expression of both novel transcripts at the RNA level was increased above canonical protein producing sequences (transcripts 1, 3–8) in SK-N-MC cells grown as spheroids compared to 2D cultures, this was most significant for novel transcript 1 (Figure 5D). The increase in ABCG1 gene expression in SK-N-MC cells from spheroids (Figure 4G,H), might then reflect expression of novel transcript 1. Further studies are required to investigate this novel ABCG1 transcript and its role in ES-CSCs.
Since ABCG1 mRNA was increased in ES-CSCs from SK-N-MC spheroids, we examined its co-expression with CD133. There was no significant difference in the percentage of cells expressing ABCG1 or CD133 in SK-N-MC cells from 3D spheroids or 2D cultures (Figure 6A,B). Consistently with the increase in ABCG1 in protein extracts from 3D spheroids (Figure 4H), the level of ABCG1 protein per cell was greater in cells from 3D spheroids compared to 2D cultures (2.5 ± 0.8-fold increase, p < 0.05, Figure 6C). Levels of CD133 per cell were also increased (6.3 ± 1 fold increase, p < 0.05, Figure 6C), consistent with the higher expression of ABCG1 protein in TC-32 CD133-positive cells compared to CD133-negative cells (Figure 6D). Expression of ABCG1 RNA was confirmed at the protein level in primary patient-derived ES cells (Figure 6E). However, there was no correlation between the percentage of CD133-positive cell lines or patient-derived cells and progeny-producing ability measured using the colony formation assay in soft agar (R2 < 0.1) or self-renewing ability from a single cell seeded on low adherent or adherent plates (R2 = 0.002; Figure 6F), suggesting in patient-derived cells CD133 may not identify cell populations with a self-renewing phenotype.
Interrogation of the publicly available GSE17618 RNA dataset revealed high expression of CD133, but not ABCG1, was associated with a more than threefold risk of an event and poor outcome (Supplementary Data S6), consistent with the hypothesis that high expression of CD133 identifies ES-CSCs that may be responsible for progression and relapse in some patients (Figure 6G).

3.5. Identification of Candidate Drugs for ES-CSC Molecular Targets

Candidate small-molecule inhibitors to 62% (13/21) of the unique molecular targets increased in 3D spheroids from TTC 466 and SK-N-MC cells and ES-CSCs (Figure 7A,B) were identified using our bespoke pipeline (Figure 7C). Of the 279 drugs directed against the targets, 66 (24%) are licensed cancer drugs and 213 (76%) are non-cancer drugs (Figure 7C and Supplementary Data S7). The majority of the drugs (202/279; 72%) target one of the two ATP-binding cassette (ABC) drug efflux proteins p-glycoprotein (n = 168, includes n = 10 targeting p-glycoprotein and additional proteins) or MRP1 (n = 13), while 21 drugs target both p-glycoprotein and MRP1. Forty eight of the 279 (17%) drugs hitting the targets have been evaluated in oncology trials, whereas the majority of drugs have been evaluated in non-cancer drug trials (n = 986; Figure 7C, Table 2, Supplementary Data S8).
POU5F1/OCT-4 interacted with 10/13 of the identified stemness and chemoresistance genes (Figure 7D), suggesting shared functional protein–protein association networks. Five of these 10 gene products (POU5F1/OCT4, c-KIT, CAV1, ITGB1 and CD44) interact (Figure 7D), regulating vital cellular processes (Supplementary Data S9). All these proteins are highly expressed in ES-CSCs and so may represent candidate prognostic biomarkers and/or therapeutic targets in ES. We are currently investigating these possibilities. Using our customised pipeline, we identified two FDA approved small-molecule inhibitors, allopurinol and phenytoin, that target POU5F1/OCT4 and are used to treat gout and control seizures in epilepsy, respectively (Table 2, https://www.dgidb.org/genes/POU5F1#_interactions, accessed on 3 January 2023). These two licensed drugs have the highest DGIdb interaction scores of the three compounds that interact with POU5F1/OCT4, inhibit a range of other cancer-relevant molecular targets and are also currently in trials as cancer therapeutics, making them attractive candidate drugs for further in vitro analysis for ES. Several inhibitors targeting multiple receptor tyrosine kinases (TKIs) implicated in the pathogenesis of sarcomas including ES were also identified [52,53], including regorafenib, which is being evaluated in combination with chemotherapy (NCT02085148/2013-003579-36; NCT04055220/REGOSTA; 2021-005061-41/INTER-EWING-1; Table 2).

4. Discussion

For the first time, we have identified expression of the membrane-associated ABC transporter protein ABCG1 on the surface of human ES cells. Levels of ABCG1 were increased in SK-N-MC cells forming 3D spheroids compared to cells in 2D culture. Expression of ABCG1 was particularly high in the outer rim of spheroids, suggesting a structural or transport function between cells within and outside the spheroid. Consistent with this hypothesis, knockdown of ABCG1 reduced the ability of ES cells to bind to each other and produce 3D spheroids, reminiscent of the ABCG1-dependent regulation across cellular and intracellular membranes [54]. ABCG1 protein was detected at lower levels in hypoxic cells surrounding the necrotic centre of spheroids, consistent with observations in mouse colon adenocarcinoma spheroids where expression of ABCG1 and hypoxia-inducible factor 1α are inversely correlated [55]. In contrast to studies in lung cancer [56] and normal haemopoietic stem cells [57], decreased ABCG1 had no direct effect on the proliferation or apoptosis of ES cells. Decreased ABCG1 is reported to promote apoptosis through increased expression of the endoplasmic reticulum (ER) stress proteins GRP78 and CHOP [58]. The ability to tolerate high levels of ER stress conferred by the EWSR1-ETS oncogene [59] may explain why knockdown of ABCG1 does not induce apoptosis in these ES cells. ABCG1 expression is higher in several cancer types compared to normal tissue [29,56], expression being associated with higher-grade tumours [60], metastasis [61] and poor response to chemotherapy [62]. In agreement with these observations, ABCG1 expression is increased in drug-resistant, self-renewing osteosarcoma cells [18], a second bone cancer characterised by recurrent disease and acquired MDR. Further studies are required to establish what effect this lipid ABC transporter protein has on plasma membrane organisation and recruitment of signalling processes that may regulate cell fate and contribute to tumour development, evasion and metastasis. Single-nucleotide polymorphisms (SNPs) of ABCG1 in the first cytoplasmic domain (within intron 2) are associated with survival of patients with non-small-cell lung cancer [63]. Mutagenesis studies have highlighted the importance of this region for effective trafficking of ABCG1 to the plasma membrane and regulation of cholesterol efflux [64,65]. However, the clinical relevance of the 17 SNPs identified within the cytoplasmic region of ABCG1 (amino acids 178–195, exon 5, www.nvbi.nlm.nih.gov/snp, accessed on 13 July 2022) has yet to be established. In ES 3D spheroids, we identified two novel ABCG1 transcripts, the most abundant transcript (transcript 1) missing exon 5 of the cytoplasmic region. The cellular functional and clinical significance of this deletion and these transcripts requires further investigation. These hypotheses and data are consistent with the premise that ABCG1 identifies ES cells that are capable of surviving chemotherapy and may be responsible for progression and relapse in some patients. We are currently investigating the expression and functional role ABCG1 mRNAs and proteins in ES.
In agreement with the premise that CD133 can identify some cells with a CSC phenotype [14,66], the number of CD133-positive cells and level of CD133 expression was higher in ES cells from 3D spheroids compared to cells grown in 2D. CD133-positive cells shared common characteristics of cancer stem-like cells including increased colony formation from a single cell. However, in patient-derived cells there was no correlation between CD133 expression and self-renewing ability from a single cell. This highlights the difference between established ES cell lines and cultures more recently derived from patient tumours and supports the premise that CD133 does not identify all ES cell populations with a self-renewing phenotype. We are currently investigating the difference in the genotype and phenotype of ES cell lines and patient-derived cells. Although CD133 is reported to be a biomarker of stem cells, this is controversial, which may arise from heterogeneity and uncertainty about its physiological role(s) [67]. In contrast to studies in metastatic melanoma [68], we found no enrichment of the canonical Wnt pathway in CD133-positive ES cells compared to CD133-negative cells [69] and the level of CD133 did not predict outcomes. The increased drug resistance of TC-32 CD133-positive cells compared to CD133-negative cells might instead be effected through greater levels of the drug efflux protein MRP1, overexpression of which induces resistance to chemotherapy [23] and high membrane expression predicts a worse clinical outcome [27]. This suggests the need for further studies on the intrinsic and acquired resistance mechanisms in ES.
Increasing evidence suggests that ES can arise in neurally derived mesenchymal stem cells (MSCs; [70]) or in cells of the neural crest [71]. The interaction between the permissive cellular environment and the EWSR1-ETS tumour-specific chimeric transcription factors leading to cellular transformation and ES is poorly understood, although is no doubt important illustrated by the high incidence of ES in Europeans compared to Africans [72,73,74]. Rewiring of the transcriptome and epigenome, in addition to rare oncogenic driver events including STAG2, TP53 and CDKN2A are also known contributing factors [75]. Interestingly, five of six cell lines produced 3D progeny with projections similar to those displayed by tissue organoids [76], not previously reported in bone cancer. In contrast, the SKES-1 cells generated a cellular core with surrounding dissociating cells reminiscent of ES spheroids from CD133-positive STA-ET 8.2 cells [15]. Heterogeneity of spheroid compactness, growth and stability over time may reflect the amount of extracellular matrix components produced by the different cell lines [77], suggesting these models may be useful tools to evaluate the efficacy of therapeutic candidates. In this study, we found that the ES-CSC phenotype was associated with the expression of ABC transporter proteins and stemness features. However, because of the bystander stemness pathways associated with the cell of origin, unravelling the functional molecular mechanisms leading to Ewing sarcomagenesis requires further investigation.
We identified POU5F1 as a promising candidate therapeutic target, interacting with 10 of the prioritised targets. The protein product of this gene (octamer-binding transcription factor 4, OCT4) regulates several functional characteristics of CSCs, including self-renewal, survival, drug resistance, epithelial–mesenchymal transition and metastasis [78]. This is consistent with emerging roles for this protein in the tumorigenesis of adult cancers and stimulation of expression by the EWSR1 proto-oncogene in several human cancers, including ES [79]. Using our bespoke pipeline, we identified two FDA-approved drugs that are off patent and target POU5F1/OCT4, allopurinol and phenytoin. The efficacy and safety of phenytoin is currently being evaluated in a phase 2 clinical trial as part of a combination maintenance therapy in patients with clinically advanced sarcoma including ES [80] and in patients with metastatic pancreatic cancer [81]. A trial combining allopurinol and mycophenolate mofetil with chemotherapy in patients with relapsed small-cell lung cancer is expected to start recruitment later this year (the CLAMP trial, NCT05049863). We are currently investigating the effect of allopurinol and phenytoin alone and in combination with additional novel drugs and standard of care chemotherapy [82] in ES.
The majority of drugs we identified target the ABC transporter proteins p-glycoprotein and/or MRP1, which contribute to MDR by reducing the level of drug within cells. Several ABC transporter proteins have been described in ES, including MRP1, p-glycoprotein, ABCG1, ABCF1, ABCA6 and ABCA7 [83]. P-glycoprotein has been most frequently studied, although it does not predict outcome [27,84,85,86]. Despite the development of third-generation inhibitors of p-glycoprotein, to date, these drugs are of limited or no clinical value, and the majority of trials have been stopped due to unacceptable toxicity. In contrast, high membrane expression of MRP1 in tumours at diagnosis predicts reduced event-free and overall survival for patients [27]. Interestingly, low levels of ABCF1 in combination with high levels of IGF2BP3 also predicts poor outcome for patients [87]. Despite the clinical failure of inhibitors to ABC transporter proteins, several targeted small molecules, including some TKIs, interact with one or more ABC transporters, suggesting inhibitors may in some cases be beneficial [88]. Increased understanding of the functional role of ABC transporter proteins, including ABCG1, in the cells responsible for progression and relapse is needed to establish which, if any, are worthwhile candidate therapeutic targets to overcome MDR.

5. Conclusions

In summary, we have evidence that ABCG1 may be a candidate biomarker that could be used to select ES patients for treatment. This hypothesis requires validation. We have identified proteins expressed by ES-CSCs that might be therapeutic targets and used our bespoke pipeline to identify repurposing drug candidates that have the potential to inhibit these targets and eradicate ES-CSCs. These drugs include FDA-approved small molecules to ABC transporter proteins and two inhibitors of POU5F1/OCT-4.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/cancers15030769/s1. Supplementary Methods. Methods relating to data described in Supplementary Data S1 are described. Supplementary Data S1. Results of experiments analysing CD133-positive and CD133-negative TC-32 and A673 populations. Supplementary Data S2 (A) Colony formation, (B) proliferation, (C) cell growth, (D) cell cycle, (E) telomere length analyses and (F) migration in A673 CD133-positive and CD133-negative cells. (G) Protein expression of MRP1 in A673 CD133-positive and -negative cells. Equal protein loading was confirmed by expression of α-tubulin. (H). Drug efflux activity of MRP1 was quantified by efflux of calcein F compared with initial fluorescence following loading with calcein AM. Results are shown as means ± SEM of three independent experiments. Calcein F efflux was compared using a non-parametric Mann–Whitney two-tailed t-test. (I) The effect of doxorubicin (3.5–224 nm) on A673 CD133-positive and -negative viable cell numbers at 48 h was quantified using the trypan blue exclusion assay. Viable cell number is presented as the percentage of cells after treatment with doxorubicin relative to the vehicle control (DMSO); results are shown as the mean ± SEM of three independent experiments. EC50 values were calculated using linear regression and compared using the extra sum of squares F test. Supplementary Data S3. (A) Colony formation, (B) proliferation, (C) cell growth, (D) cell cycle, (E) telomere length analyses and (F) migration in TC-32 CD133-positive and CD133-negative cells. (G) Protein expression of MRP1 in TC-32 CD133-positive and -negative cells. (H) Drug efflux activity of MRP1 was quantified by efflux of calcein F. (I) The effect of doxorubicin (3.5-224 nm) on TC-32 CD133-positive and -negative viable cell numbers at 48 h was quantified using the trypan blue exclusion assay. Supplementary Data S4. Expression of genes associated with differentiation, pluripotency and stemness, ABC transporter protein mRNAs and those linked with the Wnt signalling pathway in A673 and TC-32 CD133-positive and -negative cells. Supplementary Data S5. (A) Expression of genes at the RNA level associated with differentiation, pluripotency and stemness in SK-N-MC and TTC 466 cells in 2D and 3D cultures. (B) Protein expression of ABCG1 in ES cell lines (n = 6) was evaluated by Western blot. Equal protein loading was confirmed by probing the blots for GRP75. Supplementary Data S6. Prognostic significance of CD133 and ABCG1 RNA in diagnosis ES tissue. Supplementary Data S7. Complete list of drugs (n = 279) directed against the shared stemness and chemoresistance therapeutic targets. Supplementary Data S8. Drugs directed against the shared stemness and chemoresistance therapeutic targets and corresponding clinical trial identifiers. Supplementary Data S9. Statistically significant (by false discovery rate) GO terms and biological processes driven by POU5F1/OCT4, c-KIT, CAV1, ITGB1 and CD44. Refs. [18,23,31,34,89,90,91] are cited in the supplementary materials.

Author Contributions

Conceptualization, S.A.B.; methodology, E.A.R. and S.A.B.; software, E.A.R.; validation, E.A.R., P.P., D.E.L., L.A.S., G.A. and S.A.B.; formal analysis, E.A.R., P.P., D.E.L., L.A.S., G.A. and S.A.B.; investigation, E.A.R., D.E.L., L.A.S. and G.A.; resources, P.P. and S.A.B.; data curation, E.A.R., P.P., D.E.L., L.A.S., G.A. and S.A.B.; writing—original draft preparation, E.A.R. and S.A.B.; writing—review and editing, E.A.R., P.P. and S.A.B.; visualization, E.A.R. and S.A.B.; supervision, S.A.B.; project administration, S.A.B.; funding acquisition, S.A.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Yorkshire Cancer Research award L352, PPO12; Bone Cancer Research Trust award BCRT 5918; Ewing Sarcoma Research Trust award BCRT/36/13, Little Princess Trust award CCLGA 2021 09 Burchill, the Candlelighters Trust, and the University of Leeds. L.A.S. was the recipient of a University of Leeds scholarship.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki. The study was approved and performance monitored by the Integrated Research Application System (IRAS167880) and CRN Yorkshire & Humber NIHR Clinical Research Network (EDGE 79301), respectively.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The FASTQ files of sequenced ES cell lines, 2D and 3D SK-N-MC cells are available in the Research Data Leeds Repository (University of Leeds), Burchill Susan and Roundhill Elizabeth (2022): Total RNA sequencing of ES cell lines, TTC 466 and SK-N-MC 2D and 3D cultures. [Dataset]. https://doi.org/10.5518/1210.

Acknowledgments

The authors thank Kim Cass, Tom Maisey, Samantha Brownhill and Andrea Berry for technical assistance.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Balamuth, N.J.; Womer, R.B. Ewing’s sarcoma. Lancet Oncol. 2010, 11, 184–192. [Google Scholar] [CrossRef] [PubMed]
  2. Cotterill, S.J.; Ahrens, S.; Paulussen, M.; Jurgens, H.F.; Voute, P.A.; Gadner, H.; Craft, A.W. Prognostic factors in Ewing’s tumor of bone: Analysis of 975 patients from the European Intergroup Cooperative Ewing’s Sarcoma Study Group. J. Clin. Oncol. 2000, 18, 3108–3114. [Google Scholar] [CrossRef] [PubMed]
  3. Ahrens, S.; Hoffmann, C.; Jabar, S.; Braun-Munzinger, G.; Paulussen, M.; Dunst, J.; Rube, C.; Winkelmann, W.; Heinecke, A.; Gobel, U.; et al. Evaluation of prognostic factors in a tumor volume-adapted treatment strategy for localized Ewing sarcoma of bone: The CESS 86 experience. Cooperative Ewing Sarcoma Study. Med. Pediatr. Oncol. 1999, 32, 186–195. [Google Scholar] [CrossRef]
  4. Rodriguez-Galindo, C.; Navid, F.; Liu, T.; Billups, C.A.; Rao, B.N.; Krasin, M.J. Prognostic factors for local and distant control in Ewing sarcoma family of tumors. Ann. Oncol. 2007, 19, 814–820. [Google Scholar] [CrossRef] [PubMed]
  5. Stahl, M.; Ranft, A.; Paulussen, M.; Bolling, T.; Vieth, V.; Bielack, S.; Gortitz, I.; Braun-Munzinger, G.; Hardes, J.; Jurgens, H.; et al. Risk of recurrence and survival after relapse in patients with Ewing sarcoma. Pediatr. Blood Cancer 2011, 57, 549–553. [Google Scholar] [CrossRef]
  6. Barker, L.M.; Pendergrass, T.W.; Sanders, J.E.; Hawkins, D.S. Survival after recurrence of Ewing’s sarcoma family of tumors. J. Clin. Oncol. Off. J. Am. Soc. Clin. Oncol. 2005, 23, 4354–4362. [Google Scholar] [CrossRef]
  7. Huang, M.; Lucas, K. Current therapeutic approaches in metastatic and recurrent ewing sarcoma. Sarcoma 2011, 2011, 863210. [Google Scholar] [CrossRef] [Green Version]
  8. Joo, K.M.; Kim, S.Y.; Jin, X.; Song, S.Y.; Kong, D.S.; Lee, J.I.; Jeon, J.W.; Kim, M.H.; Kang, B.G.; Jung, Y.; et al. Clinical and biological implications of CD133-positive and CD133-negative cells in glioblastomas. Lab. Investig. 2008, 88, 808–815. [Google Scholar] [CrossRef] [Green Version]
  9. Ogden, A.T.; Waziri, A.E.; Lochhead, R.A.; Fusco, D.; Lopez, K.; Ellis, J.A.; Kang, J.; Assanah, M.; McKhann, G.M.; Sisti, M.B.; et al. Identification of A2B5+CD133- tumor-initiating cells in adult human gliomas. Neurosurgery 2008, 62, 505–514. [Google Scholar] [CrossRef] [Green Version]
  10. Tirino, V.; Desiderio, V.; d’Aquino, R.; De Francesco, F.; Pirozzi, G.; Graziano, A.; Galderisi, U.; Cavaliere, C.; De Rosa, A.; Papaccio, G.; et al. Detection and characterization of CD133+ cancer stem cells in human solid tumours. PLoS ONE 2008, 3, e3469. [Google Scholar] [CrossRef]
  11. Walter, D.; Satheesha, S.; Albrecht, P.; Bornhauser, B.C.; D’Alessandro, V.; Oesch, S.M.; Rehrauer, H.; Leuschner, I.; Koscielniak, E.; Gengler, C.; et al. CD133 positive embryonal rhabdomyosarcoma stem-like cell population is enriched in rhabdospheres. PLoS ONE 2011, 6, e19506. [Google Scholar] [CrossRef]
  12. Wang, J.; Sakariassen, P.O.; Tsinkalovsky, O.; Immervoll, H.; Boe, S.O.; Svendsen, A.; Prestegarden, L.; Rosland, G.; Thorsen, F.; Stuhr, L.; et al. CD133 negative glioma cells form tumors in nude rats and give rise to CD133 positive cells. Int. J. Cancer 2008, 122, 761–768. [Google Scholar] [CrossRef]
  13. Cornaz-Buros, S.; Riggi, N.; DeVito, C.; Sarre, A.; Letovanec, I.; Provero, P.; Stamenkovic, I. Targeting cancer stem-like cells as an approach to defeating cellular heterogeneity in Ewing sarcoma. Cancer Res. 2014, 74, 6610–6622. [Google Scholar] [CrossRef] [Green Version]
  14. Suva, M.L.; Riggi, N.; Stehle, J.C.; Baumer, K.; Tercier, S.; Joseph, J.M.; Suva, D.; Clement, V.; Provero, P.; Cironi, L.; et al. Identification of cancer stem cells in Ewing’s sarcoma. Cancer Res. 2009, 69, 1776–1781. [Google Scholar] [CrossRef] [Green Version]
  15. Jiang, X.; Gwye, Y.; Russell, D.; Cao, C.; Douglas, D.; Hung, L.; Kovar, H.; Triche, T.J.; Lawlor, E.R. CD133 expression in chemo-resistant Ewing sarcoma cells. BMC Cancer 2010, 10, 116. [Google Scholar] [CrossRef] [Green Version]
  16. Fujii, H.; Honoki, K.; Tsujiuchi, T.; Kido, A.; Yoshitani, K.; Takakura, Y. Sphere-forming stem-like cell populations with drug resistance in human sarcoma cell lines. Int. J. Oncol. 2009, 34, 1381–1386. [Google Scholar]
  17. Roundhill, E.A.; Chicon-Bosch, M.; Jeys, L.; Parry, M.; Rankin, K.S.; Droop, A.; Burchill, S.A. RNA sequencing and functional studies of patient-derived cells reveal that neurexin-1 and regulators of this pathway are associated with poor outcomes in Ewing sarcoma. Cell Oncol. (Dordr) 2021, 44, 1065–1085. [Google Scholar] [CrossRef]
  18. Roundhill, E.A.; Jabri, S.; Burchill, S.A. ABCG1 and Pgp identify drug resistant, self-renewing osteosarcoma cells. Cancer Lett. 2019, 453, 142–157. [Google Scholar] [CrossRef]
  19. Charafe-Jauffret, E.; Ginestier, C.; Iovino, F.; Wicinski, J.; Cervera, N.; Finetti, P.; Hur, M.H.; Diebel, M.E.; Monville, F.; Dutcher, J.; et al. Breast cancer cell lines contain functional cancer stem cells with metastatic capacity and a distinct molecular signature. Cancer Res. 2009, 69, 1302–1313. [Google Scholar] [CrossRef] [Green Version]
  20. Suzuki, E.; Chiba, T.; Zen, Y.; Miyagi, S.; Tada, M.; Kanai, F.; Imazeki, F.; Miyazaki, M.; Iwama, A.; Yokosuka, O. Aldehyde dehydrogenase 1 is associated with recurrence-free survival but not stem cell-like properties in hepatocellular carcinoma. Hepatol. Res. 2012, 42, 1100–1111. [Google Scholar] [CrossRef]
  21. Bertolini, F.; Sukhatme, V.P.; Bouche, G. Drug repurposing in oncology--patient and health systems opportunities. Nat. Rev. Clin. Oncol. 2015, 12, 732–742. [Google Scholar] [CrossRef] [PubMed]
  22. Pantziarka, P.; Verbaanderd, C.; Huys, I.; Bouche, G.; Meheus, L. Repurposing drugs in oncology: From candidate selection to clinical adoption. Semin. Cancer Biol. 2021, 68, 186–191. [Google Scholar] [CrossRef]
  23. Roundhill, E.A.; Burchill, S.A. Detection and characterisation of multi-drug resistance protein 1 (MRP-1) in human mitochondria. Br. J. Cancer 2012, 106, 1224–1233. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Aflatoonian, B.; Ruban, L.; Shamsuddin, S.; Baker, D.; Andrews, P.; Moore, H. Generation of Sheffield (Shef) human embryonic stem cell lines using a microdrop culture system. In Vitro Cell Dev. Biol. Anim. 2010, 46, 236–241. [Google Scholar] [CrossRef] [PubMed]
  25. Fischer, M.; Skowron, M.; Berthold, F. Reliable transcript quantification by real-time reverse transcriptase-polymerase chain reaction in primary neuroblastoma using normalization to averaged expression levels of the control genes HPRT1 and SDHA. J. Mol. Diagn. 2005, 7, 89–96. [Google Scholar] [CrossRef] [Green Version]
  26. Lastowska, M.; Viprey, V.; Santibanez-Koref, M.; Wappler, I.; Peters, H.; Cullinane, C.; Roberts, P.; Hall, A.G.; Tweddle, D.A.; Pearson, A.D.; et al. Identification of candidate genes involved in neuroblastoma progression by combining genomic and expression microarrays with survival data. Oncogene 2007, 26, 7432–7444. [Google Scholar] [CrossRef] [Green Version]
  27. Roundhill, E.; Burchill, S. Membrane expression of MRP-1, but not MRP-1 splicing or Pgp expression, predicts survival in patients with ESFT. Br. J. Cancer 2013, 109, 195–206. [Google Scholar] [CrossRef] [Green Version]
  28. Kodani, M.; Yang, G.; Conklin, L.M.; Travis, T.C.; Whitney, C.G.; Anderson, L.J.; Schrag, S.J.; Taylor, T.H., Jr.; Beall, B.W.; Breiman, R.F.; et al. Application of TaqMan low-density arrays for simultaneous detection of multiple respiratory pathogens. J. Clin. Microbiol. 2011, 49, 2175–2182. [Google Scholar] [CrossRef] [Green Version]
  29. Demidenko, R.; Razanauskas, D.; Daniunaite, K.; Lazutka, J.R.; Jankevicius, F.; Jarmalaite, S. Frequent down-regulation of ABC transporter genes in prostate cancer. BMC Cancer 2015, 15, 683. [Google Scholar] [CrossRef] [Green Version]
  30. Szklarczyk, D.; Gable, A.L.; Lyon, D.; Junge, A.; Wyder, S.; Huerta-Cepas, J.; Simonovic, M.; Doncheva, N.T.; Morris, J.H.; Bork, P.; et al. STRING v11: Protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019, 47, D607–D613. [Google Scholar] [CrossRef] [Green Version]
  31. Brownhill, S.C.; Taylor, C.; Burchill, S.A. Chromosome 9p21 gene copy number and prognostic significance of p16 in ESFT. Br. J. Cancer 2007, 96, 1914–1923. [Google Scholar] [CrossRef]
  32. Anders, S.; Huber, W. Differential expression analysis for sequence count data. Genome Biol. 2010, 11, R106. [Google Scholar] [CrossRef] [Green Version]
  33. Brownhill, S.; Cohen, D.; Burchill, S. Proliferation index: A continuous model to predict prognosis in patients with tumours of the Ewing’s sarcoma family. PLoS ONE 2014, 9, e104106. [Google Scholar] [CrossRef]
  34. Myatt, S.S.; Redfern, C.P.; Burchill, S.A. p38MAPK-Dependent sensitivity of Ewing’s sarcoma family of tumors to fenretinide-induced cell death. Clin. Cancer Res. 2005, 11, 3136–3148. [Google Scholar] [CrossRef] [Green Version]
  35. Wishart, D.S.; Feunang, Y.D.; Guo, A.C.; Lo, E.J.; Marcu, A.; Grant, J.R.; Sajed, T.; Johnson, D.; Li, C.; Sayeeda, Z.; et al. DrugBank 5.0: A major update to the DrugBank database for 2018. Nucleic Acids Res. 2018, 46, D1074–D1082. [Google Scholar] [CrossRef]
  36. Freshour, S.L.; Kiwala, S.; Cotto, K.C.; Coffman, A.C.; McMichael, J.F.; Song, J.J.; Griffith, M.; Griffith, O.L.; Wagner, A.H. Integration of the Drug-Gene Interaction Database (DGIdb 4.0) with open crowdsource efforts. Nucleic Acids Res. 2021, 49, D1144–D1151. [Google Scholar] [CrossRef]
  37. Pantziarka, P.; Verbaanderd, C.; Sukhatme, V.; Rica Capistrano, I.; Crispino, S.; Gyawali, B.; Rooman, I.; Van Nuffel, A.M.; Meheus, L.; Sukhatme, V.P.; et al. ReDO_DB: The repurposing drugs in oncology database. Ecancermedicalscience 2018, 12, 886. [Google Scholar] [CrossRef] [Green Version]
  38. Ochoa, D.; Hercules, A.; Carmona, M.; Suveges, D.; Gonzalez-Uriarte, A.; Malangone, C.; Miranda, A.; Fumis, L.; Carvalho-Silva, D.; Spitzer, M.; et al. Open Targets Platform: Supporting systematic drug-target identification and prioritisation. Nucleic Acids Res. 2021, 49, D1302–D1310. [Google Scholar] [CrossRef]
  39. Pinero, J.; Ramirez-Anguita, J.M.; Sauch-Pitarch, J.; Ronzano, F.; Centeno, E.; Sanz, F.; Furlong, L.I. The DisGeNET knowledge platform for disease genomics: 2019 update. Nucleic Acids Res. 2020, 48, D845–D855. [Google Scholar] [CrossRef] [Green Version]
  40. Pantziarka, P.; Capistrano, I.R.; De Potter, A.; Vandeborne, L.; Bouche, G. An Open Access Database of Licensed Cancer Drugs. Front. Pharmacol. 2021, 12, 627574. [Google Scholar] [CrossRef]
  41. Pantziarka, P.; Vandeborne, L.; Bouche, G. A Database of Drug Repurposing Clinical Trials in Oncology. Front. Pharmacol. 2021, 12, 790952. [Google Scholar] [CrossRef] [PubMed]
  42. Shi, W.; Kwon, J.; Huang, Y.; Tan, J.; Uhl, C.G.; He, R.; Zhou, C.; Liu, Y. Facile Tumor Spheroids Formation in Large Quantity with Controllable Size and High Uniformity. Sci. Rep. 2018, 8, 6837. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Zhang, Q.; Yu, S.; Lam, M.M.T.; Poon, T.C.W.; Sun, L.; Jiao, Y.; Wong, A.S.T.; Lee, L.T.O. Angiotensin II promotes ovarian cancer spheroid formation and metastasis by upregulation of lipid desaturation and suppression of endoplasmic reticulum stress. J. Exp. Clin. Cancer Res. 2019, 38, 116. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Lee, J.W.; Sung, J.S.; Park, Y.S.; Chung, S.; Kim, Y.H. Isolation of spheroid-forming single cells from gastric cancer cell lines: Enrichment of cancer stem-like cells. Biotechniques 2018, 65, 197–203. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Ivascu, A.; Kubbies, M. Rapid generation of single-tumor spheroids for high-throughput cell function and toxicity analysis. J. Biomol. Screen 2006, 11, 922–932. [Google Scholar] [CrossRef] [Green Version]
  46. Lee, C.; Lee, C.; Atakilit, A.; Siu, A.; Ramos, D.M. Differential spheroid formation by oral cancer cells. Anticancer Res. 2014, 34, 6945–6949. [Google Scholar]
  47. Ayob, A.Z.; Ramasamy, T.S. Cancer stem cells as key drivers of tumour progression. J. Biomed. Sci. 2018, 25, 20. [Google Scholar] [CrossRef] [Green Version]
  48. Ribeiro-Dantas, M.D.C.; Oliveira Imparato, D.; Dalmolin, M.G.S.; de Farias, C.B.; Brunetto, A.T.; da Cunha Jaeger, M.; Roesler, R.; Sinigaglia, M.; Siqueira Dalmolin, R.J. Reverse Engineering of Ewing Sarcoma Regulatory Network Uncovers PAX7 and RUNX3 as Master Regulators Associated with Good Prognosis. Cancers 2021, 13, 1860. [Google Scholar] [CrossRef]
  49. Cidre-Aranaz, F.; Li, J.; Holting, T.L.B.; Orth, M.F.; Imle, R.; Kutschmann, S.; Ammirati, G.; Ceranski, K.; Carreno-Gonzalez, M.J.; Kasan, M.; et al. Integrative gene network and functional analyses identify a prognostically relevant key regulator of metastasis in Ewing sarcoma. Mol. Cancer 2022, 21, 1. [Google Scholar] [CrossRef]
  50. Miller, H.E.; Gorthi, A.; Bassani, N.; Lawrence, L.A.; Iskra, B.S.; Bishop, A.J.R. Reconstruction of Ewing Sarcoma Developmental Context from Mass-Scale Transcriptomics Reveals Characteristics of EWSR1-FLI1 Permissibility. Cancers 2020, 12, 948. [Google Scholar] [CrossRef]
  51. Khoogar, R.; Li, F.; Chen, Y.; Ignatius, M.; Lawlor, E.R.; Kitagawa, K.; Huang, T.H.; Phelps, D.A.; Houghton, P.J. Single-cell RNA profiling identifies diverse cellular responses to EWSR1/FLI1 downregulation in Ewing sarcoma cells. Cell Oncol. 2022, 45, 19–40. [Google Scholar] [CrossRef] [PubMed]
  52. Grunewald, T.G.; Alonso, M.; Avnet, S.; Banito, A.; Burdach, S.; Cidre-Aranaz, F.; Di Pompo, G.; Distel, M.; Dorado-Garcia, H.; Garcia-Castro, J.; et al. Sarcoma treatment in the era of molecular medicine. EMBO Mol. Med. 2020, 12, e11131. [Google Scholar] [CrossRef] [PubMed]
  53. Wilding, C.P.; Elms, M.L.; Judson, I.; Tan, A.C.; Jones, R.L.; Huang, P.H. The landscape of tyrosine kinase inhibitors in sarcomas: Looking beyond pazopanib. Expert Rev. Anticancer Ther. 2019, 19, 971–991. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Wu, A.; Wojtowicz, K.; Savary, S.; Hamon, Y.; Trombik, T. Do ABC transporters regulate plasma membrane organization? Cell Mol. Biol. Lett. 2020, 25, 37. [Google Scholar] [CrossRef]
  55. Namba, Y.; Sogawa, C.; Okusha, Y.; Kawai, H.; Itagaki, M.; Ono, K.; Murakami, J.; Aoyama, E.; Ohyama, K.; Asaumi, J.I.; et al. Depletion of Lipid Efflux Pump ABCG1 Triggers the Intracellular Accumulation of Extracellular Vesicles and Reduces Aggregation and Tumorigenesis of Metastatic Cancer Cells. Front. Oncol. 2018, 8, 376. [Google Scholar] [CrossRef] [Green Version]
  56. Tian, C.; Huang, D.; Yu, Y.; Zhang, J.; Fang, Q.; Xie, C. ABCG1 as a potential oncogene in lung cancer. Exp. Ther. Med. 2017, 13, 3189–3194. [Google Scholar] [CrossRef] [Green Version]
  57. Yvan-Charvet, L.; Pagler, T.A.; Seimon, T.A.; Thorp, E.; Welch, C.L.; Witztum, J.L.; Tabas, I.; Tall, A.R. ABCA1 and ABCG1 protect against oxidative stress-induced macrophage apoptosis during efferocytosis. Circ. Res. 2010, 106, 1861–1869. [Google Scholar] [CrossRef] [Green Version]
  58. Xue, J.; Wei, J.; Dong, X.; Zhu, C.; Li, Y.; Song, A.; Liu, Z. ABCG1 deficiency promotes endothelial apoptosis by endoplasmic reticulum stress-dependent pathway. J. Physiol. Sci. 2013, 63, 435–444. [Google Scholar] [CrossRef]
  59. Apfelbaum, A.A.; Wrenn, E.D.; Lawlor, E.R. The importance of fusion protein activity in Ewing sarcoma and the cell intrinsic and extrinsic factors that regulate it: A review. Front. Oncol. 2022, 12, 1044707. [Google Scholar] [CrossRef]
  60. Elsnerova, K.; Mohelnikova-Duchonova, B.; Cerovska, E.; Ehrlichova, M.; Gut, I.; Rob, L.; Skapa, P.; Hruda, M.; Bartakova, A.; Bouda, J.; et al. Gene expression of membrane transporters: Importance for prognosis and progression of ovarian carcinoma. Oncol. Rep. 2016, 35, 2159–2170. [Google Scholar] [CrossRef] [Green Version]
  61. Litviakov, N.V.; Garbukov, E.; Slonimskaia, E.M.; Tsyganov, M.M.; Denisov, E.V.; Vtorushin, S.V.; Khristenko, K.; Zav’ialova, M.V.; Cherdyntseva, N.V. Correlation of metastasis-free survival in patients with breast cancer and changes in the direction of expression of multidrug resistance genes during neoadjuvant chemotherapy. Vopr. Onkol. 2013, 59, 334–340. [Google Scholar]
  62. Litviakov, N.V.; Cherdyntseva, N.V.; Tsyganov, M.M.; Slonimskaya, E.M.; Ibragimova, M.K.; Kazantseva, P.V.; Kzhyshkowska, J.; Choinzonov, E.L. Deletions of multidrug resistance gene loci in breast cancer leads to the down-regulation of its expression and predict tumor response to neoadjuvant chemotherapy. Oncotarget 2016, 7, 7829–7841. [Google Scholar] [CrossRef] [Green Version]
  63. Wang, Y.; Liu, H.; Ready, N.E.; Su, L.; Wei, Y.; Christiani, D.C.; Wei, Q. Genetic variants in ABCG1 are associated with survival of nonsmall-cell lung cancer patients. Int. J. Cancer 2016, 138, 2592–2601. [Google Scholar] [CrossRef] [Green Version]
  64. Gu, H.M.; Li, G.; Gao, X.; Berthiaume, L.G.; Zhang, D.W. Characterization of palmitoylation of ATP binding cassette transporter G1: Effect on protein trafficking and function. Biochim. Biophys. Acta 2013, 1831, 1067–1078. [Google Scholar] [CrossRef]
  65. Gu, H.M.; Wang, F.Q.; Zhang, D.W. Caveolin-1 interacts with ATP binding cassette transporter G1 (ABCG1) and regulates ABCG1-mediated cholesterol efflux. Biochim. Biophys. Acta 2014, 1841, 847–858. [Google Scholar] [CrossRef]
  66. Riggi, N.; Suva, M.L.; Stamenkovic, I. The cancer stem cell paradigm in Ewing’s sarcoma: What can we learn about these rare cells from a rare tumor? Expert Rev. Anticancer Ther. 2011, 11, 143–145. [Google Scholar] [CrossRef]
  67. Glumac, P.M.; LeBeau, A.M. The role of CD133 in cancer: A concise review. Clin. Transl. Med. 2018, 7, 18. [Google Scholar] [CrossRef]
  68. Rappa, G.; Fodstad, O.; Lorico, A. The stem cell-associated antigen CD133 (Prominin-1) is a molecular therapeutic target for metastatic melanoma. Stem Cells 2008, 26, 3008–3017. [Google Scholar] [CrossRef] [Green Version]
  69. Danieau, G.; Morice, S.; Redini, F.; Verrecchia, F.; Royer, B.B. New Insights about the Wnt/beta-Catenin Signaling Pathway in Primary Bone Tumors and Their Microenvironment: A Promising Target to Develop Therapeutic Strategies? Int. J. Mol. Sci. 2019, 20, 3751. [Google Scholar] [CrossRef] [Green Version]
  70. Tirode, F.; Laud-Duval, K.; Prieur, A.; Delorme, B.; Charbord, P.; Delattre, O. Mesenchymal stem cell features of Ewing tumors. Cancer Cell 2007, 11, 421–429. [Google Scholar] [CrossRef] [Green Version]
  71. von Levetzow, C.; Jiang, X.; Gwye, Y.; von Levetzow, G.; Hung, L.; Cooper, A.; Hsu, J.H.; Lawlor, E.R. Modeling initiation of Ewing sarcoma in human neural crest cells. PLoS ONE 2011, 6, e19305. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  72. Jawad, M.U.; Cheung, M.C.; Min, E.S.; Schneiderbauer, M.M.; Koniaris, L.G.; Scully, S.P. Ewing sarcoma demonstrates racial disparities in incidence-related and sex-related differences in outcome: An analysis of 1631 cases from the SEER database, 1973-2005. Cancer 2009, 115, 3526–3536. [Google Scholar] [CrossRef] [PubMed]
  73. Machiela, M.J.; Grunewald, T.G.P.; Surdez, D.; Reynaud, S.; Mirabeau, O.; Karlins, E.; Rubio, R.A.; Zaidi, S.; Grossetete-Lalami, S.; Ballet, S.; et al. Genome-wide association study identifies multiple new loci associated with Ewing sarcoma susceptibility. Nat. Commun. 2018, 9, 3184. [Google Scholar] [CrossRef] [PubMed]
  74. Randall, R.L.; Lessnick, S.L.; Jones, K.B.; Gouw, L.G.; Cummings, J.E.; Cannon-Albright, L.; Schiffman, J.D. Is There a Predisposition Gene for Ewing’s Sarcoma? J. Oncol. 2010, 2010, 397632. [Google Scholar] [CrossRef] [Green Version]
  75. Grunewald, T.G.P.; Cidre-Aranaz, F.; Surdez, D.; Tomazou, E.M.; de Alava, E.; Kovar, H.; Sorensen, P.H.; Delattre, O.; Dirksen, U. Ewing sarcoma. Nat. Rev. Dis. Prim. 2018, 4, 5. [Google Scholar] [CrossRef]
  76. Kumar, S.V.; Er, P.X.; Lawlor, K.T.; Motazedian, A.; Scurr, M.; Ghobrial, I.; Combes, A.N.; Zappia, L.; Oshlack, A.; Stanley, E.G.; et al. Kidney micro-organoids in suspension culture as a scalable source of human pluripotent stem cell-derived kidney cells. Development 2019, 146, dev172361. [Google Scholar] [CrossRef] [Green Version]
  77. Keller, F.; Rudiger, R.; MAthias, H. Towards optimized breast cancer 3D spheroid mono- and co-culture models for pharmacological research and screening. J. Cell. Biotechnol. 2019, 5, 89–101. [Google Scholar] [CrossRef] [Green Version]
  78. Wang, Y.J.; Herlyn, M. The emerging roles of Oct4 in tumor-initiating cells. Am. J. Physiol. Cell Physiol. 2015, 309, C709–C718. [Google Scholar] [CrossRef]
  79. Kim, J.; Kim, H.S.; Shim, J.J.; Lee, J.; Kim, A.Y.; Kim, J. Critical role of the fibroblast growth factor signalling pathway in Ewing’s sarcoma octamer-binding transcription factor 4-mediated cell proliferation and tumorigenesis. FEBS J. 2019, 286, 4443–4472. [Google Scholar] [CrossRef]
  80. Chua, V.; Chawla, S.P.; Zheng, K.; Kim, T.; Del Priore, G.; Kim, S. Phase II study of SM-88 in Ewing’s and other sarcomas. In Proceedings of the 2021 ASCO Annual Meeting I, Chicago, IL, USA, 28 May 2021. [Google Scholar]
  81. Chung, V.; Chawla, S.P.; Dong, H.; Kim, S.; Korn, R.L.; Lim, K.-H.; Noel, M.S.; Noonan, A.M.; Oberstein, P.E.; Ocean, A.J.; et al. Phase II/III study of SM-88 in patients with metastatic pancreatic cancer. In Proceedings of the 2021 Gastrointestinal Cancers Symposium, Online, 20 January 2021. [Google Scholar]
  82. Brennan, B.; Kirton, L.; Marec-Berard, P.; Martin-Broto, J.; Gelderblom, H.; Gaspar, N.; Strauss, S.; Sastre Urgelles, A.; Anderton, J.; Laurence, V.; et al. Comparison of two chemotherapy regimens in Ewing sarcoma (ES): Overall and subgroup results of the Euro Ewing 2012 randomized trial (EE2012). In Proceedings of the 2020 ASCO Annual Meeting I, Virtual, 25 May 2020; p. 11500. [Google Scholar]
  83. Serra, M.; Hattinger, C.M.; Pasello, M.; Casotti, C.; Fantoni, L.; Riganti, C.; Manara, M.C. Impact of ABC Transporters in Osteosarcoma and Ewing’s Sarcoma: Which Are Involved in Chemoresistance and Which Are Not? Cells 2021, 10, 2461. [Google Scholar] [CrossRef]
  84. Perri, T.; Fogel, M.; Mor, S.; Horev, G.; Meller, I.; Loven, D.; Issakov, J.; Kollender, Y.; Smirnov, A.; Zaizov, R.; et al. Effect of P-glycoprotein expression on outcome in the Ewing family of tumors. Pediatr. Hematol. Oncol. 2001, 18, 325–334. [Google Scholar] [CrossRef]
  85. Scotlandi, K.; Manara, M.C.; Hattinger, C.M.; Benini, S.; Perdichizzi, S.; Pasello, M.; Bacci, G.; Zanella, L.; Bertoni, F.; Picci, P.; et al. Prognostic and therapeutic relevance of HER2 expression in osteosarcoma and Ewing’s sarcoma. Eur. J. Cancer 2005, 41, 1349–1361. [Google Scholar] [CrossRef]
  86. Roessner, A.; Ueda, Y.; Bockhorn-Dworniczak, B.; Blasius, S.; Peters, A.; Wuisman, P.; Ritter, J.; Paulussen, M.; Jurgens, H.; Bocker, W. Prognostic implication of immunodetection of P glycoprotein in Ewing’s sarcoma. J. Cancer Res. Clin. Oncol. 1993, 119, 185–189. [Google Scholar] [CrossRef]
  87. Mancarella, C.; Pasello, M.; Ventura, S.; Grilli, A.; Calzolari, L.; Toracchio, L.; Lollini, P.L.; Donati, D.M.; Picci, P.; Ferrari, S.; et al. Insulin-Like Growth Factor 2 mRNA-Binding Protein 3 is a Novel Post-Transcriptional Regulator of Ewing Sarcoma Malignancy. Clin. Cancer Res. 2018, 24, 3704–3716. [Google Scholar] [CrossRef] [Green Version]
  88. Eisenmann, E.D.; Talebi, Z.; Sparreboom, A.; Baker, S.D. Boosting the oral bioavailability of anticancer drugs through intentional drug-drug interactions. Basic Clin. Pharmacol. Toxicol. 2022, 130 (Suppl. S1), 23–35. [Google Scholar] [CrossRef]
  89. Kang, Y. Analysis of cancer stem cell metastasis in xenograft animal models. Methods Mol. Biol. 2009, 568, 7–19. [Google Scholar] [CrossRef]
  90. Sampieri, K.; Fodde, R. Cancer stem cells and metastasis. Semin. Cancer Biol. 2012, 22, 187–193. [Google Scholar] [CrossRef]
  91. O’Callaghan, N.; Dhillon, V.; Thomas, P.; Fenech, M. A quantitative real-time PCR method for absolute telomere length. Biotechniques 2008, 44, 807–809. [Google Scholar] [CrossRef]
Figure 1. Self-renewing ability of ES cell lines. Single cell suspensions of ES cells were prepared in 0.3% soft agar or 1–10,000 cells seeded into each well of an ultra-low attachment plate (96 wells). Images show colonies formed at 21 days and are representative of three independent experiments.
Figure 1. Self-renewing ability of ES cell lines. Single cell suspensions of ES cells were prepared in 0.3% soft agar or 1–10,000 cells seeded into each well of an ultra-low attachment plate (96 wells). Images show colonies formed at 21 days and are representative of three independent experiments.
Cancers 15 00769 g001
Figure 2. Characterisation of CD133 cells in ES cell lines. (A) Protein expression of CD133 in ES cell lines (n = 6) was evaluated by Western blot. Equal protein loading was confirmed by probing the blots for α-tubulin. The top differentially expressed genes comparing (B) TC-32 and (C) A673 CD133-positive and -negative cells using LIMMA were not significant. The difference in ΔCt between the two groups and corresponding Q values are shown. Genes were considered significantly differentially expressed if the Q value was <0.1 and the difference in ∆Ct > 2. Genes with Ct values of >35 in both cell populations were excluded. Dark grey = ABC transporter proteins, grey = pluripotency associated genes, white = Wnt signalling pathway genes. Results of two independent experiments. SHEF-4 embryonic stem cells were included as a positive control for CD133 expression throughout.
Figure 2. Characterisation of CD133 cells in ES cell lines. (A) Protein expression of CD133 in ES cell lines (n = 6) was evaluated by Western blot. Equal protein loading was confirmed by probing the blots for α-tubulin. The top differentially expressed genes comparing (B) TC-32 and (C) A673 CD133-positive and -negative cells using LIMMA were not significant. The difference in ΔCt between the two groups and corresponding Q values are shown. Genes were considered significantly differentially expressed if the Q value was <0.1 and the difference in ∆Ct > 2. Genes with Ct values of >35 in both cell populations were excluded. Dark grey = ABC transporter proteins, grey = pluripotency associated genes, white = Wnt signalling pathway genes. Results of two independent experiments. SHEF-4 embryonic stem cells were included as a positive control for CD133 expression throughout.
Cancers 15 00769 g002aCancers 15 00769 g002b
Figure 3. Expression of pluripotency genes in SK-N-MC and TTC 466 cells in 2D and 3D spheroid cultures. (A) Summary of strategy to identify shared markers of stemness and MDR in ES cell lines. RNA (1 µg) from (B) SK-N-MC and (C) TTC 466 cells in 2D and 3D cultures were analysed by RTqPCR using the TaqMan® Human Stem Cell Pluripotency array. The level of mRNA expression is reported using the comparative Ct method, after normalisation of target Ct values to the global mean of all mRNAs. A volcano plot summarising the differentially expressed mRNAs in 2D and 3D spheroids, displaying the Q value (level of significance) and the difference in ΔCt between the cells in 2D and 3D spheroids is shown. The vertical dashed lines indicate the ΔCt thresholds of ±ΔCt of >2 (x-axis) and the horizontal line a significant Q value < 0.1 (y-axis). Black squares = mRNAs significantly decreased in SK-N-MC 3D spheroids and with a change in ΔCt > 2, Black diamonds = mRNAs significantly decreased in TTC 466 3D spheroids and with a change in ΔCt > 2, black triangles = mRNAs significantly increased in 3D spheroids and with a change in ΔCt >2, black circles = mRNAs with a ΔCt of ±<2 and Q value > 0.1. Genes with Ct values of >35 in 2D cells were excluded. The top differentially expressed genes with a Q value of <0.2 comparing (D) SK-N-MC and (E) TTC 466 cells in 2D and 3D cultures using LIMMA are shown. The difference in ΔCt between the two groups and corresponding Q values are shown. Genes were considered significantly differentially expressed if the Q value was <0.1 and the difference in ∆Ct > 2.
Figure 3. Expression of pluripotency genes in SK-N-MC and TTC 466 cells in 2D and 3D spheroid cultures. (A) Summary of strategy to identify shared markers of stemness and MDR in ES cell lines. RNA (1 µg) from (B) SK-N-MC and (C) TTC 466 cells in 2D and 3D cultures were analysed by RTqPCR using the TaqMan® Human Stem Cell Pluripotency array. The level of mRNA expression is reported using the comparative Ct method, after normalisation of target Ct values to the global mean of all mRNAs. A volcano plot summarising the differentially expressed mRNAs in 2D and 3D spheroids, displaying the Q value (level of significance) and the difference in ΔCt between the cells in 2D and 3D spheroids is shown. The vertical dashed lines indicate the ΔCt thresholds of ±ΔCt of >2 (x-axis) and the horizontal line a significant Q value < 0.1 (y-axis). Black squares = mRNAs significantly decreased in SK-N-MC 3D spheroids and with a change in ΔCt > 2, Black diamonds = mRNAs significantly decreased in TTC 466 3D spheroids and with a change in ΔCt > 2, black triangles = mRNAs significantly increased in 3D spheroids and with a change in ΔCt >2, black circles = mRNAs with a ΔCt of ±<2 and Q value > 0.1. Genes with Ct values of >35 in 2D cells were excluded. The top differentially expressed genes with a Q value of <0.2 comparing (D) SK-N-MC and (E) TTC 466 cells in 2D and 3D cultures using LIMMA are shown. The difference in ΔCt between the two groups and corresponding Q values are shown. Genes were considered significantly differentially expressed if the Q value was <0.1 and the difference in ∆Ct > 2.
Cancers 15 00769 g003aCancers 15 00769 g003b
Figure 4. Expression of ABC transporter proteins by SK-N-MC and TTC 466 cells in 2D and 3D spheroid cultures. RNA (1 µg) from (A) SK-N-MC and (D) TTC 466 cells in 2D and 3D cultures were analysed by RTqPCR using the TaqMan® Human ABC Transporter Array. The level of mRNA expression was reported using the comparative Ct method as mean ± SEM, after normalisation of target Ct values to the global mean of all mRNAs. Results are the mean of two independent experiments. * = Genes with change in ∆Ct > 2 and Q value < 0.1 comparing SK-N-MC cells in 2D and 3D cultures. Genes with Ct values of >35 in SK-N-MC cells grown in 2D were excluded. Volcano plot summarising the differentially expressed ABC transporter mRNAs in (B) SK-N-MC and (E) TTC 466 2D and 3D spheroids, displaying the Q value (level of significance) and the difference in ΔCt between the cells in 2D and 3D spheroids. The vertical dashed lines indicate the ΔCt thresholds of +/−ΔCt of >2 (x-axis) and the horizontal line a significant Q value < 0.1 (y-axis). Black square = CFTR which was significantly decreased and black triangle = ABCG1 which was significantly increased in 3D spheroids compared to 2D cultures. Black circles = ABC transporter proteins with a ΔCt of ±<2. The differentially expressed genes all have a Q value < 0.1. Mean Ct values of ABC transporter mRNAs in (C) SK-N-MC and (F) TTC 466 cells grown in 2D and 3D, the difference in ΔCt between the two groups and corresponding Q values are shown. Grey = genes with change in ∆Ct > 2 and Q value < 0.1. (G) ABC transporter and pluripotency mRNAs (n = 140) detected (Ct values < 35), highly expressed (Ct values < 25) and not detected (Ct values > 35) in 3D spheroids from TTC 466 and SK-N-MC cells. Venn diagrams show the number of mRNAs not detected, unique or shared for each cell line. (H) Validation of ABCG1 mRNA expression by RTqPCR and reported as 2−ΔΔCt in SK-N-MC cells grown in 2D and 3D culture; expression of ABCG1 is normalised to the endogenous control gene PPIA and the control cell line, SHEF-4. RNA expression was compared between 2D and 3D SK-N-MC cells, using a two-tailed t-test. The results are representative of 2 independent experiments. (I) Increased expression of ABCG1 protein in SK-N-MC cells grown as 3D spheroids compared to 2D cultures was validated by Western blot. Equal loading of each protein was confirmed by probing the Western blot for ß-actin. The results are representative of 2 independent experiments. (J) High expression of ABCG1 protein expression in the outer 50 µm region of SK-N-MC spheroids detected by IHC; nuclei are labelled with haematoxylin. Black scale bar = 100 µm. IgG control = SK-N-MC spheroid section incubated with the isotype control antibody (4 μg/mL, Negative Control Mouse IgG1, X0931 (Dako) and 20 μg/μL, Normal Rabbit Serum Control Ig mix, 086199 (Life Technologies), stained with haematoxylin.
Figure 4. Expression of ABC transporter proteins by SK-N-MC and TTC 466 cells in 2D and 3D spheroid cultures. RNA (1 µg) from (A) SK-N-MC and (D) TTC 466 cells in 2D and 3D cultures were analysed by RTqPCR using the TaqMan® Human ABC Transporter Array. The level of mRNA expression was reported using the comparative Ct method as mean ± SEM, after normalisation of target Ct values to the global mean of all mRNAs. Results are the mean of two independent experiments. * = Genes with change in ∆Ct > 2 and Q value < 0.1 comparing SK-N-MC cells in 2D and 3D cultures. Genes with Ct values of >35 in SK-N-MC cells grown in 2D were excluded. Volcano plot summarising the differentially expressed ABC transporter mRNAs in (B) SK-N-MC and (E) TTC 466 2D and 3D spheroids, displaying the Q value (level of significance) and the difference in ΔCt between the cells in 2D and 3D spheroids. The vertical dashed lines indicate the ΔCt thresholds of +/−ΔCt of >2 (x-axis) and the horizontal line a significant Q value < 0.1 (y-axis). Black square = CFTR which was significantly decreased and black triangle = ABCG1 which was significantly increased in 3D spheroids compared to 2D cultures. Black circles = ABC transporter proteins with a ΔCt of ±<2. The differentially expressed genes all have a Q value < 0.1. Mean Ct values of ABC transporter mRNAs in (C) SK-N-MC and (F) TTC 466 cells grown in 2D and 3D, the difference in ΔCt between the two groups and corresponding Q values are shown. Grey = genes with change in ∆Ct > 2 and Q value < 0.1. (G) ABC transporter and pluripotency mRNAs (n = 140) detected (Ct values < 35), highly expressed (Ct values < 25) and not detected (Ct values > 35) in 3D spheroids from TTC 466 and SK-N-MC cells. Venn diagrams show the number of mRNAs not detected, unique or shared for each cell line. (H) Validation of ABCG1 mRNA expression by RTqPCR and reported as 2−ΔΔCt in SK-N-MC cells grown in 2D and 3D culture; expression of ABCG1 is normalised to the endogenous control gene PPIA and the control cell line, SHEF-4. RNA expression was compared between 2D and 3D SK-N-MC cells, using a two-tailed t-test. The results are representative of 2 independent experiments. (I) Increased expression of ABCG1 protein in SK-N-MC cells grown as 3D spheroids compared to 2D cultures was validated by Western blot. Equal loading of each protein was confirmed by probing the Western blot for ß-actin. The results are representative of 2 independent experiments. (J) High expression of ABCG1 protein expression in the outer 50 µm region of SK-N-MC spheroids detected by IHC; nuclei are labelled with haematoxylin. Black scale bar = 100 µm. IgG control = SK-N-MC spheroid section incubated with the isotype control antibody (4 μg/mL, Negative Control Mouse IgG1, X0931 (Dako) and 20 μg/μL, Normal Rabbit Serum Control Ig mix, 086199 (Life Technologies), stained with haematoxylin.
Cancers 15 00769 g004aCancers 15 00769 g004b
Figure 5. Functional role and characterisation of ABCG1 in ES. (A) Representative dot plot of annexin V and PI labelling of SK-N-MC cells analysed by flow cytometry following knockdown of ABCG1 using SMARTpool:Accell ABCG1 siRNA and non-targeting control siRNAs. The upper left quadrant shows the mean percentage of necrotic cells, the upper right = apoptotic cells, lower left = viable cells and the lower right early apoptotic cells. The percentage of apoptotic cells was compared using a non-parametric Mann–Whitney two-tailed t-test. (B) Canonical ABCG1 RNA transcripts were downloaded from Ensembl.gov and labelled with the prefix ENST and unique transcript number. * = No protein produced from this transcript (www.ensembl.org, accessed on 27 April 2016), # = no protein produced from this transcript, retained intron (www.ensembl.org, accessed on 27 April 2016) (C) Canonical transcripts 1 to 11 were detected in SK-N-MC cells in 2D and 3D spheroids. In addition, two novel transcripts (transcript 1 and 2) were identified in 3D spheroids. Orange box = exon, grey arrows = direction of transcription, dotted line box = missing exon 5. (D) Novel transcript 1 was the most highly differentially expressed ABCG1 transcript in 3D spheroids compared to cells grown in 2D.
Figure 5. Functional role and characterisation of ABCG1 in ES. (A) Representative dot plot of annexin V and PI labelling of SK-N-MC cells analysed by flow cytometry following knockdown of ABCG1 using SMARTpool:Accell ABCG1 siRNA and non-targeting control siRNAs. The upper left quadrant shows the mean percentage of necrotic cells, the upper right = apoptotic cells, lower left = viable cells and the lower right early apoptotic cells. The percentage of apoptotic cells was compared using a non-parametric Mann–Whitney two-tailed t-test. (B) Canonical ABCG1 RNA transcripts were downloaded from Ensembl.gov and labelled with the prefix ENST and unique transcript number. * = No protein produced from this transcript (www.ensembl.org, accessed on 27 April 2016), # = no protein produced from this transcript, retained intron (www.ensembl.org, accessed on 27 April 2016) (C) Canonical transcripts 1 to 11 were detected in SK-N-MC cells in 2D and 3D spheroids. In addition, two novel transcripts (transcript 1 and 2) were identified in 3D spheroids. Orange box = exon, grey arrows = direction of transcription, dotted line box = missing exon 5. (D) Novel transcript 1 was the most highly differentially expressed ABCG1 transcript in 3D spheroids compared to cells grown in 2D.
Cancers 15 00769 g005
Figure 6. Expression profile of ABCG1, CD133 and MRP1 in ES cell lines and patient-derived cells. (A) The percentage of SK-N-MC cells grown in 2D and 3D cultures expressing ABCG1 and CD133 was quantified by flow cytometry; 10,000 events were examined for each condition. The percentage of positive cells is presented as the mean ± SEM (n = 6). ABCG1 and CD133 expression in SK-N-MC cells grown in 2D and 3D cultures was compared using a non-parametric Mann–Whitney two-tailed t-test. (B) Representative dot plot of ABCG1 and CD133 expression of SK-N-MC cells grown in 2D and disaggregated 3D cultures, analysed by flow cytometry. Quadrants represent cells with ABCG1 low CD133 high (upper left), ABCG1 high CD133 high (upper right), ABCG1 low CD133 low (lower left) and ABCG1 high CD133 low (lower right) levels of expression. Red spots = cells in 3D culture, blue spots = cells in 2D culture. (C) Protein expression of ABCG1 and CD133 in SK-N-MC cells grown in 2D and 3D culture was quantified by flow cytometry. Data are presented as the fold change in fluorescence, expressing the median target protein fluorescence in SK-N-MC cells grown in 3D relative to cells grown in 2D. Results show the mean ± SEM for two independent experiments, 3 replicates per experiment (n = 6). ABCG1 and CD133 expression in 2D and 3D cultures was compared using an unpaired two-tailed t-test. (D) ABCG1 protein expression in CD133-positive and CD133-negative TC-32 cells detected by Western blot. Equal protein loading was confirmed by probing the blots for ß-actin. Results are representative of 2 independent experiments. (E) ABCG1 protein expression in patient-derived ES cells, determined by Western blot. Equal protein loading was confirmed by probing the blots for ß-actin. Results are representative of 2 independent experiments. (F) There was no correlation between the level of CD133 expression and ability to produce progeny from a single cell. The level of CD133 expression was quantified by flow cytometry and the ability to produce progeny from a single cell evaluated in 2D adherent culture at 21 days. Correlation was examined using linear regression. Open circles = cell lines, filled circles = patient-derived cells. (G) Summary of ABCG1 and MRP1 protein expression in ES cells grown in 2D and 3D culture. Salmon-pink circle = ES cells grown in 2D, orange circle = ES cells grown in 3D (spheroids), blue circle = ABCG1 protein, green circle = CD133 protein, black line = plasma membrane.
Figure 6. Expression profile of ABCG1, CD133 and MRP1 in ES cell lines and patient-derived cells. (A) The percentage of SK-N-MC cells grown in 2D and 3D cultures expressing ABCG1 and CD133 was quantified by flow cytometry; 10,000 events were examined for each condition. The percentage of positive cells is presented as the mean ± SEM (n = 6). ABCG1 and CD133 expression in SK-N-MC cells grown in 2D and 3D cultures was compared using a non-parametric Mann–Whitney two-tailed t-test. (B) Representative dot plot of ABCG1 and CD133 expression of SK-N-MC cells grown in 2D and disaggregated 3D cultures, analysed by flow cytometry. Quadrants represent cells with ABCG1 low CD133 high (upper left), ABCG1 high CD133 high (upper right), ABCG1 low CD133 low (lower left) and ABCG1 high CD133 low (lower right) levels of expression. Red spots = cells in 3D culture, blue spots = cells in 2D culture. (C) Protein expression of ABCG1 and CD133 in SK-N-MC cells grown in 2D and 3D culture was quantified by flow cytometry. Data are presented as the fold change in fluorescence, expressing the median target protein fluorescence in SK-N-MC cells grown in 3D relative to cells grown in 2D. Results show the mean ± SEM for two independent experiments, 3 replicates per experiment (n = 6). ABCG1 and CD133 expression in 2D and 3D cultures was compared using an unpaired two-tailed t-test. (D) ABCG1 protein expression in CD133-positive and CD133-negative TC-32 cells detected by Western blot. Equal protein loading was confirmed by probing the blots for ß-actin. Results are representative of 2 independent experiments. (E) ABCG1 protein expression in patient-derived ES cells, determined by Western blot. Equal protein loading was confirmed by probing the blots for ß-actin. Results are representative of 2 independent experiments. (F) There was no correlation between the level of CD133 expression and ability to produce progeny from a single cell. The level of CD133 expression was quantified by flow cytometry and the ability to produce progeny from a single cell evaluated in 2D adherent culture at 21 days. Correlation was examined using linear regression. Open circles = cell lines, filled circles = patient-derived cells. (G) Summary of ABCG1 and MRP1 protein expression in ES cells grown in 2D and 3D culture. Salmon-pink circle = ES cells grown in 2D, orange circle = ES cells grown in 3D (spheroids), blue circle = ABCG1 protein, green circle = CD133 protein, black line = plasma membrane.
Cancers 15 00769 g006aCancers 15 00769 g006b
Figure 7. Identification of candidate shared stemness and chemoresistance therapeutic targets. (A) Summary of strategy to identify shared markers of stemness and MDR in ES preclinical models. (B) Stemness and MDR-associated genes increased in 3D spheroids and patient-derived ES-CSCs. The full and alternative gene names, the Ensembl gene abbreviation (Gene) and the data from which genes were identified (Source of the gene) are shown. The source of the gene is presented as 1 or 2, where 1 = target gene identified from transcriptome analysis of patient-derived ES-CSCs [17] and 2 = target gene identified by comparing transcriptome of SK-N-MC and TTC 466 cells grown in 3D and 2D (Figure 3 and Figure 4). (C) Pipeline to identify candidate drugs based on 21 stemness and chemoresistance genes. Green ovals = drug-repurposing data sources, pink ovals = public data sources, blue squares = number of targets or drugs at each step. (D) Predicted interaction between the 13 shared stemness and chemoresistance genes (excluding p-glycoprotein and MRP1) that have corresponding candidate drugs using STRING [30]. Five of the genes are known to bind and regulate vital biological processes, and could be candidate prognostic biomarkers and/or therapeutic targets in ES. Pink lines = known interaction that has been experimentally determined, extracted from BIND, DIP, GRID, HPRD, IntAct, MINT and PID databases. Blue lines = known interaction based on data from curated databases Biocarta, BioCyc, GO, KEGG and Reactome. Grey lines = proteins reported to be co-expressed. Green lines = interactions identified by text mining. Filled nodes = 3D structure known or predicted. Coloured nodes = first shell of interacting proteins. White nodes = second shell of interacting proteins. Parental cells = patient-derived bulk population of ES cells from which the ES-CSCs were derived [14].
Figure 7. Identification of candidate shared stemness and chemoresistance therapeutic targets. (A) Summary of strategy to identify shared markers of stemness and MDR in ES preclinical models. (B) Stemness and MDR-associated genes increased in 3D spheroids and patient-derived ES-CSCs. The full and alternative gene names, the Ensembl gene abbreviation (Gene) and the data from which genes were identified (Source of the gene) are shown. The source of the gene is presented as 1 or 2, where 1 = target gene identified from transcriptome analysis of patient-derived ES-CSCs [17] and 2 = target gene identified by comparing transcriptome of SK-N-MC and TTC 466 cells grown in 3D and 2D (Figure 3 and Figure 4). (C) Pipeline to identify candidate drugs based on 21 stemness and chemoresistance genes. Green ovals = drug-repurposing data sources, pink ovals = public data sources, blue squares = number of targets or drugs at each step. (D) Predicted interaction between the 13 shared stemness and chemoresistance genes (excluding p-glycoprotein and MRP1) that have corresponding candidate drugs using STRING [30]. Five of the genes are known to bind and regulate vital biological processes, and could be candidate prognostic biomarkers and/or therapeutic targets in ES. Pink lines = known interaction that has been experimentally determined, extracted from BIND, DIP, GRID, HPRD, IntAct, MINT and PID databases. Blue lines = known interaction based on data from curated databases Biocarta, BioCyc, GO, KEGG and Reactome. Grey lines = proteins reported to be co-expressed. Green lines = interactions identified by text mining. Filled nodes = 3D structure known or predicted. Coloured nodes = first shell of interacting proteins. White nodes = second shell of interacting proteins. Parental cells = patient-derived bulk population of ES cells from which the ES-CSCs were derived [14].
Cancers 15 00769 g007aCancers 15 00769 g007bCancers 15 00769 g007c
Table 1. Self-renewing ability of ES cell lines. A single cell suspension of each cell line was seeded into soft agar or into a 96 well ultra-low adherence plate at a density of 1–10,000 cells. The number of colonies was counted after 21 days. For studies in soft agar, the mean (±SEM) percentage colony formation is expressed as the number of colonies counted at 21 days relative to the number of cells seeded. For cells seeded onto an ultra-low adherence plate, the mean (±SEM) percentage colony formation is expressed as the number of wells containing a colony of >5 cells relative to the total number of wells with cells seeded. Results are representative of three independent experiments. Percentage colony formation was compared using ANOVA and Bonferroni’s post hoc test.
Table 1. Self-renewing ability of ES cell lines. A single cell suspension of each cell line was seeded into soft agar or into a 96 well ultra-low adherence plate at a density of 1–10,000 cells. The number of colonies was counted after 21 days. For studies in soft agar, the mean (±SEM) percentage colony formation is expressed as the number of colonies counted at 21 days relative to the number of cells seeded. For cells seeded onto an ultra-low adherence plate, the mean (±SEM) percentage colony formation is expressed as the number of wells containing a colony of >5 cells relative to the total number of wells with cells seeded. Results are representative of three independent experiments. Percentage colony formation was compared using ANOVA and Bonferroni’s post hoc test.
Cell LineClone or Spheroid Formation (%)
Clone Formation in Soft Agar from a Single CellUltra-Low Attachment Plate
1 Cell10 Cells100 Cells1000 Cells10,000 Cells
A67317 ±118 ± 588 ± 10100100100
RD-ES17 ±15 ± 145 ± 13100100100
SKES-129 ± 2 (p < 0.05)75 ± 7 (p < 0.05)98 ± 2100100100
SK-N-MC12 ± 164 ± 197 ± 3100100100
TC-327 ± 0.5 (p < 0.05)15 ± 393 ± 9100100100
TTC 46613 ± 145 ± 296 ± 6100100100
Table 2. Drugs directed against the shared stemness and chemoresistance therapeutic targets. Drugs directed against the shared ES stemness and chemoresistance targets that have been used in human clinical trials, identified using our in-house pipeline. Drug names, direct molecular targets, whether the drug has been previously used in the treatment of cancer patients and is on- or off-patent are shown. The number, trial phase and status, cancer type and detailed trial information of clinical trials are listed. NA = no phase 3 or phase 4 trials recorded for the drug. BNFC = British National Formulary- Children, a directory or all drugs approved for use in children. Y = yes, N = no. * = direct ES targets of fostamatinib (KIT, STK10, SRC, SIK1, RPS6KA1, ROCK2, RET, PTK2B, PTK2, PRKG2, PRKCD, PLK4, PLK3, PLK1, PIM3, PIK3CG, PIK3CD, PDGFRB, PDGFRA, PAK3, PAK1, NTRK3, NTRK2, NTRK1, NEK2, MST1R, MAPK14, MAP3K3, MAP2K2, LYN, LIMK1, KDR, ITK, INSRR, INSR, GSK3A, MTOR, FLT4, FLT3, FLT1, FGFR1, ERN1, ERBB4, ERBB2, EPHB4, EPHB2, EPHA3, EPHA2, EIF2AK4, EGFR, DYRK1B, DCLK3, DCLK1, CSK, CSF1R, CHEK2, CHEK1, CDK4, BTK, BRAF, AXL, AURKB, AURKA, ALK, ZAP70, WEE1, TNK2, TNK1, TIE1, TGFBR2, TEK, TAOK3, STK36, STK33, MET, JAK1, ABL1, SYK).
Table 2. Drugs directed against the shared stemness and chemoresistance therapeutic targets. Drugs directed against the shared ES stemness and chemoresistance targets that have been used in human clinical trials, identified using our in-house pipeline. Drug names, direct molecular targets, whether the drug has been previously used in the treatment of cancer patients and is on- or off-patent are shown. The number, trial phase and status, cancer type and detailed trial information of clinical trials are listed. NA = no phase 3 or phase 4 trials recorded for the drug. BNFC = British National Formulary- Children, a directory or all drugs approved for use in children. Y = yes, N = no. * = direct ES targets of fostamatinib (KIT, STK10, SRC, SIK1, RPS6KA1, ROCK2, RET, PTK2B, PTK2, PRKG2, PRKCD, PLK4, PLK3, PLK1, PIM3, PIK3CG, PIK3CD, PDGFRB, PDGFRA, PAK3, PAK1, NTRK3, NTRK2, NTRK1, NEK2, MST1R, MAPK14, MAP3K3, MAP2K2, LYN, LIMK1, KDR, ITK, INSRR, INSR, GSK3A, MTOR, FLT4, FLT3, FLT1, FGFR1, ERN1, ERBB4, ERBB2, EPHB4, EPHB2, EPHA3, EPHA2, EIF2AK4, EGFR, DYRK1B, DCLK3, DCLK1, CSK, CSF1R, CHEK2, CHEK1, CDK4, BTK, BRAF, AXL, AURKB, AURKA, ALK, ZAP70, WEE1, TNK2, TNK1, TIE1, TGFBR2, TEK, TAOK3, STK36, STK33, MET, JAK1, ABL1, SYK).
DrugEwings Sarcoma TargetsCancer DrugBNFC approvalOff-PatentTrial CountMax Trial PhaseTrial Phase (Number of Trials in Phase)Trial Status (Number of Trials in Status)Cancer Type (Number of Trials in Cancer Type)Trial Information for Phase 3 and Phase 4 trials (Trial Number, Phase, Status, Recruitment Status, Cancer Type)
CabozantinibKIT, RET, KDR, METYNN6Phase 2Phase 1 (1), Phase 2 (5)Recruiting (1), Active, not recruiting (3), Not yet recruiting (1), Other (1)Multiple sarcomas, including Ewing’s (1), Mixed paediatric, including Ewing’s (3), Ewing’s or osteosarcoma (2)NA
DasatinibKIT, MAPK14, LYN, HSPA8, EPHB4, CSK, BCR, BTK, PDGFRB, EPHA2, SRC, ABL1YYY2Phase 2Phase 1/2 (1), Phase 2 (1)Active, not recruiting (1), Completed (1)Mixed paediatric, including Ewing’s (1), Multiple sarcomas, including Ewing’s (1)NA
DecitabineDNMT3B, DNMT1YNY1Phase 1Phase 1 (1)Completed (1)Mixed, including sarcomas (1)NA
ErdafitinibKIT, KDR, PDGFRB, PDGFRA, CSF1R, RET, FGFR1YNN1Phase 2Phase 2 (1)Recruiting (1)Mixed paediatric, including Ewing’s (1)NA
ImatinibKIT, ABCB1, PDGFRB, ABL1, PDGFRA, CSF1R, NTRK1, RET, BCRYYY4Phase 2Phase 2 (4)Completed (4)Mixed, including sarcomas (1), Multiple sarcomas, including Ewing’s (1), Mixed paediatric, including Ewing’s (1), Ewing’s or DSRCT (1)NI
LenvatinibKIT, RET, PDGFRA, FGFR1, FLT4, KDR, FLT1YNN1Phase 1/2Phase 1/2 (1)Active, not recruiting (1)Mixed paediatric, including Ewing’s (1)NA
PazopanibKIT, ITK, PDGFRB, PDGFRA, FLT4, KDR, FLT1YNY1Phase 2Phase 2 (1)Completed (1)Mixed paediatric, including Ewing’s (1)NA
RegorafenibKIT, RET, ABL1, BRAF, EPHA2, NTRK1, TEK, FGFR1, PDGFRB, PDGFRA, FLT4, KDR, FLT1YNN5Phase 2Phase 1 (1), Phase 1/2 (1), Phase 2 (2), Other (1)Recruiting (4), Active, not recruiting (1)Multiple sarcomas, including Ewing’s (3), Ewing’s or osteosarcoma (1), Mixed paediatric, including Ewing’s (1)NA
SorafenibKIT, ABCB1, FLT1, RET, FGFR1, PDGFRB, FLT3, KDR, FLT4, BRAFYNUnclear2Phase 2Phase 1 (1), Phase 2 (1)Active, not recruiting (1), Completed (1)Ewing’s or DSRCT (1), Mixed paediatric, including Ewing’s (1)NA
SunitinibKIT, ABCB1, PDGFRA, CSF1R, FLT3, FLT4, KDR, FLT1, PDGFRBYNY2Phase 2Phase 1/2 (1), Phase 2 (1)Ongoing (1), Completed (1)Multiple sarcomas, including Ewing’s (1), Mixed, including sarcomas (1)NA
TemozolomideABCB1YYY23Phase 3Phase 1 (10), Phase 1/2 (4), Phase 2 (8), Phase 3 (1)Recruiting (8), Active, not recruiting (4), Completed (6), Not yet recruiting (2), Other (2), Terminated (1)Mixed paediatric, including Ewing’s (6), Ewing’s (10), Mixed, including sarcomas (4), Ewing’s or DSRCT (1), Ewing’s or RMS (2)NCT03495921 (Phase 3; Active, not recruiting)—Ewing sarcoma
AliskirenABCB1NNY1Phase 2Phase 2 (1)Recruiting (1)Multiple cancer types (1)NA
AllopurinolPOU5F1NYY1Phase 1/2Phase 1/2 (1)Not yet recruiting (1)Lung (1)NA
AtorvastatinABCB1, NR1I3, HDAC2, AHR, DPP4, HMGCRNYY16Phase 3Phase 1 (4), Phase 2 (6), Phase 2/3 (1), Phase 3 (5)Recruiting (13), Ongoing (1), Not yet recruiting (2)Multiple cancer types, Leukaemia (1), GI (4), Breast (8), Urological (3)NCT03024684 (Phase 3; Recruiting)—Hepatocellular carcinoma, NCT03819101 (Phase 3; Recruiting)—Prostate, NCT03971019 (Phase 3; Recruiting)—Breast, NCT04601116 (Phase 3; Recruiting)—Breast, NCT04026230 (Phase 3; Recruiting)—Prostate
AzithromycinABCC1NYY3Phase 2Phase 2 (3)Ongoing (3)Lymphoma (2), Breast (1)NA
CarvedilolABCB1, HIF1A, GJA1, VEGFANYY2Phase 2Phase 1 (1), Phase 2 (1)Active, not recruiting (1), Not yet recruiting (1)Urological (1), CNS (1)NA
ChlorpromazineABCB1NNY3Phase 2Phase 1 (1), Phase 1/2 (1), Phase 2 (1)Recruiting (2), Not yet recruiting (1)GI (1), CNS (2)NA
CitalopramABCB1NYY1Phase 3Phase 3 (1)Ongoing (1)CNS (1)2013-004705-59 (Phase 3; Ongoing)—Glioblastoma
ClarithromycinABCB1NYY18Phase 4Phase 1/2 (2), Phase 2 (10), Phase 3 (3), Phase 4 (1), Other (2)Recruiting (8), Active, not recruiting (8), Not yet recruiting (1), Other (1)Lymphoma (1), Other Haem-onc (15), Multiple cancer types (1), GI (1)ChiCTR2100047608 (Phase 4; Recruiting)—Multiple Myeloma, NCT02575144 (Phase 3; Active, not recruiting)—Multiple Myeloma, NCT02516696 (Phase 3; Active, not recruiting)—Multiple Myeloma, NCT04287660 (Phase 3; Recruiting)—Multiple Myeloma
ClopidogrelABCB1NNY1Phase 1Phase 1 (1)Recruiting (1)Head and Neck (1)NA
ColchicineABCB1, TUBBNNY2Phase 2Phase 1 (1), Phase 2 (1)Recruiting (2)Urological, Multiple cancer types (1), GI (1)NA
CyclosporineABCC1, ABCB1NNY2Phase 2Phase 1 (1), Phase 2 (1)Recruiting (1), Active, not recruiting (1)Leukaemia (1), Other Haem-onc (1)NA
DeferoxamineABCB1NNY4Phase 3Phase 1 (2), Phase 2 (1), Phase 3 (1)Recruiting (3), Not yet recruiting (1)Multiple cancer types (1), Breast (1), Leukaemia (1), GI (1)IRCT20200313046756N2 (Phase 3; Pending)—Acute Myeloid Leukaemia
DigoxinABCB1NYY2Phase 2Phase 1 (1), Phase 2 (1)Recruiting (2)GI (1), Multiple cancer types (1)NA
DisulfiramABCB1NNY13Phase 2/3Phase 1 (4), Phase 1/2 (2), Phase 2 (6), Phase 2/3 (1)Recruiting (7), Active, not recruiting (2), Ongoing (3), Not yet recruiting (1)Multiple cancer types, GI (1), Soft-Tissue Sarcoma, Bone Sarcoma (1), CNS (3), Other Haem-onc (1), Urological (2), Multiple cancer types, Breast (1), GI (2), Breast (2)NA
DoxycyclineABCB1NYY10Phase 2Phase 1 (1), Phase 2 (9)Recruiting (4), Active, not recruiting (4), Ongoing (1), Not yet recruiting (1)Urological (1), Lymphoma (3), Multiple cancer types (1), Breast, Gynaecological (1), Breast (1), GI (2), Head and Neck (1)NA
FenofibrateABCB1, NR1I2NYY1Phase 2Phase 2 (1)Recruiting (1)CNS (1)NA
Fostamatinib*NNN2Phase 1Phase 1 (2)Recruiting (2)Gynaecological (1), Other Haem-onc, Leukaemia (1)NA
IndomethacinABCC1, PTGS2NYY5Phase 4Phase 1 (1), Phase 1/2 (1), Phase 2 (1), Phase 3 (1), Phase 4 (1)Recruiting (2), Active, not recruiting (3)Urological (2), Breast (1), Head and Neck (2)ChiCTR2000038968 (Phase 4; Recruiting)—Prostate, NCT01265849 (Phase 3; Active, not recruiting)—Cancer of the oral cavity
ItraconazoleABCB1, ERBB2NYY15Phase 3Phase 1 (4), Phase 1/2 (2), Phase 2 (7), Phase 3 (1), Other (1)Recruiting (9), Active, not recruiting (2), Completed (1), Not yet recruiting (1), Other (1), Suspended (1)Leukaemia, Other Haem-onc (1), Lung (2), Multiple cancer types, Lymphoma, Leukaemia (1), follow-up continuing (1), Multiple cancer types (2), GI (4), Urological (1), Gynaecological (2), Skin (1)NCT03458221 (Phase 3; Not yet recruiting)—Ovarian
IvermectinABCB1NYY2Phase 2Phase 2 (2)Recruiting (1), Not yet recruiting (1)Breast (1), Multiple cancer types (1)NA
KetoconazoleABCB1, NR1I3, NR1I2, ARNYY4Phase 2Phase 1 (2), Phase 2 (2)Recruiting (1), Active, not recruiting (2), Ongoing (1)Urological (2), Breast, CNS (1), CNS (1)NA
LansoprazoleABCB1, MAPTNYY3Phase 3Phase 2 (1), Phase 3 (2)Recruiting (2), Active, not recruiting (1)Breast (1), Lymphoma (1), GI (1)NCT04874935 (Phase 3; Recruiting)—Breast, NCT03647072 (Phase 3; Recruiting)—Non-Hodgkin Lymphoma
LevetiracetamABCB1NYY2Phase 2Phase 2 (1), Other (1)Not yet recruiting (2)CNS (2)NA
LosartanABCB1NNY12Phase 3Phase 1 (4), Phase 2 (7), Phase 3 (1)Recruiting (9), Active, not recruiting (1), Not yet recruiting (2)Multiple cancer types (3), Bone Sarcoma (1), GI (7), Breast (1)CTRI/2021/05/033482 (Phase 3; Not Yet Recruiting)—Pancreatic
LovastatinABCB1, HDAC2, HMGCRNNY1OtherOther (1)Not yet recruiting (1)GI (1)NA
MaprotilineABCB1NNY1Phase 1Phase 1 (1)Not yet recruiting (1)CNS (1)NA
MefloquineABCB1NYY1Phase 1Phase 1 (1)Active, not recruiting (1)CNS (1)NA
MiconazoleABCB1, NR1I2NYY2Phase 2Phase 1 (1), Phase 2 (1)Recruiting (2)Multiple cancer types (2)NA
MidazolamGABRB3, ABCB1NYY2Phase 2Phase 2 (2)Recruiting (1), Not yet recruiting (1)Urological (2)NA
MifepristoneABCB1, NR1I2, KLK3, NR3C1, PGRNNY2Phase 3Phase 2 (1), Phase 3 (1)Active, not recruiting (1), Not yet recruiting (1)Breast (2)NCT05016349 (Phase 3; Not yet recruiting)—Breast
MiltefosineABCB1NNY2Phase 2Phase 2 (2)Recruiting (1), Active, not recruiting (1)Breast (1), Multiple cancer types (1)NA
NelfinavirABCB1NNY9Phase 3Phase 1 (2), Phase 1/2 (2), Phase 2 (3), Phase 3 (2)Recruiting (5), Active, not recruiting (2), Ongoing (2)Soft Tissue Sarcoma (1), Urological, Multiple cancer types (1), Gynaecological (3), Head and Neck (1), Other Haem-onc (1), GI (1), CNS (1)NCT03256916 (Phase 3; Recruiting)—Advanced carcinoma of the cervix, CTRI/2017/08/009265 (Phase 3; Open to Recruitment)—Advanced carcinoma of the cervix
NicardipineABCB1NNY1OtherOther (1)Not yet recruiting (1)CNS (1)NA
OmeprazoleABCB1, AHRNYY3Phase 2Phase 1 (2), Phase 2 (1)Recruiting (2), Active, not recruiting (1)Urological (1), Breast (1), GI (1)NA
PantoprazoleABCB1NNY2Phase 2Phase 1/2 (1), Phase 2 (1)Active, not recruiting (1), Not yet recruiting (1)Head and Neck (1), Urological (1)NA
PhenytoinPOU5F1, SCN8A, SCN2A, SCN3A, NR1I2, SCN1A, SCN5ANYY1Phase 2Phase 2 (1)Recruiting (1)Bone Sarcoma, Soft-Tissue Sarcoma (1)NA
PravastatinABCB1, HDAC2, HMGCRNNY2Phase 4Phase 2 (1), Phase 4 (1)Active, not recruiting (1), Not yet recruiting (1)Leukaemia (1), Breast (1)ChiCTR2000034035 (Phase 4; Pending)—Breast
PropofolGABRB3, ABCB1, SCN2A, SCN4ANYY29Phase 4Phase 2 (1), Phase 3 (2), Phase 4 (7), Other (19)Recruiting (17), Active, not recruiting (2), Ongoing (1), Not yet recruiting (9)GI (6), Breast, GI (1), Multiple cancer types (7), Breast (5), Urological (3), Lung (6), CNS (1)NCT01975064 (Phase 4; Recruiting)—Breast, colon, rectal, NCT05331911 (Phase 4; Not yet recruiting)—Liver, NCT04475705 (Phase 4; Recruiting)—Paediatric solid tumours, NCT05141877 (Phase 4; Not yet recruiting)—Brain tumour, NCT03034096 (Phase 4; Recruiting)—Adult cancer, NCT04513808 (Phase 3; Recruiting)—Oesophageal caner, ChiCTR2000040604 (Phase 4; Pending)—Non-small cell lung cancer, ACTRN12611000301965 (Phase 4; Not yet recruiting)—Breast, 2009-009114-40 (Phase 3; Ongoing)—Prostate
PropranololABCB1, EGFR, ADRB3NNY19Phase 3Phase 1 (3), Phase 1/2 (1), Phase 2 (13), Phase 2/3 (1), Phase 3 (1)Recruiting (10), Active, not recruiting (1), Not yet recruiting (7), Other (1)GI (8), Urological (2), Soft Tissue Sarcoma (2), Gynaecological (1), Skin (3), Other (1), Multiple cancer types (2)CTRI/2019/11/021924 (Phase 3; Not Yet Recruiting)—Ovarian
RifampicinABCB1, NR1I2NYY1Phase 1Phase 1 (1)Recruiting (1)Multiple cancer types, Lymphoma, Leukaemia (1)NA
RitonavirABCC1, ABCB1, NR1I2NYY1Phase 1Phase 1 (1)Recruiting (1)Leukaemia, Lymphoma (1)NA
SertralineABCB1, SLC29A4, SLC6A3NYY1Phase 1Phase 1 (1)Recruiting (1)Leukaemia (1)NA
SimvastatinABCB1, HDAC2, HMGCRNYY18Phase 4Phase 1 (4), Phase 2 (11), Phase 3 (1), Phase 4 (2)Recruiting (7), Active, not recruiting (3), Ongoing (3), Not yet recruiting (4), Suspended (1)Breast (5), Other Haem-onc (1), Lung (3), Gynaecological (1), GI (6), Multiple cancer types (2)ChiCTR2000034035 (Phase 4; Pending)—Breast, 2010-018491-24 (Phase 4; Ongoing)—Adults with bone metastasis, NCT03971019 (Phase 3; Recruiting)—Breast
SirolimusABCB1, MTORNYY28Phase 4Phase 1 (7), Phase 1/2 (4), Phase 2 (13), Phase 3 (2), Phase 4 (1), Other (1)Recruiting (17), Active, not recruiting (7), Ongoing (3), Not yet recruiting (1)Lung (2), Multiple cancer types, Lung (1), Multiple cancer types (6), Other (2), Soft-Tissue Sarcoma (3), Urological, Multiple cancer types (1), Bone Sarcoma, Soft Tissue Sarcoma (1), Gynaecological (1), Leukaemia (3), Endocrine (1), GI (3), CNS (1), Bone Sarcoma (2), Breast (1)NCT04775173 (Phase 3; Recruiting)—Haemangioendothelioma, ChiCTR1900021896 (Phase 4; Recruiting)—Liver, NCT04736589 (Phase 3; Not yet recruiting)—Breast
ValganciclovirABCB1NYY5Phase 2Phase 1/2 (2), Phase 2 (2), Other (1)Recruiting (4), Not yet recruiting (1)GI (1), Lymphoma (1), Head and Neck (1), CNS (1), Multiple cancer types (1)NA
VerapamilABCC1, ABCB1NNY1Phase 1Phase 1 (1)Active, not recruiting (1)Lymphoma (1)NA
WarfarinABCB1, AXL, NR1I2NNY1Phase 1Phase 1 (1)Recruiting (1)GI (1)NA
ZidovudineABCC1, ABCB1, TERTNYY1Phase 2Phase 2 (1)Recruiting (1)Leukaemia (1)NA
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Roundhill, E.A.; Pantziarka, P.; Liddle, D.E.; Shaw, L.A.; Albadrani, G.; Burchill, S.A. Exploiting the Stemness and Chemoresistance Transcriptome of Ewing Sarcoma to Identify Candidate Therapeutic Targets and Drug-Repurposing Candidates. Cancers 2023, 15, 769. https://doi.org/10.3390/cancers15030769

AMA Style

Roundhill EA, Pantziarka P, Liddle DE, Shaw LA, Albadrani G, Burchill SA. Exploiting the Stemness and Chemoresistance Transcriptome of Ewing Sarcoma to Identify Candidate Therapeutic Targets and Drug-Repurposing Candidates. Cancers. 2023; 15(3):769. https://doi.org/10.3390/cancers15030769

Chicago/Turabian Style

Roundhill, Elizabeth Ann, Pan Pantziarka, Danielle E. Liddle, Lucy A. Shaw, Ghadeer Albadrani, and Susan Ann Burchill. 2023. "Exploiting the Stemness and Chemoresistance Transcriptome of Ewing Sarcoma to Identify Candidate Therapeutic Targets and Drug-Repurposing Candidates" Cancers 15, no. 3: 769. https://doi.org/10.3390/cancers15030769

APA Style

Roundhill, E. A., Pantziarka, P., Liddle, D. E., Shaw, L. A., Albadrani, G., & Burchill, S. A. (2023). Exploiting the Stemness and Chemoresistance Transcriptome of Ewing Sarcoma to Identify Candidate Therapeutic Targets and Drug-Repurposing Candidates. Cancers, 15(3), 769. https://doi.org/10.3390/cancers15030769

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop