Genetic alterations of TP53 and OTX2 indicate increased risk of relapse in WNT medulloblastomas

Cohort parameters

In our cohort of 191 WNT-MBs, 125 patients were children (age 3–15 years), 36 were adolescents (16–20 years), and 30 were adults (≥ 21 years). The median age was 13 years. A female predominance was found (ratio 1.4:1 (Table 1, Online Resource Table 3). The majority of tumor samples had classic histology (n = 181; 94.8%); only nine cases presented with large cell/anaplastic (LC/A) histology (4.7%; Fig. 1, Online Resource Fig. 1 and Online Resource Table 4). Among the 120 patients with clinical data, 32 (26.7%) were clinically high-risk (HR); of these, 16 patients because of metastatic disease and the remaining because of large postoperative residual tumor or LC/A histology (eight patients each). Eighty-four patients received “standard-risk” medulloblastoma therapy with 23.4 Gy CSI with a boost to the posterior fossa or tumor bed of 54 Gy followed by maintenance chemotherapy or hyperfractionated radiotherapy with a CSI-dose of 36 Gy and boosts to posterior fossa and tumor bed according to the HIT-SIOP-PNET4 trial [27]. The remaining patients received “high-risk” medulloblastoma therapies, usually with 35.2 Gy CSI followed by a boost to the posterior fossa and maintenance chemotherapy according to the HIT’91 trial in the maintenance chemotherapy arm [17] or according to the MET-HIT2000-AB4 regimen (Table 1, Online Resource Table 5) [41].

Table 1 Cohort overview–demographical, molecular, and clinical characteristicsFig. 1figure 1

Clinical and molecular features of the WNT medulloblastoma cohort. A total of 191 patient samples with their mutational status (CTNNB1, APC, TP53) and most frequent focal and numerical chromosomal alterations are shown as well as nuclear accumulation of β-catenin and p53 protein. Survival data were available from 120 patients. Missing data are shown in gray. wt, wild-type; mut., mutation; IHC, immunohistochemistry; cnLOH, copy-neutral loss of heterozygosity; M0, no metastatic disease; M + , metastatic disease (M1-3); R0, gross-total resection (< 1.5 cm2); R + , subtotal resection; i17q, isochromosome 17q; Chr., chromosome; WCA, whole chromosomal aberrations; S, i17q and WCA data from SNP6 array instead of Molecular Inversion Probe array

Of the three major hallmarks of WNT-MBs, hotspot mutations in exon 3 of CTNNB1 were most frequent (92.2%; 176 of 191), followed by nuclear accumulation of β-catenin protein (91.5% with ≥ 5% of nuclei; Fig. 2a, b). Monosomy 6 was present in 78.6% of analyzed cases (Fig. 1). The latter hallmark was more frequent in WNT-MBs of children (87.2%) compared with adolescents (70.4%) or adults (50%).

Fig. 2figure 2

a Distribution of the 176 CTNNB1 mutations in exon 3 at or adjacent to the 4 phosphorylation sites (purple color), including 4 short in-frame deletions. b Immunohistochemical staining of nuclear β-catenin accumulation in a WNT-MB (ID149) harboring two heterozygous APC mutations; normal cells with blue nuclei serve as internal control; scalebar = 20 µM. c, d Copy-number (upper) and allele ratio (lower) plots from molecular inversion probe array of a sibling pair (ID183 (c), ID154 (d)); the copy-neutral loss of heterozygosity within chromosome arm 5q in ID183 (c) is enlarged and the APC locus is highlighted by the gray vertical line; of note, ID154 (d) has additional mosaic losses of chromosomes 10, 12 and Y. e Distribution of APC mutations in 13 CTNNB1 wild-type WNT-MBs; the homozygous R414C is the only missense mutation and likely benign according to the ACMG/AMP classification rules, but this case (ID151) has a second homozygous stop mutation (R283*) and a copy-neutral loss of heterozygosity, so that the R414C mutation will not be translated to protein and has no tumorigenic impact anyway. One CTNNB1-mutant WNT-MB (ID122) has a heterozygous R2226* APC mutation (not shown), which is likely pathogenic. f, g Kaplan–Meier progression-free (PFS) and overall survival (OS) plots of patients with CTNNB1-mutant versus APC-mutant WNT-MB; the green line shows ID191 without CTNNB1 or APC mutation, but with a WNT-activating homozygous loss of FBXW7. The color-filled regions of the Kaplan–Meier plots indicate pointwise 95%-confidence intervals for the Kaplan–Meier estimate. The numbers at risk refer to the complete cohort. p(adj), adjusted p value after correction for multiple testing (Bonferroni–Holm)

Of the 15 CTNNB1-wild-type tumors, 13 carried APC mutations. Of the remaining two tumor samples, one (ID191) had a homozygous deletion of the FBXW7 gene located on chromosome 4q31.3 (Online Resource Fig. 2), whereas from the other sufficient DNA for APC sequencing was not available. However, this case qualified for the diagnosis of WNT-MB because of the presence of both nuclear β-catenin accumulation and monosomy 6. 450k/850k methylation array data was available in 94 cases, confirming the assignment to the WNT-MB group in 87 cases; in 7 cases, no matching score was reached (Fig. 1).

TP53 mutations were found in 16.1% (30/186 cases; Table 1). Twenty six were missense mutations located in exons 5–8; only one was a missense variant in exon 4. Additionally, we found one splice site mutation and one frameshift deletion (both exon 6), and one case had a focal homozygous deletion of the TP53 locus (chromosome 17p13.1; Fig. 4a, Online Resource Fig. 3). An accumulation of p53 protein in ≥ 5% of tumor cells was present in 20.4% of cases (38 of 186; Figs. 1, 4b). TP53 mutations were found in 18 of 38 p53-positive tumor samples. The TP53 locus showed loss of heterozygosity in 18.5% (24/130 cases), in some only present in a fraction of tumor cells (mosaicism). Of the 30 TP53-mutated cases, allele status was available from 27 samples, and 11 of these had a copy loss on chromosome 17p13.1. A significant predominance of male patients was found in cases with TP53 mutations (18 of 30 (60%), p = 0.020, Chi-square test). TP53 mutations occurred in all age groups (6–23 years) without a significant association with one of the age groups. TP53 mutant and TP53 wild-type cases did not cluster separately in UMAP analysis (Online Resource Fig. 4). LC/A histology was significantly enriched in TP53-mutated WNT-MBs (4 of 30 versus 5 of 156, p = 0.020, Chi-square test).

The E17K hotspot regions of AKT1-3 were Sanger sequenced in 165 cases, and the E17K mutation was found 6 times each in AKT1 and AKT3 (12 of 165; 7.27%).

Next-generation panel sequencing (NGS panel II) was performed in 31 cases. Besides identification of APC variants and validation of CTNNB1 and TP53 mutations, recurrent mutations (n > 2) were found in DDX3X (n = 12, 38.7%), SMARCA4 (n = 5), KMT2D (n = 5), KMT2C (n = 4), PIK3CA (n = 3), and FBXW7 (n = 3; all mutations are shown in Online Resource Table 6).

Chromosomal aberrations

Except for monosomy 6, most WNT-MBs presented stable genomes. Genome-wide copy-number profiles were available from 128 cases with a mean of 2.0 whole chromosomal aberrations (WCAs) per case. Most frequent WCAs besides chromosome 6 were gains of chromosomes 8 and 19 in 16 samples each (12.5%), and losses of chromosomes 10 and 13 in 10 samples each (7.8%; Fig. 4e, Online Resource Fig. 5). TP53-mutated tumors were enriched for monosomy 6 (28/28 cases versus 74/101 TP53 wild-type tumors, p = 0.002, Chi-square test) and chromosome arm 17p loss (11/27 versus 13/88 TP53 wild-type tumors, p = 0.001, Chi-square test; Online Resource Fig. 6). An isochromosome 17q (i17q) was found in 8 of 130 samples (6.2%). Genomic identification of significant targets in cancer (GISTIC) analysis revealed several significantly altered focal aberrations, most notably CDK6 gain (chromosome 7q21.2), OTX2 gain (14q22.3), and AXIN2 gain (17q24.1; Online Resource Table 7 and Online Resource Fig. 7). Of note, OTX2 gain (n = 42; Table 1, Fig. 4f and Online Resource Fig. 8) was the most frequent focal aberration occurring both in TP53 mutant (8/27, 29.6%) and TP53 wild-type (33/79, 41.8%, p = 0.263, chi-square test) tumors. Cases with OTX2 gain versus balanced OTX2 locus did not cluster separately in UMAP analysis (Online Resource Fig. 4).

APC mutations

Of 13 APC-mutant cases, 8 had larger regions with copy-neutral loss of heterozygosity (cnLOH) spanning the APC locus up the telomeric region of chromosome 5q; one case had a focal loss at the APC locus (5q22.2). Eight of these nine cases with deletion carried a single mutation, and one carried a double mutation on the same allele (R283* and R414C) with allele variant fractions between 75 and 100%. The remaining four cases presented with two mutations with allele variant fractions of both mutations around 40–50% and 3 lacked LOH of the APC locus (1 × not done; Figs. 1, 2e). TP53 mutations and OTX2 gains both occurred at similar frequencies in APC-mutant WNT-MBs compared to CTNNB1-mutant tumors (TP53 mutation: 3/13 versus 27/172, p = 0.486; OTX2 gain: 4/12 versus 37/94, p = 0.686, Chi-square tests).

Occurrence of APC-mutant WNT-MB in siblings

Our WNT-MB cohort contained one sibling pair with a female patient, who developed a classic MB at the age of 11 with metastatic spread at diagnosis (ID183). Her brother was 19 when the tumor was detected, and histology was also classic, but without metastatic disease (ID154). Both were CTNNB1 wild-type, but APC sequencing revealed the same R213* mutation in both tumors, suggesting that it represented a germline mutation. The sister had an allele variant fraction of almost 100% and a large cnLOH on chromosome arm 5q, whereas the brother had a balanced APC locus, but an additional two base-pair frameshift deletion at amino acid L1129. Both APC mutations in this tumor had allele variant fractions of 40–50%.

Additionally, both siblings had no OTX2 gains, but were affected by TP53 alterations. The sister presented a monoallelic TP53 mutation (R282P) without loss of the TP53 locus, whereas the brother had a chromosome arm 17p loss, but wild-type TP53 (Fig. 2c, d). Both siblings are alive without relapse or secondary tumor.

Survival analyses

Staging and follow-up data from 120 patients were available with median follow-up of 7.4 years (range 0.1–19.6 years) in 109 patients alive at last follow-up. 5-year PFS was 88% [95% CI 82–95%] and 5-year OS was 93% [88–98%]. Univariable survival analysis did not show statistically significant associations of survival with gender and clinical HR stage (M + , residual disease; Fig. 3a–d), LC/A histology, nuclear β-catenin accumulation, i17q, AKT1-3 mutations, CDK6 gains, or monosomy 6 as well as clinical risk and therapy regimen (Online Resource Figs. 9–11). As in other studies, patients older than 16 had worse survival, reaching significance in OS (5-years OS 3–15 years 96% [92–100%] versus 16–20 years 80% [62–100%] versus ≥ 21 years 86% [63–100%], p = 0.030), but not in PFS (91% [85–98%] versus 79% [63–100%] versus 86% [63–100%], p = 0.286). After correction for multiple testing, the difference in age groups for OS did not reach significance anymore (p(adj) = 0.178; Fig. 3e, f). APC-mutated cases showed similar survival as CTNNB1-mutated cases (Fig. 2f, g). In contrast, TP53 mutations were associated with lower PFS (68% [49–93%] versus 93% [87–99%], p = 0.001, p(adj) = 0.003) and OS (88% [83–100%] versus 94% [89–99%], p = 0.105, p(adj) = 0.105; Fig. 4c, d), and TP53 copy loss also was associated with lower PFS (p = 0.004, p(adj) = 0.014), but not OS (Online Resource Fig. 12), as was loss of whole chromosome 10: PFS p < 0.001, p(adj) = 0.003; loss of chromosome 13: PFS p = 0.003, p(adj) = 0.017 (Online Resource Figs. 13–15). Importantly, gains of the OTX2 locus were associated with significant worse PFS (5-year PFS 72% [57–91%] versus 93% [85–100%], p = 0.017, p(adj) = 0.034) and OS (5-year OS 82% [69–98%] versus 97% [91–100%], p = 0.006, p(adj) = 0.019; Fig. 4g, h), which was independent and additive of the effect of TP53 mutations on univariable levels. No relapses or other events were observed among patients with tumors harboring neither TP53 mutations nor OTX2 gains (5-year PFS/OS 100%; Online Resource Table 8), whereas 5-year PFS/OS was 78% 65–93%/85% 73–98% in patients whose tumors either had a TP53 mutation or an OTX2 gain (n = 36, 48.7%) and was 50% [22–100]/83% [58–100%] among patients with both TP53 mutation and OTX2 gain (n = 7, p < 0.001, p(adj) = 0.001/p = 0.001, p(adj) = 0.010; Fig. 5a–c, Online Resource Fig. 16). The difference in patients whose tumors either had a TP53 mutation or an OTX2 gain for PFS seemed to be relatively stable under subsampling (p value ranging from 0 to 0.036 in 100 subsampling data sets; Online Resource Table 9).

Fig. 3figure 3

Kaplan–Meier progression-free (PFS) and overall survival (OS) plots for gender (a, b), staging with clinical high-risk factors “metastatic disease” (M +) and “subtotal resected tumor” (R + ; ≥ 1.5 cm2) versus cases with non-metastatic disease and gross-total resection (M0R0; c, d), and age (e, f). e, f Log-rank p values for individual PFS analyses are: blue versus red = 0.142; red versus green = 0.894; blue versus green = 0.284. Log-rank p values for individual OS analyses are: blue versus red = 0.020; red versus green = 0.856; blue versus green = 0.036

Fig. 4figure 4

a Distribution of TP53 mutations; colored mutations represent splice site mutation (blue), frameshift deletion (green), and nonsense mutation (red). b Immunohistochemical staining for p53 protein in a WNT-MB (ID184) with a heterozygous R248Q TP53 mutation and balanced TP53 locus; a strong nuclear staining pattern is visible focally, whereas some regions lack p53 staining (left side of the image); scalebar = 50 µM. c, d Kaplan–Meier progression-free (PFS) and overall survival (OS) plots of patients with TP53 mutant versus TP53 wild-type WNT-MB. e Summary plot of copy-number aberrations (n = 115 WNT-MBs) from molecular inversion probe array (n = 108) and SNP6 array (n = 7). The arrow indicates the OTX2 locus on chromosome arm 14q. Blue bars, gains; red bars, losses. Thickness of bars indicates frequency of alterations. f Copy-number plots from molecular inversion probe arrays of a WNT-MB case with focal OTX2 gain (upper part) and balanced OTX2 (lower part) on chromosome arm 14q. g, h Kaplan–Meier PFS and OS plots of patients with OTX2 gain versus no OTX2 gain

Fig. 5figure 5

a, b Kaplan–Meier progression-free (PFS) and overall survival (OS) plots of patients with WNT-MBs without TP53 mutation and OTX2 gain (blue line), one of these two alterations (red line), or both alterations (green line). Log-rank p values for individual PFS analyses are: blue versus red = 0.010; red versus green = 0.050; blue versus green < 0.001. Log-rank p values for individual OS analyses are: blue versus red = 0.030; red versus green = 0.100; blue versus green < 0.001. c Pie chart showing frequency of identified risk factors in 74 WNT-MB cases with PFS/OS and genomic data from Molecular Inversion Probe array. d Venn diagram showing the distribution of TP53 mutant, OTX2 gained, and chromosome 6 lost (n = 64 with at least one of these alterations; n = 10 were wt/balanced for all 3 factors). e, f Results of the multivariable Cox regression analyses for PFS and OS; *, p < 0.05; **, p < 0.01. g Results of the subsampling approach for the multivariable Cox regression analysis with variable selection. Values in brackets refer to the IQR or the range for the subsampling approach, respectively. Values outside brackets refer to the model on the complete data (n = 74). PFS: Concordance Index [range]: 0.8 [0.76; 0.86], Median Rank: 8, Frequency Rank 1: 14%, Average AIC: 65.88; OS: Concordance Index [range]: 0.82 [0.75; 0.91], Median Rank: 3, Frequency Rank 1: 35%, Average AIC: 49.92

This effect was retained in a multivariable Cox regression analysis for PFS and OS, where OTX2 gain (HR 5.98 [1.63–22.0], p = 0.007/12.9 [1.6–103], p = 0.016), TP53 mutations (HR 10.8 [2.2–53.2], p = 0.003/6.2 [1.1–34.5], p = 0.037), and monosomy 6 (HR 0.15 [0.03–0.85], p = 0.032/0.15 [0.027–0.82], p = 0.028) were identified as significant and independent prognostic markers (Fig. 5d–f). The model selection procedure applied to the 100 subsampling data sets revealed the model containing the three mentioned covariates to have the minimal average AIC as well as the minimal median rank when ranking models according to the AIC (PFS: average AIC = 65.88 OS: average AIC = 49.92; Fig. 5g and Online Resource Fig. 17).

留言 (0)

沒有登入
gif