IDHwt glioblastomas can be stratified by their transcriptional response to standard treatment, with implications for targeted therapy.

BACKGROUND: Glioblastoma (GBM) brain tumors lacking IDH1 mutations (IDHwt) have the worst prognosis of all brain neoplasms. Patients receive surgery and chemoradiotherapy but tumors almost always fatally recur. RESULTS: Using RNA sequencing data from 107 pairs of pre- and post-standard treatment locally recurrent IDHwt GBM tumors, we identify two responder subtypes based on longitudinal changes in gene expression. In two thirds of patients, a specific subset of genes is upregulated from primary to recurrence (Up responders), and in one third, the same genes are downregulated (Down responders), specifically in neoplastic cells. Characterization of the responder subtypes indicates subtype-specific adaptive treatment resistance mechanisms that are associated with distinct changes in the tumor microenvironment. In Up responders, recurrent tumors are enriched in quiescent proneural GBM stem cells and differentiated neoplastic cells, with increased interaction with the surrounding normal brain and neurotransmitter signaling, whereas Down responders commonly undergo mesenchymal transition. ChIP-sequencing data from longitudinal GBM tumors suggests that the observed transcriptional reprogramming could be driven by Polycomb-based chromatin remodeling rather than DNA methylation. CONCLUSIONS: We show that the responder subtype is cancer-cell intrinsic, recapitulated in in vitro GBM cell models, and influenced by the presence of the tumor microenvironment. Stratifying GBM tumors by responder subtype may lead to more effective treatment.

Defining the Role of Extrachromosomal DNA Amplifications in Medulloblastoma.

Circular extrachromosomal DNA (ecDNA), a common mechanism of oncogene amplification, has been identified as a major contributor to intratumoral heterogeneity and patient outcomes. In a recent publication in Nature Genetics, Chapman and colleagues further explored the role of ecDNA in the context of medulloblastoma. Using whole-genome sequencing, they found that 18% of the patients carry ecDNA amplification across a 468 medulloblastoma patient cohort. The presence of ecDNA was associated with worse survival. Single-cell FISH imaging and multiomic sequencing revealed that ecDNA copy number displayed a cell-to-cell variability within the sample, contributing to tumor heterogeneity. Furthermore, through sequencing and CRISPRi experiments, the authors uncovered frequent enhancer rewiring events on ecDNA that drive proliferation.

The Epigenetic Evolution of Glioma Is Determined by the IDH1 Mutation Status and Treatment Regimen.

Tumor adaptation or selection is thought to underlie therapy resistance in glioma. To investigate longitudinal epigenetic evolution of gliomas in response to therapeutic pressure, we performed an epigenomic analysis of 132 matched initial and recurrent tumors from patients with IDH-wildtype (IDHwt) and IDH-mutant (IDHmut) glioma. IDHwt gliomas showed a stable epigenome over time with relatively low levels of global methylation. The epigenome of IDHmut gliomas showed initial high levels of genome-wide DNA methylation that was progressively reduced to levels similar to those of IDHwt tumors. Integration of epigenomics, gene expression, and functional genomics identified HOXD13 as a master regulator of IDHmut astrocytoma evolution. Furthermore, relapse of IDHmut tumors was accompanied by histologic progression that was associated with survival, as validated in an independent cohort. Finally, the initial cell composition of the tumor microenvironment varied between IDHwt and IDHmut tumors and changed differentially following treatment, suggesting increased neoangiogenesis and T-cell infiltration upon treatment of IDHmut gliomas. This study provides one of the largest cohorts of paired longitudinal glioma samples with epigenomic, transcriptomic, and genomic profiling and suggests that treatment of IDHmut glioma is associated with epigenomic evolution toward an IDHwt-like phenotype. SIGNIFICANCE: Standard treatments are related to loss of DNA methylation in IDHmut glioma, resulting in epigenetic activation of genes associated with tumor progression and alterations in the microenvironment that resemble treatment-naïve IDHwt glioma.

Longitudinal characteristics of T2-FLAIR mismatch in IDH-mutant astrocytomas: Relation to grade, histopathology, and overall survival in the GLASS-NL cohort.

Background: The T2-FLAIR mismatch sign is defined by signal loss of the T2-weighted hyperintense area with Fluid-Attenuated Inversion Recovery (FLAIR) on magnetic resonance imaging, causing a hypointense region on FLAIR. It is a highly specific diagnostic marker for IDH-mutant astrocytoma and is postulated to be caused by intercellular microcystic change in the tumor tissue. However, not all IDH-mutant astrocytomas show this mismatch sign and some show the phenomenon in only part of the lesion. The aim of the study is to determine whether the T2-FLAIR mismatch phenomenon has any prognostic value beyond initial noninvasive molecular diagnosis. Methods: Patients initially diagnosed with histologically lower-grade (2 or 3) IDH-mutant astrocytoma and with at least 2 surgical resections were included in the GLASS-NL cohort. T2-FLAIR mismatch was determined, and the growth pattern of the recurrent tumor immediately before the second resection was annotated as invasive or expansive. The relation between the T2-FLAIR mismatch sign and tumor grade, microcystic change, overall survival (OS), and other clinical parameters was investigated both at first and second resection. Results: The T2-FLAIR mismatch sign was significantly related to Grade 2 (80% vs 51%), longer post-resection median OS (8.3 vs 5.2 years), expansive growth, and lower age at second resection. At first resection, no relation was found between the mismatch sign and OS. Microcystic change was associated with areas of T2-FLAIR mismatch. Conclusions: T2-FLAIR mismatch in IDH-mutant astrocytomas is correlated with microcystic change in the tumor tissue, favorable prognosis, and Grade 2 tumors at the time of second resection.

Quality of life after surgery for lower grade gliomas.

BACKGROUND: Few large studies have investigated quality of life (QOL) for adults diagnosed with lower grade glioma (LGG). METHODS: QOL was assessed for 320 adults with LGG (World Health Organization grade 2/3) enrolled in the International Low Grade Glioma Registry by using the Medical Outcomes Study 36-Item Short Form health survey. Data on symptoms were also collected. QOL outcomes were examined by treatment group and also compared to those from a population-based case-control study of meningioma (the Meningioma Consortium), in which 1722 meningioma cases diagnosed among residents of Connecticut, Massachusetts, California, Texas, and North Carolina from May 1, 2006 through March 14, 2013 were enrolled and frequency matched to 1622 controls by age, sex, and geography. RESULTS: The LGG sample average age is 45 years at the time of interview and 53.1% male. Almost 55% of patients had received radiation and chemotherapy (primarily temozolomide); 32.4% had received neither treatment. Two-thirds of participants with LGG report difficulty with speaking, memory, or thinking, and over one of three reports personality change or difficulty driving. After controlling for age and other comorbidities, individuals with LGG report levels of physical, emotional, and mental health functioning below those reported in a meningioma as well as a general healthy population. CONCLUSIONS: Despite being relatively young, persons with LGG report significantly reduced QOL compared to persons with nonmalignant brain tumors and to a control population, which highlights the need to better acknowledge and manage these symptoms for this group of patients diagnosed in the prime of life.

Extrachromosomal DNA in the cancerous transformation of Barrett's oesophagus.

Oncogene amplification on extrachromosomal DNA (ecDNA) drives the evolution of tumours and their resistance to treatment, and is associated with poor outcomes for patients with cancer1-6. At present, it is unclear whether ecDNA is a later manifestation of genomic instability, or whether it can be an early event in the transition from dysplasia to cancer. Here, to better understand the development of ecDNA, we analysed whole-genome sequencing (WGS) data from patients with oesophageal adenocarcinoma (EAC) or Barrett's oesophagus. These data included 206 biopsies in Barrett's oesophagus surveillance and EAC cohorts from Cambridge University. We also analysed WGS and histology data from biopsies that were collected across multiple regions at 2 time points from 80 patients in a case-control study at the Fred Hutchinson Cancer Center. In the Cambridge cohorts, the frequency of ecDNA increased between Barrett's-oesophagus-associated early-stage (24%) and late-stage (43%) EAC, suggesting that ecDNA is formed during cancer progression. In the cohort from the Fred Hutchinson Cancer Center, 33% of patients who developed EAC had at least one oesophageal biopsy with ecDNA before or at the diagnosis of EAC. In biopsies that were collected before cancer diagnosis, higher levels of ecDNA were present in samples from patients who later developed EAC than in samples from those who did not. We found that ecDNAs contained diverse collections of oncogenes and immunomodulatory genes. Furthermore, ecDNAs showed increases in copy number and structural complexity at more advanced stages of disease. Our findings show that ecDNA can develop early in the transition from high-grade dysplasia to cancer, and that ecDNAs progressively form and evolve under positive selection.

Driver Mutations Dictate the Immunologic Landscape and Response to Checkpoint Immunotherapy of Glioblastoma.

The composition of the tumor immune microenvironment (TIME) is considered a key determinant of patients' response to immunotherapy. The mechanisms underlying TIME formation and development over time are poorly understood. Glioblastoma (GBM) is a lethal primary brain cancer for which there are no curative treatments. GBMs are immunologically heterogeneous and impervious to checkpoint blockade immunotherapies. Utilizing clinically relevant genetic mouse models of GBM, we identified distinct immune landscapes associated with expression of EGFR wild-type and mutant EGFRvIII cancer driver mutations. Over time, accumulation of polymorphonuclear myeloid-derived suppressor cells (PMN-MDSC) was more pronounced in EGFRvIII-driven GBMs and was correlated with resistance to PD-1 and CTLA-4 combination checkpoint blockade immunotherapy. We determined that GBM-secreted CXCL1/2/3 and PMN-MDSC-expressed CXCR2 formed an axis regulating output of PMN-MDSCs from the bone marrow leading to systemic increase in these cells in the spleen and GBM tumor-draining lymph nodes. Pharmacologic targeting of this axis induced a systemic decrease in the numbers of PMN-MDSC, facilitated responses to PD-1 and CTLA-4 combination checkpoint blocking immunotherapy, and prolonged survival in mice bearing EGFRvIII-driven GBM. Our results uncover a relationship between cancer driver mutations, TIME composition, and sensitivity to checkpoint blockade in GBM and support the stratification of patients with GBM for checkpoint blockade therapy based on integrated genotypic and immunologic profiles.

GBMdeconvoluteR accurately infers proportions of neoplastic and immune cell populations from bulk glioblastoma transcriptomics data.

BACKGROUND: Characterizing and quantifying cell types within glioblastoma (GBM) tumors at scale will facilitate a better understanding of the association between the cellular landscape and tumor phenotypes or clinical correlates. We aimed to develop a tool that deconvolutes immune and neoplastic cells within the GBM tumor microenvironment from bulk RNA sequencing data. METHODS: We developed an IDH wild-type (IDHwt) GBM-specific single immune cell reference consisting of B cells, T-cells, NK-cells, microglia, tumor associated macrophages, monocytes, mast and DC cells. We used this alongside an existing neoplastic single cell-type reference for astrocyte-like, oligodendrocyte- and neuronal progenitor-like and mesenchymal GBM cancer cells to create both marker and gene signature matrix-based deconvolution tools. We applied single-cell resolution imaging mass cytometry (IMC) to ten IDHwt GBM samples, five paired primary and recurrent tumors, to determine which deconvolution approach performed best. RESULTS: Marker-based deconvolution using GBM-tissue specific markers was most accurate for both immune cells and cancer cells, so we packaged this approach as GBMdeconvoluteR. We applied GBMdeconvoluteR to bulk GBM RNAseq data from The Cancer Genome Atlas and recapitulated recent findings from multi-omics single cell studies with regards associations between mesenchymal GBM cancer cells and both lymphoid and myeloid cells. Furthermore, we expanded upon this to show that these associations are stronger in patients with worse prognosis. CONCLUSIONS: GBMdeconvoluteR accurately quantifies immune and neoplastic cell proportions in IDHwt GBM bulk RNA sequencing data and is accessible here: https://gbmdeconvoluter.leeds.ac.uk.

Extrachromosomal DNA (ecDNA): an origin of tumor heterogeneity, genomic remodeling, and drug resistance.

The genome of cancer cells contains circular extrachromosomal DNA (ecDNA) elements not found in normal cells. Analysis of clinical samples reveal they are common in most cancers and their presence indicates poor prognosis. They often contain enhancers and driver oncogenes that are highly expressed. The circular ecDNA topology leads to an open chromatin conformation and generates new gene regulatory interactions, including with distal enhancers. The absence of centromeres leads to random distribution of ecDNAs during cell division and genes encoded on them are transmitted in a non-mendelian manner. ecDNA can integrate into and exit from chromosomal DNA. The numbers of specific ecDNAs can change in response to treatment. This dynamic ability to remodel the cancer genome challenges long-standing fundamentals, providing new insights into tumor heterogeneity, cancer genome remodeling, and drug resistance.

Lineage-coupled clonal capture identifies clonal evolution mechanisms and vulnerabilities of BRAFV600E inhibition resistance in melanoma.

Targeted cancer therapies have revolutionized treatment but their efficacies are limited by the development of resistance driven by clonal evolution within tumors. We developed "CAPTURE", a single-cell barcoding approach to comprehensively trace clonal dynamics and capture live lineage-coupled resistant cells for in-depth multi-omics analysis and functional exploration. We demonstrate that heterogeneous clones, either preexisting or emerging from drug-tolerant persister cells, dominated resistance to vemurafenib in BRAFV600E melanoma. Further integrative studies uncovered diverse resistance mechanisms. This includes a previously unrecognized and clinically relevant mechanism, chromosome 18q21 gain, which leads to vulnerability of the cells to BCL2 inhibitor. We also identified targetable common dependencies of captured resistant clones, such as oxidative phosphorylation and E2F pathways. Our study provides new therapeutic insights into overcoming therapy resistance in BRAFV600E melanoma and presents a platform for exploring clonal evolution dynamics and vulnerabilities that can be applied to study treatment resistance in other cancers.

PDPN marks a subset of aggressive and radiation-resistant glioblastoma cells.

Treatment-resistant glioma stem cells are thought to propagate and drive growth of malignant gliomas, but their markers and our ability to target them specifically are not well understood. We demonstrate that podoplanin (PDPN) expression is an independent prognostic marker in gliomas across multiple independent patient cohorts comprising both high- and low-grade gliomas. Knockdown of PDPN radiosensitized glioma cell lines and glioma-stem-like cells (GSCs). Clonogenic assays and xenograft experiments revealed that PDPN expression was associated with radiotherapy resistance and tumor aggressiveness. We further demonstrate that knockdown of PDPN in GSCs in vivo is sufficient to improve overall survival in an intracranial xenograft mouse model. PDPN therefore identifies a subset of aggressive, treatment-resistant glioma cells responsible for radiation resistance and may serve as a novel therapeutic target.

Natural Coevolution of Tumor and Immunoenvironment in Glioblastoma.

Isocitrate dehydrogenase (IDH) wild-type glioblastoma (GBM) has a dismal prognosis. A better understanding of tumor evolution holds the key to developing more effective treatment. Here we study GBM's natural evolutionary trajectory by using rare multifocal samples. We sequenced 61,062 single cells from eight multifocal IDH wild-type primary GBMs and defined a natural evolution signature (NES) of the tumor. We show that the NES significantly associates with the activation of transcription factors that regulate brain development, including MYBL2 and FOSL2. Hypoxia is involved in inducing NES transition potentially via activation of the HIF1A-FOSL2 axis. High-NES tumor cells could recruit and polarize bone marrow-derived macrophages through activation of the FOSL2-ANXA1-FPR1/3 axis. These polarized macrophages can efficiently suppress T-cell activity and accelerate NES transition in tumor cells. Moreover, the polarized macrophages could upregulate CCL2 to induce tumor cell migration. SIGNIFICANCE: GBM progression could be induced by hypoxia via the HIF1A-FOSL2 axis. Tumor-derived ANXA1 is associated with recruitment and polarization of bone marrow-derived macrophages to suppress the immunoenvironment. The polarized macrophages promote tumor cell NES transition and migration. This article is highlighted in the In This Issue feature, p. 2711.

Extrachromosomal DNA amplifications in cancer.

Extrachromosomal DNA (ecDNA) amplification is an important driver alteration in cancer. It has been observed in most cancer types and is associated with worse patient outcome. The functional impact of ecDNA has been linked to its unique properties, such as its circular structure that is associated with altered chromatinization and epigenetic regulatory landscape, as well as its ability to randomly segregate during cell division, which fuels intercellular copy number heterogeneity. Recent investigations suggest that ecDNA is structurally more complex than previously anticipated and that it localizes to specialized nuclear bodies (hubs) and can act in trans as an enhancer for genes on other ecDNAs or chromosomes. In this Review, we synthesize what is currently known about how ecDNA is generated and how its genetic and epigenetic architecture affects proto-oncogene deregulation in cancer. We discuss how recently identified ecDNA functions may impact oncogenesis but also serve as new therapeutic vulnerabilities in cancer.

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.

Spatiotemporal dynamics of clonal selection and diversification in normal endometrial epithelium.

It has become evident that somatic mutations in cancer-associated genes accumulate in the normal endometrium, but spatiotemporal understanding of the evolution and expansion of mutant clones is limited. To elucidate the timing and mechanism of the clonal expansion of somatic mutations in cancer-associated genes in the normal endometrium, we sequence 1311 endometrial glands from 37 women. By collecting endometrial glands from different parts of the endometrium, we show that multiple glands with the same somatic mutations occupy substantial areas of the endometrium. We demonstrate that "rhizome structures", in which the basal glands run horizontally along the muscular layer and multiple vertical glands rise from the basal gland, originate from the same ancestral clone. Moreover, mutant clones detected in the vertical glands diversify by acquiring additional mutations. These results suggest that clonal expansions through the rhizome structures are involved in the mechanism by which mutant clones extend their territories. Furthermore, we show clonal expansions and copy neutral loss-of-heterozygosity events occur early in life, suggesting such events can be tolerated many years in the normal endometrium. Our results of the evolutionary dynamics of mutant clones in the human endometrium will lead to a better understanding of the mechanisms of endometrial regeneration during the menstrual cycle and the development of therapies for the prevention and treatment of endometrium-related diseases.

Live-Cell Imaging Shows Uneven Segregation of Extrachromosomal DNA Elements and Transcriptionally Active Extrachromosomal DNA Hubs in Cancer.

Oncogenic extrachromosomal DNA elements (ecDNA) play an important role in tumor evolution, but our understanding of ecDNA biology is limited. We determined the distribution of single-cell ecDNA copy number across patient tissues and cell line models and observed how cell-to-cell ecDNA frequency varies greatly. The exceptional intratumoral heterogeneity of ecDNA suggested ecDNA-specific replication and propagation mechanisms. To evaluate the transfer of ecDNA genetic material from parental to offspring cells during mitosis, we established the CRISPR-based ecTag method. ecTag leverages ecDNA-specific breakpoint sequences to tag ecDNA with fluorescent markers in living cells. Applying ecTag during mitosis revealed disjointed ecDNA inheritance patterns, enabling rapid ecDNA accumulation in individual cells. After mitosis, ecDNAs clustered into ecDNA hubs, and ecDNA hubs colocalized with RNA polymerase II, promoting transcription of cargo oncogenes. Our observations provide direct evidence for uneven segregation of ecDNA and shed new light on mechanisms through which ecDNAs contribute to oncogenesis. SIGNIFICANCE: ecDNAs are vehicles for oncogene amplification. The circular nature of ecDNA affords unique properties, such as mobility and ecDNA-specific replication and segregation behavior. We uncovered fundamental ecDNA properties by tracking ecDNAs in live cells, highlighting uneven and random segregation and ecDNA hubs that drive cargo gene transcription.See related commentary by Henssen, p. 293.This article is highlighted in the In This Issue feature, p. 275.

Single-cell multimodal glioma analyses identify epigenetic regulators of cellular plasticity and environmental stress response.

Glioma intratumoral heterogeneity enables adaptation to challenging microenvironments and contributes to therapeutic resistance. We integrated 914 single-cell DNA methylomes, 55,284 single-cell transcriptomes and bulk multi-omic profiles across 11 adult IDH mutant or IDH wild-type gliomas to delineate sources of intratumoral heterogeneity. We showed that local DNA methylation disorder is associated with cell-cell DNA methylation differences, is elevated in more aggressive tumors, links with transcriptional disruption and is altered during the environmental stress response. Glioma cells under in vitro hypoxic and irradiation stress increased local DNA methylation disorder and shifted cell states. We identified a positive association between genetic and epigenetic instability that was supported in bulk longitudinally collected DNA methylation data. Increased DNA methylation disorder associated with accelerated disease progression and recurrently selected DNA methylation changes were enriched for environmental stress response pathways. Our work identified an epigenetically facilitated adaptive stress response process and highlights the importance of epigenetic heterogeneity in shaping therapeutic outcomes.

Homozygous MTAP deletion in primary human glioblastoma is not associated with elevation of methylthioadenosine.

Homozygous deletion of methylthioadenosine phosphorylase (MTAP) in cancers such as glioblastoma represents a potentially targetable vulnerability. Homozygous MTAP-deleted cell lines in culture show elevation of MTAP's substrate metabolite, methylthioadenosine (MTA). High levels of MTA inhibit protein arginine methyltransferase 5 (PRMT5), which sensitizes MTAP-deleted cells to PRMT5 and methionine adenosyltransferase 2A (MAT2A) inhibition. While this concept has been extensively corroborated in vitro, the clinical relevance relies on exhibiting significant MTA accumulation in human glioblastoma. In this work, using comprehensive metabolomic profiling, we show that MTA secreted by MTAP-deleted cells in vitro results in high levels of extracellular MTA. We further demonstrate that homozygous MTAP-deleted primary glioblastoma tumors do not significantly accumulate MTA in vivo due to metabolism of MTA by MTAP-expressing stroma. These findings highlight metabolic discrepancies between in vitro models and primary human tumors that must be considered when developing strategies for precision therapies targeting glioblastoma with homozygous MTAP deletion.