Skip to main content
  • AACR Publications
    • Blood Cancer Discovery
    • Cancer Discovery
    • Cancer Epidemiology, Biomarkers & Prevention
    • Cancer Immunology Research
    • Cancer Prevention Research
    • Cancer Research
    • Clinical Cancer Research
    • Molecular Cancer Research
    • Molecular Cancer Therapeutics

AACR logo

  • Register
  • Log in
  • My Cart
Advertisement

Main menu

  • Home
  • About
    • The Journal
    • AACR Journals
    • Subscriptions
    • Permissions and Reprints
  • Articles
    • OnlineFirst
    • Current Issue
    • Past Issues
    • Rapid Impact Archive
    • Meeting Abstracts
    • Collections
      • COVID-19 & Cancer Resource Center
      • Metabolism Collection
      • Editors' Picks
      • "Best of" Collection
  • For Authors
    • Information for Authors
    • Author Services
    • Best of: Author Profiles
    • Spotlight on Genomic Analysis of Rare and Understudied Cancers
    • Submit
  • Alerts
    • Table of Contents
    • Editors' Picks
    • OnlineFirst
    • Citation
    • Author/Keyword
    • RSS Feeds
    • My Alert Summary & Preferences
  • News
    • Cancer Discovery News
  • COVID-19
  • Webinars
  • Search More

    Advanced Search

  • AACR Publications
    • Blood Cancer Discovery
    • Cancer Discovery
    • Cancer Epidemiology, Biomarkers & Prevention
    • Cancer Immunology Research
    • Cancer Prevention Research
    • Cancer Research
    • Clinical Cancer Research
    • Molecular Cancer Research
    • Molecular Cancer Therapeutics

User menu

  • Register
  • Log in
  • My Cart

Search

  • Advanced search
Molecular Cancer Research
Molecular Cancer Research
  • Home
  • About
    • The Journal
    • AACR Journals
    • Subscriptions
    • Permissions and Reprints
  • Articles
    • OnlineFirst
    • Current Issue
    • Past Issues
    • Rapid Impact Archive
    • Meeting Abstracts
    • Collections
      • COVID-19 & Cancer Resource Center
      • Metabolism Collection
      • Editors' Picks
      • "Best of" Collection
  • For Authors
    • Information for Authors
    • Author Services
    • Best of: Author Profiles
    • Spotlight on Genomic Analysis of Rare and Understudied Cancers
    • Submit
  • Alerts
    • Table of Contents
    • Editors' Picks
    • OnlineFirst
    • Citation
    • Author/Keyword
    • RSS Feeds
    • My Alert Summary & Preferences
  • News
    • Cancer Discovery News
  • COVID-19
  • Webinars
  • Search More

    Advanced Search

Chromatin, Epigenetics, and RNA Regulation

Epigenetic Regulation of ZBTB18 Promotes Glioblastoma Progression

Vita Fedele, Fangping Dai, Anie P. Masilamani, Dieter H. Heiland, Eva Kling, Ana M. Gätjens-Sanchez, Roberto Ferrarese, Leonardo Platania, Soroush Doostkam, Hyunsoo Kim, Sven Nelander, Astrid Weyerbrock, Marco Prinz, Andrea Califano, Antonio Iavarone, Markus Bredel and Maria S. Carro
Vita Fedele
1Department of Neurosurgery, Medical Center, University of Freiburg, Freiburg, Germany.
2Faculty of Medicine, University of Freiburg, Freiburg, Germany.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Fangping Dai
1Department of Neurosurgery, Medical Center, University of Freiburg, Freiburg, Germany.
2Faculty of Medicine, University of Freiburg, Freiburg, Germany.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Anie P. Masilamani
1Department of Neurosurgery, Medical Center, University of Freiburg, Freiburg, Germany.
2Faculty of Medicine, University of Freiburg, Freiburg, Germany.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Dieter H. Heiland
1Department of Neurosurgery, Medical Center, University of Freiburg, Freiburg, Germany.
2Faculty of Medicine, University of Freiburg, Freiburg, Germany.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Eva Kling
1Department of Neurosurgery, Medical Center, University of Freiburg, Freiburg, Germany.
2Faculty of Medicine, University of Freiburg, Freiburg, Germany.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Ana M. Gätjens-Sanchez
1Department of Neurosurgery, Medical Center, University of Freiburg, Freiburg, Germany.
2Faculty of Medicine, University of Freiburg, Freiburg, Germany.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Roberto Ferrarese
1Department of Neurosurgery, Medical Center, University of Freiburg, Freiburg, Germany.
2Faculty of Medicine, University of Freiburg, Freiburg, Germany.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Leonardo Platania
1Department of Neurosurgery, Medical Center, University of Freiburg, Freiburg, Germany.
2Faculty of Medicine, University of Freiburg, Freiburg, Germany.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Soroush Doostkam
3Institute of Neuropathology, Neurocenter, and Comprehensive Cancer Center, University of Freiburg, Freiburg, Germany.
4BIOSS Centre for Biological Signalling Studies, University of Freiburg, Freiburg, Germany.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Hyunsoo Kim
5The Jackson Laboratory for Genomic Medicine, Farmington, Connecticut.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Sven Nelander
6Department of Immunology, Genetics and Pathology and Science for Life Laboratories, University of Uppsala, Uppsala, Sweden.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Astrid Weyerbrock
1Department of Neurosurgery, Medical Center, University of Freiburg, Freiburg, Germany.
2Faculty of Medicine, University of Freiburg, Freiburg, Germany.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Marco Prinz
3Institute of Neuropathology, Neurocenter, and Comprehensive Cancer Center, University of Freiburg, Freiburg, Germany.
4BIOSS Centre for Biological Signalling Studies, University of Freiburg, Freiburg, Germany.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Andrea Califano
7Institute for Cancer Genetics, Columbia University, New York, New York.
8Department of Biomedical Informatics, Columbia University, New York, New York.
9Department of Systems Biology, Columbia University, New York, New York.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Antonio Iavarone
7Institute for Cancer Genetics, Columbia University, New York, New York.
10Department of Pathology, Columbia University Medical Center, New York, New York.
11Department of Neurology, Columbia University Medical Center, New York, New York.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Markus Bredel
12Department of Radiation Oncology, Comprehensive Cancer Center, University of Alabama at Birmingham School of Medicine, Birmingham, Alabama.
13Department of Neurosurgery, Stanford Cancer Institute, Stanford University School of Medicine, Stanford, California.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Maria S. Carro
1Department of Neurosurgery, Medical Center, University of Freiburg, Freiburg, Germany.
2Faculty of Medicine, University of Freiburg, Freiburg, Germany.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: maria.carro@uniklinik-freiburg.de
DOI: 10.1158/1541-7786.MCR-16-0494 Published August 2017
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

Abstract

Glioblastoma (GBM) comprises distinct subtypes characterized by their molecular profile. Mesenchymal identity in GBM has been associated with a comparatively unfavorable prognosis, primarily due to inherent resistance of these tumors to current therapies. The identification of molecular determinants of mesenchymal transformation could potentially allow for the discovery of new therapeutic targets. Zinc Finger and BTB Domain Containing 18 (ZBTB18/ZNF238/RP58) is a zinc finger transcriptional repressor with a crucial role in brain development and neuronal differentiation. Here, ZBTB18 is primarily silenced in the mesenchymal subtype of GBM through aberrant promoter methylation. Loss of ZBTB18 contributes to the aggressive phenotype of glioblastoma through regulation of poor prognosis–associated signatures. Restitution of ZBTB18 expression reverses the phenotype and impairs tumor-forming ability. These results indicate that ZBTB18 functions as a tumor suppressor in GBM through the regulation of genes associated with phenotypically aggressive properties.

Implications: This study characterizes the role of the putative tumor suppressor ZBTB18 and its regulation by promoter hypermethylation, which appears to be a common mechanism to silence ZBTB18 in the mesenchymal subtype of GBM and provides a new mechanistic opportunity to specifically target this tumor subclass. Mol Cancer Res; 15(8); 998–1011. ©2017 AACR.

Introduction

Glioblastoma (GBM) is the most malignant primary brain tumor, characterized by a highly invasive nature, poor prognosis, and resistance to aggressive therapies (1). Over the past decade, gene expression profiling has contributed to the identification of multiple GBM subclasses with distinct molecular and clinical characteristics (2, 3). In particular, the mesenchymal (MES) and the proneural (PN) groups appear as the most consistent subclasses reported in both studies (2, 3). The MES subtype is characterized by resistance to radiotherapy (4).

Using bioinformatics tools (5, 6), several transcription factors have been identified as master regulators of a “mesenchymal gene expression signature” (MGES) in GBM, including STAT3, CEBPB, and WWTR1 (a.k.a. TAZ; refs. 7, 8). More recently, a role of NF-κB in controlling the expression of the three master regulators and consequent mesenchymal differentiation was reported (4). Notably, the transcriptional repressor Zinc Finger and BTB Domain Containing 18 (ZBTB18; formerly ZNF238) was identified as a potential negative regulator of the MGES in GBM (8). ZBTB18 is a transcription factor that belongs to the Broad complex, Tramtrack, Bric à brac (BTB) or poxvirus and zing finger (POZ)-zinc finger (BTB/POZ-ZF) protein family and plays a crucial role in brain development and neuronal differentiation (9–12). Previous findings revealed that ZBTB18 is downregulated or lost in mouse gliomas and human GBM cell lines and have implicated ZBTB18 as a putative tumor suppressor in the brain (11). However, the mechanism of ZBTB18 downregulation in GBM remains to be defined.

Mounting evidence suggests that alteration of methylation pathways, which can induce silencing of tumor suppressor genes, is one of the earliest events in carcinogenesis that grant a predisposition to mutational changes (13, 14). The DNA repair gene O-6-methylguanine-DNA methyltransferase (MGMT) is one such gene for which promoter methylation has been shown to result in gene silencing in many cancers (15–17).

Here, we have characterized the role of ZBTB18 as a transcriptional repressor of gene signatures associated with aggressive properties and poor prognosis in GBM. Our results indicate that ZBTB18 serves as a tumor suppressor in the brain and is silenced by DNA methylation in mesenchymal GBM.

Materials and Methods

Tumor samples and culture of GBM-derived cells

GBM and cortical samples from epilepsy surgery were collected at the Department of Neurosurgery of the University Medical Center Freiburg (Freiburg, Germany) in accordance with an Institutional Review Board–approved protocol. Informed consent was obtained from all patients, in accordance with the declaration of Helsinki. Patient-derived GBM stem cells (BTSCs) were prepared from tumor specimens as previously described (18). For passaging, neurospheres were incubated in nonenzymatic cell dissociation solution (Sigma) and mechanically dissociated by pipetting. Proneural cells 3047, 3082, and 3111 were generated at the University of Uppsala (19). GBM cell lines (SNB19 and LN229) and HEK 293T cells were routinely grown in DMEM with 10% FBS. SNB19 and LN229 cells have been authenticated on 3/2/2017 by PCR-single-locus-technology (Eurofins Medigenomix). All cells were mycoplasma-free.

Classification of brain tumor stem cells

The classification of brain tumor stem cells (BTSC) was performed using 510 genes of the 840 classifier genes used by Verhaak and colleagues to classify 260 GBM samples (3) and 529 GBM tissue samples from The Cancer Genome Atlas (TCGA) with assigned subtypes for reference (Cancer Genome Atlas Research Network, 2008). The 510 genes were selected such that they reliably classify the extended set of 529 TCGA samples and were represented on the Illumina HumanHT-12v3 expression BeadChip arrays. The expression levels for these genes on the Illumina arrays and in the TCGA dataset were converted into z-scores, and the combined matrix was used to classify each BTSC sample based on a k-nearest neighbors (k = 10) and voting procedure, in which a subtype was assigned on the basis of the majority subtype among the 10 TCGA samples with highest correlation coefficients for these genes with respect to the BTSC sample. All data manipulations were performed in R (R Core Team, 2012) and MATLAB (The MathWorks, Inc.).

The classification of BTSCs was performed using a method of hierarchical clustering with a Euclidean distance metric to cluster the samples alongside the TCGA samples. Two mRNA datasets with subtype assignment were used as references: an mRNA dataset of 529 patient samples from TCGA assignment and a subset of the 810 signature genes published by Verhaak and colleagues (3). Gliomas (BT) and BTSC classification is reported in Supplementary Table S1.

ZBTB18 vector construction, lentiviral production, and infection

ZBTB18 coding sequence was PCR amplified from normal brain RNA using the following primers: BstXI-flag-hZBTB18v1-sense (TGGCCACAACCATGGACTAC AAGGACGACGATGACAAGTGTCCTAAAGGTTATGAAGACAG) and PmeI-hZBTB18-antisense (GCCTTGGTTTAAACTTATTTCCAAAGTTCTTGAGAG). The PCR product was first cloned into pDrive cloning vector (Qiagen), sequence validated, and subsequently transferred in the pCHMWS-EGFP lentiviral vector (kind gift from V. Baekelandt, University of Leuven, Leuven, Belgium). Lentiviral infections were performed as previously described (18).

Proteomic analysis

For protein analysis of the 30 kDa band detected by the ZBTB18 antibody, SNB19 cells transduced with either control or FLAG-ZBTB18 lentiviral vector were subjected to immunoprecipitation using the M2 FLAG antibody (Sigma). The immunoprecipitated proteins were analyzed by gel electrophoresis and stained by Coomassie blue staining. The 30 kDa band was cut from the gel and analyzed by MS (Agilent 6520 Q-TOF) at the Core Facility Proteomics of the Center for Biological System Analysis (ZBSA) at the University of Freiburg (Freiburg, Germany). The identified peptides were aligned using the mascot software (Matrix Science).

DNA/RNA extraction and quantitative real-time-PCR

DNA and RNA were extracted from human cortex, tumor tissue, or cell culture using the All Prep DNA/RNA Mini Kit or miRNeasy Mini Kit (Qiagen). First-strand cDNA synthesis was generated using the Superscript cDNA synthesis Kit (Invitrogen). Quantitative RT-PCR was performed using SYBR Green (Applied Biosystems) and analyzed relative to 18sRNA (housekeeping) using the ΔCt method. Primer sequences are listed in Supplementary Table S2.

Immunoblotting

Total protein extracts were prepared in RIPA buffer supplemented with protease inhibitor cocktail (Thermo Scientific), phosphatase inhibitor cocktail (Sigma), and PMSF (Sigma). The following antibodies were used: mouse anti-FLAG M2 (Sigma), rabbit anti-ZBTB18 (Abcam, ab118471), and mouse anti–β-actin (Abcam, ab7291).

Microarray expression

For microarray expression profiling, total RNA was prepared using the RNeasy Kit or the all Prep DNA/RNA/Protein mini Kit (Qiagen) and quantified using 2100 Bioanalyzer (Agilent). A total of 1.5 μg of total RNA was processed and analyzed at the Deutsches Krebsforschungszentrum (DKFZ) in Heidelberg (Germany). Hybridization was carried out on Illumina HumanHT-12v4 expression BeadChip. Microarray data were analyzed using the GSEA software (http://www.broadinstitute.org/gsea/index.jsp). Microarray gene accession number: GSE97350 (subseries GSE97347 and GSE97349).

Methylation analysis and pyrosequencing

Genomic DNA from tissues and cell lines was bisulfite modified using the EZ DNA methylation–gold Kit (The Epigenetics Company) according to the manufacturer's instructions. A pool of normal brain tissues or glioma samples were subjected to PCR to amplify ZBTB18 promoter CpG islands using a PyroMark PCR kit (Qiagen) and primers summarized in Supplementary Table S3. Non–CpG-containing primers, methylation-specific primers, and unmethylation-specific primers were designed to cover the entire region. PCR products were cloned into pDrive cloning vector (Qiagen) and submitted for sequencing with SP6 primer and T7 primer (GATC).

Pyrosequencing analysis was performed using a PyroMark Q96 instrument (Qiagen), following the manufacturer's protocol. Primers are summarized in Supplementary Table S3. The results were analyzed by PyroMark CpG software (Qiagen). The methylation index for each sample was calculated as the average value of CpG methylation in the CpG examined.

In vitro methylation and luciferase reporter assay

Different regions of the ZBTB18 promoter were amplified by PCR from brain (cortex)-derived RNA and cloned in the pGL3 vector (Promega). The activity of each promoter was measured by dual-luciferase reporter assay system (Promega) according to the manufacturer's instructions. Cells were seeded at 50% confluence in 6-well plates 24 hours prior to transfection. Cells at 80% to 90% confluence were transiently transfected with 2.5 μg of each pGL3-ZBTB18 promoter vector or pGL3 as a negative control, together with 0.5 μg of pRL-TK (Renilla luciferase reporter, Promega). Transfections were done using Lipofectamine 2000 as directed by the manufacturer (Life Technologies) in serum-free medium. All transfections were done in triplicate. After 48 hours, cells were lysed with passive lysis buffer (Promega). Luciferase activity of each sample was determined in technical triplicate using a Thermo Scientific Appliskan luminometer. All data were reported relative to luciferase activity (firefly/renilla). Twenty micrograms of pGL3-ZBTB18 promoter 4 vector was methylated with SssI methylase or with HpaII methylase (all from New England Biolabs, 2.5 U/μg DNA) in the presence of 160 μmol/L S-Adenosylmethionine (SAM; New England Biolabs), in a manufacturer-supplied buffer at 37°C for 2 hours. The unmethylated DNA was treated as above but without methylases or SAM. The plasmid DNA was extracted by phenol/chloroform, ethanol precipitated, and quantified using a Nanodrop Spectrophotometer (Peqlab). The completion of methylation reaction was controlled by digesting both methylated and unmethylated DNAs using the methylation-sensitive restriction enzyme Hpa II and the methylation-insensitive restriction enzyme McrBC.

Chromatin immunoprecipitation

Promoter analysis was performed with the MatInspector software (www.Genomatix.de). Primers were designed using the Primer 3 software (http://bioinfo.ut.ee/primer3-0.4.0/) and are listed in Supplementary Table S4. Chromatin immunoprecipitation (ChIP) was performed as previously described (20) with some modifications. SNB19 cells expressing ZBTB18 were first crosslinked with 1% formaldehyde (Polysciences) for 30 minutes at room temperature and quenched by addition of 125 mmol/L Glycine for 15 minutes. Lysates were sonicated for 15 minutes (30 seconds on/30 seconds off) using Branson digital sonifier and centrifuged at 14,000 rpm for 15 minutes at 4°C. Fifty micrograms of sonicated chromatin per IP was precleared and incubated with primary antibody (4 μg) anti-ZBTB18 (Abcam, ab118471). Chromatin–antibody complexes were eluted by two 15-minute incubations at 30°C with 250 μL Elution Buffer (1% SDS, 100 mmol/L NaHCO3). Chromatin was reverse-crosslinked by adding 20 μL of NaCl 5 mol/L and incubated at 65°C for 12 hours. DNA was extracted by phenol-chloroform after RNase and proteinase K digestion. Quantitative RT-PCR was performed using SYBR Green (Applied Biosystems).

Methylation array

For methylation analysis, DNA was prepared using the All Prep DNA/RNA/Protein Mini Kit (Qiagen) and quantified using NanoDrop2000c (Thermo Scientific). A total of 1.5 μg of DNA for each sample was processed at DKFZ by Illumina HumanMethylation450 (Illumina). Data analysis was performed using Integrative Genomics Viewer (IGV) software (Broad Institute; refs. 21, 22).

Migration, invasion, and proliferation assay

Migration and invasion assays were performed as described before (18). Images of migrating cells were taken every 24 hours. For BTSC233 and JX6 cells, laminin-coated (Invitrogen; 4 μg/mL) 60 mm dishes containing a culture insert (Ibidi) were used. Cell migration was calculated using the following formula: (premigration area − migration area)/premigration area × 100).

For invasion assay, 1 × 105 BTSC233, JX6, or SNB19 (both 2.5 × 104) cells were seeded in triplicates in the upper compartment. PDGF-BB (20 ng/mL; R&D Systems) was used as a chemoattractant. Pictures were acquired using an Axioimager 2 Microscope (Zeiss). The assays were validated in two independent experiments.

Cell proliferation was assessed using a commercially available kit for EdU detection (EdU Cell proliferation assay, base click). Cells were plated at a density of 2.0 × 104 per well in a 24-well microplates containing laminin-coated coverslips. After 24 hours of seeding, cells were incubated with EdU solution overnight following the manufacturer's instructions. Upon EdU detection, images were acquired using a fluorescent microscope (Axiovert; Zeiss).

5-AZA-2′-deoxycytidine treatment

A total of 100 mmol/L stock solutions of 5-AZA-2′-deoxycytidine (5-AZA-dC; Sigma-Aldrich) were prepared by dissolving the substances in DMSO (GIBCO) and stored at −80°C. Immediately before treatment, stock solutions were diluted in cold PBS and added to the cell culture medium. LN229 cells not expressing ZBTB18 were treated with 10 μmol/L 5-aza-dC for 2, 3, 4, 5, and 8 days. JX6 cells and BTSC161s cells were treated with equal amounts of 5-aza-dC for 3 and 6 days. Medium containing fresh drug was changed every 24 hours.

Intracranial injection and immunohistochemistry

Intracranial injections were performed in NOD/SCID mice (Charles River Laboratories) in accordance with the directive 86/609/EEC of the European Parliament, following approval by regional authorities. Experiments were performed as described before (8). Animals were monitored daily until the development of neurologic symptoms by a blinded operator. Histology was performed as previously described (23). For MIB1 staining, primary antibody MIB1 (1:50; Dako; M7240) was used.

Statistical analysis

Linear regression analyses and graphical model validation were executed using R software. Scatterplots and locally weighted least squares (LOWESS) smooths were used to confirm the suitability of linear regression analyses, and statistical significance of these relationships was assessed according to the P value for the estimated slope of the regression line. A multiple linear regression model was computed based on the ordinary least squares method. Expression analysis of the TCGA data was performed by R software. Publicly available Level 3 TCGA (https://tcga-data.nci.nih.gov/docs/publications/tcga/) data were used for analysis. Data were downloaded at the UCSC Cancer Genome Browser. Only patients with full datasets (expression, methylome, and clinical information) were included. Expression analysis was performed based on Agilent array data (TCGA GBM G4502A) for high-grade glioma and RNA-seq data (TCGA LGG HiSeqV2 PANCAN) for low-grade glioma. Both datasets were normalized and log2 transformed. Methylation data from Infinium HumanMethylation450 BeadChip for lower-grade and high-grade glioma were used for further analysis. For TCGA Expression/Methylation Analysis, normalized expression values were analyzed in tumor subtypes (mesenchymal, proneural, classical, and neural) by a one-way ANOVA model. Survival analysis was performed by Cox proportional hazards model and plotted by Kaplan–Meier survival statistics. Patients without survival data were censored. Robustness was ensured by 10-fold cross-validation. Methylation of the cg23829949 was extracted and analyzed by Wilcox regression model and one-way ANOVA. Significant level was defined as P < 0.001. Analysis was performed by survival package included in R-Software. Different tumor numbers used in the analyses reflects sample availability in different genomic datasets and their overlap.

Results

ZBTB18 is downregulated in high-grade gliomas

To study the role of ZBTB18 in gliomas, we looked at its expression in 1,161 low- and high-grade gliomas from TCGA. ZBTB18 expression was lower in GBM (WHO grade IV) compared with grade II and III gliomas (P < 0.001; Fig. 1A). Interestingly, ZBTB18 expression appeared strongly associated with the proneural subclass when we focused our analysis on high-grade gliomas only (n = 561; Fig. 1B). This finding was confirmed by cluster analysis of ZBTB18 correlating genes in high-grade gliomas (Fig. 1C). ZBTB18 was mostly associated with the proneural subtype and expressed at lower levels in the mesenchymal and classical subtype (Fig. 1C), consistent with our previous identification of ZBTB18 as a putative transcriptional repressor of the MGES described by Phillips and colleagues (2, 8). ZBTB18 protein expression analysis showed low expression of ZBTB18 in GBM-derived cells compared with normal brain samples (cortex; Fig. 1D and E), further reinforcing the notion that ZBTB18 is downregulated in high-grade gliomas. Interestingly, a lower band around 30 kDa was also detected.

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1.

ZBTB18 is highly expressed in low-grade gliomas and proneural GBMs. A, Analysis of ZBTB18 expression in low- and high-grade glioma samples from TCGA. B, Analysis of ZBTB18 expression in GBM subtypes from TCGA samples. C, Cluster analysis showing association between expression of ZBTB18-correlating genes with methylation and tumor subtype in GBM samples from TCGA. ZBTB18 appears to be mostly expressed in proneural GBMs. D, Western blot showing ZBTB18 expression in normal brain tissues and GBM-derived cells. The low band represents a shorter ZBTB18 isoform. α-Tubulin was used as loading control. E, Densitometric analysis of the Western blot displayed in D. The ZBTB18 signal was normalized to the α-tubulin signal.

ZBTB18 directly regulates MGES genes

Our previous study identified ZBTB18 as a putative negative regulator of the MGES of GBM (8). Moreover, cluster analysis shown in Fig. 1B indicates that ZBTB18 is downregulated in both mesenchymal and classical GBMs, suggesting that low ZBTB18 expression is associated with the more aggressive GBM subtypes. To validate ZBTB18-repressive function, we first analyzed the expression of mesenchymal genes previously predicted to be negatively connected to ZBTB18 in the ARACNE network (8), using SNB19 GBM cell lines transduced with FLAG-ZBTB18 or control vector (Fig. 2A; Supplementary Fig. S3A and S3B). Among 13 putative targets analyzed, 8 genes were clearly downregulated upon ZBTB18 overexpression (Fig. 2A). ChIP revealed ZBTB18 binding at ACTN1, PTRF, SERPINE1, and CD97 promoters (Fig. 2B), indicating that, at least for a subset of targets, ZBTB18 is directly involved in gene repression. To better understand the role of ZBTB18 in high-grade gliomas, we performed overexpression studies and subsequent GSEA (http://software.broadinstitute.org/gsea/index.jsp) on mesenchymal patient-derived BTSC-like cell line (BTSC233; Figs. 1D and E and 2C and D) and a patient-derived GBM xenoline (JX6; ref. 24) which was classified as classical according to the Verhaak study (Figs. 1D and E and 2E and F). Interestingly, we also detected a lower band similar to those observed in GBM-derived cells at around 30 kDa (Figs. 2C and E and 1D). Sequence analysis by MS confirmed that the peptides correspond to the N-terminal portion of ZBTB18 up to around amino acid 270 (Supplementary Fig. S1A and S1B). Western blot analysis of the transduced cells with a polyclonal anti-ZBTB18 antibody showed almost no expression of ZBTB18 in the control cells; the shorter band was weakly recognized because the antibody is directed to central region of ZBTB18 sequence (amino acids 228–498; Supplementary Fig. S1C and S1D). GSEA analysis for the gene signatures described by Phillips and colleagues (2) revealed a strong enrichment for mesenchymal genes in BTSC233 cells expressing a control GFP vector compared with BTSC233 cells expressing ZBTB18 (Fig. 2G; Supplementary Fig. S2A). Surprisingly, the same analysis using gene signatures from the Verhaak classification (3) did not show any specific signature enrichment (Supplementary Fig. S2B). Validation by qRT-PCR in BTSC233 expressing ectopic ZBTB18 confirmed the downregulation of several mesenchymal genes (Fig. 2H). The repressive function of ZBTB18 on a subset of genes was further validated in JX6 cells (Fig. 2I). GSEA showed a strong downregulation of the Phillips proliferative signature in JX6 transduced with a ZBTB18 expressing lentiviral vector compared with the control vector, but again no significant change was induced in the Verhaak signatures (Fig. 2J; Supplementary Fig. S3A). As the classification by Phillips is based on gene expression data of GBMs and grade III gliomas, with the goal of identifying survival-associated genes, we reasoned that the different results might indicate a role of ZBTB18 in the negative regulation of genes associated with poor survival (i.e., proliferative and mesenchymal genes). Interestingly, some of the validated ZBTB18 targets have been previously associated with unfavorable prognosis in glioma (25–27). Further examination of the top downregulated genes by ZBTB18 in JX6 cells by gene expression array highlighted many genes previously reported to play a role in epithelial-to-mesenchymal transition (EMT; Fig. 3A; refs. 28–34). ID1 and ID3 were also downregulated, as consistent with previous findings (Fig. 3A; ref. 35). The repressive role of ZBTB18 was validated by qRT-PCR in JX6 (Fig. 3B) and BTSC233 (Fig. 3C) cells. Many of the validated genes have been reported as part of a multi-cancer gene expression signature associated with prolonged time-to-recurrence in GBM (36, 37). As such, GSEA showed a strong loss of this multi-cancer signature (Anastassiou_cancer_mesenchymal_transition signature) upon ZBTB18 overexpression (Fig. 3D). These data further suggest that ZBTB18 downregulation in high-grade gliomas leads to re-expression of genes associated with malignant features and poor outcome. Interestingly, ZBTB18 re-expression in BTSC233 and JX6 cells led to upregulation of epithelial markers which are often repressed during EMT (refs. 38–43; Fig. 3E and F), further reinforcing the idea that ZBTB18 could play a role in suppressing an EMT-like phenotype in GBM. Consistent with our data and with its previously reported tumor-suppressive role (11), ZBTB18 reexpression in SNB19 affected cell proliferation, migration, and invasion (Supplementary Fig. S4). The same effect on cell proliferation, migration, and invasion was confirmed in BTSC233 cells (Supplementary Fig. S5A, S5C, and S5D). In JX6 cells, ZBTB18 overexpression also reduced cell proliferation and migration (Supplementary Fig. S5B and S5E), although no clear effect on invasion was observed, probably due to the higher invasive properties of those cells (data not shown).

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2.

Ectopic expression of ZBTB18 in GBM cells induces repression of poor prognosis signature genes. A, Representative qRT-PCR analysis showing gene expression changes of ZBTB18 mesenchymal targets previously predicted by ARACNE in SNB19 cells expressing ectopic ZBTB18 (n = 3 PCR replicates; error bars ± SD). Gene expression was normalized to 18sRNA. B, Representative ChIP experiment in SNB19 cells expressing ectopic ZBTB18. The panel shows ZBTB18 binding to the promoter of a subset of mesenchymal targets (n = 3 PCR replicates) expressed as percentage of the initial DNA amount in the immune-precipitated fraction. The OLR1 gene which does not contain putative ZBTB18 binding sites was used as a negative control. C, Western blot showing ectopic expression of ZBTB18 in BTSC233. The low band represents a shorter ZBTB18 isoform. D, Representative images of GFP-positive BTSC233 cells expressing control vector or FLAG-ZBTB18 construct after lentiviral infection. The scale bar represents 100 μm. E, Western blot showing ectopic expression of ZBTB18 in JX6. The shorter ZBTB18 isoform is indicated. F, Representative images of GFP-positive JX6 cells transduced with control vector or FLAG-ZBTB18 lentiviral vector. The scale bar represents 100 μm. G, GSEA enrichment plot for mesenchymal genes in the comparison of 233 cells expressing control vector vs. FLAG-ZBTB18. H and I, Validation by qPCR of selected mesenchymal genes in BTSC233 (H) or JX6 (I) expressing either control vector or FLAG-ZBTB18 construct (n = 3; error bars ± SD). *, P < 0.05; **, P < 0.01; and ***, P < 0.001. Gene expression was normalized to 18sRNA. J, GSEA enrichment plot for proliferative genes in the comparison of JX6 cells expressing control vector vs. FLAG-ZBTB18.

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3.

Ectopic expression of ZBTB18 in GBM cells affects a cancer-mesenchymal transition signature. A, List of the top 20 genes downregulated by ZBTB18 overexpression, in JX6 cells, analyzed by gene expression array. Genes selected for validation are in blue. B and C, Validation by qRT-PCR of selected genes listed in A in BTSC233 (B) or JX6 (C) expressing either control vector or FLAG-ZBTB18 construct (n = 3; error bars ± SD). *, P < 0.05; **, P < 0.01; and ***, P < 0.001. Gene expression was normalized to 18sRNA. D, GSEA enrichment plot for the Anastassiou_cancer_mesenchymal_transition signature described in ref. 37. E and F, qRT-PCR analysis of epithelial genes in BTSC233 (E) and JX6 (F) transduced with either control vector or FLAG-ZBTB18.

ZBTB18 inhibits brain tumor growth in vivo

We next addressed the role of ZBTB18 in tumor formation in vivo. Immunocompromised (NOD/SCID) mice were intracranially injected with JX6 or BTSC233 cells stably expressing either ZBTB18 or control-GFP vector. Histology analysis of the mouse brains revealed that mice injected with JX6 cells expressing control vector developed bulky tumors. Conversely, only one mouse injected with JX6 cells expressing ZBTB18 formed a very small tumor (Fig. 4A and data not shown). In accord, overall survival was significantly increased in the ZBTB18 group (Fig. 4B). The same experiment in BTSC233 confirmed the effect of ZBTB18 on survival even though the mice still developed tumors (Fig. 4C and D). Ki67/MIB1 staining showed a high level of proliferation in BTSC233 transduced with the control-GFP vector, whereas cells expressing ZBTB18 appeared to be less proliferative or confined in small satellite areas, suggesting that ZBTB18 might somehow impair proliferation or restrict it to specific tumor regions (Fig. 4E). These results suggest that expression of ZBTB18 prolongs animal survival by delaying or inhibiting tumor formation and that the extent of tumor inhibition might depend on the cell background.

Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4.

Expression of ZBTB18 affects tumor formation and survival in vivo. A, Hematoxylin and eosin (H&E) staining of representative tumors resulting from intracranial injection of JX6 cells infected with either control (left) or FLAG-ZBTB18 overexpressing vector (right) in immunocompromised mice. The scale bar represents 1 mm. B, Kaplan–Meier survival curves of mice intracranially injected with JX6 cells expressing full-length ZBTB18 or control vector (n = 8). C, H&E staining of representative mice tumors upon intracranial injection of BTSC233 cells infected with either control (left) or FLAG-ZBTB18 overexpressing vector (right) in immunocompromised mice. The scale bar represents 1 mm. D, Kaplan–Meier survival curves of mice intracranially injected in C (n = 9). E, Representative images of mice tumors resulting from intracranial injection of BTSC233 transduced with either control or FLAG-ZBTB18 stained with anti-MIB1 antibody. The scale bar represents 200 μm. A reduced proliferation or pattern of proliferating cells is observed upon ZBTB18 overexpression.

Promoter methylation silences ZBTB18 in GBM

To elucidate the mechanism by which ZBTB18 is downregulated in GBM, we examined the ZBTB18 promoter in silico (http://genome.ucsc.edu). The analysis revealed the presence of two CpG islands (Fig. 5A and B), suggesting that DNA hypermethylation could play a role in ZBTB18 transcriptional repression. To verify this hypothesis, we first cloned several promoter regions covering the two CpG islands from a pool of normal brains, GBM samples (BTs), and GBM cell lines after bisulfite modification (Fig. 5A and B). Sequence analysis of the cloned DNA fragments revealed no change in DNA methylation in the more upstream CpG island (CpG1, containing 27 CpGs; Fig. 5B). Conversely, higher methylation in CpG island 2 (CpG2, containing 9 CpGs) was detected in the pool of tumor samples (Fig. 5B). Pyrosequencing of promoter regions containing 6 CpGs located in CpG2 and 2 more downstream CpGs located in the 5′UTR revealed that, although not statistically significant, methylation of ZBTB18 CpG2 tended to be higher in the glioma samples compared with normal brain (Supplementary Fig. S6A). Furthermore, methylation of the 2 CpGs in the 5′UTR (5′UTR CpG -1 and 5′UTR CpG -2) was higher in gliomas compared with normal brain samples (P = 0.0184, unpaired t test; Supplementary Fig. S6B and S6C). The significance was even higher when only GBM samples were included in the analysis (P = 0.0017, unpaired t test; Fig. 5C).

Figure 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 5.

ZBTB18 promoter is methylated in a group of GBM. A, Schematic representation of strategy used for methylation analysis of ZBTB18 promoter. B, Detailed analysis of ZBTB18 promoter showing the region analyzed, the CpGs islands identified, and their relative methylation status in normal tumor samples (Control), brain tumors (BTs), and brain tumor–derived cells (called MBs). C, Comparison of the 5′UTR CpG-1 and -2 methylation status between controls (N = 33) and grade 4 brain tumors (n = 28). D, Expression of ZBTB18 in LN229, JX6, and BTSC161s cells after treatment with 5-Aza-2′-dC and relative control. Representative image of an experiment performed in triplicate (n = 3; error bars ± SD).

Treatment with the hypomethylating agent 5′-Aza-2′-dC (Decitabine) in GBM cell line LN229 and two patient-derived mesenchymal GBM cell lines (BT161s and JX6), all showing downregulation of ZBTB18 and hypermethylation of the ZBTB18 promoter, resulted in reexpression of ZBTB18 after 72 hours (Fig. 5D). Pyrosequencing analysis confirmed concomitant reduction of 5′UTR CpG-1 and 5′UTR CpG-2 in LN229 cells (Supplementary Fig. S6C). Together, these data support a role of promoter hypermethylation in the silencing of ZBTB18 in GBM.

To further prove that the ZBTB18 region, including 5′UTR CpGs-1 and -2, is important for ZBTB18 promoter activity, we cloned several ZBTB18 promoter regions and analyzed their activity by luciferase reporter assay. As shown in Fig. 6A and B, the ZBTB18 promoter region 3, which does not contain the core promoter and the 5′UTR CpGs-1 and -2, had no promoter activity compared with other cloned promoter regions in which the core promoter and 5′UTR CpGs-1 and -2 were included (Fig. 6A and B). Next, we investigated the effect of DNA methylation on the luciferase reporter activity controlled by the ZBTB18 promoter region with the highest activity. As displayed in Fig. 6C and D, the promoter activity was completely inhibited by the SssI methylase, an enzyme that methylates all CGs, and to a less extent by HpaII methylase, which methylates only CGs in the CCGG context (Fig. 6C and D). These results further indicate that the ZBTB18 promoter region, which includes the core promoter and 5′UTR CpGs-1 and -2, is responsible for promoter activity and sensitive to DNA methylation. This is consistent with the expected role of methylated CpGs close to the TSS in gene expression regulation (44).

Figure 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 6.

The core promoter region of ZBTB18 is essential for promoter activity and is sensitive to DNA methylation. A, Schematic representation of the ZBTB18 promoter regions cloned in the pGL3 luciferase reporter vector. B, Analysis of ZBTB18 promoter constructs activity by dual-luciferase assay. C, In vitro methylation assay of ZBTB18 promoter 4. The plasmid DNA was methylated with SssI or HpaII methylase. D, Control restriction digestion of the methylase reaction using the methylation-sensitive (HpaII) and the methylation-insensitive (McrBC) restriction enzymes.

ZBTB18 promoter methylation is a hallmark of mesenchymal GBMs

Correlation analysis between ZBTB18 expression and DNA methylation in 251 GBM samples from TCGA (HumanMethylation27 platform, http://cancergenome.nih.gov) revealed a significant inverse correlation between ZBTB18 expression and promoter methylation (P = 1.96 × 10−05, linear regression; Fig. 7A). Interestingly, the examined probe (cg23829949) mapped into the same region analyzed by pyrosequencing (CpG island 2). Consistently, in our set of patient-derived glioma samples, ZBTB18 promoter methylation of 6 CpGs located at the 3′ end of CpG island 2 (Fig. 7B) and of 5′UTR CpGs-1 and -2 (Supplementary Fig. S6B) measured by pyrosequencing was higher in a subset of samples showing low ZBTB18 expression (high-grade gliomas, red line), whereas low promoter methylation was detected in samples with higher ZBTB18 expression (low-grade gliomas, blue line; Fig. 7B). Thus, methylation of CpG island 2 and 5′UTR CpGs-1 and -2 correlates with ZBTB18 expression, at least in a subset of GBM. However, a fraction of glioma samples that did not show promoter hypermethylation still had low ZBTB18 expression, suggesting that DNA hypermethylation is not the only mechanism regulating ZBTB18 expression, or alternatively, that additional methylated regions might be involved.

Figure 7.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 7.

Methylation at the ZBTB18 promoter causes ZBTB18 repression in GBM. A, Linear regression of ZBTB18 expression on methylation status in n = 251 samples (P = 0.0000196). B, Methylation level of ZBTB18 in a series of glioma samples assessed by pyrosequencing (n = 3; error bars ± SD). ZBTB18 expression level measured by qRT-PCR is indicated by blue circles (high expression and low expression are indicated by dark blue and blue circles, respectively). Red line highlights the high methylation/low ZBTB18 expression samples, and blue line highlights the low methylation/high ZBTB18 expression samples. BTSC 145 (GBM-derived cell line) was analyzed instead of the corresponding BT 145 GBM sample which was no longer available. C, Association between ZBTB18 methylation and IDH1 mutation status. D, Cluster analysis of methylation profiles of the G-CIMP and non–G-CIMP GBMs using the 450/27 K BeadChip from TCGA. The association between expression of ZBTB18 and tumor subtypes (methylation or expression) is shown. ZBTB18 appears to be mostly expressed in G-CIMP GBMs. E, Association analysis between ZBTB18 promoter methylation and GBM subtypes. F and G, Representative experiment showing expression of mesenchymal targets in samples of brain tumors highlighted in A assessed by qRT-PCR (n = 3 PCR replicates; error bars ± SD). Red line, high methylation; blue line, low methylation. H, Kaplan–Meier estimates of time-to-tumor progression in 109 GBM patients from TCGA, with patients stratified according to low versus high (relative to the median) ZBTB18 methylation (log-rank P = 0.029).

Promoter methylation analysis of the subset of gliomas with low ZBTB18 expression/high promoter methylation (indicated in red) versus high ZBTB18 expression/low promoter methylation (indicated in blue) previously analyzed using the Infinium HumanMethylation450 BeadChip (ref. 45; Supplementary Fig. S7D) confirmed differential methylation of CpG island 2 between the two tumor groups covering CpGs located in the same region analyzed by pyrosequencing (probes cg19698993 and cg12869659). Intriguingly, tumor samples with high ZBTB18 expression and low ZBTB18 promoter methylation analyzed by DNA methylation array showed high levels of global methylation (45), suggesting that silencing of ZBTB18 by promoter methylation could be a hallmark of non G-CIMP gliomas. Consistent with this hypothesis, ZBTB18 methylation correlated with G-CIMP and IDH1 mutation status so that ZBTB18 was more methylated in IDH1 wild-type compared with IDH1-mutant gliomas (Fig. 7C P = 1.914e−10, Welch two-sample t test). Consistent with this, ZBTB18 expression was higher in G-CIMP GBMs (Fig. 7D). We then investigated the relationship between ZBTB18 methylation (HumanMethylation27 platform) and GBM subclasses in 283 TCGA GBM samples. This analysis revealed a highly significant association between ZBTB18 methylation in a region confirmed by pyrosequencing (probe cg23829949) and mesenchymal GBM subtype (Fig. 7E). Accordingly, regression analysis of ZBTB18 expression on ZBTB18 methylation revealed a strong correlation only in mesenchymal tumors but not in nonmesenchymal tumors (all GBMs: P = 0.000026; mesenchymal tumors: P = 0.0136; nonmesenchymal tumors: P = 0.13, linear regression), further highlighting that methylation-induced silencing of ZBTB18 might be particularly important in the mesenchymal subclass (Supplementary Fig. S7A–S7C). This is consistent with previous analysis indicating that the mesenchymal tumors are usually non–G-CIMP (46, 47). Expression analysis of selected mesenchymal genes (CD97, ACTN1, EMP3, and CHI3L1) revealed an inverse association with ZBTB18 promoter methylation (Fig. 7E and F), further indicating that silencing of ZBTB18 through promoter hypermethylation could play a role in mesenchymal differentiation in GBM.

Survival modeling did not show a statistically significant association with patient survival (data not shown). Instead, we observed a significant link between time-to-tumor progression and ZBTB18 methylation in a two-class model of 109 GBM patients stratified according to lower-than-median versus higher-than-median ZBTB18 methylation (log-rank P = 0.029), such that patients with high methylation demonstrated a comparatively unfavorable outcome (Fig. 7G). This association between ZBTB18 methylation and time-to-tumor progression was also evident in a continuous univariate Cox model [univariate Cox model P = 0.037, HR for tumor progression with methylated ZBTB18: 8.30, 95% confidence interval (CI), 1.13–60.80] and prevailed in a multivariate Cox model that included GBM subclass (classical, mesenchymal, neural, proneural) as a covariate (multivariate Cox model P = 0.048; HR, 7.49; 95% CI, 1.01–65.32), suggesting that ZBTB18 methylation portends a more aggressive tumor phenotype.

Discussion

Here, we describe a new role for the transcriptional repressor ZBTB18 as a negative regulator of signatures associated with poor survival in GBM, and we propose DNA methylation as a mechanism to silence ZBTB18 in the mesenchymal tumor subtype. Our finding is in accordance with a previous study reporting that ZBTB18 is lost in established human GBM cells and identifying ZBTB18 as a brain tumor-suppressor gene (11). We demonstrate that ZBTB18 is mostly expressed in low-grade gliomas and proneural GBMs but less expressed in mesenchymal GBMs. This is consistent with our previous identification of ZBTB18 as a transcription factor negatively associated with mesenchymal GBMs (8). ZBTB18 reexpression in primary BTSCs dampens the adherence poor-prognosis proliferative and mesenchymal signatures, which were identified as mutually exclusive to proneural in the subtype classification by Phillips and colleagues. Although the association between patient survival and mesenchymal subtype has not been confirmed in the previous TCGA study involving GBMs only (18), mouse models for glioma show that sequential mutations causing a shift from proneural to mesenchymal GBM are also associated with reduced survival (19, 20). Moreover, because tumors often show mixed subtype profiles (20, 21), it is possible that this might mask an association with survival. This was clearly shown by Patel and colleagues who demonstrated that pure IDH1 wild-type proneural tumors were associated with a better survival compared with highly heterogeneous proneural tumors containing cells with other subtypes (21). Similarly, tumor purity was recently shown to be an important parameter to determine a positive association of MGMT methylation with patient survival (17). We show that ZBTB18 attenuates the expression of genes associated with EMT and with time-to-tumor progression in GBM (36, 37). This is in line with our analysis that presented a strong link between ZBTB18 methylation and time-to-tumor progression of GBM but only a trend of association with patient overall survival (data not shown). This finding is in accordance with the prevailing argument that time-to-tumor progression might relate more closely to tumor repopulation, aggressiveness, and therapy resistance, which are biological properties also associated with mesenchymal differentiation of GBM. Given this link to tumor progression and the fact that the mesenchymal phenotype is more prevalent in recurrent GBM (2), ZBTB18 hypermethylation might play a role in both the mesenchymal differentiation characteristics of the de novo tumor and the progression toward a more mesenchymal phenotype. However, as reported in our previous study (8) and also confirmed by others (7), several transcription factors usually cooperate to regulate specific GBM subclasses. So, it would be interesting to study how deregulation of ZBTB18 fits in the previously described regulation of mesenchymal genes by other master regulators (i.e., STAT3, CEBPB, and TAZ). The recent report that, at the single-cell level, a GBM often consists of mosaic of subtype makes the picture even more complicated (48). Still, identifying regulators of different GBM subtype could be important since regulators of different subclasses co-existing in the tumor could be targeted.

The mechanism leading to ZBTB18 downregulation in GBM remained to be defined. Our data reveal a link between ZBTB18 promoter methylation and loss of ZBTB18 expression. However, we observed that, in some cases, low expression of ZBTB18 also occurs in the absence of promoter hypermethylation, suggesting that additional molecular mechanisms downregulating ZBTB18 are potentially operative in GBM. We further demonstrate a strong association between ZBTB18 promoter hypermethylation and the mesenchymal subtype of GBM, implying that silencing of ZBTB18 by promoter methylation is a particular hallmark of this unfavorable GBM subtype. We also report a strong correlation between ZBTB18 methylation and IDH1 wild-type, which is consistent with previous studies indicating that the majority of mesenchymal GBMs are non–G-CIMP (46, 47). Consistently, we show that key mesenchymal genes (2, 8), which are also differentially methylated in G-CIMP versus non–G-CIMP gliomas (46), are highly expressed in GBM exhibiting promoter hypermethylation/low expression of ZBTB18.

Collectively, our data identify ZBTB18 as a candidate tumor suppressor and transcriptional regulator of poor prognosis-associated signatures in GBM. We have identified promoter hypermethylation as a common mechanism to silence ZBTB18 in the mesenchymal subtype of GBM, which provides a new mechanistic opportunity to specifically target this tumor subclass.

Disclosure of Potential Conflicts of Interest

A. Califano is Cofounder and Chief Scientific Advisor at, and has an ownership interest (including patents) in, DarwinHealth, Inc. No potential conflicts of interest were disclosed by the other authors.

Authors' Contributions

Conception and design: V. Fedele, A.P. Masilamani, S. Nelander, A. Iavarone, M. Bredel, M.S. Carro

Development of methodology: V. Fedele, F. Dai, A.P. Masilamani, R. Ferrarese, A. Iavarone, M. Bredel, M.S. Carro

Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): A.P. Masilamani, D.H. Heiland, E. Kling, S. Doostkam, A. Weyerbrock, M. Prinz, M. Bredel, M.S. Carro

Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): F. Dai, D.H. Heiland, A.M. Gätjens-Sanchez, R. Ferrarese, H. Kim, S. Nelander, A. Califano, M. Bredel, M.S. Carro

Writing, review, and/or revision of the manuscript: V. Fedele, F. Dai, A.P. Masilamani, R. Ferrarese, S. Nelander, A. Weyerbrock, A. Califano, A. Iavarone, M. Bredel, M.S. Carro

Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): M. Prinz, M. Bredel

Study supervision: M.S. Carro

Other [Experimental assay (migration and proliferation assay); analysis of experimental assay; qPCR; ChIP; Cell culture)]: A.M. Gätjens-Sanchez

Other (technical support): L. Platania

Grant Support

This study was supported by Marie Curie International Reintegration Grant (MC IRG 268303), Deutsche Forschungsgemeinschaft Grant (DFG, CA 1246/2-1; both to M.S. Carro), and German Cancer Aid Grant Award (107714, M. Bredel). M. Prinz was supported by the DFG (SFB 992).

The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

Acknowledgments

The authors thank M. Oberle for mice brain histology, M. Lübbert for providing access to pyrosequencing facility, T. Feuerstein for normal brain tissues, C. Stein for technical assistance, R. Claus for help with methylation data analysis, and D.Ó. hAilín for editing the article (all University of Freiburg), V.D. Marinescu for BTSCs classification (University of Uppsala), S. Nozell for input on the article (UAB, Birmingham); Y. Gillespie (UAB, Birmingham) for providing JX6 cells; and V. Baekelandt (Katholieke Universiteit Leuven) for pLV-eGFP lentiviral vector.

Footnotes

  • Note: Supplementary data for this article are available at Molecular Cancer Research Online (http://mcr.aacrjournals.org/).

  • Received December 28, 2016.
  • Revision received April 7, 2017.
  • Accepted May 12, 2017.
  • ©2017 American Association for Cancer Research.

References

  1. 1.↵
    1. Dunn GP,
    2. Rinne ML,
    3. Wykosky J,
    4. Genovese G,
    5. Quayle SN,
    6. Dunn IF,
    7. et al.
    Emerging insights into the molecular and cellular basis of glioblastoma. Genes Dev 2012;26:756–84.
    OpenUrlAbstract/FREE Full Text
  2. 2.↵
    1. Phillips HS,
    2. Kharbanda S,
    3. Chen R,
    4. Forrest WF,
    5. Soriano RH,
    6. Wu TD,
    7. et al.
    Molecular subclasses of high-grade glioma predict prognosis, delineate a pattern of disease progression, and resemble stages in neurogenesis. Cancer Cell 2006;9:157–73.
    OpenUrlCrossRefPubMed
  3. 3.↵
    1. Verhaak RG,
    2. Hoadley KA,
    3. Purdom E,
    4. Wang V,
    5. Qi Y,
    6. Wilkerson MD,
    7. et al.
    Integrated genomic analysis identifies clinically relevant subtypes of glioblastoma characterized by abnormalities in PDGFRA, IDH1, EGFR, and NF1. Cancer Cell 2010;17:98–110.
    OpenUrlCrossRefPubMed
  4. 4.↵
    1. Bhat KP,
    2. Balasubramaniyan V,
    3. Vaillant B,
    4. Ezhilarasan R,
    5. Hummelink K,
    6. Hollingsworth F,
    7. et al.
    Mesenchymal differentiation mediated by NF-kappaB promotes radiation resistance in glioblastoma. Cancer Cell 2013;24:331–46.
    OpenUrlCrossRefPubMed
  5. 5.↵
    1. Basso K,
    2. Margolin AA,
    3. Stolovitzky G,
    4. Klein U,
    5. Dalla-Favera R,
    6. Califano A,
    7. et al.
    Reverse engineering of regulatory networks in human B cells. Nat Genet 2005;37:382–90.
    OpenUrlCrossRefPubMed
  6. 6.↵
    1. Margolin AA,
    2. Nemenman I,
    3. Basso K,
    4. Wiggins C,
    5. Stolovitzky G,
    6. Dalla Favera R,
    7. et al.
    ARACNE: an algorithm for the reconstruction of gene regulatory networks in a mammalian cellular context. BMC Bioinformatics 2006;7 Suppl 1:S7.
    OpenUrlCrossRefPubMed
  7. 7.↵
    1. Bhat KP,
    2. Salazar KL,
    3. Balasubramaniyan V,
    4. Wani K,
    5. Heathcock L,
    6. Hollingsworth F,
    7. et al.
    The transcriptional coactivator TAZ regulates mesenchymal differentiation in malignant glioma. Genes Dev 2011;25:2594–609.
    OpenUrlAbstract/FREE Full Text
  8. 8.↵
    1. Carro MS,
    2. Lim WK,
    3. Alvarez MJ,
    4. Bollo RJ,
    5. Zhao X,
    6. Snyder EY,
    7. et al.
    The transcriptional network for mesenchymal transformation of brain tumours. Nature 2010;463:318–25.
    OpenUrlCrossRefPubMed
  9. 9.↵
    1. Ohtaka-Maruyama C,
    2. Hirai S,
    3. Miwa A,
    4. Takahashi A,
    5. Okado H
    . The 5′-flanking region of the RP58 coding sequence shows prominent promoter activity in multipolar cells in the subventricular zone during corticogenesis. Neuroscience 2012;201:67–84.
    OpenUrlCrossRefPubMed
  10. 10.↵
    1. Okado H,
    2. Ohtaka-Maruyama C,
    3. Sugitani Y,
    4. Fukuda Y,
    5. Ishida R,
    6. Hirai S,
    7. et al.
    The transcriptional repressor RP58 is crucial for cell-division patterning and neuronal survival in the developing cortex. Dev Biol 2009;331:140–51.
    OpenUrlCrossRefPubMed
  11. 11.↵
    1. Tatard VM,
    2. Xiang C,
    3. Biegel JA,
    4. Dahmane N
    . ZNF238 is expressed in postmitotic brain cells and inhibits brain tumor growth. Cancer Res 2010;70:1236–46.
    OpenUrlAbstract/FREE Full Text
  12. 12.↵
    1. Xiang C,
    2. Baubet V,
    3. Pal S,
    4. Holderbaum L,
    5. Tatard V,
    6. Jiang P,
    7. et al.
    RP58/ZNF238 directly modulates proneurogenic gene levels and is required for neuronal differentiation and brain expansion. Cell Death Differ 2012;19:692–702.
    OpenUrlCrossRefPubMed
  13. 13.↵
    1. Cadieux B,
    2. Ching TT,
    3. VandenBerg SR,
    4. Costello JF
    . Genome-wide hypomethylation in human glioblastomas associated with specific copy number alteration, methylenetetrahydrofolate reductase allele status, and increased proliferation. Cancer Res 2006;66:8469–76.
    OpenUrlAbstract/FREE Full Text
  14. 14.↵
    1. Esteller M,
    2. Catasus L,
    3. Matias-Guiu X,
    4. Mutter GL,
    5. Prat J,
    6. Baylin SB,
    7. et al.
    hMLH1 promoter hypermethylation is an early event in human endometrial tumorigenesis. Am J Pathol 1999;155:1767–72.
    OpenUrlCrossRefPubMed
  15. 15.↵
    1. Esteller M,
    2. Hamilton SR,
    3. Burger PC,
    4. Baylin SB,
    5. Herman JG
    . Inactivation of the DNA repair gene O6-methylguanine-DNA methyltransferase by promoter hypermethylation is a common event in primary human neoplasia. Cancer Res 1999;59:793–7.
    OpenUrlAbstract/FREE Full Text
  16. 16.↵
    1. You JS,
    2. Jones PA
    . Cancer genetics and epigenetics: two sides of the same coin? Cancer Cell 2012;22:9–20.
    OpenUrlCrossRefPubMed
  17. 17.↵
    1. Schulze Heuling E,
    2. et al.
    Prognostic relevance of tumor purity and TERT promoter mutations on MGMT promoter methylation in glioblastoma. Mol Cancer Res 2017;15:532–40.
    OpenUrlAbstract/FREE Full Text
  18. 18.↵
    1. Ferrarese R,
    2. Harsh GR IV.,
    3. Yadav AK,
    4. Bug E,
    5. Maticzka D,
    6. Reichardt W,
    7. et al.
    Lineage-specific splicing of a brain-enriched alternative exon promotes glioblastoma progression. J Clin Invest 2014;124:2861–76.
    OpenUrlCrossRefPubMed
  19. 19.↵
    1. Xie Y,
    2. Bergström T,
    3. Jiang Y,
    4. Johansson P,
    5. Marinescu VD,
    6. Lindberg N,
    7. et al.
    The human glioblastoma cell culture resource: validated cell models representing all molecular subtypes. EBioMedicine 2015;2:1351–63.
    OpenUrl
  20. 20.↵
    1. Izzo A,
    2. Kamieniarz-Gdula K,
    3. Ramírez F,
    4. Noureen N,
    5. Kind J,
    6. Manke T,
    7. et al.
    The genomic landscape of the somatic linker histone subtypes H1.1 to H1.5 in human cells. Cell Rep 2013;3:2142–54.
    OpenUrlCrossRefPubMed
  21. 21.↵
    1. Thorvaldsdottir H,
    2. Robinson JT,
    3. Mesirov JP
    . Integrative Genomics Viewer (IGV): high-performance genomics data visualization and exploration. Brief Bioinform 2013;14:178–92.
    OpenUrlCrossRefPubMed
  22. 22.↵
    1. Robinson JT,
    2. Thorvaldsdóttir H,
    3. Winckler W,
    4. Guttman M,
    5. Lander ES,
    6. Getz G,
    7. et al.
    Integrative genomics viewer. Nat Biotechnol 2011;29:24–6.
    OpenUrlCrossRefPubMed
  23. 23.↵
    1. Goldmann T,
    2. Zeller N,
    3. Raasch J,
    4. Kierdorf K,
    5. Frenzel K,
    6. Ketscher L,
    7. et al.
    USP18 lack in microglia causes destructive interferonopathy of the mouse brain. EMBO J 2015;34:1612–29.
    OpenUrlAbstract/FREE Full Text
  24. 24.↵
    1. Han X,
    2. Zhang W,
    3. Yang X,
    4. Wheeler CG,
    5. Langford CP,
    6. Wu L,
    7. et al.
    The role of Src family kinases in growth and migration of glioma stem cells. Int J Oncol 2014;45:302–10.
    OpenUrlPubMed
  25. 25.↵
    1. Roy A,
    2. Coum A,
    3. Marinescu VD,
    4. Põlajeva J,
    5. Smits A,
    6. Nelander S,
    7. et al.
    Glioma-derived plasminogen activator inhibitor-1 (PAI-1) regulates the recruitment of LRP1 positive mast cells. Oncotarget 2015;6:23647–61.
    OpenUrl
  26. 26.↵
    1. Ernst A,
    2. Hofmann S,
    3. Ahmadi R,
    4. Becker N,
    5. Korshunov A,
    6. Engel F,
    7. et al.
    Genomic and expression profiling of glioblastoma stem cell-like spheroid cultures identifies novel tumor-relevant genes associated with survival. Clin Cancer Res 2009;15:6541–50.
    OpenUrlAbstract/FREE Full Text
  27. 27.↵
    1. Safaee M,
    2. Clark AJ,
    3. Oh MC,
    4. Ivan ME,
    5. Bloch O,
    6. Kaur G,
    7. et al.
    Overexpression of CD97 confers an invasive phenotype in glioblastoma cells and is associated with decreased survival of glioblastoma patients. PLoS One 2013;8:e62765.
    OpenUrlCrossRefPubMed
  28. 28.↵
    1. Bacigalupo ML,
    2. Manzi M,
    3. Espelt MV,
    4. Gentilini LD,
    5. Compagno D,
    6. Laderach DJ,
    7. et al.
    Galectin-1 triggers epithelial-mesenchymal transition in human hepatocellular carcinoma cells. J Cell Physiol 2015;230:1298–309.
    OpenUrlCrossRefPubMed
  29. 29.↵
    1. Tang J,
    2. Xiao L,
    3. Cui R,
    4. Li D,
    5. Zheng X,
    6. Zhu L,
    7. et al.
    CX3CL1 increases invasiveness and metastasis by promoting epithelial-to-mesenchymal transition through the TACE/TGF-alpha/EGFR pathway in hypoxic androgen-independent prostate cancer cells. Oncol Rep 2016;35:1153–62.
    OpenUrl
  30. 30.↵
    1. Sakamaki M,
    2. G Bemben M,
    3. Abe T
    . Legs and trunk muscle hypertrophy following walk training with restricted leg muscle blood flow. J Sports Sci Med 2011;10:338–40.
    OpenUrlPubMed
  31. 31.↵
    1. Anand V,
    2. Jaswal S,
    3. Singh S,
    4. Kumar S,
    5. Jena MK,
    6. Verma AK,
    7. et al.
    Functional characterization of Mammary Gland Protein-40, a chitinase-like glycoprotein expressed during mammary gland apoptosis. Apoptosis 2016;21:209–24.
    OpenUrl
  32. 32.↵
    1. Wang HX,
    2. Sharma C,
    3. Knoblich K,
    4. Granter SR,
    5. Hemler ME
    . EWI-2 negatively regulates TGF-beta signaling leading to altered melanoma growth and metastasis. Cell Res 2015;25:370–85.
    OpenUrlCrossRef
  33. 33.↵
    1. Chen X,
    2. Liu X,
    3. Lang H,
    4. Zhang S,
    5. Luo Y,
    6. Zhang J
    . S100 calcium-binding protein A6 promotes epithelial-mesenchymal transition through beta-catenin in pancreatic cancer cell line. PLoS One 2015;10:e0121319.
    OpenUrl
  34. 34.↵
    1. Taube JH,
    2. Herschkowitz JI,
    3. Komurov K,
    4. Zhou AY,
    5. Gupta S,
    6. Yang J,
    7. et al.
    Core epithelial-to-mesenchymal transition interactome gene-expression signature is associated with claudin-low and metaplastic breast cancer subtypes. Proc Natl Acad Sci U S A 2010;107:15449–54.
    OpenUrlAbstract/FREE Full Text
  35. 35.↵
    1. Diotel N,
    2. Beil T,
    3. Strähle U,
    4. Rastegar S
    . Differential expression of id genes and their potential regulator znf238 in zebrafish adult neural progenitor cells and neurons suggests distinct functions in adult neurogenesis. Gene Expr Patterns 2015;19:1–13.
    OpenUrl
  36. 36.↵
    1. Anastassiou D,
    2. Rumjantseva V,
    3. Cheng W,
    4. Huang J,
    5. Canoll PD,
    6. Yamashiro DJ,
    7. et al.
    Human cancer cells express Slug-based epithelial-mesenchymal transition gene expression signature obtained in vivo. BMC Cancer 2011;11:529.
    OpenUrlCrossRefPubMed
  37. 37.↵
    1. Cheng WY,
    2. Kandel JJ,
    3. Yamashiro DJ,
    4. Canoll P,
    5. Anastassiou D
    . A multi-cancer mesenchymal transition gene expression signature is associated with prolonged time to recurrence in glioblastoma. PLoS One 2012;7:e34705.
    OpenUrlCrossRefPubMed
  38. 38.↵
    1. Lubeseder-Martellato C,
    2. Hidalgo-Sastre A,
    3. Hartmann C,
    4. Alexandrow K,
    5. Kamyabi-Moghaddam Z,
    6. Sipos B,
    7. et al.
    Membranous CD24 drives the epithelial phenotype of pancreatic cancer. Oncotarget 2016;7:49156–68.
    OpenUrl
  39. 39.↵
    1. Vannier C,
    2. Mock K,
    3. Brabletz T,
    4. Driever W
    . Zeb1 regulates E-cadherin and Epcam (epithelial cell adhesion molecule) expression to control cell behavior in early zebrafish development. J Biol Chem 2013;288:18643–59.
    OpenUrlAbstract/FREE Full Text
  40. 40.↵
    1. Raschperger E,
    2. Neve EP,
    3. Wernerson A,
    4. Hultenby K,
    5. Pettersson RF,
    6. Majumdar A
    . The coxsackie and adenovirus receptor (CAR) is required for renal epithelial differentiation within the zebrafish pronephros. Dev Biol 2008;313:455–64.
    OpenUrlCrossRefPubMed
  41. 41.↵
    1. Li L,
    2. Ying J,
    3. Li H,
    4. Zhang Y,
    5. Shu X,
    6. Fan Y,
    7. et al.
    The human cadherin 11 is a pro-apoptotic tumor suppressor modulating cell stemness through Wnt/beta-catenin signaling and silenced in common carcinomas. Oncogene 2012;31:3901–12.
    OpenUrlCrossRefPubMed
  42. 42.↵
    1. Guaita S,
    2. Puig I,
    3. Franci C,
    4. Garrido M,
    5. Dominguez D,
    6. Batlle E,
    7. et al.
    Snail induction of epithelial to mesenchymal transition in tumor cells is accompanied by MUC1 repression and ZEB1 expression. J Biol Chem 2002;277:39209–16.
    OpenUrlAbstract/FREE Full Text
  43. 43.↵
    1. Gracz AD,
    2. Fuller MK,
    3. Wang F,
    4. Li L,
    5. Stelzner M,
    6. Dunn JC,
    7. et al.
    Brief report: CD24 and CD44 mark human intestinal epithelial cell populations with characteristics of active and facultative stem cells. Stem Cells 2013;31:2024–30.
    OpenUrlCrossRefPubMed
  44. 44.↵
    1. Lin JC,
    2. Jeong S,
    3. Liang G,
    4. Takai D,
    5. Fatemi M,
    6. Tsai YC,
    7. et al.
    Role of nucleosomal occupancy in the epigenetic silencing of the MLH1 CpG island. Cancer Cell 2007;12:432–44.
    OpenUrlCrossRefPubMed
  45. 45.↵
    1. Kling T,
    2. Ferrarese R,
    3. Ó hAilín D,
    4. Johansson P,
    5. Heiland DH,
    6. Dai F,
    7. et al.
    Integrative modeling reveals annexin A2-mediated epigenetic control of mesenchymal glioblastoma. EBioMedicine 2016;12:72–85.
    OpenUrl
  46. 46.↵
    1. Noushmehr H,
    2. Weisenberger DJ,
    3. Diefes K,
    4. Phillips HS,
    5. Pujara K,
    6. Berman BP,
    7. et al.
    Identification of a CpG island methylator phenotype that defines a distinct subgroup of glioma. Cancer Cell 2010;17:510–22.
    OpenUrlCrossRefPubMed
  47. 47.↵
    1. Ozawa T,
    2. Riester M,
    3. Cheng YK,
    4. Huse JT,
    5. Squatrito M,
    6. Helmy K,
    7. et al.
    Most human non-GCIMP glioblastoma subtypes evolve from a common proneural-like precursor glioma. Cancer Cell 2014;26:288–300.
    OpenUrlCrossRefPubMed
  48. 48.↵
    1. Patel AP,
    2. Tirosh I,
    3. Trombetta JJ,
    4. Shalek AK,
    5. Gillespie SM,
    6. Wakimoto H,
    7. et al.
    Single-cell RNA-seq highlights intratumoral heterogeneity in primary glioblastoma. Science 2014;344:1396–401.
    OpenUrlAbstract/FREE Full Text
PreviousNext
Back to top
Molecular Cancer Research: 15 (8)
August 2017
Volume 15, Issue 8
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Editorial Board (PDF)

Sign up for alerts

View this article with LENS

Open full page PDF
Article Alerts
Sign In to Email Alerts with your Email Address
Email Article

Thank you for sharing this Molecular Cancer Research article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Epigenetic Regulation of ZBTB18 Promotes Glioblastoma Progression
(Your Name) has forwarded a page to you from Molecular Cancer Research
(Your Name) thought you would be interested in this article in Molecular Cancer Research.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Epigenetic Regulation of ZBTB18 Promotes Glioblastoma Progression
Vita Fedele, Fangping Dai, Anie P. Masilamani, Dieter H. Heiland, Eva Kling, Ana M. Gätjens-Sanchez, Roberto Ferrarese, Leonardo Platania, Soroush Doostkam, Hyunsoo Kim, Sven Nelander, Astrid Weyerbrock, Marco Prinz, Andrea Califano, Antonio Iavarone, Markus Bredel and Maria S. Carro
Mol Cancer Res August 1 2017 (15) (8) 998-1011; DOI: 10.1158/1541-7786.MCR-16-0494

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Share
Epigenetic Regulation of ZBTB18 Promotes Glioblastoma Progression
Vita Fedele, Fangping Dai, Anie P. Masilamani, Dieter H. Heiland, Eva Kling, Ana M. Gätjens-Sanchez, Roberto Ferrarese, Leonardo Platania, Soroush Doostkam, Hyunsoo Kim, Sven Nelander, Astrid Weyerbrock, Marco Prinz, Andrea Califano, Antonio Iavarone, Markus Bredel and Maria S. Carro
Mol Cancer Res August 1 2017 (15) (8) 998-1011; DOI: 10.1158/1541-7786.MCR-16-0494
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Introduction
    • Materials and Methods
    • Results
    • Discussion
    • Disclosure of Potential Conflicts of Interest
    • Authors' Contributions
    • Grant Support
    • Acknowledgments
    • Footnotes
    • References
  • Figures & Data
  • Info & Metrics
  • PDF
Advertisement

Related Articles

Cited By...

More in this TOC Section

  • miR-200c Targets TDO2 in TNBC
  • A Nomogram to Predict Colorectal Cancer Prognosis
  • SWI/SNF Regulates Oncogenic Signaling in AML
Show more Chromatin, Epigenetics, and RNA Regulation
  • Home
  • Alerts
  • Feedback
  • Privacy Policy
Facebook  Twitter  LinkedIn  YouTube  RSS

Articles

  • Online First
  • Current Issue
  • Past Issues
  • Rapid Impact Archive
  • Meeting Abstracts

Information for

  • Authors
  • Subscribers
  • Advertisers
  • Librarians

About MCR

  • About the Journal
  • Editorial Board
  • Permissions
  • Submit a Manuscript
AACR logo

Copyright © 2021 by the American Association for Cancer Research.

Molecular Cancer Research
eISSN: 1557-3125
ISSN: 1541-7786

Advertisement