Genetic and epigenetic landscape of IDH-wildtype glioblastomas with FGFR3-TACC3 fusions

A subset of glioblastomas (GBMs) harbors potentially druggable oncogenic FGFR3-TACC3 (F3T3) fusions. However, their associated molecular and clinical features are poorly understood. Here we analyze the frequency of F3T3-fusion positivity, its associated genetic and methylation profiles, and its impact on survival in 906 IDH-wildtype GBM patients. We establish an F3T3 prevalence of 4.1% and delineate its associations with cancer signaling pathway alterations. F3T3-positive GBMs had lower tumor mutational and copy-number alteration burdens than F3T3-wildtype GBMs. Although F3T3 fusions were predominantly mutually exclusive with other oncogenic RTK pathway alterations, they did rarely co-occur with EGFR amplification. They were less likely to harbor TP53 alterations. By methylation profiling, they were more likely to be assigned the mesenchymal or RTK II subclass. Despite being older at diagnosis and having similar frequencies of MGMT promoter hypermethylation, patients with F3T3-positive GBMs lived about 8 months longer than those with F3T3-wildtype tumors. While consistent with IDH-wildtype GBM, F3T3-positive GBMs exhibit distinct biological features, underscoring the importance of pursuing molecular studies prior to clinical trial enrollment and targeted treatment. Electronic supplementary material The online version of this article (10.1186/s40478-020-01058-6) contains supplementary material, which is available to authorized users.


Introduction
The molecular landscape of IDH-wildtype glioblastoma (GBM) has been extensively characterized, yet it remains a disease with a dismal prognosis [8,9,21]. Oncogenic fusions have recently been recognized as molecular drivers in a subset of IDH-wildtype GBMs. Approximately 3% of IDH-wildtype GBMs have been reported to harbor activating fusions involving the tyrosine kinase domain (TKD) of the fibroblast growth factor receptor 3 (FGFR3) gene and the coiled-coil domain of the transforming acidic coiled-coil-containing protein 3 (TACC3) gene [39]. The FGFR3-TACC3 (F3T3) fusion protein is thought to promote malignant transformation by increasing downstream signaling through the MAPK pathway, activating mitochondrial biogenesis and metabolism, and recruiting endogenous TACC3 away from the mitotic spindle, leading to delayed mitotic progression and aneuploidy [20,37].
The relative rarity of F3T3-positive, IDH-wildtype GBM has hampered a full characterization of this molecular subset of GBM. While prior reports have suggested that these tumors may have recurrent histologic features and molecular profiles, including the absence of EGFR amplification and an increased frequency of CDK4 and MDM2 amplifications, extensive genomic characterization is lacking, particularly in the context of long-term clinical follow up and survival [7]. Further, since the discovery of F3T3 in GBM, DNA methylation-based tumor classification has emerged as a promising modality for improving diagnostic precision in neuropathology [10]. However, methylation profiles specific to F3T3-positive GBMs have not been reported. Further study of whether F3T3 fusions are associated with other genetic or epigenetic alterations may refine tumor subclassification efforts and impact prognosis and treatment.
The objective of this study was to describe the genomic landscape and methylation profiles of an unbiased institutional cohort of patients with F3T3-positive, IDH-wildtype GBMs and to identify potential associations with patient clinicopathologic characteristics and survival.

Design, setting, and participants
We performed a retrospective cohort study of all patients diagnosed with IDH-wildtype GBMs between January 2015 and December 2019 at Memorial Sloan Kettering Cancer Center (MSK) in New York, NY. The study was approved by the MSK Institutional Review Board and reported in accordance with the STROBE (Strengthening the Reporting of Observational Studies in Epidemiology) guidelines [18].

Eligibility criteria
All patients with F3T3-positive GBMs were included in the study. All patients with F3T3-wildtype GBMs from the same period were selected for comparison. Cases with IDH1 or IDH2 mutations curated as oncogenic by the MSK OncoKB Precision Oncology Knowledge Base were excluded [14].

Diagnostic criteria
Cases were diagnosed by board-certified MSK neuropathologists (M.K.R., T.A.B.) according to the diagnostic criteria specified in the 2016 World Health Organization Classification of Tumors of the Central Nervous System [28].

Systematic review and meta-analysis of FGFR3-TACC3 fusion prevalence
The prevalence of F3T3 positivity identified among the IDH-wildtype GBMs by MSK-IMPACT (Integrated Mutation Profiling of Actionable Cancer Targets) in this study was synthesized with previously published estimates using random effects meta-analysis according to a previously reported method [31,36].

DNA-based molecular analyses
To assess for the presence of FGFR3 fusions and other molecular alterations, tumors and matched normal blood samples were analyzed with the MSK-IMPACT nextgeneration DNA sequencing platform that targets up to 468 genes and select introns to produce data on single nucleotide variants, small insertions and deletions, copy number variation, and structural variants [16,42]. Tumor mutational burden (TMB) was defined as the number of mutations per megabase (mt/Mb).

Methylation-based molecular analyses
A subset of cases was analyzed with the Infinium Meth-ylationEPIC (850K) platform, which provides MGMT promoter hypermethylation status as well as data on > 850,000 CpG methylation sites across the genome [3,6]. Cases were assigned methylation-based classes using version 0.1.124 of the random forest-based mnp.v11b4 R package obtained from the German Cancer Research Center (DKFZ) [10]. To investigate possible differences in methylation profiles between F3T3-positive and F3T3wildtype GBMs, dimensionality reduction with principal component analysis and t-distributed stochastic neighbor embedding was performed and overlaid on a reference cohort of GBMs from the previously published DKFZ study using version 0.15 of the Rtsne package with the following non-default parameters: initial_dims = 100, max_iter = 1500, and theta = 0. Cross-reactive, sex-chromosome, and failed probes were excluded from these analyses.

Copy number alteration-based cytogenetic analyses
In addition to providing data on methylation profiles, the Infinium platform provides high resolution copy-number variant data similar to that provided by conventional whole-genome copy-number microarray. The chromosome-and arm-level cytogenetic profiles of the GBMs with F3T3 fusions were assessed using version 1.20.0 of the conumee-based MNPcnvplot function in the mnp. v11b4 package in R [26]. The overall tumor copy-number alteration burden (TCB), defined as the percentage of the analyzed genome for which the copy number was not equal to two, was also assessed using the larger MSK-IMPACT dataset. Tumor ploidy was also inferred from the tumor and matched germline MSK-IMPACT sequencing data using the previously validated FACETS algorithm [38].

Clinical data and survival analyses
Clinical charts were reviewed to extract data on age, sex, date of initial pathologic diagnosis, survival time, history of radiotherapy with or without concurrent temozolomide (TMZ), and treatment history. Reported MGMT promoter hypermethylation status, typically determined by pyrosequencing, was extracted from the electronic medical record for those patients who did not undergo methylation analysis on the Infinium platform [32]. Overall survival was defined as the time from diagnostic biopsy until the time of death due to any cause.

Statistical analyses
Differences in sample means were assessed using the unpaired Student's t test. Differences in sample medians were assessed using the unpaired Mann-Whitney-Wilcoxon test. Differences among categorical variables were assessed using Fisher's exact test. Binomial proportion confidence intervals were calculated using the Clopper-Pearson method. Clinical and pathological variables were examined in univariate and multivariate Cox proportional hazards models for associations with overall survival. The multivariate models were adjusted for potential confounders including age, sex, race, and MGMT status. Genomic data were accessed using the internal MSK cBioPortal for Cancer Genomics and statistical analyses were performed using R version 3.6.2 (R Foundation for Statistical Computing) [13,22,35]. Statistical tests were 2-sided and used a significance threshold of p < 0.05. Reported p values were not adjusted for multiple testing.

Description of the cohort
Between January 2015 and December 2019, samples from 906 patients with IDH-wildtype GBMs underwent DNA sequencing with MSK-IMPACT (Fig. 1 1 Flow diagram of specimen selection. Note: The two tumors with non-canonical fusions were excluded from further analyses because dimerization domains have not been reported in these 3′ partner genes [44] vs. 39.6% [343/867] female, p = 0.96). There was no difference in laterality or anatomic localization (p > 0.05 for both comparisons). All tumors were reviewed by a boardcertified neuropathologist and found to meet histologic criteria for GBM, WHO grade IV (i.e., highly cellular and variably anaplastic glial cells with brisk mitotic activity, microvascular proliferation, and/or necrosis) [28].

Prevalence of FGFR3-TACC3 fusions compared to previously published studies of IDH-wildtype GBMs
The 4.1% (37/906) prevalence of F3T3 fusions identified in this study was compared to prevalence estimates reported in six previously published studies [2,5,17,33,34,39]. Taken together, the previously published studies involved a total of 883 IDH-wildtype GBMs (median, 68.5; range 17-584). Meta-analytic pooling of the prevalence value from the present study with those from the six previously published studies yielded a summary F3T3-fusion prevalence of 3.7% (95% confidence interval [CI], 2.9-4.7%), with no evidence of between-study heterogeneity (Q = 8.8, τ 2 = 0, I 2 = 0%, p = 0.18) (Fig. 2). The prevalence estimates reported by the individual studies ranged from 1.3 to 11.8%. Sensitivity analysis, in which the meta-analysis was serially repeated after exclusion of each study, demonstrated that no individual study affected the overall estimate by more than 0.4%. Thus, we conclude that the prevalence of F3T3 fusions in IDHwildtype GBMs is approximately 4%.

Associations between FGFR3-TACC3 fusions and canonical drivers of glioblastoma tumorigenesis
The relative frequencies of selected mutations and copynumber alterations curated as potentially oncogenic in OncoKB [14] and/or as putative driver alterations in the Cancer Hotspots [15] database in the RTK/RAS/MAPK, PI3K, P53, cell cycle, and telomere maintenance pathways among the 37 F3T3-positive GBMs are presented in Fig. 3 and Additional file 1: Fig. S2. F3T3-positive tumors had significantly lower median TMB than F3T3-wildtype  . 4a).

Copy number alteration-based cytogenetic analysis of glioblastomas with and without FGFR3-TACC3 fusions
Among F3T3-positive tumors with analyzable data, 78.1% (25/32) had concurrent chromosome 7 gain and 10q loss, 9.4% (3/32) had chromosome 7 gain without 10q loss, 6.3% (2/32) had 10q loss without chromosome 7 gain, and 6.3% (2/32) had concurrent loss of 10q, 13q, and 14q, all consistent with the integrated diagnosis of IDH-wildtype GBM. The overall TCB was assessed using the larger MSK-IMPACT dataset that included all 904 patients. Because the F3T3 fusion protein is thought to cause aneuploidy, we hypothesized that F3T3-positive GBMs would have higher TCBs than F3T3-wildtype GBMs. Unexpectedly, F3T3-positive tumors had significantly lower TCBs. The median TCB for the F3T3-positive tumors was 15.0% (IQR, 10.4-17.5%) while that for the F3T3-wildtype tumors was 17.4% (IQR, 12.9-23.4%) (p = 0.006) (Fig. 4b). A similar finding was obtained using a subset (n = 371) of the MSK-IMPACT paired germline Fig. 6 Dimensionality reduction demonstrated that FGFR3-TACC3 fusion-positive tumors were more likely to be assigned to the mesenchymal and RTK II methylation subclasses. Note: Cases were assigned DNA methylation-based classifications according to the method described by Capper et al. [10]. F3T3-positive GBMs overwhelmingly fell into the GBM, MES and GBM, RTK II clusters data in an allele-specific copy-number analysis of tumor ploidy using the FACETS algorithm [38], which revealed that F3T3-positive tumors had slightly lower median ploidy (2. (Fig. 4c).

Treatment of patients with FGFR3-TACC3 fusions
Treatment data were available for 91.9% (34/37) of the patients with F3T3-positive GBMs, all of whom received upfront radiotherapy and TMZ (concurrently and/or as adjuvant therapy). Only a minority of patients received a therapy targeting FGF/FGFR signaling. Specifically, 5.9% (2/34) received a selective FGFR inhibitor and 5.9% (2/34) received a multikinase inhibitor/non-selective FGFR inhibitor. The small number of patients receiving targeted therapy precluded investigating its association with overall survival.

Discussion
This genomic landscape study of 906 patients with IDHwildtype GBMs treated at a tertiary cancer referral center demonstrated that 4.1% harbored potentially druggable activating F3T3 fusions with recurrent structural isoforms; characteristic associated mutational, copynumber, and methylation profiles; and clinical outcomes slightly better than patients with F3T3-wildtype tumors despite occurring in older individuals.
Our understanding of F3T3-positive GBMs has increased since the initial discovery that demonstrated F3T3 fusions comprised a rare yet recurrent oncogenic structural event resulting in an in-frame fusion protein with constitutive kinase activity and mitogenic effects, which were reversible by FGFR kinase inhibition [34,39]. Together with subsequent studies involving a total 883 IDH-wildtype GBMs, it was previously estimated that a b Fig. 7 Kaplan-Meier analysis of overall survival stratified by FGFR3-TACC3 fusion (a) and MGMT promoter hypermethylation (b) status. Note: Overall survival was defined as the time from initial diagnostic biopsy until the time of death due to any cause. Two-sided p values for difference were calculated using the log-rank test 1-12% might harbor F3T3 fusions, with 3% a commonly cited prevalence estimate [2,5,17,33,34,39]. The present study, involving a large, unbiased single-institutional cohort of 906 prospectively sequenced GBMs, revealed an F3T3-fusion prevalence of 4.1% among individuals with IDH-wildtype GBMs. Meta-analytic pooling suggested a similar prevalence, with an upper confidence interval indicating that as many as 5% of GBMs may have F3T3 fusions, slightly higher than suggested by prior consensus. Given ongoing interest in targeted FGFR3 inhibition, these data support the importance of pursuing comprehensive molecular characterization with robust fusion detection on all IDH-wildtype GBMs to identify patients for potential clinical trial enrollment and targeted treatment. The present study confirms molecular features of F3T3-positive GBMs suggested in prior reports, including the diversity of structural isoforms of the F3T3 fusion gene itself. Prior studies established that isoforms involving FGFR3 exons 17 or 18 were most common [17,34,39], in agreement with our finding that 97.3% (36/37) of cases included one of these exons. As in prior studies, we found great variability for the 3′ TACC3 exon, ranging from exons 3-12. Overall, FGFR3 exon 17 to TACC3 exon 11 and FGFR3 exon 17 to TACC3 exon 10 fusions were most common.
However, our data reveal new insights into the molecular alterations associated with F3T3 fusions in IDHwildtype GBMs, in some instances contrasting with prior smaller studies that employed more targeted sequencing methods. For example, prior reports suggested that the F3T3 fusion was mutually exclusive with EGFR, PDG-FRA, or MET amplification [23,34], and had an increased frequency of co-occurring CDK4 and MDM2 amplification [7,17]. In the present study, F3T3-positivity was indeed mutually exclusive with amplification in PDGFRA or MET. Similarly, we found that F3T3-positive GBMs were less likely to exhibit concurrent EGFR amplification (5.6% vs. 43.3%); however, our data revealed that these events were not truly mutually exclusive, as two cases with F3T3 fusions exhibited bona fide amplification in EGFR. However, it is possible that the F3T3 fusion and the EGFR amplification events were present in different subclones, as has been shown for amplification of diverse RTKs [40,41]. Further, there was no increased frequency of CDK4 or MDM2 amplification among F3T3-positive GBMs in our study.
The present study utilized a hybridization capturebased NGS assay (MSK-IMPACT) for sequencing all exons and selected introns of up to 468 genes with paired tumor-matched normal tissue analysis enabling unambiguous somatic mutation detection [16,42], allowing us to comprehensively analyze oncogenic pathway alterations.
A key negative finding gleaned from these data was a lack of significant differences in downstream PI3K, cell cycle, and telomere maintenance pathway alterations in F3T3positive and F3T3-wildtype GBMs. Also interesting was the observation that oncogenic alterations in TP53 and MDM4 were less common among F3T3-positive tumors. Combined with the lower median TMB and TCB observed in F3T3-positive GBMs, these data suggest that F3T3-positivity itself may be enough to drive oncogenesis in the relative absence of other concurrent driver alterations. The relatively genomically quiet nature of the F3T3-positive GBMs was unexpected given that experimental data suggest that the F3T3 fusion protein recruits endogenous TACC3 away from the mitotic spindle, leading to delayed mitotic progression and aneuploidy [37]; our data therefore suggest that additional mechanisms to counteract this activity (e.g., the relative rarity of concurrent TP53 and MDM4 alterations identified in this study) may be engaged in vivo.
DNA methylation-based CNS tumor classification is becoming increasingly employed as a diagnostic tool in many neuropathology practices worldwide. The present study provides insight into the F3T3-positive GBM methylome using a previously validated machine-learningbased tumor classifier developed at University Hospital Heidelberg [10]. Aside from confirming the prior observation that MGMT promoter hypermethylation does not differ between F3T3-positive and F3T3-wildtype GBMs [17], the present analysis demonstrated that the majority of F3T3-positive GBMs can be accurately classified as IDH-wildtype GBMs by the Heidelberg classifier. Further, it showed that F3T3-positive cases overwhelmingly (in 96.7% of cases) fall into the mesenchymal or RTK II subclasses, which overlap in their methylation profiles, rather than the RTK I (or any other) subclass. In contrast, GBMs with EGFR or PDGFRA amplifications were assigned the RTK I subclass in 14.3% and 50.0% of cases, respectively, consistent with prior reports [10]. While it is unclear if the role of TACC3 in regulating the epithelialmesenchymal transition may relate to the methylation characteristics of F3T3-positive GBM [24], it has been suggested that methylation subclass is reflective of tumor cell of origin as well as somatic epigenetic changes [11,19,25]. Our data show that while F3T3-positive GBMs are closely related to other IDH-wildtype GBMs, there are differences among these tumors that may be biologically significant.
The present study is among the first to provide longterm clinical follow-up data on a large cohort of patients with F3T3-positive tumors, with a median follow-up time of 55.6 months. Patients with F3T3-positive tumors were slightly older at initial diagnosis (63.6 vs. 59.4 years). Although older age at diagnosis is a well-established risk factor for poorer outcome and shorter overall survival in patients with GBM, patients with F3T3-positive GBMs in this study lived slightly longer than those without F3T3 fusions (26.7 vs. 18.5 months, p = 0.07), a difference that persisted after adjustment for potential confounders including age, sex, race, and MGMT status. Although this difference did not meet formal statistical significance, at the very least, we can conclude that patients with F3T3positive tumors have clinical features and outcomes similar to, or perhaps slightly better than, those with non-F3T3-driven, IDH-wildtype GBMs.
Of note, F3T3 fusions have also been identified in lower-grade histologic entities [17,23,27,29,43]. Sufficient long-term follow up is lacking to determine if any among these are under-sampled higher-grade gliomas, or gliomas that eventually recur, progress, or otherwise behave aggressively. Furthermore, both low-and high-grade F3T3-positive gliomas exhibit characteristic histologic features, including monomorphous oligodendroglioma-like nuclei, "chicken-wire" capillary networks, and frequent microcalcifications [7]. While these features suggest the attractive hypothesis that F3T3-positive GBMs arise from lower-grade precursor lesions, to date there has been insufficient evidence to support this idea. Rather, our study underscores that F3T3 fusions, when detected in histologic GBM, drive clinical behavior akin to other IDH-wildtype GBMs.
Importantly, given the difference in overall survival suggested by this study, our findings raise the possibility that the underlying biology of F3T3-positive GBM may be a confounding factor in interpreting outcome data in the setting of pan-GBM clinical trials. Considering ongoing research into therapeutic inhibition of tyrosine kinase signaling in general and FGFR3 signaling specifically, identification of F3T3 fusions in patients with IDH-wildtype GBMs could soon have therapeutic implications. For example, the recent phase 2 REGOMA trial showed that median overall survival was approximately 2 months longer in patients with first recurrence of GBM after surgery and radiotherapy/TMZ receiving the oral multikinase inhibitor regorafenib rather than lomustine (7.4 vs. 5.6 months, p < 0.001) [30]. Notably, that study did not enroll patients based on a molecular alteration, nor did it report results stratified by underlying somatic genetic alterations, so its relevance to patients with F3T3-positive GBMs is uncertain.
FGFR signaling has been implicated in a variety of human cancers. FGFR pathway inhibition as a therapeutic strategy remains an area of active investigation, including in ongoing and recently completed clinical trials of FGFR inhibitors in brain tumors (e.g., NCT01975701, NCT028224133, NCT02052778, and NCT01948297). Other inhibitors of FGFR3 and F3T3 have also been tested in preclinical models and in trials in other tumor types, especially urothelial carcinoma and cholangiocarcinoma, for which therapeutics such as erdafitinib, rogaratinib, infigratinib, and the monoclonal antibody vofatamab have shown promise [12]. Recent research has also demonstrated that F3T3positive tumors are characterized by mitochondrial activation, and that these tumors may be particularly susceptible to inhibitors of mitochondrial respiration and oxidative metabolism, highlighting this pathway as another potential therapeutic opportunity [20]. Aside from FGFR3, there is also interest in inhibiting TACC3, including with the potent new TACC3-targeting agent BO-264, which significantly inhibited the growth of cells harboring F3T3 fusions in a preclinical model [1].
This study has limitations. Although 906 cases were sequenced, only 137 underwent methylation array analysis. Also, an insufficient number of F3T3-positive patients in this cohort received F3T3-specific targeted therapy to evaluate the potential efficacy of such an approach, the determination of which awaits future clinical trials. Last, detailed treatment information was not available for all F3T3-wildtype patients, thus chemoradiotherapy treatment status was not adjusted for in the Cox proportional hazards regression models.
In conclusion, approximately 4% of IDH-wildtype GBMs in this study harbored potentially targetable F3T3 fusions. While their clinical and molecular characteristics were largely in keeping with other IDH-wildtype GBMs, they demonstrated characteristic associated mutational, copy-number, and methylation profiles, and patients with F3T3-positive tumors had clinical outcomes slightly better than patients with F3T3-wildtype tumors. F3T3-positivity was predominantly mutually exclusive with other RTK drivers, and F3T3-wildtype tumors were enriched for mutations in TP53. As histologic features lack specificity for identifying F3T3-positive tumors, comprehensive NGS and methylation analysis should be considered for all IDH-wildtype GBMs to identify patients for potential clinical trial enrollment and targeted treatment.