Yvonne Kahlert

Customized Oligonucleotide Microarray Gene Expression–Based Classification of Neuroblastoma Patients Outperforms Current Clinical Risk Stratification

PURPOSE To develop a gene expression-based classifier for neuroblastoma patients that reliably predicts courses of the disease. PATIENTS AND METHODS Two hundred fifty-one neuroblastoma specimens were analyzed using a customized oligonucleotide microarray comprising 10,163 probes for transcripts with differential expression in clinical subgroups of the disease. Subsequently, the prediction analysis for microarrays (PAM) was applied to a first set of patients with maximally divergent clinical courses (n = 77). The classification accuracy was estimated by a complete 10-times-repeated 10-fold cross validation, and a 144-gene predictor was constructed from this set. This classifiers predictive power was evaluated in an independent second set (n = 174) by comparing results of the gene expression-based classification with those of risk stratification systems of current trials from Germany, Japan, and the United States. RESULTS The first set of patients was accurately predicted by PAM (cross-validated accuracy, 99%). Within the second set, the PAM classifier significantly separated cohorts with distinct courses (3-year event-free survival [EFS] 0.86 +/- 0.03 [favorable; n = 115] v 0.52 +/- 0.07 [unfavorable; n = 59] and 3-year overall survival 0.99 +/- 0.01 v 0.84 +/- 0.05; both P < .0001) and separated risk groups of current neuroblastoma trials into subgroups with divergent outcome (NB2004: low-risk 3-year EFS 0.86 +/- 0.04 v 0.25 +/- 0.15, P < .0001; intermediate-risk 1.00 v 0.57 +/- 0.19, P = .018; high-risk 0.81 +/- 0.10 v 0.56 +/- 0.08, P = .06). In a multivariate Cox regression model, the PAM predictor classified patients of the second set more accurately than risk stratification of current trials from Germany, Japan, and the United States (P < .001; hazard ratio, 4.756 [95% CI, 2.544 to 8.893]). CONCLUSION Integration of gene expression-based class prediction of neuroblastoma patients may improve risk estimation of current neuroblastoma trials.

Telomerase activation by genomic rearrangements in high-risk neuroblastoma

Neuroblastoma is a malignant paediatric tumour of the sympathetic nervous system. Roughly half of these tumours regress spontaneously or are cured by limited therapy. By contrast, high-risk neuroblastomas have an unfavourable clinical course despite intensive multimodal treatment, and their molecular basis has remained largely elusive. Here we have performed whole-genome sequencing of 56 neuroblastomas (high-risk, n = 39; low-risk, n = 17) and discovered recurrent genomic rearrangements affecting a chromosomal region at 5p15.33 proximal of the telomerase reverse transcriptase gene (TERT). These rearrangements occurred only in high-risk neuroblastomas (12/39, 31%) in a mutually exclusive fashion with MYCN amplifications and ATRX mutations, which are known genetic events in this tumour type. In an extended case series (n = 217), TERT rearrangements defined a subgroup of high-risk tumours with particularly poor outcome. Despite a large structural diversity of these rearrangements, they all induced massive transcriptional upregulation of TERT. In the remaining high-risk tumours, TERT expression was also elevated in MYCN-amplified tumours, whereas alternative lengthening of telomeres was present in neuroblastomas without TERT or MYCN alterations, suggesting that telomere lengthening represents a central mechanism defining this subtype. The 5p15.33 rearrangements juxtapose the TERT coding sequence to strong enhancer elements, resulting in massive chromatin remodelling and DNA methylation of the affected region. Supporting a functional role of TERT, neuroblastoma cell lines bearing rearrangements or amplified MYCN exhibited both upregulated TERT expression and enzymatic telomerase activity. In summary, our findings show that remodelling of the genomic context abrogates transcriptional silencing of TERT in high-risk neuroblastoma and places telomerase activation in the centre of transformation in a large fraction of these tumours.

Prognostic Impact of Gene Expression–Based Classification for Neuroblastoma

PURPOSE To evaluate the impact of a predefined gene expression-based classifier for clinical risk estimation and cytotoxic treatment decision making in neuroblastoma patients. PATIENTS AND METHODS Gene expression profiles of 440 internationally collected neuroblastoma specimens were investigated by microarray analysis, 125 of which were examined prospectively. Patients were classified as either favorable or unfavorable by a 144-gene prediction analysis for microarrays (PAM) classifier established previously on a separate set of 77 patients. PAM classification results were compared with those of current prognostic markers and risk estimation strategies. RESULTS The PAM classifier reliably distinguished patients with contrasting clinical courses (favorable [n = 249] and unfavorable [n = 191]; 5-year event free survival [EFS] 0.84 +/- 0.03 v 0.38 +/- 0.04; 5-year overall survival [OS] 0.98 +/- 0.01 v 0.56 +/- 0.05, respectively; both P < .001). Moreover, patients with divergent outcome were robustly discriminated in both German and international cohorts and in prospectively analyzed samples (P <or= .001 for both EFS and OS for each). In subgroups with clinical low-, intermediate-, and high-risk of death from disease, the PAM predictor significantly separated patients with divergent outcome (low-risk 5-year OS: 1.0 v 0.75 +/- 0.10, P < .001; intermediate-risk: 1.0 v 0.82 +/- 0.08, P = .042; and high-risk: 0.81 +/- 0.08 v 0.43 +/- 0.05, P = .001). In multivariate Cox regression models based on both EFS and OS, PAM was a significant independent prognostic marker (EFS: hazard ratio [HR], 3.375; 95% CI, 2.075 to 5.492; P < .001; OS: HR, 11.119, 95% CI, 2.487 to 49.701; P < .001). The highest potential clinical impact of the classifier was observed in patients currently considered as non-high-risk (n = 289; 5-year EFS: 0.87 +/- 0.02 v 0.44 +/- 0.07; 5-year OS: 1.0 v 0.80 +/- 0.06; both P < .001). CONCLUSION Gene expression-based classification using the 144-gene PAM predictor can contribute to improved treatment stratification of neuroblastoma patients.

Mutational dynamics between primary and relapse neuroblastomas

Neuroblastoma is a malignancy of the developing sympathetic nervous system that is often lethal when relapse occurs. We here used whole-exome sequencing, mRNA expression profiling, array CGH and DNA methylation analysis to characterize 16 paired samples at diagnosis and relapse from individuals with neuroblastoma. The mutational burden significantly increased in relapsing tumors, accompanied by altered mutational signatures and reduced subclonal heterogeneity. Global allele frequencies at relapse indicated clonal mutation selection during disease progression. Promoter methylation patterns were consistent over disease course and were patient specific. Recurrent alterations at relapse included mutations in the putative CHD5 neuroblastoma tumor suppressor, chromosome 9p losses, DOCK8 mutations, inactivating mutations in PTPN14 and a relapse-specific activity pattern for the PTPN14 target YAP. Recurrent new mutations in HRAS, KRAS and genes mediating cell-cell interaction in 13 of 16 relapse tumors indicate disturbances in signaling pathways mediating mesenchymal transition. Our data shed light on genetic alteration frequency, identity and evolution in neuroblastoma.

Hox-C9 activates the intrinsic pathway of apoptosis and is associated with spontaneous regression in neuroblastoma

Neuroblastoma is an embryonal malignancy of the sympathetic nervous system. Spontaneous regression and differentiation of neuroblastoma is observed in a subset of patients, and has been suggested to represent delayed activation of physiologic molecular programs of fetal neuroblasts. Homeobox genes constitute an important family of transcription factors, which play a fundamental role in morphogenesis and cell differentiation during embryogenesis. In this study, we demonstrate that expression of the majority of the human HOX class I homeobox genes is significantly associated with clinical covariates in neuroblastoma using microarray expression data of 649 primary tumors. Moreover, a HOX gene expression-based classifier predicted neuroblastoma patient outcome independently of age, stage and MYCN amplification status. Among all HOX genes, HOXC9 expression was most prominently associated with favorable prognostic markers. Most notably, elevated HOXC9 expression was significantly associated with spontaneous regression in infant neuroblastoma. Re-expression of HOXC9 in three neuroblastoma cell lines led to a significant reduction in cell viability, and abrogated tumor growth almost completely in neuroblastoma xenografts. Neuroblastoma growth arrest was related to the induction of programmed cell death, as indicated by an increase in the sub-G1 fraction and translocation of phosphatidylserine to the outer membrane. Programmed cell death was associated with the release of cytochrome c from the mitochondria into the cytosol and activation of the intrinsic cascade of caspases, indicating that HOXC9 re-expression triggers the intrinsic apoptotic pathway. Collectively, our results show a strong prognostic impact of HOX gene expression in neuroblastoma, and may point towards a role of Hox-C9 in neuroblastoma spontaneous regression.

High ALK Receptor Tyrosine Kinase Expression Supersedes ALK Mutation as a Determining Factor of an Unfavorable Phenotype in Primary Neuroblastoma

Purpose: Genomic alterations of the anaplastic lymphoma kinase (ALK) gene have been postulated to contribute to neuroblastoma pathogenesis. This study aimed to determine the interrelation of ALK mutations, ALK expression levels, and clinical phenotype in primary neuroblastoma. Experimental Design: The genomic ALK status and global gene expression patterns were examined in 263 primary neuroblastomas. Allele-specific ALK expression was determined by cDNA cloning and sequencing. Associations of genomic ALK alterations and ALK expression levels with clinical phenotypes and transcriptomic profiles were compared. Results: Nonsynonymous point mutations of ALK were detected in 21 of 263 neuroblastomas (8%). Tumors with ALK mutations exhibited about 2-fold elevated median ALK mRNA levels in comparison with tumors with wild-type (WT) ALK. Unexpectedly, the WT allele was preferentially expressed in 12 of 21 mutated tumors. Whereas survival of patients with ALK mutated tumors was significantly worse as compared with the entire cohort of WT ALK patients, it was similarly poor in patients with WT ALK tumors in which ALK expression was as high as in ALK mutated neuroblastomas. Global gene expression patterns of tumors with ALK mutations or with high-level WT ALK expression were highly similar, and suggested that ALK may be involved in cellular proliferation in primary neuroblastoma. Conclusions: Primary neuroblastomas with mutated ALK exhibit high ALK expression levels and strongly resemble neuroblastomas with elevated WT ALK expression levels in both their clinical and molecular phenotypes. These data suggest that high levels of mutated and WT ALK mediate similar molecular functions that may contribute to a malignant phenotype in primary neuroblastoma. Clin Cancer Res; 17(15); 5082–92. ©2011 AACR.

Differential Expression of Neuronal Genes Defines Subtypes of Disseminated Neuroblastoma with Favorable and Unfavorable Outcome

Purpose: Identification of molecular characteristics of spontaneously regressing stage IVS and progressing stage IV neuroblastoma to improve discrimination of patients with metastatic disease following favorable and unfavorable clinical courses. Experimental Design: Serial analysis of gene expression profiles were generated from five stage IVS and three stage IV neuroblastoma. Differential expression of candidate genes was evaluated by real-time quantitative reverse transcription-PCR in 76 pretreatment tumor samples (stage IVS n = 27 and stage IV n = 49). Gene expression-based outcome prediction was determined by Prediction Analysis for Microarrays using 38 tumors as a training set and 38 tumors as a test set. Results: Comparison of serial analysis of gene expression profiles from stage IV and IVS neuroblastoma revealed ∼500 differentially expressed transcripts. Genes related to neuronal differentiation were observed more frequently in stage IVS tumors as determined by associating transcripts to Gene Ontology annotations. Forty-one candidate genes were evaluated by quantitative reverse transcription-PCR and 18 were confirmed to be differentially expressed (P ≤ 0.001). Classification of patients according to expression patterns of these 18 genes using Prediction Analysis for Microarrays discriminated two subgroups with significantly differing event-free survival (96 ± 6% versus 40 ± 8% at 3 years; P < 0.0001) and overall survival (100% versus 72 ± 7% at 3 years; P = 0.0003). This classifier was the only independent covariate marker in a multivariate analysis considering the variables stage, age, MYCN amplification, and gene signature. Conclusions: Spontaneously regressing and progressing metastatic neuroblastoma differ by specific gene expression patterns, indicating distinct levels of neuronal differentiation and allowing for an improved risk estimation of children with disseminated disease.

Revised risk estimation and treatment stratification of low- and intermediate-risk neuroblastoma patients by integrating clinical and molecular prognostic markers

Purpose: To optimize neuroblastoma treatment stratification, we aimed at developing a novel risk estimation system by integrating gene expression–based classification and established prognostic markers. Experimental Design: Gene expression profiles were generated from 709 neuroblastoma specimens using customized 4 × 44 K microarrays. Classification models were built using 75 tumors with contrasting courses of disease. Validation was performed in an independent test set (n = 634) by Kaplan–Meier estimates and Cox regression analyses. Results: The best-performing classifier predicted patient outcome with an accuracy of 0.95 (sensitivity, 0.93; specificity, 0.97) in the validation cohort. The highest potential clinical value of this predictor was observed for current low-risk patients [5-year event-free survival (EFS), 0.84 ± 0.02 vs. 0.29 ± 0.10; 5-year overall survival (OS), 0.99 ± 0.01 vs. 0.76 ± 0.11; both P < 0.001] and intermediate-risk patients (5-year EFS, 0.88 ± 0.06 vs. 0.41 ± 0.10; 5-year OS, 1.0 vs. 0.70 ± 0.09; both P < 0.001). In multivariate Cox regression models for low-risk/intermediate-risk patients, the classifier outperformed risk assessment of the current German trial NB2004 [EFS: hazard ratio (HR), 5.07; 95% confidence interval (CI), 3.20–8.02; OS: HR, 25.54; 95% CI, 8.40–77.66; both P < 0.001]. On the basis of these findings, we propose to integrate the classifier into a revised risk stratification system for low-risk/intermediate-risk patients. According to this system, we identified novel subgroups with poor outcome (5-year EFS, 0.19 ± 0.08; 5-year OS, 0.59 ± 0.1), for whom we propose intensified treatment, and with beneficial outcome (5-year EFS, 0.87 ± 0.05; 5-year OS, 1.0), who may benefit from treatment de-escalation. Conclusions: Combination of gene expression–based classification and established prognostic markers improves risk estimation of patients with low-risk/intermediate-risk neuroblastoma. We propose to implement our revised treatment stratification system in a prospective clinical trial. Clin Cancer Res; 21(8); 1904–15. ©2014 AACR. See related commentary by Attiyeh and Maris, p. 1782

Subclassification and Individual Survival Time Prediction from Gene Expression Data of Neuroblastoma Patients by Using CASPAR

Purpose: To predict individual survival times for neuroblastoma patients from gene expression data using the cancer survival prediction using automatic relevance determination (CASPAR) algorithm. Experimental Design: A first set of oligonucleotide microarray gene expression profiles comprising 256 neuroblastoma patients was generated. Then, CASPAR was combined with a leave-one-out cross-validation to predict individual times for both the whole cohort and subgroups of patients with unfavorable markers, including stage 4 disease (n = 67), unfavorable genetic alterations, intermediate-risk or high-risk stratification by the German neuroblastoma trial, and patients predicted as unfavorable by a recently described gene expression classifier (n = 83). Prediction accuracy of individual survival times was assessed by Kaplan-Meier analyses and time-dependent receiver operator characteristics curve analyses. Subsequently, classification results were validated in an independent cohort (n = 120). Results: CASPAR separated patients with divergent outcome in both the initial and the validation cohort [initial set, 5y-OS 0.94 ± 0.04 (predicted long survival) versus 0.38 ± 0.17 (predicted short survival), P < 0.0001; validation cohort, 5y-OS 0.94 ± 0.07 (long) versus 0.40 ± 0.13 (short), P < 0.0001]. Time-dependent receiver operator characteristics analyses showed that CASPAR-predicted individual survival times were highly accurate (initial set, mean area under the curve for first 10 years of overall survival prediction 0.92 ± 0.04; validation set, 0.81 ± 0.05). Furthermore, CASPAR significantly discriminated short (<5 years) from long survivors (>5 years) in subgroups of patients with unfavorable markers with the exception of MYCN-amplified patients (initial set). Confirmatory results with high significance were observed in the validation cohort [stage 4 disease (P = 0.0049), NB2004 intermediate-risk or high-risk stratification (P = 0.0017), and unfavorable gene expression prediction (P = 0.0017)]. Conclusions: CASPAR accurately forecasts individual survival times for neuroblastoma patients from gene expression data.

Identification of DEIN, a Novel Gene with High Expression Levels in Stage IVS Neuroblastoma

Neuroblastoma at stage IVS, defined by dissemination to specific tissues and age <1 year at diagnosis, regularly follows spontaneous regression without cytotoxic treatment. To uncover the molecular characteristics of this subtype, Serial Analysis of Gene Expression (SAGE) profiles from stage IVS and fatal stage IV tumors were compared. A SAGE tag (GCAACTTAAC) was detected that was overrepresented in stage IVS disease and that had no reliable match in current National Center for Biotechnology Information databases of SAGE profiles, thus pointing to a novel gene. The corresponding gene, which maps to chromosome 4q33-34, was identified using a modified 3′- and 5′-rapid amplification of cDNA ends PCR and was designated as DEIN (differentially expressed in neuroblastoma). The gene comprises five transcript variants and its sequence overlaps with expressed sequences of the yet uncharacterized UniGene cluster Hs.61435. DEIN exhibits nucleotide sequence conservation over a broad range of species with an overall homology of 65% between human and mouse. As none of the predicted amino acid sequences is homologous to known proteins, it remains to be determined whether DEIN represents a coding or noncoding RNA. Northern blot analysis and semiquantitative reverse transcription-PCR showed high DEIN expression in neuroblastoma, whereas expression was absent or weak in most normal adult tissues. Analysis of 121 primary neuroblastomas by real-time reverse transcription-PCR revealed a strong association with age at diagnosis <1 year and particularly with stage IVS disease (both P < 0.001). The characteristic expression pattern of DEIN suggests a specific role of this gene in the unique biology of stage IVS tumors and may help to molecularly define this special subtype of neuroblastoma. (Mol Cancer Res 2007;5(12):1276–84)