Ion channels represent a large class of drug targets, but their role in brain cancer is underexplored. Here, we identify that chloride intracellular channel 1 (CLIC1) is overexpressed in human central nervous system malignancies, including medulloblastoma, a common pediatric brain cancer. While global knockout does not overtly affect mouse development, genetic deletion of CLIC1 suppresses medulloblastoma growth in xenograft and genetically engineered mouse models. Mechanistically, CLIC1 enriches to the plasma membrane during mitosis and cooperates with potassium channel EAG2 at lipid rafts to regulate cell volume homeostasis. CLIC1 deficiency is associated with elevation of cell/nuclear volume ratio, uncoupling between RNA biosynthesis and cell size increase, and activation of the p38 MAPK pathway that suppresses proliferation. Concurrent knockdown of CLIC1/EAG2 and their evolutionarily conserved channels synergistically suppressed the growth of human medulloblastoma cells and Drosophila melanogaster brain tumors, respectively. These findings establish CLIC1 as a molecular dependency in rapidly dividing medulloblastoma cells, provide insights into the mechanism by which CLIC1 regulates tumorigenesis, and reveal that targeting CLIC1 and its functionally cooperative potassium channel is a disease-intervention strategy.
Brain tumors are the leading cause of cancer-related deaths in children and young adults. As a high-grade brain tumor at the hindbrain, medulloblastoma (MB) is one of the most common pediatric malignant brain tumors. Conventional treatment, which includes surgery, chemotherapy, and radiation therapy, is toxic and produces lifelong side effects such as reduced IQ, growth failure, lowered physical performance, and early aging (Edelstein et al., 2011; Frange et al., 2009). MB is composed of four distinct subgroups (WNT, SHH, group 3, and group 4) and further subtypes within each subgroup, which display distinct molecular profiles and clinical features (Cavalli et al., 2017; Taylor et al., 2012). Targeted therapy, which inhibits the SHH pathway receptor Smoothened (SMO), has been developed for SHH pathway–driven MB (SHH MB). After displaying initial efficacy, tumors acquire drug resistance, and relapse is almost always fatal (Rudin et al., 2009). These observations highlight the effectiveness of molecularly targeted therapy and the need to identify novel targets for developing combinatorial therapy. As mouse SHH MB arises from the cerebellar granule neuron precursors (CGNPs), loss of one allele of the SHH pathway inhibitor Ptch1 in Ptch1+/− mice results in 20% MB occurrence (Goodrich et al., 1997). Math1-Cre; SmoM2 mice display constitutive activation of the SHH pathway in CGNPs due to expression of SmoM2 (the constitutively active mutant form of Smo), driven by the CGNP-specific driver Math1-Cre. As a result, Math1-Cre; SmoM2 mice develop fully penetrant MBs (Schüller et al., 2008). Studies using the genetic mouse models of SHH MB and xenograft models of human MB can identify new disease mechanisms and therapeutic targets.
Ion channels are pore-forming, transmembrane proteins that regulate biological processes by controlling ion passage across cell membranes (Hille, 2001). The opening of ion channel pores allows the flux of ions, including potassium, chloride, calcium, or sodium, based on their electrochemical gradient. Ion channels constitute a large class of drug targets for human diseases, such as neurological and cardiovascular disorders (Clare, 2010). However, ion channel function in cancer is underexplored, and its role in pediatric brain tumors was unknown before our studies. We reported that potassium channel EAG2 is up-regulated in ∼15% of human MB across molecular subgroups. Genetic deletion of EAG2 suppressed MB growth in preclinical mouse models (Huang et al., 2012). We identified the US Food and Drug Administration–approved anti-psychotic drug thioridazine as an EAG2 blocker and demonstrated its anti-MB efficacy in mice. We treated a patient with SHH MB, which was resistant to the standard chemo- and radiation-therapy, using thioridazine. The positron emission tomography–computed tomography imaging revealed marked reduction of his tumor, demonstrating a response to the thioridazine therapy (Huang et al., 2015). Therefore, we identified MB dependency on overexpressed ion channels that can be therapeutically targeted.
Cell volume regulation is fundamental to many cellular behaviors, such as proliferation, apoptosis, and migration. Ionic flux across the plasma membrane serves as a mechanism to control intracellular osmolarity, the movement of nonprotein-bound water molecules, and cell volume increase or decrease. Mammalian cell volume undergoes stereotypical oscillations during cell cycle progression. Live imaging studies have shown that cells increase in volume at interphase, reduce volume before mitotic entry, and reach a minimal volume at metaphase, after which the cell volume increases during anaphase and telophase to facilitate cytokinesis (Habela and Sontheimer, 2007; Boucrot and Kirchhausen, 2008). Importantly, a recent study using budding yeast and human fibroblasts demonstrated that while cells scale up protein and RNA biosynthesis in accordance to cell growth and volume increase, excessive cell size beyond a certain level results in increased cytoplasm/DNA ratio (cytoplasm dilution). Enlarged cells undergo cell cycle arrest due to their inability to proportionally match macromolecule biosynthesis to their increased cell size (Neurohr et al., 2019). EAG2 localizes to the plasma membrane during late G2 phase and mitosis to promote potassium efflux and cell volume decrease to ensure mitotic entry and successful mitosis in MB (Huang et al., 2012). As cation and anion channels often function together to regulate cellular behaviors, we reason that a hitherto-unknown anionic channel functionally couples with EAG2 to maintain electroneutrality and coregulate MB malignancy through controlling cell volume homeostasis.
Chloride is a major anion and electrolyte that contributes to the osmotic pressure across the plasma membrane. Chloride intracellular channel 1 (CLIC1) belongs to a family of six CLICs (CLIC1–CLIC6). CLIC1 can exist as either a soluble cytosolic or a membrane protein, and the ion channel activity of CLIC1 has been shown by multiple groups (Littler et al., 2010). CLIC1 channel activity is increased by oxidation, low pH (Tulk et al., 2002; Warton et al., 2002), and cholesterol-rich lipid rafts (Hossain et al., 2016; Valenzuela et al., 2013). CLIC1 permeates several anions (e.g., Cl−, SCN−, F−, and HCO3−) in overexpression systems (Singh and Ashley, 2006; Valenzuela et al., 1997) and forms a chloride channel in microglia cells (Novarino et al., 2004). CLIC1 can localize to intracellular organelles such as phagosomes in macrophage and dendritic cells, where it is thought to increase phagosome chloride concentration to counterbalance the charge difference induced by proton entry (Jiang et al., 2012; Salao et al., 2016). Given its presence as a soluble protein in the cytosol, phylogenetic prediction suggests that CLIC1 may have enzymatic function (Littler et al., 2010), and soluble CLIC1 displays glutathione-dependent oxidoreductase activity in a cell-free system (Al Khamici et al., 2015). Clic1 knockout mice were viable and fertile and appeared normal, aside from showing a mild defect in platelets and macrophages (Qiu et al., 2010; Ulmasov et al., 2017).
CLIC1 overexpression has been reported in cancers of the lung, bladder, gastrointestinal tract, and glioma (Chen et al., 2007; Ding et al., 2015; Jia et al., 2016; Lu et al., 2015; Wang et al., 2011; Zhang et al., 2013). CLIC1 knockdown decreased glioblastoma (GBM) cell growth in vitro and xenograft tumor growth in mice (Setti et al., 2013). While CLIC1 knockdown inhibited cancer cell growth in most contexts, CLIC1 deficiency was reported to increase the proliferation of gastric cancer cells (Ma et al., 2012). In Xenopus laevis, CLIC1 in nontransformed cells exerted long-range suppression of mutant KRAS-induced tumorigenesis (Chernet and Levin, 2014). Therefore, different studies have implicated CLIC1 as either an oncogene or a tumor suppressor, but no direct mechanism for either role has been reported. Given what is known about CLIC1 in cancer, there are multiple knowledge gaps. First, how widespread is CLIC1 overexpression in human central nervous system (CNS) malignancies? Second, what is the in vivo role of CLIC1 in tumor initiation and progression (which cannot be fully addressed by in vitro or xenograft studies)? Third, what is the specific cell type, within the heterogeneous tumor cellular hierarchy, to which CLIC1 confers malignant traits to? Fourth, what is the mechanism by which CLIC1 regulates tumor growth? Addressing these knowledge gaps is key to developing a CLIC1-based rational strategy to combat malignant tumors. Here, we show that CLIC1 is overexpressed in multiple brain cancer types, including MB. Clic1 knockout does not overly affect mouse development but suppresses MB growth in xenograft and genetically engineered mouse models. We show that CLIC1 regulates the proliferation of rapidly dividing MB cells by functionally coupling with potassium channels to control cell volume homeostasis and mitosis.
CLIC1 is overexpressed in multiple human brain cancer types, including MB
To search for anion channels that regulate brain tumor malignancy, we performed bioinformatic interrogation of Gene Expression Omnibus (GEO) datasets. CLIC1 expression is up-regulated in multiple types of primary brain tumors, including MB, GBM, pilocytic astrocytoma (grade I glioma), ependymoma, atypical teratoid rhabdoid tumor (AT/RT), and primitive neuroectodermal tumor (PNET; Fig. 1 A). Within MB, SHH MB, and WNT MB displayed further CLIC1 up-regulation compared with group 3 and group 4 MB (Fig. 1 B). In a human MB tissue microarray, which included primary MB subgrouped based on DNA methylation (Cavalli et al., 2017) and normal cerebellum and cortex tissues as controls, we detected high CLIC1 protein expression in ∼58% (11 of 19) of the patient tumors, while CLIC1 was undetectable in the nontumor controls (Fig. 1 C). CLIC1 overexpression is recapitulated in the tumors of Ptch1+/− and Math1-Cre; SmoM2 mouse models of SHH MB, as well as group 3 and 4 xenograft MB tumors (Fig. 1 D). Quantitative PCR (qPCR) analysis showed high Clic1 expression in MB and early postnatal mouse cerebella compared with cerebella in adult mice (Fig. 1 E). Recent genomic studies showed that genetic alterations in primary MB tumors poorly overlap with those in metastatic tumors, highlighting the need to validate whether a candidate therapeutic target remains in the tumor after metastasis (Wu et al., 2012). We found that high Clic1 and CLIC1 expression is maintained in the metastatic MB tumors in Ptch1+/−; Math1-SB11; T2Onc2 mice and group 3 and 4 MB xenografts, respectively (Fig. 1 F). Collectively, CLIC1 is overexpressed in a wide spectrum of human brain tumors, including MB, and CLIC1 overexpression is recapitulated in mouse MB models amenable to further functional investigation.
CLIC1 knockdown suppresses the growth of MB tumors in multiple xenograft models
To demonstrate the translatability of targeting CLIC1 in human tumors, we investigated the function of CLIC1 using human MB cell lines and xenograft models. shRNA-mediated CLIC1 knockdown decreased CLIC1 protein levels (Fig. 2 A and,Fig. 4 A) and reduced the viability and clonogenic potential of serum-cultured ONS76, DAOY, and MED8A cells grown as monolayers (Fig. 2, B and C). CLIC1 knockdown also suppressed the growth of SOX2+; NESTIN+ spheres of D425 and MB002 human MB cells cultured using defined growth factors (leukemia inhibitory factor, epidermal growth factor [EGF], basic fibroblast growth factor [bFGF]) in serum-free media (Fig. 2, B and C; and Fig. S1 A). CLIC1 knockdown did not impact the growth of human astrocytes or human embryonic kidney 293T cells (Fig. 2 B), suggesting that CLIC1 deficiency reduces cell viability in a tumor-specific manner. Importantly, CLIC1 knockdown suppressed the in vivo growth of MB tumors and extended the survival of tumor-bearing mice in multiple xenograft models (Fig. 2, D and E). These human cell line and xenograft studies demonstrate a cell-autonomous function of CLIC1 in MB across different molecular subgroups. Interestingly, in time-matched or end-point tumors, CLIC1 expression levels were comparable or appeared higher in CLIC1 shRNA-treated tumors (Fig. S1, B and C). Similarly, the histology and percentages of mitotic cells were comparable between these tumors (Fig. S1, B–D). We postulate that the xenograft tumors that eventually displayed in vivo growth were likely from “escaper cells” with insufficient or no perturbation of CLIC1 expression.
Clic1 knockout does not overtly affect mouse development but inhibits the proliferation of tumor cells, thereby extending the survival of SHH MB–bearing mice
Next, we asked what role Clic1 plays in mouse development and in genetically engineered mouse models of MB. We detected Clic1 in mouse SHH MB cells in vivo at various cell cycle phases (Fig. 3 A). We used CRISPR-Cas9 to generate a Clic1 knockout mouse line by introducing a 1,458-bp deletion between introns 1 and 5, which results in the absence of Clic1 protein (Fig. S2, A and B). Clic1−/− mice were born at Mendelian ratio and displayed normal survival within the observed time frame of 120 d, as well as normal body weight and fertility (Fig. S2, C–F). Adult Clic1+/+ and Clic1−/− mice displayed comparable brain morphology (Fig. S2 G). As CGNPs are the cell of origin of mouse SHH MB, we compared the proliferative capacity of CGNPs of Clic1+/+ and Clic1−/− mice at postnatal day 7 (P7), a stage when CGNPs are rapidly dividing at the external granule neuron layer. The percentages of cycling and mitotic CGNPs are comparable between Clic1+/+ and Clic1−/− mice (Fig. S2, H and I). Collectively, these data show that Clic1 is dispensable for normal mouse development.
We next studied the Math1-Cre; SmoM2 mouse model of SHH MB with and without genetic deletion of Clic1. Remarkably, in contrast to the lack of overt phenotype during development, Clic1 knockout reduced MB tumor burden, suppressed tumor-associated hydrocephalus and cranium bulging, and significantly extended the survival of tumor-bearing mice (Fig. 3, B and C). MB in Math1-Cre; SmoM2; Clic1−/− mice displayed reduced proliferation and increased apoptosis compared with the tumors in Math1-Cre; SmoM2 mice (Fig. 3 D). A lineage-tracing study has shown that mouse SHH MB harbors a rare population of slow-cycling, Sox2+ brain tumor–initiating cells. The growth of SHH MB mirrors a dysregulated neurogenic program in that the Sox2+ cells produce Dcx+, transient, rapidly cycling tumor progenitor cells that subsequently generate postmitotic NeuN+ cells. Dcx+ and NeuN+ cells together comprise the tumor bulk (Vanner et al., 2014). We found that Math1-Cre; SmoM2; Clic1−/− tumors displayed significantly decreased mitosis in Dcx+ cells (Fig. 3 E), whereas the number or percentage of mitotic Sox2+ cells was not altered by Clic1 knockout (Fig. 3 F). Furthermore, we isolated and cultured Sox2+ cells from primary mouse MB tumors and determined that Clic1 knockout did not overtly alter the overall growth or cell diameters (Fig. 3 G). Taken together, these data show that Clic1 regulates rapidly dividing tumor cells, but not slow-cycling Sox2+ cells, in SHH MB.
CLIC1 regulates cell volume homeostasis for cell cycle progression
We next investigated the mechanism by which CLIC1 regulates MB cell cycle progression. To achieve high resolution and define the subcellular localization of CLIC1, we performed CLIC1 immunostaining followed by stimulated emission depletion (STED) superresolution imaging of human ONS76 MB cells. Differing from the potassium channels EAG2 or KCNT2 that only localizes at the plasma membrane during mitosis (Huang et al., 2012, 2015), CLIC1 is present at the plasma membrane during both interphase and mitosis (Fig. 4 A). Furthermore, we performed total internal reflection fluorescence (TIRF) imaging to unequivocally visualize the ion channel proteins localized on the plasma membrane. Lending support to our STED imaging data, TIRF imaging showed that CLIC1 is present on the plasma membrane during both interphase and mitosis, whereas EAG2 and KCNT2 displayed mitosis-specific localization on the plasma membrane (Fig. 4 B). To determine CLIC1 ion channel activity, we performed perforated whole-cell patch clamp recording. We added the antibiotic gramicidin to the recording solution, which forms pores on the plasma membrane that are exclusively permeable to monovalent cations and small uncharged molecules. Gramicidin perforated-patch recording has the advantage of not altering intracellular chloride concentration and is suitable for the study of anionic channels when intracellular chloride is a variable of interest (Kyrozis and Reichling, 1995). Mitotic cells were identified based on their mitotic rounding morphology and the presence of condensed chromosomes aligned at the metaphase plate (Fig. 4 C). Currents were elicited by voltage steps from −80 mV to 80 mV in 20-mV increments. Noticeably, mitotic MB cells displayed a significantly higher level of total currents compared with cells at interphase, while CLIC1 knockdown decreased the total currents during both interphase and mitosis (Fig. 4, C and D). It is important to note that we recorded total currents carried by all endogenous ion channels, including but not limited to CLIC1. The recording solutions (extracellular and internal) included potassium and chloride as the electrolytes, and the 5 mM tetraethylammonium chloride in the extracellular solution may block a substantial portion of, but not all, potassium channels. Therefore, the recorded whole-cell currents represent an aggregate of both CLIC1-mediated chloride and potassium channel–mediated potassium currents. The reduction of total currents in CLIC1 knockdown cells suggests that CLIC1 channel activity is prominent at both interphase and mitosis.
As we reported that EAG2-mediated potassium efflux is necessary for cell volume reduction and mitotic entry (Huang et al., 2012), CLIC1-mediated chloride efflux may cooperate with potassium to control cell volume and cell cycle progression. Indeed, CLIC1 knockdown significantly suppressed mitosis (Fig. 4, E and F), accompanied by an increase in cell volume, which was determined by both confocal imaging–based 3D cell volume reconstruction (Fig. 4 G) and Coulter Counter–based cell diameter measurement (Fig. 4 H). Since CLIC1, EAG2, and KCNT2 channels are present on the plasma membrane during mitosis, we asked whether these three ion channels cooperatively regulate cell volume. While ectopically expressing CLIC1, EAG2, or KCNT2 alone was sufficient to reduce the volume of COS7 cells, coexpressing two (CLIC1 and EAG2, CLIC1, and KCNT2) or three ion channels all together synergistically decreased cell volume (Fig. 4 I).
To further investigate the relationship between cell volume and cell cycle progression, we manipulated MB cell volume using various chemical compounds, including capsazepine, an activator of epithelial sodium channels, to increase sodium influx, niflumic acid to inhibit general chloride efflux, and two cholesterol synthesis inhibitors, atorvastin and lovastatin, to suppress lipid raft formation. These pharmacological treatments resulted in a dose-dependent increase of MB cell volume (Fig. S3 A). Notably, by analyzing these same cell populations, we detected robust positive Pearson’s correlation between cell volume and the percentage of MB cells at G1/G0 phase and robust negative Pearson’s correlation between cell volume and the percentage of MB cells at G2/M phase (Fig. S3 B). Given that we manipulated MB cell volume using distinct pharmacological or genetic approaches, which all resulted in defects in mitotic entry, we suggest MB cell cycle progression critically depends on dynamic cell volume controls mediated by ion channel activity. Plasma membrane–localized CLIC1 conducts chloride current, CLIC1 channel activity increases during mitosis, and its deficiency results in cell swelling and defective cell cycle progression.
CLIC1 deficiency activates p38 MAPK pathway to suppress MB cell growth and is associated with uncoupling between RNA biosynthesis and cell size increase
EAG2-deficient MB cells display abnormal cell volume increase, which activates the osmolarity-sensing p38 MAPK pathway to prevent mitotic entry (Huang et al., 2012). Since CLIC1 knockdown resulted in mitotic defects and cell swelling, we asked whether the p38 MAPK pathway is similarly activated. Indeed, CLIC1 knockdown induced the phosphorylation of p38 and its upstream kinases MKK3/6, whereas the phosphorylation level of ERK1/2 appeared unaltered (Fig. 5 A). We analyzed RNA sequencing (RNA-seq) data of nontargeting and CLIC1 shRNA-treated MB cells and identified enrichment of p38 MAPK pathway target genes upon CLIC1 knockdown (Fig. 5, A and B; and Fig. S3 C). Furthermore, we treated three MB cell lines (ONS76, MED8A, and DAOY) with nontargeting shRNA and CLIC1 shRNA and suppressed p38 MAPK activity using two structurally different kinase inhibitors (SB 203580 and SB 239063). While p38 inhibitors did not affect growth of the nontargeting shRNA treated cells, adding either p38 inhibitor profoundly rescued the growth defect caused by CLIC1 knockdown in all three MB cell lines (Fig. 5 C and,Fig. S3 E). Interestingly, while p38 inhibition did not affect cell volume in control or nontargeting shRNA-treated MB cells, it partially suppressed the cell volume increase phenotype upon CLIC1 knockdown (Fig. S3 F). It is possible that inhibiting p38 activity in the presence of CLIC1 knockdown reduces cell volume through bypassing cell cycle arrest and allowing cell division. Alternatively, p38 activation may exacerbate the cell volume increase phenotype in addition to the cell swelling caused by deficient CLIC1 chloride efflux. Taken together, these data establish that CLIC1 deficiency activates the p38 MAPK pathway, which exacerbates cell swelling, and that activation of the cell volume–sensing p38 MAPK pathway is a shared mechanism that deficiency in either EAG2 or CLIC1 engages to suppress MB cell growth.
Cells scale up protein and nucleic acid biosynthesis in accordance to cell volume increase. However, aberrant cell size increase past a certain threshold results in an elevated cytoplasm/nucleus ratio, uncouples macromolecule biosynthesis, and impairs the cell cycle (Neurohr et al., 2019). We asked whether CLIC1 regulates cell volume to prevent cytoplasm dilution during cell cycle progression. To this end, we quantified the nuclear and overall cell volume. Strikingly, while the nuclear volume was not altered, CLIC1-deficient cells displayed significantly increased cell/nuclear volume ratio, suggesting that CLIC1 specifically regulates the volume of the cytoplasm (Fig. 5 D and,Fig. S3 D). We performed RNA-seq and bioinformatic analysis to determine the CLIC1-dependent gene expression and signaling pathways. Gene set enrichment analysis (GSEA) identified pathways that regulate cell cycle, cell size, and cell growth (Fig. 5 E). More importantly, pathways of ribosome, RNA polymerase, Spt-Ada-Gcn5 acetyltransferase (SAGA) complex, DNA repair, proteasome, and mitochondria, all of which have been functionally linked to the oversized cells and deficient biosynthesis (Neurohr et al., 2019), are among the most altered after CLIC1 knockdown (Fig. 5 E and,Fig. S4 A). Lastly, we compared the global RNA quantity of the control and CLIC1 knockdown cells. Despite an ∼1.8-fold increase in cell volume (Fig. 4 G), CLIC1 knockdown MB cells express significantly decreased total RNA compared with control cells (93.6 ± 0.6% vs. 100 ± 0.8%, P = 0.004; Fig. 5 F), suggesting that their total RNA was diluted almost twofold relative to the control. These results suggest that CLIC1 regulates the cytoplasm/nuclear volume ratio and prevents cytoplasm dilution and indicate a role for CLIC1 in scaling RNA levels to MB cell size during cell cycle progression.
CLIC1 and EAG2 are synthetic vulnerabilities in MB cells
CLIC1 undergoes palmitoylation, a posttranslational modification that target proteins to lipid rafts (Fang et al., 2016). Cholesterol promotes CLIC1 membrane insertion and ion channel activity (Hossain et al., 2016). Since lipid rafts are cholesterol-rich microdomains on the plasma membrane, we hypothesize that lipid rafts may represent designated membrane sites where CLIC1 and EAG2 colocalize and functionally cooperate in MB cells. To study the dynamics of EAG2 and CLIC1 localization during mitosis, we first quantified EAG2 puncta and CLIC1 puncta localized at lipid rafts, and colocalization of EAG2 and CLIC1 (which we define as being within 0.2 µm apart) at lipid rafts on the plasma membrane and nonlipid raft plasma membrane, during both interphase and mitosis. While the extent of either EAG2 or CLIC1 localization at lipid rafts did not differ between interphase and mitosis (Fig. 6 A and,Fig. S4 B), colocalization of EAG2 and CLIC1 on lipid rafts significantly increased at mitosis, supporting the notion that these two ion channels increase physical proximity at lipid rafts to promote potassium and chloride efflux for mitosis. We did not detect appreciable CLIC1 and KCNT2 colocalization at lipid rafts during mitosis (Fig. S4 C).
Next, we asked whether CLIC1 and EAG2 channels represent synthetic vulnerabilities in MB cells. To this end, we tested multiple concentrations of lentivirus-expressing shRNAs against CLIC1 or EAG2, either alone or in combination. We identified a specific lentiviral concentration at which single-ion-channel knockdown did not elicit discernible phenotype. However, combinatorial knockdown of CLIC1 and EAG2 induced synergistic and significant suppression of the growth of human MB cells (Fig. 6 B and,Fig. S4 D). CLIC1 and EAG2 also displayed synergistic regulation of cell volume (Fig. 6 C), supporting a mechanism whereby CLIC1 and EAG2 cooperate to regulate MB cell growth by converging on cell volume control. Either single or double knockdown activated the p38 MAPK pathway (Fig. S4 E). Interestingly, we noted that p38 activation was not enhanced by comparing double- and single-channel knockdown. Since CLIC1 or EAG2 deficiency not only activated the p38 pathway but also resulted in aberrant cell volume increase and reduction of global mRNA density, a probable explanation for this is that CLIC1 and EAG2 synergize to regulate MB cell cycle progression via both p38-dependent and p38-independent mechanisms. Furthermore, these data suggest that different subgroups of MB tumors may evolve a common mechanism by which CLIC1 cooperates with potassium channels for cell cycle control. In support of this notion, our bioinformatics analysis revealed that 10 potassium channels, including KCNT2, display the highest positive correlation with CLIC1 across human MB subgroups (Fig. S5 A).
Drosophila melanogaster Eag and Clic channels cooperate to regulate brain tumor growth in vivo
Drosophila represents a powerful genetically tractable model to study brain tumor biology. During Drosophila larval brain development, neuroblasts (neural stem cells) proliferate to self-renew and generate neural progenitor cells, neurons, and glia. Timely cell cycle entry and exit is key to brain development, and aberrant cell cycle progression leads to brain tumorigenesis (Homem and Knoblich, 2012). insc-Gal4 is a genetic driver active in the proliferative neuroblast and progenitor cells. insc-Gal4-driven expression of the basic helix-loop-helix oncogene Dpn results in uncontrolled neuroblast self-renewal (Zhu et al., 2012), which mimics neural stem cell– or progenitor-derived human brain tumors, including MB. Using this model, we previously identified the potassium channel Eag as a critical regulator of the growth of Drosophila brain tumors (Huang et al., 2015). Here, we first asked whether the Drosophila chloride intracellular channel (Clic) plays an evolutionarily conserved role in governing the proliferation of brain tumor cells. A previous study has identified two hypomorphic alleles of Clic, ClicEY04209 and ClicG0472, in which transposon insertions abrogated transcript expression (Bhandari et al., 2012). As Clic is located on the X chromosome, we generated male brain tumor–bearing larvae with either one of the two Clic mutations. We used the insc-Gal4 to express mCD8-GFP that labels cell membranes and performed confocal imaging to quantify the volume of tumor tissues. We also performed immunostaining to quantify the numbers of neuroblasts and neural progenitor cells (Dpn+) and glia cells (Repo+). Both hypomorphic mutations reduced the number of Dpn+; pHis3+ mitotic neuroblast and neural progenitor cells and decreased the volume of mCD8-GFP+ tumor tissue (Fig. 6 D). As a control, global Clic deficiency in the ClicEY04209 or ClicG0472 mutants (without brain tumors) did not result in consistent alteration of the numbers of Dpn+ or Repo+ cells of the brain lobes and ventral nerve cord (VNC; Fig. S5 B), highlighting a heightened dependency on Clic in brain tumors compared with normal brains.
Having established the essential role of Clic, we then modeled our combinatorial RNA interference (RNAi) experiments in human MB cells and asked whether Clic and Eag functionally cooperate to regulate Drosophila brain tumor growth. We selected specific RNAi lines such that single Clic or Eag RNAi was unable to cause tumor suppression, likely due to insufficient gene knockdown. Notably, we observed robust synergism when we combined the RNAi of both Clic and Eag (Fig. 6 E). As a control, either single or double ion channel knockdown did not elicit difference in the number of Dpn+ neuroblast and progenitor cells or the volume of insc-Gal4–driven mCD8-GFP–labeled neuroblast lineage in the nontumor-bearing larvae (Fig. S5 C). Collectively, our Drosophila data demonstrate a dependency on Clic in oncogenically transformed cells and establish Clic and potassium channel cooperation as an evolutionarily conserved mechanism that fuels malignant brain tumor growth in vivo.
An integrative cross-species platform to identify ion channels as dependencies of brain tumor
MB is one of the most prevalent childhood brain cancers. As a large class of cell surface proteins and drug targets for a variety of human diseases, ion channels in brain tumors are underexplored. By bioinformatic interrogation of large Cancer Genome Atlas datasets, we found that CLIC1 is overexpressed in multiple human CNS cancer types, including MB (Fig. 1). Using xenograft models, we demonstrated that CLIC1 knockdown suppressed the growth of MBs across molecular subgroups (Fig. 2). By generating Clic1 knockout mice and studying genetically engineered mouse models, we showed that Clic1 is dispensable for normal mouse development, whereas Clic1 knockout suppressed tumor growth in a mouse model of SHH MB (Fig. 3). Mechanistically, we discovered that CLIC1-dependent ion channel activity increases at mitosis (Fig. 4). We showed that CLIC1 deficiency activates the cell volume–sensing p38 MAPK pathway to suppress MB cell growth and uncovered that CLIC1 regulates MB growth through controlling the molecular machineries that converge on regulation of cell size and RNA biosynthesis to prevent cytoplasm dilution (Fig. 5). Furthermore, we identified lipid rafts as the plasma membrane microdomains where CLIC1 and EAG2 enrich during mitosis (Fig. 6) and demonstrated that combinatorial targeting of CLIC1/EAG2 and Clic/Eag suppresses the growth of human MB cells in vitro and Drosophila brain tumors in vivo, respectively (Fig. 6). Our studies highlight the utility of a cross-species platform to identify evolutionarily conserved tumor dependencies, establish CLIC1 and EAG2 channels as synthetic vulnerabilities in MB, and reveal the mechanism of action by which CLIC1 regulates tumorigenesis (Fig. 6 F).
CLIC1 in cancer
We discovered previously unreported CLIC1 overexpression in multiple types of primary brain tumors, including MB, pilocytic astrocytoma, ependymoma, AT/RT, and PNET. This widespread alteration of CLIC1 expression in human tumors indicates that CLIC1 may regulate fundamental aspects of malignant cells. We propose that CLIC1 functionally couples with potassium channels to control cell volume oscillation to facilitate tumor cell cycle progression. Clic1−/− mice are viable and fertile, without overt phenotypes, highlighting the value of CLIC1 as a cancer therapeutic target. We note that CLIC1 may play additional roles in the cytosol or at intracellular organelles to regulate tumor growth. CLIC1 may also function in nontumoral cells, such as macrophage or microglia, to influence tumor progression. It is intriguing that CLIC1 knockdown resulted in stronger growth suppression in human MB xenograft tumors compared with genetically deleting Clic1 in the Math1-Cre; SmoM2 mouse SHH MB. We consider several possible reasons. First, constitutive global Clic1 knockout in mice may allow compensatory mechanisms to reduce tumor dependency on Clic1. Second, MB is a highly heterogeneous disease. The Math1-Cre; SmoM2 mouse model of SHH MB may not represent the human MB tumors we used in the xenograft studies. Third, MB tumors that arise in mouse and human may develop species-specific dependency on Clic1 and CLIC1, respectively.
Functionally coupled ion channels as regulators for tumor growth
Ion channels cooperate to regulate biological processes. For example, activating the calcium-permeating TRPV4 channel results in increased activity of the calcium-activated chloride channel TMEM16A, which controls chloride flux to regulate the volume of choroid plexus epithelial cells. This cooperation between TRPV4 and TMEM16A is facilitated by physical interaction (Takayama et al., 2014). TRPV4 and calcium-activate potassium channel SK3 cocluster in the caveolae of endothelial cells to regulate vasodilation (Goedicke-Fritz et al., 2015). Superresolution imaging demonstrated various classes of ion channels intimately associate with each other to form signaling macromolecular complexes in neurons (Zhang et al., 2016). We speculate that a network of ion channels temporally and spatially enriches to membrane microdomains and functionally couple to regulate cancer cell behavior, including, but not limited to, mitosis. Identifying the components of such ion channel networks will reveal ample new therapeutic opportunities.
Therapeutically targeting CLIC1 in cancer
We identified that the US Food and Drug Administration–approved antipsychotic drug thioridazine is an EAG2 channel blocker, which displayed therapeutic efficacy in treating MB in mice and a patient (Huang et al., 2015). IAA-94, which inhibits CLIC1, decreased the capacity for colorectal cancer cells to regulate volume during cell migration (Wang et al., 2012). However, the high concentration required for IAA-94 to function and its effect on CLIC1–CLIC6 preclude its use as a CLIC1-specific agent for cancer treatment. A specific and potent CLIC1 inhibitor is needed for further therapeutic development. As the structure of soluble CLIC1 is known, rational design of CLIC1-interacting small molecules or peptides may offer new avenues to blocking CLIC1 function. Alternatively, high-throughput screening of agents that either antagonize CLIC1 channel activity or prevent its membrane trafficking may lead to the development of CLIC1-specific targeted therapy. A combination therapy targeting CLIC1 and its functionally coupled potassium channels (e.g., EAG2 in MB) should maximize the tumor-suppressive effect. Leveraging tumor dependency on overexpressed ion channels such as CLIC1 represents a novel approach to effectively treat devastating cancers such as MB.
Materials and methods
The following fly stocks were used: insc-Gal4 (Bloomington stock number 8751), Gal4 inserted in the inscuteable promoter to drive gene expression in the neuroblast lineage; UAS-Dpn (Yuh Nung Jan laboratory, University of California, San Francisco, San Francisco, CA), UAS-driving expression of Dpn; ClicEY04209 (Bloomington stock number 15725), hypomorphic allele of Clic; ClicG0472 (Bloomington stock number 12267), hypomorphic allele of Clic; RNAi-Clic (Vienna Drosophila Resource Center stock number 105975), UAS-driving expression of RNAi against Clic; RNAi-control (Bloomington stock number 35786), UAS-driving RNAi expression against a scrambled sequence; RNAi-Eag (Vienna Drosophila Resource Center stock numbers 31678 and 9126), UAS-driving RNAi expression against Eag; UAS-mCD8-GFP (Bloomington stock number 32189), UAS-driving expression of mCD8-GFP.
Drosophila immunostaining, confocal microscopy, and image quantifications
All larvae and flies were cultured at 25°C, aside from those involving the Clic hypomorphic mutants, which were cultured at room temperature. Randomly selected third-instar larvae were dissected in PBS, fixed in 4% formaldehyde for 20 min at room temperature, and incubated with primary antibody overnight at 4°C and secondary antibody for 2 h at room temperature. Primary antibodies were chicken anti-GFP (1:1,000; Aves Labs), rabbit anti-Dpn (1:500; Yuh Nung Jan laboratory), mouse anti-Repo (1:100; Developmental Studies Hybridoma Bank) and rat (1:200; Abcam) or rabbit (1:500; Abcam) anti-phosphorylated Histone H3 (Ser10). Secondary antibodies conjugated to Alexa Fluor 488, 555, or 647 (Cell Signaling Technologies) were used at 1:500. Images were acquired with a 20× objective on a Leica SP8 confocal microscope or a Quorum spinning disk confocal equipped with a Leica DMi8 microscope and Photometric Prime 95B camera, with a Z-step size of 5 µm. Tissue volume (insc-Gal4-driven mCD8-GFP) or Dpn and Repo was determined by Imaris software (v9.2) after 3D reconstruction of z-stack confocal images using the Surfaces function (with surface detail smoothing = 1 µm, manual thresholding based on absolute intensity, background subtraction, and filtering based on signal intensity mean). Mitotic cells (pHis3+), Dpn, and Repo cells in the nontumor brains were quantified using the Spots quantification function on Imaris (with an estimated diameter = 4 µm, manual thresholding, background subtraction, filtering, and manual deletion of nonspecific spots). To quantify the number of Dpn+; pHis3+ double-positive cells, a MATLAB extension on Imaris (Find Spots close to Surfaces, with a distance threshold determined by the average of manual distance measurements between the spots and surfaces of interest) was used. To assess the survival of adult Drosophila, the resulting progeny classes from each cross were scored. Crosses and scoring were performed three independent times.
Mouse lines, immunohistochemistry, and quantification
Ptch1+/−, Math1-Cre, SmoM2 mice were previously described and cited in the text. Clic1−/− mice were generated at The Centre for Phenogenomics by injecting Cas9 mRNA and two guide RNAs with spacer sequences 5′-TTAGATAGAGGCTTAGTGAA-3′ and 5′-CCAGTACCCAATTAAAAAGC-3′. This resulted in a 1,458-bp deletion in chromosome 17 (35,052,089–35,053,546), which spans from intron 1–2 to intron 4–5 in the Clic1 gene. The common name for this allele is Clic1_em2_del. The official Nomenclature is C57BL / 6N-Clic1em2Tcp. All procedures were performed in compliance with the Animals for Research Act of Ontario and the Guidelines of the Canadian Council on Animal Care. The Centre for Phenogenomics Animal Care Committee reviewed and approved our protocol 19-0288H. All strains were backcrossed at least eight generations, and littermates were used as controls for all experiments.
Immunohistochemistry, immunofluorescence, and H&E staining were performed on paraffin-embedded tissue sections. The primary antibodies were mouse anti-CLIC1 (1:50; Santa Cruz Biotechnologies), rabbit anti-Sox2 (1:100; Abcam), rabbit anti-Dcx (1:100; Abcam), mouse anti-phospho-Histone H3 (1:1,000; Cell Signaling Technologies), rabbit anti-cleaved Caspase 3 (1:200; Cell Signaling Technologies), rabbit anti-Pax6 (1:100; BioLegend), and rabbit anti-Ki67 (1:100; Abcam). Images were acquired with a Leica SP8 Lightning Confocal (DMI6000) microscope. The numbers of Pax6+, Ki67+, Sox2+, pHis3+, or cleaved-caspase3+ cells were quantified using the “Spots” function on Imaris software with an estimated diameter of 2.5–4.5 µm, manual thresholding, background subtraction, and filtering. Dcx volume was determined through the “Surface” function in Imaris after 3D reconstruction of z-stack confocal images (with smoothing level of 0.5 µm, manual thresholding, background subtraction, and filtering based on signal intensity mean). To quantify the number of Dcx+; pHis3+ double-positive cells, a MATLAB extension on Imaris: Find Spots close to Surfaces with a distance threshold determined by the average of manual distance measurements between the spots and surfaces of interest, was used.
The group 3 xenograft MB used for CLIC1 immunohistochemistry in Fig. 1 is Med-411FH, which is a MYC amplified group 3 xenograft obtained from the Brain Tumor Resource Laboratory, and profiled using genome-wide methylation arrays at our institution and reported (Brabetz et al., 2018). The group 4 xenograft MB used for CLIC1 immunohistochemistry in Fig. 1 is DMB006, which was previously reported (Pei et al., 2016).
Mouse Sox2+ MB cells were isolated from Math1-cre; SmoM2 and Math1-cre; SmoM2; Clic1−/− mice and cultured as previously described (Huang et al., 2010). Briefly, primary MB tumors were mechanically dissociated by repetitive pipetting using ice-cold PBS without Mg2+ and Ca2+, followed by treatment using 50% Accutase diluted in PBS. The dissociated cells were cultured on a plate surface coated with poly-L-ornithine and laminin using Neurocult NS-A Basal media supplemented with 2 mmol/liter L-glutamine, N2 supplement, B27 supplement, 75 µg/ml BSA, 2 µg/ml Heparin, 10 ng/ml bFGF, and 10 ng/ml human EGF.
Human MB tissue microarray
Primary MB samples subgrouped using DNA methylation were used to assemble the tissue microarray (Cavalli et al., 2017). Three representative cores per sample were used along with normal cerebellum and cortex tissues as controls. Tissue samples (sections) were deparaffinized and rehydrated through an alcohol gradient to water for antigen retrieval in 10 mM citrate buffer, pH 6.0, in a rice cooker. The mouse anti-CLIC1 antibody (1:50; Santa Cruz Biotechnologies) was incubated overnight at 4°C, and the anti-rabbit HRP step was performed with 1-h incubation at room temperature. The 3′-diaminobenzidine color reaction was developed at room temperature and stopped after 5 min. The slides were then counterstained with hematoxylin. Tissue sections were imaged using a 3DHistech Pannoramic 250 Flash II Slide Scanner.
Cell culture, clonogenic and cell viability assays, immunofluorescence staining, and drug inhibition experiments
DAOY and MED8A human MB cell lines were cultured using DMEM with 10% FBS. The ONS76 human MB cell line was cultured using RPMI with 10% FBS and supplemented with sodium pyruvate. D425 and MB002 human MB cell lines were cultured using 50:50 neural basal medium and DMEM/F12 supplemented with nonessential amino acids, sodium pyruvate, Hepes, glutamine, B27, heparin, 25 ng/ml human leukemia inhibitory factor (LIF), 25 ng/ml bFGF, and 25 ng/ml human EGF without addition of serum. All cell lines were regularly checked for mycoplasma infections and treated with Plasmocin (Invivogen) when infection was noted. No cell lines are listed in the database of commonly misidentified cell lines maintained by the International Cell Line Authentication Committee and National Center for Biotechnology Information BioSample.
To determine the clonogenic potential of brain tumor cells, cells were infected with lentivirus for 4 d and replated at clonal density (150 cells/ml of culture medium) into 60-mm plates. Cell colonies were stained with 0.05% crystal violet 7–10 d after seeding. Representative results from three independent experiments were shown. Cell viabilities were determined using CellTiter 96 AQueous One Solution Cell Proliferation Assay (Promega) after plating (500–1,000 cells per well). For the nonadherent D425 or MB002 cell lines, the cells were plated in serial dilutions (ranging from 2,000 to 3 cells per well) on nonadherent 96-well plates in six biological replicates. 1 wk after plating, the CellTiter 96 AQueous One Solution Cell Proliferation Assay was performed to measure cell viabilities on a spectrophotometer.
To perform the double knockdown (CLIC1 and EAG2), MB cells were plated in 96-well plates at a density of 500–1,000 cells per well. Multiple concentrations of the lentivirus, which encode shRNA targeting CLIC1 or EAG2, either alone or in combinations were added to the wells. Each well received an equal virus titer (in the single knockdown conditions, control shRNA was added equalize viral titer). The relative viabilities were determined using CellTiter 96 AQueous One Solution Cell Proliferation Assay (Promega) after lentivirus treatment as described (in three biological replicates).
To perform immunofluorescence staining, cells were grown on 1-mm coverslips (or 1.5 mm for superresolution and TIRF imaging), fixed in 4% formaldehyde at 4°C overnight, and incubated with primary antibody overnight at 4°C and secondary antibody for 1 h at room temperature. For nonadherent lines, the staining was performed in 0.5-ml tubes with constant nutation. The primary antibodies were chicken anti-GFP (1:2,000; Living Colors), rabbit anti-EAG2 (1:1,000; Alomone), rabbit anti-KCNT2 (1:500; Novus), mouse anti-CLIC1 (1:250; Developmental Studies Hybridoma Bank), rat anti-phospho-Histone H3 (1:10,000; Abcam), rabbit anti-Sox2 (1:1,000; Abcam), and rabbit anti-Nestin (1:500; Abcam). Secondary antibodies conjugated to Alexa Fluor 488, 555, and 647 (Cell Signaling Technologies) were used at 1:1,000, and those conjugated to Alexa Fluor 405, 594, and 532 (Thermo Fisher) were used at 1:2,000. For STED superresolution staining, images were acquired using a Leica SP8 Lightning Confocal (DMI6000) with STED 3X microscope using Hybrid Detector (HyD) detectors and processed using Huygens Professional deconvolution software. The concentrations for the primary and secondary antibodies were increased 5× and 3×, respectively. DAPI or TO-PRO3 was used to visualize DNA. For TIRF imaging, cells were grown and stained as previously described but remained unmounted. TIRF images were acquired using a Zeiss Elyra PS1 (Axio Observer Z1) microscope with an Andor iXon3 885 (SIM) detector and processed using Imaris software.
Lipid rafts were labeled using the Molecular Probes Vybrant Lipid Raft Labeling Kits–Alexa Fluor 555, followed by staining with a goat anti-cholera toxin B (1:200; Fitzgerald Antibodies) to amplify the signal. Confocal and STED superresolution images were acquired with a Leica SP8 Lightning Confocal (DMI6000) with STED 3X microscope using HyD detectors. For superresolution imaging, a 93× (glycerol immersion) objective was used using constant-wave STED 594- and 660-nm depletion lasers with gating. Image deconvolution for the lipid rafts experiment was performed using the Leica Lightning deconvolution software. Image deconvolution for superresolution images was performed using Huygens Profession Software with a STED function.
Analysis and quantification for the lipid rafts experiment was performed using Imaris in a multistep process. To define the localization of EAG2 and CLIC1, we first quantified the immunofluorescence signals defined as spots using Imaris software on lipid rafts (lipid raft label+). The Spots function was used to create a 3D reconstruction and to quantify the number of EAG2, CLIC1, or KCNT2, with background subtraction and filtering using intensity mean. The Surfaces function was used to create a 3D reconstruction of the lipid rafts or membrane GFP. Two MATLAB extensions were used with manual thresholding: Colocalize Spots to quantify the total EAG2 + CLIC1 and KCNT2 + CLIC1 colocalization and Find Spots Close to Surface to determine the number of each EAG2, CLIC1, KCNT2, EAG2 + CLIC1, or KCNT2 + CLIC1 that are on either lipid rafts or membrane GFP. A colocalization channel was built using the “Coloc” function on Imaris to determine the are covered by both lipid rafts (lipid raft label+) and nonlipid raft plasma membrane (membrane GFP+; lipid raft label−), and this was 3D-reconstructed using the Surfaces function. Then MATLAB: Find Spots Close to Surface was used to determine the number of EAG2, CLIC1, KCNT2, EAG2 + CLIC1, or KCNT2 + CLIC1 that are both on lipid rafts and membrane GFP. The relative number of colocalized EAG2 + CLIC1 and KCNT2 + CLIC1 on lipid rafts was determined using the following equation: EAG2 + CLIC1 spots close to lipid rafts + mGFP surface/total EAG2 + CLIC1 spots close to mGFP surface.
To determine cell volume, we expressed soluble GFP in MB cells and performed serial confocal imaging. We acquired the image stacks with a z-step size at 1 µm per optical slice using a 63× 1.4-NA oil lens throughout the entire thickness of the cells. All subsequent image processing was conducted with Imaris 9.2 software and volume was measured using the Surfaces function as previously described above. For cell diameter measurements, cells were plated at a density of 60,000 cells per well in 6-well plates and treated as described for knockdown or pharmacological experiments. Cell diameter measurements were taken 3–4 d after infection or drug treatment. Cells were resuspended using trypsin and incubated in isotonic solution, and diameter was measured using the Multisizer 4 Coulter Counter (Beckman-Coulter) using standard protocols and a threshold of 12–60 µm. All experiments were repeated three times.
p38 MAPK inhibitors SB203580 and SB239063 were obtained from Tocris (catalog numbers 1202 and 1962) and dissolved in DMSO as 10 mM stock. Indicated concentrations of SB203580 and SB239063 were added on the second day post lentivirus infection (following the replacing virus-containing medium with virus-free new medium). Cell viability assays were performed at day 5–10 after infection. Cell diameter measurements were obtained using the Multisizer 4 Coulter Counter (Beckman-Coulter) using standard protocols and a threshold of 12–30 µm.
Drug treatment, cell diameter measurement, and cell cycle analysis by propidium iodide staining
Atorvastatin, lovastatin, capsazepine, and niflumic acid were obtained from Tocris (catalog numbers 3776, 1530, 0464, and 4112, respectively) and dissolved in DMSO to the following stock concentrations: atorvastatin, 10 mM; lovastatin, 10 mM; capsazepine, 30 mM; and niflumic acid, 50 mM. ONS76 MB cells were plated and treated with the indicated concentration of drugs and collected 3 d after for cell diameter measurement and propidium iodide staining. Cell diameter measurements were obtained using the Multisizer 4 Coulter Counter (Beckman-Coulter) using standard protocols and a threshold of 12–30 µm.
For cell cycle analysis, resuspended MB cells were fixed overnight in 80% ethanol at 4°C. Fixed cells were incubated in 2 mg/ml bovine RNase IIA (Sigma) for 5 min at room temperature. Cells were stained in 0.1 mg/ml propidium iodide solution (Life Technologies) in HBSS (Gibco) containing 0.6% Triton X-100 (Bio Basic) for 30 min at room temperature, protected from light. Stained cells were resuspended in HBSS containing 10 mM Hepes·NaOH, pH 7.2, 2% FBS (Gibco), and 10 mM NaN3, then filtered through cell-strainer caps (Falcon) into FACS tubes immediately before flow cytometry measurement. Data were collected on a Becton Dickinson (BD) LSRFortessa cell analyzer, with acquisition software BD FACSDiva v.8.0. Subsequent cell cycle analysis was performed on FlowJo 10 using the Watson Pragmatic algorithm cell cycle platform.
Total protein was extracted using a lysis buffer containing 50 mM Hepes, pH 7.4, 150 mM NaCl, 1% NP-40, 1 mM dithiothreitol, 1 mM EDTA supplemented with Complete Protease Inhibitor Cocktail (Roche). Protein lysate samples, 25ug each, were resolved on 4–12% Bis-Tris gels using 3-(N-morpholino)propanesulfonic acid (MOPS) buffer (Invitrogen) and transferred onto a polyvinylidene fluoride membrane (Amersham Biosciences) using a buffer containing 25 mM Tris base, pH 8.3, 192 mM glycine, and 20% methanol for 2 h at 200 V at 4°C. Western blot assays were performed using the appropriate dilution of primary antibody in TBS supplemented with 0.1% Tween-20 and 0.5% milk or 0.5% BSA. Immunoreactive bands were visualized using HRP-conjugated secondary antibodies (Cell Signaling Technologies), followed by chemiluminescence with ECL-plus Western Blotting Detection System (Amersham). Chemiluminescence was imaged and analyzed using Molecular Imager@VersaDoc MP4000 system (Bio-Rad). The primary antibodies used were mouse anti-CLIC1 (1:500; Santa Cruz Biotechnologies), mouse-anti-αTubulin (1:10,000; Sigma), rabbit-anti-phospho-AKT (Thr308; 1:1,000; Cell Signaling Technologies), rabbit-anti-phospho-p38 (Thr180/Tyr182; 1:1,000; Cell Signaling Technologies), rabbit-anti-phospho-MKK3 (Ser189)/MKK6 (Ser207; 1:1,000; Cell Signaling Technologies), rabbit anti-phospho-ERK1/2 (Thr202/Tyr204; 1:1,000; Cell Signaling Technologies), and rabbit-anti-p38 (1:1,000; Cell Signaling Technologies). Average fold change of protein level from three independent experiments were indicated below the corresponding protein bands.
RNA extraction, reverse transcription, and real-time RT-PCR
For mouse cerebellum and MB tissues, total RNAs were extracted using TRIzol (Thermo Fisher) and cleaned up with a RNeasy Plus Mini kit (Qiagen). 150 ng to 1 µg RNA was reverse transcribed to cDNA using the SensiFAST cDNA synthesis kit (BIOLINE). All DNA and RNA concentrations were measured by a NanoDrop 1000 Spectrophotometer. Real-time detection and quantification of cDNAs were performed with the Viia7 Cycler (Applied Biosystems). qPCR was performed in a 15-µl reaction mixture using SYBR Green qPCR Supermixes (Applied Biosystems). 40 cycles of amplification were performed according to manufacturer’s instructions. Fluorescence data were collected at annealing stages and real-time analysis performed with Viia7 System Software. Ct values were determined with automatically set baseline and manually adjusted fluorescence threshold. Gene expressions were normalized with that of a control housekeeping gene and analyzed using the ΔΔCt method. All experiments were repeated three times.
The following primers were used (h, human; m, mouse): hCLIC1, 5′-TTCAAACCCAGCACTCAATG-3′ (forward), 5′-TCAGCACTGGTTTCATCCAC-3′ (reverse); hGAPDH, 5′-CTCCTGCACCACCAACTGCT-3′ (forward), 5′-GGGCCATCCACAGTCTTCTG-3′ (reverse); mClic1, 5′-ACTCAAACCCAGCCCTCAAT-3′ (forward), 5′-CGCTGGTTTCATCCACTTCT-3′ (reverse); and mGapdh, 5′-TTCACCACCATGGAGAAGGC-3′ (forward), 5′-GGCATGGACTGTGGTCATGA-3′ (reverse).
Total RNAs were extracted from DAOY and ONS76 cells (control and CLIC1-knockdown treated) using TRIzol (Thermo Fisher) and cleaned up with an RNeasy Plus Mini kit (Qiagen). Total RNA (genomic DNA–depleted) samples were sent to The Hospital for Sick Children TCAG Sequencing Facility for poly(A) library prep and sequenced on an Illumina Hiseq 2500 platform with a target of 30 million reads per sample.
Lentivirus-mediated shRNA study
Human pLKO.1 lentiviral shRNA target gene set against CLIC1 and pLKO.1-TRC-control vector were obtained from Open Biosystems. Virus infections were performed within antibiotics-free culture medium for 24 h. The CLIC1 shRNA mature antisense sequences are 5′-TTCAGCACTGGTTTCATCCAC-3′ (#1) and 5′-TCAACGGTGGTAACATTGAAG-3′ (#2). The EAG2 shRNA mature antisense sequence is 5′-ATATTATTCTTGAGTCGCAGC-3′.
Xenograft and in vivo bioluminescence imaging
Tumor cells with firefly luciferase–expressing reporter
MB cells were transduced with a lentiviral vector containing firefly luciferase under the control of the spleen focus-forming virus promoter. Cells were screened for infection efficiency by treatment with luciferin (D-luciferin potassium salt; Gold Biotechnology) in vitro and examination by the Xenogen IVIS Lumina System. More than 95% of cells were infected.
Surgical procedure for implantation of tumor cells
6–8-wk-old female NOD scid gamma/J#5557 immunodeficient mice were used for xenograft experiments. The mice were housed under aseptic conditions, which included filtered air and sterilized food, water, bedding, and cages. Age-matched mice with the same sex (females) were randomly assigned to experimental groups. There were no specific inclusion or exclusion criteria. Mice were anesthetized using gaseous isoflurane and immobilized in a stereotaxic head frame. The skull of the mouse was exposed, and a small opening was made using a sterile dental drill (Precision Guide) at 2 mm lateral and 3 mm posterior to lambda. At this location, 2 × 105 MB002, D425, or ONS76 cells in 2–3 µl culture media were slowly injected (over 2 min) 2 mm deep to the surface of the skull using a 26G Hamilton syringe. All procedures were performed under sterile conditions.
In vivo bioluminescence monitoring
In vivo bioluminescence imaging was performed using the Xenogen IVIS Lumina System coupled to LivingImage software for data acquisition. Mice were anesthetized using gaseous isoflurane and imaged 10 min after intraperitoneal injection of luciferin. Signal intensity was quantified within a region of interest over the head defined by the LivingImage software.
ONS76 cells were seeded on coverslips for 48–72 h and voltage-clamped in perforated whole-cell configurations at room temperature using an Axopatch 200B amplifier (Molecular Devices). Cells were placed in an electrophysiology recording chamber and bathed with (in mM) 130 N-methyl-d-glucamine chloride, 2 MgCl2, 2 CaCl2, 10 Hepes, 10 glucose, and 5 tetraethylammonium chloride, pH 7.4, with NMDG. The internal (pipette) solution consisted of (in mM): 145 KCl, 1 MgCl2, and 10 Hepes, pH 7.25, with KOH. 2.5 µg/ml Gramicidin was added to the internal solution before each recording experiment. Patch electrodes (borosilicate glass) with a resistance of ∼4 MΩ were used. Pipette and whole-cell capacitance were compensated. The voltage protocol consisted of 600-ms pulses from −80 mV to +80 mV (20-mV voltage steps). The data were filtered at 2 kHz, digitized at 10 kHz, and analyzed using pClamp10 (Molecular Devices). Current-voltage curves were generated by plotting peak current amplitude at different voltages. Quantified data were graphed using GraphPad Prism 5.
Data availability statement
The RNA-seq raw data of this study are available from the Sequence Read Archive repository (accession numbers PRJNA517338 and PRJNA592377).
Code availability statement
The code for bioinformatics analysis can be acquired from https://github.com/SiyiWanggou/CLIC1_in_Medulloblastoma.
Datasets and data processing
GSE50161 is gene expression data of nontumoral brain tissue, MB, GBM, pilocytic astrocytoma, and ependymoma. The gene expression profiles were generated by Affymetrix HG-U133plus2 chips.
GSE66354 is gene expression data of nontumoral brain tissue, MB, GBM, pilocytic astrocytoma, ependymoma, and AT/RT. The gene expression profiles were generated by Affymetrix HG-U133plus2 chips.
GSE74195 is gene expression data of PNET, MB, ependymoma, and normal brain tissue from the cerebellum. The gene expression profiles were generated by HG-U133plus2Chips.
GSE21140 is gene expression data of MB from different molecular subgroups. The gene expression profiles were generated by Affymetrix Human Exon 1.0 ST Array platform.
GSE37418 is gene expression data of MB from different molecular subgroups.
GSE85217 includes a total of 763 MB samples, representing four molecular subgroups.
GSE37382 is a dataset that includes MB samples from molecular subgroups of SHH, WNT, and group 4.
For GEO microarray datasets, background subtraction, normalization, and expression summarization of raw files were performed using the Robust Multi-array Average algorithm in R and Bioconductor packages affy and beadarray.
Statistics for CLIC1 expression comparison in high-throughput datasets
Equality of variances was tested by Leven’s test. For equal variances assumed data, one-way ANOVA was used for mean comparison and least significant difference test was used for stepwise multiple comparison of P values. For equal variances not assumed data, we used Welch’s test for mean comparison and Tamhane’s T2 test for stepwise multiple comparison of P values.
Determining CLIC1 correlated ion channel genes in MB
The GSE85217 dataset was processed for detecting CLIC1 correlated ion channel genes in MB. We used the R package psych to calculate matrix correlations and probabilities of all the gene pairs in the dataset. Gene pairs with r value >0.3 and false discovery rate value <0.05 were identified as significant. Finally, we regained a matrix containing 763 MB patients and 48 CLIC1 significantly correlated ion channel genes. The hierarchical clustering with agglomerative average linkage targeting the new matrixes was performed with the R/ConsensusClusterPlus package and the heatmap were generated via R/Heatmap function by Z score.
RNA-seq analysis of CLIC1 knockdown in DAOY and ONS76 cell lines
The adapters of fastq files were trimmed by TrimGalore. The clean paired-end reads were mapped to GRCh38 genome reference using HISAT2. StringTie was used for transcripts assembly and quantification. Differential expression was determined by DEseq2. Differentially expressed genes with an adjusted P value (false discovery rate) < 0.1 and fold change >2 (or <0.5) were accepted as significant.
For differentially expressed genes in RNA-seq, we performed enrichment analysis targeting on pathways in Gene Ontology and Kyoto Encyclopedia of Genes and Genomes. GSEA was used for enrichment analysis via the R package clusterprofiler. Pathways with a false discovery rate <0.1 were considered as enriched. The GSEA enrichment plots were generated by R package enrichplot.
Total RNA estimation
The total RNA from RNA-seq in each DAOY sample was determined by combining the read counts of all genes. The estimated total RNA was generated by dividing total RNA-seq read counts with average read counts of reference genes. Housekeeping genes, including ABCB7, ABCD3, ABCF1, ABCF2, MFSD3, MFSD5, SLC25A28, SLC25A38, SLC25A39, SLC30A5, SLC35B1, SLC39A7, and SLC46A3, were used as reference genes.
No statistical methods were used to predetermine sample sizes. The statistical analyses were done afterwards without interim data analysis. No data points were excluded. A two-tailed Student’s t test was performed for comparison between two groups of samples. Two-way ANOVA analyses were used to assess significance of multiple data points. The Kaplan–Meier estimator and GraphPad Prism software were used to generate survival curves. Differences between survival curves were calculated using a log-rank test. The data meet the assumptions of the tests. The variance has been tested in each group of the data, and the variance is similar among genotypes. Data distribution was assumed to be normal, but this was not formally tested. All data were collected and processed randomly. Each experiment has been successfully reproduced at least three times and was performed on multiple days. All data are expressed as mean ± SEM. We considered a P value <0.05 to be statistically significant.
Online supplemental material
Fig. S1 shows characterizations of human MB cell lines and xenograft tumors. Fig. S2 shows that Clic1 knockout does not affect mouse development, survival, or body weight. Fig. S3 shows the analyses of cell volume, cell cycle, and the p38 MAPK pathway of human MB cells. Fig. S4 shows the subcellular localization of CLIC1 and KCNT2 and identification of CLIC1-correlated ion channels in human MB tumors. Fig. S5 shows that Clic mutations or RNAi knockdown of Clic, Eag, or both does not elicit overt phenotype in Drosophila larval neuroblast, progenitor, or glia cells.
We are grateful to Lily Jan and Yuh Nung Jan for their support during the beginning phase of this project and Susan Younger for her mentoring on Drosophila genetics. The Dpn antibody was a gift from the Jan laboratory. We thank Paul Paroutis and Kimberly Lau at the SickKids Imaging Facility for their help with confocal imaging and image analysis. We thank Chi-chung Hui and Peter Dirks for constructive suggestions on the project and their laboratories for technical support. We thank all Huang laboratory members for critical reading of the manuscript. X. Huang acknowledges love and support from Lucas Huang and Liam Huang.
This work is supported by the SickKids Foundation, the Arthur and Sonia Labatt Brain Tumour Research Centre, Garron Family Cancer Centre, b.r.a.i.n.child, Meagan’s Walk, the Natural Sciences and Engineering Research Council of Canada (Discovery Grant RGPIN-2016-05644), the Ontario Institute for Cancer Research (Translational Research Initiative), the US Department of Defense Peer Reviewed Cancer Research Program (Career Development Award CA160264), the Canadian Institutes of Health Research (Project Grants 377122 and 406569), and the Sontag Foundation (Distinguished Scientist Award to X. Huang). X. Huang is a Catalyst Scholar of The Hospital for Sick Children and Canada Research Chair in Cancer Biophysics.
Author contributions: Conception of project: X. Huang; Experimental design, data acquisition and interpretation: M.A. Francisco, S. Wanggou, J.J. Fan, W. Dong, X. Chen, A. Momin, N. Abeysundara, H-K. Min, J. Chan, R. McAdam, M. Sia, R.J. Pusong, S. Liu, N. Patel, V. Ramaswamy, N. Kijima, L-Y. Wang, Y. Song, R. Kafri, M.D. Taylor, X. Li, and X. Huang; Manuscript writing: M.A. Francisco, J.J. Fan, and X. Huang.
Disclosures: Dr. Ramaswamy reported personal fees from Astra-Zeneca outside the submitted work. No other disclosures were reported.