Skip to main content

Global activation of oncogenic pathways underlies therapy resistance in diffuse midline glioma

Abstract

Diffuse midline gliomas (DMGs) are aggressive pediatric brain tumors with dismal prognosis due to therapy-resistant tumor growth and invasion. We performed the first integrated histologic/genomic/proteomic analysis of 21 foci from three pontine DMG cases with supratentorial dissemination. Histone H3.3-K27M was the driver mutation, usually at high variant allele fraction due to recurrent chromosome 1q copy number gain, in combination with germline variants in ATM, FANCM and MYCN genes. Both previously reported and novel recurrent copy number variations and somatic pathogenic mutations in chromatin remodeling, DNA damage response and PI3K/MAPK growth pathways were variably detected, either in multiple or isolated foci. Proteomic analysis showed global upregulation of histone H3, lack of H3-K27 trimethylation, and further impairment of polycomb repressive complex 2 by ASXL1 downregulation. Activation of oncogenic pathways resulted from combined upregulation of N-MYC, SOX2, p65/p50 NF-κB and STAT3 transcription factors, EGFR, FGFR2, PDGFRα/β receptor tyrosine kinases, and downregulation of PHLPP1/2, PTEN and p16/INK4A tumor suppressors. Upregulation of SMAD4, PAI-1, CD44, and c-SRC in multiple foci most likely contributed to invasiveness. This integrated comprehensive analysis revealed a complex spatiotemporal evolution in diffuse intrisic pontine glioma, recommending pontine and cerebellar biopsies for accurate populational genetic characterization, and delineated common signaling pathways and potential therapeutic targets. It also revealed an unsuspected activation of a multitude of oncogenic pathways, including cancer cell reprogramming, explaining the resistance of DMG to current therapies.

Introduction

Diffuse midline glioma with histone H3-K27M mutation (DMG/K27M) is a newly defined entity in the World Health Organization (WHO) group of grade IV diffuse gliomas, usually seen in the pediatric population [1]. Most of these tumors occur in the pons, as diffuse intrinsic pontine gliomas (DIPG), and a smaller proportion occur in the thalamus, cerebellum or spinal cord [2]. The prognosis is dismal, with a median patient survival of 1 year, very similar to glioblastoma. In contrast to glioblastoma, even subtotal surgical resection is not possible due to the tumor location in vital brain areas. The treatment therefore relies on radiotherapy, usually in combination with chemotherapy and/or novel immunotherapies [3,4,5,6]. After initial therapy, a relatively brief period of tumor regression may be seen, after which the remaining tumor cells regrow and invade at distance from the primary focus, establishing secondary foci resistant to any therapy.

The discovery of histone H3-K27M mutations in a large subset of pediatric DMGs [7, 8] was underscored by a significantly poorer survival than in the non-mutated subset [9], although this survival difference has not been confirmed in a recent large study [10]. The p.K27M mutation occurs predominantly in H3F3A, encoding histone H3.3, and rarely in HIST1H3B/C genes, encoding H3.1 and usually affecting younger patients with better prognosis [11, 12]. The genetic landscaping has revealed mutations in other genes, such as TP53, ATRX and ACVR1, the last usually occurring in conjunction with H3.1 mutations in approximately 20% of DIPG cases [11]. Receptor tyrosine kinase (RTK) gene gain/amplification and pathogenic mutations were also found, especially of PDGFRA and MET [9, 13,14,15,16]. Most of these studies relied on pontine biopsies, and only limited studies addressed the spatiotemporal evolution of DIPG, mainly within infratentorial foci [13, 17].

We present here the spatiotemporal evolution of DIPG with H3F3A p.K27M driver mutation in three patients with complex infratentorial and supratentorial involvement. The genetic landscaping delineated known and novel recurrent chromosomal alterations and gene mutations. Importantly, all three patients harbored germline mutations in genes previously shown to increase the risk for cancer development. The proteomic analysis, a field previously unexplored in DIPG, interrogated oncogenic and cell invasion pathways. Our integrated results delineated novel migration patterns during DIPG progression and showed activation of many oncogenic pathways, deriving new considerations for candidate therapeutic targets.

Material and methods

Autopsy, histology and tumor burden quantification

The autopsies were performed as previously described [18], in accordance to hospital regulations. The parents consented the patients’ autopsies for diagnosis and research. The nursing team ensured prompt post-mortem transportation for autopsy. A new standardized sampling and tumor burden quantification protocol was applied. The latest MRI studies were available for all patients and used to guide autopsy sampling. The fresh brain was weighed and compared to normal standards [19]. Fresh and pilot overnight formalin-fixed paraffin-embedded (FFPE) sections were harvested from primary and secondary tumor foci, as well as from apparently normal brain and pituitary gland. The brains were subsequently fixed in 20% formalin for 10–12 days, and extensively sampled to include the following structures: whole multiple levels of the cervical spinal cord and medulla, and complete axial sections of the pons and midbrain. Bilateral samples were also harvested for paired structures: cerebellar hemispheres, thalami, basal ganglia and hippocampi. Additional samples included frontal and occipital cortex, anterior body and splenium of the corpus callosum, and fornix. Care was taken to include sections from the walls of all of the ventricles. FFPE sections were stained with hematoxylin-eosin (H&E). For patient F10, the initial biopsy was available for comparison. Images were acquired at various magnifications with Aperio Scanscope CS2 whole slide image system (Leica Biosystems, San Diego, CA) or with Nikon Eclipse Ci microscope equipped with Nikon Digital Sight DS-Fi2 camera (Nikon Instruments Inc., Melville, NY), as previously described [20, 21]. The histologic tumor burden was quantified on a 0-to-4 scale, with 0 representing no involvement, 1 – invasion by single cells, 2 – moderate cellularity, 3 – high cellularity, and 4 – massive involvement with either tumor necrosis or microvascular proliferation. Numerical data were represented graphically by using GraphPad Prism (Version 8.3.0, GraphPad Software, La Jolla, CA).

Immunohistochemistry (IHC)

IHC was performed on selected sections, as described [20, 22]. Primary antibodies were: histone H3-K27M (Millipore/Sigma, Burlington, MA), IDH1-R132H (DIA-H09, Dianova, Hamburg, Germany), p53 (DO-7), vimentin (V9), Ki-67 (30–9) (Roche/Ventana, Tucson, AZ), Olig-2 (387 M-15), GFAP (EP672Y) (Ventana/Cell Marque, Rocklin, CA).

Transmission electron microscopy

Freshly collected autopsy samples were processed as previously described [20, 23].

Next generation sequencing (NGS)

Nucleic acids were extracted from fresh frozen or FFPE samples, as described [22]. FFPE microdissection was performed on selected samples either to separate morphologically different neoplastic populations or tumor cells from normal cells. Matched normal tissue was obtained either from pituitary gland (F5 and F10), or uninvolved frontal cerebral cortex (F12). All samples were sequenced by using a customized 295-gene DNA library (SureSelect XT-HS, Agilent, Santa Clara, CA), as described [22, 23]. For all F5 and F12 samples, and for F10 ganglion-like cells (GCs), NGS was also performed at Tempus Labs (Chicago, IL), by using the xT 596-gene panel, as described [22, 24]. The two F10 biopsies were additionally sequenced at University of California/San Francisco. Variant analysis and interpretation were performed as described [22,23,24]. Copy number variation (CNV) analysis was performed at Tempus Labs [25]. Loss of heterozygosity (LOH) refers to alterations with loss of one allele. Tumor mutation burden (TMB) represents the number of single nucleotide protein-altering mutations per million base pairs.

Proteomic analysis

Fresh frozen tissue lysis and Western blotting (WB) were performed as described [26,27,28]. Briefly, 25–50 mg frozen tissue was homogenized on ice with 500 μl Tissue Protein Extraction Reagent (T-PER) reagent (Pierce Biotechnology, Rockford, IL) containing Halt Protease Inhibitor Single-use Cocktail, EDTA-Free (Pierce Biotechnology). Following centrifugation at 10,000×g for 5 min, the supernatant was collected, protein concentration was measured and the samples were denatured and resolved by SDS-polyacrylamide gel electrophoresis. Multiple primary antibodies were tested for the same protein, and the ones showing the highest specificity were selected for WB analysis (Suppl. Table S1). Positive controls for specificity testing included previously tested glioblastoma cell lines [28], and patient samples with known mutations. The WBs for each antibody were repeated at least twice, with similar results. The densitometric analysis was performed by scanning the X-ray films with optimal exposures on a ChemiDoc™ Touch imager (Bio-Rad, Hercules, CA). The bands were further quantified by using the Image Lab 6.0 software (Bio-Rad). Individual protein values were normalized to the corresponding actin values, except for phosphoprotein values that were normalized to the corresponding unphosphorylated protein values. Minus values were manually adjusted as zero. Results were expressed as fold-increase or fold-decrease compared to normal control.

Results

Clinical overview of three DIPG patients

In order to understand disease progression and resistance to therapy in DMG/K27M, we performed an integrative comparative analysis of three autopsies of 5-, 10- and 12-year-old White females (F5, F10, F12) with normal psychomotor development, normal weight, no co-morbidities, and residing within the same geographic area (Table 1).

Table 1 Patient clinical data

F5 initially complained of daily headaches and recurrent ear and sinus infections, followed by behavioral changes manifested by extreme anger, and right eye strabismus. Magnetic resonance imaging (MRI) showed a 5.5 × 5.0 × 3.9 cm, non-enhancing, pontine mass with secondary hydrocephalus of 4th and 3rd ventricles that was radiologically diagnosed as DIPG without biopsy. She initially received 6 weeks of radiotherapy and veliparib (PARP1/2 inhibitor), followed over several months by 8 carboplatin and HDAC inhibitor infusions by convection enhanced delivery with significant tumor regression. Fourteen months post-diagnosis, the patient presented with dysphagia, and MRI showed tumor regrowth and extension to right cerebellum (Suppl. Fig. S1A). F5 expired 20 months post-diagnosis and the definitive diagnosis of DMG/K27M, WHO grade IV, was rendered at autopsy.

F10 presented with chronic headaches, hearing loss and subsequent left facial droop. MRI showed a 5.3 × 4.4 × 3.9 cm, rim-enhancing, pontocerebellar mass, with 4th ventricle compression, tonsillar herniation and 3rd and lateral ventricles dilatation (Suppl. Fig. S1B). A biopsy rendered the diagnosis DMG/K27M, WHO grade IV. The patient received 6 weeks of initial brainstem radiotherapy, followed by K27M-peptide vaccine therapy (7 cycles; 3 weeks/month) until week 24, when MRI showed progression by involvement of supratentorial periventricular brain structures (Suppl. Fig. S1C). A repeat frontal corpus callosum biopsy was performed at 10 months post-diagnosis, confirming secondary tumor involvement. F10 further received frontal area proton beam therapy (3 cycles) and nivolumab (3 cycles). At this time, she presented chronic left 6th and 7th cranial nerve palsy, and abnormalities of coordination and gait. She further received CUDC907, a dual HDAC and PI3K-p110α inhibitor (1 cycle), but her status declined abruptly and she expired 13.5 months post-diagnosis.

F12 presented with headache, left-sided facial weakness and balance disturbance. MRI showed a large pontine mass that was biopsied and diagnosed as glioblastoma, WHO grade IV. She received 6 weeks of radiotherapy and concurrent erlotinib and bevacizumab chemotherapy. An MRI performed 11.5 months post-diagnosis showed tumor progression with bilateral cerebellum and right frontal periventricular involvement (Suppl. Fig. S1D). She expired 13 months post-diagnosis, and was retrospectively diagnosed at autopsy with DMG/K27M, WHO grade IV.

Patterns of histologic invasion in DIPG

Brain gross examination showed diffuse enlargement of the pons in all cases, and regional foliar expansion and effacement in the right cerebellar hemisphere (Fig. 1a). The F10 fresh brain had increased weight (Suppl. Table S2), with mild gyri flattening, consistent with cerebral edema. All cases presented with shunted hydrocephalus and F5 had right cerebellar tonsil herniation (Suppl. Table S2). Brain sectioning showed liquefying necrosis in all tumor foci noted radiologically (Fig. 1b, black arrows), except for F10 pons and F12 right frontal focus that appeared fleshy, red, and indurated (Fig. 1b, blue arrows).

Fig. 1
figure1

Histologic patterns of tumor cell invasion in DIPG. a. Brain gross appearance. Red and green arrows indicate the primary pontine and the secondary cerebellar foci, respectively. R, right; L, left. b. Sections of indicated brain structures with gross visible tumor foci. Black and blue arrows indicate liquefying necrosis and fleshy tumor, respectively. c. Semiquantitative analysis of histological tumor involvement of the specified color-coded structures. MVP, microvascular proliferation; 4 V, 4th ventricle; PF, posterior fossa; CN, cranial nerve; 3 V, 3rd ventricle; CC, corpus callosum; Temp, temporal; Occip, occipital; WM, white matter; Pyr, pyramidal tracts; IC, internal capsule. The panels and graphs from (b) and (c), respectively, correspond to the brains from (a)

Tumor burden digital quantification (Fig. 1c; Suppl. Fig. S2) delineated distinct cell dispersal patterns: (1) centrifugal parenchymal infiltration from the primary pontine focus in F5, (2) massive seeding of ventricles and leptomeninges with minimal neuropil infiltration in F10, and a mixed pattern in F12 that showed the highest tumor burden.

At least two morphologically distinct populations were apparent in each case (Fig. 2; Suppl. Fig. S3). F5 and F10 showed areas of small neoplastic cells, with elongated or irregular nuclei and scant cytoplasm (Fig. 2, F5: panels 1–2, 4–6; F10: panels 10, 13–14), and areas of larger neoplastic cells, with abundant eosinophilic cytoplasm, resembling reactive astrocytes (Fig. 2, F5: panel 3; F10: panels 15, 17–18, Suppl. Fig. S3B). The initial F10 pontine biopsy showed similar morphology to the postmortem pontine sections except for the presence of a myxoid background (Fig. 2, compare panels 7 to 10). A third unexpected morphology of large, dysplastic ganglion-like cells (GCs), as in ganglioglioma, previously reported in DMG/K27M [29], with H3-K27M expression but lacking BRAF mutations, was present in F10 pons focally (Fig. 2, panels 11–12, Suppl. Fig. S4), demonstrating the metaplastic potential of these neoplastic cells. An F10 incidental autopsy finding was the presence of striated muscle fibers surviving within the tumor at the site of the previous biopsy, most likely carried by the surgical procedure (Fig. 2, panel 9). F12 showed predominantly glioblastoma-like pleomorphic cells (Fig. 2, panels 19, 22, Suppl. Fig. S3C), and focally in pons, elongated bipolar neoplastic cells (Fig. 2, panel 20, Suppl. Fig. S3D). Common features for all tumors were pseudopalisading necrosis and microvascular proliferation. Leptomeningeal invasion was developed in F10 and F12 (Fig. 2, panels 13&24, red arrows). Mitotic activity was brisk in F5 sections (Fig. 2, panels 1, 5, arrowheads), but scant in F10 and F12 postmortem sections.

Fig. 2
figure2

Morphologic spectrum of the 3 DIPGs across infratentorial and supratentorial foci. H&E and IHC with Ki-67, histone H3-K7M, PTEN and p53 antibodies of selected sections indicated by numbers and anatomic locations. The columns correspond to each case, with two central columns shown for F10. L, left; R, right; Pyr decuss, pyramidal tract decussation; LGB, lateral geniculate body; GP, globus pallidus; periV, periventricular; CC, corpus callosum. Blue arrowheads and black, red, green and yellow arrows indicate mitotic figures, microvascular proliferation, leptomeningeal spread, subventricular ependymal cell proliferation, and foamy macrophages, respectively. Among other features, note various neoplastic cell morphologies, supratentorial involvement for all tumors, brisk mitotic activity in multiple F5 foci and in the F10 biopsy, viable striated muscle fibers in F10, most likely at prior biopsy site (panel 9), labeling with H3-K27M antibody of F10 bi-nucleated ganglion cells (blue contour in panel 12), and PTEN expression detected in endothelial cells but not in tumor cells, in F12

All neoplastic cells expressed H3-K27M mutant by IHC (Fig. 2; Suppl. Fig. S5, Suppl. Table S3). IDH1-R132H was negative in all cases, p53 expression was high in F12 (Fig. 2, panel 24), but not in F5 or F10 (Suppl. Fig. S5), and PTEN expression loss was noted in F12 neoplastic cells (Fig. 2, panel 23). Olig-2 was expressed in F5 but not in F10 (Suppl. Table S3). GFAP labeled F5 and F10, but not F12 neoplastic cells, whereas vimentin was mostly negative in F5 and F12 neoplastic cells, and positive in the F10 small neoplastic cells (Suppl. Fig. S5).

NGS revealed germline and progressive somatic alterations in DIPG

NGS was performed in the primary pontine and secondary infratentorial and supratentorial foci, two F10 biopsies, matching normal tissue, totaling 21 tumor and 3 normal samples. All the patients had confirmed or potentially deleterious germline mutations in cancer-associated genes (Fig. 3; Suppl. Table S4). F5 harbored a heterozygous germline MYCN in-frame deletion of three residues p.264_DDE_266, mapping to the conserved central acidic protein region that is a CKII phosphorylation target [30]. Hypothetically, this deletion may generate a consensus p53-binding motif that may lead to p53 inactivation similarly to viral oncogenes. Interestingly, MYCN normal allele CN gain occurred in the right cerebellar focus (Fig. 3; Suppl. Table S5). F10 presented heterozygous germline ATM splice-site mutation predicted to truncate the carboxyl (C)-terminus and classified as germline pathogenic mutation for ataxia-telangiectasia (AT) syndrome (ClinVar, multiple submitters). This germline mutation was accompanied by normal allele somatic inactivation through various mechanisms (see below), consistent with a tumor suppressor role for ATM. F12 harbored a heterozygous germline FANCM nonsense mutation, p.Q1701*, predicted to truncate FANCM C-terminus, and shown to confer increased risk for breast cancer development [31]. Paradoxically, this mutant allele was lost in all the tumor foci, by chromosome 14q CN loss (Fig. 3; Suppl. Table S5).

Fig. 3
figure3

Spatio-temporal genomic profiling. a. Pathway color-coded representation of gene mutations and selected CNVs: chromatin remodeling – pink; DDR – purple; transcription factors and modulators – green; growth factor receptors and mediators – yellow; other – blue. White boxes lack mutations, and numbers in boxes indicate the number of mutant variants for the respective gene. FFPE microdissected tumor foci are showed in blue font and show decreased TBM for the F5 samples, indicating limited tumor material available for NGS that most likely lowered the mutation and CNV detection rate in these samples. b. Color-coded CNV analysis. c. VAF graphic representation: germline mutations as dots, somatic mutations as bars. FFPE-microdissected samples are shown in blue font and show decreased VAF for the F5 samples

H3F3A p.K27M was the only somatic aberration present in all tumor foci from all cases (Fig. 3a). Its variant allele fraction (VAF) was higher than that of other heterozygous gene alterations in F5 and F10, most likely due to CN gain of the mutant H3F3A allele by a recurrent 1q gain (Fig. 3b-c; Suppl. Table S5). Somatic mutations in three other chromatin remodeling genes, ASXL1, BCOR and KDM6A, were detected in isolated foci from F10 and F12 (Fig. 3a). Although ASXL1 mutations are usually seen in non-glial neoplasms [32], rare mutations have been described in DMG/K27M [33]. Interestingly, the pathogenic protein-truncating ASXL1 mutation was solely detected in the F10 pontine biopsy and pontine autopsy GC population (Fig. 3a), suggesting common derivation. BCOR, is a component of the noncanonical polycomb repressive complex (PRC) 1.1 [34], and BCOR mutations have been described focally in DIPG with H3.3-K27M mutation [13]. In F12, a BCOR frameshift loss of function mutation was detected in the pontine focus and also in the frontal focus at subclonal VAF. KDM6A is a di- and trimethyl H3-K27 demethylase that modulates the recruitment of PRC1 and the monoubiquitination of histone H2A [35]; KDM6A alterations have not been previously described in DIPG. An inactivating KDM6A splice variant was detected at subclonal VAF in F12 frontal focus (Fig. 3; Suppl. Table S4).

Somatic pathogenic mutations in the DNA damage response (DDR) pathway were present in all three cases (Fig. 3; Suppl. Table S4). In F12, the TP53 p.R248W pathogenic mutation with or without LOH was detected in all tumor foci. In F5, as many as five clonal or subclonal TP53 mutations were detected in pontine and right cerebellar foci. In F10, a gain-of function truncating mutation in exon 6 of PPM1D that encodes the wild-type p53-induced phosphatase 1D (WIP1) was detected at clonal VAF in the cerebellar and supratentorial foci. PPM1D/WIP1 is a phosphatase that dephosphorylates and inactivates many DDR mediators, such as ATM, p53, CHK2 and H2AX, and truncating exon 6 mutations were previously reported in DIPG as mutually exclusive with TP53 mutations [36,37,38]. In addition, F10 harbored second-hit inactivation of the normal ATM allele by different mechanisms in different foci: (1) LOH by chromosome 11q loss containing the ATM locus, for the pontine biopsy and autopsy GC populations, (2) p.L2877F somatic mutation at clonal VAF, for the cerebellar and supratentorial foci, and (3) p.G2342S somatic mutation, for the left pontine population. The L2877F somatic mutation was previously reported as pathogenic mainly in lung adenocarcinomas [39], and a germline G2342S variant of uncertain significance (VUS) was reported in the AT syndrome (ClinVar/Invitae).

Somatic mutations in transcription factors or modulators involved in cell differentiation and cell cycle were detected in F5 and F12 (Fig. 3a). Most of the alterations were focal, due to homozygous CN loss in DAXX, BTG1, CUX1, RUNX1 (Suppl. Table S5), MYCN amplification, CAMTA1 frameshift or TCF7L2 missense mutation (Suppl. Table S4). A somatic TCF12 p.L352* nonsense mutation present in all F12 tumor foci at heterozygous VAF suggested early co-occurrence with the H3F3A p.K27M mutation, as well as selective advantage (Fig. 3c). This mutation truncates half of the protein, including the C-terminal DNA-binding bHLH domain that is required for transcriptional activity [40]. Although TCF12 mutations were not previously reported in DIPG, they were shown in anaplastic oligodendroglioma to correlate with more aggressive tumor types [40].

Somatic alterations of the RTK/PI3K/MAPK/mTOR pathways were also seen in all cases, although very focal in F10 (Fig. 3; Suppl. Tables S4-S5). PDGFRA and KIT amplification was detected in F12 cerebellar and frontal foci. An oncogenic PDGFRA p.S564L missense mutation mapping to the juxtamembrane region, previously described in colorectal carcinoma [41], was present at low VAF in F12 cerebellar tumor. PIK3CA p.E545G and p.V344G pathogenic mutations mapping to hot spots in the C2 and helical domains of PI3K-p110α catalytic subunit, respectively, were detected in almost all F5 and F12 foci, respectively. A PIK3CA p.E542K mutation was detected at very low VAF in the F10 GC population. Pathogenic PIK3R1 p.E558_M563 in frame deletion and p.E451_Y452 in-frame duplication mapping to the inter-SH2 region of PI3K-p85 regulatory subunit were detected in a left-sided F5 infratentorial subpopulation and in F10 frontal biopsy, respectively. Although not identical, similar oncogenic mutations have been reported in endometrial cancer [39, 42]. A KRAS gain-of-function mutation was detected focally in an F5 midbrain subpopulation. An NF1 R1590W somatic missense mutation mapping to the Sec14-PH domain and previously described in syndromic NF1 [43] was detected focally in F10 pons. A TSC1 pathogenic frameshift mutation was detected throughout F5, suggestive of early occurrence and selective advantage. An STK11 homozygous CN loss was noted in F12 pons.

Focal somatic pathogenic alterations occurred in B2M and GRIN2A genes encoding transmembrane proteins (Fig. 3a), the latter gene being mutated in metastatic glioblastoma [24] and frequently, in melanoma [44]. Other VUS were detected in F12: EGFR p.A613T at low VAF in frontal focus, and STAT4 p.E194* at heterozygous VAF in pontine and frontal foci (Suppl. Table S4).

The TMB in various tumor foci from all three cases showed little variation, with values between 1.3 and 1.7 mutations/MB (Fig. 3a), ranking at the lower end for various cancer types, as previously described [45].

The CNV analysis showed a combination of CN gains and losses in all cases, consisting of recurrent, case-specific and isolated CNVs (Fig. 3b; color-coded Suppl. Table S5). Recurrent CNVs, called if detected in at least two cases, were either (1) gains of 1p-distal-end; 1q, including the H3F3A locus; 4q12, mapping to PDGFRA and KIT loci; 7q, containing SERPINE1 that encodes plasminogen activator inhibitor 1 (PAI-1); 8q-distal-end; focal 14q11.2, encompassing TRA (T cell receptor alpha) locus; 17p13, including CRK proto-oncogene involved in glioblastoma invasiveness [46, 47]; 17q25; 21q11-q21, or (2) losses of 10q-distal-end, encompassing ECHS1 encoding enoyl-CoA hydratase-1, whose deficiency has been shown to activate mTOR [48]; 16q, including PHLPP2 tumor suppressor. Entire chromosome 14 loss was detected in F10 pontine biopsy and GC samples, and 14q loss, in all F12 samples. Focal 10q25-q26 CN gain was detected in three F5 samples and the pontine F12 sample, and contained high or low gain depending on FGFR2 presence or absence, respectively. Many of these recurrent CNVs are novel, compared to a previous report describing CNVs in DIPG [49]. Case-specific CNVs were noted in some or all samples from a given case, e.g. F5 and F12 contained as many as 7 or 6 different case-specific CNVs, respectively, each detected in at least two samples. F12 showed the most extensive and complex case-specific CNVs, including complex alterations on chromosome 7 in all samples, with additional LOH in the cerebellar and frontal foci, and loss of chromosome 18 in all samples, with focal high CN gain of distal 18p, containing TYMS (thymidylate synthetase), in the cerebellar and frontal foci. The 17p13 chromosomal region containing CRK and TP53, showed low CN gain in F5 and F10, sometimes without involving TP53, and an interesting evolution pattern in F12: neutral LOH in pons, CRK neutral LOH and TP53 heterozygous CN loss in cerebellum, and entire region heterozygous CN loss in the frontal focus, suggesting a role for CRK in DIPG invasion. Except for the F5 midbrain sample, all individual samples contained isolated CNVs; e.g. chromosome 10 loss, a recurrent CNV in glioblastoma, was present only in F12 cerebellar sample.

Activation of epigenetic, DDR, cell growth and migration pathways in DIPG

To profile pathway activation in DIPG, a semiquantitative proteomic analysis was undertaken in grossly apparent tumor foci and control normal cortex (Fig. 4; Suppl. Fig. S6). Global epigenetic changes included strong histone H3 upregulation, absent H3-K27 trimethylation, undetectable EZH2 expression, except for the F12 frontal focus that showed overexpression, and ASXL1 downregulated levels, without mutations in these samples. Focal H3-K27 acetylation was detected in F5 and F12, but not F10.

Fig. 4
figure4

Proteomic profiling of DIPG foci. a-b. WB analysis with indicated antibodies of total protein lysates (50 μg proteins) from tumor (T) and normal (N) fresh frozen autopsy samples. Primary pontine (P) foci are indicated in red font and secondary cerebellar (C) and frontal (Fr) foci, in green font. Sample laterality: R, right; L, left. Supplemental Fig. S6 shows the positive control for H3 K27Me3 antibody from an adult glioblastoma autopsy case with EZH2 overexpression. Note lack of H3 K27 methylation in the presence of EZH2 overexpression for F12 right frontal (RFr) tumor focus. c. Heat map of semiquantitative WB analysis, as shown quantified and normalized to loading controls in Supplemental Fig. S6. Shades of red and blue indicate expression level increase or decrease, respectively, as compared to normal controls

DDR involves an early step of histone variant H2A.X phosphorylation by ATM [50]; phospho-Ser139 H2A.X was minimal in F5, but was relatively high in F12 foci. Surprisingly, it was also detected with lower intensity in all F10 tumor foci with bi-allelic ATM inactivation, suggesting compensation by a related kinase. PARP1 levels were increased in F5 foci (F5 received an initial course with PARP1/2 inhibitors), strongly increased in all F12 foci, and minimally increased in some F10 foci. Elevated p53 levels were present in F5 and F12 tumor foci that harbored TP53 missense mutations known to stabilize the mutant protein [51]. N-MYC transcription factor has been implicated in pathogenesis of high-grade pediatric gliomas, including DIPG without H3-K27M mutation [11]. It showed high levels in F5 right pons and cerebellum, correlating with the detected MYCN amplification in the latter focus, but not in the normal sample, suggesting that the germline MYCN alteration does not influence protein expression. N-MYC was also strongly upregulated in all F12 foci and in F10 right cerebellar tumor that did not harbor gene amplification, suggesting increased transcription, protein stability, or both. SOX2, a transcription factor essential for maintaining undifferentiated neural stem cells [52], was upregulated in all three DIPGs, although patchy in F10, and with high overexpression in F12 cerebellar and frontal foci.

The levels of many oncoproteins were upregulated in DIPGs, either in generalized or more focal patterns. The heterodimeric NF-κB transcription factor proteins p50 and p65 were upregulated in all three cases in many foci. Similarly, STAT3 transcription factor was globally upregulated. Among the RTKs tested, EGFR showed upregulation in all cases, but with mild to moderate phosphorylation only in F10 and F12 foci. PDGFRα was strongly upregulated in F12 cerebellar and frontal foci, reflecting the PDGFRA amplification detected in these samples (see Fig. 3), and PDGFRβ was mildly or moderately upregulated in multiple foci. FGFR2 was mildly to moderately upregulated in foci from all three cases and correlated with CN gain in F5 right cerebellar focus (Suppl. Table S5). MET was mildly upregulated in F10 cerebellar foci.

Known mediators of glioma invasion pathways showed also increased expression levels. SMAD4, the common transducer of TGF-β and bone morphogenic protein pathways, and its downstream target PAI-1 [53], both involved in glioma invasion and angiogenesis [54, 55], were upregulated in nearly all foci. In contrast to PAI-1 upregulation that may be partly due to recurrent 7q CN gain, SMAD4 showed CN loss in all F12 foci due to chromosome 18 loss (Suppl. Table S5). Its increase is most likely controlled by epigenetic transcriptional regulation and protein turnover mechanisms. CD44, a hyaluronic acid receptor that associates and promotes c-SRC activation [56] and is involved in glioma cell invasion [57, 58], was upregulated in almost all foci and correlated with c-SRC activation. In contrast, FAK showed focal strong activation in the F12 cerebellar focus.

As expected from the overall upregulation of RTKs, the ERK/MAPK and PI3K/AKT pathways were activated in virtually all foci. The phosphorylated ERK1–2 species in F12 pontine and cerebellar foci are not shown or quantified, but were present as strong bands shifted upward, most likely due to tissue degradation, and, as expected, did not react with other MAPK family members, including JNK (not shown). The expression levels of selected tumor suppressors showed global or more focal decrease. PHLPP1/2 proteins that are PI3K/AKT [27, 59] and NF-κB pathway inhibitors [26], were globally downregulated, and their respective genes showed monoallelic loss by case-specific or recurrent chromosome 18 and 16q loss, respectively (Suppl. Table S5). PTEN tumor suppressor was downregulated in most F5, F10 and F12 foci. However, contamination by normal neurons containing high PTEN expression levels is hard to avoid and cannot be excluded in cerebellar and pontine samples. The two cyclin D-CDK4/6 complex cell cycle inhibitors showed different expression trends: p16 INK4A, encoded by CDKN2A, was downregulated in all foci, except for the F12 pontine focus that showed also PTEN retention, whereas p27 Kip1 was not significantly decreased except for the F10 frontal focus. Conversely, cyclin D1 expression levels were mildly increased in almost all foci.

Discussion

DMG/K27M are therapy-resistant aggressive tumors with extremely poor survival, occurring mainly in the pediatric population [2]. Tremendous therapeutic efforts to achieve tumor growth control [3,4,5,6] are coupled to clinicopathologic efforts for prognostic subgrouping [9,10,11,12]. A better understanding of the natural evolution of the disease, especially of the unrelenting dispersal of the tumor cells, has also been attempted [13, 17], but is hampered by the lack of pathological material from autopsy. To tackle this drawback, we performed the first integrated histopathologic/genomic/proteomic analysis of a large number of matched infratentorial and supratentorial tumor foci from 3 DIPG autopsies.

Patterns of invasion in DIPG

The semiquantitative histologic analysis of tumor spread suggested the presence of a common ponto-cerebellar migration pattern (Fig. 5a, red), and two secondary patterns: centrifugal parenchymal dispersal (Fig. 5a, blue-pattern-1), and CSF spread with anterograde seeding of supratentorial foci (Fig. 5a, green-pattern-2). The latter pathways may present relatively separated or may coexist, and it appears that the CSF spread correlated with shorter survival. The genomic/proteomic analysis of infratentorial and supratentorial foci confirmed and expanded these histologic invasion patterns. We have previously shown that the tumor mutation signature shows little spatiotemporal variability, being reliable for tracking the origin of invasive populations, whereas the CNV composition shows high variability in distinct tumor foci [24]. A common mutation set characterized all tumor foci in each DIPG case, and additional genetic alterations allowed tracking the spatiotemporal selection and migration of neoplastic populations. The CNVs reflected more precisely the marked heterogeneity of the various subsets of neoplastic cells and their patterns of invasion (Fig. 5b). For F5 that exhibited histologic centrifugal invasion, the pontine population showed a TSC1 pathogenic mutation at subclonal VAF, and this mutant subpopulation was selected and subsequently migrated to all other tumor areas. However, it appears that many different recurrent or case specific CNVs occurred in this subpopulation, tracking distinct migratory subpopulations either towards the cerebellum and midbrain or cervical spinal cord (Fig. 5b, red and orange lines), or in centrifugal paths towards the spinal cord and supratentorial foci (Fig. 5b, light blue lines). Interestingly, the spinal cord and supratentorial foci shared four additional CNVs that were not detected in the brainstem or cerebellar foci (Fig. 5b, yellow circles), and the pons and supratentorial foci shared another two similar CNVs (Fig. 5b, blue circles). We speculate that these CNVs may confer a selective migration advantage to these cell populations, and listed putative invasion-promoting genes in these chromosomal regions for further confirmatory studies (Suppl. Table S5). To add to the complexity, a pattern of centrifugal migration may originate from the cerebellum and circumvent the pons, perhaps occurring later in the evolution of the tumor (Fig. 5b, green lines). For the second histologic invasion pattern, with extensive leptomeningeal and ventricular dissemination, the genetic analysis revealed many alterations common only to the cerebellar and supratentorial foci, indicating that the supratentorial seeding originates from a cerebellar rather than pontine neoplastic population (Fig. 5a, green circle). This finding challenges the conclusion from a previous study based on two autopsy cases with supratentorial foci without histological characterization stating that DIPG is a homogenous tumor, and therefore, the pontine population present in the initial biopsy is representative for the entire tumor [17]. Especially for tumors with the second invasion pattern and/or cerebellar involvement at presentation, our data showed significant spatial genetic divergence between the pontine and cerebellar or supratentorial populations. Moreover, for F10, for whom tissue was available from 2 biopsies and autopsy, neither the initial pontine biopsy nor the pontine autopsy populations were representative for the major neoplastic population that disseminated throughout the brain.

Fig. 5
figure5

Models of invasion, pathway activation and therapeutic targeting in DIPG. a. Histologic invasion pathways include the common ponto-cerebellar dissemination pathway, shown in red, and the more specific centrifugal and CSF/leptomeningeal pathways shown in blue and green, respectively. The green circle indicates the cerebellar tumor as the most common origin of supratentorial CSF seeding. b. Spatial color-coded CNV-tracking of invasive neoplastic populations within the F5 tumor that showed histologic centrifugal migration pattern. Contiguous neoplastic populations are indicated with line tracings, and isolated foci sharing the same CNV composition, with circles. Cerebellar bars indicate bilateral spread of the respective population. The corresponding color-coded CNVs are indicated in bold for gains, and regular font, for losses. c. Schematic model of epigenetic changes targeting histone H3 K27 residue. The histone H2A/H2B/H3/H4 nucleosome is represented as a blue barrel, the DNA, as a red thread, and the N-terminal histone H3 “tail”, as a black curved line. The methyl and acetyl groups are shown in orange and green, respectively. The stop sign shows blockage of enzymatic or transcriptional activity. d. Pathway activation and drug targeting model in DIPG. The oncogenic H3 K27M driver protein and mediators with dual tumor suppressor or oncogenic roles depending on context are shown in red and purple, respectively. Otherwise, the oncogenic and the tumor suppressor proteins are shown on left and right, respectively. The pathways are shown in bold font, and red and blue indicate activation or suppression, respectively. The candidate targets for therapy are underlined

Oncogenic pathway profiling and candidate therapies in DIPG

Histone H3.3-K27M mutation was ubiquitous, consistent with its role as driver oncogenic mutation, as previously reported [17, 60]. Moreover, the mutant allele had CN gain, an observation recently reported by others [61]. A therapeutic vaccine approach against the mutant peptide [4] has been used for F10, but did not appear to extend survival. H3-K27 residue is target to multiple posttranslational modifications, being methylated by the PRC2 catalytic components EZH1/2, and acetylated by the CBP/p300 acetyltransferase complex [62] (Fig. 5c). PRC2 is recruited in the presence of ASXL1 to unmethylated CpG islands to maintain silent genes in repressed state [63, 64]. Although an isolated ASXL1 mutation was detected in F10, the proteomic analysis showed global ASXL1 expression downregulation, suggesting an additional common mechanism contributing to PRC2 inactivation in DIPG. H3-K27M mutant inhibits EZH2 catalytic activity in a dominant-negative manner [65], resulting in global H3-K27 hypomethylation and transcription de-repression [66,67,68]. The unmethylated K27 may be acetylated to further activate transcription, and we found its acetylation in F5 and F12, unrelated to prior anti-HDAC therapy. Additional chromatin remodeling genes were altered more focally, either by mutation or by CNV. TYMS acquired high CN gain in F12 cerebellar and frontal foci and may represent a 5-fluorouracil/capecitabine therapy target [69, 70].

DDR pathway mutations, either germline or somatic, were present either in the common genetic signature, in F10 and F12, or more focally, in F5 (Fig. 5d). TP53, PPM1D and ATM somatic mutations have been previously described in DIPG [11, 33, 38]. The detection of ATM germline mutation in F10 adds DMG/K27M to the list of malignancies for which AT syndrome carriers are at risk. The multiple mechanisms for second ATM allele somatic inactivation suggested independent spatial evolution from an initial pauci-mutated ponto-cerebellar population. The acquisition of a mutation in a second DDR modifier, PPM1D, conferred most likely a selective advantage to the population in the cerebellum that further populated the supratentorial compartment. PPM1D/WIP1 activation is expected to result in further ATM and p53 inactivation [71] with suppression of DDR and cell cycle check points (Fig. 5d). Surprisingly, effects previously attributed to ATM loss, such as defective H2AX-Ser139 phosphorylation and EZH2 expression increase [50, 72], were not observed in our proteomic analysis, perhaps through compensation by other mechanisms [73, 74]. This observation, and the resistance of this tumor to radiotherapy, cautions the current development of ATM pharmacologic inhibitors designed to synergize with radiation-induced DNA damage in cancer cell killing [75]. PARP1 is implicated in DDR and maintenance of genome stability [76] and PARP inhibitors confer synthetic lethality in tumors with BRCA1/BRCA2 alterations and DNA repair deficiency [77]. F5 was initially treated with PARP inhibitors without response, but did not show recurrent DDR pathway alterations as F10, for whom other therapies were tried. FANCM, a component of the Fanconi anemia complex, was considered a tumor suppressor involved in DDR, until very recent studies showed a novel role in telomere maintenance, specifically in cancer cells with alternative lengthening of telomeres (ALT) [78,79,80]. Depletion of FANCM in these cells induced reduced cell viability, and hence FANCM inhibition has been proposed as therapeutic strategy specifically in cancer cells with ALT [78,79,80]. In F12, the FANCM p.Q1701* germline mutation, deleting 2 C-terminal domains that may be involved in the ALT phenotype [78], is accompanied by paradoxical mutant allele loss in the tumor, suggesting perhaps a requirement for one functional FANCM allele. Although the FANCM p.Q1701* mutation is well established as genetic risk factor for breast cancer development [31], there are no studies addressing the FANCM somatic events in patient tumors. Therefore, the possibility of distinct roles of FANCM in normal and ALT-exhibiting cancer cells, as well as the significance of the p.Q1701* and other truncating FANCM mutations in cancer remain to be elucidated.

The detection of germline mutations in cancer-related genes in all three patients was unexpected, knowing that, despite abundant literature documenting somatic mutations, germline mutations in DMG/K27M have not been similarly explored. A pathogenic germline mutation in MUTYH, involved in Lynch syndrome and encoding a DNA base excision repair enzyme, was reported in a pediatric case of spinal cord DMG/K27M [81]. These findings warrant studies on larger cohorts to assess the germline mutation incidence, which might currently be underestimated in DMG/K27M.

The RTK/PI3K/MAPK/mTOR pathway was consistently activated in all cases by proteomic analysis, although only F5 showed PIK3CA and additional TSC1 mutations in the common signature. Focal populations or subpopulations in all three cases acquired a variety of oncogenic genomic alterations during the development of the tumor, including in PDGFRA, KIT, FGFR2, KRAS, NF1, PIK3CA, PIK3R1 and STK11 genes. Except for PDGFRA amplification that correlated with PDGFRα overexpression, the proteomic analysis showed more global upregulation or activation of RTKs, including EGFR, PDGFRβ and FGFR2. RTK inhibitors are widely used in the treatment of various cancers, but response requires tumor dependency to the respective RTK signaling, and residual disease develops as a result of tumor cell reprogramming [82]. EGFR inhibitors are successfully used in solid cancers harboring EGFR mutations or overexpression [82], and were used in F12 that showed EGFR overexpression and mild activation. However, mechanisms of rapid EGFR resistance may involve transcriptional upregulation of other RTKs, such as PDGFRβ and FGFR2/3 [83, 84], and we noted high PDGFRβ levels and FGFR2 upregulation without CN gain in F12 foci, explaining therapy resistance. Signaling through RTKs and PI3K/PTEN/PHLPP pathway induces increased cell invasion [85,86,87,88], and contributions from upregulated TGFβ/PAI-1 and CD44/c-SRC pathways [54,55,56,57,58] most likely fueled invasion in DIPG. All these pathways are targetable and efforts are underway to test drug inhibitors in glioma cells [28, 54, 89]. However, due to the multitude of activated pathways, and regardless of the epigenetic or therapy-induced mechanism of upregulation of these RTKs, transcription regulators or invasiveness mediators, we believe that more general strategies, either immune or epigenetic, aiming shutdown of aberrant transcriptional programs, such as use of BRD4 inhibitors [90] (Fig. 5c-d), might prove beneficial in DIPG.

Conclusions

This study contributed several novel genetic alterations, and described a populational heterogeneity not previously appreciated in DMG. The integrated analysis of infratentorial and supratentorial foci allowed outlining migratory pathways and revealed a complex spatiotemporal evolution in DMG that recommends, when applicable, harvesting of both pontine and cerebellar biopsies for accurate populational genetic characterization. The unique proteomic exploration from this study revealed the activation of a multitude of oncogenic and invasion pathways, including positive feedback loops triggered by chemotherapy. In addition to tumor heterogeneity, this global activation that includes cancer cell reprogramming following initial miscellaneous therapeutic attempts explains most likely DMG therapy resistance. This integrated analysis also uncovered previously unknown targetable pathways, strongly supporting the identification of rational combination therapies in DIPG.

Availability of data and materials

Supporting data for this manuscript are available in the Supplemental Material and upon request to the corresponding author.

Abbreviations

ALT:

Alternative lengthening of telomeres

AT:

Ataxia telangiectasia

ATM:

Ataxia telangiectasia mutated

C-:

Carboxyl

CN:

Copy number

CNS:

Central nervous system

CSF:

Cerebrospinal fluid

CT:

Computed tomography

DDR:

DNA damage response

GCs:

Ganglion-like cells

DIPG:

Diffuse intrinsic pontine glioma

DMG:

Diffuse midline glioma

EGFR:

Epidermal growth factor receptor

ERK/MAPK:

Extracellular signal-regulated kinase/mitogen-activated protein kinase

FAK:

Focal adhesion kinase

FANCM:

Fanconi anemia complementation group M

FFPE:

Formalin-fixed paraffin-embedded

FGFR:

Fibroblast growth factor receptor

FLAIR:

Fluid attenuated inversion recovery

GFAP:

Glial fibrillary acidic protein

H3F3A gene:

H3.3 histone A/H3 histone family member 3A

H&E:

Hematoxylin eosin

IHC:

Immunohistochemistry

LOH:

Loss of heterozygosity

MRI:

Magnetic resonance imaging

N-:

Amino

NGS:

Next generation sequencing

PAI-1:

Plasminogen activator inhibitor 1, encoded by SERPINE1 gene

PARP:

Poly (ADP-Ribose) Polymerase

PDGFR:

Platelet-derived growth factor receptor

PI3K p110α/PIK3CA :

Phosphatidylinositol 3-OH kinase catalytic subunit alpha

PI3K p85/PIK3R1 :

Phosphatidylinositol 3-OH kinase regulatory subunit

PRC:

Polycomb repressive complex

RTK:

Receptor tyrosine kinase

SNP:

Single nucleotide polymorphism

TGF-β:

Transforming growth factor β

TMB:

Tumor mutation burden

VAF:

Variant allele fraction

VUS:

Variant of unknown significance

W:

Weighted (on MRI sequences T1 and T2)

WB:

Western blot

WHO:

World Health Organization

References

  1. 1.

    Louis DN, Ohgaki H, Wiestler OD, Caveneee WK (2016) WHO classification of Tumors of the central nervous system. IARC, Lyon

    Google Scholar 

  2. 2.

    Chiang JC, Ellison DW (2017) Molecular pathology of paediatric central nervous system tumours. J Pathol 241:159–172

    CAS  PubMed  Google Scholar 

  3. 3.

    Benitez-Ribas D, Cabezon R, Florez-Grau G, Molero MC, Puerta P, Guillen A et al (2018) Immune response generated with the Administration of Autologous Dendritic Cells Pulsed with an Allogenic Tumoral cell-lines lysate in patients with newly diagnosed diffuse intrinsic Pontine Glioma. Front Oncol 8:127

    PubMed  PubMed Central  Google Scholar 

  4. 4.

    Chheda ZS, Kohanbash G, Okada K, Jahan N, Sidney J, Pecoraro M et al (2018) Novel and shared neoantigen derived from histone 3 variant H3.3K27M mutation for glioma T cell therapy. J Exp Med 215:141–157

    CAS  PubMed  PubMed Central  Google Scholar 

  5. 5.

    Mathew RK, Rutka JT (2018) Diffuse intrinsic Pontine Glioma : clinical features, molecular genetics, and novel targeted therapeutics. J Korean Neurosurg Soc 61:343–351

    CAS  PubMed  PubMed Central  Google Scholar 

  6. 6.

    Williams MJ, Singleton WG, Lowis SP, Malik K, Kurian KM (2017) Therapeutic targeting of histone modifications in adult and Pediatric high-grade Glioma. Front Oncol 7:45

    PubMed  PubMed Central  Google Scholar 

  7. 7.

    Schwartzentruber J, Korshunov A, Liu XY, Jones DT, Pfaff E, Jacob K et al (2012) Driver mutations in histone H3.3 and chromatin remodelling genes in paediatric glioblastoma. Nature 482:226–231

    CAS  PubMed  Google Scholar 

  8. 8.

    Wu G, Broniscer A, McEachron TA, Lu C, Paugh BS, Becksfort J et al (2012) Somatic histone H3 alterations in pediatric diffuse intrinsic pontine gliomas and non-brainstem glioblastomas. Nat Genet 44:251–253

    CAS  PubMed  PubMed Central  Google Scholar 

  9. 9.

    Khuong-Quang DA, Buczkowicz P, Rakopoulos P, Liu XY, Fontebasso AM, Bouffet E et al (2012) K27M mutation in histone H3.3 defines clinically and biologically distinct subgroups of pediatric diffuse intrinsic pontine gliomas. Acta Neuropathol 124:439–447

    CAS  PubMed  PubMed Central  Google Scholar 

  10. 10.

    Hoffman LM, Veldhuijzen van Zanten SEM, Colditz N, Baugh J, Chaney B, Hoffmann M et al (2018) Clinical, radiologic, pathologic, and molecular characteristics of long-term survivors of diffuse intrinsic Pontine Glioma (DIPG): a collaborative report from the international and European Society for Pediatric Oncology DIPG registries. J Clin Oncol 36:1963–1972

    CAS  PubMed  PubMed Central  Google Scholar 

  11. 11.

    Buczkowicz P, Hoeman C, Rakopoulos P, Pajovic S, Letourneau L, Dzamba M et al (2014) Genomic analysis of diffuse intrinsic pontine gliomas identifies three molecular subgroups and recurrent activating ACVR1 mutations. Nat Genet 46:451–456

    CAS  PubMed  PubMed Central  Google Scholar 

  12. 12.

    Castel D, Philippe C, Calmon R, Le Dret L, Truffaux N, Boddaert N et al (2015) Histone H3F3A and HIST1H3B K27M mutations define two subgroups of diffuse intrinsic pontine gliomas with different prognosis and phenotypes. Acta Neuropathol 130:815–827

    CAS  PubMed  PubMed Central  Google Scholar 

  13. 13.

    Hoffman LM, DeWire M, Ryall S, Buczkowicz P, Leach J, Miles L et al (2016) Spatial genomic heterogeneity in diffuse intrinsic pontine and midline high-grade glioma: implications for diagnostic biopsy and targeted therapeutics. Acta Neuropathol Commun 4:1

    PubMed  PubMed Central  Google Scholar 

  14. 14.

    Koschmann C, Zamler D, MacKay A, Robinson D, Wu YM, Doherty R et al (2016) Characterizing and targeting PDGFRA alterations in pediatric high-grade glioma. Oncotarget 7:65696–65706

    PubMed  PubMed Central  Google Scholar 

  15. 15.

    Paugh BS, Broniscer A, Qu C, Miller CP, Zhang J, Tatevossian RG et al (2011) Genome-wide analyses identify recurrent amplifications of receptor tyrosine kinases and cell-cycle regulatory genes in diffuse intrinsic pontine glioma. J Clin Oncol 29:3999–4006

    CAS  PubMed  PubMed Central  Google Scholar 

  16. 16.

    Paugh BS, Zhu X, Qu C, Endersby R, Diaz AK, Zhang J et al (2013) Novel oncogenic PDGFRA mutations in pediatric high-grade gliomas. Cancer Res 73:6219–6229

    CAS  PubMed  PubMed Central  Google Scholar 

  17. 17.

    Nikbakht H, Panditharatna E, Mikael LG, Li R, Gayden T, Osmond M et al (2016) Spatial and temporal homogeneity of driver mutations in diffuse intrinsic pontine glioma. Nat Commun 7:11185

    CAS  PubMed  PubMed Central  Google Scholar 

  18. 18.

    Georgescu MM, Pinho Mda C, Richardson TE, Torrealba J, Buja LM, Milewicz DM et al (2015) The defining pathology of the new clinical and histopathologic entity ACTA2-related cerebrovascular disease. Acta Neuropathol Commun 3:81

    PubMed  PubMed Central  Google Scholar 

  19. 19.

    Dekaban AS (1978) Changes in brain weights during the span of human life: relation of brain weights to body heights and body weights. Ann Neurol 4:345–356

    CAS  PubMed  Google Scholar 

  20. 20.

    Georgescu MM, Olar A, Mobley BC, Faust PL, Raisanen JM (2018) Epithelial differentiation with microlumen formation in meningioma: diagnostic utility of NHERF1/EBP50 immunohistochemistry. Oncotarget 9:28652–28665

    PubMed  PubMed Central  Google Scholar 

  21. 21.

    Georgescu MM, Yell P, Mobley BC, Shang P, Georgescu T, Wang SH et al (2015) NHERF1/EBP50 is an organizer of polarity structures and a diagnostic marker in ependymoma. Acta Neuropathol Commun 3:11

    PubMed  PubMed Central  Google Scholar 

  22. 22.

    Georgescu MM, Li Y, Islam MZ, Notarianni C, Sun H, Olar A et al (2019) Mutations of the MAPK/TSC/mTOR pathway characterize periventricular glioblastoma with epithelioid SEGA-like morphology-morphological and therapeutic implications. Oncotarget 10:4038–4052

    PubMed  PubMed Central  Google Scholar 

  23. 23.

    Georgescu MM, Nanda A, Li Y, Mobley BC, Faust PL, Raisanen JM et al (2020) Mutation status and epithelial differentiation stratify recurrence risk in Chordoid meningioma-a Multicenter study with high prognostic relevance. Cancers (Basel) 12:225

    CAS  Google Scholar 

  24. 24.

    Georgescu MM, Olar A (2020) Genetic and histologic spatiotemporal evolution of recurrent, multifocal, multicentric and metastatic glioblastoma. Acta Neuropathol Commun 8:10

    CAS  PubMed  PubMed Central  Google Scholar 

  25. 25.

    Beaubier N, Bontrager M, Huether R, Igartua C, Lau D, Tell R et al (2019) Integrated genomic profiling expands clinical options for patients with cancer. Nat Biotechnol 37:1351–1360

    CAS  PubMed  Google Scholar 

  26. 26.

    Agarwal NK, Zhu X, Gagea M, White CL 3rd, Cote G, Georgescu MM (2014) PHLPP2 suppresses the NF-kappaB pathway by inactivating IKKbeta kinase. Oncotarget 5:815–823

    PubMed  PubMed Central  Google Scholar 

  27. 27.

    Molina JR, Agarwal NK, Morales FC, Hayashi Y, Aldape KD, Cote G et al (2012) PTEN, NHERF1 and PHLPP form a tumor suppressor network that is disabled in glioblastoma. Oncogene 31:1264–1274

    CAS  PubMed  Google Scholar 

  28. 28.

    Zhu X, Morales FC, Agarwal NK, Dogruluk T, Gagea M, Georgescu MM (2013) Moesin is a glioma progression marker that induces proliferation and Wnt/beta-catenin pathway activation via interaction with CD44. Cancer Res 73:1142–1155

    CAS  PubMed  Google Scholar 

  29. 29.

    Solomon DA, Wood MD, Tihan T, Bollen AW, Gupta N, Phillips JJ et al (2016) Diffuse midline Gliomas with histone H3-K27M mutation: a series of 47 cases assessing the Spectrum of morphologic variation and associated genetic alterations. Brain Pathol 26:569–580

    CAS  PubMed  Google Scholar 

  30. 30.

    Hagiwara T, Nakaya K, Nakamura Y, Nakajima H, Nishimura S, Taya Y (1992) Specific phosphorylation of the acidic central region of the N-myc protein by casein kinase II. Eur J Biochem 209:945–950

    CAS  PubMed  Google Scholar 

  31. 31.

    Kiiski JI, Pelttari LM, Khan S, Freysteinsdottir ES, Reynisdottir I, Hart SN et al (2014) Exome sequencing identifies FANCM as a susceptibility gene for triple-negative breast cancer. Proc Natl Acad Sci U S A 111:15172–15177

    CAS  PubMed  PubMed Central  Google Scholar 

  32. 32.

    Katoh M (2013) Functional and cancer genomics of ASXL family members. Br J Cancer 109:299–306

    CAS  PubMed  PubMed Central  Google Scholar 

  33. 33.

    Mackay A, Burford A, Carvalho D, Izquierdo E, Fazal-Salom J, Taylor KR et al (2017) Integrated molecular meta-analysis of 1,000 Pediatric high-grade and diffuse intrinsic Pontine Glioma. Cancer Cell 32:520–37.e5

    CAS  PubMed  PubMed Central  Google Scholar 

  34. 34.

    Astolfi A, Fiore M, Melchionda F, Indio V, Bertuccio SN, Pession A (2019) BCOR involvement in cancer. Epigenomics 11:835–855

    CAS  PubMed  PubMed Central  Google Scholar 

  35. 35.

    Lee MG, Villa R, Trojer P, Norman J, Yan KP, Reinberg D et al (2007) Demethylation of H3K27 regulates polycomb recruitment and H2A ubiquitination. Science 318:447–450

    CAS  PubMed  Google Scholar 

  36. 36.

    Shreeram S, Demidov ON, Hee WK, Yamaguchi H, Onishi N, Kek C et al (2006) Wip1 phosphatase modulates ATM-dependent signaling pathways. Mol Cell 23:757–764

    CAS  PubMed  Google Scholar 

  37. 37.

    Wu G, Diaz AK, Paugh BS, Rankin SL, Ju B, Li Y et al (2014) The genomic landscape of diffuse intrinsic pontine glioma and pediatric non-brainstem high-grade glioma. Nat Genet 46:444–450

    CAS  PubMed  PubMed Central  Google Scholar 

  38. 38.

    Zhang L, Chen LH, Wan H, Yang R, Wang Z, Feng J et al (2014) Exome sequencing identifies somatic gain-of-function PPM1D mutations in brainstem gliomas. Nat Genet 46:726–730

    CAS  PubMed  PubMed Central  Google Scholar 

  39. 39.

    Zehir A, Benayed R, Shah RH, Syed A, Middha S, Kim HR et al (2017) Mutational landscape of metastatic cancer revealed from prospective clinical sequencing of 10,000 patients. Nat Med 23:703–713

    CAS  PubMed  PubMed Central  Google Scholar 

  40. 40.

    Labreche K, Simeonova I, Kamoun A, Gleize V, Chubb D, Letouze E et al (2015) TCF12 is mutated in anaplastic oligodendroglioma. Nat Commun 6:7207

    CAS  PubMed  PubMed Central  Google Scholar 

  41. 41.

    Kovaleva V, Geissler AL, Lutz L, Fritsch R, Makowiec F, Wiesemann S et al (2016) Spatio-temporal mutation profiles of case-matched colorectal carcinomas and their metastases reveal unique de novo mutations in metachronous lung metastases by targeted next generation sequencing. Mol Cancer 15:63

    PubMed  PubMed Central  Google Scholar 

  42. 42.

    Urick ME, Rudd ML, Godwin AK, Sgroi D, Merino M, Bell DW (2011) PIK3R1 (p85alpha) is somatically mutated at high frequency in primary endometrial cancer. Cancer Res 71:4061–4067

    CAS  PubMed  PubMed Central  Google Scholar 

  43. 43.

    Upadhyaya M, Maynard J, Osborn M, Harper PS (1997) Six novel mutations in the neurofibromatosis type 1 (NF1) gene. Hum Mutat 10:248–250

    CAS  PubMed  Google Scholar 

  44. 44.

    Wei X, Walia V, Lin JC, Teer JK, Prickett TD, Gartner J et al (2011) Exome sequencing identifies GRIN2A as frequently mutated in melanoma. Nat Genet 43:442–446

    CAS  PubMed  PubMed Central  Google Scholar 

  45. 45.

    Jones C, Baker SJ (2014) Unique genetic and epigenetic mechanisms driving paediatric diffuse high-grade glioma. Nat Rev Cancer 14. 651-61

  46. 46.

    Kumar S, Lu B, Dixit U, Hossain S, Liu Y, Li J et al (2015) Reciprocal regulation of Abl kinase by Crk Y251 and Abi1 controls invasive phenotypes in glioblastoma. Oncotarget 6:37792–37807

    PubMed  PubMed Central  Google Scholar 

  47. 47.

    Takino T, Nakada M, Miyamori H, Yamashita J, Yamada KM, Sato H (2003) CrkI adapter protein modulates cell migration and invasion in glioblastoma. Cancer Res 63:2335–2337

    CAS  PubMed  Google Scholar 

  48. 48.

    Zhang YK, Qu YY, Lin Y, Wu XH, Chen HZ, Wang X et al (2017) Enoyl-CoA hydratase-1 regulates mTOR signaling and apoptosis by sensing nutrients. Nat Commun 8:464

    PubMed  PubMed Central  Google Scholar 

  49. 49.

    Zarghooni M, Bartels U, Lee E, Buczkowicz P, Morrison A, Huang A et al (2010) Whole-genome profiling of pediatric diffuse intrinsic pontine gliomas highlights platelet-derived growth factor receptor alpha and poly (ADP-ribose) polymerase as potential therapeutic targets. J Clin Oncol 28:1337–1344

    CAS  PubMed  Google Scholar 

  50. 50.

    Burma S, Chen BP, Murphy M, Kurimasa A, Chen DJ (2001) ATM phosphorylates histone H2AX in response to DNA double-strand breaks. J Biol Chem 276:42462–42467

    CAS  PubMed  Google Scholar 

  51. 51.

    Zambetti GP, Levine AJ (1993) A comparison of the biological activities of wild-type and mutant p53. FASEB J 7:855–865

    CAS  PubMed  Google Scholar 

  52. 52.

    Graham V, Khudyakov J, Ellis P, Pevny L (2003) SOX2 functions to maintain neural progenitor identity. Neuron 39:749–765

    CAS  PubMed  Google Scholar 

  53. 53.

    Hua X, Miller ZA, Wu G, Shi Y, Lodish HF (1999) Specificity in transforming growth factor beta-induced transcription of the plasminogen activator inhibitor-1 gene: interactions of promoter DNA, transcription factor muE3, and Smad proteins. Proc Natl Acad Sci U S A 96:13130–13135

    CAS  PubMed  PubMed Central  Google Scholar 

  54. 54.

    Hjelmeland MD, Hjelmeland AB, Sathornsumetee S, Reese ED, Herbstreith MH, Laping NJ et al (2004) SB-431542, a small molecule transforming growth factor-beta-receptor antagonist, inhibits human glioma cell line proliferation and motility. Mol Cancer Ther 3:737–745

    CAS  PubMed  Google Scholar 

  55. 55.

    Hjortland GO, Bjornland K, Pettersen S, Garman-Vik SS, Emilsen E, Nesland JM et al (2003) Modulation of glioma cell invasion and motility by adenoviral gene transfer of PAI-1. Clin Exp Metastasis 20:301–309

    CAS  PubMed  Google Scholar 

  56. 56.

    Bourguignon LY, Zhu H, Shao L, Chen YW (2001) CD44 interaction with c-Src kinase promotes cortactin-mediated cytoskeleton function and hyaluronic acid-dependent ovarian tumor cell migration. J Biol Chem 276:7327–7336

    CAS  PubMed  Google Scholar 

  57. 57.

    Merzak A, Koocheckpour S, Pilkington GJ (1994) CD44 mediates human glioma cell adhesion and invasion in vitro. Cancer Res 54:3988–3992

    CAS  PubMed  Google Scholar 

  58. 58.

    Radotra B, McCormick D (1997) Glioma invasion in vitro is mediated by CD44-hyaluronan interactions. J Pathol 181:434–438

    CAS  PubMed  Google Scholar 

  59. 59.

    Gao T, Brognard J, Newton AC (2008) The phosphatase PHLPP controls the cellular levels of protein kinase C. J Biol Chem 283:6300–6311

    CAS  PubMed  Google Scholar 

  60. 60.

    Kallappagoudar S, Yadav RK, Lowe BR, Partridge JF (2015) Histone H3 mutations--a special role for H3.3 in tumorigenesis? Chromosoma 124:177–189

    CAS  PubMed  PubMed Central  Google Scholar 

  61. 61.

    Maeda S, Ohka F, Okuno Y, Aoki K, Motomura K, Takeuchi K et al (2020) H3F3A mutant allele specific imbalance in an aggressive subtype of diffuse midline glioma, H3 K27M-mutant. Acta Neuropathol Commun 8:8

    CAS  PubMed  PubMed Central  Google Scholar 

  62. 62.

    Audia JE, Campbell RM (2016) Histone modifications and Cancer. Cold Spring Harb Perspect Biol 8:a019521

    PubMed  PubMed Central  Google Scholar 

  63. 63.

    Abdel-Wahab O, Adli M, LaFave LM, Gao J, Hricik T, Shih AH et al (2012) ASXL1 mutations promote myeloid transformation through loss of PRC2-mediated gene repression. Cancer Cell 22:180–193

    CAS  PubMed  PubMed Central  Google Scholar 

  64. 64.

    Comet I, Riising EM, Leblanc B, Helin K (2016) Maintaining cell identity: PRC2-mediated regulation of transcription and cancer. Nat Rev Cancer 16:803–810

    CAS  PubMed  Google Scholar 

  65. 65.

    Justin N, Zhang Y, Tarricone C, Martin SR, Chen S, Underwood E et al (2016) Structural basis of oncogenic histone H3K27M inhibition of human polycomb repressive complex 2. Nat Commun 7:11316

    CAS  PubMed  PubMed Central  Google Scholar 

  66. 66.

    Bender S, Tang Y, Lindroth AM, Hovestadt V, Jones DT, Kool M et al (2013) Reduced H3K27me3 and DNA hypomethylation are major drivers of gene expression in K27M mutant pediatric high-grade gliomas. Cancer Cell 24:660–672

    CAS  PubMed  Google Scholar 

  67. 67.

    Chan KM, Fang D, Gan H, Hashizume R, Yu C, Schroeder M et al (2013) The histone H3.3K27M mutation in pediatric glioma reprograms H3K27 methylation and gene expression. Genes Dev 27:985–990

    CAS  PubMed  PubMed Central  Google Scholar 

  68. 68.

    Stafford JM, Lee CH, Voigt P, Descostes N, Saldana-Meyer R, Yu JR et al (2018) Multiple modes of PRC2 inhibition elicit global chromatin alterations in H3K27M pediatric glioma. Sci Adv 4:eaau5935

    CAS  PubMed  PubMed Central  Google Scholar 

  69. 69.

    Kornmann M, Hebart H, Danenberg K, Goeb R, Staib L, Kron M et al (2012) Response prediction in metastasised colorectal cancer using intratumoural thymidylate synthase: results of a randomised multicentre trial. Eur J Cancer 48:1443–1451

    CAS  PubMed  Google Scholar 

  70. 70.

    Rao S, Beckman RA, Riazi S, Yabar CS, Boca SM, Marshall JL et al (2017) Quantification and expert evaluation of evidence for chemopredictive biomarkers to personalize cancer treatment. Oncotarget 8:37923–37934

    PubMed  Google Scholar 

  71. 71.

    Shreeram S, Hee WK, Demidov ON, Kek C, Yamaguchi H, Fornace AJ Jr et al (2006) Regulation of ATM/p53-dependent suppression of myc-induced lymphomas by Wip1 phosphatase. J Exp Med 203:2793–2799

    CAS  PubMed  PubMed Central  Google Scholar 

  72. 72.

    Li J, Hart RP, Mallimo EM, Swerdel MR, Kusnecov AW, Herrup K (2013) EZH2-mediated H3K27 trimethylation mediates neurodegeneration in ataxia-telangiectasia. Nat Neurosci 16:1745–1753

    PubMed  PubMed Central  Google Scholar 

  73. 73.

    Furuta T, Takemura H, Liao ZY, Aune GJ, Redon C, Sedelnikova OA et al (2003) Phosphorylation of histone H2AX and activation of Mre11, Rad50, and Nbs1 in response to replication-dependent DNA double-strand breaks induced by mammalian DNA topoisomerase I cleavage complexes. J Biol Chem 278:20303–20312

    CAS  PubMed  Google Scholar 

  74. 74.

    Park EJ, Chan DW, Park JH, Oettinger MA, Kwon J (2003) DNA-PK is activated by nucleosomes and phosphorylates H2AX within the nucleosomes in an acetylation-dependent manner. Nucleic Acids Res 31:6819–6827

    CAS  PubMed  PubMed Central  Google Scholar 

  75. 75.

    Durant ST, Zheng L, Wang Y, Chen K, Zhang L, Zhang T et al (2018) The brain-penetrant clinical ATM inhibitor AZD1390 radiosensitizes and improves survival of preclinical brain tumor models. Sci Adv 4:eaat1719

    PubMed  PubMed Central  Google Scholar 

  76. 76.

    Ray Chaudhuri A, Nussenzweig A (2017) The multifaceted roles of PARP1 in DNA repair and chromatin remodelling. Nat Rev Mol Cell Biol 18:610–621

    CAS  PubMed  PubMed Central  Google Scholar 

  77. 77.

    Bryant HE, Schultz N, Thomas HD, Parker KM, Flower D, Lopez E et al (2005) Specific killing of BRCA2-deficient tumours with inhibitors of poly (ADP-ribose) polymerase. Nature 434:913–917

    CAS  PubMed  Google Scholar 

  78. 78.

    Lu R, O'Rourke JJ, Sobinoff AP, Allen JAM, Nelson CB, Tomlinson CG et al (2019) The FANCM-BLM-TOP3A-RMI complex suppresses alternative lengthening of telomeres (ALT). Nat Commun 10:2252

    PubMed  PubMed Central  Google Scholar 

  79. 79.

    Pan X, Drosopoulos WC, Sethi L, Madireddy A, Schildkraut CL, Zhang D (2017) FANCM, BRCA1, and BLM cooperatively resolve the replication stress at the ALT telomeres. Proc Natl Acad Sci U S A 114:E5940–E59E9

    CAS  PubMed  PubMed Central  Google Scholar 

  80. 80.

    Silva B, Pentz R, Figueira AM, Arora R, Lee YW, Hodson C et al (2019) FANCM limits ALT activity by restricting telomeric replication stress induced by deregulated BLM and R-loops. Nat Commun 10:2253

    PubMed  PubMed Central  Google Scholar 

  81. 81.

    Kline CN, Joseph NM, Grenert JP, van Ziffle J, Yeh I, Bastian BC et al (2016) Inactivating MUTYH germline mutations in pediatric patients with high-grade midline gliomas. Neuro-Oncology 18:752–753

    PubMed  PubMed Central  Google Scholar 

  82. 82.

    Kleczko EK, Heasley LE (2018) Mechanisms of rapid cancer cell reprogramming initiated by targeted receptor tyrosine kinase inhibitors and inherent therapeutic vulnerabilities. Mol Cancer 17:60

    PubMed  PubMed Central  Google Scholar 

  83. 83.

    Akhavan D, Pourzia AL, Nourian AA, Williams KJ, Nathanson D, Babic I et al (2013) De-repression of PDGFRbeta transcription promotes acquired resistance to EGFR tyrosine kinase inhibitors in glioblastoma patients. Cancer Discov 3:534–547

    CAS  PubMed  PubMed Central  Google Scholar 

  84. 84.

    Ware KE, Marshall ME, Heasley LR, Marek L, Hinz TK, Hercule P et al (2010) Rapidly acquired resistance to EGFR tyrosine kinase inhibitors in NSCLC cell lines through de-repression of FGFR2 and FGFR3 expression. PLoS One 5:e14117

    CAS  PubMed  PubMed Central  Google Scholar 

  85. 85.

    Keller S, Schmidt MHH (2017) EGFR and EGFRvIII promote angiogenesis and cell invasion in Glioblastoma: combination therapies for an effective treatment. Int J Mol Sci 18:1295

    PubMed Central  Google Scholar 

  86. 86.

    Molina JR, Hayashi Y, Stephens C, Georgescu MM (2010) Invasive glioblastoma cells acquire stemness and increased Akt activation. Neoplasia 12:453–463

    CAS  PubMed  PubMed Central  Google Scholar 

  87. 87.

    Pennacchietti S, Michieli P, Galluzzo M, Mazzone M, Giordano S, Comoglio PM (2003) Hypoxia promotes invasive growth by transcriptional activation of the met protooncogene. Cancer Cell 3:347–361

    PubMed  Google Scholar 

  88. 88.

    Takahashi Y, Morales FC, Kreimann EL, Georgescu MM (2006) PTEN tumor suppressor associates with NHERF proteins to attenuate PDGF receptor signaling. EMBO J 25:910–920

    CAS  PubMed  PubMed Central  Google Scholar 

  89. 89.

    Kubala MH, DeClerck YA (2019) The plasminogen activator inhibitor-1 paradox in cancer: a mechanistic understanding. Cancer Metastasis Rev 38:483–492

    PubMed  Google Scholar 

  90. 90.

    Rhyasen GW, Hattersley MM, Yao Y, Dulak A, Wang W, Petteruti P et al (2016) AZD5153: a novel bivalent BET Bromodomain inhibitor highly active against hematologic malignancies. Mol Cancer Ther 15:2563–2574

    CAS  PubMed  Google Scholar 

Download references

Acknowledgements

Special acknowledgements go to the patients’ parents and their commitment for research and progress in pediatric brain cancer. The authors thank Michael Caldwell and Jonathan Ball from Tempus for help with NGS, and Corey Gemelli from LSU/Shreveport for excellent technical assistance with IHC. We are especially grateful to Drs. Charles Timmons and Kathreena Kurian, as well as to Nikki Dauzat from STEPScenla for logistic assistance.

Funding

This work was supported by awards from Feist-Weiller Cancer Center and NeuroMarkers PLLC [NM2020–1] to M.-M.G.

Author information

Affiliations

Authors

Contributions

MMG performed conception and design of the study, development of methodology, acquisition of data, analysis and interpretation of data, writing and review of the manuscript, study supervision and provided material support. MZI performed acquisition and analysis of data and review of the manuscript. YL, MC, JT, CN and CK performed acquisition of data. DB performed acquisition of data and review of the manuscript. The author(s) read and approved the final manuscript.

Corresponding author

Correspondence to M.-M. Georgescu.

Ethics declarations

Ethics approval and consent to participate

The autopsies were performed in accordance to hospital ethical guidelines and regulations. The parents consented the patients’ autopsies for diagnosis, research and publication.

Consent for publication

The parents consented the patients’ autopsies for diagnosis, research and publication.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1

: Supplemental Table S1. Primary antibodies for WB. Supplemental Table S2. Autopsy gross parameters. Supplemental Table S3. IHC results. Supplemental Table S4. Mutations in DIPG. Supplemental Table S5. CNVs. Supplemental Figure S1. Radiologic characteristics of DIPG cases. A. F5 MRI: axial (left and center) and coronal (right) T2W-FLAIR images showing hyperintense diffusely infiltrating pontine mass extending posteriorly into the right cerebellum and rostrally into the midbrain and subthalamus. Images were acquired 2.5 months pre-mortem. B-C. F10 MRI at initial diagnosis (B) and 10 months later (C) showing rim-enhancing pontocerebellar mass with initial extension into the left cerebellum (B) and subsequent extension into the right cerebellum and frontal periventricular region (C). Yellow arrows show post-contrast enhancing tumor foci. D. F12 axial and sagittal T2W-FLAIR images acquired 1-month pre-mortem showing the frontal and cerebellar secondary foci with yellow arrows. Supplemental Figure S2. Example of histologic examination for H&E quantified analysis. F12 H&E sections are shown at the levels schematically illustrated on a normal brain coronal section. Higher magnification of areas marked with small blue squares are shown in adjacent pictures with blue borders. Note massive infiltration of periventricular areas, white matter tracts, septum pellucidum, hippocampus, brainstem and cerebellum (red arrow shows the invasive tumor front in the left cerebellum) by neoplastic cells with pleomorphic nuclei resembling glioblastoma. LV, lateral ventricle; R, right; L, left. Supplemental Figure S3. Variable ultrastructural morphology of neoplastic cells in DIPG. A. Variable nuclear morphology of F5 neoplastic cells, showing elongated/fibrillary forms (yellow arrows) admixed with larger forms (blue arrows). The red arrow indicates a myelinated axon “hugged” by a fibrillary neoplastic cell. B. F10 larger neoplastic astrocyte with abundant cytoplasm (blue arrow) dissecting the myelinated axons (red arrow) of the corpus callosum. C. F12 neoplastic astrocytes with nuclear pleomorphism (blue arrows) within the frontal white matter (myelinated axons shown with red arrows and oligodendroglia shown with green arrow). D. F12 large pontine fibrillary neoplastic astrocyte with prominent nucleolus. Supplemental Figure S4. Ganglioglioma-like pontine focus in DIPG. F10 H&E of a focal proliferation of large binucleated dysplastic ganglion cells, interspersed with small astrocytic neoplastic cells, conferring a typical ganglioglioma appearance to a limited pontine area of the tumor. IHC and NGS following microdissection of this area revealed histone H3 K27M mutation in all the neoplastic cells, consistent with DMG/K27M. Supplemental Figure S5. IHC profiles in DIPG. IHC with indicated antibodies of the pontine tumor areas from F5, F10 biopsy and autopsy, and F12. The F10 autopsy focus of dysplastic binucleated ganglion cells (red arrows) is shown. These cells express GFAP but not vimentin, whereas the small neoplastic cells express both types of intermediate filaments. Supplemental Figure S6. Quantitative protein expression analysis in DIPG. A-B. Bar graphs representing the indicated protein levels from various foci of the 3 DIPG autopsies. The quantification of the WB bands was performed by densitometric analysis, as described in Material and Methods. Individual densitometric values were normalized to the corresponding actin values. In (B), the phosphoprotein values were normalized to the corresponding total unphosphorylated protein values, as indicated. Results are expressed as fold-increase or fold-decrease in comparison to normal control. The WBs for each antibody were repeated at least twice, with similar results. The labeling is as in Fig. 4: T, tumor; N, normal; primary pontine (P) foci indicated in red; secondary cerebellar (C) and frontal (Fr) foci, in green. Sample laterality: R, right; L, left. For the tri-methylated histone H3 K27 residue (H3K27Me3) antibody in (A) that showed undetectable expression in both normal and tumor samples, the WB was repeated in the presence of a positive control from an autopsy case of glioblastoma, IDH-mutant, WHO grade IV, known to harbor increased levels of H3 K27 tri-methylation.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Georgescu, M., Islam, M.Z., Li, Y. et al. Global activation of oncogenic pathways underlies therapy resistance in diffuse midline glioma. acta neuropathol commun 8, 111 (2020). https://doi.org/10.1186/s40478-020-00992-9

Download citation

Keywords

  • Diffuse midline glioma, Autopsy, Next generation sequencing (NGS), Copy number variation (CNV), Proteomics