- Open Access
Characterizing temporal genomic heterogeneity in pediatric high-grade gliomas
Acta Neuropathologica Communications volume 5, Article number: 78 (2017)
Pediatric high-grade gliomas (pHGGs) are aggressive neoplasms representing approximately 20% of brain tumors in children. Current therapies offer limited disease control, and patients have a poor prognosis. Empiric use of targeted therapy, especially at progression, is increasingly practiced despite a paucity of data regarding temporal and therapy-driven genomic evolution in pHGGs. To study the genetic landscape of pHGGs at recurrence, we performed whole exome and methylation analyses on matched primary and recurrent pHGGs from 16 patients. Tumor mutational profiles identified three distinct subgroups. Group 1 (n = 7) harbored known hotspot mutations in Histone 3 (H3) (K27M or G34V) or IDH1 (H3/IDH1 mutants) and co-occurring TP53 or ACVR1 mutations in tumor pairs across the disease course. Group 2 (n = 7), H3/IDH1 wildtype tumor pairs, harbored novel mutations in chromatin modifiers (ZMYND11, EP300 n = 2), all associated with TP53 alterations, or had BRAF V600E mutations (n = 2) conserved across tumor pairs. Group 3 included 2 tumors with NF1 germline mutations. Pairs from primary and relapsed pHGG samples clustered within the same DNA methylation subgroup. ATRX mutations were clonal and retained in H3G34V and H3/IDH1 wildtype tumors, while different genetic alterations in this gene were observed at diagnosis and recurrence in IDH1 mutant tumors. Mutations in putative drug targets (EGFR, ERBB2, PDGFRA, PI3K) were not always shared between primary and recurrence samples, indicating evolution during progression. Our findings indicate that specific key driver mutations in pHGGs are conserved at recurrence and are prime targets for therapeutic development and clinical trials (e.g. H3 post-translational modifications, IDH1, BRAF V600E). Other actionable mutations are acquired or lost, indicating that re-biopsy at recurrence will provide better guidance for effective targeted therapy of pHGGs.
Genetic and epigenetic molecular profiling techniques have revolutionized our understanding of the etiology and biology of pediatric high-grade gliomas (pHGGs) (reviewed in ). Unfortunately, this has not yet led to an improvement in outcome for children with this disease  despite the use of agents that target pathways identified through these biological advances. Novel agents for the treatment of pHGGs are first tested in the relapse setting, and target genomic alterations typically present in therapy-naïve diagnostic tumor samples or models. However, there is limited data on the relevance of genomic aberrations at diagnosis on disease progression after multimodal therapy, making the effectiveness of this approach questionable. An improved understanding of temporal and therapy-driven evolution of recurrent pHGGs is therefore needed, especially in the context of hemispheric HGGs that show increased genetic heterogeneity [5, 12, 13, 19, 37, 50, 51].
Clonal evolution is a dynamic process that has been reported in many cancer types [3, 28, 39, 48], even without exposure to therapy . Morrissy et al., have recently demonstrated poor overlap in genetic events between primary and post-treatment medulloblastoma both in murine models and human samples . This included a marked divergence in actionable genes between diagnosis and recurrence, despite conservation of molecular subgroup affiliation [28, 36, 47]. Whole exome sequencing (WES) of 23 initial and recurrent gliomas in adults by Johnson et al., revealed variable genetic relatedness across pairs; in 10 cases, most mutations from diagnosis were not conserved in the recurrent sample, including the BRAF V600E hotspot mutation . In adult glioblastoma multiforme (GBM), a longitudinal study of the genetic landscape of 114 untreated and recurrent paired tumors revealed a switch in expression-based subtypes in 63% of cases. Enrichment of a hypermutated phenotype in recurrent disease exposed to temozolomide (TMZ) was also identified, suggesting the occurrence of therapy-induced mutagenesis . Moreover, an analysis of tumor phylogeny revealed that dominant clones at recurrence were infrequently direct descendants of dominant clones from diagnosis . We have previously shown that disease-defining somatic mutations in oncohistones [K27M in Histone 3 (H3) variants (H3F3A, HIST1H3B)] are spatially stable in diffuse intrinsic pontine glioma (DIPG), and co-occur with highly conserved partners throughout geographically distinct tumor sites [18, 30]. However, limited data on disease recurrence are available for supratentorial pHGGs. This is of major therapeutic interest as hemispheric pHGGs show more genetic variability at diagnosis than midline tumors, the vast majority of which are defined by H3K27M mutations (> 90%) [14, 51]. In the current study, we characterize the temporal genomic heterogeneity in pHGGs by assessing the mutational profile and methylome of paired primary and recurrent tumors with emphasis on supratentorial pHGGs.
Materials and methods
Institutional review board approval was obtained to perform this retrospective study at Cincinnati Children’s Hospital Medical Center (CCHMC, Study ID: 2014-6849) and Nationwide Children’s Hospital (NCH: IRB15-00143). The patient cohort was chosen based on the availability of material from both the primary and recurrent tumor for each case with a confirmed HGG diagnosis Two neuropathologists (CF and JK) independently reviewed tumor samples. Patient tumor samples were acquired from diagnosis as well as recurrence or autopsy and preserved either as fresh-frozen or formalin fixed paraffin embedded (FFPE) tissue. Blood or other matched normal tissue was obtained when available for germline analysis. To ensure adequate tumor content, hematoxylin and eosin (H&E) slides were reviewed from each frozen specimen, the initial cut of each FFPE block, and an additional cut of FFPE block after scrolls were obtained for DNA extraction. All patient tumor and matched blood samples were collected after informed consent was provided by patients or legal guardians through institutional review board approved protocols at the respective institutions.
DNA extraction was carried out from frozen tissue using the Qiagen AllPrep DNA/RNA/miRNA Universal Kit following the manufacturer’s instructions. DNA from FFPE scrolls or core punches were isolated by suspending the paraffin scrolls in deparaffinization solution (Qiagen) followed by DNA extraction using the QIAamp DNA FFPE Tissue Kit. DNA quantification was conducted using the Quant-iT Picogreen dsDNA assay kit (Thermo Fisher Scientific). Droplet digital PCR (ddPCR) assays for H3K27M mutations were performed as previously described .
Whole Exome Sequencing (WES) analysis
The Nextera Rapid Capture Exome kit (Illumina) was used to prepare 36 libraries, and the Agilent SureSelect Reagent Exome kit (Agilent) was used to prepare 6 libraries according to the manufacturer’s instructions. Genomic DNA was extracted from frozen tissue and FFPE blocks representing tumor or normal tissue and from monocytes. Sequencing was performed on the Illumina HiSeq 2000 using rapid-run mode with 100 bp paired-end reads. Adaptor sequences were removed, and reads trimmed for quality using the FASTX-Toolkit (http://hannonlab.cshl.edu/fastx_toolkit/). An in-house program was used to ensure the presence of exclusively paired-reads. We next aligned the reads using Burrows-Wheeler Aligner (BWA) 0.7.7 to GRC37/hg19 as a reference genome. Indel realignment was performed using the Genome Analysis Toolkit (GATK) 29 (https://software.broadinstitute.org/gatk/). Duplicate reads were marked using Picard (http://broadinstitute.github.io/picard/), and excluded from further analyses. The average coverage for all the samples was 69X. Single Nucleotide Variants (SNVs) and short indels were called using our in-house pipeline that exploits three different variant callers: FreeBayes 1.1.0 (https://arxiv.org/abs/1207.3907), SAMtools 1.3.1 (http://samtools.sourceforge.net/) and GATK HaplotypeCaller 3.7 . Thresholds were set for calling a true variant to two out of three variant callers. Next, variants were filtered for quality so at least 10% of reads supported each variant call. ANNOVAR  and in-house programs were used to annotate variants that affect protein-coding sequence. Variants were screened to assess whether they had previously been observed in public datasets including the 1000 Genomes Project data set (November 2011), the National Heart, Lung and Blood Institute (NHLBI) Grand Opportunity (GO) exomes as well as in over 3000 exomes previously sequenced at our center (including cancer and non-cancer samples).
Somatic and putative somatic mutation identification from Whole Exome Sequencing
Protein coding variants were identified as nonsynonymous missense, frameshifts, stopgains, indels, splice variants, and present with a minimum of 5% mutant allele frequency (MAF) or greater. A MAF of > 5% was used in order to reduce removal of low frequency mutations in genes of interest (H3F3A, TP53, ATRX, ZMYND11, LZTR1). All other reported mutations have a MAF > 15%. To be considered somatic, zero mutant variant reads were present in the normal sample. In some cases where < 2 variant reads were present in the normal, the sequencing alignments were manually checked to verify as sequencing artifacts. Additionally, all somatic and putative somatic variants were manually checked for alignment and sequencing artifacts. To further remove germline SNPs where normal was not available, all samples were filtered based on a MAF (mutant allele frequency) < 0.0005 in the 1000 Genomes Project (November 2011), EVS, and < 0.00005 from ExAC databases. Finally, in cases with no associated normal, variants were only reported when present in COSMIC or high/medium functional impact as assessed using MutationAssessor (http://mutationassessor.org/r3/). Targeted validations on selected variants were performed using Sanger sequencing or high-depth sequencing on an Illumina Miseq as previously described .
Allelic Imbalance and CNV analysis
Allelic Imbalance was assessed in whole exome sequencing data using the ExomeAI program  for 15 pHGG tumor pairs using the default parameters of the program. HGG1 was not analyzed due to low WES coverage of the recurrence tumor. Copy number variations (CNVs) were analyzed in 8 pHGG tumor-normal pairs using an in-house program (CNAXX; unpublished) we developed that takes both coverage (normalized average) and the deviation of B allele frequency from 50% into account (adapted from methods used in FishingCNV  and ExomeAI ). Different CNV events (amplification, deletion and copy neutral LOH) were called based upon the status of the normalized coverage and the B allelic imbalance as we described previously in . We assessed the CNV events at both the chromosomal arm level and the union of the segments called by each of the two features used by the program (i.e. B allele frequency and normalized average coverage). Paired t-tests were performed using the statistical program R version 3.1.1.
DNA methylation analysis
Illumina 450 K methylation chips were used on 28 samples, and profiling data was analyzed as previously described [14, 32]. In some cases, only the tumor from the primary (HGG1, HGG3, HGG11) or recurrence (HGG12) were included in the final analysis due to insufficient or poor quality FFPE DNA from the respective tumor pair. The raw data were subject to quality control and preprocessing utilizing the R package minfi, and normalized for technical variation between the Infinium I and II probes using the SWAN method. We removed probes on sex chromosomes (chr X, Y), those containing SNPs (dbSNP: https://www.ncbi.nlm.nih.gov/projects/SNP/) as well as non-specific probes that bind to multiple genomic locations. Unsupervised hierarchical clustering was performed using average linkage, and Pearson rank correlation distance on the top 3000 most variable probes selected based on standard deviation of beta values (β-values).
Immunohistochemistry for MMR proteins
Immunohistochemistry (IHC) staining for expression of MMR proteins (MLH1, MSH2, MSH6, PMS2) was performed on slides cut from FFPE blocks of pHGG samples using conventional methods .
Paired primary and recurrent pHGG tumor samples were available from 16 patients. Median age at diagnosis was 15 years (range: 4–29 years). Among the 16 cases, one patient with a germline NF1 mutation was initially diagnosed with a WHO grade II glioma with pilocytic features (HGG16) before recurring as a HGG. Median time to progression was 13 months (range: 4–45 months). Fourteen patients received radiation therapy at diagnosis, and 10 patients (63%) had received temozolomide (TMZ) prior to last progression. Five tumors (31%) were midline [2 pontine tumors (HGG2, HGG4), 1 thalamic tumor (HGG8), 2 spinal tumors (HGG1, HGG3)], and 11 tumors (69%) were hemispheric. Two patients (HGG15 and HGG16) were previously known to have germline neurofibromatosis type 1 (NF1). All clinicopathological features are shown in Table 1.
Core oncogenic mutations are shared in primary and recurrent pHGGs
To determine the temporal stability of the mutational landscape of pHGGs, WES was performed on all 16 tumor pairs. Matched germline DNA was available for 9 patients. Samples were sequenced with an average coverage of 69X (range 3.5–200.4X) (Additional file 1: Table S1). Tumor pairs were stratified into three distinct molecular groups based on the identified mutation patterns (Fig. 1a, Additional file 2: Table S2). The first group was comprised of H3/IDH1 mutant tumors (7/16 pairs, 44%), where both the primary and recurrence harbored shared epigenetic driver mutations in H3F3A K27M (n = 3, HGG1, HGG2, HGG3), HIST1H3B K27M (n = 1, HGG4), H3F3A G34V (n = 1, HGG5), or IDH1 R132H/S (n = 2, HGG6, HGG7). Consistent with previous studies, all H3K27M mutant tumors were in the brain midline (spine and pons), while H3F3A G34 V and IDH1 mutant tumors were hemispheric [37, 41]. In addition, all six H3F3A and IDH1 mutant tumors had co-occurring TP53 mutations while the HIST1H3B (Histone 3.1) K27M mutant DIPG had a known associated ACVR1 R258G mutation (HGG4) [14, 42, 51]. H3K27M/G34V and IDH1 mutations were as expected clonal and conserved across primary and recurrent tumor pairs. Similarly, ACVR1 and TP53 mutations, which have been shown to be obligate partners in H3 and IDH mutagenesis, were also conserved throughout disease progression.
The second group had no identifiable H3 or IDH1 mutations, termed H3/IDH1 wildtype (7/16 pairs, 44%). This is a heterogeneous group of supratentorial, mostly hemispheric, tumors (6/7, 86%) with clonal mutations identified in TP53 (4/7, 57%), ATRX (3/7, 43%), BRAF V600E (2/7, 29%), BCOR (1/7, 14%), CDKN2A (1/7, 14%), ZMYND11 (1/7, 14%), and EP300 (1/7, 14%) in tumor pairs (Fig. 1a). HGG8, the only H3/IDH1 thalamic wildtype tumor, had TP53 missense, ATRX missense, CDKN2A nonsense, and interestingly, EP300 missense mutations that were retained at recurrence. One sample (HGG11) had a hypermutated phenotype with 151 protein coding somatic mutations at diagnosis and 670 at relapse. The H3/IDH1 wildtype group was enriched for multiple mutations in receptor tyrosine kinase (RTK), phosphoinositide 3-kinase (PI3K) pathway, histone modifiers, integrin, and cadherin genes (Fig. 1a, b). Except for rare mutations, most of these genetic alterations were not present across tumor pairs, and seemed to occur more frequently at disease recurrence (Fig. 1a, b).
Lastly, the third molecular group was composed of patients with NF1 germline truncating mutations (2/16, 12%) with associated HGG (HGG15, HGG16). There were few SNVs identified in the primary tumors from these 2 cases. The primary tumor of HGG15 had no detectable driver SNVs, however, it acquired an ATRX frameshift mutation as well as a RAD50 missense mutation at recurrence (Fig. 1b). In HGG16, both primary and recurrent samples harbored PPM1D nonsense (L513X) and ATRX missense mutations. Previous studies have shown that PPM1D mutations affect p53 function and are mutually exclusive of TP53 mutations . This mutation pattern in the NF1 mutant group is reminiscent of G34R/V-H3F3A mutated HGGs which co-occur with ATRX and TP53 alterations.
Other mutational patterns at diagnosis and recurrence
In all groups of pHGG tumors, core oncogenic driver mutations in H3 variants, IDH1, TP53, ACVR1, BRAF V600E, and PPM1D were conserved at recurrence. As expected, the TP53 pathway (TP53 or PPM1D) was the most frequently altered (11/16, 69%) across all tumor subgroups (Fig. 1a). ATRX mutations were also frequent (8/16, 50%) and, as previously described, were enriched in supratentorial samples (7/10, 70%) [1, 9, 14, 26]. In the H3/IDH1 mutant group, ATRX mutations were observed in 3 hemispheric tumor pairs: one pair with H3G34V (HGG5), and two with IDH1 R132H/S (HGG6, HGG7) mutations. Additionally, ATRX mutations were identified in three H3/IDH1 wildtype tumor pairs, two of which co-occurred with alterations in chromatin modifiers (HGG8, HGG9), and in both samples from patients with germline NF1. We observed the same ATRX mutation at recurrence in all tumor pairs of the H3/IDH1 wildtype subgroup. In IDH1 mutant tumors, HGG7 harbored two ATRX alterations including a frameshift mutation exclusive to the primary tumor, while HGG6 acquired an additional ATRX missense mutation in the recurrence, in keeping with previous findings in IDH1 mutagenesis . While taking into consideration the limitations imposed by our relatively small cohort, our findings indicate that in H3/IDH1 wildtype pHGGs, the same ATRX mutation, when present at diagnosis, is retained at disease recurrence despite gross total resection of the primary tumor in most cases. Interestingly, while loss of BRAF V600E at recurrence has been reported in adult gliomas , the two cases in our study (HGG12, HGG13) retained this mutation at recurrence.
Novel epigenetic alterations in H3/IDH1 wildtype pHGGs
Two tumor pairs in the H3/IDH1 wildtype group showed potentially novel epigenetic drivers that converge to affect the same histone mark directly affected in pHGGs carrying H3.3 G34R/V or H3K27M mutations. Both primary and recurrent samples from HGG9, an adolescent patient with a parietal brain tumor, shared a ZMYND11 frameshift mutation that possibly abrogates expression of this protein, and they also harbored concurrent TP53 and ATRX missense mutations. These characteristics are typical of pHGGs that harbor H3.3G34R/V mutations. ZMYND11/BS69 has been shown to specifically recognize H3K36me3 to regulate and repress transcription [17, 49]. Conversely, H3.3G34R/V mutant nucleosomes affect the trimethylation of K36 on the H3.3 mutant nucleosomes , impairing the recognition of H3K36me3 ZMYND11/BS69 and its action on modulating translation . In HGG8, both primary and recurrent thalamic tumors shared a mutation in EP300 (D1339N), concurrent with TP53 and ATRX missense mutations and a CDKN2A nonsense mutation. The EP300 gene encodes a histone acetyltransferase (HAT), with acetylation activity at H3K27, and has been shown to regulate transcription through multiple mechanisms of chromatin remodeling . The D1339N mutation has been identified in multiple tumor types, but not in HGG, and is presumed to affect the enzymatic activity of the protein [8, 25]. It is worth noting that HGG8 harboring mutant EP300 was a midline tumor, a location where H3K27M mutations known to affect H3K27 acetylation usually account for the vast majority of driver mutations. Taken together, these data suggest that both these novel epigenetic mutations are relevant in the setting of pHGGs, and may possibly mimic the effect of genetic drivers in H3-mutant tumors.
Activating mutations of the RTK-PI3K pathway are not always conserved
Mutations in the RTK genes involving EGFR, ERBB2, ERBB4, FLT1 and EPHA/B, were identified in 5/7 (71%) of the H3/IDH1 wildtype tumor pairs. Many of these were exclusive to the primary or recurrent sample (Fig. 1b). In the H3/IDH1 mutant tumor HGG5, a PDGFRA Y288C missense mutation was acquired at recurrence. Similar to previous reports, PI3K mutations were temporally heterogeneous , where 1/5 tumor (20%) exhibited a shared mutation in a PI3K regulatory subunit PIK3R3 (HGG12). Additionally, subclonal mutations in the PI3K catalytic subunit were private to the primary tumor only in one case (HGG3) and were acquired at recurrence in two cases (HGG4, HGG13). Prior to our study, testing of a different tumor sampling from HGG3 was performed using a clinical genomic panel. This analysis identified a low frequency PIK3CA H1047R hotspot mutation that was not found in either the WES analysis of a different primary tumor tissue block or targeted high-depth sequencing of multiple samplings of the recurrent tumor (Additional file 3: Figure S1). In the hypermutated HGG11 tumor pair, the primary tumor harbored a missense MTOR mutation, while at recurrence the tumor acquired a PI3K catalytic subunit PIK3CD passenger mutation (Fig. 1b).
Mutational burden, allelic imbalance and copy number variations
Analysis of the mutational burden showed no statistically significant difference in the number of mutations between primary and recurrent tumors across all groups (paired t-test, p = 0.24) (Fig. 2a, Additional file 4: Table S3, Additional file 5: Figure S2). It is worth noting that within the limitations of sample size, we observed a trend towards an increase in the mutational burden at recurrence that did not reach statistical significance despite the use of TMZ as adjuvant therapy in 10/16 (63%) pHGGs. In HGG11, we observed a marked increase in the number of somatic mutations in the primary (n = 151) and at recurrence (n = 670) compared to all other tumor samples, indicating a hypermutated phenotype. We identified and validated a germline MLH1 splice missense mutation, and also performed immunohistochemistry on MMR proteins (MLH1, MSH6, MSH2, PMS2) on the primary HGG11 tumor (Additional file 6: Figure S3). Although IHC results did not show loss of any MMR proteins, we hypothesize that the splice mutation that translated extra inframe amino acids (data not shown), resulted in a dysfunctional yet nuclear-localized MLH1 protein. This may explain MMR IHC nuclear positivity in the setting of mismatch repair deficiency resulting in hypermutation. Interestingly, the mutation burden in that case dramatically increased at recurrence, which may be attributable to the combined effects of radiation and TMZ treatment .
To further assess chromosomal alterations in all of the primary-recurrent tumor pairs, we used WES data to analyze the state of allelic imbalance using ExomeAI . Copy Number Variations (CNVs) were analyzed in eight tumor pairs with available matched normal. We calculated the numbers of allelic imbalance regions as shared or specific to the primary or recurrent tumor (Fig. 2b, Additional file 4: Table S3), regardless of the size of each region. Similar to mutation counts, there was no significant difference in the number of regions of allelic imbalance between the primary and recurrent tumors across all subgroups (paired t-test, p = 0.071). One tumor pair, HGG9, was particularly remarkable as there was an increased number of allelic imbalance regions in the recurrent tumor compared to the primary. Assessment of copy number variations confirmed genome-wide loss of heterozygosity events at recurrence resulting in a copy number neutral genome (Additional file 7: Figure S4), compatible with radiotherapy-induced chromosomal alterations [22, 54]. Both NF1 germline cases also showed an increase in the number of regions of allelic imbalance. In both NF1 germline recurrent tumors, CNV analysis showed loss of heterozygosity in the whole q arm of chromosome 17, containing the NF1 gene locus (Additional files 8 and 9: Tables S4 and S5). Our limited CNV analysis did not show any focal deletions or amplifications in genes of interest (PDGFRA, MYC, CDKN2A/B, CDK6) that have been previously implicated in pHGG . In 5 of 6 tumors with TP53 mutations, large copy number alterations in the p arm of chromosome 7 were present that included the TP53 gene loci (Additional file 8: Table S4).
DNA methylation subgroup is maintained at recurrence in pHGG
We performed unsupervised clustering of DNA methylation in tumor pairs integrated into a larger pHGG in-house data set (Fig. 3). As expected, tumors in the H3/IDH1 mutant group clustered within their respective methylation cluster, while H3/IDH1 wildtype and NF1 germline groups clustered with H3/IDH1 wildtype HGGs. In all cases, the recurrent tumors clustered within the same methylation subgroup, similar to findings in other brain tumors including DIPG and medulloblastoma [30, 31].
In this work, we performed whole exome sequencing on 16 primary and recurrent pHGG pairs including two pHGGs from patients with germline NF1 mutations, and provide insight into the temporal genomic evolution of these tumors. A direct comparison of the mutational landscape of paired samples reveals that oncogenic driver mutations are typically conserved. The identification of these mutations in both the primary and recurrent tumors suggests that these mutations are early initiating events in tumorigenesis, are stable, and unaffected by treatment. This is in contrast to adult GBM where cancer driver mutations can be subclonal in the primary and recurrent tumors .
In our dataset, 10 of 16 patients were treated with TMZ, and despite a trend towards an increase in the number of mutations at recurrence, there was no statistically significant increase in mutational burden. This is in contrast to adult GBM, where an increase in mutational burden is observed with TMZ treatment . Although our sample size is small, our findings may reflect a different biological process in response to TMZ in children compared to adults and warrants further evaluation. We observed one H3/IDH1 wildtype primary tumor (HGG11) with an increased number of somatic mutations compared to other primary tumors, and was identified to harbor a germline MLH1 splice missense mutation. Immunohistochemical analysis did not show loss of the MMR proteins, however, we hypothesize that the missense splice mutation likely resulted in the translation of a dysfunctional MLH1 protein product to cause mismatch repair deficiency (MMRD) and hypermutation. After treatment with radiation and TMZ, this tumor acquired an increased number of somatic mutations compared to the primary tumor, suggesting that treatment further exacerbated the hypermutated phenotype. Several controversial and contradictory studies have variably reported the presence of microsatellite instability which results in mismatch repair deficiency in pediatric HGG and adults [10, 44], highlighting the need for further studies. Future genetic testing for MMRD in pediatric HGG patients could steer treatment towards immunotherapy, as immune checkpoint blockade has shown clinical benefits in MMRD colorectal cancers as well as children with high-grade glioma [4, 23].
Similar to findings in adult IDH1-mutant gliomas , we identify heterogeneous ATRX alterations among IDH1 mutant pHGG tumor pairs. While IDH1 mutant tumors are more common in adult GBM and occur in up to 98% of secondary GBMs, they make up less than 10% of all pediatric HGGs [2, 52]. In contrast to IDH1-mutant gliomas, ATRX mutations associated with H3G34V, ZMYND11, EP300, or BRAF V600E were stable across the disease course in our study. Additionally, the BRAF V600E mutation was present in both primary and relapse samples in two children in our study which is in contrast to adult studies where it was identified either at diagnosis or at recurrence .
H3/IDH1 wildtype pHGGs have previously been shown to be a diverse group of tumors with mutations in many cancer pathways [35, 37, 51], but have not been directly linked to any particular epigenetic driver as is the case with H3 and IDH1 mutant tumors. Our data reflect the heterogeneity of tumors in the H3/IDH1 wildtype group while also identifying two novel pHGG epigenetic cancer drivers (ZMYND11 and EP300) in this group. ZMYND11 has recently been described as an epigenetic regulator that specifically interacts with H3K36me3 to regulate transcription. Wen et al. have reported that H3 G34R/V mutations impair binding of ZMYND11 to an H3.3K36me3 peptide, suggesting that H3.3 G34R/V and ZMYND11 mutations alter H3K36me3 levels in similar fashions . To the best of our knowledge, ZMYND11 mutations have not been previously described in pHGGs. The tumor harboring this mutation (HGG9) was located in the right parietal lobe and carried partner mutations in ATRX and TP53, further supporting its similarity to hemispheric H3.3 G34R/V mutated tumors. In addition, inactivating mutations identified in the HAT gene EP300 have been implicated in a wide array of cancer types including diffuse large B cell lymphoma , head and neck, esophageal, colorectal, medulloblastoma and non-small cell lung carcinoma [7, 15]. We also report a specific EP300 hotspot D1399N mutation (HGG8) which has not been previously identified in HGGs. Structural analysis of EP300 has shown that the D1399 residue has effects on the conformation of the HAT domain, specifically the L1 loop . This is also an inactivating mutation which abolishes autoacetylation required for HAT activity, thus affecting post-translational modification of K27 on H3 variants . Interestingly, EP300 D1399Y mutations alter its interaction with transcription factor AP-2alpha indirectly leading to the transactivation of Myc . Moreover, the tumor harboring the EP300 mutation was located in the thalamus which is a neuroanatomical structure in the brain midline where the majority of HGGs harbor H3K27M mutations. This novel epigenetic mutation may reproduce some of the effects of K27M in a wildtype H3K27 tumor. In our study, the tumor with the EP300 D1399N mutation had increased Myc expression (data not shown), suggesting that this particular EP300 mutation may also play a role in Myc-related oncogenesis similar to K27M mutagenesis. Although interesting, these findings need further testing and functional validation in relevant disease models. The two HGGs from patients with germline NF1 did not show a high mutational burden at diagnosis or at recurrence, and no clear associated driver mutation. Interestingly, a tendency towards increased copy number alteration was observed in both pairs at recurrence. These findings also need further validation on a larger sample set.
Somatic mutations in RTKs are common in adult GBM [5, 6] and are generally found at low frequencies in pHGGs . Similar to our previous report , the H3/IDH1 wildtype group in this study seemed enriched with RTK mutations (5/7, 71%). One striking finding in this molecular group was the discovery of EGFR missense mutations in the primary occurrence of HGG10 (T790M and E709A), which were lost in the recurrence. A shared EGFR R222C missense mutation was present in both the primary and recurrent tumors, indicating that alteration of the RTK pathway is nonetheless conserved in the recurrent tumor. The EGFR T790M mutation has been implicated in acquired resistance to most EGFR tyrosine kinase inhibitors [21, 27]. This may, in part, explain tumor progression in this patient despite treatment with lapatinib (Novartis, East Hanover, NJ), and highlights the importance of identifying resistance-promoting mutations in the clinical setting. We also identified three tumors with targetable RTK lesions (PDGFRA, ERBB2, ERBB4) that were exclusive to the recurrent tumor (HGG5, HG8, HGG11), indicating that genomic data from tumor tissue at recurrence may provide better guidance for therapeutic choices. Conversely, one case harbored a low level subclonal PIK3CA mutation that was discovered by a clinical genomics panel in the primary tumor, but was not identified by WES in different primary tumor blocks from the same case, nor in multiple samplings of the recurrent tumor. Excluding the subclonal nature of this mutation, and confirming its maintenance at recurrence are important therapeutic considerations before embarking on targeted treatment, especially with single agents such as rapamycin used in this patient.
In conclusion, this study further highlights the molecular distinction between pediatric and adult HGGs, especially in therapy-induced tumor evolution. We show that genes with driver mutations (H3, TP53, PPMID, ZMYND11, EP300) as well as some targetable mutations (e.g. IDH1, BRAF V600E) are conserved. Importantly, we demonstrate that some actionable mutations are unstable (PI3K, EGFR), indicating that re-biopsy is warranted in order to optimize personalized therapy. The presence of subclonal targetable alterations concurrently with driver mutations supports the use of combination therapy approaches to address disease biology and evolution with the aim of improving patient outcomes.
isocitrate dehydrogenase 1
tumor protein 53
activin a receptor, type I
zinc finger MYND domain-containing protein 11
histone acetyltransferase p300
alpha-thalassemia/mental retardation syndrome, Nondeletion type, x-linked
epidermal growth factor receptor
Erb-B2 receptor tyrosine kinase 2
platelet derived growth factor receptor alpha
Bai H, Harmanci AS, Erson-Omay EZ, Li J, Coskun S, Simon M et al (2016) Integrated genomic characterization of IDH1-mutant glioma malignant progression. Nat Genet 48:59–66. doi: 10.1038/ng.3457
Balss J, Meyer J, Mueller W, Korshunov A, Hartmann C, von Deimling A (2008) Analysis of the IDH1 codon 132 mutation in brain tumors. Acta Neuropathol 116:597–602. doi: 10.1007/s00401-008-0455-2
Beltran H, Prandi D, Mosquera JM, Benelli M, Puca L, Cyrta J et al (2016) Divergent clonal evolution of castration-resistant neuroendocrine prostate cancer. Nat Med 22:298–305. doi: 10.1038/nm.4045
Bouffet E, Larouche V, Campbell BB, Merico D, de Borja R, Aronson M et al (2016) Immune checkpoint inhibition for Hypermutant Glioblastoma Multiforme resulting from Germline Biallelic mismatch repair deficiency. J Clin Oncol 34:2206–2211. doi: 10.1200/jco.2016.66.6552
Brennan CW, Verhaak RG, McKenna A, Campos B, Noushmehr H, Salama SR et al (2013) The somatic genomic landscape of glioblastoma. Cell 155:462–477. doi: 10.1016/j.cell.2013.09.034
Cancer Genome Atlas Research N (2008) Comprehensive genomic characterization defines human glioblastoma genes and core pathways. Nature 455:1061–1068. doi: 10.1038/nature07385
Cerami E, Gao J, Dogrusoz U, Gross BE, Sumer SO, Aksoy BA et al (2012) The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. Cancer Discov 2:401–404. doi: 10.1158/2159-8290.CD-12-0095
Delvecchio M, Gaucher J, Aguilar-Gurrieri C, Ortega E, Panne D (2013) Structure of the p300 catalytic core and implications for chromatin targeting and HAT regulation. Nat Struct Mol Biol 20:1040–1046. doi: 10.1038/nsmb.2642
Eckel-Passow JE, Lachance DH, Molinaro AM, Walsh KM, Decker PA, Sicotte H et al (2015) Glioma groups based on 1p/19q, IDH, and TERT promoter mutations in tumors. N Engl J Med 372:2499–2508. doi: 10.1056/NEJMoa1407279
Eckert A, Kloor M, Giersch A, Ahmadi R, Herold-Mende C, Hampl JA et al (2007) Microsatellite instability in pediatric and adult high-grade gliomas. Brain Pathol 17:146–150. doi: 10.1111/j.1750-3639.2007.00049.x
Eirew P, Steif A, Khattra J, Ha G, Yap D, Farahani H et al (2015) Dynamics of genomic clones in breast cancer patient xenografts at single-cell resolution. Nature 518:422–426. doi: 10.1038/nature13952
Fontebasso A, Schwartzentruber J, Khuong-Quang D-A, Liu X-Y, Sturm D, Korshunov A et al (2013) Mutations in SETD2 and genes affecting histone H3K36 methylation target hemispheric high-grade gliomas. Acta Neuropathol 125:659–669. doi: 10.1007/s00401-013-1095-8
Fontebasso AM, Gayden T, Nikbakht H, Neirinck M, Papillon-Cavanagh S, Majewski J et al (2014) Epigenetic dysregulation: a novel pathway of oncogenesis in pediatric brain tumors. Acta Neuropathol 128:615–627. doi: 10.1007/s00401-014-1325-8
Fontebasso AM, Papillon-Cavanagh S, Schwartzentruber J, Nikbakht H, Gerges N, Fiset PO et al (2014) Recurrent somatic mutations in ACVR1 in pediatric midline high-grade astrocytoma. Nat Genet 46:462–466. doi: 10.1038/ng.2950
Gao J, Aksoy BA, Dogrusoz U, Dresdner G, Gross B, Sumer SO et al (2013) Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Sci Signal 6:pl1. doi: 10.1126/scisignal.2004088
Gaubatz S, Imhof A, Dosch R, Werner O, Mitchell P, Buettner R et al (1995) Transcriptional activation by Myc is under negative control by the transcription factor AP-2. EMBO J 14:1508–1519
Guo R, Zheng L, Park JW, Lv R, Chen H, Jiao F et al (2014) BS69/ZMYND11 reads and connects histone H3.3 lysine 36 trimethylation-decorated chromatin to regulated pre-mRNA processing. Mol Cell 56:298–310. doi: 10.1016/j.molcel.2014.08.022
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. doi: 10.1186/s40478-015-0269-0
Johnson BE, Mazor T, Hong C, Barnes M, Aihara K, McLean CY et al (2014) Mutational analysis reveals the origin and therapy-driven evolution of recurrent glioma. Science 343:189–193. doi: 10.1126/science.1239947
Jones C, Karajannis MA, Jones DTW, Kieran MW, Monje M, Baker SJ et al (2017) Pediatric high-grade glioma: biologically and clinically in need of new thinking. Neuro-Oncology 19:153–161. doi: 10.1093/neuonc/now101
Kobayashi S, Boggon TJ, Dayaram T, Janne PA, Kocher O, Meyerson M et al (2005) EGFR mutation and resistance of non-small-cell lung cancer to gefitinib. N Engl J Med 352:786–792. doi: 10.1056/NEJMoa044238
Korbel JO, Campbell PJ (2013) Criteria for inference of chromothripsis in cancer genomes. Cell 152:1226–1236. doi: 10.1016/j.cell.2013.02.023
Le DT, Uram JN, Wang H, Bartlett BR, Kemberling H, Eyring AD et al (2015) PD-1 blockade in tumors with mismatch-repair deficiency. N Engl J Med 372:2509–2520. doi: 10.1056/NEJMoa1500596
Lewis PW, Muller MM, Koletsky MS, Cordero F, Lin S, Banaszynski LA et al (2013) Inhibition of PRC2 activity by a gain-of-function H3 mutation found in pediatric glioblastoma. Science 340:857–861. doi: 10.1126/science.1232245
Liu X, Wang L, Zhao K, Thompson PR, Hwang Y, Marmorstein R et al (2008) The structural basis of protein acetylation by the p300/CBP transcriptional coactivator. Nature 451:846–850. doi: 10.1038/nature06546
Liu XY, Gerges N, Korshunov A, Sabha N, Khuong-Quang DA, Fontebasso AM et al (2012) Frequent ATRX mutations and loss of expression in adult diffuse astrocytic tumors carrying IDH1/IDH2 and TP53 mutations. Acta Neuropathol 124:615–625. doi: 10.1007/s00401-012-1031-3
Lynch TJ, Bell DW, Sordella R, Gurubhagavatula S, Okimoto RA, Brannigan BW et al (2004) Activating mutations in the epidermal growth factor receptor underlying responsiveness of non-small-cell lung cancer to gefitinib. N Engl J Med 350:2129–2139. doi: 10.1056/NEJMoa040938
Morrissy AS, Garzia L, Shih DJ, Zuyderduyn S, Huang X, Skowron P et al (2016) Divergent clonal selection dominates medulloblastoma at recurrence. Nature 529:351–357. doi: 10.1038/nature16478
Nadaf J, Majewski J, Fahiminiya S (2015) ExomeAI: detection of recurrent allelic imbalance in tumors using whole-exome sequencing data. Bioinformatics 31:429–431. doi: 10.1093/bioinformatics/btu665
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. doi: 10.1038/ncomms11185
Northcott PA, Buchhalter I, Morrissy AS, Hovestadt V, Weischenfeldt J, Ehrenberger T et al (2017) The whole-genome landscape of medulloblastoma subtypes. Nature 547:311–317. doi: 10.1038/nature22973
Papillon-Cavanagh S, Lu C, Gayden T, Mikael LG, Bechet D, Karamboulas C et al (2017) Impaired H3K36 methylation defines a subset of head and neck squamous cell carcinomas. Nat Genet 49:180–185. doi: 10.1038/ng.3757
Pasini D, Malatesta M, Jung HR, Walfridsson J, Willer A, Olsson L et al (2010) Characterization of an antagonistic switch between histone H3 lysine 27 methylation and acetylation in the transcriptional regulation of Polycomb group target genes. Nucleic Acids Res 38:4958–4969. doi: 10.1093/nar/gkq244
Pasqualucci L, Dominguez-Sola D, Chiarenza A, Fabbri G, Grunn A, Trifonov V et al (2011) Inactivating mutations of acetyltransferase genes in B-cell lymphoma. Nature 471:189–195. doi: 10.1038/nature09730
Paugh BS, Qu C, Jones C, Liu Z, Adamowicz-Brice M, Zhang J et al (2010) Integrated molecular genetic profiling of pediatric high-grade gliomas reveals key differences with the adult disease. J Clin Oncol 28:3061–3068. doi: 10.1200/JCO.2009.26.7252
Ramaswamy V, Remke M, Bouffet E, Faria CC, Perreault S, Cho YJ et al (2013) Recurrence patterns across medulloblastoma subgroups: an integrated clinical and molecular analysis. Lancet Oncol 14:1200–1207. doi: 10.1016/S1470-2045(13)70449-2
Schwartzentruber J, Korshunov A, Liu X-Y, Jones DTW, Pfaff E, Jacob K et al (2012) Driver mutations in histone H3.3 and chromatin remodelling genes in paediatric glioblastoma. Nature 482:226–231. doi: 10.1038/nature10833
Shi Y, Majewski J (2013) FishingCNV: a graphical software package for detecting rare copy number variations in exome-sequencing data. Bioinformatics 29:1461–1462. doi: 10.1093/bioinformatics/btt151
Siravegna G, Mussolin B, Buscarino M, Corti G, Cassingena A, Crisafulli G et al (2015) Clonal evolution and resistance to EGFR blockade in the blood of colorectal cancer patients. Nat Med 21:795–801. doi: 10.1038/nm.3870
Staedtke V, Bai RY, Laterra J (2016) Investigational new drugs for brain cancer. Expert Opin Investig Drugs 25:937–956. doi: 10.1080/13543784.2016.1182497
Sturm D, Witt H, Hovestadt V, Khuong-Quang D-A, Jones David TW, Konermann C et al (2012) Hotspot mutations in H3F3A and IDH1 define distinct epigenetic and biological subgroups of Glioblastoma. Cancer Cell 22:425–437. doi: 10.1016/j.ccr.2012.08.024
Taylor KR, Mackay A, Truffaux N, Butterfield YS, Morozova O, Philippe C et al (2014) Recurrent activating ACVR1 mutations in diffuse intrinsic pontine glioma. Nat Genet 46:457–461. doi: 10.1038/ng.2925
Van der Auwera GA, Carneiro MO, Hartl C, Poplin R, Del Angel G, Levy-Moonshine A et al (2013) From FastQ data to high confidence variant calls: the genome analysis toolkit best practices pipeline. Curr Protoc Bioinformatics 43:11.10.1–11.1033. doi: 10.1002/0471250953.bi1110s43
Viana-Pereira M, Lee A, Popov S, Bax DA, Al-Sarraj S, Bridges LR et al (2011) Microsatellite instability in pediatric high grade glioma is associated with genomic profile and differential target gene inactivation. PLoS One 6:e20588. doi: 10.1371/journal.pone.0020588
Wang J, Cazzato E, Ladewig E, Frattini V, Rosenbloom DI, Zairis S et al (2016) Clonal evolution of glioblastoma under therapy. Nat Genet 48:768–776. doi: 10.1038/ng.3590
Wang K, Li M, Hakonarson H (2010) ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res 38:e164. doi: 10.1093/nar/gkq603
Wang X, Dubuc AM, Ramaswamy V, Mack S, Gendoo DM, Remke M et al (2015) Medulloblastoma subgroups remain stable across primary and metastatic compartments. Acta Neuropathol 129:449–457. doi: 10.1007/s00401-015-1389-0
Wang Y, Waters J, Leung ML, Unruh A, Roh W, Shi X et al (2014) Clonal evolution in breast cancer revealed by single nucleus genome sequencing. Nature 512:155–160. doi: 10.1038/nature13600
Wen H, Li Y, Xi Y, Jiang S, Stratton S, Peng D et al (2014) ZMYND11 links histone H3.3K36me3 to transcription elongation and tumour suppression. Nature 508:263–268. doi: 10.1038/nature13045
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. doi: 10.1038/ng.1102
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. doi: 10.1038/ng.2938
Yan H, Parsons DW, Jin G, McLendon R, Rasheed BA, Yuan W et al (2009) IDH1 and IDH2 mutations in gliomas. N Engl J Med 360:765–773. doi: 10.1056/NEJMoa0808710
Yuan L, Chi Y, Chen W, Chen X, Wei P, Sheng W et al (2015) Immunohistochemistry and microsatellite instability analysis in molecular subtyping of colorectal carcinoma based on mismatch repair competency. Int J Clin Exp Med 8:20988–21000
Zhang CZ, Leibowitz ML, Pellman D (2013) Chromothripsis and beyond: rapid genome evolution from complex chromosomal rearrangements. Genes Dev 27:2513–2530. doi: 10.1101/gad.229559.113
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. doi: 10.1038/ng.2995
We would like to thank all the patients and families who have supported this project. We would also like to thank Dr. Blake Gilks for his expert consultation of the immunohistochemical interpretation. This research was supported by: Cancer Free KIDS (RS), US National Institutes of Health (NIH) grant P01CA196539 (NJ, JM), the Canadian Institutes for Health Research (CIHR) grant MOP 286756 to NJ; the Fonds de recherche du Québec Santé (FRQS) to CK. JM holds a Canada Research Chair (tier 2). Computational infrastructure was provided by Compute Canada. This work was performed within the context of the I-CHANGE consortium and supported by funding from Genome Canada, Genome Quebec, the Institute for Cancer Research of the CIHR, McGill University and the Montreal Children’s Hospital Foundation. NJ is a member of the Penny Cole Laboratory and the recipient of a Chercheur Boursier, Chaire de Recherche Award from the FRQS. MKM is funded by a CIHR Banting postdoctoral fellowship, NDJ is a recipient of an FRQS studentship.
Ethics approval and consent to participate
All procedures performed in studies involving human participants were in accordance with the ethical standards of the institutional and/or national research committee and with the 1964 Helsinki declaration and its later amendments or comparable ethical standards.
Consent for publication
Informed consent was obtained from all individual participants included in this study.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Tissue location and type, tumor content, WES coverage and digital droplet PCR results in 16 pairs of pHGG analyzed in this study. (XLSX 13 kb)
Somatic and putative somatic mutation information for 16 pairs of pHGGs analyzed in this study. (XLSX 301 kb)
IGV views a subclonal low frequency PIK3CA mutation in HGG3 from a clinical sequencing panel, WES, and targeted sequencing. (PDF 2380 kb)
Number of Single Nucleotide Variants (SNVs) and regions of Allelic Imbalance (AI) present in tumors as shared, primary only, or recurrence only in the pHGG tumor pairs analyzed in this study. (XLSX 13 kb)
Percentages of SNVs and regions of Allelic Imbalance as shared, primary only and recurrence only. (PDF 908 kb)
Immunohistochemical staining for the MMR panel (MLH1, MSH2, MSH6 and PMS2) in the HGG11 primary tumor. (PDF 23521 kb)
Genome-wide view of copy number variations in HGG9 primary and recurrence tumors calculated from Whole Exome Sequencing data. (PDF 2757 kb)
Chromosomal location of AI segments in 15 pairs of pHGG analyzed in this study. (XLSX 61 kb)
Copy number variation (CNV) segments in primary and recurrence tumors from 8 of 16 pairs of pHGG with matched normal tissue available. (XLSX 71 kb)