Repurposing Clemastine to Target Glioblastoma Cell Stemness

Simple Summary Brain tumor-initiating cells (BTICs) drive tumor progression and resistance to treatments, posing formidable challenges to advancing effective treatments against glioblastoma (GBM). We postulated that inducing BTIC differentiation can serve as a solution to diminishing their malignant features. In this study, we found that clemastine, an over-the-counter oral medication for allergy relief, attenuated the propagation and promoted the differentiation of BTICs, and we uncovered the indispensable role of EBP (Emopamil-binding protein) in maintaining the BTIC population. Taken together, our study implicates specific pathways in the perpetuation of BTICs, and identifies a non-oncology drug with a well-established safety profile that can be repurposed to mitigate the malignant properties of BTICs in GBM. Abstract Brain tumor-initiating cells (BTICs) and tumor cell plasticity promote glioblastoma (GBM) progression. Here, we demonstrate that clemastine, an over-the-counter drug for treating hay fever and allergy symptoms, effectively attenuated the stemness and suppressed the propagation of primary BTIC cultures bearing PDGFRA amplification. These effects on BTICs were accompanied by altered gene expression profiling indicative of their more differentiated states, resonating with the activity of clemastine in promoting the differentiation of normal oligodendrocyte progenitor cells (OPCs) into mature oligodendrocytes. Functional assays for pharmacological targets of clemastine revealed that the Emopamil Binding Protein (EBP), an enzyme in the cholesterol biosynthesis pathway, is essential for BTIC propagation and a target that mediates the suppressive effects of clemastine. Finally, we showed that a neural stem cell-derived mouse glioma model displaying predominantly proneural features was similarly susceptible to clemastine treatment. Collectively, these results identify pathways essential for maintaining the stemness and progenitor features of GBMs, uncover BTIC dependency on EBP, and suggest that non-oncology, low-toxicity drugs with OPC differentiation-promoting activity can be repurposed to target GBM stemness and aid in their treatment.


Introduction
Glioblastoma (GBM) tumors demonstrate striking aggressiveness and therapeutic resistance, which are driven by brain tumor-initiating cells (BTICs), a subpopulation of GBM cells that exhibit phenotypic plasticity and stem/progenitor-like features [1][2][3][4][5][6].The pathological ramifications of GBM plasticity and stemness are highlighted by the heterogeneous nature of GBM with the presence of various subtypes, including proneural (PN), classical (CL), and mesenchymal (MES) subtypes [7][8][9][10][11][12].Among the GBM subtypes, PN tumors are characterized by stem/progenitor-like signatures, a failure to respond to more aggressive chemo-and radio-therapy [9,13], and the ability to undergo the process of proneural-mesenchymal transition (PMT) [14,15].Facilitated and driven by GBM plasticity and stemness, these tumors can transition from one subtype to another in the course of tumor progression and recurrence, and give rise to the most aggressive tumors with the worst prognosis, the MES subtype GBMs [9,13].In addition, a previous study has provided evidence to suggest that the PN-like precursor cells can serve broadly as progenitors for GBMs and give rise to various subtypes of tumor cells [16].Collectively, these findings support targeting GBM cell stemness as a strategy for mitigating GBM progression.
Numerous studies have provided critical insights into the mechanism of stemness maintenance in GBM [17][18][19][20][21], and offered promising strategies, including epigenetic and metabolic approaches, for targeting this property of GBM cells [22][23][24][25][26][27].One alternative strategy for targeting stemness and plasticity of cancer cells is to promote their differentiation, an approach that shows great promise in blood cancers [28].Although employing this strategy for solid tumors has proven more complicated, recent studies have supported the feasibility of attenuating tumor cells' stemness/plasticity features and mitigating tumor progression in solid tumors via inducing their differentiation [29][30][31].In gliomas, IDH1 mutant-induced differentiation blockage has been exploited for therapeutic purposes [32,33], and metabolism-based strategies for directing GBM differentiation have been proposed [34,35].In addition to the benefit of attenuating stemness, plasticity, and self-renewal capacity, differentiation of tumor cells in solid tumors, including gliomas, can potentially promote tumor cell senescence and sensitize gliomas to immunotherapy, suggesting its multifaceted benefits from a therapeutic perspective [36][37][38][39].
Therapeutic differentiation of oligodendrocyte progenitor cells (OPCs) to become myelinating oligodendrocytes has been extensively pursued for treating multiple sclerosis (MS), a demyelinating disease that affects millions of patients worldwide.Several drugs are shown to effectively differentiate OPCs to myelinating oligodendrocytes in various in vivo disease models [40][41][42].Notably, the PN subtype GBMs display molecular signatures reminiscent of OPCs, both at the individual tumor level and at the single cell level, suggesting their OPC origin and progenitor features [9,12].We hypothesized that diminishing the progenitor features of these GBM cells via induction of differentiation can suppress their stemness and tumorigenicity.We identified clemastine, a low-toxicity, non-oncology drug for alleviating allergy symptoms with high potency in promoting OPC differentiation [40][41][42], as an agent that can be used to inhibit the propagation of BTIC cultures.Using BTIC cultures bearing PDGFRA amplification to represent PN subtype GBMs, we found that the tumor suppressive effects of clemastine were accompanied by diminished stemness of tumor cells and altered molecular signatures indicative of more differentiated states.Corroborating these findings, loss-of-function assays of putative pharmacological targets of clemastine identified the Emopamil Binding Protein (EBP), a crucial enzyme in the cholesterol biosynthesis pathway, as an essential protein for maintaining BTIC proliferation and stemness.Finally, clemastine suppressed the in vivo tumorigenicity of a mouse glioma cell line resembling the PN subtype of GBMs.This study identifies new candidate proteins in GBM cells that can be therapeutically targeted (i.e., EBP), and suggests that clemastine and other agents capable of inducing OPC differentiation can be leveraged for attenuating the stemness/plasticity of GBM cells, particularly those bearing PN features.Our findings support the principle of repurposing non-oncology drugs for inducing differentiation and thus mitigating GBM stemness.
shRNA-mediated gene knockdown cells were generated by transducing cell lines of interest with lentivirus, containing pLKO.1 vector with non-targeting control shRNA or shRNA sequences targeting either EBP or CHRM1 (purchased from Duke Functional Genomics Shared Resource, shRNA sequences listed in Table S1) at an MOI of 1 or 3. CRISPR/Cas9-mediated gene knockout cells were generated by transducing BTIC#102 cells with HRH1 and CHRM3 CRISPR/Cas9 knockout lentivirus.Luciferase-expressing cells were generated by transducing C266 cells with CMV-Firefly luciferase lentivirus (Cellomics Technology, #PLV-10003-50, Halethorpe, MD, USA) at an MOI of 1.All transduction reactions were supplemented with 6 µg/mL polybrene (MilliporeSigma, #TR-1003, Darmstadt, Germany).Virus-containing media were removed after 24 h of transduction and replaced with the appropriate media.Transduced cells were selected with puromycin (Millipore-Sigma, #P8833, Darmstadt, Germany) at 1 µg/mL for 5 days and subsequently maintained in 0.3 µg/mL puromycin.

Chemicals, Other Reagents and Plasmids
Chemicals and other molecular biology reagents that were used are listed in Table S2.Briefly, stock solutions of all remyelinating agents were prepared in dimethyl sulfoxide (DMSO), with clemastine and benztropine prepared at 10 mM and other remyelinating agents at 20 mM, and stored at −20 • C. The 200 µM working solutions were then prepared by diluting the 10 mM stock solutions with human or mouse standard NSC media and stored at −20 • C. CW3388 stock solutions were prepared in dimethyl sulfoxide (DMSO) at 2 mM and stored at 4 • C. Lathosterol stock solutions were prepared in ethanol at 12.5 mM, water-soluble cholesterol stock solutions were prepared in PBS at 4 mM, and stock solutions of H1R inhibitors were prepared in DMSO at 10 mM.The pcDNA3.1 V5-His A plasmid was a kind gift from Dr. Vidyalakshmi Chandramohan, and the pcDNA3.1+-C-(K)-DYKplasmid expressing Flag-EBP was purchased from GenScript (#OHu18817, Piscataway, NJ, USA).Primers for site-directed mutagenesis of the EBP gene were designed using the QuikChange Primer Design Program and listed in Table S1.Three EBP mutant constructs (E80K, R147H, W196S) were constructed from pcDNA3.1+-C-(K)-DYK-EBPplasmids using the QuikChange Lightning Multi Site-Directed Mutagenesis Kit (Agilent Technologies, #210519, Santa Clara, CA, USA) following the manufacturer's protocol.All plasmid sequences were verified by Sanger sequencing provided by Genewiz (South Plainfield, NJ, USA).

Plasmid Electroporation
Cells were transfected with pcDNA3.1 V5-His A, pcDNA3.1+-C-(K)-DYKplasmids encoding either wild-type EBP or one of the three mutant EBPs (E80K, R147H, W196S), EGFP-hGal3, or PM-GFP plasmids using Neon Transfection System 10 µL Kit (Thermo Fisher Scientific, #MPK1025, Waltham, MA, USA) following the manufacturer's protocol.In brief, cells were prepared and resuspended in buffer R at a concentration of 100,000 cells/10 µL buffer R.Then, 10 µL of cell suspension was added to each plasmid solution containing 0.5 µg or 1 µg of plasmid DNA.The cell-plasmid DNA mixtures were then subjected to electroporation using the Neon Transfection System and transferred into laminin-coated 24 well plates containing 500 µL of prewarmed standard NSC media in each well.Transfected cells were selected with G418 sulfate (ThermoFisher Scientific, #10131035, Waltham, MA, USA) at 250 µg/mL for 5 days and subsequently maintained in 100 µg/mL G418 sulfate.

BTIC Proliferation Assay, Cell Cycle Analysis and ELDA
Cell proliferation assays were performed by seeding the cell lines of interest in their respective media, with either vehicle or drugs of interest and without any antibiotics, at a density of 1500-3000 cells/well and a final volume of 200 µL/well in laminin-coated 96-well clear flat bottom plates.Cells were incubated overnight in a humidified incubator at 37 • C and 5% CO 2 , and then placed in the IncuCyte S3 Live-cell analysis system (Sartorius, Göttingen, Germany) and scanned every 4 h using the whole well, phase contrast acquisition mode (4× objective) for 7-14 days.Phase area confluence was obtained and calculated by the IncuCyte system, normalized to day 0 to generate relative phase object area (fold change), and then presented as mean ± S.E.M.
For cell cycle analyses, cells were re-suspended in 500 µL of ice-cold PBS and triturated to obtain a single cell suspension.Ethanol fixation was performed by first adding 3 mL of ice-cold 70% ethanol to the cell suspension in a drop-wise manner while vortexing, and then incubating the cell-ethanol mix at −20 • C for 1 h.The fixed cells were washed with ice-cold PBS three times, resuspended in 500 µL FxCycle™ PI/RNase Staining Solution (ThermoFisher Scientific, #F10797, Waltham, MA, USA), and incubated at room temperature for 30 min while protected from light.Samples were analyzed with BD Fortessa X-20 with the BD FACSDiVa software (BD Biosciences, San Jose, CA, USA).ELDAs with or without pretreatment were performed as previously described [43,44].

In Vivo Drug Treatment
Orthotopic intracranial tumor implantation and in vivo bioluminescent imaging were performed as previously described [43,44].Briefly, C266 mouse glioma cells (~100,000 cells, expressing luciferase) were mixed 1:3 with methylcellulose and injected into the right caudate nucleus of female athymic nude mice (Jackson Labs, strain #:007850, Bar Harbor, ME, USA), and tumor treatment was initiated 4 days post-implantation.Mice were treated with vehicle control or clemastine at 30 mg/kg via intraperitoneal (i.p.) injection daily for five days per week (15% DMSO in PBS was used as the vehicle control).In vivo drug response was monitored weekly by bioluminescent imaging of mice and analyzed using the Living Image software (PerkinElmer, Waltham, MA, USA).Mouse body weight was documented at least every two days.All animal experiments were performed in accordance with protocols approved by the Duke University Institutional Animal Care and Use Committee (IACUC), Protocol #A133-19-06 (approved 27 June 2019).

Reverse Transcription and Real-Time Quantitative PCR (RT-qPCR)
BTIC cultures were plated at a density of 80,000-120,000 cells/well in laminin-coated 6 well plates in standard NSC media and incubated overnight to allow them to adhere.On the next day, cells were treated with either vehicle or drugs.Cells were passaged every 5-7 days into fresh media containing respective treatments, and then harvested on day 9-14 for total RNA extraction using the Quick-RNA Miniprep kit (Genesee Scientific, #11-328, El Cajon, CA, USA) following the manufacturer's protocol.The concentration and the A260/A280 ratio of each RNA sample was measured with a Nanodrop Lite Spectrophotometer (Thermo Fisher Scientific, #ND-LITE, Waltham, MA, USA).RNA samples with concentration > 30 ng/µL and A260/280 ≥ 2.0 were then reverse-transcribed to cDNA using the EcoDry cDNA synthesis kit (Takara Bio, #639548, Shiga, Japan) following the manufacturer's protocol.Real-time quantitative PCR was performed using gene-specific qPCR primers listed in Table S1 and the KAPA SYBR Fast qPCR Master Mix (2×) kit (Kapa Biosystems, #KK4602, Wilmington, MA, USA), and ran on the CFX96 Real-Time PCR Detection System (BIO-RAD, Hercules, CA, USA).The results were analyzed using the CFX Maestro Software (BIO-RAD, Hercules, CA, USA).The housekeeping genes ACTB, B2M, and GAPDH were used as internal controls for gene expression normalization.

Protein Extraction and Immunoblotting
Total cellular protein was extracted from cell pellets using 1% SDS lysis buffer.Briefly, frozen cell pellets were resuspended with 80-120 µL 95 • C preheated SDS lysis buffer (1% SDS, 50 mM NaF, 1 mM Na 3 VO 4 in PBS, pH 7.4) and triturated until the suspensions were less gluey.The cell suspensions were then sonicated with Bioruptor Standard (Diagenode, #UCD-200, Denville, NJ, USA) for 5 min at high intensity to achieve complete cell lysis, and heated at 95 • C for 5 min.Cell lysates were centrifuged to remove residual cell debris, and supernatants, which were the protein extracts, were collected into new Eppendorf tubes.Protein concentrations were determined using the Pierce BCA Protein Assay kit (Thermo Fisher Scientific, #23225, Waltham, MA, USA) following the manufacturer's protocol.The residual protein extracts were mixed with 4× Laemmli sample buffer (BIO-RAD, #1610747, Hercules, CA, USA) at a ratio of 1:3, heated at 95 • C for 5 min, cooled on ice for 5 min, and stored at −20 500 µL of NuPAGE Antioxidant (Thermo Fisher Scientific, #NP0005, Waltham, MA, USA).After gel electrophoresis, proteins were transferred from gels to nitrocellulose membranes (BIO-RAD, #1620115, Hercules, CA, USA) by semi-dry transfer methods using the Trans-Blot Turbo Transfer System (BIO-RAD, #1704150, Hercules, CA, USA).The membranes were then blocked with 5% non-fat dry milk (Genesee Scientific, #20-241, El Cajon, CA, USA) in TBST (0.1% Tween 20 in TBS) for 1 h at room temperature, and incubated overnight at 4 • C with primary antibodies diluted in antibody dilution buffer (5% BSA, 0.02% NaN 3 in TBST) according to the manufacturer's suggested dilutions (detailed information for antibodies and their dilutions were listed in Table S3).On the next day, the membranes were incubated with appropriate horseradish peroxidase-conjugated secondary antibodies (Cell Signaling Technology, anti-rabbit IgG, #7074S; anti-mouse IgG, #7076S, Danvers, MA, USA) diluted in TBST for 1 h at room temperature.Then, the chemiluminescent signals were enhanced by incubating the membranes with SuperSignal West Pico PLUS Chemiluminescent Substrate (Thermo Fisher Scientific, #34580, Waltham, MA, USA) following the manufacturer's protocol.The membranes were imaged with ChemiDoc MP System (BIO-RAD, #1708280, Hercules, CA, USA) and analyzed with Image Lab Software (BIO-RAD, Hercules, CA, USA).To re-probe the membranes, membranes were stripped using a Restore Western Blot Stripping Buffer (Thermo Fisher Scientific, #21059, Waltham, MA, USA) and incubated with other primary antibodies following the manufacturer's protocol.
Images from PDGFRA and GALC immunofluorescence staining of the vehicle or clemastine-treated BTIC#102 cells were analyzed using the CellProfiler software (version 4.2.1, Broad Institute, Cambridge, MA, USA) [47].For quantification of mean fluorescence intensity of PDGFRA and GALC, "Cell" objects were first identified using the DAPI signal for nuclear regions and the PDGFRA/GALC signal to outline cell borders.The mean fluorescence intensity of PDGFRA and GALC per "Cell" object was measured, and the mean of mean fluorescence intensity per "Cell" object per field of view was then calculated.For quantification of percent positivity of PDGFRA and GALC, PDGFRA + and GALC + Cancers 2023, 15, 4619 7 of 22 "Cell" objects were first identified using the mean fluorescence intensity per "Cell" object as thresholds, which were set using isotype-stained images as negative controls.The percentage of PDGFRA and GALC positive cells were then calculated as number of positive Cell object total number of Cell Object per field of view.To quantify the frequency distribution of PDGFRA + and/or GALC + cells, GALC positivity was determined in every PDGFRA − and PDGFRA + cell from all fields of view regardless of treatment status, and the contingency table was generated based on the calculated frequency distribution.

RNA-Seq Data and Pathway Analysis
Next generation sequencing for mRNA-seq, including library construction and sequencing, was provided by Novogene Corporation Inc. (Sacramento, CA, USA).NovaSeq 6000 was used for PE150 sequencing.mRNA-seq data were analyzed using the Galaxy web platform via the public server at usegalaxy.org[48].The workflow of RNA-seq analysis was adapted from a previously described procedure [49].Briefly, raw data were trimmed by Trim Galore, aligned by Hisat2, re-assembled by StringTie, and differentially expressed genes (DEGs) were analyzed by Deseq2.All parameters were set at default.For pathway enrichment analysis, differentially expressed genes (DEGs, sorted by adjusted p-value < 0.05) were imported into ShinyGO [50] (http://bioinformatics.sdstate.edu/go/,accessed on 14 October 2021).For pathway analysis, KEGG and Panther pathways were used, and genes with adjusted p-values > 0 were input as the background gene list.For GSEA, all genes or genes with FDR ≤ 0.25 from human or mouse mRNA-seq data, respectively, were ranked by fold change and subjected to GSEA preranked analyses following the established protocol [51] using the GSEA software (version 4.1.0,Broad Institute, Cambridge, MA, USA).For pathway impact (two-evidence) analysis, all genes with adjusted p-values ≤ 0.05 were selected and subjected to analysis using ROntoTools [52], and pathways were ranked by combined FDR p-value (pComb.fdr).

Statistical Analysis
All in vitro proliferation assays, cell cycle analyses, ELDA, quantitative RT-PCR, immunoblotting, immunofluorescence staining, and in vivo clemastine treatment experiments were repeated in at least two independent experiments.The number of independent samples was noted in the figure legends.Data were presented as mean ± standard error of the mean (S.E.M.), except for gene expression data, which were presented as geometric mean ± geometric standard deviation (S.D.), and ELDA data, which were shown as trend lines ± two-sided 95% confidence intervals.Mean values between two groups were compared by student t-tests (with Welch's correction when variances were deemed significantly different by F tests) or non-parametric Mann-Whitney tests.Mean values between 3 or more groups were compared by one-way or two-way ANOVA followed by Dunnett's or Tukey's multiple comparisons tests, respectively.Mean values between two or multiple groups over time were compared to the control group by two-way repeated measures ANOVA followed by Sidak's or Dunnett's multiple comparisons tests, respectively.Mean values of multiple groups were compared between all possible group pairings by two-way repeated measures ANOVA, followed by Tukey's multiple comparisons tests.Geometric mean values of gene expression fold change were log 2 -transformed and analyzed by one-sample t-tests compared to zero.Survival analyses were performed using Log-rank (Mantel-Cox) tests.Contingency tables were analyzed using one-sided chi-square tests.All tests were twosided, if not otherwise specified, and deemed statistically significant when p-values < 0.05.Normality was tested before conducting any parametric test using Shapiro-Wilk normality tests.For proliferation assays, significance was calculated using the data from the last timepoints unless stated otherwise.All statistical analyses were performed using GraphPad Prism software (version 9.3.1,San Diego, CA, USA) except for ELDA (analyzed by ELDA online software: https://bioinf.wehi.edu.au/software/elda,accessed on 17 April 2021), Gliovis gene expression data (analyzed by Gliovis portal: http://gliovis.bioinfo.cnio.es/,accessed on 18 October 2021), RNA-seq, GSEA, and pathway analyses.

Clemastine Suppresses the Propagation of Patient-Derived BTIC Cultures
Several drugs with OPC-differentiating activity have been shown to effectively stimulate normal OPC differentiation and subsequent remyelination in multiple in vivo disease models [40][41][42].These findings prompted us to test the effects of these drugs on GBM cells bearing OPC features.The PN subtype of GBM, typified by OPC-like transcriptomes, is characterized by IDH1 mutations or PDGFRA amplification, while EGFR amplification associates with the CL subtype [9].Therefore, we first identified GBM lines with predominantly EGFR amplification (EGFR + , representing CL GBMs), both PDGFRA and EGFR amplifications (in agreement with the heterogeneous nature of GBMs that are observed), as well as lines that bear predominantly PDGFRA amplification (PDGFRA + , representing PN GBMs) (Figure S1A).We then cultured two PDGFRA + lines (BTIC#102 and BTIC#148, with ~25× and ~3× PDGFRA amplification, respectively) in serum-free neural stem cell medium to maintain their stemness properties as BTIC cultures, and determined the effects of a small panel of OPC-differentiating drugs on BTIC propagation.We found that all of these drugs suppressed the propagation of BTICs in a dose-dependent manner (Figure S1B,C).Among them, tamoxifen, benztropine, and clemastine were previously found to promote OPC differentiation by affecting the cholesterol biosynthesis pathway [53].Of note, clemastine is a safe, over-the-counter (OTC) allergy relief medicine that can cross the blood-brain barrier (BBB) and induce remyelination [40,41,54].It has shown promising therapeutic benefits via differentiation induction in demyelinating models [40][41][42] and, most remarkably, in clinical trials for MS patients [54,55].These findings prompted us further to characterize the inhibitory effects of clemastine.We found that the effects of clemastine on BTICs were accompanied by cell morphology alterations (Figures 1A,B and S2A,B) and suppressed cell proliferation (Figures 1C and S2C).Notably, following 10-day clemastine pre-treatments, the antiproliferative effects of clemastine persisted even after clemastine was removed from the media (Figures 1D and S2D).Collectively, these results suggest that clemastine induced lasting alterations in PDGFRA + BTICs and suppressed the propagation of these tumor cells.

2023, 15, x FOR PEER REVIEW 9
implicate clemastine as a BTIC-suppressive drug and provide a basis for repurposing mastine, and potentially other remyelinating agents, for targeting GBMs.While we had focused on PN-like BTICs that are believed to serve broadly as progenitors of GBMs [16] and transition between one subtype to another (e.g., PN to MES GBMs) [14,15], the intratumoral heterogeneity of GBMs consisting of various subtypes [7][8][9][10][11][12] prompted us to examine the effects of clemastine in primary BTIC cultures resembling other subtypes.We used two primary EGFR + BTIC cultures (representing the CL subtype: BTIC#127 and BTIC#095) (Figure S1A) and two primary BTIC cultures that were shown to bear MES transcriptomic signatures (GSC-1123 and JK-046 [45,46]) for this purpose.We found that BTIC cultures resembling the CL and MES subtypes were also susceptible to clemastine treatment, as indicated by reduced cell proliferation (Figure 1E), suggesting that the antiproliferative effects of this drug impact BTICs broadly and are independent of tumor subtypes.Importantly, unlike the tumor cells tested, the proliferation of normal human astrocytes was not affected by the treatment (Figure 1F).Collectively, these results implicate clemastine as a BTIC-suppressive drug and provide a basis for repurposing clemastine, and potentially other remyelinating agents, for targeting GBMs.

Clemastine Attenuates the Stemness and Progenitor Cell Features of PDGFRA + BTICs
Clemastine has been shown to promote the differentiation of OPCs to mature oligodendrocytes in treating MS [53,56].This prompted us to investigate whether the effects of clemastine on BTICs involved such a differentiation induction mechanism.We used PDGFRA + BTIC cultures as the prototypic models with the rationale that the OPC-like status of these tumor cells can facilitate assessing their progenitor versus differentiated states.We treated the PDGFRA + BTIC cultures with clemastine and evaluated the expression of key markers indicative of neural stem cell (NSC) and OPC properties.In agreement, clemastine treatment resulted in the decreased expression of marker genes for NSCs (NES [57], SOX2 [58]), OPCs (OLIG2 [59], PDGFRA [60], CSPG4 [61]), and negative regulators of OPC differentiation (e.g., NOTCH1, CSPG5 [62]) in the two PDGFRA + BTIC cultures (Figures 2A and S3A).Intriguingly, MOG, an oligodendrocyte maturation marker gene [63], displayed upregulated expression in clemastine-treated cells (Figure 2A).The reduced expression of these molecular markers of NSCs and OPCs was confirmed at their protein levels by immunoblots (Figures 2B and S3B) and immunofluorescent staining (Figures 2C and S3C).To further investigate cell state changes following clemastine treatment, we utilized BTIC#102, the line with a high level of PDGFRA amplification and PDGFRA expression (i.e., more OPC-like), to facilitate studying the differentiating effects.
The diminished progenitor states and potential differentiation of tumor cells were corroborated by the induction of GALC protein expression, another oligodendrocyte maturation marker [63], in parallel with the reduced level of PDGFRA protein (Figure 2D,E).The increased presence of GALC + and reduced presence of PDGFRA + cells (Figure 2F), together with the mutually exclusive staining pattern of these proteins (Figure 2G), suggest these tumor cells indeed became more differentiated upon clemastine treatment.The altered expression of marker genes associated with progenitor or differentiated cells prompted us to examine the changes in the global gene expression (mRNA-seq) of the two PDGFRA + BTIC cultures in response to clemastine.Consistent with our observations mentioned above, transcriptomic profiling from the mRNA-seq data revealed an overall downregulated expression of most dominant marker genes for OPCs [64] in response to clemastine treatment in both BTIC cultures (Figures 2H and S3D).mRNA-seq analyses also revealed upregulated expression of a fraction of marker genes for newly formed oligodendrocyte (NFO) and myelinating oligodendrocytes (MO) [64] in response to clemastine treatment (Figures 2I and S3E).Gene Set Enrichment Analysis (GSEA) confirmed the clemastine-induced downregulation of genes associated with OPCs (Figure S3F,G) [65], and downregulation of genes in the protein translation machinery (Figure S3H,I), which serves a critical role in cancer cell stemness/plasticity and as a promising therapeutic target [66,67].

EBP, A Pharmacological Target of Clemastine, Is Essential for BTIC Propagation
To assess the functional consequence of the altered gene expression profiles described above, we determined the effects of clemastine on the stem-like properties of BTICs, as measured by their self-renewal capability [68,69], via the Extreme Limiting Dilution Assay (ELDA) [70].We found that clemastine treatment diminished the self-renewal capability of BTICs (Figure 2J).More notably, reminiscent of the lasting effects of clemastine on BTIC propagation, BTICs that were exposed to clemastine and then maintained in clemastine-free media also displayed a diminished self-renewal capacity, suggesting a persistent effect of clemastine on BTIC stemness (Figure 2K).Together with the altered gene expression profiles, these results suggest that clemastine attenuates the stemness/progenitor properties of BTICs by promoting differentiation.

EBP, A Pharmacological Target of Clemastine, Is Essential for BTIC Propagation
Pharmacologically, clemastine is known to act as an antagonist for histamine H1 receptor (H1R, encoded by HRH1) and muscarinic receptors (M1R-M5R, encoded by CHRM1-CHRM5, respectively) [71], and to target EBP, an enzyme in the cholesterol biosynthesis pathway.The well-established antagonist activity of clemastine against H1R and the inhibitory effects of M1R and M3R signaling on oligodendrocyte differentiation [72,73] prompted us to determine the expression profiles of these three receptor genes (HRH1, CHRM1, and CHRM3) in gliomas.We found that HRH1 and CHRM3 were expressed at higher levels in GBM (grade IV) in comparison to lower-grade gliomas (grade II and III), while CHRM1 displayed the opposite trend (Figure S4A).This result led us to experimentally assess the roles of H1R and M3R in BTICs.First, two second-generation H1R-antihistamines, fexofenadine (Allegra) and cetirizine (Zyrtec) [74], were found to have no suppressive effects on BTICs, suggesting these cells were not affected by the H1R signaling blockade (Figure S4B).Furthermore, CRISPR/Cas9-mediated individual gene knockout of HRH1 or CHRM3 (Figure S4C) did not affect the propagation of BTICs or the susceptibility of BTICs to clemastine (Figure S4D).As the tumor-supporting roles of each gene may be affected by the variable genetic makeups specific for each tumor model, these results did not conclusively rule out the involvement of these receptors in mediating the effects of clemastine.Nevertheless, these findings led us to investigate other potential targets of clemastine.Interestingly, the pharmacological inhibition of EBP has been shown to effectively induce the differentiation of OPCs to become oligodendrocytes [53,56], which prompted us to examine the roles of EBP in GBM (Figure S5A).
First, analysis of mRNA-seq data from The Cancer Genome Atlas (TCGA) [75] via the Gliovis data portal [76] revealed that, distinctively, from genes encoding other enzymes in the cholesterol biosynthesis pathway, EBP displayed an upregulated expression in GBMs in comparison to non-tumor tissues, as well as to lower grade gliomas (grade II and III) (Figure S5B).Notably, two additional properties distinguish EBP from other enzymes in the cholesterol biosynthesis pathway: EBP mutations have been linked to Conradi-Hünermann-Happle (CHH) syndrome [77], and it is known to be a protein target of numerous structurally diverse pharmacological agents [78][79][80].In fact, EBP is the common pharmacological target shared by seven out of the eight OPC-differentiating drugs included in our initial test [53,56,78,79] (Figure S1B,C).These unique features, together with the therapeutic potential of targeting EBP, led us further to assess the pathological functions of this protein in GBM.
We used the patient-derived BTIC cultures to determine the roles of EBP in GBM.The knockdown of EBP in BTICs resulted in significant suppression of BTIC propagation, which was accompanied by altered cell morphology (Figures 3A,B and S6A,B) and reduced expression levels of featured OPC marker genes, including PDGFRA, CSPG4, and OLIG2 (Figure S6C).Conversely, echoing the finding from the loss-of-function assays, an overexpression of exogenous EBP promoted BTIC propagation (Figures 3C and S6D,E).Several pathogenic, missense mutations of EBP have been identified in CHH patients with various degrees of reduction in enzymatic activities [77,81,82].Therefore, we assessed if enzymatic activity is required for the promoting effect of EBP by overexpressing three of the EBP mutants, EBP-E80K, EBP-R147H, EBP-W196S, or the wild-type EBP in BTIC#102 cells (Figure S6F,G).Of note, levels of these mutant EBP proteins were lower than the wild-type EBP protein (Figure S6F,G), likely due to their various stabilities.Nevertheless, among the three mutants, which had comparable expression levels, EBP-E80K, the mutant known to have minimal enzymatic activity retained, had lower potency in promoting the growth of the GBM culture compared to EBP-R147H, a mutant that possesses partial enzymatic activity [77,83] (Figure 3D).Taken together, these results suggest that EBP is essential for BTIC propagation and that the GBM-promoting effects of EBP likely require its enzymatic activity.
Cancers 2023, 15, x FOR PEER REVIEW 13 of 24 suppressive effects of clemastine on BTIC propagation (Figures 3E and S6H).Furthermore, the presence of the immediate metabolic product of EBP, lathosterol, also partially countered the suppressive effects of clemastine on BTIC propagation.Of note, the addition of exogenous lathosterol alone at a higher dose moderately suppressed tumor cell propagation (Figures 3F and S6I), suggesting that an optimal dose of this metabolite and maintaining the homeostasis of the cholesterol pathway are essential for the optimal growth of tumor cells.Finally, we found that a stable overexpression of EBP in primary BTICs conferred partial protection from clemastine (Figure 3G).Collectively, the above results suggest that BTIC propagation depends on EBP, a pharmacological target of clemastine, and the homeostasis of the cholesterol pathway, and provide a potential mechanism underlying the suppressive effects of clemastine.Several additional lines of evidence corroborate the findings mentioned above and support that suppressing EBP partially mediates the effects of clemastine on BTICs.First, the addition of exogenous cholesterol to the culture media partially rescued the suppressive effects of clemastine on BTIC propagation (Figures 3E and S6H).Furthermore, the presence of the immediate metabolic product of EBP, lathosterol, also partially countered the suppressive effects of clemastine on BTIC propagation.Of note, the addition of exogenous lathosterol alone at a higher dose moderately suppressed tumor cell propagation (Figures 3F and S6I), suggesting that an optimal dose of this metabolite and maintaining the homeostasis of the cholesterol pathway are essential for the optimal growth of tumor cells.Finally, we found that a stable overexpression of EBP in primary BTICs conferred partial protection from clemastine (Figure 3G).Collectively, the above results suggest that BTIC propagation depends on EBP, a pharmacological target of clemastine, and the homeostasis of the cholesterol pathway, and provide a potential mechanism underlying the suppressive effects of clemastine.

Clemastine Suppresses Tumorigenicity in A Mouse Glioma Model Representative of PN GBMs
We postulated that a GBM model with a better-defined genetic makeup could provide advantages to further characterize the tumor suppressive effects of clemastine.For this purpose, we generated mouse NSC lines from a conditional Trp53 knockout mouse model previously described [84,85].NSCs derived from these models underwent adenoviral Crerecombinase mediated Trp53 deletion followed by the overexpression of exogenous Pdgfb, an oncogene that has been shown to induce gliomas that resemble PN GBMs [86][87][88][89][90].The genetically modified NSC lines were then orthotopically implanted into mice to generate PDGFB-driven, aggressive glioma tumors with high penetrance and short latency (Figure S7A).Gene expression profiling (mRNA-seq) of the tumor tissues and subtyping analysis utilizing the previously defined subtype-specific gene signatures [90] confirmed their predominant PN features, and the detectable expression of genes associated with other subtypes recapitulated the heterogeneous nature of GBMs (Figure S7B).A glioma cell line derived from these tumors (named C266) was used to investigate its dependency on the EBP protein and susceptibility to clemastine, benztropine, and a small-molecule drug (CW3388) known to inhibit mouse EBP [56].Results from these experiments using the mouse glioma cell line corroborate findings from human BTICs, as elaborated in the following: First, the knockdown of Ebp resulted in attenuated growth (Figures 4A and S7C,D) and reduced self-renewal capacity (Figure 4B) of the mouse glioma cell line.Second, treating C266 cells with CW3388 resulted in decreased cell proliferation (Figure 4C).Similarly, the mouse glioma line also displayed suppressed propagation in response to clemastine and benztropine treatment, recapitulating the findings from human BTICs (Figures 4D and S7E).Third, echoing the effect of EBP's abundance on human BTIC's susceptibility to clemastine, Ebp knockdown rendered mouse glioma cells more sensitive to clemastine and CW3388 (Figures 4E and S7F), suggesting that clemastine acted in the same manner as CW3388.Taking advantage of the aggressive nature of this tumor cell line in orthotopic mouse models, we tested the suppressive effects of clemastine in vivo.We found that the clemastine treatment led to a delayed progression of tumors in vivo, as evidenced by the slower tumor progression (Figure 4F), and a small subset of mice (20%) that had extended survival with measurable tumor signal two months post-implantation (at which point, mice were terminated) (Figure 4G).The reasons for this delayed yet extended benefit remained unclear, and we speculate that the treatment likely altered the progenitor properties of tumor cells instead of acutely killing them in vivo.Finally, we noted that mice subjected to this treatment regimen displayed no different body weight compared to the control mice, suggesting its minimal toxicity (Figure S7G).
clemastine treatment led to a delayed progression of tumors in vivo, as evidenced by th slower tumor progression (Figure 4F), and a small subset of mice (20%) that had extende survival with measurable tumor signal two months post-implantation (at which poin mice were terminated) (Figure 4G).The reasons for this delayed yet extended benefit re mained unclear, and we speculate that the treatment likely altered the progenitor proper ties of tumor cells instead of acutely killing them in vivo.Finally, we noted that mice sub jected to this treatment regimen displayed no different body weight compared to the con trol mice, suggesting its minimal toxicity (Figure S7G).

Clemastine Treatment Broadly Alters Multiple Signaling Pathways in Glioma Cells
While clemastine has multiple pharmacological targets and its effects on cell signaling pathways are expected to be broad and dependent on the genetic composition of GBM cells, we postulated that identifying pathways that were perturbed by clemastine and/or CW3388 in the mouse glioma cell line could provide additional insights into the molecular mechanisms of clemastine.Therefore, we performed mRNA-seq to identify differentially expressed genes in clemastine or CW3388-treated C266 cells (Table S4).Several findings were noted.First, among the 10,292 genes identified, 18.3% (1888 genes) of them were affected by clemastine.In comparison, among those genes affected by CW3388 (2759 genes), the number increased to 45.7% (1261 genes), suggesting genes affected by clemastine were enriched in the CW3388-responsive gene set (p < 0.0001).Additionally, among genes affected by both agents, an overwhelming majority of them (1210 out of 1261; 96.0%) displayed the same directional changes (i.e., up-or down-regulation) in response to both agents, further supporting that clemastine and CW3388 act on overlapping genes/pathways in the tumor cells.Second, GSEA of differentially expressed genes presented a positive enrichment of genes corresponding to oligodendrocytes [91] in response to either clemastine or CW3388, in agreement with their expected effects on a tumor cell's progenitor identity (Figure S8).Third, in clemastine-treated cells, KEGG pathway analysis identified pathways in cancers, metabolism, Wnt, and PI3K-AKT signaling pathways (Table S5).Notably, similar analyses comparing the CW3388-treated versus the vehicle-treated cells identified an overlapping set of KEGG pathways (Table S5).Panther pathway analysis pinpointed smaller numbers of altered pathways using the same criteria (FDR 0.05), and revealed that a majority of the CW3388-affected pathways were also identified as being perturbed by clemastine (Table S5).These results suggest that the impacts of clemastine overlap with and exceed those of CW3388, in agreement with its expected pharmacological activity targeting multiple proteins, including EBP.
We employed signaling pathway impact analysis, which incorporates an overrepresentation analysis and functional class scoring and was shown to have superior specificity and sensitivity in identifying altered pathways [52], to analyze differentially expressed genes of patient-derived BTIC cultures in response to clemastine (Table S6).This analysis identified numerous KEGG pathways that were affected in the BTIC cultures (196 and 154 were identified in BTIC#102 and BTIC#148, respectively).Although the pathways varied between these models, likely due to the various genetic background and different basal gene expressions in matched control cells, 142 common pathways were identified in both BTIC cultures, including cell adhesion, cholesterol metabolism, and the PI3K-AKT signaling pathway (Table S7).Collectively, these pathway analyses suggest the pleiotropic consequences of clemastine treatment at the molecular level, and also highlight the substantial alterations in cellular processes and oncogenic signaling pathways in tumor cells, including cell metabolism, the Wnt signaling pathway, and the PI3K-AKT pathway.While such assays could not distinguish pathway alterations that contributed to the tumor suppression from those that resulted from cell adaptation/population evolution, these results resonate with our findings in the mouse glioma cell line, and yield additional evidence to support the tumor suppressive effects of clemastine.

Discussion
BTICs are key drivers of GBM's resistance to treatments, progression, and recurrence.In this study, we showed that clemastine, an over-the-counter drug used for alleviating allergies, suppresses the propagation of patient-derived BTICs bearing PDGFRA amplification (as a surrogate representing the PN GBMs), as well as those resembling CL or MES subtypes.We further used the PN/OPC-like BTICs to demonstrate that this drug attenuates the stemness/progenitor features of tumor cells.Finally, we showed that pharmacological targets of clemastine, as highlighted by EBP, are essential for the proliferation of these tumor cells.
The suppressive effects of clemastine on the BTICs described in the current study corroborate the findings from previous research on clemastine.First, it was shown that clemastine could compromise the lysosomal membrane integrity in glioma cells to suppress glioma cell tumorigenicity [92].This study aligns with our findings in supporting that this drug can act on multiple pathways and cellular processes to exert tumor suppressive effects, and provide a foundation for investigating the detailed mechanisms in further studies.Second, the induction of normal OPC differentiation into myelinating oligodendrocytes by clemastine has yielded promising therapeutic benefits in clinical trials in MS patients, validating the in vivo therapeutic efficacy of clemastine and its permeability across the blood-brain barrier [54,55].
We note that, unlike in normal cells, genetic and epigenetic alterations in tumor cells most likely preclude the possibility of accomplishing terminal differentiation in the latter.Nevertheless, the suppressive effects of clemastine on GBM stemness could presumptively subvert the resistance mechanisms (e.g., heightened DNA repair machinery [93] and metabolic plasticity [94]), address the heterogeneity issue attributed to BTIC multipotency, and provide therapeutic benefits.Specifically, we propose to fully exploit the therapeutic potential of clemastine by combinatorial strategies with therapeutics such as (i) standard chemo-and radiation therapy; and (ii) drugs that target other subtypes of GBMs (e.g., anti-mesenchymal drug PTC596 [95], anti-classical drugs gefitinib [96] and nimotuzumab [97]).These strategies can potentially mitigate the emergence or selection of the evolved populations, such as those driven by EGFR oncogenic signaling [15].
The fact that an OTC medication, such as clemastine, demonstrates effectiveness in differentiating and attenuating stemness in BTICs offers an intriguing opportunity for GBM treatment development.In comparison with traditional chemo-and radiation therapy, the differentiation-based mechanism of clemastine is less dependent on proliferation status, and clemastine boasts the established BBB permeability and an incomparable safety profile with only minor adverse effects such as sedation and xerostomia (oral dryness).In support, many other OTC medications have gained sufficient experimental evidence and entered clinical trials for the investigation of GBM treatment (Table S8 [98][99][100][101][102][103]), and clemastine represents a promising therapeutic agent worthy of further investigation and awaiting to join the list.
The loss-of-function of EBP, a protein that is pharmacologically inhibited by clemastine [53,56], recapitulated the suppressive effects of clemastine, as indicated by both attenuated BTIC stemness and the accompanied differential gene expression.Taken together with the dependency of BTICs on EBP, our study implicates sterol metabolism in GBM stemness maintenance.Prior studies support that inhibitors for certain enzymes in the cholesterol biosynthesis pathway, including EBP, can promote remyelination, linking this pathway and the involved intermediate metabolites to promoting normal OPC differentiation [53,56].In addition, lipid biosynthesis and the cholesterol pathway have been demonstrated to be promising therapeutic targets for treating GBMs [104][105][106][107][108]. For instance, intermediate metabolites in the cholesterol pathway were found to exert the tumor suppressive effects via activating endogenous liver X receptors (LXRs) in GBM [104,107] and another type of brain tumor, diffuse intrinsic pontine glioma (DIPG) [109].Thus, our results are consistent with the previous literature that highlights the tumor cell's dependency on the homeostasis of the cholesterol pathway.
Three outstanding issues and limitations remain to be further investigated and addressed.First, although EBP inhibition likely contributed to the suppressive effects of clemastine, we cannot rule out the potential promiscuity of clemastine and that additional pathways (e.g., histamine/muscarinic receptor signaling) might participate in mediating clemastine's effects, as suggested by the remaining sensitivity of shEbp mouse glioma cells to clemastine (Figures 4E and S7F).The underlying mechanism of clemastine-mediated EBP inhibition also remains to be studied.This issue is complicated by EBP's unique multi-drug binding capacity and the possibility that EBP lowers the intracellular concentration of drugs that inhibit its activity [83,110,111], which provide an alternative explanation for the findings in which the knockdown of EBP conferred the tumor cell's sensitivity (and vice versa) to its inhibitors (such as clemastine and CW3388).Second, the essentiality of EBP in BTICs needs to be assessed and validated in vivo by performing limiting dilution experiments using the isogenic Ebp-knockdown mouse glioma cell line in our mouse models.Third, illustrating the molecular mechanisms of the tumor suppressive effects of EBP inhibition will be critical for assessing the potential of this enzyme as a therapeutic target.It is possible that EBP inhibition alters the intermediate metabolite composition in tumor cells (cell-autonomous) and/or in the tumor microenvironments (e.g., non-tumor cells).Alternatively, the upregulated expression of EBP in GBMs, a trend that is opposite to several other enzymes in the cholesterol biosynthesis pathway, as well as the complicated pharmacological property of this protein [78,79], raise the possibility that its inhibition affects cellular processes other than simply disrupting the biosynthesis pathway of cholesterol.

Conclusions
In summary, findings from this study support the feasibility of inducing differentiation for targeting the stemness of gliomas bearing OPC features, which include not only PN-like GBM, but also oligodendroglioma and DIPG [112,113].We provide evidence to support exploiting clemastine, or other non-oncology, low-toxicity drugs with similar activity in targeting EBP and/or inducing differentiation [78,79,114], for targeting GBM stemness.In addition, the study also nominates EBP as a therapeutic target in GBM.During the preparation of this manuscript, the FDA has granted orphan drug designation (ODD) for DSP-0390, a newly developed inhibitor of EBP (https://www.targetedonc.com/view/fdagrants-orphan-drug-designation-to-dsp-0390-in-brain-cancer,accessed on 19 July 2022), for which a phase I clinical trial for treating recurrent high-grade gliomas has been initiated (NCT05023551) [115].The results presented in this study provide evidence to support such actions.We propose that further studies assessing the underlying pathogenic mechanism of EBP and the therapeutic efficacy of EBP inhibitors on various subtypes of GBM cells in combination with other treatment approaches are warranted.

Supplementary Materials:
The following supporting information can be downloaded at: https:// www.mdpi.com/article/10.3390/cancers15184619/s1, Figure S1: OPC-differentiating agents inhibit the proliferation of patient-derived BTIC cultures bearing PDGFRA amplification; Figure S2: Clemastine suppresses the proliferation of patient-derived BTIC cultures bearing PDGFRA amplification; Figure S3: Clemastine attenuates the stemness and progenitor cell features of PDGFRA + BTIC cells; Figure S4: Assessing the roles of potential pharmacological targets of clemastine in BTICs; Figure S5: EBP displays distinct expression patterns in gliomas in comparison to other genes in the cholesterol biosynthesis pathway; Figure S6: EBP is essential for BTIC propagation; Figure S7: Loss of Ebp impairs the growth of C266 mouse glioma cells derived from p53-null mouse NSCs with PDGFB overexpression; Figure S8: Both clemastine and CW3388 treatments upregulate the oligodendrocyte gene set in C266 mouse glioma cells; Table S1: List of primers and shRNA sequences; Table S2: List of plasmids, chemical compounds, reagents, kits, and equipment; Table S3: List of antibodies used for immunoblotting and immunofluorescence staining; Table S4: Differentially expressed genes of clemastine or CW3388-treated C266 mouse glioma cells; Table S5: Clemastine affects multiple signaling pathways in C266 mouse glioma cells; Table S6: Differentially expressed genes of clemastinetreated patient-derived BTIC cultures; Table S7: Clemastine alters multiple signaling pathways in patient-derived BTIC cultures; Table S8: List of over-the-counter (OTC) medications repurposed for GBM treatment and entered clinical trials; File S1: Original images for blots and gels.Informed Consent Statement: Informed consent was obtained from all subjects involved in the study.

Figure 2 .
Figure 2. Clemastine attenuates the stemness and progenitor cell features of PDGFRA + BTICs.(A) mRNA expression levels of genes associated with NSCs, OPCs, oligodendrocytes, and suppressors of OPC differentiation in BTIC#102 cells treated with vehicle or clemastine (4 µM) for 9-15 days as assessed by quantitative RT-PCR.n = 7 per gene except MOG, CSPG5, and NOTCH1, n = 6, NCAN, n = 5, and NES, SOX2, n = 3. (B,C) Protein levels of NSC and OPC markers in BTIC#102 cells treated with clemastine at indicated doses for 13-14 days were assessed by (B) immunoblot assays (representative results from three independent experiments) or (C) immunofluorescence (IF) staining (representative images from two independent experiments; 20×, scale bar: 50 µm).(D-G) IF staining of PDGFRA and GALC in BTIC#102 cells treated with vehicle or clemastine (Clem; at 6 µM) for 26 days.(D) Representative IF images (20×, scale bar: 50 µm), (E) quantification of mean PDGFRA and GALC fluorescence intensities, and (F) percentages of PDGFRA + and GALC + cells of the vehicle and clemastine-treated groups.(E,F) Each data point was calculated from a different field of view.n = 22 fields of view for all vehicle groups, and n = 32 for all clemastine-treated groups.(G) The contingency table of PDGFRA -/+ and/or GALC -/+ cell counts to assess the mutual exclusivity of PDGFRA and GALC positivity.A total of 5804 cells were identified from 54 fields of view.(H,I) MA plots summarizing differential mRNA expression of top 40 (H) OPC-specific genes or (I) oligodendrocytes-specific genes between clemastine versus vehicle-treated (15-day treatment) BTIC#102 cells.Larger symbols indicate genes with adjusted p-values < 0.05.(J) Extreme limiting dilution analysis (ELDA) of BTIC#102 cells treated with clemastine (Clem) at indicated doses for 14 days to assess their renewal capacity.p-values indicate significance between vehicle and clemastine-treated groups.(K) ELDA of BTIC#102 (left panel) and BTIC#148 (right panel) cells with or without 15-day clemastine (4 µM) pre-treatments (pre-Tx) and/or subsequent two-week clemastine (4 µM) treatments (cont.Tx).p-values indicate significance between the pre-Tx -cont.Tx -group and other groups.n = 4 per condition."-": no clemastine in the media; "+" with clemastine in the media.Data are represented as (A) geometric mean ± geometric S.D. of fold change relative to vehicle-treated groups or (E,F) mean ± S.E.M. (E,F) Significance was calculated using unpaired t tests (F left panel with Welch's correction) and represented as * p < 0.05, ** p < 0.01, *** p < 0.001, ns: not significant.

Figure 2 .
Figure 2. Clemastine attenuates the stemness and progenitor cell features of PDGFRA + BTICs.(A) mRNA expression levels of genes associated with NSCs, OPCs, oligodendrocytes, and suppressors of OPC differentiation in BTIC#102 cells treated with vehicle or clemastine (4 µM) for 9-15 days as assessed by quantitative RT-PCR.n = 7 per gene except MOG, CSPG5, and NOTCH1, n = 6, NCAN, n = 5, and NES, SOX2, n = 3. (B,C) Protein levels of NSC and OPC markers in BTIC#102 cells treated with clemastine at indicated doses for 13-14 days were assessed by (B) immunoblot assays (representative results from three independent experiments) or (C) immunofluorescence (IF) staining (representative images from two independent experiments; 20×, scale bar: 50 µm).(D-G) IF staining of PDGFRA and GALC in BTIC#102 cells treated with vehicle or clemastine (Clem; at 6 µM) for 26 days.(D) Representative IF images (20×, scale bar: 50 µm), (E) quantification of mean PDGFRA and GALC fluorescence intensities, and (F) percentages of PDGFRA + and GALC + cells of the vehicle and clemastine-treated groups.(E,F) Each data point was calculated from a different field of view.n = 22 fields of view for all vehicle groups, and n = 32 for all clemastine-treated groups.(G) The contingency table of PDGFRA −/+ and/or GALC −/+ cell counts to assess the mutual exclusivity of PDGFRA and GALC positivity.A total of 5804 cells were identified from 54 fields of view.(H,I) MA plots summarizing differential mRNA expression of top 40 (H) OPC-specific genes or (I) oligodendrocytes-specific genes between clemastine versus vehicle-treated (15-day treatment) BTIC#102 cells.Larger symbols indicate genes with adjusted p-values < 0.05.(J) Extreme limiting dilution analysis (ELDA) of BTIC#102 cells treated with clemastine (Clem) at indicated doses for 14 days to assess their renewal capacity.p-values indicate significance between vehicle and clemastine-treated groups.(K) ELDA of BTIC#102 (left panel) and BTIC#148 (right panel) cells with or without 15-day clemastine (4 µM) pre-treatments (pre-Tx) and/or subsequent two-week clemastine (4 µM) treatments (cont.Tx).p-values indicate significance between the pre-Tx − cont.Tx − group and other groups.n = 4 per condition."-": no clemastine in the media; "+" with clemastine in the media.Data are represented as (A) geometric mean ± geometric S.D. of fold change relative to vehicle-treated groups or (E,F) mean ± S.E.M. (E,F) Significance was calculated using unpaired t tests (F left panel with Welch's correction) and represented as * p < 0.05, ** p < 0.01, *** p < 0.001, ns: not significant.

Figure 4 .Figure 4 .
Figure 4. Genetic knockdown and small-molecule inhibition of Ebp suppress the growth of C26 mouse glioma cells in vitro and in vivo.(A) Proliferation of non-targeting shRNA control (shCtl) o Ebp shRNA (shEbp) C266 cells 12 days after lentiviral transduction.p-values indicate significanc between the two groups at indicated timepoints.n = 18 per condition.(B) ELDA of shCtl and shEb C266 cells to assess their self-renewal capacity.(C) Proliferation of C266 cells treated with CW338 at indicated doses.n = 6 per condition except vehicle, n = 12.(D) Proliferation of C266 cells treate with clemastine (Clem) at indicated doses.n = 6 per condition except vehicle, n = 12.(E) Quantifica tion of relative cell proliferation of shCtl and shEbp C266 cells treated with clemastine (left panel) o