Modeling epigenetic lesions that cause gliomas.

Epigenetic lesions that disrupt regulatory elements represent potential cancer drivers. However, we lack experimental models for validating their tumorigenic impact. Here, we model aberrations arising in isocitrate dehydrogenase-mutant gliomas, which exhibit DNA hypermethylation. We focus on a CTCF insulator near the PDGFRA oncogene that is recurrently disrupted by methylation in these tumors. We demonstrate that disruption of the syntenic insulator in mouse oligodendrocyte progenitor cells (OPCs) allows an OPC-specific enhancer to contact and induce Pdgfra, thereby increasing proliferation. We show that a second lesion, methylation-dependent silencing of the Cdkn2a tumor suppressor, cooperates with insulator loss in OPCs. Coordinate inactivation of the Pdgfra insulator and Cdkn2a drives gliomagenesis in vivo. Despite locus synteny, the insulator is CpG-rich only in humans, a feature that may confer human glioma risk but complicates mouse modeling. Our study demonstrates the capacity of recurrent epigenetic lesions to drive OPC proliferation in vitro and gliomagenesis in vivo.

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.

Inhibitory CD161 receptor identified in glioma-infiltrating T cells by single-cell analysis.

T cells are critical effectors of cancer immunotherapies, but little is known about their gene expression programs in diffuse gliomas. Here, we leverage single-cell RNA sequencing (RNA-seq) to chart the gene expression and clonal landscape of tumor-infiltrating T cells across 31 patients with isocitrate dehydrogenase (IDH) wild-type glioblastoma and IDH mutant glioma. We identify potential effectors of anti-tumor immunity in subsets of T cells that co-express cytotoxic programs and several natural killer (NK) cell genes. Analysis of clonally expanded tumor-infiltrating T cells further identifies the NK gene KLRB1 (encoding CD161) as a candidate inhibitory receptor. Accordingly, genetic inactivation of KLRB1 or antibody-mediated CD161 blockade enhances T cell-mediated killing of glioma cells in vitro and their anti-tumor function in vivo. KLRB1 and its associated transcriptional program are also expressed by substantial T cell populations in other human cancers. Our work provides an atlas of T cells in gliomas and highlights CD161 and other NK cell receptors as immunotherapy targets.

Integrated Proteogenomic Characterization across Major Histological Types of Pediatric Brain Cancer.

We report a comprehensive proteogenomics analysis, including whole-genome sequencing, RNA sequencing, and proteomics and phosphoproteomics profiling, of 218 tumors across 7 histological types of childhood brain cancer: low-grade glioma (n = 93), ependymoma (32), high-grade glioma (25), medulloblastoma (22), ganglioglioma (18), craniopharyngioma (16), and atypical teratoid rhabdoid tumor (12). Proteomics data identify common biological themes that span histological boundaries, suggesting that treatments used for one histological type may be applied effectively to other tumors sharing similar proteomics features. Immune landscape characterization reveals diverse tumor microenvironments across and within diagnoses. Proteomics data further reveal functional effects of somatic mutations and copy number variations (CNVs) not evident in transcriptomics data. Kinase-substrate association and co-expression network analysis identify important biological mechanisms of tumorigenesis. This is the first large-scale proteogenomics analysis across traditional histological boundaries to uncover foundational pediatric brain tumor biology and inform rational treatment selection.

Histone H3.3G34-Mutant Interneuron Progenitors Co-opt PDGFRA for Gliomagenesis.

Histone H3.3 glycine 34 to arginine/valine (G34R/V) mutations drive deadly gliomas and show exquisite regional and temporal specificity, suggesting a developmental context permissive to their effects. Here we show that 50% of G34R/V tumors (n = 95) bear activating PDGFRA mutations that display strong selection pressure at recurrence. Although considered gliomas, G34R/V tumors actually arise in GSX2/DLX-expressing interneuron progenitors, where G34R/V mutations impair neuronal differentiation. The lineage of origin may facilitate PDGFRA co-option through a chromatin loop connecting PDGFRA to GSX2 regulatory elements, promoting PDGFRA overexpression and mutation. At the single-cell level, G34R/V tumors harbor dual neuronal/astroglial identity and lack oligodendroglial programs, actively repressed by GSX2/DLX-mediated cell fate specification. G34R/V may become dispensable for tumor maintenance, whereas mutant-PDGFRA is potently oncogenic. Collectively, our results open novel research avenues in deadly tumors. G34R/V gliomas are neuronal malignancies where interneuron progenitors are stalled in differentiation by G34R/V mutations and malignant gliogenesis is promoted by co-option of a potentially targetable pathway, PDGFRA signaling.

Rapid Generation of Somatic Mouse Mosaics with Locus-Specific, Stably Integrated Transgenic Elements.

In situ transgenesis methods such as viruses and electroporation can rapidly create somatic transgenic mice but lack control over copy number, zygosity, and locus specificity. Here we establish mosaic analysis by dual recombinase-mediated cassette exchange (MADR), which permits stable labeling of mutant cells expressing transgenic elements from precisely defined chromosomal loci. We provide a toolkit of MADR elements for combination labeling, inducible and reversible transgene manipulation, VCre recombinase expression, and transgenesis of human cells. Further, we demonstrate the versatility of MADR by creating glioma models with mixed reporter-identified zygosity or with "personalized" driver mutations from pediatric glioma. MADR is extensible to thousands of existing mouse lines, providing a flexible platform to democratize the generation of somatic mosaic mice. VIDEO ABSTRACT.

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.

SMARCA4 vulnerability in H3K27M midline glioma: A silver bullet for a lethal disease.

To investigate epigenetic dependencies and identify therapeutic vulnerabilities, Mo et al.1 and Panditharatna et al.2 performed CRISPR screens and show that deadly H3K27M gliomas are dependent on mammalian BAF (SWI/SNF) chromatin remodeling complex.

Promoter and enhancer RNAs regulate chromatin reorganization and activation of miR-10b/HOXD locus, and neoplastic transformation in glioma.

miR-10b is silenced in normal neuroglial cells of the brain but commonly activated in glioma, where it assumes an essential tumor-promoting role. We demonstrate that the entire miR-10b-hosting HOXD locus is activated in glioma via the cis-acting mechanism involving 3D chromatin reorganization and CTCF-cohesin-mediated looping. This mechanism requires two interacting lncRNAs, HOXD-AS2 and LINC01116, one associated with HOXD3/HOXD4/miR-10b promoter and another with the remote enhancer. Knockdown of either lncRNA in glioma cells alters CTCF and cohesin binding, abolishes chromatin looping, inhibits the expression of all genes within HOXD locus, and leads to glioma cell death. Conversely, in cortical astrocytes, enhancer activation is sufficient for HOXD/miR-10b locus reorganization, gene derepression, and neoplastic cell transformation. LINC01116 RNA is essential for this process. Our results demonstrate the interplay of two lncRNAs in the chromatin folding and concordant regulation of miR-10b and multiple HOXD genes normally silenced in astrocytes and triggering the neoplastic glial transformation.

GD2-CAR T cell therapy for H3K27M-mutated diffuse midline gliomas.

Diffuse intrinsic pontine glioma (DIPG) and other H3K27M-mutated diffuse midline gliomas (DMGs) are universally lethal paediatric tumours of the central nervous system1. We have previously shown that the disialoganglioside GD2 is highly expressed on H3K27M-mutated glioma cells and have demonstrated promising preclinical efficacy of GD2-directed chimeric antigen receptor (CAR) T cells2, providing the rationale for a first-in-human phase I clinical trial (NCT04196413). Because CAR T cell-induced brainstem inflammation can result in obstructive hydrocephalus, increased intracranial pressure and dangerous tissue shifts, neurocritical care precautions were incorporated. Here we present the clinical experience from the first four patients with H3K27M-mutated DIPG or spinal cord DMG treated with GD2-CAR T cells at dose level 1 (1 × 106 GD2-CAR T cells per kg administered intravenously). Patients who exhibited clinical benefit were eligible for subsequent GD2-CAR T cell infusions administered intracerebroventricularly3. Toxicity was largely related to the location of the tumour and was reversible with intensive supportive care. On-target, off-tumour toxicity was not observed. Three of four patients exhibited clinical and radiographic improvement. Pro-inflammatory cytokine levels were increased in the plasma and cerebrospinal fluid. Transcriptomic analyses of 65,598 single cells from CAR T cell products and cerebrospinal fluid elucidate heterogeneity in response between participants and administration routes. These early results underscore the promise of this therapeutic approach for patients with H3K27M-mutated DIPG or spinal cord DMG.

High-Resolution In Vivo Identification of miRNA Targets by Halo-Enhanced Ago2 Pull-Down.

The identification of microRNA (miRNA) targets by Ago2 crosslinking-immunoprecipitation (CLIP) methods has provided major insights into the biology of this important class of non-coding RNAs. However, these methods are technically challenging and not easily applicable to an in vivo setting. To overcome these limitations and facilitate the investigation of miRNA functions in vivo, we have developed a method based on a genetically engineered mouse harboring a conditional Halo-Ago2 allele expressed from the endogenous Ago2 locus. By using a resin conjugated to the HaloTag ligand, Ago2-miRNA-mRNA complexes can be purified from cells and tissues expressing the endogenous Halo-Ago2 allele. We demonstrate the reproducibility and sensitivity of this method in mouse embryonic stem cells, developing embryos, adult tissues, and autochthonous mouse models of human brain and lung cancers. This method and the datasets we have generated will facilitate the characterization of miRNA-mRNA networks in vivo under physiological and pathological conditions.

Proline Hydroxylation Primes Protein Kinases for Autophosphorylation and Activation.

Activation of dual-specificity tyrosine-phosphorylation-regulated kinases 1A and 1B (DYRK1A and DYRK1B) requires prolyl hydroxylation by PHD1 prolyl hydroxylase. Prolyl hydroxylation of DYRK1 initiates a cascade of events leading to the release of molecular constraints on von Hippel-Lindau (VHL) ubiquitin ligase tumor suppressor function. However, the proline residue of DYRK1 targeted by hydroxylation and the role of prolyl hydroxylation in tyrosine autophosphorylation of DYRK1 are unknown. We found that a highly conserved proline in the CMGC insert of the DYRK1 kinase domain is hydroxylated by PHD1, and this event precedes tyrosine autophosphorylation. Mutation of the hydroxylation acceptor proline precludes tyrosine autophosphorylation and folding of DYRK1, resulting in a kinase unable to preserve VHL function and lacking glioma suppression activity. The consensus proline sequence is shared by most CMGC kinases, and prolyl hydroxylation is essential for catalytic activation. Thus, formation of prolyl-hydroxylated intermediates is a novel mechanism of kinase maturation and likely a general mechanism of regulation of CMGC kinases in eukaryotes.

Accurate estimation of pathway activity in single cells for clustering and differential analysis.

Inferring which and how biological pathways and gene sets change is a key question in many studies that utilize single-cell RNA sequencing. Typically, these questions are addressed by quantifying the enrichment of known gene sets in lists of genes derived from global analysis. Here we offer SiPSiC, a new method to infer pathway activity in every single cell. This allows more sensitive differential analysis and utilization of pathway scores to cluster cells and compute UMAP or other similar projections. We apply our method to COVID-19, lung adenocarcinoma and glioma data sets, and demonstrate its utility. SiPSiC analysis results are consistent with findings reported in previous studies in many cases, but SiPSiC also reveals the differential activity of novel pathways, enabling us to suggest new mechanisms underlying the pathophysiology of these diseases and demonstrating SiPSiC's high accuracy and sensitivity in detecting biological function and traits. In addition, we demonstrate how it can be used to better classify cells based on activity of biological pathways instead of single genes and its ability to overcome patient-specific artifacts.

The Power of Human Cancer Genetics as Revealed by Low-Grade Gliomas.

The human brain contains a vast number of cells and shows extraordinary cellular diversity to facilitate the many cognitive and automatic commands governing our bodily functions. This complexity arises partly from large-scale structural variations in the genome, evolutionary processes to increase brain size, function, and cognition. Not surprisingly given recent technical advances, low-grade gliomas (LGGs), which arise from the glia (the most abundant cell type in the brain), have undergone a recent revolution in their classification and therapy, especially in the pediatric setting. Next-generation sequencing has uncovered previously unappreciated diverse LGG entities, unraveling genetic subgroups and multiple molecular alterations and altered pathways, including many amenable to therapeutic targeting. In this article we review these novel entities, in which oncogenic processes show striking age-related neuroanatomical specificity (highlighting their close interplay with development); the opportunities they provide for targeted therapies, some of which are already practiced at the bedside; and the challenges of implementing molecular pathology in the clinic.

Malignant glioma: lessons from genomics, mouse models, and stem cells.

Eighty percent of malignant tumors that develop in the central nervous system are malignant gliomas, which are essentially incurable. Here, we discuss how recent sequencing studies are identifying unexpected drivers of gliomagenesis, including mutations in isocitrate dehydrogenase 1 and the NF-κB pathway, and how genome-wide analyses are reshaping the classification schemes for tumors and enhancing prognostic value of molecular markers. We discuss the controversies surrounding glioma stem cells and explore how the integration of new molecular data allows for the generation of more informative animal models to advance our knowledge of glioma's origin, progression, and treatment.

A vaccine targeting mutant IDH1 in newly diagnosed glioma.

Mutated isocitrate dehydrogenase 1 (IDH1) defines a molecularly distinct subtype of diffuse glioma1-3. The most common IDH1 mutation in gliomas affects codon 132 and encodes IDH1(R132H), which harbours a shared clonal neoepitope that is presented on major histocompatibility complex (MHC) class II4,5. An IDH1(R132H)-specific peptide vaccine (IDH1-vac) induces specific therapeutic T helper cell responses that are effective against IDH1(R132H)+ tumours in syngeneic MHC-humanized mice4,6-8. Here we describe a multicentre, single-arm, open-label, first-in-humans phase I trial that we carried out in 33 patients with newly diagnosed World Health Organization grade 3 and 4 IDH1(R132H)+ astrocytomas (Neurooncology Working Group of the German Cancer Society trial 16 (NOA16), ClinicalTrials.gov identifier NCT02454634). The trial met its primary safety endpoint, with vaccine-related adverse events restricted to grade 1. Vaccine-induced immune responses were observed in 93.3% of patients across multiple MHC alleles. Three-year progression-free and death-free rates were 0.63 and 0.84, respectively. Patients with immune responses showed a two-year progression-free rate of 0.82. Two patients without an immune response showed tumour progression within two years of first diagnosis. A mutation-specificity score that incorporates the duration and level of vaccine-induced IDH1(R132H)-specific T cell responses was associated with intratumoral presentation of the IDH1(R132H) neoantigen in pre-treatment tumour tissue. There was a high frequency of pseudoprogression, which indicates intratumoral inflammatory reactions. Pseudoprogression was associated with increased vaccine-induced peripheral T cell responses. Combined single-cell RNA and T cell receptor sequencing showed that tumour-infiltrating CD40LG+ and CXCL13+ T helper cell clusters in a patient with pseudoprogression were dominated by a single IDH1(R132H)-reactive T cell receptor.

Histone Variant and Cell Context Determine H3K27M Reprogramming of the Enhancer Landscape and Oncogenic State.

Development of effective targeted cancer therapies is fundamentally limited by our molecular understanding of disease pathogenesis. Diffuse intrinsic pontine glioma (DIPG) is a fatal malignancy of the childhood pons characterized by a unique substitution to methionine in histone H3 at lysine 27 (H3K27M) that results in globally altered epigenetic marks and oncogenic transcription. Through primary DIPG tumor characterization and isogenic oncohistone expression, we show that the same H3K27M mutation displays distinct modes of oncogenic reprogramming and establishes distinct enhancer architecture depending upon both the variant of histone H3 and the cell context in which the mutation occurs. Compared with non-malignant pediatric pontine tissue, we identify and functionally validate both shared and variant-specific pathophysiology. Altogether, we provide a powerful resource of epigenomic data in 25 primary DIPG samples and 5 rare normal pediatric pontine tissue samples, revealing clinically relevant functional distinctions previously unidentified in DIPG.

Rapid detection of IDH mutations in gliomas by intraoperative mass spectrometry.

The development and performance of two mass spectrometry (MS) workflows for the intraoperative diagnosis of isocitrate dehydrogenase (IDH) mutations in glioma is implemented by independent teams at Mayo Clinic, Jacksonville, and Huashan Hospital, Shanghai. The infiltrative nature of gliomas makes rapid diagnosis necessary to guide the extent of surgical resection of central nervous system (CNS) tumors. The combination of tissue biopsy and MS analysis used here satisfies this requirement. The key feature of both described methods is the use of tandem MS to measure the oncometabolite 2-hydroxyglutarate (2HG) relative to endogenous glutamate (Glu) to characterize the presence of mutant tumor. The experiments i) provide IDH mutation status for individual patients and ii) demonstrate a strong correlation of 2HG signals with tumor infiltration. The measured ratio of 2HG to Glu correlates with IDH-mutant (IDH-mut) glioma (P < 0.0001) in the tumor core data of both teams. Despite using different ionization methods and different mass spectrometers, comparable performance in determining IDH mutations from core tumor biopsies was achieved with sensitivities, specificities, and accuracies all at 100%. None of the 31 patients at Mayo Clinic or the 74 patients at Huashan Hospital were misclassified when analyzing tumor core biopsies. Robustness of the methodology was evaluated by postoperative re-examination of samples. Both teams noted the presence of high concentrations of 2HG at surgical margins, supporting future use of intraoperative MS to monitor for clean surgical margins. The power of MS diagnostics is shown in resolving contradictory clinical features, e.g., in distinguishing gliosis from IDH-mut glioma.

Automethylation of PRC2 promotes H3K27 methylation and is impaired in H3K27M pediatric glioma.

The histone methyltransferase activity of PRC2 is central to the formation of H3K27me3-decorated facultative heterochromatin and gene silencing. In addition, PRC2 has been shown to automethylate its core subunits, EZH1/EZH2 and SUZ12. Here, we identify the lysine residues at which EZH1/EZH2 are automethylated with EZH2-K510 and EZH2-K514 being the major such sites in vivo. Automethylated EZH2/PRC2 exhibits a higher level of histone methyltransferase activity and is required for attaining proper cellular levels of H3K27me3. While occurring independently of PRC2 recruitment to chromatin, automethylation promotes PRC2 accessibility to the histone H3 tail. Intriguingly, EZH2 automethylation is significantly reduced in diffuse intrinsic pontine glioma (DIPG) cells that carry a lysine-to-methionine substitution in histone H3 (H3K27M), but not in cells that carry either EZH2 or EED mutants that abrogate PRC2 allosteric activation, indicating that H3K27M impairs the intrinsic activity of PRC2. Our study demonstrates a PRC2 self-regulatory mechanism through its EZH1/2-mediated automethylation activity.