EZH2 Inhibition Sensitizes IDH1R132H-Mutant Gliomas to Histone Deacetylase Inhibitor

Isocitrate Dehydrogenase-1 (IDH1) is commonly mutated in lower-grade diffuse gliomas. The IDH1R132H mutation is an important diagnostic tool for tumor diagnosis and prognosis; however, its role in glioma development, and its impact on response to therapy, is not fully understood. We developed a murine model of proneural IDH1R132H-mutated glioma that shows elevated production of 2-hydroxyglutarate (2-HG) and increased trimethylation of lysine residue K27 on histone H3 (H3K27me3) compared to IDH1 wild-type tumors. We found that using Tazemetostat to inhibit the methyltransferase for H3K27, Enhancer of Zeste 2 (EZH2), reduced H3K27me3 levels and increased acetylation on H3K27. We also found that, although the histone deacetylase inhibitor (HDACi) Panobinostat was less cytotoxic in IDH1R132H-mutated cells (either isolated from murine glioma or oligodendrocyte progenitor cells infected in vitro with a retrovirus expressing IDH1R132H) compared to IDH1-wild-type cells, combination treatment with Tazemetostat is synergistic in both mutant and wild-type models. These findings indicate a novel therapeutic strategy for IDH1-mutated gliomas that targets the specific epigenetic alteration in these tumors.


Introduction
Diffuse gliomas are the most common type of primary brain tumor and include several biologically distinct tumor types with different molecular and genetic features.The current WHO classification system for CNS tumors incorporates common genetic alterations, and diffuse gliomas are classified based on the presence or absence of mutations in the metabolic enzyme Isocitrate Dehydrogenase (IDH), the most common being IDH1-R132H [1].Glioblastoma (GBM), the most common and malignant primary brain tumor, is IDH-wild-type, whereas IDH-mutant gliomas are classified as oligodendroglioma or astrocytoma.The majority of IDH-mutant gliomas have a proneural transcriptional phenotype; the gene expression profiles of these cells resemble oligodendrocyte progenitor cells (OPCs) or neural progenitor cells (NPCs) [2][3][4].Therefore, OPCs are considered a likely cell of origin for IDH-mutant gliomas [5][6][7][8].While the cell of origin in gliomagenesis is still disputed, the proliferation ability, distribution, and abundance in the adult brain further implicate OPCs in this process [9][10][11][12][13][14].Our lab and others have established the tumorigenic potential of OPCs residing in the adult white matter of mice and rats [6,7,10,[15][16][17][18].Although the Cells 2024, 13, 219 2 of 14 transforming effects of mutations in IDH have been extensively investigated, the specific effects of mutant IDH1 in the cellular transformation of OPCs have not been established.The connection between tumor genotype and phenotype suggests that the transforming effects of IDH mutation are highly dependent on the cellular context and that OPCs are uniquely susceptible, although the mechanism is not understood.
The ability of OPCs to acquire the expression profile of myelinating oligodendrocytes requires the repression of several transcripts related to other lineages, which is initiated by deacetylation and methylation of specific lysine residues on histone H3 [19][20][21].Enhancer of Zeste 2 (EZH2) is the main writer for H3K27 di-and trimethylation [22] and has a role in OPC lineage specification and proliferation as well as restriction of neuronal and astrocytic lineage commitment [23].Thus, the intrinsic functional necessity for epigenetic regulation and maintenance of myelination in the adult brain imparts OPCs with an innate susceptibility to metabolic stress.

Materials and Methods
Information about antibodies used and concentrations can be found in Supplementary Table S3.
Retrovirus construction: PDGF-IRES-Cre was generated by cloning human PDGF-B and Cre into a pQXIX vector (Clontech, Mountain View, CA, USA) as previously described [6].IDH1R132H-IRES-CRE and PDGF-IRES-GFP retroviruses were generated by cloning the human IDH1R132H sequence and Cre or PDGF-B and GFP into the same vector.A VSVG pseudotyped retrovirus was generated as previously described [6].Briefly, 293GP2 cells were incubated with plasmids expressing VSVG, and each viral plasmid was mixed with CaCl 2 and Hepes-buffered saline.Retroviruses were harvested and filtered and then resuspended in PBS.
Intracerebral stereotactic injections: Brain surgery procedures were performed in adult mixed background transgenic mice harboring floxed Trp53 and stop-floxed mCherry luciferase (gift from Thomas Ludwig).For subcortical white matter targeting, we used the coordinates of 2.0 mm lateral, 2.0 mm rostral, and 2.0 mm deep, with bregma as the reference point.Using a stereotaxis platform, 1 µL of the specified retrovirus was injected using a Hamilton syringe (flow rate 0.2 µL/min) at a depth of 2 mm.For serial transplantation experiments, 50 × 10 3 primary tumor cells were resuspended in Opti-MEM (Invitrogen, Waltham, MA, USA) and injected at the same coordinate.All injections were performed in mice between six and eight weeks of age.Overexpression of the PDGF-BB mitogen, combined with the deletion of the tumor-suppressor gene, forms brain tumors with the histopathologic and molecular features of proneural glioblastomas.Tumor progression was monitored using Xenogen IVIS Spectrum imager 10-15 min after 100 µL intraperitoneal injection of 30 mg/mL luciferin (Caliper Life Sciences, Hopkinton, MA, USA) with 1-min exposure.
MALDI imaging mass spectrometry by MS/MS: At the tumor end-stage, the whole brains of mice with IDH1R132H and IDH1WT glioma conditions were harvested and coronally bisected at the retroviral infusion site and posterior to the tumor.The anterior Cells 2024, 13, 219 3 of 14 section was used for H&E and IDH1R132H immunohistochemical staining.The coronal tumor-bearing slab was sliced at 12 microns and used for acquiring two-dimensional ion density images for glutamic acid (m/z 146-28) and 2-hydroxyglutaric acid (m/z 147-129) on a MALDI LTQ XL linear ion trap mass spectrometer at a spatial resolution of 50 microns according to methods previously described [40].Solution standards of 2-hydroxyglutaric acid concentrations were used for generating a standard calibration curve (intensity ratio of 2-HG/ internal standard vs. concentration of 2-HG (pmol/µL) using weighted linear regression (1/X 2 ).
Primary tumor cell lines: Primary tumor cell lines were derived as described previously [6,43,44].Once mice showed evidence of neurological features consistent with tumor morbidity, mice were sacrificed by cervical dislocation and decapitation.Brains were acutely harvested and maintained in prewarmed PBS (37 • C).Tumors were surgically resected, minced in 2.5% Trypsin, and serially triturated with a syringe using 18 G and 21 G needles.Triturated tumor tissue was filtered (70 µm) then neutralized in 50% FBS and centrifuged at 1500 rpm for 5 min.The supernatant was discarded, and the tissue was cultured in defined tissue culture media containing high glucose DMEM (5.4 g/L D-glucose, L-glutamine), 0.5% FBS, N2 supplement, PDGF-AA (10 ng/mL), bFGF (10 ng/mL), and antibiotic-antimycotic solution.
Drug treatment cell viability assay: Cells were plated in black-walled 96-well plates and allowed to settle for 24 h.For combined treatment studies, cells were pretreated with Tazemetostat (5 µM) or DMSO for 24 h and treated by doubling serial Panobinostat dilutions (1-128 nM) with Tazemetostat or DMSO for 72 h.CellTiter-Glo (Promega, Madison, WI, USA) was added to cells following the manufacturer's protocol.The luminescence was read on a GloMax microplate reader (Promega) using the provided manufacturer's protocol.
Drug combination treatments: For Western blot and flow cytometry studies, cells were plated and allowed to settle for 24 h.Cells were pretreated with Tazemetostat (5 µM) or DMSO for 24 h.Cells were then co-treated with the addition of Panobinostat (20 µM) or DMSO for an additional 24 h, providing four groups: control (DMSO alone), Tazemetostat (and DMSO), Panobinostat (and DMSO), and Tazemetostat and Panobinostat.
Whole-cell lysis: Cells were washed with ice-cold PBS and then scraped in Cell Extraction Buffer (Invitrogen).Cells were moved to Eppendorf tubes and kept on ice for lysis.Tubes were sonicated every 10 min for 30 min.Samples were spun at 13,000 rpm for 10 min at 4 • C, and the lysate was transferred to new tubes and stored at −80 • C until use.
Subcellular protein fractionation from cultured cells: Cells were lysed in hypotonic buffer (10 mM HEPES, pH 7.9, 1.5 mM MgCl 2 ,10 mM KCl) supplemented with 0.5 mM dithiothreitol (DTT), 1 mM phenylmethylsulfonyl fluoride (PMSF), 5 mM sodium butyrate, phosphatase inhibitor cocktail, and protease inhibitor cocktail freshly prepared (at 4 • C for 15 min), followed by 0.5% NP40 treatment (vortex for 10 s) to disrupt the cell membranes.Lysates were then centrifuged at 1500× g for 10 min at 4 • C to separate the cytoplasmic components (supernatant) from the nuclei-enriched fractions (pellets).Volumes (0.11 µL/µL) of 10× cytoplasmic extraction buffer (0.3 M HEPES, 1.4 M KCl, and 30 mM MgCl 2 ) were added to the supernatant and sonicated for 30 s ON/OFF for 5 min at high power in Bioruptor (Diagenode, Denville, NJ, USA).After centrifugation at 16,000× g for 10 min at 4 • C, the soluble fraction was collected as a cytoplasmic extract, and the pellet as the nuclear fraction.
Differential centrifugation subcellular protein fractionation from brain homogenate: Tissue samples were homogenized in 50 mM HEPES, 125 mM NaCl, 0.1 M sucrose with fresh protease inhibitor cocktail.Samples were spun for 10 min at 1000× g at 4 • C; the supernatant was collected as the cytoplasmic fraction, and the pellet as the nuclear fraction.
Histone extraction and purification after fractionation.Histones were extracted from the nuclear fraction using the acid extraction method with all steps performed at 4 • C [45].The nuclear pellet was incubated for 1.5 h to overnight in 0.4 N H 2 SO 4 .After centrifugation at 16,000× g for 10 min, nuclear debris was removed, and acid-soluble histones were then precipitated using trichloroacetic acid for at least 30 min.Samples were spun for 10 min at 16,000× g, washed twice in acetone, and then resuspended in water.
Immunofluorescence and immunohistochemistry: Cells for immunofluorescence were washed with PBS and then fixed with 4% paraformaldehyde (PFA) for 10 min at room temperature, and then the membrane was permeabilized and blocked with 0.5% (vol/vol) Triton X-100 (Fisher Scientific, Hampton, NH, USA) and 5% normal goat serum for 1 h.Primary antibodies were applied overnight at 4 • C followed by incubation of appropriate secondary antibodies conjugated with fluorophores.Images were captured using a Nikon TE2000 epifluorescence microscope equipped with MetaMorph software (Version.V.7.6.5.0,Molecular Devices, San Jose, CA, USA).Quantification of the immunofluorescent intensity was performed using ImageJ (version 1.53).Tumor-bearing mice were anesthetized and perfused with 4% PFA.After tissue processing and paraffin embedding, sections of 5 µm were cut by the Columbia University Medical Center Molecular Pathology Shared Resource facility.For immunohistochemistry, sections were de-paraffinized, immersed in 10 mM citrate buffer, pH 6.0, for 10 min in the microwave at 650 W, followed by blocking with 10% normal goat serum, before overnight incubation of primary antibodies at 4 • C. Appropriate secondary antibodies conjugated with fluorophores were used the following day to complete the staining.DAPI (4 ′ ,6-diamidino-2phenylindole) was used as a nuclear counterstain.
Microscopy: Stained tissue sections and fluorescent reporters of labeled tumor-derived cells were imaged using a Nikon TE2000 epifluorescence microscope equipped with Meta-Morph software (Version.V.7.6.5.0,Molecular Devices).Micrographs were processed and merged using MetaMorph and ImageJ.
Cellular dissociation for flow cytometry: Slice culture and brain tumors were dissociated using filtered Carica Papain extract diluted 1:20 in 1 N NaOH PBS (supplemented with 2 mg/mL L-cysteine) for 30 min at 37 • C with shaking and then centrifuged in excess PBS at 400× g for 5 min.The slices were manually dissociated with a pipette; then sucrose was added (final 15%), and the samples were spun at 1000 rpm for 5 min to remove debris.The remaining bottom ~30% was washed in PBS and then processed for flow cytometry.
Flow cytometry: Cells were collected and fixed with 4% formaldehyde (methanolfree) for 15 min at room temperature.Cells were centrifuged in excess PBS 400× g for 5 min for all washes.The cells were permeabilized by adding 100% ice-cold methanol drop-wise while gently vortexing cells to a final concentration of 90% methanol.Cells were incubated on ice for 30 min and then washed in excess PBS to remove methanol.Cells were resuspended in diluted primary antibodies and incubated for 1 h at room temperature.They were washed in excess PBS and then resuspended in diluted fluorochrome-conjugated secondary antibodies and incubated for 30 min at room temperature.Cells were washed again in excess PBS and then resuspended in PBS for flow cytometry reading.The BD LSRFortessa was used for optical measurements, and analysis was performed using FlowJo version 9 (Ashland, OR, USA: Becton, Dickinson and Company; 2019).
RNA extraction and sequencing: Total mRNA was extracted using the Qiagen RNeasy Mini kit (Qiagen, Hilden, Germany).The Columbia Sulzberger Genome Center sequenced the RNA samples (mm10, Illumina iGenomes, NIH HPC Biowulf cluster, Bethesda, MD, USA, http://hpc.nih.gov), and Tophat2 was used to map the reads.Differential expression analysis was performed using Cuffdiff2.We used expression data from the murine brain transcriptome database [46] to create cell-type-specific gene sets with the top 200 enriched genes, based on fold change above average expression, found in astrocytes, oligodendrocytes, neurons, OPCs, newly formed oligodendrocytes, myelinating oligodendrocytes, microglia, and endothelia.We performed Gene Set Enrichment Analysis [47] on the RNA sequencing data using these cell-type-specific gene sets.
Statistical analysis and graphing: The coefficient of drug interaction (CDI) was used to determine if drug-drug and drug-mutation interactions were synergistic (CDI < 1) or antagonistic (CDI > 1) [48,49].The CDI was calculated using CDI = AB/(A × B), where AB is the OD value ratio of the combination group to the control group and A or B is the OD value ratio of a single agent to the control group.All graphs with statistical analysis were made in Graphpad Prism version 8.0.0 for Mac (GraphPad Software, San Diego, California USA, www.graphpad.com).All volcano plots were made and all statistical analyses were performed using R [50].Heatmaps were made using R or Multi-experiment viewer (MeV) version 10.2 [51].Western blots were analyzed and quantified using Image Studio Lite analysis software V. 5.2 (LI-COR Biosciences, Lincoln, NE, USA).

Retroviral Delivery of IDH1R132H Mutation with PDGF and p53 Deletion Induces Glioma Formation with Increased 2HG
To determine the influence of IDH1 mutation on gliomagenesis induced by loss of p53 and overexpression of PDGF, mice with floxed p53 were stereotactically injected to target the glial progenitors of the subcortical white matter with a retrovirus expressing PDGFB and Cre-recombinase, alone or in combination with an IDH1R132H-expressing retrovirus.These two models allow us to compare tumorigenesis and phenotype between wild-type and IDH1R132H OPC transformation while controlling for other genetic alterations.Tumor growth was monitored by luciferase imaging, but no significant difference was seen in the growth kinetics.Both IDH1-mutant and wild-type murine gliomas reached the tumor morbidity end stage with a median survival of 40 days post-injection (Figure 1A), suggesting the IDH1 mutation did not facilitate or inhibit gliomagenesis in this co-delivery model.Tumors collected at morbidity showed large diffusely infiltrating lesions expanding from the site of injection (Figure 1B).The histology showed that both the wild-type and IDH1R132H tumors resembled malignant high-grade human gliomas.The IDH1-mutant tumors recapitulated the proneural glioma while retaining the mutation (Figure 1B and Supplementary Figure S1A).Furthermore, the mutation did not significantly affect the proliferation index of the fully formed gliomas (Supplementary Figure S1B), consistent with the observation that the IDH1R132H and wild-type glioma models have a similar disease progression.
IDH1 normally converts isocitrate to α-KG; however, the R132H mutation changes the enzymatic properties, causing a higher affinity for α-KG rather than isocitrate, thus producing 2HG [52,53].Low levels of 2HG are seen under physiological conditions; however, IDH1R132H produces elevated levels which have metabolic and epigenetic influences [54].We used matrix-assisted laser desorption ionization imaging mass spectroscopy (MALDI MS) to determine if IDH1R132H-expressing tumors had increased 2HG compared to normal brain and wild-type IDH1 tumor tissue [40].IDH1R132H-mutant gliomas showed increased 2HG in the tumor region (Figure 1C,D).We also analyzed α-KG and glutamate levels but found no significant difference in mutant tumors compared to the wild-type or normal brain (Supplementary Figure S1C,D).and glutamate levels but found no significant difference in mutant tumors compared to the wild-type or normal brain (Supplementary Figure S1C,D).

ID1R132H Mutation Alters the Transcriptional Phenotype of Mouse Gliomas
We performed transcriptional profiling by RNA sequencing on the IDH1R132H-mutant and wild-type tumors to determine if these tumors would recapitulate the molecular profile of human proneural gliomas.We compared the expression profile of the wild-type (N = 5) and IDH1R132H (N= 6) murine tumors to the TCGA database using Gene Set Enrichment Analysis (GSEA), with focused analysis on the Verhaak subtype classifier genes: Proneural (PN), Classical (CL), Mesenchymal (MES), Neural-high (NA) and Neural-low (NL) [4,6,47].We found that both the IDH1R132H-mutant tumors and the wild-type tumors were highly enriched in the PN molecular signature and depleted in the CL and MES signatures (Figure 2A, raw data can be found in Supplementary Table S1 and Supplementary Table S2).
We further analyzed the transcriptional differences between the IDH-mutant and IDH-WT mouse gliomas by performing GSEA using the Barres murine brain transcriptome database [46] to create cell-type-specific gene sets for astrocytes, oligodendrocytes, neurons, OPCs, newly formed oligodendrocytes, myelinating oligodendrocytes, microglia, and endothelia for analysis.Analysis of the gene sets showed that the IDH1R132H tumors were enriched in OPC genes (Figure 2B), while the IDH-wild-type tumors were significantly enriched in neuronal, astrocytic, myelinating oligodendroglial, and microglial genes (Figure 2C).These findings suggested that the IDH1R132H mutation fostered an OPC-restricted transcriptional phenotype during gliomagenesis by inhibiting the expression of other neural lineage genes.

ID1R132H Mutation Alters the Transcriptional Phenotype of Mouse Gliomas
We performed transcriptional profiling by RNA sequencing on the IDH1R132Hmutant and wild-type tumors to determine if these tumors would recapitulate the molecular profile of human proneural gliomas.We compared the expression profile of the wild-type (N = 5) and IDH1R132H (N= 6) murine tumors to the TCGA database using Gene Set Enrichment Analysis (GSEA), with focused analysis on the Verhaak subtype classifier genes: Proneural (PN), Classical (CL), Mesenchymal (MES), Neural-high (NA) and Neural-low (NL) [4,6,47].We found that both the IDH1R132H-mutant tumors and the wild-type tumors were highly enriched in the PN molecular signature and depleted in the CL and MES signatures (Figure 2A, raw data can be found in Supplementary Tables S1 and S2).
We further analyzed the transcriptional differences between the IDH-mutant and IDH-WT mouse gliomas by performing GSEA using the Barres murine brain transcriptome database [46] to create cell-type-specific gene sets for astrocytes, oligodendrocytes, neurons, OPCs, newly formed oligodendrocytes, myelinating oligodendrocytes, microglia, and endothelia for analysis.Analysis of the gene sets showed that the IDH1R132H tumors were enriched in OPC genes (Figure 2B), while the IDH-wild-type tumors were significantly enriched in neuronal, astrocytic, myelinating oligodendroglial, and microglial genes (Figure 2C).These findings suggested that the IDH1R132H mutation fostered an OPCrestricted transcriptional phenotype during gliomagenesis by inhibiting the expression of other neural lineage genes.B,C).GSEA of significant differentially expressed genes using cell-specific brain transcriptome gene sets.Gene sets were made with the top 200 enriched genes based on fold change above average for each cell type: astrocytes, oligodendrocytes, neurons, OPCs, newly formed oligodendrocytes, myelinating oligodendrocytes, microglia, and endothelia.Only showing OPC lineage (B) and significant GSEA (neuron, astrocytes, myelinating oligodendrocytes, and microglia) (C).

Tazemetostat Sensitizes IDH1R132H Cells to Panobinostat, and Combined Treatment Increases H3K27ac
Based on the increased 2HG levels in the IDH1R132H-mutant tumors as well as the transcriptomic data, we wanted to examine if H3K27me3 levels differed between IDH1 mutant and wild-type mouse glioma cells.Western blot analysis of histone extracts from IDH1R132H cells revealed significantly higher H3K27me3 levels compared to IDH1 wildtype cells (Figure 3A,B).
Based on the differences in post-translational histone modifications, we also compared the viability of our cells after treatment with epigenome-modifying drugs.We found that the wild-type and IDH1R132H mouse glioma cells were equally insensitive to treatment after 72 h with EZH2 inhibitor Tazemetostat and did not reach IC50 even at 64 µM (Supplementary Figure S1E), although the treatment was effective in reducing H3K27me3 at doses as low as 0.03 µM (Supplementary Figure S1F).
Notably, Tazemetostat (5 µM) significantly reduced H3K27me3 in both IDH1R132H and wild-type cells by 48 h (Figure 3A,B) and significantly increased H3K27ac in IDH1R132H cells by 48 h, and in wild-type cells by 116 h (Supplementary Figure S1G), revealing that Tazemetostat has reciprocal effects on methylation and acetylation at H3K27 in both cell lines.
We used a histone deacetylase inhibitor, Panobinostat (20 µM for 24 h), to determine if co-treatment with Tazemetostat would synergistically increase H3K27ac.We acutely treated glioma cells with low drug doses and assessed the presence of post-translational modifications (Figure 3A).Most notably, co-treatment in IDH1WT and IDH1R132H glioma cells revealed significantly higher levels of H3K27ac compared to vehicle-treated control cells (Figure 3C).We further demonstrated the synergistic effects of combined treatment on the H3K27 marks via flow cytometry by comparing across treatment conditions the numbers of cells above a specified immunofluorescence threshold for each histone mark (Supplementary Figure S2C,D).
Based on these findings, we asked whether a combination treatment with Tazemetostat and Panobinostat would be more effective in inducing cytotoxicity in glioma cells.Based on a 72-h viability assay, the IDH1R132H (IC50 = 47.64)glioma cells were less sensitive to Panobinostat treatment than wild-type (IC50 = 20.04)glioma cells.However, 24-h pretreatment with Tazemetostat followed by co-treatment increased the Panobinostat sensitivity of IDH1R132H cells to levels similar to those seen with wild-type cells (WT IC50 = 17.65 and IDH1R132H IC50 = 22.95) (Figure 3D and Supplementary Figure S2A,B).We used the coefficient of drug interaction (CDI) to determine the nature of drug coaction, in which a CDI value of <1 indicates synergy and a CDI > 1 indicates antagonism.We determined that IDH1R132H mutation and Panobinostat are antagonistic (CDI = 1.52) (Figure 3F).These data show that there is synergy between Tazemetostat and Panobinostat in both cell lines (WT CDI = 0.61 and IDH1R132H CDI = 0.64) (Figure 3E).The glioma cells harboring the IDH1R132H mutation were resistant to Panobinostat, but this is rectified by EZH2 inhibition.This suggests that the synergistic mechanism of action of the combined treatment is to block methylation, allowing for more acetylation than either single treatment alone.
treatment is to block methylation, allowing for more acetylation than either single treatment alone.

Combined Treatment in OPCs Recapitulates the Metabolic and Epigenetic Alterations Found in IDH1-Mutated Gliomas
Given that the IDH1R132H mutation is observed in the early stages of glioma formation, we hypothesize that this mutation and the metabolic alterations that it causes have their primary effect on the glial progenitor cells that give rise to gliomas.To test this hypothesis OPCs were isolated from mCherry-luciferase stop-flox p53flox/flox mice and then infected in vitro with retroviruses expressing cre alone (X-cre) to delete p53 or cre and IDH1R132H (IDH1R132H-cre).

Combined Treatment in OPCs Recapitulates the Metabolic and Epigenetic Alterations Found in IDH1-Mutated Gliomas
Given that the IDH1R132H mutation is observed in the early stages of glioma formation, we hypothesize that this mutation and the metabolic alterations that it causes have their primary effect on the glial progenitor cells that give rise to gliomas.To test this hypothesis OPCs were isolated from mCherry-luciferase stop-flox p53flox/flox mice and then infected in vitro with retroviruses expressing cre alone (X-cre) to delete p53 or cre and IDH1R132H (IDH1R132H-cre).
To examine the dynamics of H3K27 regulation in the early stages of gliomagenesis, we assessed the epigenetic changes that EZH2i and HDACi combination treatment had on early-passage (p < 6) p53-deleted OPCs (Figure 4A).IDH1R132H-cre OPCs had higher baseline trimethylation than wild-type X-cre OPCs, similar to what we observed in glioma cells derived from our fully formed IDH1R132H murine tumor model.Tazemetostat monotherapy decreased trimethylation for both retroviral cell conditions (Figure 4B).Both IDH1R132H-cre and X-cre OPCs also had a marked increase in H3K27ac after Tazemetostat and Panobinostat co-treatment, with a 66.2-and 55-fold change for wild-type and IDH1R132H, respectively (Figure 4C).We further demonstrated these findings using flow cytometry (Supplementary Figure S3C,D).
To examine the dynamics of H3K27 regulation in the early stages of gliomagenesis, we assessed the epigenetic changes that EZH2i and HDACi combination treatment had on early-passage (p < 6) p53-deleted OPCs (Figure 4A).IDH1R132H-cre OPCs had higher baseline trimethylation than wild-type X-cre OPCs, similar to what we observed in glioma cells derived from our fully formed IDH1R132H murine tumor model.Tazemetostat monotherapy decreased trimethylation for both retroviral cell conditions (Figure 4B).Both IDH1R132H-cre and X-cre OPCs also had a marked increase in H3K27ac after Tazemetostat and Panobinostat co-treatment, with a 66.2-and 55-fold change for wild-type and IDH1R132H, respectively (Figure 4C).We further demonstrated these findings using flow cytometry (Supplementary Figure S3C,D).To determine if Tazemetostat and Panobinostat have a synergistic effect on retrovirally transformed OPCs, we repeated the combination treatment survival study in wildtype-and IDH1R132H-mutated p53-deleted OPCs.The IDH1R132H-mutated OPCs were less sensitive to Panobinostat (Figure 4D and Supplementary Figure S3A,B), further To determine if Tazemetostat and Panobinostat have a synergistic effect on retrovirally transformed OPCs, we repeated the combination treatment survival study in wild-typeand IDH1R132H-mutated p53-deleted OPCs.The IDH1R132H-mutated OPCs were less sensitive to Panobinostat (Figure 4D and Supplementary Figure S3A,B), further suggesting mutated IDH1 protects cells from Panobinostat-triggered cell death, despite the higher levels of histone acetylation (WT IC50 = 15.72 vs. IDH1R132H IC50 = 23.65).Wild-type and IDH1R132H-mutant OPCs were more sensitive to the co-treatment than Panobinostat alone (WT IC50 = 11.5 and IDH1R132H IC50 = 13.28, Figure 4D and Supplementary Figure S3A,B).Both the wild-type (CDI = 0.42) and IDH1-mutant (CDI = 0.49) OPCs showed synergy be-tween Tazemetostat and Panobinostat treatment (Figure 4E).We also found that IDH1R132H is antagonistic to Panobinostat treatment (CDI = 1.3) (Figure 4F), further suggesting that IDH1R132H protects against the cytotoxic effects of Panobinostat in OPCs.These data demonstrate that the IDH1R132H mutation makes retrovirally transformed OPCs more resistant to Panobinostat but that combined treatment with Tazemetostat sensitizes cells to Panobinostat.
Together, these results demonstrate that Tazemetostat and Panobinostat co-treatment is synergistic in increasing H3K27ac and reducing cell viability more than either monotherapy in glioma cells and OPCs.Furthermore, the similarities between OPCs and glioma cells suggest that the epigenetic vulnerabilities are inherited from the cell of origin.

Discussion
In this study, we designed an approach to compare the epigenetic and transcriptional phenotype of IDH1R132H and IDH1 wild-type mouse glioma models and how early stages of OPC gliomagenesis may provide insight into mechanisms of glial lineage specification that can be therapeutically exploited.We find that while IDH1WT and IDH1R132H tumors had comparable growth dynamics, IDH1R132H recapitulated superphysiological levels of intratumoral alpha-2 hydroxyglutarate and acquired a different transcriptional phenotype characterized by similar levels of OPC-related genes and significantly lower levels of expression for genes associated with neurons, astrocytes, oligodendroglia and microglia.This suggests that IDH1R132H mutation and the associated metabolic and epigenetic alterations restricted cell transformation to an OPC-like phenotype during gliomagenesis.
Although Tazemetostat monotherapy treatment was not cytotoxic in our murine glioma model, it decreased H3K27me3 and increased H3K27ac, which are histone modifications necessary for OPC lineage restriction [55].To exploit the cross-talk between acetylation and trimethylation on H3K27, we treated tumor-derived cells with both Tazemetostat and the pan-HDAC inhibitor Panobinostat.We hypothesized that combining an HDAC inhibitor and an EZH2 inhibitor would be synergistic by targeting epigenetic mechanisms (i.e., histone deacetylase and histone methyltransferase), whereby EZH2 inhibitors remove H3K27me3 to allow for acetylation of H3K27, which could be beneficial to proneural glioma patient prognosis [56].Others have found that the combination of EZH2 inhibition with an HDACi is synergistic in other types of cancer [57][58][59][60]; however, the effect has previously not been tested in IDH1R132H glioma.Both tumor-derived glioma cells and transformed OPCs show similar responses to EZH2 and HDAC inhibition, suggesting that this mechanism is retained throughout the process of OPC gliomagenesis.Using both Tazemetostat and Panobinostat increased H3K27ac more than either monotherapy, and the co-treatment was more cytotoxic to the cells than monotherapy specifically in IDH1-mutated cells.
We found that IDH1R132H in combination with other genetic alterations transformed OPCs to proneural glioma.A detailed understanding of epigenetic alterations provides a therapeutic target to leverage.Understanding the cell of origin and how genetic alterations play out within that cellular context can help in understanding therapeutic effects and developing new therapeutic strategies.We found that IDH1R132H-mutant murine tumors had increased H3K27me3 which rendered those cells resistant to Panobinostat treatment.Treating cells with an EZH2 inhibitor removed H3K27me3 to allow for Panobinostat to be more effective.It is important to note that our model retains both endogenous copies of wild-type IDH1 with IDH1-mutant overexpression, which we acknowledge does not recapitulate physiological gene ratios in human IDH1-mutant tumors.We also looked at global epigenetic changes and not the transcriptional outcome of these changes.Future studies should determine the transcriptional effects of using these drugs in the more clinically relevant setting of heterozygous IDH1R132H mutation.These results show that although an epigenome-modifying drug might not be effective as a monotherapy [61], it can still be beneficial in combination treatments for glioma.Epigenome-modifying drugs can be more effective in combination by targeting antagonistic modifications (such as methylation and acetylation).

Figure 1 .
Figure 1.IDH1R132H-mutant tumor model does not influence murine survival with increased 2HG.(A) Kaplan-Meier curve comparing the survival of WT and R132H mouse groups; log-rank (Mantel-Cox) test p-value = 0.16 and Gehan-Breslow Wilcoxon test p-value = 0.46.(B) WT and R132H tumors are histologically indistinguishable using H&E.Upper panel (1.5×) with 1 mm scale bar and lower panel (10×) with 100 μm scale bar.(C) Left panel of H&E staining to show tumor margins.Right panels of MALDI MS two-dimensional ion density images for 2-hydroxyglutaric acid collected at a spatial resolution of 200 μm.(D) Quantification of average concentration (nmol/mg) of 2-HG across tumor (TT) and non-tumor (NT) regions.One-way ANOVA with Tukey's correction for multiple comparison was used.**** p-value < 0.0001.

Figure 1 .
Figure 1.IDH1R132H-mutant tumor model does not influence murine survival with increased 2HG.(A) Kaplan-Meier curve comparing the survival of WT and R132H mouse groups; log-rank (Mantel-Cox) test p-value = 0.16 and Gehan-Breslow Wilcoxon test p-value = 0.46.(B) WT and R132H tumors are histologically indistinguishable using H&E.Upper panel (1.5×) with 1 mm scale bar and lower panel (10×) with 100 µm scale bar.(C) Left panel of H&E staining to show tumor margins.Right panels of MALDI MS two-dimensional ion density images for 2-hydroxyglutaric acid collected at a spatial resolution of 200 µm.(D) Quantification of average concentration (nmol/mg) of 2-HG across tumor (TT) and non-tumor (NT) regions.One-way ANOVA with Tukey's correction for multiple comparison was used.**** p-value < 0.0001.

Figure 2 .
Figure 2. RNA sequencing reveals that wild-type and IDH1R132H murine tumors are proneural and mutant gliomas are more lineage-restricted.(A) Heatmap of Gene Set Enrichment Analysis (GSEA) used to compare the RNA sequence profiles of IDH1R132H-mutant and WT murine gliomas to the

Figure 3 .
Figure 3. Tazemetostat and Panobinostat co-treatment is synergistic in glioma cells.(A) Representative Western blot of wild-type and IDH1R132H glioma cells pretreated with 5 μM Tazemetostat (EZH2i) or DMSO for 24 h then retreated with the addition of 20 μM Panobinostat (HDACi) or DMSO for 24 h.Cells were fractionated for the cytoplasmic fraction (CF), and the nuclear fraction (NF) was acid-extracted.Quantification of (B) H3K27me3 and (C) H3K27ac normalized to total H3.All bar graphs are from three independent experiments and shown as fold change to control.(D) Area under the curve analysis from three independent viability experiments normalized to wildtype cells.Paired T-test performed for statistical analysis.(E) Coefficient of drug interaction from three independent experiments looking at Panobinostat and Tazemetostat synergy in wild-type and IDH1R132H cells.(F) Coefficient of drug interaction from three independent experiments looking at Panobinostat and IDH1R132H.All bar graphs are from at least three independent experiments and shown as fold change to control.Student's T-test used for statistical analysis * p-value ≤ 0.05, ** p-value ≤ 0.01.

Figure 3 .
Figure 3. Tazemetostat and Panobinostat co-treatment is synergistic in glioma cells.(A) Representative Western blot of wild-type and IDH1R132H glioma cells pretreated with 5 µM Tazemetostat (EZH2i) or DMSO for 24 h then retreated with the addition of 20 µM Panobinostat (HDACi) or DMSO for 24 h.Cells were fractionated for the cytoplasmic fraction (CF), and the nuclear fraction (NF) was acid-extracted.Quantification of (B) H3K27me3 and (C) H3K27ac normalized to total H3.All bar graphs are from three independent experiments and shown as fold change to control.(D) Area under the curve analysis from three independent viability experiments normalized to wild-type cells.Paired T-test performed for statistical analysis.(E) Coefficient of drug interaction from three independent experiments looking at Panobinostat and Tazemetostat synergy in wild-type and IDH1R132H cells.(F) Coefficient of drug interaction from three independent experiments looking at Panobinostat and IDH1R132H.All bar graphs are from at least three independent experiments and shown as fold change to control.Student's T-test used for statistical analysis * p-value ≤ 0.05, ** p-value ≤ 0.01.

Figure 4 .
Figure 4. Tazemetostat and Panobinostat co-treatment is synergistic in OPCs.(A) Representative Western blot of X-cre and IDH1R132H-cre OPCs pretreated with 5 μM Tazemetostat (EZH2i) or DMSO for 24 h then retreated with the addition of 20 μM Panobinostat (HDACi) or DMSO for 24 h.Cells were fractionated for the cytoplasmic fraction (CF), and the nuclear fraction (NF) was acidextracted.Quantification of (B) H3K27me3 (C) and H3K27ac normalized to total H3.(D) Area under the curve analysis from three independent viability experiments normalized to wild-type cells.Paired T-test performed for statistical analysis.(E) Coefficient of drug interaction from three independent experiments looking at Tazemetostat and Panobinostat synergy in X-cre and IDH1R132Hcre infected OPCs.(F) Coefficient of drug interaction from three independent experiments looking at Panobinostat and IDH1R132H.All bar graphs are from at least three independent experiments and shown as fold change to control.Student's T-test used for statistical analysis.* p-value ≤ 0.05, ** p-value ≤ 0.01 and *** p-value ≤ 0.001.

Figure 4 .
Figure 4. Tazemetostat and Panobinostat co-treatment is synergistic in OPCs.(A) Representative Western blot of X-cre and IDH1R132H-cre OPCs pretreated with 5 µM Tazemetostat (EZH2i) or DMSO for 24 h then retreated with the addition of 20 µM Panobinostat (HDACi) or DMSO for 24 h.Cells were fractionated for the cytoplasmic fraction (CF), and the nuclear fraction (NF) was acidextracted.Quantification of (B) H3K27me3 (C) and H3K27ac normalized to total H3.(D) Area under the curve analysis from three independent viability experiments normalized to wild-type cells.Paired T-test performed for statistical analysis.(E) Coefficient of drug interaction from three independent experiments looking at Tazemetostat and Panobinostat synergy in X-cre and IDH1R132H-cre infected OPCs.(F) Coefficient of drug interaction from three independent experiments looking at Panobinostat and IDH1R132H.All bar graphs are from at least three independent experiments and shown as fold change to control.Student's T-test used for statistical analysis.* p-value ≤ 0.05, ** p-value ≤ 0.01 and *** p-value ≤ 0.001.