Clinically relevant molecular hallmarks of PFA ependymomas display intratumoral heterogeneity and correlate with tumor morphology.

Posterior fossa type A (PF-EPN-A, PFA) ependymoma are aggressive tumors that mainly affect children and have a poor prognosis. Histopathology shows significant intratumoral heterogeneity, ranging from loose tissue to often sharply demarcated, extremely cell-dense tumor areas. To determine molecular differences in morphologically different areas and to understand their clinical significance, we analyzed 113 PF-EPN-A samples, including 40 corresponding relapse samples. Cell-dense areas ranged from 0 to 100% of the tumor area and displayed a higher proportion of proliferating tumor cells (p < 0.01). Clinically, cell density was associated with poor progression-free and overall survival (pPFS = 0.0026, pOS < 0.01). Molecularly, tumor areas with low and high cell density showed diverging DNA methylation profiles regarding their similarity to distinct previously discovered PF-EPN-A subtypes in 9/21 cases. Prognostically relevant chromosomal changes at 1q and 6q showed spatial heterogeneity within single tumors and were significantly enriched in cell-dense tumor areas as shown by single-cell RNA (scRNA)-sequencing as well as copy number profiling and fluorescence in situ hybridization (FISH) analyses of different tumor areas. Finally, spatial transcriptomics revealed cell-dense areas of different tumors to be more similar than various different areas of the same tumor. High-density areas distinctly overexpressed genes encoding histone proteins, WNT5A, TGFB1, or IGF2. Relapsing tumors displayed a higher proportion of cell-dense areas (p = 0.036), a change in PF-EPN-A methylation subtypes (13/32 patients), and novel chromosome 1q gains and 6q losses (12/32 cases) compared to corresponding primary tumors. Our data suggest that PF-EPN-A ependymomas habor a previously unrecognized intratumoral heterogeneity with clinical implications, which has to be accounted for when selecting diagnostic material, inter alia, by histological evaluation of the proportion of cell-dense areas.

A road map for the treatment of pediatric diffuse midline glioma.

Recent clinical trials for H3K27-altered diffuse midline gliomas (DMGs) have shown much promise. We present a consensus roadmap and identify three major barriers: (1) refinement of experimental models to include immune and brain-specific components; (2) collaboration among researchers, clinicians, and industry to integrate patient-derived data through sharing, transparency, and regulatory considerations; and (3) streamlining clinical efforts including biopsy, CNS-drug delivery, endpoint determination, and response monitoring. We highlight the importance of comprehensive collaboration to advance the understanding, diagnostics, and therapeutics for DMGs.

Diffuse midline glioma invasion and metastasis rely on cell-autonomous signaling.

BACKGROUND: Diffuse midline gliomas (DMG) are pediatric tumors with negligible 2-year survival after diagnosis characterized by their ability to infiltrate the central nervous system. In the hope of controlling the local growth and slowing the disease, all patients receive radiotherapy. However, distant progression occurs frequently in DMG patients. Current clues as to what causes tumor infiltration circle mainly around the tumor microenvironment, but there are currently no known determinants to predict the degree of invasiveness. METHODS: In this study, we use patient-derived glioma stem cells (GSCs) to create patient-specific 3D avatars to model interindividual invasion and elucidate the cellular supporting mechanisms. RESULTS: We show that GSC models in 3D mirror the invasive behavior of the parental tumors, thus proving the ability of DMG to infiltrate as an autonomous characteristic of tumor cells. Furthermore, we distinguished 2 modes of migration, mesenchymal and ameboid-like, and associated the ameboid-like modality with GSCs derived from the most invasive tumors. Using transcriptomics of both organoids and primary tumors, we further characterized the invasive ameboid-like tumors as oligodendrocyte progenitor-like, with highly contractile cytoskeleton and reduced adhesion ability driven by crucial over-expression of bone morphogenetic pathway 7 (BMP7). Finally, we deciphered MEK, ERK, and Rho/ROCK kinases activated downstream of the BMP7 stimulation as actionable targets controlling tumor cell motility. CONCLUSIONS: Our findings identify 2 new therapeutic avenues. First, patient-derived GSCs represent a predictive tool for patient stratification in order to adapt irradiation strategies. Second, autocrine and short-range BMP7-related signaling becomes a druggable target to prevent DMG spread and metastasis.

An ERK5-PFKFB3 axis regulates glycolysis and represents a therapeutic vulnerability in pediatric diffuse midline glioma.

Metabolic reprogramming in pediatric diffuse midline glioma is driven by gene expression changes induced by the hallmark histone mutation H3K27M, which results in aberrantly permissive activation of oncogenic signaling pathways. Previous studies of diffuse midline glioma with altered H3K27 (DMG-H3K27a) have shown that the RAS pathway, specifically through its downstream kinase, extracellular-signal-related kinase 5 (ERK5), is critical for tumor growth. Further downstream effectors of ERK5 and their role in DMG-H3K27a metabolic reprogramming have not been explored. We establish that ERK5 is a critical regulator of cell proliferation and glycolysis in DMG-H3K27a. We demonstrate that ERK5 mediates glycolysis through activation of transcription factor MEF2A, which subsequently modulates expression of glycolytic enzyme PFKFB3. We show that in vitro and mouse models of DMG-H3K27a are sensitive to the loss of PFKFB3. Multi-targeted drug therapy against the ERK5-PFKFB3 axis, such as with small-molecule inhibitors, may represent a promising therapeutic approach in patients with pediatric diffuse midline glioma.

TIM-3 blockade in diffuse intrinsic pontine glioma models promotes tumor regression and antitumor immune memory.

Diffuse intrinsic pontine glioma (DIPG) is an aggressive brain stem tumor and the leading cause of pediatric cancer-related death. To date, these tumors remain incurable, underscoring the need for efficacious therapies. In this study, we demonstrate that the immune checkpoint TIM-3 (HAVCR2) is highly expressed in both tumor cells and microenvironmental cells, mainly microglia and macrophages, in DIPG. We show that inhibition of TIM-3 in syngeneic models of DIPG prolongs survival and produces long-term survivors free of disease that harbor immune memory. This antitumor effect is driven by the direct effect of TIM-3 inhibition in tumor cells, the coordinated action of several immune cell populations, and the secretion of chemokines/cytokines that create a proinflammatory tumor microenvironment favoring a potent antitumor immune response. This work uncovers TIM-3 as a bona fide target in DIPG and supports its clinical translation.

Mapping pediatric brain tumors to their origins in the developing cerebellum.

BACKGROUND: Distinguishing the cellular origins of childhood brain tumors is key for understanding tumor initiation and identifying lineage-restricted, tumor-specific therapeutic targets. Previous strategies to map the cell-of-origin typically involved comparing human tumors to murine embryonal tissues, which is potentially limited due to species-specific differences. The aim of this study was to unravel the cellular origins of the 3 most common pediatric brain tumors, ependymoma, pilocytic astrocytoma, and medulloblastoma, using a developing human cerebellar atlas. METHODS: We used a single-nucleus atlas of the normal developing human cerebellum consisting of 176 645 cells as a reference for an in-depth comparison to 4416 bulk and single-cell transcriptome tumor datasets, using gene set variation analysis, correlation, and single-cell matching techniques. RESULTS: We find that the astroglial cerebellar lineage is potentially the origin for posterior fossa ependymomas. We propose that infratentorial pilocytic astrocytomas originate from the oligodendrocyte lineage and MHC II genes are specifically enriched in these tumors. We confirm that SHH and Group 3/4 medulloblastomas originate from the granule cell and unipolar brush cell lineages. Radiation-induced gliomas stem from cerebellar glial lineages and demonstrate distinct origins from the primary medulloblastoma. We identify tumor genes that are expressed in the cerebellar lineage of origin, and genes that are tumor specific; both gene sets represent promising therapeutic targets for future study. CONCLUSION: Based on our results, individual cells within a tumor may resemble different cell types along a restricted developmental lineage. Therefore, we suggest that tumors can arise from multiple cellular states along the cerebellar "lineage of origin."

A comprehensive genomic study of 390 H3F3A-mutant pediatric and adult diffuse high-grade gliomas, CNS WHO grade 4.

Malignant brain tumors, known as H3K27-altered diffuse midline glioma (DMG) and H3G34-mutant diffuse hemispheric glioma (DHG), can affect individuals of all ages and are classified as CNS WHO grade 4. We comprehensively characterized 390 H3F3A-mutant diffuse gliomas (201 females, 189 males) arising in pediatric patients (under 20 years old) and adults (20 years and older) evaluated by the CGP program at Foundation Medicine between 2013 and 2020. We assessed information from pathology reports, histopathology review, and clinical data. The cohort included 304 H3K27M-mutant DMG (156 females, 148 males) and 86 H3G34-mutant DHG (45 females, 41 males). Median patient age was 20 years (1-74 years). The frequency of H3K27M-mutant DMG was similar in both pediatric and adult patients in our cohort-48.6% of the patients were over 20 years old, 31.5% over 30, and 18% over 40 at initial diagnosis. FGFR1 hotspot point mutations (N546K and K656E) were exclusively identified in H3K27M-mutant DMG tumors (64/304, 21%; p = 0.0001); these tend to occur in older patients (median age: 32.5 years) and mainly arose in the diencephalon. H3K27M-mutant DMG had higher rates of mutations in NF1 (31.0 vs 8.1%; p = 0.0001) and PIK3CA/PIK3R1 (27.9% vs 15.1%; p = 0.016) compared to H3G34-mutant DHG. However, H3G34-mutant DHG had higher rates of targetable alterations in cell-cycle pathway genes (CDK4 and CDK6 amplification; CDKN2A/B deletion) (27.0 vs 9.0%). Potentially targetable PDGFRA alterations were identified in ~ 20% of both H3G34-mutant DHG and H3K27M-mutant DMG. Overall, in the present study H3K27M-mutant DMG occurred at similar rates in both adult and patient patients. Through our analysis, we were able to identify molecular features characteristic of DMG and DHG. By identifying the recurrent co-mutations including actionable FGFR1 point mutations found in nearly one-third of H3K27M-mutant DMG in young adults, our findings can inform clinical translational studies, patient diagnosis, and clinical trial design.

Why haven't we solved intracranial pediatric ependymoma? Current questions and barriers to treatment advances.

Pediatric intracranial ependymoma has seen a recent exponential expansion of biological findings, rapidly dividing the diagnosis into several subgroups, each with specific molecular and clinical characteristics. While such subdivision may complicate clinical conclusions from historical trials, this knowledge also provides an opportunity for interrogating the major clinical and biological questions preventing near-term translation into effective therapy for children with ependymoma. In this article, we briefly review some of the most critical clinical questions facing both patient management and the construct of future trials in childhood ependymoma, as well as explore some of the current barriers to efficient translation of preclinical discovery to the clinic.

Common molecular features of H3K27M DMGs and PFA ependymomas map to hindbrain developmental pathways.

Globally decreased histone 3, lysine 27 tri-methylation (H3K27me3) is a hallmark of H3K27-altered diffuse midline gliomas (DMGs) and group-A posterior fossa ependymomas (PFAs). H3K27-altered DMGs are largely characterized by lysine-to-methionine mutations in histone 3 at position 27 (H3K27M). Most PFAs overexpress EZH inhibitory protein (EZHIP), which possesses a region of similarity to the mutant H3K27M. Both H3K27M and EZHIP inhibit the function of the polycomb repressive complex 2 (PRC2) responsible for H3K27me3 deposition. These tumors often arise in neighboring regions of the brainstem and posterior fossa. In rare cases PFAs harbor H3K27M mutations, and DMGs overexpress EZHIP. These findings together raise the possibility that certain cell populations in the developing hindbrain/posterior fossa are especially sensitive to modulation of H3K27me3 states. We identified shared molecular features by comparing genomic, bulk transcriptomic, chromatin-based profiles, and single-cell RNA-sequencing (scRNA-seq) data from the two tumor classes. Our approach demonstrated that 1q gain, a key biomarker in PFAs, is prognostic in H3.1K27M, but not H3.3K27M gliomas. Conversely, Activin A Receptor Type 1 (ACVR1), which is associated with mutations in H3.1K27M gliomas, is overexpressed in a subset of PFAs with poor outcome. Despite diffuse H3K27me3 reduction, previous work shows that both tumors maintain genomic H3K27me3 deposition at select sites. We demonstrate heterogeneity in shared patterns of residual H3K27me3 for both tumors that largely segregated with inferred anatomic tumor origins and progenitor populations of tumor cells. In contrast, analysis of genes linked to H3K27 acetylation (H3K27ac)-marked enhancers showed higher expression in astrocytic-like tumor cells. Finally, common H3K27me3-marked genes mapped closely to expression patterns in the human developing hindbrain. Overall, our data demonstrate developmentally relevant molecular similarities between PFAs and H3K27M DMGs and support the overall hypothesis that deregulated mechanisms of hindbrain development are central to the biology of both tumors.

Neuronal-Activity Dependent Mechanisms of Small Cell Lung Cancer Progression.

Neural activity is increasingly recognized as a critical regulator of cancer growth. In the brain, neuronal activity robustly influences glioma growth both through paracrine mechanisms and through electrochemical integration of malignant cells into neural circuitry via neuron-to-glioma synapses, while perisynaptic neurotransmitter signaling drives breast cancer brain metastasis growth. Outside of the CNS, innervation of tumors such as prostate, breast, pancreatic and gastrointestinal cancers by peripheral nerves similarly regulates cancer progression. However, the extent to which the nervous system regulates lung cancer progression, either in the lung or when metastatic to brain, is largely unexplored. Small cell lung cancer (SCLC) is a lethal high-grade neuroendocrine tumor that exhibits a strong propensity to metastasize to the brain. Here we demonstrate that, similar to glioma, metastatic SCLC cells in the brain co-opt neuronal activity-regulated mechanisms to stimulate growth and progression. Optogenetic stimulation of cortical neuronal activity drives proliferation and invasion of SCLC brain metastases. In the brain, SCLC cells exhibit electrical currents and consequent calcium transients in response to neuronal activity, and direct SCLC cell membrane depolarization is sufficient to promote the growth of SCLC tumors. In the lung, vagus nerve transection markedly inhibits primary lung tumor formation, progression and metastasis, highlighting a critical role for innervation in overall SCLC initiation and progression. Taken together, these studies illustrate that neuronal activity plays a crucial role in dictating SCLC pathogenesis in both primary and metastatic sites.

Activation of Hedgehog signaling by the oncogenic RELA fusion reveals a primary cilia-dependent vulnerability in supratentorial ependymoma.

BACKGROUND: Supratentorial RELA fusion (ST-RELA) ependymomas (EPNs) are resistant tumors without an approved chemotherapeutic treatment. Unfortunately, the molecular mechanisms that lead to chemoresistance traits of ST-RELA remain elusive. The aim of this study was to assess RELA fusion-dependent signaling modules, specifically the role of the Hedgehog (Hh) pathway as a novel targetable vulnerability in ST-RELA. METHODS: Gene expression was analyzed in EPN from patient cohorts, by microarray, RNA-seq, qRT-PCR, and scRNA-seq. Inhibitors against Smoothened (SMO) (Sonidegib) and Aurora kinase A (AURKA) (Alisertib) were evaluated. Protein expression, primary cilia formation, and drug effects were assessed by immunoblot, immunofluorescence, and immunohistochemistry. RESULTS: Hh components were selectively overexpressed in EPNs induced by the RELA fusion. Single-cell analysis showed that the Hh signature was primarily confined to undifferentiated, stem-like cell subpopulations. Sonidegib exhibited potent growth-inhibitory effects on ST-RELA cells, suggesting a key role in active Hh signaling; importantly, the effect of Sonidegib was reversed by primary cilia loss. We, thus, tested the effect of AURKA inhibition by Alisertib, to induce cilia stabilization/reassembly. Strikingly, Alisertib rescued ciliogenesis and synergized with Sonidegib in killing ST-RELA cells. Using a xenograft model, we show that cilia loss is a mechanism for acquiring resistance to the inhibitory effect of Sonidegib. However, Alisertib fails to rescue cilia and highlights the need for other strategies to promote cilia reassembly, for treating ST-RELA tumors. CONCLUSION: Our study reveals a crucial role for the Hh pathway in ST-RELA tumor growth, and suggests that rescue of primary cilia represents a vulnerability of the ST-RELA EPNs.

Therapeutic targeting of prenatal pontine ID1 signaling in diffuse midline glioma.

BACKGROUND: Diffuse midline gliomas (DMG) are highly invasive brain tumors with rare survival beyond two years past diagnosis and limited understanding of the mechanism behind tumor invasion. Previous reports demonstrate upregulation of the protein ID1 with H3K27M and ACVR1 mutations in DMG, but this has not been confirmed in human tumors or therapeutically targeted. METHODS: Whole exome, RNA, and ChIP-sequencing was performed on the ID1 locus in DMG tissue. Scratch-assay migration and transwell invasion assays of cultured cells were performed following shRNA-mediated ID1-knockdown. In vitro and in vivo genetic and pharmacologic [cannabidiol (CBD)] inhibition of ID1 on DMG tumor growth was assessed. Patient-reported CBD dosing information was collected. RESULTS: Increased ID1 expression in human DMG and in utero electroporation (IUE) murine tumors is associated with H3K27M mutation and brainstem location. ChIP-sequencing indicates ID1 regulatory regions are epigenetically active in human H3K27M-DMG tumors and prenatal pontine cells. Higher ID1-expressing astrocyte-like DMG cells share a transcriptional program with oligo/astrocyte-precursor cells (OAPCs) from the developing human brain and demonstrate upregulation of the migration regulatory protein SPARCL1. Genetic and pharmacologic (CBD) suppression of ID1 decreases tumor cell invasion/migration and tumor growth in H3.3/H3.1K27M PPK-IUE and human DIPGXIIIP* in vivo models of pHGG. The effect of CBD on cell proliferation appears to be non-ID1 mediated. Finally, we collected patient-reported CBD treatment data, finding that a clinical trial to standardize dosing may be beneficial. CONCLUSIONS: H3K27M-mediated re-activation of ID1 in DMG results in a SPARCL1+ migratory transcriptional program that is therapeutically targetable with CBD.

Characterizing the biology of primary brain tumors and their microenvironment via single-cell profiling methods.

Genomic and transcriptional heterogeneity is prevalent among the most common and aggressive primary brain tumors in children and adults. Over the past 20 years, advances in bioengineering, biochemistry and bioinformatics have enabled the development of an array of techniques to study tumor biology at single-cell resolution. The application of these techniques to study primary brain tumors has helped advance our understanding of their intra-tumoral heterogeneity and uncover new insights regarding their co-option of developmental programs and signaling from their microenvironment to promote tumor proliferation and invasion. These insights are currently being harnessed to develop new therapeutic approaches. Here we provide an overview of current single-cell techniques and discuss relevant biology and therapeutic insights uncovered by their application to primary brain tumors in children and adults.

The landscape of tumor cell states and spatial organization in H3-K27M mutant diffuse midline glioma across age and location.

Histone 3 lysine27-to-methionine (H3-K27M) mutations most frequently occur in diffuse midline gliomas (DMGs) of the childhood pons but are also increasingly recognized in adults. Their potential heterogeneity at different ages and midline locations is vastly understudied. Here, through dissecting the single-cell transcriptomic, epigenomic and spatial architectures of a comprehensive cohort of patient H3-K27M DMGs, we delineate how age and anatomical location shape glioma cell-intrinsic and -extrinsic features in light of the shared driver mutation. We show that stem-like oligodendroglial precursor-like cells, present across all clinico-anatomical groups, display varying levels of maturation dependent on location. We reveal a previously underappreciated relationship between mesenchymal cancer cell states and age, linked to age-dependent differences in the immune microenvironment. Further, we resolve the spatial organization of H3-K27M DMG cell populations and identify a mitotic oligodendroglial-lineage niche. Collectively, our study provides a powerful framework for rational modeling and therapeutic interventions.

Immunotherapy approaches for the treatment of diffuse midline gliomas.

Diffuse midline gliomas (DMG) are a highly aggressive and universally fatal subgroup of pediatric tumors responsible for the majority of childhood brain tumor deaths. Median overall survival is less than 12 months with a 90% mortality rate at 2 years from diagnosis. Research into the underlying tumor biology and numerous clinical trials have done little to change the invariably poor prognosis. Continued development of novel, efficacious therapeutic options for DMGs remains a critically important area of active investigation. Given that DMGs are not amenable to surgical resection, have only limited response to radiation, and are refractory to traditional chemotherapy, immunotherapy has emerged as a promising alternative treatment modality. This review summarizes the various immunotherapy-based treatments for DMG as well as their specific limitations. We explore the use of cell-based therapies, oncolytic virotherapy or immunovirotherapy, immune checkpoint inhibition, and immunomodulatory vaccination strategies, and highlight the recent clinical success of anti-GD2 CAR-T therapy in diffuse intrinsic pontine glioma (DIPG) patients. Finally, we address the challenges faced in translating preclinical and early phase clinical trial data into effective standardized treatment for DMG patients.

A druggable addiction to de novo pyrimidine biosynthesis in diffuse midline glioma.

Diffuse midline glioma (DMG) is a uniformly fatal pediatric cancer driven by oncohistones that do not readily lend themselves to drug development. To identify druggable targets for DMG, we conducted a genome-wide CRISPR screen that reveals a DMG selective dependency on the de novo pathway for pyrimidine biosynthesis. This metabolic vulnerability reflects an elevated rate of uridine/uracil degradation that depletes DMG cells of substrates for the alternate salvage pyrimidine biosynthesis pathway. A clinical stage inhibitor of DHODH (rate-limiting enzyme in the de novo pathway) diminishes uridine-5'-phosphate (UMP) pools, generates DNA damage, and induces apoptosis through suppression of replication forks-an "on-target" effect, as shown by uridine rescue. Matrix-assisted laser desorption/ionization (MALDI) mass spectroscopy imaging demonstrates that this DHODH inhibitor (BAY2402234) accumulates in the brain at therapeutically relevant concentrations, suppresses de novo pyrimidine biosynthesis in vivo, and prolongs survival of mice bearing intracranial DMG xenografts, highlighting BAY2402234 as a promising therapy against DMGs.

An affinity for brainstem microglia in pediatric high-grade gliomas of brainstem origin.

Background: High-grade gliomas (HGG) in children have a devastating prognosis and occur in a remarkable spatiotemporal pattern. Diffuse midline gliomas (DMG), including diffuse intrinsic pontine gliomas (DIPG), typically occur in mid-childhood, while cortical HGGs are more frequent in older children and adults. The mechanisms behind this pattern are not clear. Methods: We used mouse organotypic slice cultures and glial cell cultures to test the impact of the microenvironment on human DIPG cells. Comparing the expression between brainstem and cortical microglia identified differentially expressed secreted proteins. The impact of some of these proteins on DIPGs was tested. Results: DIPGs, pediatric HGGs of brainstem origin, survive and divide more in organotypic slice cultures originating in the brainstem as compared to the cortex. Moreover, brainstem microglia are better able to support tumors of brainstem origin. A comparison between the two microglial populations revealed differentially expressed genes. One such gene, interleukin-33 (IL33), is highly expressed in the pons of young mice and its DIPG receptor is upregulated in this context. Consistent with this observation, the expression levels of IL33 and its receptor, IL1RL1, are higher in DIPG biopsies compared to low-grade cortical gliomas. Furthermore, IL33 can enhance proliferation and clonability of HGGs of brainstem origin, while blocking IL33 in brainstem organotypic slice cultures reduced the proliferation of these tumor cells. Conclusions: Crosstalk between DIPGs and the brainstem microenvironment, in particular microglia, through IL33 and other secreted factors, modulates spatiotemporal patterning of this HGG and could prove to be an important future therapeutic target.

VRK1 as a synthetic lethal target in VRK2 promoter-methylated cancers of the nervous system.

Collateral lethality occurs when loss of a gene/protein renders cancer cells dependent on its remaining paralog. Combining genome-scale CRISPR/Cas9 loss-of-function screens with RNA sequencing in over 900 cancer cell lines, we found that cancers of nervous system lineage, including adult and pediatric gliomas and neuroblastomas, required the nuclear kinase vaccinia-related kinase 1 (VRK1) for their survival in vivo. VRK1 dependency was inversely correlated with expression of its paralog VRK2. VRK2 knockout sensitized cells to VRK1 loss, and conversely, VRK2 overexpression increased cell fitness in the setting of VRK1 loss. DNA methylation of the VRK2 promoter was associated with low VRK2 expression in human neuroblastomas and adult and pediatric gliomas. Mechanistically, depletion of VRK1 reduced barrier-to-autointegration factor phosphorylation during mitosis, resulting in DNA damage and apoptosis. Together, these studies identify VRK1 as a synthetic lethal target in VRK2 promoter-methylated adult and pediatric gliomas and neuroblastomas.

Single-cell epigenetic analysis reveals principles of chromatin states in H3.3-K27M gliomas.

Cancer cells are highly heterogeneous at the transcriptional level and epigenetic state. Methods to study epigenetic heterogeneity are limited in throughput and information obtained per cell. Here, we adapted cytometry by time-of-flight (CyTOF) to analyze a wide panel of histone modifications in primary tumor-derived lines of diffused intrinsic pontine glioma (DIPG). DIPG is a lethal glioma, driven by a histone H3 lysine 27 mutation (H3-K27M). We identified two epigenetically distinct subpopulations in DIPG, reflecting inherent heterogeneity in expression of the mutant histone. These two subpopulations are robust across tumor lines derived from different patients and show differential proliferation capacity and expression of stem cell and differentiation markers. Moreover, we demonstrate the use of these high-dimensional data to elucidate potential interactions between histone modifications and epigenetic alterations during the cell cycle. Our work establishes new concepts for the analysis of epigenetic heterogeneity in cancer that could be applied to diverse biological systems.