Glioma progression is shaped by genetic evolution and microenvironment interactions.

The factors driving therapy resistance in diffuse glioma remain poorly understood. To identify treatment-associated cellular and genetic changes, we analyzed RNA and/or DNA sequencing data from the temporally separated tumor pairs of 304 adult patients with isocitrate dehydrogenase (IDH)-wild-type and IDH-mutant glioma. Tumors recurred in distinct manners that were dependent on IDH mutation status and attributable to changes in histological feature composition, somatic alterations, and microenvironment interactions. Hypermutation and acquired CDKN2A deletions were associated with an increase in proliferating neoplastic cells at recurrence in both glioma subtypes, reflecting active tumor growth. IDH-wild-type tumors were more invasive at recurrence, and their neoplastic cells exhibited increased expression of neuronal signaling programs that reflected a possible role for neuronal interactions in promoting glioma progression. Mesenchymal transition was associated with the presence of a myeloid cell state defined by specific ligand-receptor interactions with neoplastic cells. Collectively, these recurrence-associated phenotypes represent potential targets to alter disease progression.

Synthetic lethality-mediated precision oncology via the tumor transcriptome.

Precision oncology has made significant advances, mainly by targeting actionable mutations in cancer driver genes. Aiming to expand treatment opportunities, recent studies have begun to explore the utility of tumor transcriptome to guide patient treatment. Here, we introduce SELECT (synthetic lethality and rescue-mediated precision oncology via the transcriptome), a precision oncology framework harnessing genetic interactions to predict patient response to cancer therapy from the tumor transcriptome. SELECT is tested on a broad collection of 35 published targeted and immunotherapy clinical trials from 10 different cancer types. It is predictive of patients' response in 80% of these clinical trials and in the recent multi-arm WINTHER trial. The predictive signatures and the code are made publicly available for academic use, laying a basis for future prospective clinical studies.

A Hematogenous Route for Medulloblastoma Leptomeningeal Metastases.

While the preponderance of morbidity and mortality in medulloblastoma patients are due to metastatic disease, most research focuses on the primary tumor due to a dearth of metastatic tissue samples and model systems. Medulloblastoma metastases are found almost exclusively on the leptomeningeal surface of the brain and spinal cord; dissemination is therefore thought to occur through shedding of primary tumor cells into the cerebrospinal fluid followed by distal re-implantation on the leptomeninges. We present evidence for medulloblastoma circulating tumor cells (CTCs) in therapy-naive patients and demonstrate in vivo, through flank xenografting and parabiosis, that medulloblastoma CTCs can spread through the blood to the leptomeningeal space to form leptomeningeal metastases. Medulloblastoma leptomeningeal metastases express high levels of the chemokine CCL2, and expression of CCL2 in medulloblastoma in vivo is sufficient to drive leptomeningeal dissemination. Hematogenous dissemination of medulloblastoma offers a new opportunity to diagnose and treat lethal disseminated medulloblastoma.

Haploinsufficiency of phosphodiesterase 10A activates PI3K/AKT signaling independent of PTEN to induce an aggressive glioma phenotype.

Glioblastoma is universally fatal and characterized by frequent chromosomal copy number alterations harboring oncogenes and tumor suppressors. In this study, we analyzed exome-wide human glioblastoma copy number data and found that cytoband 6q27 is an independent poor prognostic marker in multiple data sets. We then combined CRISPR-Cas9 data, human spatial transcriptomic data, and human and mouse RNA sequencing data to nominate PDE10A as a potential haploinsufficient tumor suppressor in the 6q27 region. Mouse glioblastoma modeling using the RCAS/tv-a system confirmed that Pde10a suppression induced an aggressive glioma phenotype in vivo and resistance to temozolomide and radiation therapy in vitro. Cell culture analysis showed that decreased Pde10a expression led to increased PI3K/AKT signaling in a Pten-independent manner, a response blocked by selective PI3K inhibitors. Single-nucleus RNA sequencing from our mouse gliomas in vivo, in combination with cell culture validation, further showed that Pde10a suppression was associated with a proneural-to-mesenchymal transition that exhibited increased cell adhesion and decreased cell migration. Our results indicate that glioblastoma patients harboring PDE10A loss have worse outcomes and potentially increased sensitivity to PI3K inhibition.

Molecular Profiling Reveals Biologically Discrete Subsets and Pathways of Progression in Diffuse Glioma.

Therapy development for adult diffuse glioma is hindered by incomplete knowledge of somatic glioma driving alterations and suboptimal disease classification. We defined the complete set of genes associated with 1,122 diffuse grade II-III-IV gliomas from The Cancer Genome Atlas and used molecular profiles to improve disease classification, identify molecular correlations, and provide insights into the progression from low- to high-grade disease. Whole-genome sequencing data analysis determined that ATRX but not TERT promoter mutations are associated with increased telomere length. Recent advances in glioma classification based on IDH mutation and 1p/19q co-deletion status were recapitulated through analysis of DNA methylation profiles, which identified clinically relevant molecular subsets. A subtype of IDH mutant glioma was associated with DNA demethylation and poor outcome; a group of IDH-wild-type diffuse glioma showed molecular similarity to pilocytic astrocytoma and relatively favorable survival. Understanding of cohesive disease groups may aid improved clinical outcomes.

Identification of causal genetic drivers of human disease through systems-level analysis of regulatory networks.

Identification of driver mutations in human diseases is often limited by cohort size and availability of appropriate statistical models. We propose a framework for the systematic discovery of genetic alterations that are causal determinants of disease, by prioritizing genes upstream of functional disease drivers, within regulatory networks inferred de novo from experimental data. We tested this framework by identifying the genetic determinants of the mesenchymal subtype of glioblastoma. Our analysis uncovered KLHL9 deletions as upstream activators of two previously established master regulators of the subtype, C/EBPß and C/EBPδ. Rescue of KLHL9 expression induced proteasomal degradation of C/EBP proteins, abrogated the mesenchymal signature, and reduced tumor viability in vitro and in vivo. Deletions of KLHL9 were confirmed in > 50% of mesenchymal cases in an independent cohort, thus representing the most frequent genetic determinant of the subtype. The method generalized to study other human diseases, including breast cancer and Alzheimer's disease.

The somatic genomic landscape of glioblastoma.

We describe the landscape of somatic genomic alterations based on multidimensional and comprehensive characterization of more than 500 glioblastoma tumors (GBMs). We identify several novel mutated genes as well as complex rearrangements of signature receptors, including EGFR and PDGFRA. TERT promoter mutations are shown to correlate with elevated mRNA expression, supporting a role in telomerase reactivation. Correlative analyses confirm that the survival advantage of the proneural subtype is conferred by the G-CIMP phenotype, and MGMT DNA methylation may be a predictive biomarker for treatment response only in classical subtype GBM. Integrative analysis of genomic and proteomic profiles challenges the notion of therapeutic inhibition of a pathway as an alternative to inhibition of the target itself. These data will facilitate the discovery of therapeutic and diagnostic target candidates, the validation of research and clinical observations and the generation of unanticipated hypotheses that can advance our molecular understanding of this lethal cancer.

PKM2 phosphorylates histone H3 and promotes gene transcription and tumorigenesis.

Tumor-specific pyruvate kinase M2 (PKM2) is essential for the Warburg effect. In addition to its well-established role in aerobic glycolysis, PKM2 directly regulates gene transcription. However, the mechanism underlying this nonmetabolic function of PKM2 remains elusive. We show here that PKM2 directly binds to histone H3 and phosphorylates histone H3 at T11 upon EGF receptor activation. This phosphorylation is required for the dissociation of HDAC3 from the CCND1 and MYC promoter regions and subsequent acetylation of histone H3 at K9. PKM2-dependent histone H3 modifications are instrumental in EGF-induced expression of cyclin D1 and c-Myc, tumor cell proliferation, cell-cycle progression, and brain tumorigenesis. In addition, levels of histone H3 T11 phosphorylation correlate with nuclear PKM2 expression levels, glioma malignancy grades, and prognosis. These findings highlight the role of PKM2 as a protein kinase in its nonmetabolic functions of histone modification, which is essential for its epigenetic regulation of gene expression and tumorigenesis.

A clinically applicable integrative molecular classification of meningiomas.

Meningiomas are the most common primary intracranial tumour in adults1. Patients with symptoms are generally treated with surgery as there are no effective medical therapies. The World Health Organization histopathological grade of the tumour and the extent of resection at surgery (Simpson grade) are associated with the recurrence of disease; however, they do not accurately reflect the clinical behaviour of all meningiomas2. Molecular classifications of meningioma that reliably reflect tumour behaviour and inform on therapies are required. Here we introduce four consensus molecular groups of meningioma by combining DNA somatic copy-number aberrations, DNA somatic point mutations, DNA methylation and messenger RNA abundance in a unified analysis. These molecular groups more accurately predicted clinical outcomes compared with existing classification schemes. Each molecular group showed distinctive and prototypical biology (immunogenic, benign NF2 wild-type, hypermetabolic and proliferative) that informed therapeutic options. Proteogenomic characterization reinforced the robustness of the newly defined molecular groups and uncovered highly abundant and group-specific protein targets that we validated using immunohistochemistry. Single-cell RNA sequencing revealed inter-individual variations in meningioma as well as variations in intrinsic expression programs in neoplastic cells that mirrored the biology of the molecular groups identified.

Single-cell methylation sequencing data reveal succinct metastatic migration histories and tumor progression models.

Recent studies exploring the impact of methylation in tumor evolution suggest that although the methylation status of many of the CpG sites are preserved across distinct lineages, others are altered as the cancer progresses. Because changes in methylation status of a CpG site may be retained in mitosis, they could be used to infer the progression history of a tumor via single-cell lineage tree reconstruction. In this work, we introduce the first principled distance-based computational method, Sgootr, for inferring a tumor's single-cell methylation lineage tree and for jointly identifying lineage-informative CpG sites that harbor changes in methylation status that are retained along the lineage. We apply Sgootr on single-cell bisulfite-treated whole-genome sequencing data of multiregionally sampled tumor cells from nine metastatic colorectal cancer patients, as well as multiregionally sampled single-cell reduced-representation bisulfite sequencing data from a glioblastoma patient. We show that the tumor lineages constructed reveal a simple model underlying tumor progression and metastatic seeding. A comparison of Sgootr against alternative approaches shows that Sgootr can construct lineage trees with fewer migration events and with more in concordance with the sequential-progression model of tumor evolution, with a running time a fraction of that used in prior studies. Lineage-informative CpG sites identified by Sgootr are in inter-CpG island (CGI) regions, as opposed to intra-CGIs, which have been the main regions of interest in genomic methylation-related analyses.

Longitudinal molecular trajectories of diffuse glioma in adults.

The evolutionary processes that drive universal therapeutic resistance in adult patients with diffuse glioma remain unclear1,2. Here we analysed temporally separated DNA-sequencing data and matched clinical annotation from 222 adult patients with glioma. By analysing mutations and copy numbers across the three major subtypes of diffuse glioma, we found that driver genes detected at the initial stage of disease were retained at recurrence, whereas there was little evidence of recurrence-specific gene alterations. Treatment with alkylating agents resulted in a hypermutator phenotype at different rates across the glioma subtypes, and hypermutation was not associated with differences in overall survival. Acquired aneuploidy was frequently detected in recurrent gliomas and was characterized by IDH mutation but without co-deletion of chromosome arms 1p/19q, and further converged with acquired alterations in the cell cycle and poor outcomes. The clonal architecture of each tumour remained similar over time, but the presence of subclonal selection was associated with decreased survival. Finally, there were no differences in the levels of immunoediting between initial and recurrent gliomas. Collectively, our results suggest that the strongest selective pressures occur during early glioma development and that current therapies shape this evolution in a largely stochastic manner.

Metabologenomic characterization uncovers a clinically aggressive IDH mutant glioma subtype.

Mutations in the pivotal metabolic isocitrate dehydrogenase (IDH) enzymes are recognized to drive the molecular footprint of diffuse gliomas, and patients with IDH mutant gliomas have overall favorable outcomes compared to patients with IDH wild-type tumors. However, survival still varies widely among patients with IDH mutated tumors. Here, we aimed to characterize molecular signatures that explain the range of IDH mutant gliomas. By integrating matched epigenome-wide methylome, transcriptome, and global metabolome data in 154 patients with gliomas, we identified a group of IDH mutant gliomas with globally altered metabolism that resembled IDH wild-type tumors. IDH-mutant gliomas with altered metabolism have significantly shorter overall survival from their IDH mutant counterparts that is not fully accounted for by recognized molecular prognostic markers of CDKN2A/B loss and glioma CpG Island Methylator Phenotype (GCIMP) status. IDH-mutant tumors with dysregulated metabolism harbored distinct epigenetic alterations that converged to drive proliferative and stem-like transcriptional profiles, providing a window to target novel dependencies in gliomas.

Central nervous system embryonal tumors with EWSR1-PLAGL1 rearrangements reclassified as INI-1 deficient tumors at relapse.

PURPOSE: Central nervous system (CNS) embryonal tumors are a diverse group of malignant tumors typically affecting pediatric patients that recently have been better defined, and this paper describes evolution of a unique type of embryonal tumor at relapse. METHODS: Two pediatric patients with CNS embryonal tumors with EWSR1-PLAGL1 rearrangements treated at Arkansas Children's Hospital with histopathologic and molecular data are described. RESULTS: These two patients at diagnosis were classified as CNS embryonal tumors with EWSR1-PLAGL1 rearrangements based on histologic appearance and molecular data. At relapse both patient's disease was reclassified as atypical teratoid rhabdoid tumor (ATRT) based on loss of INI-1, presence of SMARCB1 alterations, and methylation profiling results. CONCLUSION: CNS embryonal tumors with EWSR1-PLAGL1 rearrangements acquire or include a population of cells with SMARCB1 alterations that are the component that predominate at relapse, suggesting treatment aimed at this disease component at diagnosis should be considered.

Virtual multi-institutional tumor board: a strategy for personalized diagnoses and management of rare CNS tumors.

PURPOSE: Multidisciplinary tumor boards (MTBs) integrate clinical, molecular, and radiological information and facilitate coordination of neuro-oncology care. During the COVID-19 pandemic, our MTB transitioned to a virtual and multi-institutional format. We hypothesized that this expansion would allow expert review of challenging neuro-oncology cases and contribute to the care of patients with limited access to specialized centers. METHODS: We retrospectively reviewed records from virtual MTBs held between 04/2020-03/2021. Data collected included measures of potential clinical impact, including referrals to observational or therapeutic studies, referrals for specialized neuropathology analysis, and whether molecular findings led to a change in diagnosis and/or guided management suggestions. RESULTS: During 25 meetings, 32 presenters discussed 44 cases. Approximately half (n = 20; 48%) involved a rare central nervous system (CNS) tumor. In 21% (n = 9) the diagnosis was changed or refined based on molecular profiling obtained at the NIH and in 36% (n = 15) molecular findings guided management. Clinical trial suggestions were offered to 31% (n = 13), enrollment in the observational NCI Natural History Study to 21% (n = 9), neuropathology review and molecular testing at the NIH to 17% (n = 7), and all received management suggestions. CONCLUSION: Virtual multi-institutional MTBs enable remote expert review of CNS tumors. We propose them as a strategy to facilitate expert opinions from specialized centers, especially for rare CNS tumors, helping mitigate geographic barriers to patient care and serving as a pre-screening tool for studies. Advanced molecular testing is key to obtaining a precise diagnosis, discovering potentially actionable targets, and guiding management.

Diffuse Pediatric-type High-grade Glioma Arising in an Ovarian Mature Cystic Teratoma.

Immature neuroectodermal tissue can be found in the ovary as part of an immature teratoma or as part of a teratoma with malignant neuroectodermal transformation. Such lesions may closely resemble central nervous system tumors, but their biologic similarity is unclear. We describe an 18-yr-old female who presented with abdominal pain caused by an ovarian mass with widespread metastases. Histology showed a primitive, high-grade tumor arising in the background of a mature teratoma. The tumor was SOX10 positive, with focal expression of GFAP, S100, NSE, and synaptophysin. Molecular analysis demonstrated co-amplification of PDGFRA and KIT , alterations common in high-grade gliomas. By whole-genome methylation profiling, it clustered into the "diffuse pediatric-type high-grade glioma, RTK1 subtype, subclass c" group. Despite progressing through 2 lines of chemotherapy with widespread metastatic disease, she achieved an excellent response to chemotherapy directed toward aggressive germ cell tumors. This case emphasizes the importance of immunohistochemical, genomic, and epigenetic analyses to accurately classify these exceedingly rare tumors and determine the optimal therapy.

Therapeutic targeting of ependymoma as informed by oncogenic enhancer profiling.

Genomic sequencing has driven precision-based oncology therapy; however, the genetic drivers of many malignancies remain unknown or non-targetable, so alternative approaches to the identification of therapeutic leads are necessary. Ependymomas are chemotherapy-resistant brain tumours, which, despite genomic sequencing, lack effective molecular targets. Intracranial ependymomas are segregated on the basis of anatomical location (supratentorial region or posterior fossa) and further divided into distinct molecular subgroups that reflect differences in the age of onset, gender predominance and response to therapy. The most common and aggressive subgroup, posterior fossa ependymoma group A (PF-EPN-A), occurs in young children and appears to lack recurrent somatic mutations. Conversely, posterior fossa ependymoma group B (PF-EPN-B) tumours display frequent large-scale copy number gains and losses but have favourable clinical outcomes. More than 70% of supratentorial ependymomas are defined by highly recurrent gene fusions in the NF-κB subunit gene RELA (ST-EPN-RELA), and a smaller number involve fusion of the gene encoding the transcriptional activator YAP1 (ST-EPN-YAP1). Subependymomas, a distinct histologic variant, can also be found within the supratetorial and posterior fossa compartments, and account for the majority of tumours in the molecular subgroups ST-EPN-SE and PF-EPN-SE. Here we describe mapping of active chromatin landscapes in 42 primary ependymomas in two non-overlapping primary ependymoma cohorts, with the goal of identifying essential super-enhancer-associated genes on which tumour cells depend. Enhancer regions revealed putative oncogenes, molecular targets and pathways; inhibition of these targets with small molecule inhibitors or short hairpin RNA diminished the proliferation of patient-derived neurospheres and increased survival in mouse models of ependymomas. Through profiling of transcriptional enhancers, our study provides a framework for target and drug discovery in other cancers that lack known genetic drivers and are therefore difficult to treat.

Mitochondria-Translocated PGK1 Functions as a Protein Kinase to Coordinate Glycolysis and the TCA Cycle in Tumorigenesis.

It is unclear how the Warburg effect that exemplifies enhanced glycolysis in the cytosol is coordinated with suppressed mitochondrial pyruvate metabolism. We demonstrate here that hypoxia, EGFR activation, and expression of K-Ras G12V and B-Raf V600E induce mitochondrial translocation of phosphoglycerate kinase 1 (PGK1); this is mediated by ERK-dependent PGK1 S203 phosphorylation and subsequent PIN1-mediated cis-trans isomerization. Mitochondrial PGK1 acts as a protein kinase to phosphorylate pyruvate dehydrogenase kinase 1 (PDHK1) at T338, which activates PDHK1 to phosphorylate and inhibit the pyruvate dehydrogenase (PDH) complex. This reduces mitochondrial pyruvate utilization, suppresses reactive oxygen species production, increases lactate production, and promotes brain tumorigenesis. Furthermore, PGK1 S203 and PDHK1 T338 phosphorylation levels correlate with PDH S293 inactivating phosphorylation levels and poor prognosis in glioblastoma patients. This work highlights that PGK1 acts as a protein kinase in coordinating glycolysis and the tricarboxylic acid (TCA) cycle, which is instrumental in cancer metabolism and tumorigenesis.