Search

Supplementary Material and Methods

Supplementary Table S1-S4

Isocitrate dehydrogenase-mutant low-grade gliomas (IDHmut-LGG) grow slowly but frequently undergo malignant transformation, which eventually leads to premature death. Chemotherapy and radiotherapy treatments prolong survival, but can also induce genetic (or epigenetic) alterations involved in transformation. Here, we developed a mathematical model of tumor progression based on serial tumor volume data and treatment history of 276 IDHmut-LGGs classified by chromosome 1p/19q codeletion (IDHmut/1p19qcodel and IDHmut/1p19qnoncodel) and performed genome-wide mutational analyses, including targeted sequencing and longitudinal whole-exome sequencing data. These analyses showed that tumor mutational burden correlated positively with malignant transformation rate, and chemotherapy and radiotherapy significantly suppressed tumor growth but increased malignant transformation rate per cell by 1.8 to 2.8 times compared with before treatment. This model revealed that prompt adjuvant chemoradiotherapy prolonged malignant transformation-free survival in small IDHmut-LGGs (≤ 50 cm3). Furthermore, optimal treatment differed according to genetic alterations for large IDHmut-LGGs (> 50 cm3); adjuvant therapies delayed malignant transformation in IDHmut/1p19qnoncodel but often accelerated it in IDHmut/1p19qcodel. Notably, PI3K mutation was not associated with malignant transformation but increased net postoperative proliferation rate and decreased malignant transformation-free survival, prompting the need for adjuvant therapy in IDHmut/1p19qcodel. Overall, this model uncovered therapeutic strategies that could prevent malignant transformation and, consequently, improve overall survival in patients with IDHmut-LGGs.Significance:A mathematical model successfully estimates malignant transformation-free survival and reveals a link between genetic alterations and progression, identifying precision medicine approaches for optimal treatment of IDH-mutant low-grade gliomas.

Supplementary Figures S1-S8

Glioblastoma (GBM), the most common primary brain tumor, is the most aggressive human cancers, with a median survival rate of only 14.6 months. Temozolomide (TMZ) is the frontline chemotherapeutic drug in GBM. Drug resistance is the predominant obstacle in TMZ therapy. Drug resistance occurs via multiple pathways such as DNA mismatch repair and base excision repair systems, by which glioma cells acquire chemoresistance to some extent (5% and 95%, respectively). Histone3 Lysin27 residue-acetylation (H3K27ac) status regulates cis-regulatory elements, which increases the likelihood of gene transcription. Histone deacetylase (HDAC) complex deacetylate lysine residues on core histones, leading to a decrease in gene transcription. In cis-regulatory element regions, complexes with HDAC repress histones by H3K27ac deacetylation. The cis-regulating and three-dimensional transcriptional mechanism is called "super-enhancer". RET finger protein (RFP) is a protein that is expressed in many kinds of cancer. RFP forms a protein complex with HDAC1. The disruption of the RFP-HDAC1 complex has resulted in increased drug sensitivity in other cancers. We conclude that the downregulation of RFP or the disruption of the RFP/HDAC1 complex leads to an increase in TMZ efficacy in glioblastoma by changing histone modifications which lead to changes in cell division, cell cycle and apoptosis.

Isocitrate dehydrogenase-mutant low-grade gliomas (IDHmut-LGG) grow slowly but frequently undergo malignant transformation, which eventually leads to premature death. Chemotherapy and radiotherapy treatments prolong survival, but can also induce genetic (or epigenetic) alterations involved in transformation. Here, we developed a mathematical model of tumor progression based on serial tumor volume data and treatment history of 276 IDHmut-LGGs classified by chromosome 1p/19q codeletion (IDHmut/1p19qcodel and IDHmut/1p19qnoncodel) and performed genome-wide mutational analyses, including targeted sequencing and longitudinal whole-exome sequencing data. These analyses showed that tumor mutational burden correlated positively with malignant transformation rate, and chemotherapy and radiotherapy significantly suppressed tumor growth but increased malignant transformation rate per cell by 1.8 to 2.8 times compared with before treatment. This model revealed that prompt adjuvant chemoradiotherapy prolonged malignant transformation-free survival in small IDHmut-LGGs (≤ 50 cm3). Furthermore, optimal treatment differed according to genetic alterations for large IDHmut-LGGs (> 50 cm3); adjuvant therapies delayed malignant transformation in IDHmut/1p19qnoncodel but often accelerated it in IDHmut/1p19qcodel. Notably, PI3K mutation was not associated with malignant transformation but increased net postoperative proliferation rate and decreased malignant transformation-free survival, prompting the need for adjuvant therapy in IDHmut/1p19qcodel. Overall, this model uncovered therapeutic strategies that could prevent malignant transformation and, consequently, improve overall survival in patients with IDHmut-LGGs. Significance: A mathematical model successfully estimates malignant transformation-free survival and reveals a link between genetic alterations and progression, identifying precision medicine approaches for optimal treatment of IDH-mutant low-grade gliomas.

Glioblastoma is one of the most malignant forms of cancer, for which no effective targeted therapy has been found. Although The Cancer Genome Atlas has provided a list of fusion genes in glioblastoma, their role in progression of glioblastoma remains largely unknown. To search for novel fusion genes, we obtained RNA-seq data from TGS-01 human glioma-initiating cells, and identified a novel fusion gene (HMGA2-EGFR), encoding a protein comprising the N-terminal region of the high-mobility group AT-hook protein 2 (HMGA2) fused to the C-terminal region of epidermal growth factor receptor (EGFR), which retained the transmembrane and kinase domains of the EGFR. This fusion gene product showed transforming potential and a high tumor-forming capacity in cell culture and in vivo. Mechanistically, HMGA2-EGFR constitutively induced a higher level of phosphorylated STAT5B than EGFRvIII, an in-frame exon deletion product of the EGFR gene that is commonly found in primary glioblastoma. Forced expression of HMGA2-EGFR enhanced orthotopic tumor formation of the U87MG human glioma cell line. Furthermore, the EGFR kinase inhibitor erlotinib blocked sphere formation of TGS-01 cells in culture and inhibited tumor formation in vivo. These findings suggest that, in addition to gene amplification and in-frame exon deletion, EGFR signaling can also be activated by gene fusion, suggesting a possible avenue for treatment of glioblastoma.

World Health Organization grade II and III gliomas most frequently occur in the central nervous system (CNS) in adults. Gliomas are not circumscribed; tumor edges are irregular and consist of tumor cells, normal brain tissue, and hyperplastic reactive glial cells. Therefore, the tumors are not fully resectable, resulting in recurrence, malignant progression, and eventual death. Approximately 69–80% of grade II and III gliomas harbor mutations in the isocitrate dehydrogenase 1 gene (IDH1), of which 83–90% are found to be the IDH1-R132H mutation. Detection of the IDH1-R132H mutation should help in the differential diagnosis of grade II and III gliomas from other types of CNS tumors and help determine the boundary between the tumor and normal brain tissue. In this study, we established a highly sensitive antibody-based device, referred to as the immuno-wall, to detect the IDH1-R132H mutation in gliomas. The immuno-wall causes an immunoreaction in microchannels fabricated using a photo-polymerizing polymer. This microdevice enables the analysis of the IDH1 status with a small sample within 15 min with substantially high sensitivity. Our results suggested that 10% content of the IDH1-R132H mutation in a sample of 0.33 μl volume, with 500 ng protein, or from 500 cells is theoretically sufficient for the analysis. The immuno-wall device will enable the rapid and highly sensitive detection of the IDH1-R132H mutation in routine clinical practice.

RET finger protein (RFP) forms a complex with histone deacetylase 1, resulting in aberrant deacetylation of H3K27ac and dysregulation of cis-regulatory elements. We evaluated the modulatory effects of RFP knockdown on cis-regulatory elements, gene expression, and chemosensitivity to temozolomide both in glioblastoma cells and in an intracranial glioblastoma model. The combination of RFP knockdown and temozolomide treatment markedly suppressed the glioblastoma cell growth due to oxidative stress and aberrant cell cycle and increased survival time in mice with glioblastoma. ChIP-seq and RNA-seq revealed that RFP knockdown increased or decreased activity of numerous cis-regulatory elements that lie adjacent to genes that control functions such as apoptosis, mitosis, DNA replication, and cell cycle: FOXO1, TBP2, and PARPBP. This study suggests that RFP contributes to chemoresistance via aberrant deacetylation of histone H3 at K27, whereas dysregulation of RFP-associated cis-regulatory elements in glioma and RFP knockdown combined with temozolomide is an effective treatment strategy for lethal glioma.

The authors recently found that 80% of lower-grade gliomas (LGGs) harbored a mutation in IDH1. Intraoperative detection of the mutated IDH1 helps not only differentiate LGGs from other type of brain tumors, but determine the resection border. In the current study, the authors have applied an automated genetic typing involving a quenching probe to detect the mutated IDH1. If tumor cells with the mutated IDH1 contained 10% or more in the mixture of normal and tumor cells, the device could detect it sensitively. The intraoperative assessment of IDH1 mutation is useful in brain tumor surgeries.

Comprising more than 80% of malignant brain tumors, glioma has proven to be a daunting cause of mortality in a vast majority of the human population. Progressive and extensive research on malignant glioma has substantially enhanced our understanding of glioma cell biology and molecular pathology. Subtypes of glioma such as astrocytoma and oligodendroglioma are currently grouped together into one pathological class, where they show many differences in histology and molecular etiology. This indicates that it may be beneficial to consider a new and radical subclassification. Thus, we summarize recent developments in glioblastoma multiforme (GBM) subtypes, immunohistochemical analyses useful for diagnoses and the biological evaluation and therapeutic implications of gliomas in this review.

RET finger protein (RFP, aka TRIM27) plays pivotal roles in acquirement of chemoresistance, via the formation of the complex with nuclear transcription factor Y (NF-Y) and histone deacetylase1 (HDAC1), which are associated with deacetylation of histone H3K27 and specific enhancer activity. Although the chemoresistance driven by RFP may result from aberrant deacetylation of histone H3K27 and dysregulation of enhancers, the details of the mechanism has not been studied so far. The aim of this study was to comprehensively investigate the effects of RFP on super-enhancers and gene expression. We also evaluated the effects of RFP knock-down (KD) on growth of GBM treated with temozolomide (TMZ) both in vitro and vivo. We found that the combination of RFP-KD and TMZ markedly suppressed growth of GBM cells and extended the survival time of mice bearing intracranial GBM compared to TMZ alone. Chromatin immunoprecipitation sequencing (ChIP-seq) using anti-H3K27ac antibody revealed that RFP-KD GBM cells aberrantly acquired several super-enhancers that were different from ones in RFP normal GBM cells. Whole-transcriptome sequencing (RNA-seq) confirmed overexpression of the corresponding genes in RFP-KD GBM cells. This study suggests that RFP contributes to chemoresistance via aberrant deacetylation of histone H3K27 and dysregulation of RFP-associated super-enhancers in GBM and that the combination of targeting RFP and TMZ might become an effective novel treatment strategy for lethal GBM.

Glioblastoma is one of the most malignant forms of cancer, for which no effective targeted therapy has been found. Although The Cancer Genome Atlas has provided a list of fusion genes in glioblastoma, their role in progression of glioblastoma remains largely unknown. To search for novel fusion genes, we obtained RNA-seq data from TGS-01 human glioma-initiating cells, and identified a novel fusion gene (HMGA2-EGFR), encoding a protein comprising the N-terminal region of the high-mobility group AT-hook protein 2 (HMGA2) fused to the C-terminal region of epidermal growth factor receptor (EGFR), which retained the transmembrane and kinase domains of the EGFR. This fusion gene product showed transforming potential and a high tumor-forming capacity in cell culture and in vivo. Mechanistically, HMGA2-EGFR constitutively induced a higher level of phosphorylated STAT5B than EGFRvIII, an in-frame exon deletion product of the EGFR gene that is commonly found in primary glioblastoma. Forced expression of HMGA2-EGFR enhanced orthotopic tumor formation of the U87MG human glioma cell line. Furthermore, the EGFR kinase inhibitor erlotinib blocked sphere formation of TGS-01 cells in culture and inhibited tumor formation in vivo. These findings suggest that, in addition to gene amplification and in-frame exon deletion, EGFR signaling can also be activated by gene fusion, suggesting a possible avenue for treatment of glioblastoma.

BACKGROUND In WHO grade II diffuse gliomas (low-grade gliomas, hereafter called LGGs), chemotherapy and radiotherapy contribute to prolonged survival but could induce somatic mutations. The optimal timing of treatment in LGGs remain poorly understood. To delineate this, we designed a mathematical model for tumor growth and investigate the association among the treatment, malignant transformation (MT), and the accumulation of somatic mutations revealed by whole exome sequencing (WES) in LGGs. METHODS Totally, 290 patients with LGGs between 1990 and 2018 were analyzed. We assessed the statuses of IDH mutation and 1p19q co-deletion in all tumors. Among all, 114 patients (39%) underwent MT during follow-up periods (mean: 82.6 months). Tumor volume was evaluated with FLAIR and/or T2-weighted MRI. MT was evaluated with contrast-enhanced MRI and/or pathological diagnosis. To investigate the number of somatic mutations in a cohort of LGGs and their patient matched recurrence, WES was performed on 88 serial samples collected at least two time-points from 39 patients. RESULTS Oligodendroglioma, IDH-mutant and 1p/19q-codeleted (OD) showed longer transformation-free survival compared to other subtypes. An exponential model was chosen to estimate growth rate in LGGs, since the exponential model provided a better fit to our data as compared to a linear model. The growth rate significantly decreased in the middle of chemotherapy and after radiotherapy. By contrast, these treatments increased the number of somatic mutations identified by WES and the rate of MT in each subtype. The increasing number of mutations in recurrent tumors showed strong correlation with the rise in MT rate. Based on the growth rate and the risk of MT, optimal timing of treatments could be calculated for each genetic subtype. CONCLUSIONS The mathematical model and WES analysis delineates the optimal timing of treatments in each subtype, which will help to decide the treatment for LGGs.

RET finger protein (RFP) plays a pivotal role in the acquisition of chemoresistance via formation of a complex with nuclear transcription factor Y and histone deacetylase 1, producing specific enhancer activity. We hypothesized that chemoresistance mediated by RFP may result from aberrant deacetylation of H3K27 and dysregulation of novel cis-regulatory (active) enhancers. Therefore, we investigated the effects of RFP on active enhancers and gene expression and evaluated the effects of RFP depletion on the growth of glioma cells treated with temozolomide both in vitro and in vivo. We found that the combination of RFP depletion and TMZ treatment markedly suppressed the growth of glioma cells and extended the survival time of intracranial tumor-bearing mice compared to that of TMZ alone. Chromatin immunoprecipitation and whole-transcriptome sequencing revealed that RFP depletion weakened a significant number of enhancers, and diminished RNA with functions related to mitosis, DNA-dependent DNA replication and cell cycle. Further, the transcriptomes of FOXO1 and TBP2 were significantly increased, while that of PARP-binding protein (PARPBP) was decreased, resulting in induction of reactive oxygen species and cell death. This study suggests that RFP contributes to chemoresistance via aberrant deacetylation of histone H3K27 and dysregulation of RFP-associated active-enhancers in glioma and that the combination of targeting RFP and TMZ has potential as an effective novel treatment strategy for lethal glioma.

Grade II and III gliomas are generally slowly progressing brain cancers, many of which eventually transform into more aggressive tumors. Despite recent findings of frequent mutations in IDH1 and other genes, knowledge about their pathogenesis is still incomplete. Here, combining two large sets of high-throughput sequencing data, we delineate the entire picture of genetic alterations and affected pathways in these glioma types, with sensitive detection of driver genes. Grade II and III gliomas comprise three distinct subtypes characterized by discrete sets of mutations and distinct clinical behaviors. Mutations showed significant positive and negative correlations and a chronological hierarchy, as inferred from different allelic burdens among coexisting mutations, suggesting that there is functional interplay between the mutations that drive clonal selection. Extensive serial and multi-regional sampling analyses further supported this finding and also identified a high degree of temporal and spatial heterogeneity generated during tumor expansion and relapse, which is likely shaped by the complex but ordered processes of multiple clonal selection and evolutionary events.

Isocitrate dehydrogenase 1 (IDH1), which localizes to the cytosol and peroxisomes, catalyzes the oxidative decarboxylation of isocitrate to α-ketoglutarate (α-KG) and in parallel converts NADP(+) to NADPH. IDH1 mutations are frequently detected in grades 2-4 gliomas and in acute myeloid leukemias (AML). Mutations of IDH1 have been identified at codon 132, with arginine being replaced with histidine in most cases. Mutant IDH1 gains novel enzyme activity converting α-KG to D-2-hydroxyglutarate (2-HG) which acts as a competitive inhibitor of α-KG. As a result, the activity of α-KG-dependent enzyme is reduced. Based on these findings, 2-HG has been proposed to be an oncometabolite. In this study, we established HEK293 and U87 cells that stably expressed IDH1-WT and IDH1-R132H and investigated the effect of glutaminase inhibition on cell proliferation with 6-diazo-5-oxo-L-norleucine (DON). We found that cell proliferation was suppressed in IDH1-R132H cells. The addition of α-KG restored cell proliferation. The metabolic features of 33 gliomas with wild type IDH1 (IDH1-WT) and with IDH1-R132H mutation were examined by global metabolome analysis using capillary electrophoresis time-of-flight mass spectrometry (CE-TOFMS). We showed that the 2-HG levels were highly elevated in gliomas with IDH1-R132H mutation. Intriguingly, in gliomas with IDH1-R132H, glutamine and glutamate levels were significantly reduced which implies replenishment of α-KG by glutaminolysis. Based on these results, we concluded that glutaminolysis is activated in gliomas with IDH1-R132H mutation and that development of novel therapeutic approaches targeting activated glutaminolysis is warranted.