Reliable assessment of telomere maintenance mechanisms in neuroblastoma.

BACKGROUND: Telomere maintenance mechanisms (TMM) are a hallmark of high-risk neuroblastoma, and are conferred by activation of telomerase or alternative lengthening of telomeres (ALT). However, detection of TMM is not yet part of the clinical routine, and consensus on TMM detection, especially on ALT assessment, remains to be achieved. METHODS: Whole genome sequencing (WGS) data of 68 primary neuroblastoma samples were analyzed. Telomere length was calculated from WGS data or by telomere restriction fragment analysis (n = 39). ALT was assessed by C-circle assay (CCA, n = 67) and detection of ALT-associated PML nuclear bodies (APB) by combined fluorescence in situ hybridization and immunofluorescence staining (n = 68). RNA sequencing was performed (n = 64) to determine expression of TERT and telomeric long non-coding RNA (TERRA). Telomerase activity was examined by telomerase repeat amplification protocol (TRAP, n = 15). RESULTS: Tumors were considered as telomerase-positive if they harbored a TERT rearrangement, MYCN amplification or high TERT expression (45.6%, 31/68), and ALT-positive if they were positive for APB and CCA (19.1%, 13/68). If all these markers were absent, tumors were considered TMM-negative (25.0%, 17/68). According to these criteria, the majority of samples were classified unambiguously (89.7%, 61/68). Assessment of additional ALT-associated parameters clarified the TMM status of the remaining seven cases with high likelihood: ALT-positive tumors had higher TERRA expression, longer telomeres, more telomere insertions, a characteristic pattern of telomere variant repeats, and were associated with ATRX mutations. CONCLUSIONS: We here propose a workflow to reliably detect TMM in neuroblastoma. We show that unambiguous classification is feasible following a stepwise approach that determines both, activation of telomerase and ALT. The workflow proposed in this study can be used in clinical routine and provides a framework to systematically and reliably determine telomere maintenance mechanisms for risk stratification and treatment allocation of neuroblastoma patients.

Telomerase-targeting compounds Imetelstat and 6-thio-dG act synergistically with chemotherapy in high-risk neuroblastoma models.

BACKGROUND: The majority of high-risk neuroblastomas harbor telomerase activity, and telomerase-interacting compounds, such as 6-thio-2'-deoxyguanosine (6-thio-dG), have been found to impair the growth of telomerase-positive neuroblastoma cell lines. It has remained unclear, however, how such drugs can be combined with other compounds used in current treatment concepts for neuroblastoma patients. METHODS: Growth-inhibitory effects of varying concentrations of 6-thio-dG in combination with etoposide, doxorubicin or ceritinib were determined in eight telomerase-positive neuroblastoma cell lines with distinct genetic backgrounds. Tumor growth inhibition of subcutaneous xenografts from three different cell lines was assessed upon treatment with 6-thio-dG, the competitive telomerase inhibitor imetelstat, etoposide, or combinations of these compounds. RESULTS: Robust synergistic anti-tumor effects were observed for combinations of 6-thio-dG and etoposide or doxorubicin, but not for 6-thio-dG and ceritinib, in telomerase-positive neuroblastoma cell lines in vitro. Treatment of mouse xenografts with combinations of 6-thio-dG and etoposide significantly attenuated tumor growth and improved mouse survival over etoposide alone in two of three cell line models. Treatment of xenograft tumors by imetelstat monotherapy decreased telomerase activity by roughly 50% and significantly improved survival over control in all three models, whereas treatment with imetelstat plus etoposide led to enhanced survival over etoposide monotherapy in one model. Mechanistically, the synergistic effect was found to be due to both increased apoptosis and cell cycle arrest. CONCLUSION: Our study indicates that telomerase is an actionable target in telomerase-positive neuroblastoma, and demonstrates that combination therapies including telomerase-interacting compounds may improve the efficacy of established cytotoxic drugs. Targeting telomerase may thus represent a therapeutic option in high-risk neuroblastoma patients.

Chromothripsis followed by circular recombination drives oncogene amplification in human cancer.

The mechanisms behind the evolution of complex genomic amplifications in cancer have remained largely unclear. Using whole-genome sequencing data of the pediatric tumor neuroblastoma, we here identified a type of amplification, termed 'seismic amplification', that is characterized by multiple rearrangements and discontinuous copy number levels. Overall, seismic amplifications occurred in 9.9% (274 of 2,756) of cases across 38 cancer types, and were associated with massively increased copy numbers and elevated oncogene expression. Reconstruction of the development of seismic amplification showed a stepwise evolution, starting with a chromothripsis event, followed by formation of circular extrachromosomal DNA that subsequently underwent repetitive rounds of circular recombination. The resulting amplicons persisted as extrachromosomal DNA circles or had reintegrated into the genome in overt tumors. Together, our data indicate that the sequential occurrence of chromothripsis and circular recombination drives oncogene amplification and overexpression in a substantial fraction of human malignancies.

Telomerase Is a Prognostic Marker of Poor Outcome and a Therapeutic Target in Neuroblastoma.

PURPOSE: Telomere maintenance is a hallmark of high-risk neuroblastoma; however, the contribution of telomerase and alternative lengthening of telomeres (ALT) to clinical phenotypes has remained unclear. We aimed to determine the clinical relevance of telomerase activation versus ALT as biomarkers in pretreatment neuroblastoma and to assess the potential value of telomerase as a therapeutic target. MATERIALS AND METHODS: The genomic status of TERT and MYCN was assessed in 457 pretreatment neuroblastomas by fluorescence in situ hybridization. ALT was examined in 273 of 457 tumors by detection of ALT-associated promyelocytic leukemia nuclear bodies, and TERT expression was determined by microarrays in 223 of these. Cytotoxic effects of telomerase-interacting compounds were analyzed in neuroblastoma cell lines in vitro and in vivo. RESULTS: We detected TERT rearrangements in 46 of 457 cases (10.1%), MYCN amplification in 93 of 457 cases (20.4%), and elevated TERT expression in tumors lacking TERT or MYCN alterations in 10 of 223 cases (4.5%). ALT activation was found in 49 of 273 cases (17.9%). All these alterations occurred almost mutually exclusively and were associated with unfavorable prognostic variables and adverse outcome. The presence of activated telomerase (ie, TERT rearrangements, MYCN amplification, or high TERT expression without these alterations) was associated with poorest overall survival and was an independent prognostic marker in multivariable analyses. We also found that the telomerase-interacting compound 6-thio-2'-deoxyguanosine effectively inhibited viability and proliferation of neuroblastoma cells bearing activated telomerase. Similarly, tumor growth was strongly impaired upon 6-thio-2'-deoxyguanosine treatment in telomerase-positive neuroblastoma xenografts in mice. CONCLUSION: Our data suggest telomerase activation and ALT define distinct neuroblastoma subgroups with adverse outcome and that telomerase may represent a promising therapeutic target in many high-risk neuroblastomas.

A mechanistic classification of clinical phenotypes in neuroblastoma.

Neuroblastoma is a pediatric tumor of the sympathetic nervous system. Its clinical course ranges from spontaneous tumor regression to fatal progression. To investigate the molecular features of the divergent tumor subtypes, we performed genome sequencing on 416 pretreatment neuroblastomas and assessed telomere maintenance mechanisms in 208 of these tumors. We found that patients whose tumors lacked telomere maintenance mechanisms had an excellent prognosis, whereas the prognosis of patients whose tumors harbored telomere maintenance mechanisms was substantially worse. Survival rates were lowest for neuroblastoma patients whose tumors harbored telomere maintenance mechanisms in combination with RAS and/or p53 pathway mutations. Spontaneous tumor regression occurred both in the presence and absence of these mutations in patients with telomere maintenance-negative tumors. On the basis of these data, we propose a mechanistic classification of neuroblastoma that may benefit the clinical management of patients.

Molecular Classification Substitutes for the Prognostic Variables Stage, Age, and MYCN Status in Neuroblastoma Risk Assessment.

BACKGROUND: Current risk stratification systems for neuroblastoma patients consider clinical, histopathological, and genetic variables, and additional prognostic markers have been proposed in recent years. We here sought to select highly informative covariates in a multistep strategy based on consecutive Cox regression models, resulting in a risk score that integrates hazard ratios of prognostic variables. METHODS: A cohort of 695 neuroblastoma patients was divided into a discovery set (n=75) for multigene predictor generation, a training set (n=411) for risk score development, and a validation set (n=209). Relevant prognostic variables were identified by stepwise multivariable L1-penalized least absolute shrinkage and selection operator (LASSO) Cox regression, followed by backward selection in multivariable Cox regression, and then integrated into a novel risk score. RESULTS: The variables stage, age, MYCN status, and two multigene predictors, NB-th24 and NB-th44, were selected as independent prognostic markers by LASSO Cox regression analysis. Following backward selection, only the multigene predictors were retained in the final model. Integration of these classifiers in a risk scoring system distinguished three patient subgroups that differed substantially in their outcome. The scoring system discriminated patients with diverging outcome in the validation cohort (5-year event-free survival, 84.9±3.4 vs 63.6±14.5 vs 31.0±5.4; P<.001), and its prognostic value was validated by multivariable analysis. CONCLUSION: We here propose a translational strategy for developing risk assessment systems based on hazard ratios of relevant prognostic variables. Our final neuroblastoma risk score comprised two multigene predictors only, supporting the notion that molecular properties of the tumor cells strongly impact clinical courses of neuroblastoma patients.

Transcription factor activating protein 2 beta (TFAP2B) mediates noradrenergic neuronal differentiation in neuroblastoma.

Neuroblastoma is an embryonal pediatric tumor that originates from the developing sympathetic nervous system and shows a broad range of clinical behavior, ranging from fatal progression to differentiation into benign ganglioneuroma. In experimental neuroblastoma systems, retinoic acid (RA) effectively induces neuronal differentiation, and RA treatment has been therefore integrated in current therapies. However, the molecular mechanisms underlying differentiation are still poorly understood. We here investigated the role of transcription factor activating protein 2 beta (TFAP2B), a key factor in sympathetic nervous system development, in neuroblastoma pathogenesis and differentiation. Microarray analyses of primary neuroblastomas (n = 649) demonstrated that low TFAP2B expression was significantly associated with unfavorable prognostic markers as well as adverse patient outcome. We also found that low TFAP2B expression was strongly associated with CpG methylation of the TFAP2B locus in primary neuroblastomas (n = 105) and demethylation with 5-aza-2'-deoxycytidine resulted in induction of TFAP2B expression in vitro, suggesting that TFAP2B is silenced by genomic methylation. Tetracycline inducible re-expression of TFAP2B in IMR-32 and SH-EP neuroblastoma cells significantly impaired proliferation and cell cycle progression. In IMR-32 cells, TFAP2B induced neuronal differentiation, which was accompanied by up-regulation of the catecholamine biosynthesizing enzyme genes DBH and TH, and down-regulation of MYCN and REST, a master repressor of neuronal genes. By contrast, knockdown of TFAP2B by lentiviral transduction of shRNAs abrogated RA-induced neuronal differentiation of SH-SY5Y and SK-N-BE(2)c neuroblastoma cells almost completely. Taken together, our results suggest that TFAP2B is playing a vital role in retaining RA responsiveness and mediating noradrenergic neuronal differentiation in neuroblastoma.

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.

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.

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.

FOXP1 inhibits cell growth and attenuates tumorigenicity of neuroblastoma.

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 < 0.001 each). Moreover, FOXP1 expression predicts patient outcome accurately and independently from well-established prognostic markers. Array-based CGH analysis of 159 neuroblastomas revealed that heterozygous loss of the FOXP1 locus was a rare event (n = 4), but if present, was associated with low FOXP1 expression. By contrast, DNA methylation analysis in 47 neuroblastomas indicated that hypermethylation is not regularly involved in FOXP1 gene silencing. Re-expression of FoxP1 significantly impaired cell proliferation, viability and colony formation in soft agar. Furthermore, induction of FOXP1 expression led to cell cycle arrest and apoptotic cell death of neuroblastoma cells. CONCLUSIONS: Our results suggest that down-regulation of FOXP1 expression is a common event in high-risk neuroblastoma pathogenesis and may contribute to tumor progression and unfavorable patient outcome.

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.

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

Co-regulated expression of HAND2 and DEIN by a bidirectional promoter with asymmetrical activity in neuroblastoma.

BACKGROUND: HAND2, a key regulator for the development of the sympathetic nervous system, is located on chromosome 4q33 in a head-to-head orientation with DEIN, a recently identified novel gene with stage specific expression in primary neuroblastoma (NB). Both genes are expressed in primary NB as well as most NB cell lines and are separated by a genomic sequence of 228 bp. The similar expression profile of both genes suggests a common transcriptional regulation mediated by a bidirectional promoter. RESULTS: Northern Blot analysis of DEIN and HAND2 in 20 primary NBs indicated concurrent expression levels of the two genes, which was confirmed by microarray analysis of 236 primary NBs (Pearson's correlation coefficient r = 0.65). While DEIN expression in the latter cohort was associated with stage 4S (p = 0.02), HAND2 expression was not associated with tumor stage. In contrast, both HAND2 and DEIN transcript levels were highly associated with age at diagnosis <12 months (p = 0.001). The intergenic region shows substantial homology in different species (89%, 72% and 53% identity between human and mouse, chicken and zebrafish, respectively) and contains many highly conserved putative transcription factor binding sites. Using luciferase reporter gene constructs, asymmetrical bidirectional promoter activity was found in four NB cell lines: In DEIN orientation, an average 3.4 fold increase in activity was observed as compared to the promoterless vector, whereas an average 15.4 fold activation was detected in HAND2 orientation. The presence of two highly conserved putative regulatory elements, one of which was shown to enhance HAND2 expression in branchial arches previously, displayed weak repressor activity for both genes. CONCLUSION: HAND2 and DEIN represent a gene pair that is tightly linked by a bidirectional promoter in an evolutionary highly conserved manner. Expression of both genes in NB is co-regulated by asymmetrical activity of this promoter and modulated by the activity of two cis-regulatory elements acting as weak repressors. The concurrent quantitative and tissue specific expression of HAND2 and DEIN suggests a functional link between both genes.

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.

Classification of neuroblastoma patients by published gene-expression markers reveals a low sensitivity for unfavorable courses of MYCN non-amplified disease.

Currently, Pubmed lists 385 marker genes for neuroblastoma outcome. Using a customized neuroblastoma-microarray, we evaluated the prognostic impact of the gene-expression pattern of 349 of these candidates (90.6%) in 127 neuroblastoma patients with divergent outcome. By significance analysis of microarrays (SAM) and both uncorrected and Bonferroni-corrected ANOVA, 166/349 (47.5%), 218/349 (62.5%) and 128/349 (36.4%) candidates showed significant differential expression between patients with contrasting outcome. By Prediction Analysis for Microarrays (PAM), a 38-gene-classifier was derived from all markers, which classified patients outcome with an overall accuracy of 78.5%. However, patients with unfavorable outcome of MYCN non-amplified disease were largely misclassified (accuracy: 35%), suggesting that these courses are not identified by current marker genes.

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.

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 classifier's 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.

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 approximately 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

Characterization of a complex genomic alteration on chromosome 2p that leads to four alternatively spliced fusion transcripts in the neuroblastoma cell lines IMR-5, IMR-5/75 and IMR-32.

Genetic aberrations in neuroblastoma (NB) have been extensively characterized over the last years. Alterations of the short arm of chromosome 2 (2p) have been of particular interest, since amplification of the MYCN oncogene on 2p24 is associated with an adverse outcome in NB patients. Here, we report on the characterization of a novel genomic rearrangement involving genetic material from 2p13 and 2p24 in NB cell lines that was discovered based on a serial analysis of gene expression (SAGE) profile of the MYCN-amplified NB cell line IMR-5. By analysis of a highly expressed SAGE tag not matching a Unigene cluster we identified four alternatively spliced corresponding transcripts, each of which consisted of the first 14 exons of the anthrax toxin receptor 1 gene (2p13.1) and varying combinations of exons of an unidentified gene located 1.3 Mb telomeric of MYCN (2p24.3) that was termed novel neuroblastoma gene 1. By Southern Blotting, Fluorescent In Situ Hybridization and Long Distance Inverse-PCR we disclosed that these transcripts result from a genomic alteration including material from distinct regions of chromosome 2p and four genomic breakpoints that are joined by short sequences of unknown origin. Furthermore, we show that this rearrangement lies within the homogeneous staining regions (HSR) in IMR-32 cells and is prevalent in both IMR-32 cells and their sub-clones IMR-5 and IMR-5/75, but not in a panel of 70 primary NB tumors. Our work is the first study discovering a fusion transcript based on a SAGE profile and for the first time precisely describes the DNA sequence of amplified breakpoint regions in NB.