- Open Access
An 18 gene expression-based score classifier predicts the clinical outcome in stage 4 neuroblastoma
Journal of Translational Medicine volume 14, Article number: 142 (2016)
The prognosis of children with metastatic stage 4 neuroblastoma (NB) has remained poor in the past decade.
Patients and methods
Using microarray analyses of 342 primary tumors, we here developed and validated an easy to use gene expression-based risk score including 18 genes, which can robustly predict the outcome of stage 4 patients.
This classifier was a significant predictor of overall survival in two independent validation cohorts [cohort 1 (n = 214): P = 6.3 × 10−5; cohort 2 (n = 27): P = 3.1 × 10−2]. The prognostic value of the risk score was validated by multivariate analysis including the established markers age and MYCN status (P = 0.027). In the pooled validation cohorts (n = 241), integration of the risk score with the age and/or MYCN status identified subgroups with significantly differing overall survival (ranging from 35 to 100 %).
Together, the 18-gene risk score classifier can identify patients with stage 4 NB with favorable outcome and may therefore improve risk assessment and treatment stratification of NB patients with disseminated disease.
Neuroblastoma (NB) is the most frequent solid tumor of early childhood with a remarkable variation in clinical presentation ranging from favorable localized tumors that can spontaneously regress to metastatic disease with unfavorable outcome . Within the cohort of patients with disseminated disease the International Neuroblastoma Staging System (INSS) separates unfavorable stage 4 NB, which comprises about 45–50 % of the cases and is defined as a primary tumor with dissemination to distant lymph nodes, bone, bone marrow, liver, skin, or other organs , from favorable stage 4s (special) disease.
A recent review study on 11,037 children with NB from Australia, Europe, Japan, North America has shown that, during the period between 1974 and 2002 the event-free survival of stage 1, 2, 3 and 4 s patients has consistently increased while it hardly changed for stage 4 patients . Another paper reviewed the clinical and survival data of 2216 children with NB enrolled in the Italian Neuroblastoma Registry over a 27-year period (1979–2005). From 1992 to 2005, the overall survival (OS) of patients with stage 3 significantly improved from 67.3 to 88.5 %, whereas the OS of stage 4 patients increased only by 3 % (26–29 %) . These findings demonstrate that no substantial progress in survival has been made for stage 4 patients.
To date, age and MYCN status remain the most important markers of outcome in patients with stage 4 NB. Patients ≥18 months of age with stage 4 NB and those with MYCN-amplified stage 4 disease are defined as “high-risk” . They are usually sensitive to dose-intensive chemotherapy: a majority of patients achieve remission after induction chemotherapy, surgery and radiotherapy, but most patients relapse even with consolidation therapy. Despite intensive multimodal treatment, these high-risk NB patients therefore have an OS of less than 40 %  and discrimination of ultra-poor outcome patients from those with a more favorable prognosis remains poor with current classification systems.
Gene expression profiling by means of microarrays [5, 6] has been shown to be useful in classifying tumors and predicting patient outcome in various types of cancer [7–9]. As such, numerous prognostic gene signatures have been developed to classify NB patients [10, 11]. In spite of the robustness of the published signatures in predicting NB outcome, so far none has been introduced into clinical risk stratification systems. This is probably due to various reasons such as: (i) most of the gene classifiers are built on heterogeneous patient cohorts without differentiating among INSS stage and other clinical and genetic markers, (ii) gene expression profiles can vary according to the microarray platform, (iii) analytic strategy used. Only two prognostic studies have been focused on the well-defined molecular and clinical subgroup of patients with metastatic NB lacking MYCN amplification [12, 13].
To overcome these limitations, we developed a robust and reproducible 18-gene expression based risk-scoring system able to predict OS of children with stage 4 NB by using a different and innovative strategy based on microarrays independent validations.
Patients and methods
The analytic strategy is shown in Fig. 1. Our study has selected 520 unique probes for 426 genes as putative clinical markers based on the following sources:
203 markers by Medline search through the PubMed database (1990–2014) by using the keywords “Neuronal marker” and “Neuroendocrine marker”. 115 articles were selected on the basis of the best available evidence for the specific question;
112 genes Oberthuer et al. 
55 genes Asgharzadeh et al. 
59 genes Vermeulen et al. 
32 genes Fardin et al. 
14 genes Asgharzadeh et al. ;
101 genes/probes obtained by analysis of publicly available microarray of gene expression data (GSE9169) on different NB cell lines treated with retinoic acid (RA). This latter data were used to select genes involved in morphologic differentiation to ganglioneuromatous histopathology which are recognized as a positive prognostic sign in NB.
Identification of the optimal gene set to predict overall survival
The optimal outcome predictor of stage 4 NB was built by using normalized gene expression array data of two independent sets of NB patients (n = 142):
“Seeger dataset” including 102 stage 4 samples downloaded from the website Oncogenomics (http://home.ccr.cancer.gov/oncology/oncogenomics/). Affymetrix HG-U133A and HG-U133B array (GSE16254);
“Versteeg dataset” including 40 stage 4 samples downloaded from the website R2: microarray analysis and visualization platform (http://hgserver1.amc.nl/cgi-bin/r2/main.cgi). Affymetrix HG-U133 Plus 2.0 array (GSE16476).
Development of the gene expression–based prognostic model is described in Additional file 1.
Development of the 18-gene stage 4 NB (Stage4NB) risk score
We adopted a previously developed strategy using the Cox regression coefficient of each gene among 18-gene set from the Seeger cohort [17, 18]. The risk score for each patient was derived by multiplying the expression level of a gene by its corresponding coefficient (risk score = sum of Cox coefficient of Gene Gi X expression value of Gene Gi). Patients of Seeger dataset were dichotomized into both a high-risk group and a low-risk group, using the 50th percentile (median) cutoff of the risk score as the threshold value.
Validation of the 18-gene stage 4 NB (Stage4NB) risk score score in independent test sets
Both the coefficient and the threshold value derived from the Seeger cohort were directly applied to the gene expression data from the exploration data set (Seeger cohort) and independent test sets of microarray experiments termed the “validation cohort 1” and “validation cohort 2” comprising gene expression microarray data of 214 (GSE45547) and 27 (GSE79910) stage 4 NB patients collected at University of Cologne Children’s Hospital in Germany (n = 20) and Gaslini Children Hospital in Italy (n = 7). Table 1 summarizes the features of patient cohorts. The latter set of 27 tumors has been newly analyzed for this study.
RNA preparation was performed essentially as described previously . Collection of patient data and samples was by the clinicians responsible for patient care, with written informed consent obtained from all children’s parents or legal guardians, and with approval by local university ethical committee (Ethical Committee of the University Federico II, C. Romano, Napoli, Italy). Subsequently, single-color gene-expression profiles were generated using customized 4 × 44 K oligonucleotide microarrays produced by Agilent Technologies (Palo Alto, CA, USA) as described previously . More details are in Additional file 1.
Kaplan–Meier estimates for OS were calculated and compared by log-rank test. Only death from disease was considered as an event. Cox regression models were applied using a stepwise selection procedure recommended by Collett  to analyze the prognostic value of potentially prognostic factors. Support Vector Machines (SVM)-based area under receiver operating characteristic (ROC) curve method as implemented in Gene Expression Model Selector (GEMS) software utilizing ten-fold cross-validation and linear polynomial kernel for SVM  was used to evaluate the performance of the published gene signatures and the 18-gene optimal outcome predictor. The gene network and gene ontology (GO) analysis was performed by the website GeneMANIA .
Selection of 18 genes for predicting overall survival in stage 4 NB
The analytic approach based on a re-analysis of public data and the application of Cox regression and ROC curve method (Fig. 1; Additional file 1: Fig. S1 and S2) allowed us to identify a signature as predictor of overall survival (OS) composed of 20 probes for 18 unique genes for patients with stage 4 NB (Tables 1, 2). The predictor score that maximized the area under the curve (AUC) in the Seeger dataset (Additional file 1: Fig. S3) contained 15 genes (AUC = 0.94) whereas that in Versteeg dataset (Additional file 1: Fig. S3) contained 4 genes (AUC = 0.96). One gene (FOXP1) was shared between the two analyses. The distribution of the 18 genes according to the original source is reported in Additional file 1: Table S4. To evaluate the prognostic ability of each signature included in this study, we used SVM-based AUC analysis. The 18-gene signature predicted the OS of better than other published gene signatures [12–16] in both the Seeger (AUC = 94.63) and Versteeg (AUC = 88.33) datasets (Additional file 1: Table S5).
Development and validation of the 18-gene Stage4NB risk score in independent test sets
We generated an outcome predictor of a minimal size and maximum accuracy using 18 unique genes and a model that is based on the relative contributions of each gene. The risk score for each patient was calculated using the regression coefficient of each gene in the 18-gene signature (Table 2). Patients in the Seeger training set were dichotomized according to their 18-gene Stage4NB risk score, and OS was significantly worse in the patient group with a high-risk score (P = 1.1 × 10−12; Fig. 2a). The 5-year survival in low- and high-risk groups in the Seeger cohort was 88 ± 5 and 20 ± 6 %, respectively. Then, gene expression data from validation cohort 1 were analyzed using the 18-gene Stage4NB risk score. With direct application of the Cox regression coefficient from the Seeger training set and the 50th percentile cutoff threshold, OS in the two patient groups differed significantly in validation cohort 1 (P = 6.3 × 10−5; Fig. 2b). Finally, the 18-gene Stage4NB risk score was further validated using an independent cohort of 27 tumors profiled by the same array platform (4 × 44 K oligonucleotide microarrays, Agilent). Again, the Stage4NB risk score separated patient subgroups with a more favorable and an unfavorable outcome (P = 3.1 × 10−2; Fig. 2c). The 5-year OS of patients classified to be favorable or unfavorable was 82 ± 7 and 43 ± 4 % for validation cohort 1, respectively, and 73 ± 16 and 31 ± 12 % for validation cohort 2, respectively. In the pooled validation cohorts, the 18-gene Stage4NB risk score showed high capability in identifying patients at different risk levels (P = 1.1 × 10−5; 5 years OS: 80 ± 6, 42 ± 4 %; Fig. 2d). The same cutoff threshold predicted very well the patients’ event-free survival in validation cohort 1 whereas its prediction ability was less marked in validation cohort 2 (Additional file 1: Fig. S4). One possible explanation could be that the classifier is more suitable for the identification of patients at different risk levels of death as it has been built by using only OS data which is in line with the purpose of this research work.
The 18-gene Stage4NB risk score is an independent prognostic factor
To test whether the risk score is an independent prognostic factor in the pooled validation cohorts, we performed multivariate Cox regression analysis of available prognostic factors, including the risk score, age at diagnosis and MYCN status. The risk score turned out to be a significant prognostic marker in addition to age and the MYCN status (Table 3). We therefore combined these three prognostic markers in a novel risk stratification system for stage 4 NB. Kaplan–Meier analysis of the pooled validation cohort according to this system showed that age ≥18 months and/or MYCN amplification identify subsets with significantly differing survival among favorably and unfavorably classified patients (Fig. 3a). The combination of these three independent prognostic variables thus provided an accurate tool to identify subgroups of patients with a substantially distinct risk to die from disease, ranging from a 5-year survival of 35 ± 4 % to 100 % (P = 2.23 × 10−11), respectively (Fig. 3b). The 18-gene Stage4NB risk score also showed a high capability in identifying patients at different risk levels in patient groups with age <18 months (P = 0.004), age ≥18 months (P = 0.04), age <18 months and MYCN non-amplified, and age ≥18 months (P = 0.089) and MYCN non-amplified (P = 0.087) (Additional file 1: Fig. S5). By contrast, our risk score was less accurate in predicting outcome of children with MYCN-amplified tumors (P = 0.26). In this subgroup, the vast majority of patients (79/81) was classified as high-risk (Additional file 1: Fig. S5). This observation might be due to the fact that the 18-gene risk score consisted mainly of genes involved in neuronal differentiation, which are generally repressed in MYCN over-expressing cells  (Additional file 1: Table S6; Fig. S6).
Gene network and GO analysis
All 18 genes were down-regulated in the subgroup of patients with a high risk score as compared to those with a low risk score in all analyzed datasets (Additional file 1: Fig. S6). To evaluate the potential functional relevance of the signature genes, we performed gene network and GO analyses. The results showed that “activation of protein kinase A (PKA) activity” was the most enriched biological term (Additional file 1: Fig. S7; Table S6) which is highly involved in neuronal differentiation .
By applying a multistep exploration and validation strategy, we identified and validated a risk score-based classifier using the expression patterns of 18 genes that is able to identify two subsets of stage 4 NB patients with different OS. Children with a high risk score uniformly had a poor outcome, with 20–43 % OS at 5 years after diagnosis, whereas those with a low risk score had 73–88 % OS. Moreover, several lines of evidence strongly support that the 18-gene Stage4NB risk score is an independent and significant predictor of prognosis. First, the risk score was a significant predictive factor for OS in multivariate analysis including MYCN status and age at diagnosis. Second, the combination of the three independent prognostic variables Stage4NB risk score, age and MYCN provided an accurate tool to identify a specific subgroup of patients with favorable outcome (5-year rate of survival of 100 %). Taken together, these results strongly support the notion that the 18-gene risk score identifies groups of patients at different risk levels of death in stage 4 patients and that may represent an accurate tool to improve risk estimation of patients who are currently believed to be at high risk to die from disease.
Recent research has elucidated the biology of NB allowing more accurate stratification, which has permitted to develop appropriate treatments for children with localized tumors reducing cytotoxic therapy and increasing the survival rate . Major challenges still remain for children with metastatic (stage 4) NB older than 18 months or those whose tumors are MYCN-amplified, with 5-year survival rates of only 30–40 % . Our findings suggest that the 18-gene Stage4NB risk score may be capable to improve the current risk stratification system of high-risk patients, as it is able to identify a subset of patients with favorable outcome who may require less intensive therapies. The reduction of treatment intensities for patients with a more favorable outcome may substantially decrease the risk of serious side effects such as sepsis, primary hypothyroidism, growth hormone deficiency, deafness and cardiovascular problems. In our cohort, 22 patients classified as low-risk by our risk score had been treated according to a high-risk protocol. The favorable outcome of these patients indicates that they might have had a similar outcome with less intensive treatment. Thus, we propose that reduction of cytotoxic intensities in such patients should be evaluated prospectively.
Our model included genes that may be interesting for further research based on either their chromosomal location, their known function, or their possible role as drug targets (ADCY1, AKR1C1 and SCNA3). Not surprisingly, this set of 18 predictive genes contains numerous genes that have been reported to have a role in the neuronal differentiation which if arrested contributes to early event in NB pathogenesis as also demonstrated by our recent work on genetic susceptibility to NB . For instance, ADCY1 GNAI1 and PRKACB genes are associated with the cAMP-mediated signaling which plays a crucial role in initiating differentiation in transformed and embryonic cells of neuronal and glial origin . cAMP-stimulating agents also induce differentiation in human and mouse NB cells . ARHGEF10L gene is a member of the Rho family of guanine nucleotide exchange factors (GEF) that activate Rho GTPases. Interestingly, frequent mutations of RAC-RHO pathway genes regulating neuritogenesis have been found in NBs stage 3 and 4 . Further genes reported to have a role in neuronal differentiation are HOXC6 , SOX4 , FOXP1 , GFRA3 , and PTPRH . Our recent study shows the biological role of FOXP1 in contributing to NB progression and unfavorable patient outcome . This is in line with the evidence that high risk neuroblastomas are characterized by low expression of genes involved in neuronal differentiation [15, 35–37]. Importantly, the gene network and GO analysis showed that “PKA activity”, which includes ADCY1 and PRKACB genes, was the most enriched biological term. ADCY1 encodes a form of adenylate cyclase whereas PRKACB encodes a catalytic subunit of PKA . Both genes show lower expression values in patients classified to be at high risk. Recently, adenylyl cyclases have emerged as potential drug target in diverse diseases  whereas PKA signaling pathway is known to antagonize Hedgehog signaling . Interestingly, the activation of PKA pathway by forskolin (ADCY1 activator) has been associated with a reduction of cell proliferation and an induction of apoptosis by inhibition of Hedgehog signal in NB cell lines . Moreover, a recent study demonstrated that the neuropeptide pituitary adenylyl cyclase activating polypeptide (PACAP), another ADCY1 activator, inhibits proliferation of primary medulloblastoma derived tumorsphere cultures by PKA activation and inhibition of Hedgehog signal . Together, these data support the idea that regulation of PKA signaling by ADCY1 activation might be an additional therapeutic strategy for stage 4 NB. The correlation of high AKR1C1 expression with cancer is supported by two recently studies [42, 43], possibly due to the ability of AKR1C1 to act as tumor suppressor gene. Particularly in stromal fibroblasts and carcinoma cells, high AKR1C1 expression correlates with favorable tumor characteristics and longer survival in primary breast cancer patients . Recent advances in the molecular biology of esophageal squamous cell carcinoma (ESCC) have shown that in ESCC patients, high AKR1C1 expression increase the sensitive of ESCC cells to ethyl-3,4-dihydroxybenzoate (EDHB) providing potential guidance for the chemoprevention of ESCC . AKR1C1 shows lower expression values in high risk NB patients and these data supported that therapeutic modulation of the AKR1C1 expression, could be an attractive therapeutic possibility for patients classified to be at high risk.
In conclusion, we have established and validated a robust prognostic scoring system in the largest stage 4 NB population to date (n = 342). Our study provides further evidence that gene expression-based classification works well in NB. We demonstrated an excellent performance of our classifier on independent data sets, involving stage 4 patients from different countries and using two validation cohorts. The gene-based risk score shows high performance in risk estimation of stage 4 NB patients alone and when integrated with the currently used variables age and MYCN status. Together, our findings encourage large prospective studies on the clinical value of the scoring system, which may ultimately improve risk assessment and treatment stratification of children with metastatic NB.
In the last decade, survival rate of children with advanced stage 4 neuroblastoma, the most common solid extracranial tumor of infancy, has improved little whereas remarkable progresses have been recorded for children with localized or 4s disease. Gene expression–based classification has been demonstrated to precisely predict neuroblastoma outcome; however, no such classifier is used in clinical practice to date. Here we present a risk estimation for children affected by advanced stage 4 neuroblastoma which integrates both a highly accurate gene expression–based classifier and established prognostic markers. According to this system, we identified novel subgroups of patients with favourable prognosis among high-risk patients. Our study may help clinicians in choosing a more appropriate therapy to reduce side effects for those children with low-risk profile or a more intensive treatments for those children with high-risk profile.
international neuroblastoma staging system
Maris JM. Recent advances in neuroblastoma. N Engl J Med. 2010;362:2202–11.
Brodeur GM, Pritchard J, Berthold F, Carlsen NL, Castel V, Castelberry RP, De Bernardi B, Evans AE, Favrot M, Hedborg F. Revisions of the international criteria for neuroblastoma diagnosis, staging, and response to treatment. J Clin Oncol. 1993;11:1466–77.
Moroz V, Machin D, Faldum A, Hero B, Iehara T, Mosseri V, Ladenstein R, De Bernardi B, Rubie H, Berthold F, et al. Changes over three decades in outcome and the prognostic influence of age-at-diagnosis in young patients with neuroblastoma: a report from the International Neuroblastoma Risk Group Project. Eur J Cancer. 2011;47:561–71.
Haupt R, Garaventa A, Gambini C, Parodi S, Cangemi G, Casale F, Viscardi E, Bianchi M, Prete A, Jenkner A, et al. Improved survival of children with neuroblastoma between 1979 and 2005: a report of the Italian Neuroblastoma Registry. J Clin Oncol. 2010;28:2331–8.
Hoheisel JD. Microarray technology: beyond transcript profiling and genotype analysis. Nat Rev Genet. 2006;7:200–10.
Schena M, Shalon D, Davis RW, Brown PO. Quantitative monitoring of gene expression patterns with a complementary DNA microarray. Science. 1995;270:467–70.
Ramaswamy S, Ross KN, Lander ES, Golub TR. A molecular signature of metastasis in primary solid tumors. Nat Genet. 2003;33:49–54.
Golub TR, Slonim DK, Tamayo P, Huard C, Gaasenbeek M, Mesirov JP, Coller H, Loh ML, Downing JR, Caligiuri MA, et al. Molecular classification of cancer: class discovery and class prediction by gene expression monitoring. Science. 1999;286:531–7.
Ludwig JA, Weinstein JN. Biomarkers in cancer staging, prognosis and treatment selection. Nat Rev Cancer. 2005;5:845–56.
Capasso M. Neuroblastoma and whole genome searches. In: ÖnerÖzdemir, editors. Contemporary pediatrics. InTech; 2012. pp. 297–314.
Capasso M, Diskin SJ. Genetics and genomics of neuroblastoma. Cancer Treat Res. 2010;155:65–84.
Asgharzadeh S, Pique-Regi R, Sposto R, Wang H, Yang Y, Shimada H, Matthay K, Buckley J, Ortega A, Seeger RC. Prognostic significance of gene expression profiles of metastatic neuroblastomas lacking MYCN gene amplification. J Natl Cancer Inst. 2006;98:1193–203.
Asgharzadeh S, Salo JA, Ji L, Oberthuer A, Fischer M, Berthold F, Hadjidaniel M, Liu CW, Metelitsa LS, Pique-Regi R, et al. Clinical significance of tumor-associated inflammatory cells in metastatic neuroblastoma. J Clin Oncol. 2012;30:3525–32.
Oberthuer A, Berthold F, Warnat P, Hero B, Kahlert Y, Spitz R, Ernestus K, Konig R, Haas S, Eils R, et al. Customized oligonucleotide microarray gene expression-based classification of neuroblastoma patients outperforms current clinical risk stratification. J Clin Oncol. 2006;24:5070–8.
Vermeulen J, De Preter K, Naranjo A, Vercruysse L, Van Roy N, Hellemans J, Swerts K, Bravo S, Scaruffi P, Tonini GP, et al. Predicting outcomes for children with neuroblastoma using a multigene-expression signature: a retrospective SIOPEN/COG/GPOH study. Lancet Oncol. 2009;10:663–71.
Fardin P, Barla A, Mosci S, Rosasco L, Verri A, Versteeg R, Caron HN, Molenaar JJ, Ora I, Eva A, et al. A biology-driven approach identifies the hypoxia gene signature as a predictor of the outcome of neuroblastoma patients. Mol Cancer. 2010;9:185.
Chen HY, Yu SL, Chen CH, Chang GC, Chen CY, Yuan A, Cheng CL, Wang CH, Terng HJ, Kao SF, et al. A five-gene signature and clinical outcome in non-small-cell lung cancer. N Engl J Med. 2007;356:11–20.
Kim SM, Leem SH, Chu IS, Park YY, Kim SC, Kim SB, Park ES, Lim JY, Heo J, Kim YJ, et al. Sixty-five gene-based risk score classifier predicts overall survival in hepatocellular carcinoma. Hepatology. 2012;55:1443–52.
Oberthuer A, Hero B, Berthold F, Juraeva D, Faldum A, Kahlert Y, Asgharzadeh S, Seeger R, Scaruffi P, Tonini GP, et al. Prognostic impact of gene expression-based classification for neuroblastoma. J Clin Oncol. 2010;28:3506–15.
Oberthuer A, Juraeva D, Li L, Kahlert Y, Westermann F, Eils R, Berthold F, Shi L, Wolfinger RD, Fischer M, Brors B. Comparison of performance of one-color and two-color gene-expression analyses in predicting clinical endpoints of neuroblastoma patients. Pharmacogenomics J. 2010;10:258–66.
Collett D. Strategy for model selection. In: Collett D, editor. Modelling survival data in medical research. London: Chapman & Hall; 1994. p. 78–83.
Statnikov A, Tsamardinos I, Dosbayev Y, Aliferis CF. GEMS: a system for automated cancer diagnosis and biomarker discovery from microarray gene expression data. Int J Med Inform. 2005;74:491–503.
Warde-Farley D, Donaldson SL, Comes O, Zuberi K, Badrawi R, Chao P, Franz M, Grouios C, Kazi F, Lopes CT, et al. The GeneMANIA prediction server: biological network integration for gene prioritization and predicting gene function. Nucleic Acids Res. 2010;38:W214–20.
Valentijn LJ, Koster J, Haneveld F, Aissa RA, van Sluis P, Broekmans ME, Molenaar JJ, van Nes J, Versteeg R. Functional MYCN signature predicts outcome of neuroblastoma irrespective of MYCN amplification. Proc Natl Acad Sci USA. 2012;109:19190–5.
Adzhubei IA, Schmidt S, Peshkin L, Ramensky VE, Gerasimova A, Bork P, Kondrashov AS, Sunyaev SR. A method and server for predicting damaging missense mutations. Nat Methods. 2010;7:248–9.
Capasso M, Diskin S, Cimmino F, Acierno G, Totaro F, Petrosino G, Pezone L, Diamond M, McDaniel L, Hakonarson H, et al. Common genetic variants in NEFL influence gene expression and neuroblastoma risk. Cancer Res. 2014;74:6913–24.
Prasad KN, Cole WC, Yan XD, Nahreini P, Kumar B, Hanson A, Prasad JE. Defects in cAMP-pathway may initiate carcinogenesis in dividing nerve cells: a review. Apoptosis. 2003;8:579–86.
Molenaar JJ, Koster J, Zwijnenburg DA, van Sluis P, Valentijn LJ, van der Ploeg I, Hamdi M, van Nes J, Westerman BA, van Arkel J, et al. Sequencing of neuroblastoma identifies chromothripsis and defects in neuritogenesis genes. Nature. 2012;483:589–93.
Zhang X, Hamada J, Nishimoto A, Takahashi Y, Murai T, Tada M, Moriuchi T. HOXC6 and HOXC11 increase transcription of S100beta gene in BrdU-induced in vitro differentiation of GOTO neuroblastoma cells into Schwannian cells. J Cell Mol Med. 2007;11:299–306.
Potzner MR, Tsarovina K, Binder E, Penzo-Mendez A, Lefebvre V, Rohrer H, Wegner M, Sock E. Sequential requirement of Sox4 and Sox11 during development of the sympathetic nervous system. Development. 2010;137:775–84.
Dasen JS, De Camilli A, Wang B, Tucker PW, Jessell TM. Hox repertoires for motor neuron diversity and connectivity gated by a single accessory factor, FoxP1. Cell. 2008;134:304–16.
Masure S, Geerts H, Cik M, Hoefnagel E, Den Van, Kieboom G, Tuytelaars A, Harris S, Lesage AS, Leysen JE, Der Van, Helm L, et al. Enovin, a member of the glial cell-line-derived neurotrophic factor (GDNF) family with growth promoting activity on neuronal cells. Existence and tissue-specific expression of different splice variants. Eur J Biochem. 1999;266:892–902.
Johnson KG, Van Vactor D. Receptor protein tyrosine phosphatases in nervous system development. Physiol Rev. 2003;83:1–24.
Ackermann S, Kocak H, Hero B, Ehemann V, Kahlert Y, Oberthuer A, Roels F, Theissen J, Odenthal M, Berthold F, Fischer M. FOXP1 inhibits cell growth and attenuates tumorigenicity of neuroblastoma. BMC Cancer. 2014;14:840.
Benard J, Raguenez G, Kauffmann A, Valent A, Ripoche H, Joulin V, Job B, Danglot G, Cantais S, Robert T, et al. MYCN-non-amplified metastatic neuroblastoma with good prognosis and spontaneous regression: a molecular portrait of stage 4S. Mol Oncol. 2008;2:261–71.
Fischer M, Oberthuer A, Brors B, Kahlert Y, Skowron M, Voth H, Warnat P, Ernestus K, Hero B, Berthold F. Differential expression of neuronal genes defines subtypes of disseminated neuroblastoma with favorable and unfavorable outcome. Clin Cancer Res. 2006;12:5118–28.
Fredlund E, Ringner M, Maris JM, Pahlman S. High Myc pathway activity and low stage of neuronal differentiation associate with poor outcome in neuroblastoma. Proc Natl Acad Sci USA. 2008;105:14094–9.
Pierre S, Eschenhagen T, Geisslinger G, Scholich K. Capturing adenylyl cyclases as potential drug targets. Nat Rev Drug Discov. 2009;8:321–35.
Hammerschmidt M, Bitgood MJ, McMahon AP. Protein kinase A is a common negative regulator of Hedgehog signaling in the vertebrate embryo. Genes Dev. 1996;10:647–58.
Yamanaka H, Oue T, Uehara S, Fukuzawa M. Forskolin, a Hedgehog signal inhibitor, inhibits cell proliferation and induces apoptosis in pediatric tumor cell lines. Mol Med Rep. 2010;3:133–9.
Cohen JR, Resnick DZ, Niewiadomski P, Dong H, Liau LM, Waschek JA. Pituitary adenylyl cyclase activating polypeptide inhibits gli1 gene expression and proliferation in primary medulloblastoma derived tumorsphere cultures. BMC Cancer. 2010;10:676.
Li W, Hou G, Zhou D, Lou X, Xu Y, Liu S, Zhao X. The roles of AKR1C1 and AKR1C2 in ethyl-3,4-dihydroxybenzoateinduced esophageal squamous cell carcinoma cell death. Oncotarget. 2016.
Wenners A, Hartmann F, Jochens A, Roemer AM, Alkatout I, Klapper W, van Mackelenbergh M, Mundhenke C, Jonat W, Bauer M. Stromal markers AKR1C1 and AKR1C2 are prognostic factors in primary human breast cancer. Int J Clin Oncol. 2015.
Study concept and design was by MC, MF, and DF. GP, VAL, and PP were responsible for data acquisition. Quality control of data was carried out by MC, LL, SV and FC. Data analysis and interpretation was performed by DF, AI, AO, GPT. DF, VAL, GP, and MC performed the statistical analysis. The manuscript was prepared, edited, and reviewed by DF, MC, and MF. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
This study was supported by grants from Associazione Italiana per la Ricerca sul Cancro (Grant no. 10537 to MC), MIUR- FIRB Ricerca in Futuro (Grant no. RBFR08DWQ3 to MC), Fondazione Italiana per la Lotta al Neuroblastoma (to MC) and Associazione Oncologia Pediatrica e Neuroblastoma (to MC), the German Cancer Aid (Grant no. 110122 to MF), the German Ministry of Science and Education (BMBF) as part of the e:Med initiative (Grant nos. 01ZX1303A and 01ZX1307D to MF) and the Fördergesellschaft Kinderkrebs-Neuroblastom-Forschung e.V. (to MF).