Glioblastoma evolution and heterogeneity from a 3D whole-tumor perspective.

Treatment failure for the lethal brain tumor glioblastoma (GBM) is attributed to intratumoral heterogeneity and tumor evolution. We utilized 3D neuronavigation during surgical resection to acquire samples representing the whole tumor mapped by 3D spatial coordinates. Integrative tissue and single-cell analysis revealed sources of genomic, epigenomic, and microenvironmental intratumoral heterogeneity and their spatial patterning. By distinguishing tumor-wide molecular features from those with regional specificity, we inferred GBM evolutionary trajectories from neurodevelopmental lineage origins and initiating events such as chromothripsis to emergence of genetic subclones and spatially restricted activation of differential tumor and microenvironmental programs in the core, periphery, and contrast-enhancing regions. Our work depicts GBM evolution and heterogeneity from a 3D whole-tumor perspective, highlights potential therapeutic targets that might circumvent heterogeneity-related failures, and establishes an interactive platform enabling 360° visualization and analysis of 3D spatial patterns for user-selected genes, programs, and other features across whole GBM tumors.

Integrative analysis of 111 reference human epigenomes.

The reference human genome sequence set the stage for studies of genetic variation and its association with human disease, but epigenomic studies lack a similar reference. To address this need, the NIH Roadmap Epigenomics Consortium generated the largest collection so far of human epigenomes for primary cells and tissues. Here we describe the integrative analysis of 111 reference human epigenomes generated as part of the programme, profiled for histone modification patterns, DNA accessibility, DNA methylation and RNA expression. We establish global maps of regulatory elements, define regulatory modules of coordinated activity, and their likely activators and repressors. We show that disease- and trait-associated genetic variants are enriched in tissue-specific epigenomic marks, revealing biologically relevant cell types for diverse human traits, and providing a resource for interpreting the molecular basis of human disease. Our results demonstrate the central role of epigenomic information for understanding gene regulation, cellular differentiation and human disease.

Clonal expansion and epigenetic reprogramming following deletion or amplification of mutant IDH1.

IDH1 mutation is the earliest genetic alteration in low-grade gliomas (LGGs), but its role in tumor recurrence is unclear. Mutant IDH1 drives overproduction of the oncometabolite d-2-hydroxyglutarate (2HG) and a CpG island (CGI) hypermethylation phenotype (G-CIMP). To investigate the role of mutant IDH1 at recurrence, we performed a longitudinal analysis of 50 IDH1 mutant LGGs. We discovered six cases with copy number alterations (CNAs) at the IDH1 locus at recurrence. Deletion or amplification of IDH1 was followed by clonal expansion and recurrence at a higher grade. Successful cultures derived from IDH1 mutant, but not IDH1 wild type, gliomas systematically deleted IDH1 in vitro and in vivo, further suggestive of selection against the heterozygous mutant state as tumors progress. Tumors and cultures with IDH1 CNA had decreased 2HG, maintenance of G-CIMP, and DNA methylation reprogramming outside CGI. Thus, while IDH1 mutation initiates gliomagenesis, in some patients mutant IDH1 and 2HG are not required for later clonal expansions.

Estimating absolute methylation levels at single-CpG resolution from methylation enrichment and restriction enzyme sequencing methods.

Recent advancements in sequencing-based DNA methylation profiling methods provide an unprecedented opportunity to map complete DNA methylomes. These include whole-genome bisulfite sequencing (WGBS, MethylC-seq, or BS-seq), reduced-representation bisulfite sequencing (RRBS), and enrichment-based methods such as MeDIP-seq, MBD-seq, and MRE-seq. These methods yield largely comparable results but differ significantly in extent of genomic CpG coverage, resolution, quantitative accuracy, and cost, at least while using current algorithms to interrogate the data. None of these existing methods provides single-CpG resolution, comprehensive genome-wide coverage, and cost feasibility for a typical laboratory. We introduce methylCRF, a novel conditional random fields-based algorithm that integrates methylated DNA immunoprecipitation (MeDIP-seq) and methylation-sensitive restriction enzyme (MRE-seq) sequencing data to predict DNA methylation levels at single-CpG resolution. Our method is a combined computational and experimental strategy to produce DNA methylomes of all 28 million CpGs in the human genome for a fraction (<10%) of the cost of whole-genome bisulfite sequencing methods. methylCRF was benchmarked for accuracy against Infinium arrays, RRBS, WGBS sequencing, and locus-specific bisulfite sequencing performed on the same human embryonic stem cell line. methylCRF transformation of MeDIP-seq/MRE-seq was equivalent to a biological replicate of WGBS in quantification, coverage, and resolution. We used conventional bisulfite conversion, PCR, cloning, and sequencing to validate loci where our predictions do not agree with whole-genome bisulfite data, and in 11 out of 12 cases, methylCRF predictions of methylation level agree better with validated results than does whole-genome bisulfite sequencing. Therefore, methylCRF transformation of MeDIP-seq/MRE-seq data provides an accurate, inexpensive, and widely accessible strategy to create full DNA methylomes.

Functional DNA methylation differences between tissues, cell types, and across individuals discovered using the M&M algorithm.

DNA methylation plays key roles in diverse biological processes such as X chromosome inactivation, transposable element repression, genomic imprinting, and tissue-specific gene expression. Sequencing-based DNA methylation profiling provides an unprecedented opportunity to map and compare complete DNA methylomes. This includes one of the most widely applied technologies for measuring DNA methylation: methylated DNA immunoprecipitation followed by sequencing (MeDIP-seq), coupled with a complementary method, methylation-sensitive restriction enzyme sequencing (MRE-seq). A computational approach that integrates data from these two different but complementary assays and predicts methylation differences between samples has been unavailable. Here, we present a novel integrative statistical framework M&M (for integration of MeDIP-seq and MRE-seq) that dynamically scales, normalizes, and combines MeDIP-seq and MRE-seq data to detect differentially methylated regions. Using sample-matched whole-genome bisulfite sequencing (WGBS) as a gold standard, we demonstrate superior accuracy and reproducibility of M&M compared to existing analytical methods for MeDIP-seq data alone. M&M leverages the complementary nature of MeDIP-seq and MRE-seq data to allow rapid comparative analysis between whole methylomes at a fraction of the cost of WGBS. Comprehensive analysis of nineteen human DNA methylomes with M&M reveals distinct DNA methylation patterns among different tissue types, cell types, and individuals, potentially underscoring divergent epigenetic regulation at different scales of phenotypic diversity. We find that differential DNA methylation at enhancer elements, with concurrent changes in histone modifications and transcription factor binding, is common at the cell, tissue, and individual levels, whereas promoter methylation is more prominent in reinforcing fundamental tissue identities.

Conserved role of intragenic DNA methylation in regulating alternative promoters.

Although it is known that the methylation of DNA in 5' promoters suppresses gene expression, the role of DNA methylation in gene bodies is unclear. In mammals, tissue- and cell type-specific methylation is present in a small percentage of 5' CpG island (CGI) promoters, whereas a far greater proportion occurs across gene bodies, coinciding with highly conserved sequences. Tissue-specific intragenic methylation might reduce, or, paradoxically, enhance transcription elongation efficiency. Capped analysis of gene expression (CAGE) experiments also indicate that transcription commonly initiates within and between genes. To investigate the role of intragenic methylation, we generated a map of DNA methylation from the human brain encompassing 24.7 million of the 28 million CpG sites. From the dense, high-resolution coverage of CpG islands, the majority of methylated CpG islands were shown to be in intragenic and intergenic regions, whereas less than 3% of CpG islands in 5' promoters were methylated. The CpG islands in all three locations overlapped with RNA markers of transcription initiation, and unmethylated CpG islands also overlapped significantly with trimethylation of H3K4, a histone modification enriched at promoters. The general and CpG-island-specific patterns of methylation are conserved in mouse tissues. An in-depth investigation of the human SHANK3 locus and its mouse homologue demonstrated that this tissue-specific DNA methylation regulates intragenic promoter activity in vitro and in vivo. These methylation-regulated, alternative transcripts are expressed in a tissue- and cell type-specific manner, and are expressed differentially within a single cell type from distinct brain regions. These results support a major role for intragenic methylation in regulating cell context-specific alternative promoters in gene bodies.

Evolution of DNA repair defects during malignant progression of low-grade gliomas after temozolomide treatment.

Temozolomide (TMZ) increases the overall survival of patients with glioblastoma (GBM), but its role in the clinical management of diffuse low-grade gliomas (LGG) is still being defined. DNA hypermethylation of the O (6) -methylguanine-DNA methyltransferase (MGMT) promoter is associated with an improved response to TMZ treatment, while inactivation of the DNA mismatch repair (MMR) pathway is associated with therapeutic resistance and TMZ-induced mutagenesis. We previously demonstrated that TMZ treatment of LGG induces driver mutations in the RB and AKT-mTOR pathways, which may drive malignant progression to secondary GBM. To better understand the mechanisms underlying TMZ-induced mutagenesis and malignant progression, we explored the evolution of MGMT methylation and genetic alterations affecting MMR genes in a cohort of 34 treatment-naïve LGGs and their recurrences. Recurrences with TMZ-associated hypermutation had increased MGMT methylation compared to their untreated initial tumors and higher overall MGMT methylation compared to TMZ-treated non-hypermutated recurrences. A TMZ-associated mutation in one or more MMR genes was observed in five out of six TMZ-treated hypermutated recurrences. In two cases, pre-existing heterozygous deletions encompassing MGMT, or an MMR gene, were followed by TMZ-associated mutations in one of the genes of interest. These results suggest that tumor cells with methylated MGMT may undergo positive selection during TMZ treatment in the context of MMR deficiency.

Epigenome analyses using BAC microarrays identify evolutionary conservation of tissue-specific methylation of SHANK3.

CpG islands are present in one-half of all human and mouse genes and typically overlap with promoters or exons. We developed a method for high-resolution analysis of the methylation status of CpG islands genome-wide, using arrays of BAC clones and the methylation-sensitive restriction enzyme NotI. Here we demonstrate the accuracy and specificity of the method. By computationally mapping all NotI sites, methylation events can be defined with single-nucleotide precision throughout the genome. We also demonstrate the unique expandability of the array method using a different methylation-sensitive restriction enzyme, BssHII. We identified and validated new CpG island loci that are methylated in a tissue-specific manner in normal human tissues. The methylation status of the CpG islands is associated with gene expression for several genes, including SHANK3, which encodes a structural protein in neuronal postsynaptic densities. Defects in SHANK3 seem to underlie human 22q13 deletion syndrome. Furthermore, these patterns for SHANK3 are conserved in mice and rats.

Whole tumor analysis reveals early origin of the TERT promoter mutation and intercellular heterogeneity in TERT expression.

BACKGROUND: The TERT promoter mutation (TPM) is acquired in most IDH-wildtype glioblastomas (GBM) and IDH-mutant oligodendrogliomas (OD) enabling tumor cell immortality. Previous studies on TPM clonality show conflicting results. This study was performed to determine whether TPM is clonal on a tumor-wide scale. METHODS: We investigated TPM clonality in relation to presumed early events in 19 IDH-wildtype GBM and 10 IDH-mutant OD using 3-dimensional comprehensive tumor sampling. We performed Sanger sequencing on 264 tumor samples and deep amplicon sequencing on 187 tumor samples. We obtained tumor purity and copy number estimates from whole exome sequencing. TERT expression was assessed by RNA-seq and RNAscope. RESULTS: We detected TPM in 100% of tumor samples with quantifiable tumor purity (219 samples). Variant allele frequencies (VAF) of TPM correlate positively with chromosome 10 loss in GBM (R = 0.85), IDH1 mutation in OD (R = 0.87), and with tumor purity (R = 0.91 for GBM; R = 0.90 for OD). In comparison, oncogene amplification was tumor-wide for MDM4- and most EGFR-amplified cases but heterogeneous for MYCN and PDGFRA, and strikingly high in low-purity samples. TPM VAF was moderately correlated with TERT expression (R = 0.52 for GBM; R = 0.65 for OD). TERT expression was detected in a subset of cells, solely in TPM-positive samples, including samples equivocal for tumor. CONCLUSIONS: On a tumor-wide scale, TPM is among the earliest events in glioma evolution. Intercellular heterogeneity of TERT expression, however, suggests dynamic regulation during tumor growth. TERT expression may be a tumor cell-specific biomarker.

Challenges in the discovery of tumor-specific alternative splicing-derived cell-surface antigens in glioma.

Despite advancements in cancer immunotherapy, solid tumors remain formidable challenges. In glioma, profound inter- and intra-tumoral heterogeneity of antigen landscape hampers therapeutic development. Therefore, it is critical to consider alternative sources to expand the repertoire of targetable (neo-)antigens and improve therapeutic outcomes. Accumulating evidence suggests that tumor-specific alternative splicing (AS) could be an untapped reservoir of antigens. In this study, we investigated tumor-specific AS events in glioma, focusing on those predicted to generate major histocompatibility complex (MHC)-presentation-independent, cell-surface antigens that could be targeted by antibodies and chimeric antigen receptor-T cells. We systematically analyzed bulk RNA-sequencing datasets comparing 429 tumor samples (from The Cancer Genome Atlas) and 9166 normal tissue samples (from the Genotype-Tissue Expression project), and identified 13 AS events in 7 genes predicted to be expressed in more than 10% of the patients, including PTPRZ1 and BCAN, which were corroborated by an external RNA-sequencing dataset. Subsequently, we validated our predictions and elucidated the complexity of the isoforms using full-length transcript amplicon sequencing on patient-derived glioblastoma cells. However, analyses of the RNA-sequencing datasets of spatially mapped and longitudinally collected clinical tumor samples unveiled remarkable spatiotemporal heterogeneity of the candidate AS events. Furthermore, proteomics analysis did not reveal any peptide spectra matching the putative antigens. Our investigation illustrated the diverse characteristics of the tumor-specific AS events and the challenges of antigen exploration due to their notable spatiotemporal heterogeneity and elusive nature at the protein levels. Redirecting future efforts toward intracellular, MHC-presented antigens could offer a more viable avenue.

Integrated genomic and epigenomic analyses pinpoint biallelic gene inactivation in tumors.

Aberrant methylation of CpG islands and genomic deletion are two predominant mechanisms of gene inactivation in tumorigenesis, but the extent to which they interact is largely unknown. The lack of an integrated approach to study these mechanisms has limited the understanding of tumor genomes and cancer genes. Restriction landmark genomic scanning (RLGS; ref. 1) is useful for global analysis of aberrant methylation of CpG islands, but has not been amenable to alignment with deletion maps because the identity of most RLGS fragments is unknown. Here, we determined the nucleotide sequence and exact chromosomal position of RLGS fragments throughout the genome using the whole chromosome of origin of the fragments and in silico restriction digestion of the human genome sequence. To study the interaction of these gene-inactivation mechanisms in primary brain tumors, we integrated RLGS-based methylation analysis with high-resolution deletion maps from microarray-based comparative genomic hybridization (array CGH; ref. 3). Certain subsets of gene-associated CpG islands were preferentially affected by convergent methylation and deletion, including genes that exhibit tumor-suppressor activity, such as CISH1 (encoding SOCS1; ref. 4), as well as genes such as COE3 that have been missed by traditional non-integrated approaches. Our results show that most aberrant methylation events are focal and independent of deletions, and the rare convergence of these mechanisms can pinpoint biallelic gene inactivation without the use of positional cloning.

Challenges in the discovery of tumor-specific alternative splicing-derived cell-surface antigens in glioma.

Background: Despite advancements in cancer immunotherapy, solid tumors remain formidable challenges. In glioma, profound inter-and intra-tumoral heterogeneity of antigen landscape hampers therapeutic development. Therefore, it is critical to consider alternative sources to expand the repertoire of targetable (neo-)antigens and improve therapeutic outcomes. Accumulating evidence suggests that tumor-specific alternative splicing (AS) could be an untapped reservoir of neoantigens. Results: In this study, we investigated tumor-specific AS events in glioma, focusing on those predicted to generate major histocompatibility complex (MHC)-presentation-independent, cell-surface neoantigens that could be targeted by antibodies and chimeric antigen receptor (CAR)-T cells. We systematically analyzed bulk RNA-sequencing datasets comparing 429 tumor samples (from The Cancer Genome Atlas [TCGA]) and 9,166 normal tissue samples (from the Genotype-Tissue Expression project [GTEx]), and identified 13 AS events in 7 genes predicted to be expressed in more than 10% of the patients, including PTPRZ1 and BCAN , which were corroborated by an external RNA-sequencing dataset. Subsequently, we validated our predictions and elucidated the complexity of the isoforms using full-length transcript amplicon sequencing on patient-derived glioblastoma cells. However, analyses of the RNA-sequencing datasets of spatially mapped and longitudinally collected clinical tumor samples unveiled remarkable spatiotemporal heterogeneity of the candidate AS events. Furthermore, proteomics analysis did not reveal any peptide spectra matching the putative neoantigens. Conclusions: Our investigation illustrated the diverse characteristics of the tumor-specific AS events and the challenges of antigen exploration due to their notable spatiotemporal heterogeneity and elusive nature at the protein levels. Redirecting future efforts toward intracellular, MHC-presented antigens could offer a more viable avenue.

Tumor-wide RNA splicing aberrations generate immunogenic public neoantigens.

T-cell-mediated immunotherapies are limited by the extent to which cancer-specific antigens are homogenously expressed throughout a tumor. We reasoned that recurrent splicing aberrations in cancer represent a potential source of tumor-wide and public neoantigens, and to test this possibility, we developed a novel pipeline for identifying neojunctions expressed uniformly within a tumor across diverse cancer types. Our analyses revealed multiple neojunctions that recur across patients and either exhibited intratumor heterogeneity or, in some cases, were tumor-wide. We identified CD8+ T-cell clones specific for neoantigens derived from tumor-wide and conserved neojunctions in GNAS and RPL22 , respectively. TCR-engineered CD8 + T-cells targeting these mutations conferred neoantigen-specific tumor cell eradication. Furthermore, we revealed that cancer-specific dysregulation in splicing factor expression leads to recurrent neojunction expression. Together, these data reveal that a subset of neojunctions are both intratumorally conserved and public, providing the molecular basis for novel T-cell-based immunotherapies that address intratumoral heterogeneity.

Novel SOX10 indel mutations drive schwannomas through impaired transactivation of myelination gene programs.

BACKGROUND: Schwannomas are common peripheral nerve sheath tumors that can cause severe morbidity given their stereotypic intracranial and paraspinal locations. Similar to many solid tumors, schwannomas and other nerve sheath tumors are primarily thought to arise due to aberrant hyperactivation of the RAS growth factor signaling pathway. Here, we sought to further define the molecular pathogenesis of schwannomas. METHODS: We performed comprehensive genomic profiling on a cohort of 96 human schwannomas, as well as DNA methylation profiling on a subset. Functional studies including RNA sequencing, chromatin immunoprecipitation-DNA sequencing, electrophoretic mobility shift assay, and luciferase reporter assays were performed in a fetal glial cell model following transduction with wildtype and tumor-derived mutant isoforms of SOX10. RESULTS: We identified that nearly one-third of sporadic schwannomas lack alterations in known nerve sheath tumor genes and instead harbor novel recurrent in-frame insertion/deletion mutations in SOX10, which encodes a transcription factor responsible for controlling Schwann cell differentiation and myelination. SOX10 indel mutations were highly enriched in schwannomas arising from nonvestibular cranial nerves (eg facial, trigeminal, vagus) and were absent from vestibular nerve schwannomas driven by NF2 mutation. Functional studies revealed these SOX10 indel mutations have retained DNA binding capacity but impaired transactivation of glial differentiation and myelination gene programs. CONCLUSIONS: We thus speculate that SOX10 indel mutations drive a unique subtype of schwannomas by impeding proper differentiation of immature Schwann cells.

Imaging biomarkers of TERT or GABPB1 silencing in TERT-positive glioblastoma.

BACKGROUND: TERT promoter mutations are observed in 80% of wild-type IDH glioblastoma (GBM). Moreover, the upstream TERT transcription factor GABPB1 was recently identified as a cancer-specific therapeutic target for tumors harboring a TERT promoter mutation. In that context, noninvasive imaging biomarkers are needed for the detection of TERT modulation. METHODS: Multiple GBM models were investigated as cells and in vivo tumors and the impact of TERT silencing, either directly or by targeting GABPB1, was determined using 1H and hyperpolarized 13C magnetic resonance spectroscopy (MRS). Changes in associated metabolic enzymes were also investigated. RESULTS: 1H-MRS revealed that lactate and glutathione (GSH) were the most significantly altered metabolites when either TERT or GABPB1 was silenced, and lactate and GSH levels were correlated with cellular TERT expression. Consistent with the drop in lactate, 13C-MRS showed that hyperpolarized [1-13C]lactate production from [1-13C]pyruvate was also reduced when TERT was silenced. Mechanistically, the reduction in GSH was associated with a reduction in pentose phosphate pathway flux, reduced activity of glucose-6-phosphate dehydrogenase, and reduced NADPH. The drop in lactate and hyperpolarized lactate were associated with reductions in glycolytic flux, NADH, and expression/activity of GLUT1, monocarboxylate transporters, and lactate dehydrogenase A. CONCLUSIONS: Our study indicates that MRS-detectable GSH, lactate, and lactate production could serve as metabolic biomarkers of response to emerging TERT-targeted therapies for GBM with activating TERT promoter mutations. Importantly these biomarkers are readily translatable to the clinic, and thus could ultimately improve GBM patient management.

Conserved features of TERT promoter duplications reveal an activation mechanism that mimics hotspot mutations in cancer.

Mutations in the TERT promoter represent the genetic underpinnings of tumor cell immortality. Beyond the two most common point mutations, which selectively recruit the ETS factor GABP to activate TERT, the significance of other variants is unknown. In seven cancer types, we identify duplications of wildtype sequence within the core promoter region of TERT that have strikingly similar features including an ETS motif, the duplication length and insertion site. The duplications recruit a GABP tetramer by virtue of the native ETS motif and its precisely spaced duplicated counterpart, activate the promoter and are clonal in a TERT expressing multifocal glioblastoma. We conclude that recurrent TERT promoter duplications are functionally and mechanistically equivalent to the hotspot mutations that confer tumor cell immortality. The shared mechanism of these divergent somatic genetic alterations suggests a strong selective pressure for recruitment of the GABP tetramer to activate TERT.

Human placental cytotrophoblast epigenome dynamics over gestation and alterations in placental disease.

The human placenta and its specialized cytotrophoblasts rapidly develop, have a compressed lifespan, govern pregnancy outcomes, and program the offspring's health. Understanding the molecular underpinnings of these behaviors informs development and disease. Profiling the extraembryonic epigenome and transcriptome during the 2nd and 3rd trimesters revealed H3K9 trimethylation overlapping deeply DNA hypomethylated domains with reduced gene expression and compartment-specific patterns that illuminated their functions. Cytotrophoblast DNA methylation increased, and several key histone modifications decreased across the genome as pregnancy advanced. Cytotrophoblasts from severe preeclampsia had substantially increased H3K27 acetylation globally and at genes that are normally downregulated at term but upregulated in this syndrome. In addition, some cases had an immature pattern of H3K27ac peaks, and others showed evidence of accelerated aging, suggesting subtype-specific alterations in severe preeclampsia. Thus, the cytotrophoblast epigenome dramatically reprograms during pregnancy, placental disease is associated with failures in this process, and H3K27 hyperacetylation is a feature of severe preeclampsia.

Temozolomide-induced hypermutation is associated with distant recurrence and reduced survival after high-grade transformation of low-grade IDH-mutant gliomas.

BACKGROUND: Chemotherapy improves overall survival after surgery and radiotherapy for newly diagnosed high-risk IDH-mutant low-grade gliomas (LGGs), but a proportion of patients treated with temozolomide (TMZ) will develop recurrent tumors with TMZ-induced hypermutation. We aimed to determine the prevalence of TMZ-induced hypermutation at recurrence and prognostic implications. METHODS: We sequenced recurrent tumors from 82 patients with initially low-grade IDH-mutant gliomas who underwent reoperation and correlated hypermutation status with grade at recurrence and subsequent clinical outcomes. RESULTS: Hypermutation was associated with high-grade disease at the time of reoperation (OR 12.0 95% CI 2.5-115.5, P = .002) and was identified at transformation in 57% of recurrent LGGs previously exposed to TMZ. After anaplastic (grade III) transformation, hypermutation was associated with shorter survival on univariate and multivariate analysis (HR 3.4, 95% CI 1.2-9.9, P = .024), controlling for tumor grade, subtype, age, and prior radiotherapy. The effect of hypermutation on survival after transformation was validated in an independent, published dataset. Hypermutated (HM) tumors were more likely to develop discontiguous foci of disease in the brain and spine (P = .003). To estimate the overall incidence of high-grade transformation among low-grade IDH-mutant tumors, data from a phase II trial of TMZ for LGG were analyzed. Eight-year transformation-free survival was 53.8% (95% CI 42.8-69.2), and 61% of analyzed transformed cases were HM. CONCLUSIONS: TMZ-induced hypermutation is a common event in transformed LGG previously treated with TMZ and is associated with worse prognosis and development of discontiguous disease after recurrence. These findings impact tumor classification at recurrence, prognostication, and clinical trial design.

Epigenome scans and cancer genome sequencing converge on WNK2, a kinase-independent suppressor of cell growth.

Human cancer genome and epigenome projects aim to identify new cancer genes and targets for therapy that have been overlooked by conventional approaches. Here we integrated large-scale genomics and epigenomics of 31 human infiltrative gliomas and identified low-frequency deletion and highly recurrent epigenetic silencing of WNK2, encoding a putative serine/threonine kinase. Prior cancer genome sequencing projects also identified point mutations in WNK1-4, suggesting that WNK family genes may have a role in cancers. We observed consistent gene silencing in tumors with dense aberrant methylation across 1.3 kb of the CpG island but more variable expression when the 5'-most region remained unmethylated. This primary tumor data fit well with WNK2 promoter analysis, which showed strong promoter activity in the 5'-most region, equivalent to the simian virus 40 promoter, but no activity in the 3' region. WT WNK2 exhibited autophosphorylation and protein kinase activity that was enhanced in cells exposed to hypertonic conditions, similar to WNK1. WNK2 inhibited up to 78% of colony formation by glioma cells but in an unexpectedly kinase-independent manner. The WNK2 silencing by epigenetic mechanisms was significantly associated (P < 0.01) with a known genetic signature of chemosensitive oligodendroglial tumors, 1p and 19q deletion, in two small but independent tumor sets. Taken together, the epigenetic silencing, occasional deletion and point mutation, and functional assessment suggest that aberrations of WNK2 may contribute to unregulated tumor cell growth. Thus, our integrated genetic and epigenetic approach might be useful to identify genes that are widely relevant to cancer, even when genetic alterations of the locus are infrequent.

MGMT promoter methylation level in newly diagnosed low-grade glioma is a predictor of hypermutation at recurrence.

BACKGROUND: Emerging data suggest that a subset of patients with diffuse isocitrate dehydrogenase (IDH)-mutant low-grade glioma (LGG) who receive adjuvant temozolomide (TMZ) recur with hypermutation in association with malignant progression to higher-grade tumors. It is currently unclear why some TMZ-treated LGG patients recur with hypermutation while others do not. MGMT encodes O6-methylguanine-DNA methyltransferase, a DNA repair protein that removes cytotoxic and potentially mutagenic lesions induced by TMZ. Here, we hypothesize that epigenetic silencing of MGMT by promoter methylation facilitates TMZ-induced mutagenesis in LGG patients and contributes to development of hypermutation at recurrence. METHODS: We utilize a quantitative deep sequencing assay to characterize MGMT promoter methylation in 109 surgical tissue specimens from initial tumors and post-treatment recurrences of 37 TMZ-treated LGG patients. We utilize methylation arrays to validate our sequencing assay, RNA sequencing to assess the relationship between methylation and gene expression, and exome sequencing to determine hypermutation status. RESULTS: Methylation level at the MGMT promoter is significantly higher in initial tumors of patients that develop hypermutation at recurrence relative to initial tumors of patients that do not (45.7% vs 34.8%, P = 0.027). Methylation level in initial tumors can predict hypermutation at recurrence in univariate models and multivariate models that incorporate patient age and molecular subtype. CONCLUSIONS: These findings reveal a mechanistic basis for observed differences in patient susceptibility to TMZ-driven hypermutation. Furthermore, they establish MGMT promoter methylation level as a potential biomarker to inform clinical management of LGG patients, including monitoring and treatment decisions, by predicting risk of hypermutation at recurrence.