Nanne K. Kloosterhof

Intrinsic Gene Expression Profiles of Gliomas Are a Better Predictor of Survival than Histology

Gliomas are the most common primary brain tumors with heterogeneous morphology and variable prognosis. Treatment decisions in patients rely mainly on histologic classification and clinical parameters. However, differences between histologic subclasses and grades are subtle, and classifying gliomas is subject to a large interobserver variability. To improve current classification standards, we have performed gene expression profiling on a large cohort of glioma samples of all histologic subtypes and grades. We identified seven distinct molecular subgroups that correlate with survival. These include two favorable prognostic subgroups (median survival, >4.7 years), two with intermediate prognosis (median survival, 1-4 years), two with poor prognosis (median survival, <1 year), and one control group. The intrinsic molecular subtypes of glioma are different from histologic subgroups and correlate better to patient survival. The prognostic value of molecular subgroups was validated on five independent sample cohorts (The Cancer Genome Atlas, Repository for Molecular Brain Neoplasia Data, GSE12907, GSE4271, and Li and colleagues). The power of intrinsic subtyping is shown by its ability to identify a subset of prognostically favorable tumors within an external data set that contains only histologically confirmed glioblastomas (GBM). Specific genetic changes (epidermal growth factor receptor amplification, IDH1 mutation, and 1p/19q loss of heterozygosity) segregate in distinct molecular subgroups. We identified a subgroup with molecular features associated with secondary GBM, suggesting that different genetic changes drive gene expression profiles. Finally, we assessed response to treatment in molecular subgroups. Our data provide compelling evidence that expression profiling is a more accurate and objective method to classify gliomas than histologic classification. Molecular classification therefore may aid diagnosis and can guide clinical decision making.

Isocitrate dehydrogenase-1 mutations: a fundamentally new understanding of diffuse glioma?

The discovery of somatic mutations in the gene encoding isocitrate dehydrogenase-1 (IDH1) in glioblastomas was remarkable because the enzyme was not previously identified with any known oncogenic pathway. IDH1 is mutated in up to 75% of grade II and grade III diffuse gliomas. Apart from acute myeloid leukaemia, other tumour types do not carry IDH1 mutations. Mutations in a homologous gene, IDH2, have also been identified, although they are much rarer. Although TP53 mutations and 1p/19q codeletions are mutually exclusive in gliomas, in both of these genotypes IDH1 mutations are common. IDH1 and IDH2 mutations are early events in the development of gliomas. Moreover, IDH1 and IDH2 mutations are a major prognostic marker for overall and progression-free survival in grade II-IV gliomas. Mutated IDH1 has an altered catalytic activity that results in the accumulation of 2-hydroxyglutarate. Molecularly, IDH1 and IDH2 mutations are heterozygous, affect only a single codon, and rarely occur together. Because IDH1 does not belong to a traditional oncogenic pathway and is specifically and commonly mutated in gliomas, the altered enzymatic activity of IDH1 may provide a fundamentally new understanding of diffuse glioma.

Cytogenetic Prognostication Within Medulloblastoma Subgroups

PURPOSE Medulloblastoma comprises four distinct molecular subgroups: WNT, SHH, Group 3, and Group 4. Current medulloblastoma protocols stratify patients based on clinical features: patient age, metastatic stage, extent of resection, and histologic variant. Stark prognostic and genetic differences among the four subgroups suggest that subgroup-specific molecular biomarkers could improve patient prognostication. PATIENTS AND METHODS Molecular biomarkers were identified from a discovery set of 673 medulloblastomas from 43 cities around the world. Combined risk stratification models were designed based on clinical and cytogenetic biomarkers identified by multivariable Cox proportional hazards analyses. Identified biomarkers were tested using fluorescent in situ hybridization (FISH) on a nonoverlapping medulloblastoma tissue microarray (n = 453), with subsequent validation of the risk stratification models. RESULTS Subgroup information improves the predictive accuracy of a multivariable survival model compared with clinical biomarkers alone. Most previously published cytogenetic biomarkers are only prognostic within a single medulloblastoma subgroup. Profiling six FISH biomarkers (GLI2, MYC, chromosome 11 [chr11], chr14, 17p, and 17q) on formalin-fixed paraffin-embedded tissues, we can reliably and reproducibly identify very low-risk and very high-risk patients within SHH, Group 3, and Group 4 medulloblastomas. CONCLUSION Combining subgroup and cytogenetic biomarkers with established clinical biomarkers substantially improves patient prognostication, even in the context of heterogeneous clinical therapies. The prognostic significance of most molecular biomarkers is restricted to a specific subgroup. We have identified a small panel of cytogenetic biomarkers that reliably identifies very high-risk and very low-risk groups of patients, making it an excellent tool for selecting patients for therapy intensification and therapy de-escalation in future clinical trials.

IDH1 R132H decreases proliferation of glioma cell lines in vitro and in vivo

A high percentage of grade II and III gliomas have mutations in the gene encoding isocitrate dehydrogenase (IDH1). This mutation is always a heterozygous point mutation that affects the amino acid arginine at position 132 and results in loss of its native enzymatic activity and gain of alternative enzymatic activity (producing D‐2‐hydroxyglutarate). The objective of this study was to investigate the cellular effects of R132H mutations in IDH1.

Segregation of Non-p.R132H Mutations in IDH1 in Distinct Molecular Subtypes of Glioma

Mutations in the gene encoding the isocitrate dehydrogenase 1 gene (IDH1) occur at a high frequency (up to 80%) in many different subtypes of glioma. In this study, we have screened for IDH1 mutations in a cohort of 496 gliomas. IDH1 mutations were most frequently observed in low grade gliomas with c.395G>A (p.R132H) representing >90% of all IDH1 mutations. Interestingly, non‐p.R132H mutations segregate in distinct histological and molecular subtypes of glioma. Histologically, they occur sporadically in classic oligodendrogliomas and at significantly higher frequency in other grade II and III gliomas. Genetically, non‐p.R132H mutations occur in tumors with TP53 mutation, are virtually absent in tumors with loss of heterozygosity on 1p and 19q and accumulate in distinct (gene‐expression profiling based) intrinsic molecular subtypes. The IDH1 mutation type does not affect patient survival. Our results were validated on an independent sample cohort, indicating that the IDH1 mutation spectrum may aid glioma subtype classification. Functional differences between p.R132H and non‐p.R132H mutated IDH1 may explain the segregation in distinct glioma subtypes.

Integrated genomic profiling identifies candidate genes implicated in glioma-genesis and a novel LEO1-SLC12A1 fusion gene.

We performed genotyping and exon‐level expression profiling on 21 glioblastomas (GBMs) and 19 oligodendrogliomas (ODs) to identify genes involved in glioma initiation and/or progression. Low‐copy number amplifications (2.5 < n < 7) and high‐copy number amplifications (n > 7) were more frequently observed in GBMs; ODs generally have more heterozygous deletions per tumor. Four high‐copy amplicons were identified in more than one sample and resulted in overexpression of the known oncogenes EGFR, MDM2, and CDK4. In the fourth amplicon, RBBP5, a member of the RB pathway, may act as a novel oncogene in GBMs. Not all hCNAs contain known genes, which may suggest that other transcriptional and/or regulatory elements are the target for amplification. Regions with most frequent allelic loss, both in ODs and GBMs, resulted in a reduced expression of known tumor suppressor genes. We identified a homozygous deletion spanning the Pragmin gene in one sample, but direct sequencing of all coding exons in 20 other glioma samples failed to detect additional genetic changes. Finally, we screened for fusion genes by identifying aberrant 5′‐3′ expression of genes that lie over regions of a copy number change. A fusion gene between exon 11 of LEO1 and exon 10 of SLC12A1 was identified. Our data show that integrated genomic profiling can identify genes involved in tumor initiation, and/or progression and can be used as an approach to identify novel fusion genes.

The CASPR2 cell adhesion molecule functions as a tumor suppressor gene in glioma

Genomic translocations have been implicated in cancer. In this study, we performed a screen for genetic translocations in gliomas based on exon-level expression profiles. We identified a translocation in the contactin-associated protein-like 2 (CASPR2) gene, encoding a cell adhesion molecule. CASPR2 mRNA was fused to an expressed sequence tag that likely is part of the nuclear receptor coactivator 1 gene. Despite high mRNA expression levels, no CASPR2 fusion protein was detected. In a set of 25 glioblastomas and 22 oligodendrogliomas, mutation analysis identified two additional samples with genetic alterations in the CASPR2 gene and all three identified genetic alterations are likely to reduce CASPR2 protein expression levels. Methylation of the CASPR2 gene was also observed in gliomas and glioma cell lines. CASPR2-overexpressing cells showed decreased proliferation rates, likely because of an increase in apoptosis. Moreover, high CASPR2 mRNA expression level is positively correlated with survival and is an independent prognostic factor. These results indicate that CASPR2 acts as a tumor suppressor gene in glioma.

Molecular subtypes of glioma identified by genome‐wide methylation profiling

Recent studies have indicated a prognostic role for genome‐wide methylation in gliomas: Tumors that show an overall increase in DNA methylation at CpG sites (CIMP+; CpG island methylator phenotype) have a more favorable prognosis than CIMP− gliomas. Here, we have determined whether methylation profiling can identify more and clinically relevant molecular subtypes of glioma by performing genome‐wide methylation profiling on 138 glial brain tumors of all histological diagnosis. Hopach (Hierarchical ordered partitioning and collapsing hybrid) clustering using the 1,000 most variable CpGs identified three distinct glioma subtypes (C+1p19q, C+wt, and C−) and one adult brain subtype. All “C+1p19q” and “C+wt” tumors were CIMP+ whereas most (50/54) “C−” tumors were CIMP−. The C− subtype gliomas contained many glioblastomas and all pilocytic astrocytomas. 1p19q LOH was frequent in the C+1p19q subtype. Other genetic changes (IDH1 mutation and EGFR amplification) and gene‐expression based molecular subtypes also segregated in distinct methylation subtypes, demonstrating that these subtypes are also genetically distinct. Each subtype was associated with its own prognosis: median survival for C−, C+1p19q, and C+wt tumors was 1.18, 5.00, and 2.62 years, respectively. The prognostic value of these methylation subtypes was validated on an external dataset from the TCGA. Analysis of recurrences of 14 primary tumors samples indicates that shifts between some C+wt and C+1p/19q tumors can occur between the primary and recurrent tumor, but CIMP status remained stable. Our data demonstrate that methylation profiling identifies at least three prognostically relevant subtypes of glioma that can aid diagnosis and potentially guide treatment for patients.

Mutation specific functions of EGFR result in a mutation-specific downstream pathway activation

BACKGROUND Epidermal growth factor receptor (EGFR) is frequently mutated in various types of cancer. Although all oncogenic mutations are considered activating, different tumour types have different mutation spectra. It is possible that functional differences underlie this tumour-type specific mutation spectrum. METHODS We have determined whether specific mutations in EGFR (EGFR, EGFRvIII and EGFR-L858R) have differences in binding partners, differences in downstream pathway activation (gene expression and phosphoproteins), and have functional consequences on cellular growth and migration. RESULTS Using biotin pulldown and subsequent mass spectrometry we were able to detect mutation specific binding partners for EGFR. Differential binding was confirmed using a proximity ligation assay and/or Western Blot for the dedicator of cytokinesis 4 (DOCK4), UDP-glucose glycoprotein glucosyltransferase 1 (UGGT1), MYC binding protein 2 (MYCBP2) and Smoothelin (SMTN). We also demonstrate that each mutation induces the expression of a specific set of genes, and that each mutation is associated with specific phosphorylation patterns. Finally, we demonstrate using stably expressing cell lines that EGFRvIII and EGFL858R display reduced growth and migration compared to EGFR wildtype expressing cells. CONCLUSION Our results indicate that there are distinct functional differences between different EGFR mutations. The functional differences between different mutations argue for the development of mutation specific targeted therapies.

Trp53 inactivation leads to earlier phaeochromocytoma formation in pten knockout mice

Phaeochromocytomas (PCCs) are benign neuroendocrine tumours of the adrenal medulla. Approximately 10% of PCC patients develop metastases, but this frequency is much higher in specific subtypes of patients. The reliable diagnosis of malignant PCC can only be made after identification of a metastasis. To study the effect of Trp53 inactivation on PCC pathogenesis in Pten KO mice, we investigated the adrenals of a large cohort of mice with conditional monoallelic and biallelic inactivation of Trp53 and Pten. The adrenal weights were determined for all mice, and in a proportion of these mice, immunohistochemistry for tyrosine hydroxylase and dopamine β-hydroxylase was performed on the adrenals and corresponding lungs. Finally, comparative genomic hybridization (CGH) was performed. The histological and immunohistochemical results confirmed that the adrenal tumours were PCCs. Inactivation of one or both alleles of Trp53 resulted in earlier tumour occurrence in the Pten(loxP/loxP) mice as well as in the Pten(loxP/+) mice. In addition, lung metastases were found in up to 67% of mice. The CGH results showed that the most frequent genomic alterations were loss of chromosome 19 (86%) and gain of chromosome 15 (71%). In this study, we have shown that Pten/Trp53 KO mice showed metastatic PCC at high frequency and primary tumours occurred at younger ages in mice with Trp53 inactivation. Therefore, the present model appears to be a suitable model that might allow the preclinical study of new therapeutics for these tumours.