Marcel Kool

Integrative Genomic Analysis of Medulloblastoma Identifies a Molecular Subgroup That Drives Poor Clinical Outcome

PURPOSE Medulloblastomas are heterogeneous tumors that collectively represent the most common malignant brain tumor in children. To understand the molecular characteristics underlying their heterogeneity and to identify whether such characteristics represent risk factors for patients with this disease, we performed an integrated genomic analysis of a large series of primary tumors. PATIENTS AND METHODS We profiled the mRNA transcriptome of 194 medulloblastomas and performed high-density single nucleotide polymorphism array and miRNA analysis on 115 and 98 of these, respectively. Non-negative matrix factorization-based clustering of mRNA expression data was used to identify molecular subgroups of medulloblastoma; DNA copy number, miRNA profiles, and clinical outcomes were analyzed for each. We additionally validated our findings in three previously published independent medulloblastoma data sets. RESULTS Identified are six molecular subgroups of medulloblastoma, each with a unique combination of numerical and structural chromosomal aberrations that globally influence mRNA and miRNA expression. We reveal the relative contribution of each subgroup to clinical outcome as a whole and show that a previously unidentified molecular subgroup, characterized genetically by c-MYC copy number gains and transcriptionally by enrichment of photoreceptor pathways and increased miR-183∼96∼182 expression, is associated with significantly lower rates of event-free and overall survivals. CONCLUSION Our results detail the complex genomic heterogeneity of medulloblastomas and identify a previously unrecognized molecular subgroup with poor clinical outcome for which more effective therapeutic strategies should be developed.

FSTL5 Is a Marker of Poor Prognosis in Non-WNT/Non-SHH Medulloblastoma

PURPOSE Integrated genomics approaches have revealed at least four distinct biologic variants of medulloblastoma: WNT (wingless), SHH (sonic hedgehog), group C, and group D. Because of the remarkable clinical heterogeneity of group D tumors and the dismal prognosis of group C patients, it is vital to identify molecular biomarkers that will allow early and effective treatment stratification in these non-WNT/non-SHH tumors. PATIENTS AND METHODS We combined transcriptome and DNA copy-number analyses for 64 primary medulloblastomas. Bioinformatic tools were used to discover marker genes of molecular variants. Differentially expressed transcripts were evaluated for prognostic value in the screening cohort. The prognostic power of follistatin-like 5 (FSTL5) immunopositivity was tested for 235 nonoverlapping medulloblastoma samples on two independent tissue microarrays. RESULTS Comprehensive analyses of transcriptomic and genetic alterations delineate four distinct variants of medulloblastoma. Stable subgroup separation was achieved by using the 300 transcripts that varied the most. Distinct expression patterns of FSTL5 in each molecular subgroup were confirmed by quantitative real-time polymerase chain reaction. Immunopositivity of FSTL5 identified a large cohort of patients (84 of 235 patients; 36%) at high risk for relapse and death. Importantly, more than 50% of non-WNT/non-SHH tumors displayed FSTL5 negativity, delineating a large patient cohort with a good prognosis who would otherwise be considered intermediate or high-risk on the basis of current molecular subgrouping. CONCLUSION FSTL5 expression denoted a dismal prognosis both within and across medulloblastoma subgroups. The addition of FSTL5 immunohistochemistry to existing molecular stratification schemes constitutes a reliable and cost-effective tool for prognostication in future clinical trials of medulloblastoma.

Genomic imbalances in rhabdomyosarcoma cell lines affect expression of genes frequently altered in primary tumors: an approach to identify candidate genes involved in tumor development.

Rhabdomyosarcomas (RMS) are the most common pediatric soft tissue sarcomas. They resemble developing skeletal muscle and are histologically divided into two main subtypes; alveolar and embryonal RMS. Characteristic genomic aberrations, including the PAX3‐ and PAX7‐FOXO1 fusion genes in alveolar cases, have led to increased understanding of their molecular biology. Here, we determined the effect of genomic copy number on gene expression levels through array comparative genomic hybridization (CGH) analysis of 13 RMS cell lines, confirmed by multiplex ligation‐dependent probe amplification copy number analyses, combined with their corresponding expression profiles. Genes altered at the transcriptional level by genomic imbalances were identified and the effect on expression was proportional to the level of genomic imbalance. Extrapolating to a public expression profiling dataset for 132 primary RMS identified features common to the cell lines and primary samples and associations with subtypes and fusion gene status. Genes identified such as CDK4 and MYCN are known to be amplified, overexpressed, and involved in RMS tumorigenesis. Of the many genes identified, those with likely functional relevance included CENPF, DTL, MYC, EYA2, and FGFR1. Copy number and expression of FGFR1 was validated in additional primary material and found amplified in 6 out of 196 cases and overexpressed relative to skeletal muscle and myoblasts, with significantly higher expression levels in the embryonal compared with alveolar subtypes. This illustrates the ability to identify genes of potential significance in tumor development through combining genomic and transcriptomic profiles from representative cell lines with publicly available expression profiling data from primary tumors.

Evaluation of the similarity of gene expression data estimated with SAGE and Affymetrix GeneChips

BackgroundSerial Analysis of Gene Expression (SAGE) and microarrays have found awidespread application, but much ambiguity exists regarding the evaluation of these technologies. Cross-platform utilization of gene expression data from the SAGE and microarray technology could reduce the need for duplicate experiments and facilitate a more extensive exchange of data within the research community. This requires a measure for the correspondence of the different gene expression platforms. To date, a number of cross-platform evaluations (including a few studies using SAGE and Affymetrix GeneChips) have been conducted showing a variable, but overall low, concordance. This study evaluates these overall measures and introduces the between-ratio difference as a concordance measure pergene.ResultsIn this study, gene expression measurements of Unigene clusters represented by both Affymetrix GeneChips HG-U133A and SAGE were compared using two independent RNA samples. After matching of the data sets the final comparison contains a small data set of 1094 unique Unigene clusters, which is unbiased with respect to expression level. Different overall correlation approaches, like Up/Down classification, contingency tables and correlation coefficients were used to compare both platforms. In addition, we introduce a novel approach to compare two platforms based on the calculation of differences between expression ratios observed in each platform for each individual transcript. This approach results in a concordance measure per gene (with statistical probability value), as opposed to the commonly used overall concordance measures between platforms.ConclusionWe can conclude that intra-platform correlations are generally good, but that overall agreement between the two platforms is modest. This might be due to the binomially distributed sampling variation in SAGE tag counts, SAGE annotation errors and the intensity variation between probe sets of a single gene in Affymetrix GeneChips. We cannot identify or advice which platform performs better since both have their (dis)-advantages. Therefore it is strongly recommended to perform follow-up studies of interesting genes using additional techniques. The newly introduced between-ratio difference is a filtering-independent measure for between-platform concordance. Moreover, the between-ratio difference per gene can be used to detect transcripts with similar regulation on both platforms.

Molecular diagnostics of CNS embryonal tumors

Tremendous progress has recently been made in both molecular subgrouping, and the establishment of prognostic biomarkers for embryonal brain tumors, particularly medulloblastoma. Several prognostic biomarkers that were initially identified in retrospective cohorts of medulloblastoma, including MYC and MYCN amplification, nuclear β-catenin accumulation, and chromosome 17 aberrations have now been validated in clinical trials. Moreover, molecular subgroups based on distinct transcriptome profiles have been consistently reported from various groups on different platforms demonstrating that the concept of distinct medulloblastoma subgroups is very robust. Well-described subgroups of medulloblastomas include tumors showing wingless signaling pathway (Wnt) activation, and another characterized by sonic hedgehog pathway activity. Two or more additional subgroups were consistently reported to contain the vast majority of high-risk tumors, including most tumors with metastatic disease at diagnosis and/or large cell/anaplastic histology. Several years ago, atypical teratoid rhabdoid tumor (AT/RT) was recognized as a separate entity based on its distinct biology and particularly aggressive clinical behavior. These tumors may occur supra or infratentorially and are usually found to have genetic alterations of SMARCB1 (INI1/hSNF5), a tumor suppressor gene located on chromosome 22q. Subsequent loss of SMARCB1 protein expression comprises a relatively specific and sensitive diagnostic marker for AT/RT. For CNS primitive neuroectodermal tumors (CNS PNETs), a consistent finding has been that they are molecularly distinct from medulloblastoma. Furthermore, a distinct fraction of CNS PNETs with particularly poor prognosis only occurring in young children was delineated, which was previously labeled ependymoblastoma or embryonal tumor with abundant neuropil and true rosettes (ETANTR) and which is morphologically characterized by the presence of multilayered “ependymoblastic” rosettes. This group of tumors shows a unique cytogenetic abnormality not seen in other brain tumors: focal amplification of a micro-RNA cluster at chromosome 19q13.42, which has never been found to be amplified in other CNS PNETs, medulloblastoma or AT/RT. In summary, these consistent findings have significantly contributed to our ability to sub-classify embryonal brain tumors into clinically and biologically meaningful strata and, for some of the subgroups, have led to the identification of specific targets for future development of molecularly targeted therapies.

Silencing of SPRY1 Triggers Complete Regression of Rhabdomyosarcoma Tumors Carrying a Mutated RAS Gene

RAS oncogenes are among the most frequently mutated genes in human cancer, but effective strategies for therapeutic inhibition of the RAS pathway have been elusive. Sprouty1 (SPRY1) is an upstream antagonist of RAS that is activated by extracellular signal-related kinase (ERK), providing a negative feedback loop for RAS signaling, and other evidence suggests that SPRY1 may have a tumor suppressor function. Studies of RAS status in the human childhood tumor rhabdomyosarcoma (RMS) indicated mutations in approximately half of the tumors of the embryonal rhabdomyosarcoma subtype (ERMS) but not the alveolar subtype (ARMS). ERMS tumors also showed overexpression of SPRY1, which was indeed upregulated by mutant RAS. However, we found that, in the presence of mutant RAS, the function of SPRY1 was changed from an antagonist to an agonist of RAS signaling. Thus, SPRY1 supported formation of activated ERK and mitogen-activated protein/ERK kinase and was essential for ERMS cell proliferation and survival. Conversely, silencing of SPRY1 in ERMS cells (but not ARMS cells) abolished their tumorigenicity in mice. Moreover, silencing of SPRY1 caused regression of established ERMS tumors (but not ARMS tumors) formed in xenograft settings. Our findings argue that SPRY1 inhibition can offer a therapeutic strategy to treat childhood RMS and possibly other tumors carrying oncogenic RAS mutations.

OTX2 sustains a bivalent-like state of OTX2-bound promoters in medulloblastoma by maintaining their H3K27me3 levels

Recent studies showed frequent mutations in histone H3 lysine 27 (H3K27) demethylases in medulloblastomas of Group 3 and Group 4, suggesting a role for H3K27 methylation in these tumors. Indeed, trimethylated H3K27 (H3K27me3) levels were shown to be higher in Group 3 and 4 tumors compared to WNT and SHH medulloblastomas, also in tumors without detectable mutations in demethylases. Here, we report that polycomb genes, required for H3K27 methylation, are consistently upregulated in Group 3 and 4 tumors. These tumors show high expression of the homeobox transcription factor OTX2. Silencing of OTX2 in D425 medulloblastoma cells resulted in downregulation of polycomb genes such as EZH2, EED, SUZ12 and RBBP4 and upregulation of H3K27 demethylases KDM6A, KDM6B,JARID2 and KDM7A. This was accompanied by decreased H3K27me3 and increased H3K27me1 levels in promoter regions. Strikingly, the decrease of H3K27me3 was most prominent in promoters that bind OTX2. OTX2-bound promoters showed high levels of the H3K4me3 and H3K9ac activation marks and intermediate levels of the H3K27me3 inactivation mark, reminiscent of a bivalent modification. After silencing of OTX2, H3K27me3 levels strongly dropped, but H3K4me3 and H3K9ac levels remained high. OTX2-bound bivalent genes showed high expression levels in D425, but the expression of most of these genes did not change after OTX2 silencing and loss of the H3K27me3 mark. Maintaining promoters in a bivalent state by sustaining H3K27 trimethylation therefore seems to be an important function of OTX2 in medulloblastoma, while other transcription factors might regulate the actual expression levels of these genes.

High levels of PROM1 (CD133) transcript are a potential predictor of poor prognosis in medulloblastoma

The surface marker PROM1 is considered one of the most important markers of tumor-initiating cells, and its expression is believed to be an adverse prognostic factor in gliomas and in other malignancies. To date, to our knowledge, no specific studies of its expression in medulloblastoma series have been performed. The aims of our study were to evaluate the expression profile of the PROM1 gene in medulloblastoma and to assess its possible role as a prognostic factor. The PROM1 gene expression was evaluated by quantitative- polymerase chain reaction on 45 medulloblastoma samples by using specific dye-labeled probe systems. A significantly higher expression of PROM1 was found both in patients with poorer prognosis (P= .007) and in those with metastasis (P= .03). Kaplan-Meier analysis showed that both overall survival (OS) and progression-free survival (PFS) were shorter in patients with higher PROM1 mRNA levels than in patients with lower expression, even when the desmoplastic cases were excluded (P= .0004 and P= .002, for OS and PFS for all cases, respectively; P= .002 and P= .008 for OS and PFS for nondesmoplastic cases, respectively). Cox regression model demonstrated that PROM1 expression is an independent prognostic factor (hazard ratio, 4.56; P= .008). The result was validated on an independent cohort of 42 cases by microarray-based analysis (P= .019). This work suggests that high mRNA levels of PROM1 are associated with poor outcome in pediatric medulloblastoma. Furthermore, high PROM1 expression levels seem to increase the likelihood of metastases. Such results need to be confirmed in larger prospective series to possibly incorporate PROM1 gene expression into risk classification systems to be used in the clinical setting.

Mismatch repair deficiency: a temozolomide resistance factor in medulloblastoma cell lines that is uncommon in primary medulloblastoma tumours

Background:Tumours are responsive to temozolomide (TMZ) if they are deficient in O6-methylguanine-DNA methyltransferase (MGMT), and mismatch repair (MMR) proficient.Methods:The effect of TMZ on medulloblastoma (MB) cell killing was analysed with clonogenic survival assays. Expression of DNA repair genes and enzymes was investigated using microarrays, western blot, and immunohistochemistry. DNA sequencing and promoter methylation analysis were employed to investigate the cause of loss of the expression of MMR gene MLH1.Results:Temozolomide exhibited potent cytotoxic activity in D425Med (MGMT deficient, MLH1 proficient; IC50=1.7 μM), moderate activity against D341Med (MGMT proficient, MLH1 deficient), and DAOY MB cells (MGMT proficient, MLH1 proficient). MGMT inhibitor O6-benzylguanine sensitised DAOY, but not D341Med cells to TMZ. Of 12 MB cell lines, D341Med, D283Med, and 1580WÜ cells exhibited MMR deficiency due to MLH1 promoter hypermethylation. DNA sequencing of these cells provided no evidence for somatic genetic alterations in MLH1. Expression analyses of MMR and MGMT in MB revealed that all patient specimens (n=74; expression array, n=61; immunostaining, n=13) are most likely MMR proficient, whereas some tumours had low MGMT expression levels (according to expression array) or were totally MGMT deficient (3 out of 13 according to immunohistochemistry).Conclusion:A subset of MB may respond to TMZ as some patient specimens are MGMT deficient, and tumours appear to be MMR proficient.

WNT activation by lithium abrogates TP53 mutation associated radiation resistance in medulloblastoma

TP53 mutations confer subgroup specific poor survival for children with medulloblastoma. We hypothesized that WNT activation which is associated with improved survival for such children abrogates TP53 related radioresistance and can be used to sensitize TP53 mutant tumors for radiation. We examined the subgroup-specific role of TP53 mutations in a cohort of 314 patients treated with radiation. TP53 wild-type or mutant human medulloblastoma cell-lines and normal neural stem cells were used to test radioresistance of TP53 mutations and the radiosensitizing effect of WNT activation on tumors and the developing brain. Children with WNT/TP53 mutant medulloblastoma had higher 5-year survival than those with SHH/TP53 mutant tumours (100% and 36.6% ± 8.7%, respectively (p < 0.001)). Introduction of TP53 mutation into medulloblastoma cells induced radioresistance (survival fractions at 2Gy (SF2) of 89% ± 2% vs. 57.4% ± 1.8% (p < 0.01)). In contrast, β-catenin mutation sensitized TP53 mutant cells to radiation (p < 0.05). Lithium, an activator of the WNT pathway, sensitized TP53 mutant medulloblastoma to radiation (SF2 of 43.5% ± 1.5% in lithium treated cells vs. 56.6 ± 3% (p < 0.01)) accompanied by increased number of γH2AX foci. Normal neural stem cells were protected from lithium induced radiation damage (SF2 of 33% ± 8% for lithium treated cells vs. 27% ± 3% for untreated controls (p = 0.05). Poor survival of patients with TP53 mutant medulloblastoma may be related to radiation resistance. Since constitutive activation of the WNT pathway by lithium sensitizes TP53 mutant medulloblastoma cells and protect normal neural stem cells from radiation, this oral drug may represent an attractive novel therapy for high-risk medulloblastomas.