Search

Background Genomic alterations of the anaplastic lymphoma kinase gene (ALK) occur recurrently in neuroblastoma, a pediatric malignancy of the sympathetic nervous system. However, information on their development over time has remained sparse. Methods ALK alterations were assessed in neuroblastomas at diagnosis and/or relapse from a total of 943 patients, covering all stages of disease. Longitudinal information on diagnostic and relapsed samples from individual patients was available in 101 and 102 cases for mutation and amplification status, respectively. Results At diagnosis, ALK point mutations occurred in 10.5% of all cases, with highest frequencies in stage 4 patients ALK alteration frequency increased by 70%, both in high-risk and non-high-risk cases. The increase was most likely due to de novo mutations, frequently leading to R1275Q substitutions, which are sensitive to pharmacological ALK inhibition. By contrast, the frequency of ALK amplifications did not change over the course of the disease. ALK amplifications, but not mutations, were associated with poor patient outcome. Conclusions The considerably increased frequency of ALK mutations at relapse and their high prevalence in young stage 4 patients suggest surveying the genomic ALK status regularly in these patient cohorts, and to evaluate ALK-targeted treatment also in intermediate-risk patients.

Supplementary Data 1. This supplementary data comprises a step-by-step description on how our gene-expression based classification models were generated and validated

XLS file - 1949K

XLS file - 61K

Supplementary Figure 1. Highlights the EFS and OS for the cohort of neuroblastoma patients with MYCN-amplified disease (n=114, Fig 1a) and for th subcohort of patients >18 months of age with stage 4, MYCN non-amplified disease (n=102, fig. 1b). F, favorable; UF, unfavorable.

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

Supplementary Table 1. Clinical co-variates for the 709 patients who participated in the study. This supplementary table comprises detailed information on clinical co-variates for all 709 neuroblastoma patients who participated in the study.

XLS file - 10K

XLS file - 25K

Supplementary Table 2. External performance validation of the top five classifiers. This supplementary table comprises the external performance metrics of the top five classifiers. Indicated are the values for classification accuracy, sensitivity , specificity and Matthew's correlation coefficient (MCC) in the prediction of 325 patients of the test set who fulfilled the criteria for classifier training (Favorable, n=187; Unfavorable, n=138).

XLS file - 50K

Supplementary Table 4. Transcripts that contribute to the SVM_th10 classifier. This supplementary table highlights the transcripts that constitute the SVM_th10 classifier. Indicated for each feature are the Oligo-ID, the Gene symbols, the Entrez gene IDs, the RefSeq IDs, the Gene IDs and the Transcript IDs.

Legends to Supplementaries. This documents comprises the legends to the supplementary figures and tables

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.

Supplementary Data 2. This supplementary data comprises the r-algorithm scripts required to perform the both model selection and validation as described in the step-by-step protocol

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

Extrachromosomal circularization of DNA is an important genomic feature in cancer. The structure, composition and genome-wide frequency of extrachromosomal circular DNA, however, have not yet been extensively profiled. Here, we combined genomic and transcriptomic approaches to describe the landscape of extrachromosomal circular DNA in neuroblastoma, a tumor arising in childhood from primitive cells of the sympathetic nervous system. Our analysis identifies and characterizes a wide catalog of somatically acquired and undescribed extrachromosomal circular DNAs. Moreover, we find that extrachromosomal circular DNAs are an unanticipated major source of somatic rearrangements, contributing to oncogenic remodeling through chimeric circularization and reintegration of circular DNA into the linear genome. Cancer-causing lesions can emerge out of circle-derived rearrangements and are associated with adverse clinical outcome. It is highly probable that circle-derived rearrangements represent an ongoing mutagenic process. Thus, extrachromosomal circular DNAs represent a multi-hit mutagenic process, with important functional and clinical implications for the origins of genomic remodeling in cancer.

Background Long-term survival of high-risk neuroblastoma patients is still below 50% despite intensive multimodal treatment. This trial aimed to address whether the addition of two topotecan-containing chemotherapy courses compared to standard induction therapy improves event-free survival (EFS) of these patients. Patients and methods An open-label, multicenter, prospective randomized controlled trial was carried out at 58 hospitals in Germany and Switzerland. Patients aged 1–21 years with stage 4 neuroblastoma and patients aged 6 months to 21 years with MYCN-amplified tumors were eligible. The primary endpoint was EFS. Patients were randomly assigned to standard induction therapy with six chemotherapy courses or to experimental induction chemotherapy starting with two additional courses of topotecan, cyclophosphamide, and etoposide followed by standard induction chemotherapy (eight courses in total). After induction chemotherapy, all patients received high-dose chemotherapy with autologous hematopoietic stem cell rescue and isotretinoin for consolidation. Radiotherapy was applied to patients with active tumors at the end of induction chemotherapy. Results Of 536 patients enrolled in the trial, 422 were randomly assigned to the control arm (n = 211) and the experimental arm (n = 211); the median follow-up time was 3.32 years (interquartile range 1.65–5.92). At data lock, the 3-year EFS of experimental and control patients was 34% and 32% [95% confidence Interval (CI) 28% to 40% and 26% to 38%; P = 0.258], respectively. Similarly, the 3-year overall survival of the patients did not differ [54% and 48% (95% CI 46% to 62% and 40% to 56%), respectively; P = 0.558]. The response to induction chemotherapy was not different between the arms. The median number of non-fatal toxicities per patient was higher in the experimental group while the median number of toxicities per chemotherapy course was not different. Conclusion While the burden for the patients was increased by prolonging the induction chemotherapy and the toxicity, the addition of two topotecan-containing chemotherapy courses did not improve the EFS of high-risk neuroblastoma patients and thus cannot be recommended. Clinical Trials.gov number NCT number 03042429.

Background Neuroblastoma is characterized by substantial clinical heterogeneity. Despite intensive treatment, the survival rates of high-risk neuroblastoma patients are still disappointingly low. Somatic chromosomal copy number aberrations have been shown to be associated with patient outcome, particularly in low- and intermediate-risk neuroblastoma patients. To improve outcome prediction in high-risk neuroblastoma, we aimed to design a prognostic classification method based on copy number aberrations. Methods In an international collaboration, normalized high-resolution DNA copy number data (arrayCGH and SNP arrays) from 556 high-risk neuroblastomas obtained at diagnosis were collected from nine collaborative groups and segmented using the same method. We applied logistic and Cox proportional hazard regression to identify genomic aberrations associated with poor outcome. Results In this study, we identified two types of copy number aberrations that are associated with extremely poor outcome. Distal 6q losses were detected in 5.9% of patients and were associated with a 10-year survival probability of only 3.4% (95% confidence interval [CI] = 0.5% to 23.3%, two-sided P = .002). Amplifications of regions not encompassing the MYCN locus were detected in 18.1% of patients and were associated with a 10-year survival probability of only 5.8% (95% CI = 1.5% to 22.2%, two-sided P < .001). Conclusions Using a unique large copy number data set of high-risk neuroblastoma cases, we identified a small subset of high-risk neuroblastoma patients with extremely low survival probability that might be eligible for inclusion in clinical trials of new therapeutics. The amplicons may also nominate alternative treatments that target the amplified genes.

Neuroblastoma is a malignant paediatric tumour of the sympathetic nervous system1. 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 elusive2–4. 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 type1,2,5. 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.

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.

The prognostic relevance of chromosome 17 gain in neuroblastoma is still discussed. This investigation specifies the frequency, type, size, and transcriptional relevance in a large patient cohort. Primary tumor material of 202 patients was analyzed using high-resolution oligonucleotide array-based comparative genomic hybridization (aCGH) and correlated with clinical and survival data. A subset (n = 145) was correlated for differentially expressed genes (DEG) by microarray analysis. Chromosome 17 aCGH analysis showed numerical gain in 94/202 patients (47%), partial gain in 93/202 patients (46%), and no gain in 15/202 patients (7%). The frequency of partial gain was higher in stage 4 neuroblastoma (stage 1 15%; stage 2 12%; stage 3 16%; stage 4S 7%; and stage 4 50%). Overall survival (OS) was superior in patients with numerical gain compared with patients with partial gain or no gain (5-y-OS: 0.95 ± 0.02 vs. 0.63 ± 0.05 vs. 0.60 ± 0.13; P < 0.001). Gene expression analysis demonstrated 95/130 DEGs between tumors with numerical or partial chromosome/no gain. Only one DEG (CCKBR) was detected comparing tumors with partial gain and those with no gain. In patients with partial gain, the distribution of breakpoints did not correlate with stage and 11q status, but with MYCN amplification and 1p status. The "best" breakpoints in cases with partial 17q gain were at 42.5 Mb for event-free and 26.6 Mb for OS. Numerical gain of chromosome 17 is associated with a better prognosis than partial and no gain. The group of tumors with partial gain was similar to the group without gain with respect to stage distribution, outcome, and gene expression profile.

Motivation: Recurrent DNA breakpoints in cancer genomes indicate the presence of critical functional elements for tumor development. Identifying them can help determine new therapeutic targets. High-dimensional DNA microarray experiments like arrayCGH afford the identification of DNA copy number breakpoints with high precision, offering a solid basis for computational estimation of recurrent breakpoint locations. Results: We introduce a method for identification of recurrent breakpoints (consensus breakpoints) from copy number aberration datasets. The method is based on weighted kernel counting of breakpoints around genomic locations. Counts larger than expected by chance are considered significant. We show that the consensus breakpoints facilitate consensus segmentation of the samples. We apply our method to three arrayCGH datasets and show that by using consensus segmentation we achieve significant dimension reduction, which is useful for the task of prediction of tumor phenotype based on copy number data. We use our approach for classification of neuroblastoma tumors from different age groups and confirm the recent recommendation for the choice of age cut-off for differential treatment of 18 months. We also investigate the (epi)genetic properties at consensus breakpoint locations for seven datasets and show enrichment in overlap with important functional genomic regions. Availability: Implementation in R of our approach can be found at http://www.mpi-inf.mpg.de/ ∼laura/FeatureGrouping.html. Contact: laura@mpi-inf.mpg.de. Supplementary information: Supplementary data are available at Bioinformatics online.

The TP53 tumor suppressor pathway is abrogated by TP53 mutations in the majority of human cancers. Increased levels of wild-type TP53 in aggressive neuroblastomas appear paradox but are tolerated by tumor cells due to co-activation of the TP53 ubiquitin ligase, MDM2. The role of the MDM2 antagonist, p14(ARF), in controlling the TP53-MDM2 balance in neuroblastoma is unresolved. In the present study, we show that conditional p14(ARF) expression substantially suppresses viability, clonogenicity and anchorage-independent growth in p14(ARF)-deficient or MYCN-amplified neuroblastoma cell lines. Furthermore, ectopic 14(ARF) expression induced accumulation of cells in the G1 phase and apoptosis, which was paralleled by accumulation of TP53 and its targets. Comparative genomic hybridization analysis of 193 primary neuroblastomas detected one homozygous deletion of CDKN2A (encoding both p14(ARF) and p16(INK4A)) and heterozygous loss of CDKN2A in 22% of tumors. Co-expression analysis of p14(ARF) and its transactivator, E2F1, in a set of 68 primary tumors revealed only a weak correlation, suggesting that further regulatory mechanisms govern p14(ARF) expression in neuroblastomas. Intriguingly, analyses utilizing chromatin immunoprecipitation revealed different histone mark-defined epigenetic activity states of p14(ARF) in neuroblastoma cell lines that correlated with endogenous p14(ARF) expression but not with episomal p14(ARF) promoter reporter activity, indicating that the native chromatin context serves to epigenetically repress p14(ARF) in neuroblastoma cells. Collectively, the data pinpoint p14(ARF) as a critical factor for efficient TP53 response in neuroblastoma cells and assign p14(ARF) as a neuroblastoma suppressor candidate that is impaired by genomic loss and epigenetic repression.

Poor prognosis in neuroblastoma is associated with genetic amplification of MYCN. MYCN is itself a target of let-7, a tumour suppressor family of microRNAs implicated in numerous cancers. LIN28B, an inhibitor of let-7 biogenesis, is overexpressed in neuroblastoma and has been reported to regulate MYCN. Here we show, however, that LIN28B is dispensable in MYCN-amplified neuroblastoma cell lines, despite de-repression of let-7. We further demonstrate that MYCN messenger RNA levels in amplified disease are exceptionally high and sufficient to sponge let-7, which reconciles the dispensability of LIN28B. We found that genetic loss of let-7 is common in neuroblastoma, inversely associated with MYCN amplification, and independently associated with poor outcomes, providing a rationale for chromosomal loss patterns in neuroblastoma. We propose that let-7 disruption by LIN28B, MYCN sponging, or genetic loss is a unifying mechanism of neuroblastoma development with broad implications for cancer pathogenesis.

International audience; About 50% of children with neuroblastoma (NB) show a metastatic disease and have a poor prognosis. However, disease progression is greatly variable and depends on patients' age and MYCN oncogene amplification. To investigate the role of patients' age in tumor aggressiveness, we performed array-CGH and gene expression profiles of three groups (G) of metastatic NBs: G1, stage 4S patients and MYCN single copy (MYCN−) tumors; G2, stage 4 patients, ≤18 months of age, MYCN− tumors and favorable outcome and G3, Stage 4 patients, ≥19 months with unfavorable outcome. G1 was characterized by numerical aberrations prevalently; on the contrary, all G3 tumors had structural rearrangements, whereas G2 showed an intermediate pattern. The average of numerical alterations decreased significantly from G1 to G2 to G3 (p < 0.01). Contrarily, the number of structural aberrations increased from G1 to G2 to G3 (p < 2.35 E−05). Noteworthy, G3/MYCN− NBs were characterized by several complex intrachromosome rearrangements. Expression analysis of the three groups showed significant differences in genes of Rho and Ras signaling pathways, development and adhesion, cell cycle regulation and telomerase activity. Accumulation of structural alterations increased with patients' age and was associated with a more aggressive disease. Abnormal expression of genes involved in cell cycle and telomerase in G3 may be responsible for the genomic instability in this cohort of patients. The higher DNA instability observed in G3/MYCN− NBs than in MYCN-amplified G3 may also explain why patients ≥19 months have a poor outcome independently by MYCN status.

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 ≤ .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.

Background Gene expression profiling is being widely applied in cancer research to identify biomarkers for clinical endpoint prediction. Since RNA-seq provides a powerful tool for transcriptome-based applications beyond the limitations of microarrays, we sought to systematically evaluate the performance of RNA-seq-based and microarray-based classifiers in this MAQC-III/SEQC study for clinical endpoint prediction using neuroblastoma as a model. Results We generate gene expression profiles from 498 primary neuroblastomas using both RNA-seq and 44 k microarrays. Characterization of the neuroblastoma transcriptome by RNA-seq reveals that more than 48,000 genes and 200,000 transcripts are being expressed in this malignancy. We also find that RNA-seq provides much more detailed information on specific transcript expression patterns in clinico-genetic neuroblastoma subgroups than microarrays. To systematically compare the power of RNA-seq and microarray-based models in predicting clinical endpoints, we divide the cohort randomly into training and validation sets and develop 360 predictive models on six clinical endpoints of varying predictability. Evaluation of factors potentially affecting model performances reveals that prediction accuracies are most strongly influenced by the nature of the clinical endpoint, whereas technological platforms (RNA-seq vs. microarrays), RNA-seq data analysis pipelines, and feature levels (gene vs. transcript vs. exon-junction level) do not significantly affect performances of the models. Conclusions We demonstrate that RNA-seq outperforms microarrays in determining the transcriptomic characteristics of cancer, while RNA-seq and microarray-based models perform similarly in clinical endpoint prediction. Our findings may be valuable to guide future studies on the development of gene expression-based predictive models and their implementation in clinical practice. Electronic supplementary material The online version of this article (doi:10.1186/s13059-015-0694-1) contains supplementary material, which is available to authorized users.

Background Segmental genomic copy number alterations, such as loss of 11q or 3p and gain of 17q, are well established markers of poor outcome in neuroblastoma, and have been suggested to comprise tumor suppressor genes or oncogenes, respectively. The gene forkhead box P1 (FOXP1) maps to chromosome 3p14.1, a tumor suppressor locus deleted in many human cancers including neuroblastoma. FoxP1 belongs to a family of winged-helix transcription factors that are involved in processes of cellular proliferation, differentiation and neoplastic transformation. Methods Microarray expression profiles of 476 neuroblastoma specimens were generated and genes differentially expressed between favorable and unfavorable neuroblastoma were identified. FOXP1 expression was correlated to clinical markers and patient outcome. To determine whether hypermethylation is involved in silencing of FOXP1, methylation analysis of the 5′ region of FOXP1 in 47 neuroblastomas was performed. Furthermore, FOXP1 was re-expressed in three neuroblastoma cell lines to study the effect of FOXP1 on growth characteristics of neuroblastoma cells. Results Low expression of FOXP1 is associated with markers of unfavorable prognosis like stage 4, age >18 months and MYCN amplification and unfavorable gene expression-based classification (P

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

The prognostic relevance of chromosome 17 gain in neuroblastoma is still discussed. This investigation specifies the frequency, type, size, and transcriptional relevance in a large patient cohort. Primary tumor material of 202 patients was analyzed using high-resolution oligonucleotide array-based comparative genomic hybridization (aCGH) and correlated with clinical and survival data. A subset (n = 145) was correlated for differentially expressed genes (DEG) by microarray analysis. Chromosome 17 aCGH analysis showed numerical gain in 94/202 patients (47%), partial gain in 93/202 patients (46%), and no gain in 15/202 patients (7%). The frequency of partial gain was higher in stage 4 neuroblastoma (stage 1 15%; stage 2 12%; stage 3 16%; stage 4S 7%; and stage 4 50%). Overall survival (OS) was superior in patients with numerical gain compared with patients with partial gain or no gain (5-y-OS: 0.95 ± 0.02 vs. 0.63 ± 0.05 vs. 0.60 ± 0.13; P0.001). Gene expression analysis demonstrated 95/130 DEGs between tumors with numerical or partial chromosome/no gain. Only one DEG (CCKBR) was detected comparing tumors with partial gain and those with no gain. In patients with partial gain, the distribution of breakpoints did not correlate with stage and 11q status, but with MYCN amplification and 1p status. The "best" breakpoints in cases with partial 17q gain were at 42.5 Mb for event-free and 26.6 Mb for OS. Numerical gain of chromosome 17 is associated with a better prognosis than partial and no gain. The group of tumors with partial gain was similar to the group without gain with respect to stage distribution, outcome, and gene expression profile.

The let-7 microRNA family are known tumor suppressors often deregulated in cancer, yet the underlying mechanisms of let-7 disruption remain poorly understood. Neuroblastoma, a neural crest derived tumor, is defined in part by poor prognosis associated with genetic amplification of MYCN, itself a let-7 target. The let-7 biogenesis inhibitor LIN28B has recently been implicated as a critical regulator of MYCN, but through CRISPR-mediated gene disruption we show that LIN28B is dispensable for both MYCN protein expression and growth of MYCN-amplified neuroblastoma cell lines despite robust de-repression of let-7, prompting us to explore additional mechanisms for let-7 disruption. Consequently, we have found a novel non-coding role for amplified MYCN mRNA as a potent let-7 sponge that through exceptionally high expression defines a sub-class of self-sponging amplified-competing-endogenous-RNA (aceRNA) and reconciling the dispensability of LIN28B in neuroblastoma cell lines. Furthermore, by analyzing a large cohort of tumor samples from patients, we observe frequent genomic loss of let-7 that inversely associates with MYCN-amplification, providing a functional explanation for the known MYCN-amplification-independent pattern of chromosome 3p and 11q loss, which harbor let-7g and let-7a2, respectively. We thus propose a model whereby let-7 disruption by genetic loss, LIN28B expression, or aceRNA sponging is a unifying mechanism of neuroblastoma pathogenesis. Indeed, our data show that the majority of neuroblastomas have at least one let-7 disruption event and that genetic loss in non-MYCN-amplified tumors marks decreased survival, further underscoring its importance. The inverse selective relationship between allelic loss and sponging of let-7 from highly expressed or amplified oncogenes may have broad implications for oncogenesis. Citation Format: John T. Powers, Kaloyan Tsanov, Frederik Roels, Catherine Spina, Richard Ebright, Marc Seligson, Yvanka de Soysa, Patrick Cahan, Daniel Pearson, Jessica Theißen, Ho-Chou Tu, Grace LaPier, Jihan Osborne, Samantha Ross, James Collins, Frank Berthold, George Daley. Multiple mechanisms disrupt let-7 miRNA biogenesis and function in neuroblastoma. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr LB-165.

Several microRNA (miRNA) loci are found within genomic regions frequently deleted in primary neuroblastoma, including miR-885-5p at 3p25.3. In this study, we demonstrate that miR-885-5p is downregulated on loss of 3p25.3 region in neuroblastoma. Experimentally enforced miR-885-5p expression in neuroblastoma cell lines inhibits proliferation triggering cell cycle arrest, senescence and/or apoptosis. miR-885-5p leads to the accumulation of p53 protein and activates the p53 pathway, resulting in upregulation of p53 targets. Enforced miR-885-5p expression consistently leads to downregulation of cyclin-dependent kinase (CDK2) and mini-chromosome maintenance protein (MCM5). Both genes are targeted by miR-885-5p via predicted binding sites within the 3′-untranslated regions (UTRs) of CDK2 and MCM5. Transcript profiling after miR-885-5p introduction in neuroblastoma cells reveals alterations in expression of multiple genes, including several p53 target genes and a number of factors involved in p53 pathway activity. Taken together, these data provide evidence that miR-885-5p has a tumor suppressive role in neuroblastoma interfering with cell cycle progression and cell survival.

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.

Imbalances in chromosome 11q occur in approximately 30% of primary neuroblastoma and are associated with poor outcome. It has been suggested that 11q loss constitutes a distinct clinico-genetic neuroblastoma subgroup by affecting expression levels of corresponding genes. This study analysed the relationship of 11q loss, clinical phenotype and global transcriptomic profiles in four clinico-genetic subgroups (11q alteration/favourable outcome, n=7; 11q alteration/unfavourable outcome, n=14; no 11q alteration/favourable outcome, n=81; no 11q alteration/unfavourable outcome, n=8; tumours with MYCN amplification and/or 1p loss were excluded). Unsupervised and supervised comparisons of gene expression profiles consistently showed significantly different mRNA patterns between favourable and unfavourable neuroblastomas, both in the subgroups with and without 11q loss. In contrast, favourable tumours with and without 11q loss showed highly similar transcriptomic profiles. Disproportionate downregulation of 11q genes was observed only in unfavourable tumours with 11q loss. The diverging molecular profiles were neither caused by considerable differences in the size of the deleted regions nor by differential methylation patterns of 11q genes. Together, this study shows that neuroblastoma with 11q loss comprises two biological subgroups that differ both in their clinical phenotype and gene expression patterns, indicating that 11q loss is not a primary determinant of neuroblastoma tumour behaviour.

For many decades, neuroblastoma has remained a challenging disease for both clinicians and researchers. Now, techniques that efficiently specify both comprehensive genetic and gene-expression alterations of neuroblastoma tumors have provided molecular markers that indicate tumor behavior and patient outcome with very high accuracy. Once the anticipated value of these markers has been confirmed in ongoing studies, patients may profit from more accurate risk assessment by integrating these markers into clinical routine. Moreover, disclosing further tumor-initiating events, such as the recently revealed oncogenic mutations of ALK, will further promote the elucidation of the genetic etiology of the disease. Together with recent information on altered signaling pathways in aggressively growing tumors, this knowledge will help to establish therapeutic strategies specifically targeting molecular key factors of neuroblastoma tumor progression.

Purpose: MYCN amplification is an important therapy-stratifying marker in neuroblastoma. Fluorescence in situ hybridization with signal detection on the single-cell level allows a critical judgement of MYCN intratumoral heterogeneity. Experimental Design: The MYCN status was investigated by fluorescence in situ hybridization at diagnosis and relapse. Heterogeneity was defined as the simultaneous presence of amplified cells (≥5 cells per slide) and nonamplified cells within one tumor or sequential change of the amplification status during the course of the disease. Likewise, heterogeneity can be detected between primary tumor and metastasis. Results: From 1,341 patients analyzed, 1,071 showed no amplification, 250 showed homogeneous amplification, and 20 patients showed MYCN heterogeneity. Of the patients with heterogeneity, 12 of 20 had clusters of MYCN amplifications, 3 of 20 had amplified single cells, 3 of 20 showed MYCN amplifications in the bone marrow but not in the primary tumor, and 2 of 20 acquired MYCN amplification during the course of the disease. All stage 4 patients were treated according to high-risk protocols; 7 of 8 later progressed. Four patients with localized disease were treated according to high-risk protocol because of MYCN-amplified clusters; 1 of 4 later progressed. One patient treated with mild chemotherapy experienced progression. Seven patients with localized/4S disease underwent no chemotherapy: 4 of 5 patients with MYCN heterogeneity at diagnosis remained disease-free, and 1 of 5 experienced local progression. Two patients had normal MYCN status at diagnosis but acquired MYCN amplification during the course of the disease. Conclusion: MYCN heterogeneity is rare. Our results suggest that small amounts of MYCN-amplified cells are not correlated to adverse outcomes. More patients with heterogeneity are warranted to clarify the role of MYCN heterogeneity for risk classification.

Neuroblastoma (NB), the most common solid tumor of childhood, is characterized by a remarkable heterogeneity of patients' courses and survival rates of older patients with metastatic disease have remained poor. Numerous prognostic factors including amplification of the MYCN oncogene have been described. However, despite advanced clinical risk stratification, current trials still fail to determine the best treatment strategy for a substantial number of patients. On the route to clinical application of improved array-based classifiers to predict the risk profile of individual patients, significant progress has been achieved over the past years. We here present data on a cohort of 138 primary NB using a novel approach including information for all human coding exons described to date (Affymetrix ExonST Array). Using a classifier trained on 100 patient samples and then used to predict the outcome of the remaining 38 patients, we were able to achieve prediction accuracies >80% in the independent test set using support vector machine (SVM) learning algorithms, which is superior to the current clinical risk stratification. Interestingly, the histone demethylase JARID1C, which had not been linked to cancer formation previously, was up-regulated in relapsing NB and found to have prognostic value independent of MYCN amplification. JARID1C was also highly expressed in all NB cell lines investigated. Down-regulation of JARID1C using siRNA inhibited cell proliferation in vitro. Taken together, exon level analysis appears to be a highly promising tool both for prediction of outcome and for providing new insights into tumor biology.

The let-7 microRNA family are known tumor suppressors often deregulated in cancer, yet the underlying mechanisms of let-7 disruption remain poorly understood. Neuroblastoma, a neural crest derived tumor, is defined in part by poor prognosis associated with genetic amplification of MYCN, itself a let-7 target. The let-7 biogenesis inhibitor LIN28B has recently been implicated as a critical regulator of MYCN, but we have employed siRNA and CRISPR-mediated gene disruption to show that LIN28B is dispensable for both MYCN protein expression and growth of MYCN-amplified neuroblastoma cell lines despite robust de-repression of let-7, which prompted us to explore additional mechanisms for let-7 disruption. Consequently, we have found that amplified MYCN mRNA is a potent let-7 sponge that through exceptionally high expression defines a sub-class of self-sponging amplified-competing-endogenous-RNA (aceRNA), which reconciles the dispensability of LIN28B in NB cell lines. In addition, by analyzing a large cohort of tumor samples from patients, we observe frequent genomic loss of let-7 that inversely associates with MYCN-amplification, providing a functional explanation for the known MYCN-amplification-independent pattern of chromosome 3p and 11q loss, which harbor let-7g and let-7a2, respectively. We thus propose a model whereby let-7 disruption by genetic loss, LIN28B expression, or aceRNA sponging is a unifying mechanism of neuroblastoma pathogenesis. Indeed, our data show that the majority of neuroblastomas have at least one let-7 disruption event and that genetic loss in non-MYCN-amplified tumors marks decreased survival, further underscoring its importance. The inverse selective relationship between allelic loss and sponging of let-7 from highly expressed or amplified oncogenes may have broad implications for oncogenesis. Note: This abstract was not presented at the meeting. Citation Format: John T. Powers, Kaloyan M. Tsanov, Frederik Roles, Richard Ebright, Marc Seligson, Yvanka de Soysa, Patrick Cahan, Jessica Theissen, Grace S. LaPier, Dan S. Pearson, Frank Berthold, George Q. Daley. Multiple distinct mechanisms disrupt let-7 miRNA biogenesis and function in neuroblastoma. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr LB-290. doi:10.1158/1538-7445.AM2015-LB-290

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.

International audience; We aimed to identify novel molecular prognostic markers to better predict relapse risk estimate for children with high-risk (HR) metastatic neuroblastoma (NB). We performed genome- and/or transcriptome-wide analyses of 129 stage 4 HR NBs. Children older than 1 year of age were categorized as “short survivors” (dead of disease within 5 years from diagnosis) and “long survivors” (alive with an overall survival time ≥ 5 years). We reported that patients with less than three segmental copy number aberrations in their tumor represent a molecularly defined subgroup with a high survival probability within the current HR group of patients. The complex genomic pattern is a prognostic marker independent of NB-associated chromosomal aberrations, i.e., MYCN amplification, 1p and 11q losses, and 17q gain. Integrative analysis of genomic and expression signatures demonstrated that fatal outcome is mainly associated with loss of cell cycle control and deregulation of Rho guanosine triphosphates (GTPases) functioning in neuritogenesis. Tumors with MYCN amplification show a lower chromosome instability compared to MYCN single-copy NBs (P = .0008), dominated by 17q gain and 1p loss. Moreover, our results suggest that the MYCN amplification mainly drives disruption of neuronal differentiation and reduction of cell adhesion process involved in tumor invasion and metastasis. Further validation studies are warranted to establish this as a risk stratification for patients.