CASP8 SNP D302H (rs1045485) Is Associated with Worse Survival in MYCN-Amplified Neuroblastoma Patients

Abstract

Background

Neuroblastoma is a pediatric cancer that exhibits a wide clinical spectrum ranging from spontaneous regression in low-risk patients to fatal disease in high-risk patients. The identification of single nucleotide polymorphisms (SNPs) may help explain the heterogeneity of neuroblastoma and assist in identifying patients at higher risk for poor survival. SNPs in the TP53 pathway are of special importance, as several studies have reported associations between TP53 pathway SNPs and cancer. Of note, less than 2% of neuroblastoma tumors have a TP53 mutation at diagnosis.

Patients and Methods

We selected 21 of the most frequently studied SNPs in the TP53 pathway and evaluated their association with outcome in 500 neuroblastoma patients using TaqMan allelic discrimination assays.

Results and Conclusion

We investigated the impact of 21 SNPs on overall survival, event-free survival, age at diagnosis, MYCN status, and stage of the disease in 500 neuroblastoma patients. A missense SNP in exon 10 of the CASP8 gene SNP D302H was associated with worse overall and event-free survival in patients with MYCN-amplified neuroblastoma tumors.

About the Authors:

Ali Rihani

Affiliation: Center for Medical Genetics, Ghent University, Ghent, Belgium

Bram De Wilde

Affiliation: Center for Medical Genetics, Ghent University, Ghent, Belgium

Fjoralba Zeka

Affiliation: Center for Medical Genetics, Ghent University, Ghent, Belgium

Geneviève Laureys

Affiliation: Department of Pediatric Hematology, Oncology and Stem Cell Transplantation, Ghent University Hospital, Ghent, Belgium

Nadine Francotte

Affiliation: Département de pédiatrie, hémato-oncologie, SUHOPL- CHC (Service Universitaire d′HématoOncologie Pédiatrique Centre Hospitalier Chrétien) Espérance, St Nicolas Belgium

Gian Paolo Tonini

Affiliation: Neuroblastoma Laboratory, Onco/Hematology Laboratory, Department of Women's and Children's Health, University of Padua, Pediatric Research Institute, Fondazione Città della Speranza, Padua, Italy

Simona Coco

Affiliation: Lung Cancer Unit, IRCCS (Istituto di Ricovero e Cura a Carattere Scientifico), Azienda Ospedaliera Universitaria San Martino – IST (Istituto Nazionale per la Ricerca sul Cancro), Genoa, Italy

Rogier Versteeg

Affiliation: Departement of Human Genetics, Academic Medical Center, University of Amsterdam, Amsterdam, the Netherlands

Rosa Noguera

Affiliation: Department of Pathology, Medical School, University of Valencia, Valencia, Spain

Johannes H. Schulte

Affiliations German Cancer Consortium (DKTK), Germany, Translational Neuro-Oncology, West German Cancer Center, University Hospital Essen, University Duisburg-Essen, Essen, Germany, German Cancer Research Center (DKFZ), Heidelberg, Germany, Department of Pediatric Oncology and Haematology, University Children's Hospital Essen, Essen, Germany, Centre for Medical Biotechnology, University Duisburg-Essen, Essen, Germany

Angelika Eggert

Affiliation: Department of Pediatric Oncology and Haematology, University Children's Hospital Essen, Essen, Germany

Raymond L. Stallings

Affiliation: Department of Molecular and Cellular Therapeutics, Royal College of Surgeons in Ireland, and National Children′s Research Centre, Dublin, Ireland

Frank Speleman

Affiliation: Center for Medical Genetics, Ghent University, Ghent, Belgium

Jo Vandesompele

Affiliation: Center for Medical Genetics, Ghent University, Ghent, Belgium

Tom Van Maerken

* E-mail: [email protected]

Affiliation: Center for Medical Genetics, Ghent University, Ghent, Belgium

Abstract

Background

Neuroblastoma is a pediatric cancer that exhibits a wide clinical spectrum ranging from spontaneous regression in low-risk patients to fatal disease in high-risk patients. The identification of single nucleotide polymorphisms (SNPs) may help explain the heterogeneity of neuroblastoma and assist in identifying patients at higher risk for poor survival. SNPs in the TP53 pathway are of special importance, as several studies have reported associations between TP53 pathway SNPs and cancer. Of note, less than 2% of neuroblastoma tumors have a TP53 mutation at diagnosis.

Patients and Methods

We selected 21 of the most frequently studied SNPs in the TP53 pathway and evaluated their association with outcome in 500 neuroblastoma patients using TaqMan allelic discrimination assays.

Results and Conclusion

We investigated the impact of 21 SNPs on overall survival, event-free survival, age at diagnosis, MYCN status, and stage of the disease in 500 neuroblastoma patients. A missense SNP in exon 10 of the CASP8 gene SNP D302H was associated with worse overall and event-free survival in patients with MYCN-amplified neuroblastoma tumors.

Figures

Table 3

Figure 2

Table 4

Table 1

Figure 1

Table 2

Table 3

Figure 2

Table 4

Table 1

Figure 1

Table 2

Citation: Rihani A, De Wilde B, Zeka F, Laureys G, Francotte N, Tonini GP, et al. (2014) CASP8 SNP D302H (rs1045485) Is Associated with Worse Survival in MYCN-Amplified Neuroblastoma Patients. PLoS ONE 9(12): e114696. https://doi.org/10.1371/journal.pone.0114696

Editor: Roberto Amendola, ENEA, Italy

Received: June 26, 2014; Accepted: November 13, 2014; Published: December 11, 2014

Copyright: © 2014 Rihani et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: The authors confirm that all data underlying the findings are fully available without restriction. All relevant data are within the paper and its Supporting Information files.

Funding: AR is supported by a PhD fellowship from the Ghent University research fund (BOF; 01D02210). TVM is a postdoctoral fellow of the FWO (12B9912N). RN is supported by FIS (contract PI10/15) and RTICC (contracts RD06/0020/0102; RD12/0036/0020), Instituto Carlos III Madrid & ERDF, Spain. Italian Neuroblastoma Foundation supported the Italian samples collection and preparation. This work was supported by the BOF, the FWO, the Childhood Cancer Fund, and the National Cancer Plan of the Belgian State (Action 29). The study sponsors have no role in the study design, in the collection, analysis and interpretation of data; in the writing of the manuscript; or in the decision to submit the manuscript for publication.

Competing interests: The authors have declared that no competing interests exist.

Background

Neuroblastoma (NB) is a complex and heterogeneous pediatric cancer that can manifest with tumors that regress spontaneously, and also with tumors that metastasize and acquire resistance to therapy, leading to severe illness and death [1]. Currently, the treatment protocols depend on the risk stratification at diagnosis, which in turn is based on the clinico-genetic features of the patient and the tumor. The risk stratification system has been improved over the past decades [2]; however, it remains imperative to further exploit the underlying biology of clinico-genetic features of the disease for the sake of tailoring therapies to defined patient populations. Inter-individual variability in response to therapy can result from single nucleotide polymorphisms (SNPs) in critical genes involved in cell cycle control and induction of apoptosis [3]. Some SNPs have proven their usefulness in serving as prognostic markers in development of targeted therapies, or defining patient populations in terms of potential response to conventional therapies. One of the most critical genes that acts as a guardian of the genome and also mediates response to therapy is TP53 [4]. TP53 is mutated in more than 50% of human cancers [5]; however, less than 2% of NB tumors harbor a TP53 mutation at diagnosis [6]. This suggests that other players in the TP53 pathway such as SNPs may play a role in NB.

Genetic polymorphisms in TP53 pathway components have been targeted for personalizing current therapies and developing new treatment modalities [7]. However, only a limited number of studies have focused on the association between SNPs in TP53 pathway genes and NB. Therefore, we selected 21 SNPs in fifteen TP53 pathway associated genes (Table 1), which have been reported to be involved in susceptibility to several cancer types. SNP selection was primarily based on a study that identified potential SNPs involved in TP53 stress response using publicly available genotypes and drug response data from the NCI60 human tumor cell lines [8]. NCI60 is a well-characterized panel of 60 human cancer cell lines from different cancer entities with publicly available anti-cancer drug screening data [9] and genotype data from 100,000 SNPs in cancer related genes [10], [11]. We expanded this selection to include SNPs in TP53 pathway associated genes that have been reported elsewhere to be involved in cancer susceptibility [4], [12]–[15]. The 21 SNPs are found in genes involved in DNA damage response, apoptotic response, cell cycle regulation, direct regulation of TP53, or the TP53 gene itself (Table 1). We genotyped these SNPs in 500 NB cases and investigated the association between these SNPs and overall and event-free survival, age at diagnosis, MYCN status, and the stage of the disease of NB patients.

[Figure omitted. See PDF.]

Table 1. 21 SNPs in 15 p53 pathway genes.

https://doi.org/10.1371/journal.pone.0114696.t001

Materials and Methods

Patients

The study group comprised 500 NB patients who were evaluated according to the International NB Staging System (INSS) [16]. All patients included in this study provided written informed consent for anonymous use of the samples for research related to neuroblastoma biology. Approval for this study was granted by the medical ethical committee of Ghent University Hospital. Registration number: B670201111331.

Statistical methods

Comparisons of Kaplan-Meier survival curves between different genotype groups were done using the log-rank test. Survival time was defined as the time from diagnosis until death from the disease or until the date of the last follow-up of living patients. Patient samples were divided into two groups based on SNP status (wild-type versus homozygous or heterozygous for the minor allele). Comparisons of age at diagnosis between the different groups were done using the Mann-Whitney U test. Association of the SNPs with MYCN status or stage of the disease was done using Pearson's Chi-square test. Statistical tests were two-sided and results with corrected p-values (q-values) <0.05 were considered statistically significant. All statistical analyses were conducted using SPSS version 20 software. All statistical analyses were followed by Benjamini-Hochberg multiple testing correction in R version 3.0.2 using “stats” statistical package.

SNP selection and genotyping

SNPs in potentially important TP53 pathway genes were selected from previously reported studies [4], [8], [12]–[15]. The protein products of these genes function at different levels of the TP53 pathway. MDM2 SNPs (285, 309, 344, and 354) were genotyped by Sanger sequencing and the remaining 17 SNPs in TP53 pathway genes were genotyped using the TaqMan allelic discrimination assay. NB tumor DNA was first amplified using illustra GenomiPhi V2 DNA amplification kit (GE Healthcare Life Sciences, Buckinghamshire UK) according to the manufacturer's protocol. Genotyping was carried out using 10 ng of amplified genomic DNA per reaction, with 2.5 µL of 2X TaqMan Genotyping Master Mix and 160 nM primers (Life Technologies, Merelbeke, Belgium) in a 5 µL reaction volume. All genotypes were determined by endpoint calling on an ABI 7900 real-time PCR instrument (Life Technologies, Merelbeke, Belgium). PCR cycling conditions were 50°C for 2 min, 95°C for 10 min, followed by 40 cycles of 92°C for 15 s and 60°C for 1 min. Sanger sequencing was carried out by PCR amplification and subsequent sequencing. The primers used for PCR amplification were tagged with universal M13 sequencing tags (underlined): forward primer CACGACGTTGTAAAACGACTGGCTTTGCGGAGGTT and reverse primer CAGGAAACAGCTATGACCTCGGAACGTGTCTGAACTT. The PCR was performed in a 25 µl reaction mix using 20 ng of DNA, 1x KAPA2G buffer, 2.5 mM MgCl₂, 0.2 mM dNTP mix, 0.5 µM forward and reverse primers, and 1 U KAPA2G Robust HotStart DNA polymerase. PCR conditions were as follows: initial denaturation at 94°C for 4 min, denaturation at 94°C for 20 s, primer annealing at 68°C for 15 s, and extension at 72°C for 1 min for 12 cycles. This was followed by another 34 cycles: denaturation at 94°C for 30 s, primer annealing at 55°C for 30 s, and extension at 72°C for 30 s. Finally, an extension step at 72°C for 4 min was performed, followed by cooling at 15°C for 1 min. The size of the amplified PCR product was verified using Caliper Labchip GX (PerkinElmer, Waltham, MA, USA). The sequencing was done using Sanger DNA sequencing service (Genewiz, South Plainfield, NJ, USA). The generated data was analyzed using Seqscape software v.2.7 (Life Technologies, Ghent, Belgium). 18 SNPs were genotyped in 31 NB cell lines (amplified and non-amplified DNA) to ensure the quality of the DNA amplification procedure.

Reverse transcription quantitative polymerase chain reaction (RT-qPCR)

RT-qPCR was performed on duplicate samples using 5 ng of template cDNA, 2.5 µL of 2x Soo Advanced Reaction Mix (Roche), and 0.25 µL of a 5 µM solution of each primer in a 5 µL total reaction volume. The PCR conditions were as follows: 10 min at 95°C, followed by 45 cycles of denaturation (10 s at 95°C), and elongation (45 s at 60°C) using LightCycler480 (Roche). Primers for caspase-8 and the reference genes (Alu-Sq, HMBS, HPRT1, SDHA, and UBC) were designed using an in-house web tool (www.primerXL.org). The sequence of the primers can be found in Table S1 in S1 File. Relative expression of caspase 8 was quantified using Biogazelle's qbase+ qPCR data-analysis software version 2.y [17]

Results

CASP8 SNP D302H and survival in NB patients

We used whole genome amplified DNA from tumors of 500 NB patients to genotype 21 SNPs in fifteen TP53 pathway genes (Table 1) using hydrolysis probe based qPCR genotyping assays for 17 SNPs, and Sanger sequencing for four MDM2 SNPs (285, 309, 344, and 354). The clinical features of the patients are shown in Table 2. The study population consisted mostly of European patients from Caucasian origin. We first validated the approach by genotyping the SNPs on whole genome amplified and non-amplified DNA from 31 NB cell lines, and results showed that all SNPs were genotyped correctly (Table S2 in S1 File). In patients, CASP8 SNP D302H was the only SNP that showed an association with worse overall (OS) (p = 0.0006; multiple testing corrected p-value, q-value = 0.049) and event-free survival (EFS) (p = 0.0002; multiple testing corrected p-value, q-value = 0.042) in NB patients with MYCN amplification (Fig. 1, Table 3). The stratified survival analyses of all the other SNPs are shown in table S3 in S1 File. To confirm that there is no difference in the SNP status between the tumor and the germline of the patients, we collected 41 representative blood samples from the patients in our cohort and genotyped them using the same TaqMan assay for the CASP8 SNP D302H. Our results show that all the 41 patients have exactly the same SNP status in their blood as in their tumor (data not shown). This additional analysis makes it highly unlikely that different SNP status exists in the tumor compared to the blood of the patients. We performed also multivariate Cox proportional hazard analysis and investigated the effect of CASP8 SNP D302H on the survival of neuroblastoma patients taking into account the different clinical covariates (stage of the disease, MYCN status, and the age at diagnosis). Our results show that CASP8 SNP D302H is an independent prognostic factor (Table S4 and S5 in S1 File)

[Figure omitted. See PDF.]

Figure 1. Overall and event-free survival of NB patients by CASP8 SNP D302H.

Comparison of Kaplan-Meier survival curves between different genotypes of CASP8 SNP D302H. Overall survival of MYCN-amplified NB patients (A), event-free survival of MYCN-amplified NB patients (B). Raw P-values were calculated by the log rank test. q-values are adjusted p-values after Benjamini-Hochberg multiple testing correction.

https://doi.org/10.1371/journal.pone.0114696.g001

[Figure omitted. See PDF.]

Table 2. Clinical characteristics of the patients.

https://doi.org/10.1371/journal.pone.0114696.t002

[Figure omitted. See PDF.]

Table 3. Survival of NB patients by CASP8 D302H status.

https://doi.org/10.1371/journal.pone.0114696.t003

Impact of SNPs on age at diagnosis and association with the stage of the disease and MYCN status

We investigated the correlation between 21 SNPs and NB onset in a non-stratified analysis (all subgroups) and also in an analysis stratified by MYCN status or by clinical stage (localized vs. stage 4). CASP8 SNP D302H was not correlated to age at diagnosis (Table S6 in S1 File). The same results were obtained for the remaining SNPs (data not shown). We also investigated whether CASP8 SNP D302H or any of the 21 SNPs is associated with the MYCN status or stage of the disease. CASP8 SNP D302H was not associated with MYCN status or stage of the disease (Table S7 in S1 File). Similar results were obtained for the remaining SNPs (data not shown).

CASP8 SNP D302H and expression of caspase 8

Since decreased expression of caspase 8 is a recurrent event in NB [18], we tested whether CASP8 SNP D302H could affect the expression of caspase 8. We measured caspase 8 mRNA expression in two independent cohorts of NB tumor samples. Cohort 1 included 162 samples treated according to the International Society of Pediatric Oncology Europe Neuroblastoma Group (SIOPEN, https://www.siopen-r-net.org/). Cohort 2 included 161 samples from the Neuroblastoma Research Consortium (NRC), which is a collaboration between several European NB research groups. All samples from cohorts 1 and 2 were derived from patients that had been included in the entire cohort of 500 patients analyzed in Table 1. We found that CASP8 SNP D302H had no effect on caspase 8 mRNA expression levels, regardless of MYCN status or the stage of the disease (Table 4, Fig. 2).

[Figure omitted. See PDF.]

Figure 2. Expression of caspase-8 by CASP8 SNP D302H.

Comparison of caspase-8 expression between the different genotype groups of CASP8 SNP D302H. Cohort 1 (A), Cohort 2 (B). P-values were calculated by the Mann-Whitney U test. q-values are adjusted p-values after Benjamini-Hochberg multiple testing correction.

https://doi.org/10.1371/journal.pone.0114696.g002

[Figure omitted. See PDF.]

Table 4. Distribution of caspase-8 expression by CASP8 D302H.

https://doi.org/10.1371/journal.pone.0114696.t004

Discussion

NB is a heterogeneous disease with a broad clinical spectrum ranging from spontaneous regression and excellent outcome in low-risk patients to widely disseminated disease and poor survival in high-risk patients. The phenomenon of spontaneous NB regression has been attributed to successful delayed apoptosis of tumor cells, whereas progression towards uncontrolled proliferation and tumor formation largely reflects defects in the apoptotic process [19]. Evasion of growth suppressors such as TP53 is a hallmark of tumor formation and development [20], and functional inactivation of the TP53 pathway could be achieved either via mutations in the TP53 gene itself or anywhere along the TP53 regulatory network. Because TP53 is rarely mutated in NB at the time of diagnosis (<2%) [21], other components of the TP53 pathway could be responsible for circumvention of this fundamental antitumor barrier. Numerous reports have suggested that germline alterations in TP53 and TP53 pathway genes are associated with cancer risk and clinical outcomes [4]. Several SNPs that predispose to NB have been identified by genome-wide association studies (GWAS). The SNPs most highly associated with increased NB susceptibility and poor outcome were found in HACE1 (rs4336470) [22], LIN28 (rs17065417) [22], [23], BARD1 (rs6435862, rs3768716) [24], LMO1 (rs110419) [25], DUSP12 (rs1027702), DDX1 (rs2619046), IL31RA (rs10055201), and HSD17B12 (rs11037575) [26]. However, despite the power of GWAS for identifying new predisposing SNPs, these studies still lack the ability to identify all SNPs and there is a high probability of false negatives [27]. Rather than choosing a GWAS approach with complex data analysis, we decided to conduct a focused study to examine particular SNPs in the TP53 pathway that are thought to influence susceptibility to cancer [4], [8], [12]–[15]. Our study represents the largest targeted SNPs study in a large cohort of 500 NB patients. We investigated the association of 21 SNPs in fifteen TP53 pathway associated genes with the survival, age at diagnosis, stage of the disease, and MYCN status of 500 NB patients. We did multiple testing correction using Benjamini-Hochberg method and considered only the q-values (corrected p-values) less than 0.05 as significant. Our results showed that one SNP, CASP8 SNP D302H (rs1045485) was associated with worse overall and event-free survival in NB patients with MYCN amplification. Of note, we analyzed 41 representative blood samples from the patients in our cohort and observed that all these patients have the same SNP status in their blood as in their tumor (data not shown), indicating that CASP8 SNP D302H is a germline risk factor for NB patients with MYCN amplification. CASP8 SNP D302H is a missense variant located in exon 10 of the CASP8 gene. CASP8 encodes for inactive procaspase 8 that is later activated by dimerization to form caspase 8. Caspase 8 plays an important role in execution of the extrinsic cell death pathway (death receptor-mediated killing) and is involved in several other cell signaling pathways [28]. Down-regulation of caspase 8 is a well-characterized apoptotic defect in NB tumors [18]. CASP8 is located on 2q33, a region frequently associated with loss of heterozygosity (LOH) in NB [29]. Down-regulation of caspase 8 can be achieved via deletion or epigenetic silencing, preferentially in MYCN-amplified NB tumors [30]. A recent study has shown that loss of caspase 8 expression in a Th-MYCN/caspase-8 deleted mouse model significantly enhanced development of advanced NB and metastasis to the bone marrow [31]. Although these two events (MYCN amplification and reduced caspase 8 expression) are closely related in NB, it has been shown that MYCN overexpression does not affect the methylation status or expression of CASP8 [32]. We investigated the effect of CASP8 SNP D302H on caspase 8 mRNA expression in two independent cohorts. Our results showed that CASP8 SNP D302H had no effect on caspase 8 expression; being a non-synonymous SNP, it may therefore have a possible role in the activation of procaspase 8 or in the interaction of caspase 8 with elements that form the death-inducing signaling complex (DISC). CASP8 SNP D302H encodes an amino acid residue that is localized on the external surface of procaspase 8 and that could affect the processing of procaspase 8 or caspase-8 interaction with other apoptosis regulators. However, until this is confirmed, the effect of this SNP on caspase 8 activity remains elusive. The well-known MDM2 SNP309 was also included in our targeted analysis of 21 SNPs. Previous studies reported an association of this SNP with neuroblastoma [33], [34], and other cancers [13]; however, we have recently shown that this SNP is not associated with neuroblastoma. Ethnicity of the studied populations, statistical power, or genotyping data quality could be a reason for this discrepancy [35].

Conclusions

In conclusion our study shows that SNP D302H, a missense SNP in exon 10 of CASP8 gene, is associated with worse overall survival and event-free survival in NB patients with MYCN amplification.

Supporting Information

[Figure omitted. See PDF.]

S1 File.

Supporting tables. S1 Table, RT-qPCR primers. S2 Table, Comparison of genotyping data on amplified DNA (A) and non-amplified DNA (NA) for 18 SNPs in 31 NB cell lines. S3 Table, Stratified survival analysis of all the SNP. S4 Table, multivariate cox proportional hazards analysis (overall survival data). S5 Table, multivariate proportional hazards analysis (event-free survival data). S6 Table, CASP8 SNPD302H and age at diagnosis. S7 Table, CASP8 SNPD302H association with Stage 4 and MYCN status.

https://doi.org/10.1371/journal.pone.0114696.s001

(DOCX)

Acknowledgments

We thank Nurten Yigit for technical assistance. AR is supported by a PhD fellowship from the Ghent University research fund (BOF; 01D02210). TVM is a postdoctoral fellow of the FWO (12B9912N). RN is supported by FIS (contract PI10/15) and RTICC (contracts RD06/0020/0102; RD12/0036/0020), Instituto Carlos III Madrid & ERDF, Spain. Italian Neuroblastoma Foundation supported the Italian samples collection and preparation. This work was supported by the BOF, the FWO, the Childhood Cancer Fund, and the National Cancer Plan of the Belgian State (Action 29).

Author Contributions

Conceived and designed the experiments: AR JV FS TVM. Performed the experiments: AR. Analyzed the data: AR JV TVM. Contributed reagents/materials/analysis tools: BDW FZ GL NF GPT SC RV RN JS AE RS. Wrote the paper: AR TVM FS JV. Revised the manuscript critically for important intellectual content: GL NF GPT SC RV RN JS AE RS.

Citation: Rihani A, De Wilde B, Zeka F, Laureys G, Francotte N, Tonini GP, et al. (2014) CASP8 SNP D302H (rs1045485) Is Associated with Worse Survival in MYCN-Amplified Neuroblastoma Patients. PLoS ONE 9(12): e114696. https://doi.org/10.1371/journal.pone.0114696

References

1. Maris JM (2010) Recent advances in neuroblastoma. N Engl J Med 362:2202–2211

2. Cohn SL, Pearson ADJ, London WB, Monclair T, Ambros PF, et al. (2009) The International Neuroblastoma Risk Group (INRG) classification system: an INRG Task Force report. J Clin Oncol 27:289–297

3. Dai Z, Papp AC, Wang D, Hampel H, Sadee W (2008) Genotyping panel for assessing response to cancer chemotherapy. BMC Med Genomics 1:24

4. Grochola LF, Zeron-Medina J, Mériaux S, Bond GL (2010) Single-nucleotide polymorphisms in the p53 signaling pathway. Cold Spring Harb Perspect Biol 2:a001032

5. Levine AJ (1997) p53, the cellular gatekeeper for growth and division. Cell 88:323–331.

6. Tweddle DA, Malcolm AJ, Bown N, Pearson AD, Lunec J (2001) Evidence for the development of p53 mutations after cytotoxic therapy in a neuroblastoma cell line. Cancer Res 61:8–13.

7. Vazquez A, Bond EE, Levine AJ, Bond GL (2008) The genetics of the p53 pathway, apoptosis and cancer therapy. Nat Rev Drug Discov 7:979–987

8. Vazquez A, Grochola LF, Bond EE, Levine AJ, Taubert H, et al. (2010) Chemosensitivity Profiles Identify Polymorphisms in the p53 Network Genes 14-3-3 and CD44 That Affect Sarcoma Incidence and Survival. Cancer Res 70:172–180

9. Shoemaker RH (2006) The NCI60 human tumour cell line anticancer drug screen. Nat Rev Cancer 6:813–823

10. Garraway LA, Widlund HR, Rubin MA, Getz G, Berger AJ, et al. (2005) Integrative genomic analyses identify MITF as a lineage survival oncogene amplified in malignant melanoma. Nature 436:117–122

11. Ikediobi ON, Davies H, Bignell G, Edkins S, Stevens C, et al. (2006) Mutation analysis of 24 known cancer genes in the NCI-60 cell line set. Molecular Cancer Therapeutics 5:2606–2612

12. Paulin FEM, O'Neill M, McGregor G, Cassidy A, Ashfield A, et al. (2008) MDM2 SNP309 is associated with high grade node positive breast tumours and is in linkage disequilibrium with a novel MDM2 intron 1 polymorphism. BMC Cancer 8:281

13. Bond GL, Hu W, Bond EE, Robins H, Lutzker SG, et al. (2004) A single nucleotide polymorphism in the MDM2 promoter attenuates the p53 tumor suppressor pathway and accelerates tumor formation in humans. Cell 119:591–602

14. Boersma BJ, Howe TM, Goodman JE, Yfantis HG, Lee DH, et al. (2006) Association of breast cancer outcome with status of p53 and MDM2 SNP309. JNCI Journal of the National Cancer Institute 98:911–919

15. Wynendaele J, Bohnke A, Leucci E, Nielsen SJ, Lambertz I, et al. (2010) An Illegitimate microRNA Target Site within the 3′ UTR of MDM4 Affects Ovarian Cancer Progression and Chemosensitivity. Cancer Res 70:9641–9649

16. Brodeur GM, Pritchard J, Berthold F, Carlsen NL, Castel V, et al. (1994) Revisions of the international criteria for neuroblastoma diagnosis, staging and response to treatment. Prog Clin Biol Res 385:363–369.

17. Hellemans J, Mortier G, De Paepe A, Speleman F, Vandesompele J (2007) qBase relative quantification framework and software for management and automated analysis of real-time quantitative PCR data. Genome Biol 8:R19

18. Kamimatsuse A, Matsuura K, Moriya S, Fukuba I, Yamaoka H, et al. (2009) Detection of CpG island hypermethylation of caspase-8 in neuroblastoma using an oligonucleotide array. Pediatr Blood Cancer 52:777–783

19. Park JR, Eggert A, Caron H (2010) Neuroblastoma: biology, prognosis, and treatment. Hematol Oncol Clin North Am 24:65–86

20. Hanahan D, Weinberg RA (2011) Hallmarks of Cancer: The Next Generation. Cell 144:646–674

21. Van Maerken T, Vandesompele J, Rihani A, De Paepe A, Speleman F (2009) Escape from p53-mediated tumor surveillance in neuroblastoma: switching off the p14(ARF)-MDM2-p53 axis. Cell Death and Differentiation 16:1563–1572

22. Diskin SJ, Capasso M, Schnepp RW, Cole KA, Attiyeh EF, et al. (2012) Common variation at 6q16 within HACE1 and LIN28B influences susceptibility to neuroblastoma. Nat Genet 44:1126–1130

23. Maris JM, Mosse YP, Bradfield JP, Hou C, Monni S, et al. (2008) Chromosome 6p22 locus associated with clinically aggressive neuroblastoma. N Engl J Med 358:2585–2593

24. Capasso M, Devoto M, Hou C, Asgharzadeh S, Glessner JT, et al. (2009) Common variations in BARD1 influence susceptibility to high-risk neuroblastoma. Nat Genet 41:718–723

25. Wang K, Diskin SJ, Zhang H, Attiyeh EF, Winter C, et al. (2011) Integrative genomics identifies LMO1 as a neuroblastoma oncogene. Nature 469:216–220

26. Nguyen LB, Diskin SJ, Capasso M, Wang K, Diamond MA, et al. (2011) Phenotype restricted genome-wide association study using a gene-centric approach identifies three low-risk neuroblastoma susceptibility Loci. PLoS Genet 7:e1002026

27. Moore JH, Asselbergs FW, Williams SM (2010) Bioinformatics challenges for genome-wide association studies. Bioinformatics 26:445–455

28. Stupack DG (2013) Caspase-8 as a therapeutic target in cancer. Cancer Lett 332:133–140

29. Takita J, Yang HW, Chen YY, Hanada R, Yamamoto K, et al. (2001) Allelic imbalance on chromosome 2q and alterations of the caspase 8 gene in neuroblastoma. Oncogene 20:4424–4432.

30. Decock A, Ongenaert M, Vandesompele J, Speleman F (2011) Neuroblastoma epigenetics: from candidate gene approaches to genome-wide screenings. Epigenetics 6:962–970

31. Teitz T, Inoue M, Valentine MB, Zhu K, Rehg JE, et al. (2013) Th-MYCN mice with caspase-8 deficiency develop advanced neuroblastoma with bone marrow metastasis. Cancer Res 73:4086–4097

32. van Noesel MM, Pieters R, Voûte PA, Versteeg R (2003) The N-myc paradox: N-myc overexpression in neuroblastomas is associated with sensitivity as well as resistance to apoptosis. Cancer Lett 197:165–172.

33. Perfumo C, Parodi S, Mazzocco K, Defferrari R, Inga A, et al. (2009) MDM2SNP309 genotype influences survival of metastatic but not of localized neuroblastoma. Pediatr Blood Cancer 53:576–583

34. Perfumo C, Parodi S, Mazzocco K, Defferrari R, Inga A, et al. (2008) Impact of MDM2 SNP309 genotype on progression and survival of stage 4 neuroblastoma. Eur J Cancer 44:2634–2639

35. Rihani A, Van Maerken T, De Wilde B, Zeka F, Laureys G, et al. (2014) Lack of association between MDM2 promoter SNP309 and clinical outcome in patients with neuroblastoma. Pediatr Blood Cancer. doi:https://doi.org/10.1002/pbc.24927.

36. Stacey SN, Sulem P, Jonasdottir A, Masson G, Gudmundsson J, et al. (2011) A germline variant in the TP53 polyadenylation signal confers cancer susceptibility. Nat Genet 43:1098–1103