Introduction: Low-grade and secondary high-grade gliomas frequently contain a mutation in the IDH1 metabolic enzyme that is hypothesized to drive tumorigenesis by inhibiting many of the chromatin-regulating enzymes that regulate DNA structure. Much of the previous research has focused on DNA methylation leaving histone modifications relatively under-studied. Histone deacetylase inhibitors are promising anti-cancer agents and have already been used in clinical trials, however a clear understanding of their mechanism or gene targets is lacking. In this study, the authors genetically dissect patient-derived IDH1 mutant cultures to determine which HDAC enzymes drive growth in IDH1 mutant gliomas and what are the down-stream gene targets. This information will allow clinicians to design more specific HDAC inhibitors as well as monitor for treatment response. Methods: A panel of patient-derived gliomasphere cell lines (2 IDH1 mutant lines, 3 IDH1 wildtype lines) were subjected to a drug-screen of 110 epigenetic modifying drugs from 34 different epigenetic classes. The effect of LBH (Panobinostat) on gene expression and chromatin structure was tested on patient-derived IDH1 mutant lines. The role of each of the highly expressed HDAC enzymes was molecularly dissected using lentiviral RNA interference knock-down vectors and a patient-derived IDH1 mutant in vitro model of glioblastoma (HK252). These results were then confirmed in an in vivo xenotransplant model (BT-142). Results: The IDH1 mutation leads to gene down-regulation, DNA hypermethylation, increased DNA accessibility and H3K27 hypo-acetylation in two distinct IDH1 mutant over-expression models. The drug screen identified histone deacetylase inhibitors (HDACi) and Panobinostat (LBH) more specifically as the most selective compounds to inhibit growth in IDH1 mutant glioma lines. Of the eleven annotated HDAC enzymes (HDAC1-11) only six are expressed in IDH1 mutant glioma tissue samples and patient-derived gliomasphere lines (HDAC1, HDAC2, HDAC3, HDAC4, HDAC6, HDAC9). Lentiviral knock-down experiments revealed that HDAC1 and HDAC6 are the most consistently essential for growth both in vitro and in vivo and target very different gene modules. Knock-down of HDAC1 or HDAC6 in vivo led to a more circumscribed less invasive tumor. Conclusions: The gene dysregulation induced by the IDH1 mutation is wide-spread and only partially reversible by direct IDH1 inhibition. This study reports a previously under-described phenomenon of histone deacetylation induced by the IDH1 mutant enzyme and identifies HDAC1 and HDAC6 as important and drug-targetable enzymes that are necessary for growth and invasiveness in IDH1 mutant gliomas.

The pathological changes in epigenetics and gene regulation that accompany the progression of low-grade to high-grade gliomas are under-studied. The authors use a large set of paired atac-seq and RNA-seq data from surgically resected glioma specimens to infer gene regulatory relationships in glioma. Thirty-eight glioma patient samples underwent atac-seq sequencing and 16 samples underwent additional RNA-seq analysis. Using an atac-seq/RNA-seq correlation matrix, atac-seq peaks were paired with genes based on high correlation values (|r2| > 0.6). Samples clustered by IDH1 status but not by grade. Surprisingly there was a trend for IDH1 mutant samples to have more peaks. The majority of peaks are positively correlated with survival and positively correlated with gene expression. Constructing a model of the top six atac-seq peaks created a highly accurate survival prediction model (r2 = 0.68). Four of these peaks were still significant after controlling for age, grade, pathology, IDH1 status and gender. Grade II, III, and IV (primary) samples have similar transcription factors and gene modules. However, grade IV (recurrent) samples have strikingly few peaks. Patient-derived glioma cultures showed decreased peak counts following radiation indicating that this may be radiation-induced. This study supports the notion that IDH1 mutant and IDH1 wildtype gliomas have different epigenetic landscapes and that accessible chromatin sites mapped by atac-seq peaks tend to be positively correlated with expression. The data in this study leads to a new model of treatment response wherein glioma cells respond to radiation therapy by closing open regions of DNA.

Introduction: The pathological changes in epigenetics and gene regulation that accompany the progression of low-grade to high-grade gliomas are under-studied. The authors use a large set of paired atac-seq and RNA-seq data from surgically resected glioma specimens to infer gene regulatory relationships in glioma.Methods: Thirty-eight glioma patient samples underwent atac-seq sequencing and sixteen samples underwent additional RNA-seq analysis. Using an atac-seq/RNA-seq correlation matrix, atac-seq peaks were paired with genes based on high correlation values (|r2|>0.6). Results: Samples clustered by IDH1 status but not by grade. Surprisingly there was a trend for IDH1 mutant samples to have more peaks. The majority of peaks are positively correlated with survival and positively correlated with gene expression. Constructing a model of the top six atac-seq peaks created a highly accurate survival prediction model (r2=0.68). Four of these peaks were still significant after controlling for age, grade, pathology, IDH1 status and gender. Grade II, III, and IV (primary) samples have similar transcription factors and gene modules. However, grade IV (recurrent) samples have strikingly few peaks. Patient-derived glioma cultures showed decreased peak counts following radiation indicating that this may be radiation-induced.Conclusions: This study supports the notion that IDH1 mutant and IDH1 wildtype gliomas have different epigenetic landscapes and that accessible chromatin sites mapped by atac-seq peaks tend to be positively correlated with expression. The data in this study leads to a new model of treatment response wherein glioma cells respond to radiation therapy by closing open regions of DNA.