Genome-wide promoter methylation analysis in neuroblastoma identifies prognostic methylation biomarkers
- Anneleen Decock†1,
- Maté Ongenaert†1,
- Jasmien Hoebeeck1, 2,
- Katleen De Preter1,
- Gert Van Peer1,
- Wim Van Criekinge3, 4, 5,
- Ruth Ladenstein6,
- Johannes H Schulte7,
- Rosa Noguera8,
- Raymond L Stallings9, 10,
- An Van Damme11,
- Geneviève Laureys12,
- Joëlle Vermeulen13,
- Tom Van Maerken1, 14,
- Frank Speleman1 and
- Jo Vandesompele1Email author
© Decock et al.; licensee BioMed Central Ltd. 2012
Received: 30 May 2012
Accepted: 3 October 2012
Published: 3 October 2012
Accurate outcome prediction in neuroblastoma, which is necessary to enable the optimal choice of risk-related therapy, remains a challenge. To improve neuroblastoma patient stratification, this study aimed to identify prognostic tumor DNA methylation biomarkers.
To identify genes silenced by promoter methylation, we first applied two independent genome-wide methylation screening methodologies to eight neuroblastoma cell lines. Specifically, we used re-expression profiling upon 5-aza-2'-deoxycytidine (DAC) treatment and massively parallel sequencing after capturing with a methyl-CpG-binding domain (MBD-seq). Putative methylation markers were selected from DAC-upregulated genes through a literature search and an upfront methylation-specific PCR on 20 primary neuroblastoma tumors, as well as through MBD- seq in combination with publicly available neuroblastoma tumor gene expression data. This yielded 43 candidate biomarkers that were subsequently tested by high-throughput methylation-specific PCR on an independent cohort of 89 primary neuroblastoma tumors that had been selected for risk classification and survival. Based on this analysis, methylation of KRT19, FAS, PRPH, CNR1, QPCT, HIST1H3C, ACSS3 and GRB10 was found to be associated with at least one of the classical risk factors, namely age, stage or MYCN status. Importantly, HIST1H3C and GNAS methylation was associated with overall and/or event-free survival.
This study combines two genome-wide methylation discovery methodologies and is the most extensive validation study in neuroblastoma performed thus far. We identified several novel prognostic DNA methylation markers and provide a basis for the development of a DNA methylation-based prognostic classifier in neuroblastoma.
Neuroblastoma (NB) is a neuroectodermal tumor that originates from precursor cells of the sympathetic nervous system and represents the most common extra-cranial solid tumor of early childhood. NB displays a highly variable clinical course, ranging from spontaneous regression to life-threatening disease .
Despite advances in multimodal anti-cancer therapies, survival rates for children with aggressive NB remain disappointingly low. Survival rates vary widely, depending on clinical features, such as age at diagnosis and tumor stage, as well as biological characteristics of the tumor. Amongst the latter, MYCN amplification has been used for many years as a genetic marker for therapy stratification . More recently, a subset of high-risk tumors with non-amplified MYCN and 11q deletions was identified, while absence of segmental aberrations upon genome-wide DNA copy number analysis was found to be associated with excellent survival [2, 3]. In order to facilitate the comparison of risk-based clinical trials, a new consensus approach for pretreatment risk classification has been designed including genetic parameters [1, 4]. Despite this progress, additional markers for therapeutic stratification are warranted in order to avoid under- or overtreatment and to improve selection of ultra-high-risk patients for new experimental therapies. Recently, prognostic mRNA and microRNA (miRNA) signatures were developed to accommodate this need [5–7]. Here, we propose that the use of DNA methylation markers is a new and promising method for prognostic classification.
DNA methylation is the addition of a methyl group to carbon 5 of the cytosine within the dinucleotide CpG. Dense clusters of CpG dinucleotides, termed CpG islands, are often present in gene promoters and methylation of those regions typically results in transcriptional silencing of the gene. As such, abnormal DNA methylation in cancer cells leads to aberrant expression patterns . In NB, the most described epigenetic alterations are DNA methylation of CASP8  and RASSF1A , both associated with risk factors, such as MYCN amplification (MNA), age at diagnosis and tumor stage [11–15]. Recently, a few genome-wide methylation screening methodologies have been applied in NB, including re-expression analysis after treatment with 5-aza-2'-deoxycytidine (DAC), DNA methylation promoter arrays after capturing with methylated DNA immunoprecipitation (MeDIP) and methylation microarrays. These studies indicate that aberrant DNA methylation makes an important contribution toward NB tumor biology by downregulating specific genes and show the potential of using DNA methylation in future patient therapy stratification protocols [16–18]. Furthermore, the power of DNA methylation as a non-invasive, sensitive and specific biomarker has been demonstrated by measuring DNA methylation of RASSF1A in serum of primary NB patients  (for a detailed review see ). In order to improve the outcome prediction of NB patients, this study aims at establishing robust DNA methylation biomarkers that can identify patients with unfavorable prognosis.
Discovery and integrated analysis: genome-wide methylation screening for selection of candidate biomarkers
The second genome-wide DNA methylation screening methodology we applied, to the same eight NB cell lines, was MBD-seq: massively parallel sequencing of methylation-enriched DNA fragments, whereby the enrichment is based on the capture of methylated sheared DNA using the high affinity of the methyl-CpG-binding domain () of the protein MBD2 towards methylated cytosines. Sequencing yielded 4.4 to 8.6 million paired-end reads, depending on the cell line, and after peak calling 70,816 to 112,412 peaks were detected, representing genomic regions methylated in the corresponding cell line. Between 7,612 and 11,178 of these peaks (around 10% of all identified peaks) are located in promoter regions of annotated genes (-1,500 bp to +1,000 bp around the transcriptional start site (TSS)). These 'methylation peaks' were visualized in the Integrative Genomic Viewer , showing that promoter regions that are well known to be heavily methylated in NB were confirmed - for example, the protocadherin β (PCDHB) family cluster (Additional file 2) [23, 24]. In some regions (for example, in the promoter regions of HIST1H3C and ACSS3) it was also possible to distinguish different DNA methylation profiles between MNA (IMR-32, LAN-2 and N206) and MYCN single copy (SH-SY5Y, SK-N-AS, CLB-GA and SJNB-1) NB cell lines (Additional file 2). Using the R/BioC package DESeq , 510 regions were identified as differentially methylated between MYCN amplified and single copy cell lines, of which 95 are in close proximity to an annotated TSS (-1,500 bp to +1,000 bp). Also, some miRNAs appeared to be methylated in their promoter region.
After peak calling, we also performed gene set enrichment analysis , using a custom, ranked list of genes with at least one MBD peak present in a region -1,500 bp to +500 bp around its TSS, in order to explore whether promoter regions that are enriched after MBD capture are often re-expressed as well upon DAC treatment. This analysis clearly showed a high enrichment score for each cell line (enrichment scores from 0.32 to 0.36; FDR q-value <0.01), demonstrating that a large portion of methylated regions (captured by MBD) are indeed reactivated upon DAC treatment. The overlap between the two genome-wide datasets can be further explored by intersecting them. In total, 183 genes are both re-activated upon DAC treatment (at least 1 log2 difference after and before treatment) and have an MBD peak in their promoter regions (-1,500 bp to +1,000 bp around the TSS) in at least 2 of the 8 investigated NB cell lines. Of these 183 genes, 46 are both re-expressed and methylated in 3 cell lines, 9 in 4 cell lines and 5 in at least 5 cell lines.
As we feared that only using cell lines in the selection phase of potential prognostic DNA methylation biomarkers would lead to the identification of methylated markers not necessarily related with prognosis, six publicly available mRNA expression studies [27–34] were included in the analysis. In these studies, which comprise mRNA expression data of 380 primary NB tumors, identifying differentially expressed probes (genes) between prognostic groups would allow us to pinpoint potential prognostic methylated promoter regions in our methylome maps. Finally, a rank-based scoring system was used to prioritize genes that show methylation, re-expression after DAC treatment and differential expression (related with risk) across the prognostic groups. This score scheme uses the individual ranks of each analysis. In brief, DAC reactivation is ranked according to FDR rate (as determined by RankProd analysis), MDB-seq data are ranked according to peak P-values and expression data are ranked according to FDR (determined by RankProd analysis). Each data source is given the same weight and a combined rank is calculated (for details, see Materials and methods). This scoring system combined all generated data and allowed us to select 43 top-ranking and thus strong prognostic methylation candidate genes without the need to use rather artificial threshold values for the different datasets.
Validation: determining the prognostic power of DNA-methylation biomarkers
For these 43 genes, 48 MSP assays were designed and tested on 3 NB cell lines (IMR-32, SK-N-AS and SH-SY5Y) and the HCT-116 DKO cell lines, along with an independent cohort of 89 primary NB samples. Within the 89 primary NB sample set all three prognostic groups (LR-SURV, HR-DOD and high-risk survivors (HR-SURV); for details see Material and methods) were approximately equally represented. The complete matrix with all MSP results of all samples and a global overview of the MSP results per assay can be found in Additional file 3. Over 60% of the designed assays indeed detected methylation for the respective marker in at least 10% of the selected NB tumors. Ten MSP assays (COL6A3, miR-1225, miR-3177, PCDHA6, PLXNC1, ANKRD43, ADRB2, APOE, miR-671 and QPCT) revealed methylation in at least 75% of the patient samples, and the MSP assays for KCND2, PRPH, KRT19 (assay 83159) and TNFRSF10D were methylated in 50% to 75% of the patient samples. We could also detect DNA methylation in the promoter region of miR-1225, miR-3177, miR-671 and miR-663, methylated in 99%, 99%, 79% and 4% of the patient samples, respectively.
Several individual markers are differentially methylated between the prognostic groups and neuroblastoma risk factors
MYCN single copy
Age at diagnosis > 12 months
Age at diagnosis < 12 months
Age at diagnosis > 18 months
Age at diagnosis < 18 months
Statistics (Fisher's exact P-value)
Age cutoff 12 months
Age cutoff 18 months
Some individual MSP assays were also associated with one or more NB risk factors (stage, MYCN status and age at diagnosis), and are thus potential prognostic biomarkers in NB (Table 1). In this analysis, the age at diagnosis was tested using two different age cutoffs. The 12 months cutoff was chosen as it was used for therapy stratification and as a criterion in the sample selection. The more recently established cutoff of 18 months [1, 35, 36] was also taken into account. Newly discovered methylated markers are FAS, PRPH, CNR1, QPCT, HIST1H3C, ACSS3 and GRB10, methylation of which is associated with at least one of the NB risk factors. Table 1 further indicates that the difference in the methylation status of HIST1H3C and ACSS3 between MYCN single copy and MNA NB cell lines as detected by MBD-seq is reflected in the MSP results of the primary tumors as well, as HIST1H3C and ACSS3 are almost exclusively methylated in MNA samples.
Several individual DNA methylation markers are associated with survival
HR-DOD versus LR-SURV (P-value OS)
HR-DOD versus LR-SURV (P-value EFS)
HR-DOD versus LR-SURV and HR-SURV (P-value OS)
HR-DOD versus LR-SURV and HR-SURV (P-value EFS)
HR-DOD versus LR-SURV and HR-SURV (P-value OS (stratum))
Age < 12 months 0.035
Stage 4 0.033
MYCN = 0 0.033
HR-DOD versus LR-SURV and HR-SURV (P-value EFS (stratum))
Age < 12 months 0.014
MYCN = 0 0.001
mRNA expression profiling: determining transcriptional silencing of DNA methylation biomarkers
mRNA expression level of several markers associates with neuroblastoma risk factors, prognostic groups and survival
Age cutoff 12 months
Age cutoff 18 months
Overall survival status
Survival analysis using Cox proportional hazards further shows that low mRNA expression levels of CNR1 (hazard ratio (HR) 0.768; 95% CI 0.619 to 0.953; P = 0.028), GRB10 (HR 0.613; 95% CI 0.433 to 0.866; P = 0.015) and PRPH (HR 0.714; 95% CI 0.566 to 0.922; P = 0.015) were significantly associated with poor survival. After dichotomization of the mRNA expression data, using the median relative mRNA expression value as a cutoff, Kaplan-Meier survival curves were plotted (log-rank test; Additional file 5).
Thus far, most of the studies analyzing DNA methylation patterns in NB have been candidate gene-based, with the methylation status of the promoter region for only a limited number of genes being tested. These candidate genes were selected based either on prior knowledge of NB tumor biology or on the fact of being methylated in other tumor types. As a consequence, only few DNA methylation biomarkers, such as KRT19, TNFRSF10D, CASP8, ZMYND10 and RASSF1A, were previously related with NB risk factors or survival [11, 13–15, 18, 37–41]. In order to identify new DNA methylation biomarkers in NB, we applied a multilevel experimental approach. In the discovery phase we established a genome-wide methylome map of eight NB cell lines. These cell lines were profiled using gene expression microarrays before and after DAC treatment, and using MBD capture followed by next-generation sequencing (NGS). The combination of both methodologies enabled the identification of regions that are both methylated and undergo re-expression upon DAC treatment. So far, only MeDIP chips were used in whole promoter profiling studies on NB , making this study the first one using NGS for unbiased and more sensitive assessment of genome-wide DNA methylation patterns in NB. Our results emphasize the potential of this epigenetic sequencing technique, as it enables the investigation of the methylome or epigenome of a sample in great detail at a feasible cost.
Integration of these methylome maps with genome-wide gene expression profiles led to a selection of 43 candidate biomarkers that were tested on 89 primary NB patient samples. All samples were assigned to one of three discrete prognostic patient groups (low-risk survivors (LR-SURV), high-risk deceased (HR-DOD) and high-risk survivors (HR-SURV)). While most NB methylation studies did not discriminate between HR-SURV and HR-DOD patients, we believe this is an important clinical question, as both prognostic groups are currently considered high-risk and uniformly treated, making the present study unique in its concept. As we make use of amplified bisulfite-converted DNA, only limited amounts (100 to 200 ng) of tumor DNA are required to test over 100 MSP assays. The MBD-seq results greatly help in designing the assays in the most informative regions, which is important as the assay location is critically important, again confirmed in this study for a number of genes for which multiple assays were designed (for example, TGFBI and KRT19). The combination of the number of samples and assays used in this study further makes it the most comprehensive methylation study in NB. Furthermore, the high-throughput validation pipeline allows fast and accurate follow-up validation of potential candidate DNA methylation biomarkers for large numbers of patients. Indeed, PCR-based detection methodologies are robust and can thus be used in a wide range of laboratory settings for a low price without the need of special equipment other than for qPCR and (microfluidic) electrophoresis, both present in most molecular laboratories. The presented DNA methylation screening and validation methodology can thus easily be adapted by (cancer) researchers addressing similar questions in other research fields.
In this study, several novel biomarkers were established in addition to known DNA methylation biomarkers in NB, such as KRT19, TGFBI, TNFRSF10D and TNFRSF10A [14, 18, 37, 42, 43]. Interestingly, some of these novel genes were previously reported to be important in NB biology (without reference to their epigenetically altered status) or were described as epigenetic biomarkers in other tumor entities, such as FAS, which encodes a member of the tumor necrosis factor receptor (TNFR) superfamily [44–50]. Several other novel methylation biomarkers were also shown to be differentially methylated between HR and LR patients, and many of these were associated with NB risk factors or with survival. However, discriminating HR-DOD and HR-SURV patients is challenging. While only a few individual MSP designs (HIST1H3C, KRT19 and ACSS3) were moderately discriminatory between these two HR subgroups, the combined methylation data analysis of CNR1, ACSS3, HIST1H3C and PRPH indicates the potential of DNA methylation biomarkers in stratifying HR NB patients. In this study, the difficulty of identifying individual biomarkers that differentiate between HR-DOD and HR-SURV patients may be explained by the fact that NB cell lines were used in the discovery phase, thus enriching for genes discriminating between HR and LR patients as NB cell lines can be considered models for aggressive HR tumors. To accommodate this, we plan to perform a large-scale discovery using MBD capture followed by NGS on primary NB tumors equally distributed over the three prognostic groups used here.
PRPH is one of the novel biomarkers identified and is differentially methylated across the prognostic patient groups. This gene encodes the cytoskeletal protein peripherin found in neurons of the peripheral nervous system, and is probably associated with maturation of the neuronal phenotype and hence serves as a differentiation marker for tumors derived from the neural crest . In our study, PRPH methylation was mainly detected in more advanced tumor stages. Since promoter methylation may cause transcriptional silencing of the gene and advanced NB tumor stages are less differentiated , this is in line with the idea that high levels of peripherin contribute to more differentiated tumor stages. As demonstrated in this study, this idea is further strengthened by the fact that PRPH mRNA expression levels gradually decreased with increasing aggressiveness of the tumor. As whole genome sequence analysis recently showed that genes involved in neuritogenesis are recurrently affected in high-stage NB , the identification of PRPH methylation opens new research perspectives regarding NB therapy.
Next to protein-coding genes, some MSP assays were designed in the promoter region of miRNAs. Aberrant miRNA expression contributes majorly to NB tumor biology and has been extensively studied during the past few years. Most of these studies used miRNA microarrays or high-throughput RT-qPCR to analyze the miRNA expression patterns in primary NB tumor samples [54–56]. Although a broad deregulation of the miRNA expression profile in NB has been described, miRNA promoter hypermethylation is relatively unexplored. Up until now, the only miRNA for which the promoter region is known to be methylated in NB is miR-200b . Interestingly, miR-1225, miR-3177 and miR-671 were found to be methylated in their promoter region in more than 75% of the NB tumors in our study. Currently, little is known about the putative function of these miRNAs, as they are not well described or not described at all in the literature [58, 59].
Another interesting finding is that MYCN single copy and MNA samples show differential promoter methylation of HIST1H3C and ACSS3. Currently, little is known about the association between MYCN and DNA methylation of certain genes in NB, nor about the underlying molecular mechanisms. Previously, Teitz et al. [9, 60] showed that DNA methylation of CASP8 is almost exclusively associated with MNA in both NB cell lines and primary tumors. They further noticed that CASP8 was hemi-methylated (only one allele) in stage 1, 2 and 3 NB, which may indicate that complete methylation of CASP8 may be coupled to another event, such as amplification of the MYCN gene. While this suggests that MNA is functionally linked to complete methylation of both CASP8 alleles, it is not clear if these two events occur concurrently, or if one event leads to the other. Obviously, genes differentially methylated between MNA and MYCN single copy samples need to be further functionally characterized, as this may lead to new insights into NB biology.
Although international collaboration in the field of NB has invested tremendous effort in optimizing patient stratification and therapy protocols, OS rates remain low. This study shows that DNA methylation biomarkers have the potential to refine current risk assessment schemes. In contrast to most NB methylation studies that are candidate gene-based, we applied two genome-wide detection methodologies to discover hypermethylated regions in NB: re-expression analysis after demethylating DAC treatment and NGS after MBD capture. Furthermore, we present a high-throughput and semi-automated MSP pipeline, which was used to test the candidate DNA methylation markers on a large patient tumor cohort. We have identified novel aberrant promoter hypermethylation of protein coding genes and miRNAs in NB. Some of these DNA methylation biomarkers are associated with NB risk factors and/or survival, emphasizing the prognostic value of these markers and their potential to be used in a DNA methylation-based prognostic classifier in NB. The use of such a DNA methylation signature, discriminating HR patients, is demonstrated here by the combined methylation data analysis of CNR1, ACSS3, HIST1H3C and PRPH. Furthermore, some DNA methylation biomarkers showed low levels of mRNA expression in patient groups with high methylation levels. This suggests that promoter methylation may contribute to transcriptional silencing of these genes, which may be important in the pathogenesis of NB. Encouraged by these results, we will now extensively further validate these DNA methylation biomarkers and refine the methylome map of different prognostic NB patient groups.
Materials and methods
Neuroblastoma cell lines and primary tumors
In total, 33 well-characterized NB cell lines, authenticated using array comparative genomic hybridization and short tandem repeat genotyping, were included in this study (Additional files 2 and 4). DNA was isolated using the QIAamp DNA Mini Kit (Qiagen, Venlo, The Netherlands). In addition, 109 primary tumor samples of NB patients were collected prior to therapy at the Ghent University Hospital (Ghent, Belgium), the University Children's Hospital Essen (Essen, Germany), Our Lady's Children's Hospital Dublin (Dublin, Ireland) or the Hospital Clínico Universitario (Valencia, Spain). Informed consent was obtained from each patient's guardian and the study was approved by the ethical committee of the Ghent University Hospital (approval number B67020109912). Clinical characteristics of the patients are shown in Additional files 1 and 3. All NB patient samples were assigned to one of three defined risk groups based on risk parameters (tumor stage, MYCN status and age at diagnosis) and disease outcome. First, HR patients that died of disease (HR-DOD) as defined by stage 2/3, MNA, DOD; stage 4, age at diagnosis <12 months, MNA, DOD; or stage 4, age at diagnosis >12 months, DOD (n = 39). Second, HR patients alive (HR-SURV) after follow-up time >1,000 days (n = 30). Third, LR patients alive (LR-SURV) defined by stage 1/2, MYCN single copy, follow-up time >1,000 days; stage 3, MYCN single copy, age <12 months, follow-up time >1,000 days (status at last known follow-up is alive; n = 40). The clinical data of the 366 primary NB tumors (SIOPEN/GPOH cDNA library ), used to test the mRNA expression levels of the most promising DNA methylation biomarkers, can be found in Additional file 5.
Microarray after re-expression analysis
Eight NB cell lines (CHP-902R, CLB-GA, IMR-32, LAN-2, N206, SH-SY5Y, SK-N-AS and SJNB-1) were grown in the presence of 3 µM DAC (Sigma, Bornem, Belgium)) for 3 days, as previously described, and untreated controls were also prepared . After harvesting, RNA was extracted with the RNeasy Mini kit (Qiagen), accompanied by RNase free DNase treatment on column (Qiagen). After RNA quality check on the Experion (Bio-Rad, Nazareth, Belgium), sample preparation, hybridization to Affymetrix Human Genome U133 Plus 2.0 oligonucleotide chips and scanning were carried out according to the manufacturer's protocol at the VIB MicroArray Facility. Standard quality metrics (simpleaffy BioC package  - boxplots, visual inspection of the slides, 5'-3' degradation plot) demonstrated that the oligonucleotide chip data were of good quality. The BioC affy package was used to normalize (gc-RMA normalization) the expression levels and to obtain present/absent (expression/no expression) MAS 5.0 calls for each probe set. For all cell lines and for each probe set, the number of reactivation events was counted (absent in untreated cells and present in treated cells). Expression data (before and after DAC treatment) have been deposited into the Gene Expression Omnibus [GEO: GSE31229], according to the MIAME guidelines.
DNA samples (1 µg DNA) of the eight NB cell lines were sheared (Covaris S2) to an average length of 200 bp. Fragment distribution was determined by the Agilent 2100 Bioanalyzer and the concentration was determined using the Quant-iT PicoGreen dsDNA HS Assay Kit (Invitrogen, Ghent, Belgium). Starting from 200 ng sheared DNA, the MethylCollector Kit (ActiveMotif, La Hulpe, Belgium) was used to enrich for methylated fragments. Library preparation for multiplex Illumina sequencing was done by combining the DNA Sample Prep Master Mix Set 1 (New England Biolabs, Frankfurt am Main, Germany) and the Multiplexing Sample Preparation Oligo Kit (Illumina). Size selection of the library was done on a 2% agarose gel. Fragments of around 300 bp (±50 bp) were excised and purified. Illumina library amplification (21 cycles) was performed and concentration was determined. Paired-end sequencing was used for high confidence mapping of captured fragments (2 × 45 bp sequencing - Illumina GAIIx, NXTGNT). Paired-end reads were mapped on the human reference genome (GRCh37) using Bowtie 0.12.7 and peaks were called using MACS 1.4beta. For differential methylation analysis, PCR duplicates were removed and sequence tags counted by using the BioC packages ShortRead and rtracklayer [63, 64]. Sequence tag counts per sample were used to compose a count matrix that could be processed by the BioC package DESeq . Sequencing data (raw sequence files, WIG files for visualization of the mapping results and the BED peak files as determined by MACS) have been deposited into GEO [GEO:GSE31353].
Selection of candidate biomarkers
Initial high-throughput MSP
In total, 212 MSP assays (Additional file 1) were designed in the promoter region of 120 corresponding genes re-expressed after DAC treatment, and tested on both the DAC-treated and untreated NB cell lines, 9 LR-SURV patients and 11 HR-DOD patients (Additional file 1). A total of 500 to 1,000 ng DNA of these samples was bisulfite-treated (EZ DNA Methylation Kit, Zymo Research, Irvine, CA, USA), eluted in 30 µl elution buffer and then tested on the BioTrove OpenArray (Life Technologies, Ghent, Belgium). Beta actin (ACTB) was used as a control and to normalize samples. The in vitro methylated HCT-116 DKO cell line (treated with SssI, Zymo Research) was used as a positive control. The methylation status for each MSP assay was determined, and called methylated if the melting temperature (Tm) of the amplicon was within a specific interval as defined by the positive control sample. These methylation calls were further analyzed by determining specificity and sensitivity of the HR-DOD samples versus LR-SURV samples.
Publicly available mRNA expression studies
Six publicly available mRNA expression studies [27–34] [GEO:GSE19274, GEO:GSE16237, GEO:GSE14880, GEO:GSE12460, GEO:GSE13136, GEO:GSE3960] were analyzed using RankProd analysis (BioC package ), to identify differentially expressed probes between prognostic groups (high-risk versus low-risk, high-stage versus low-stage, and MNA versus MYCN single copy).
Each analysis score of a promoter region (for example, RankProd FDR value and P-value for differential expression between risk groups, and P-values of the peak after MBD-seq) was ranked and given a score, ranging from tan(1) to 0 according to their rank. These individual scores were then summed and 43 top-ranking genes were selected for further analysis.
MSP assays were designed to only amplify the bisulfite-converted target region of interest and do not anneal to genomic DNA. As each primer contains at least two CpG sites, this means that a PCR product will only be generated if the template is methylated. We choose not to design the according U primers (that would amplify the non-methylated bisulfite-converted DNA) as we do not assess methylation in a quantitative way. After in silico assay evaluation, 48 selected MSP primers (including the ACTB control; Additional file 3) were empirically validated on the Roche LightCycler 480 (LC480) using the in vitro methylated HCT-116 DKO (positive control), the HCT-116 DKO (negative control) and NB cell lines. Based on melting curve and amplicon size analysis, all assays were considered amplicon specific. The MSP assays were tested on 89 samples, selected from the previously described patient groups (31 LR-SURV patients, 28 HR-DOD patients and 30 HR-SURV patients; Additional file 3). A no template control (NTC) sample was loaded as well. For all samples, 500 to 1,000 ng DNA was bisulfite-treated (EZ DNA Methylation Kit, Zymo Research) and eluted in 40 µl elution buffer. Prior to MSP, bisulfite-treated DNA (BT-DNA) was amplified using the EpiTect Whole Bisulfitome Kit (Qiagen), starting from 100 ng BT-DNA. After amplification, the yield was determined by the Qubit 2.0 fluorometer in combination with the Quant-iT PicoGreen dsDNA BR Assay Kit (Invitrogen). The MSP was performed on the LC480 and plates were prepared using the Tecan freedom Evo robot, using a design that assures that all samples were tested for the same assay in the same run . MSP amplifications were performed in 10 µl containing 5 µl LC480 SYBR Green I Master Mix (2×; Roche, Vilvoorde, Belgium), 1 mg/ml bovine serum albumin (Roche), 1 mM MgCl2 (Roche), 125 nM forward and reverse primer (IDT, Leuven, Belgium), sample (20 ng amplified BT-DNA) and nuclease-free water (Sigma). MSP conditions were as follows: activation for 10 minutes at 95°C, 45 amplification cycles (10 s at 95°C, 30 s at 60°C and 5 s at 72°C), followed by melting curve analysis (5 s at 95°C - melting curve from 60 to 95°C) and cool down to 45°C. Afterwards, the size of the amplicons was determined using the Caliper LabChip GX. A MSP assay was considered methylated if (1) its Cq value <35 (calculated by the LC480 software using the second derivative maximum method), (2) its melting temperature (Tm) differed no more than 2°C from that of the positive control sample, and (3) the amplicon length differed no more than 10 bp from the band size of the positive control sample. In addition, the band height, as determined by the LabChip GX software, was required to be higher than 20.
mRNA expression profiling
The mRNA expression levels of CNR1, GRB10, KRT19, PRPH and QPCT were profiled on the NB SIOPEN/GPOH cDNA library generated from 366 primary NB tumor samples (Additional file 5) . For each DNA methylation marker a qPCR mRNA assay was designed and validated in silico and in vitro (Additional file 5) . PCR plates were prepared as described in the previous section and RT-qPCR was performed on the LC480 as described in . Relative gene expression levels were then normalized using the geometric mean of five reference sequences (HPRT1, SDHA, UBC, HMBS and AluSq) . For HIST1H3C and ACSS3, a qPCR mRNA assay (Additional file 4) was designed and tested on 31 NB cell lines on which the corresponding MSP assay was tested as well. Here, qPCR amplifications were performed in 5 µl containing 2.5 µl SsoAdvanced SYBR Green Supermix (2×; Bio-Rad), 0.25 µl forward and reverse primer (5 µM each) and 2 µl cDNA sample (corresponding to 5 ng cDNA). Relative gene expression levels were normalized using the geometric mean of the reference sequences SDHA, UBC and AluSq. All RT-qPCR data analysis was done in qbasePLUS version 2.0 (Biogazelle, Ghent, Belgium) . Logged and normalized qPCR data can be found in Additional file 4 and 5.
Statistical analyses were performed using IBM SPSS software version 19.0. All statistical tests were two-sided and P-values <0.05 were considered statistically significant. Differential methylation across the prognostic groups was determined by the Chi square test. The relationship between the methylation status and NB risk factors was determined using Fisher's exact test. Univariate survival analysis was performed with the Kaplan-Meier method and log-rank statistics (Mantel-Cox) to determine the impact of methylation status on EFS and OS. EFS was defined as the time between initial diagnosis and relapse or death of disease, or time between diagnosis and last follow-up if no event had occurred. OS is the time to disease-related death or last follow-up. Hierarchical clustering and leave-one-out decision tree analysis were performed using R 2.13.0 (rpart package). The relationship between logged mRNA expression levels and the prognostic groups, OS status and NB risk factors was determined using the nonparametric Kruskal-Wallis test or Mann-Whitney test. Hazard ratios between logged mRNA expression data and survival were estimated using the Cox proportional hazard model. Kaplan-Meier curves were created by dichotomizing the logged mRNA expression data, using the median mRNA expression value as a cutoff. For HIST1H3C, the relationship between logged mRNA expression levels and the methylation status of the gene, and the MYCN status, was determined using the Mann-Whitney test. For all the above mentioned statistical tests, multiple hypothesis testing correction was performed (Benjamini-Hochberg method by using the R function p.adjust).
false discovery rate
Gene Expression Omnibus
high-risk deceased patients
massively parallel sequencing after capturing with an MBD
methylated DNA immunoprecipitation
transcriptional start site.
The authors thank Els De Smet, Nurten Yigit, Justine Nuytens, Johan Vandersmissen and Jean-Pierre Renard (MDxHealth) for their kind help and NXTGNT for the excellent sequencing services. AD, KDP and TVM are supported by a grant of the Research Foundation Flanders (FWO), and GVP by a grant of the Ghent University Special Research Fund (BOF). RLS was a recipient of grants from Science Foundation Ireland (07/IN.1/B1776), Children's Medical and Research Foundation, and National Institutes of Health (NIH; 5R01CA127496). RN was a recipient of grants (RD06/0020/0102) from Thematic Network of Cooperative Research on Cancer (RTICC), Institute of Health Carlos III (ISCIII) and European Regional Development Fund (ERDF), and Spanish Health Research Fund (FIS; 2010/PI015). This study was further supported by Emmanuel Van der Schueren Foundation (scientific partner of the Flemish League Against Cancer (VLK)), Fournier-Majoie Foundation for Innovation (FFMI), UGent GOA (grant number 01G01910), Institute for the Promotion of Innovation by Science and Technology in Flanders (IWT) and FWO (grant number G0530.12N).
- Øra I, Eggert A: Progress in treatment and risk stratification of neuroblastoma: impact on future clinical and basic research. Semin Cancer Biol. 2011, 21: 217-228. 10.1016/j.semcancer.2011.07.002.PubMedView ArticleGoogle Scholar
- Schleiermacher G, Michon J, Ribeiro A, Pierron G, Mosseri V, Rubie H, Munzer C, Bénard J, Auger N, Combaret V, Janoueix-Lerosey I, Pearson A, Tweddle DA, Bown N, Gerrard M, Wheeler K, Noguera R, Villamon E, Cañete A, Castel V, Marques B, de Lacerda A, Tonini GP, Mazzocco K, Defferrari R, de Bernardi B, di Cataldo A, van Roy N, Brichard B, Ladenstein R, et al: Segmental chromosomal alterations lead to a higher risk of relapse in infants with MYCN-non-amplified localised unresectable/disseminated neuroblastoma (a SIOPEN collaborative study). Br J Cancer. 2011, 105: 1940-1948. 10.1038/bjc.2011.472.PubMedPubMed CentralView ArticleGoogle Scholar
- Vandesompele J, Baudis M, De Preter K, Van Roy N, Ambros P, Bown N, Brinkschmidt C, Christiansen H, Combaret V, Lastowska M, Nicholson J, O'Meara A, Plantaz D, Stallings R, Brichard B, Van den Broecke C, De Bie S, De Paepe A, Laureys G, Speleman F: Unequivocal delineation of clinicogenetic subgroups and development of a new model for improved outcome prediction in neuroblastoma. J Clin Oncol. 2005, 23: 2280-2299. 10.1200/JCO.2005.06.104.PubMedView ArticleGoogle Scholar
- Cohn SL, Pearson ADJ, London WB, Monclair T, Ambros PF, Brodeur GM, Faldum A, Hero B, Iehara T, Machin D, Mosseri V, Simon T, Garaventa A, Castel V, Matthay KK: The International Neuroblastoma Risk Group (INRG) classification system: an INRG Task Force report. J Clin Oncol. 2009, 27: 289-297. 10.1200/JCO.2008.16.6785.PubMedPubMed CentralView ArticleGoogle Scholar
- Schulte JH, Schowe B, Mestdagh P, Kaderali L, Kalaghatgi P, Schlierf S, Vermeulen J, Brockmeyer B, Pajtler K, Thor T, de Preter K, Speleman F, Morik K, Eggert A, Vandesompele J, Schramm A: Accurate prediction of neuroblastoma outcome based on miRNA expression profiles. Int J Cancer. 2010, 127: 2374-2385. 10.1002/ijc.25436.PubMedView ArticleGoogle Scholar
- Vermeulen J, De Preter K, Naranjo A, Vercruysse L, Van Roy N, Hellemans J, Swerts K, Bravo S, Scaruffi P, Tonini GP, De Bernardi B, Noguera R, Piqueras M, Cañete A, Castel V, Janoueix-Lerosey I, Delattre O, Schleiermacher G, Michon J, Combaret V, Fischer M, Oberthuer A, Ambros PF, Beiske K, Bénard J, Marques B, Rubie H, Kohler J, Pätschger U, Ladenstein R, et al: Predicting outcomes for children with neuroblastoma using a multigene-expression signature: a retrospective SIOPEN/COG/GPOH study. Lancet Oncol. 2009, 10: 663-671. 10.1016/S1470-2045(09)70154-8.PubMedPubMed CentralView ArticleGoogle Scholar
- De Preter K, Mestdagh P, Vermeulen J, Zeka F, Naranjo A, Bray I, Castel V, Chen C, Drozynska E, Eggert A, Hogarty MD, Izycka E, London WB, Noguera R, Piqueras M, Bryan K, Schowe B, van Sluis P, Molenaar JJ, Schramm A, Schulte JH, Stallings RL, Versteeg R, Laureys G, Van Roy N, Speleman F, Vandesompele J: miRNA expression profiling enables risk stratification in archived and fresh neuroblastoma tumor samples. Clin Cancer Res. 2011, 17: 7684-7692. 10.1158/1078-0432.CCR-11-0610.PubMedPubMed CentralView ArticleGoogle Scholar
- Portela A, Esteller M: Epigenetic modifications and human disease. Nat Biotechnol. 2010, 28: 1057-1068. 10.1038/nbt.1685.PubMedView ArticleGoogle Scholar
- Teitz T, Wei T, Valentine MB, Vanin EF, Grenet J, Valentine VA, Behm FG, Look AT, Lahti JM, Kidd VJ: Caspase 8 is deleted or silenced preferentially in childhood neuroblastomas with amplification of MYCN. Nat Med. 2000, 6: 529-535. 10.1038/75007.PubMedView ArticleGoogle Scholar
- Astuti D, Agathanggelou A, Honorio S, Dallol A, Martinsson T, Kogner P, Cummins C, Neumann HP, Voutilainen R, Dahia P, Eng C, Maher ER, Latif F: RASSF1A promoter region CpG island hypermethylation in phaeochromocytomas and neuroblastoma tumours. Oncogene. 2001, 20: 7573-7577. 10.1038/sj.onc.1204968.PubMedView ArticleGoogle Scholar
- Michalowski MB, de Fraipont F, Plantaz D, Michelland S, Combaret V, Favrot M-C: Methylation of tumor-suppressor genes in neuroblastoma: The RASSF1A gene is almost always methylated in primary tumors. Pediatr Blood Cancer. 2008, 50: 29-32. 10.1002/pbc.21279.PubMedView ArticleGoogle Scholar
- Yang Q, Liu S, Tian Y, Hasan C, Kersey D, Salwen HR, Chlenski A, Perlman EJ, Cohn SL: Methylation-associated silencing of the heat shock protein 47 gene in human neuroblastoma. Cancer Res. 2004, 64: 4531-4538. 10.1158/0008-5472.CAN-04-0956.PubMedView ArticleGoogle Scholar
- Grau E, Martinez F, Orellana C, Canete A, Yañez Y, Oltra S, Noguera R, Hernandez M, Bermudez JD, Castel V: Hypermethylation of apoptotic genes as independent prognostic factor in neuroblastoma disease. Mol Carcinog. 2010, 50: 153-162.PubMedView ArticleGoogle Scholar
- Banelli B, Gelvi I, Di Vinci A, Scaruffi P, Casciano I, Allemanni G, Bonassi S, Tonini GP, Romani M: Distinct CpG methylation profiles characterize different clinical groups of neuroblastic tumors. Oncogene. 2005, 24: 5619-5628. 10.1038/sj.onc.1208722.PubMedView ArticleGoogle Scholar
- Misawa A, Tanaka S, Yagyu S, Tsuchiya K, Iehara T, Sugimoto T, Hosoi H: RASSF1A hypermethylation in pretreatment serum DNA of neuroblastoma patients: a prognostic marker. Br J Cancer. 2009, 100: 399-404. 10.1038/sj.bjc.6604887.PubMedPubMed CentralView ArticleGoogle Scholar
- Murphy DM, Buckley PG, Bryan K, Das S, Alcock L, Foley NH, Prenter S, Bray I, Watters KM, Higgins D, Stallings RL: Global MYCN transcription factor binding analysis in neuroblastoma reveals association with distinct E-box motifs and regions of DNA hypermethylation. PLoS One. 2009, 4: e8154-10.1371/journal.pone.0008154.PubMedPubMed CentralView ArticleGoogle Scholar
- Buckley PG, Das S, Bryan K, Watters KM, Alcock L, Koster J, Versteeg R, Stallings RL: Genome-wide DNA methylation analysis of neuroblastic tumors reveals clinically relevant epigenetic events and large-scale epigenomic alterations localized to telomeric regions. Int J Cancer. 2010, 128: 2296-2305.View ArticleGoogle Scholar
- Caren H, Djos A, Nethander M, Sjoberg R-M, Kogner P, Enstrom C, Nilsson S, Martinsson T: Identification of epigenetically regulated genes that predict patient outcome in neuroblastoma. BMC Cancer. 2011, 11: 66-10.1186/1471-2407-11-66.PubMedPubMed CentralView ArticleGoogle Scholar
- Decock A, Ongenaert M, Vandesompele J, Speleman F: Neuroblastoma epigenetics: from candidate gene approaches to genome-wide screenings. Epigenetics. 2011, 6: 962-970. 10.4161/epi.6.8.16516.PubMedView ArticleGoogle Scholar
- Hoque MO, Kim MS, Ostrow KL, Liu J, Wisman GB a, Park HL, Poeta ML, Jeronimo C, Henrique R, Lendvai A, Schuuring E, Begum S, Rosenbaum E, Ongenaert M, Yamashita K, Califano J, Westra W, van der Zee AGJ, Van Criekinge W, Sidransky D: Genome-wide promoter analysis uncovers portions of the cancer methylome. Cancer Res. 2008, 68: 2661-2670. 10.1158/0008-5472.CAN-07-5913.PubMedPubMed CentralView ArticleGoogle Scholar
- Ongenaert M, Wisman GB a, Volders HH, Koning AJ, Zee AGVD, van Criekinge W, Schuuring E: Discovery of DNA methylation markers in cervical cancer using relaxation ranking. BMC Med Genomics. 2008, 1: 57-10.1186/1755-8794-1-57.PubMedPubMed CentralView ArticleGoogle Scholar
- Huss M: Introduction into the analysis of high-throughput-sequencing based epigenome data. Brief Bioinform. 2010, 11: 512-523. 10.1093/bib/bbq014.PubMedView ArticleGoogle Scholar
- Abe M, Ohira M, Kaneda A, Yagi Y, Yamamoto S, Kitano Y, Takato T, Nakagawara A, Ushijima T: CpG island methylator phenotype is a strong determinant of poor prognosis in neuroblastomas. Cancer Res. 2005, 65: 828-834.PubMedGoogle Scholar
- Banelli B, Brigati C, Di Vinci A, Casciano I, Forlani A, Borzì L, Allemanni G, Romani M: A pyrosequencing assay for the quantitative methylation analysis of the PCDHB gene cluster, the major factor in neuroblastoma methylator phenotype. Lab Invest. 2011, 92: 458-465.PubMedView ArticleGoogle Scholar
- Anders S, Huber W: Differential expression analysis for sequence count data. Genome Biol. 2010, 11: R106-10.1186/gb-2010-11-10-r106.PubMedPubMed CentralView ArticleGoogle Scholar
- Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP: Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005, 102: 15545-15550. 10.1073/pnas.0506580102.PubMedPubMed CentralView ArticleGoogle Scholar
- Cole K a, Huggins J, Laquaglia M, Hulderman CE, Russell MR, Bosse K, Diskin SJ, Attiyeh EF, Sennett R, Norris G, Laudenslager M, Wood AC, Mayes P a, Jagannathan J, Winter C, Mosse YP, Maris JM: RNAi screen of the protein kinome identifies checkpoint kinase 1 (CHK1) as a therapeutic target in neuroblastoma. Proc Natl Acad Sci USA. 2011, 108: 3336-3341. 10.1073/pnas.1012351108.PubMedPubMed CentralView ArticleGoogle Scholar
- Ohtaki M, Otani K, Hiyama K, Kamei N, Satoh K, Hiyama E: A robust method for estimating gene expression states using Affymetrix microarray probe level data. BMC Bioinformatics. 2010, 11: 183-10.1186/1471-2105-11-183.PubMedPubMed CentralView ArticleGoogle Scholar
- Bourdeaut F, Janoueix-Lerosey I, Lucchesi C, Paris R, Ribeiro A, de Pontual L, Amiel J, Lyonnet S, Pierron G, Michon J, Peuchmaur M, Delattre O: Cholinergic switch associated with morphological differentiation in neuroblastoma. J Pathol. 2009, 219: 463-472. 10.1002/path.2614.PubMedView ArticleGoogle Scholar
- Janoueix-Lerosey I, Lequin D, Brugières L, Ribeiro A, de Pontual L, Combaret V, Raynal V, Puisieux A, Schleiermacher G, Pierron G, Valteau-Couanet D, Frebourg T, Michon J, Lyonnet S, Amiel J, Delattre O: Somatic and germline activating mutations of the ALK kinase receptor in neuroblastoma. Nature. 2008, 455: 967-970. 10.1038/nature07398.PubMedView ArticleGoogle Scholar
- Łastowska M, Viprey V, Santibanez-Koref M, Wappler I, Peters H, Cullinane C, Roberts P, Hall AG, Tweddle DA, Pearson ADJ, Lewis I, Burchill SA, Jackson MS: Identification of candidate genes involved in neuroblastoma progression by combining genomic and expression microarrays with survival data. Oncogene. 2007, 26: 7432-7444. 10.1038/sj.onc.1210552.PubMedView ArticleGoogle Scholar
- Wang Q, Diskin S, Rappaport E, Attiyeh E, Mosse Y, Shue D, Seiser E, Jagannathan J, Shusterman S, Bansal M, Khazi D, Winter C, Okawa E, Grant G, Cnaan A, Zhao H, Cheung N-K, Gerald W, London W, Matthay KK, Brodeur GM, Maris JM: Integrative genomics identifies distinct molecular classes of neuroblastoma and shows that multiple genes are targeted by regional alterations in DNA copy number. Cancer Res. 2006, 66: 6050-6062. 10.1158/0008-5472.CAN-05-4618.PubMedView ArticleGoogle Scholar
- Balamuth NJ, Wood A, Wang Q, Jagannathan J, Mayes P, Zhang Z, Chen Z, Rappaport E, Courtright J, Pawel B, Weber B, Wooster R, Sekyere EO, Marshall GM, Maris JM: Serial transcriptome analysis and cross-species integration identifies centromere-associated protein E as a novel neuroblastoma target. Cancer Res. 2010, 70: 2749-2758. 10.1158/0008-5472.CAN-09-3844.PubMedPubMed CentralView ArticleGoogle Scholar
- Wang K, Diskin SJ, Zhang H, Attiyeh EF, Winter C, Hou C, Schnepp RW, Diamond M, Bosse K, Mayes PA, Glessner J, Kim C, Frackelton E, Garris M, Wang Q, Glaberson W, Chiavacci R, Nguyen L, Jagannathan J, Saeki N, Sasaki H, Grant SFA, Iolascon A, Mosse YP, Cole KA, Li H, Devoto M, McGrady PW, London WB, Capasso M, et al: Integrative genomics identifies LMO1 as a neuroblastoma oncogene. Nature. 2011, 469: 216-220. 10.1038/nature09609.PubMedPubMed CentralView ArticleGoogle Scholar
- Schmidt ML, Lal A, Seeger RC, Maris JM, Shimada H, O'Leary M, Gerbing RB, Matthay KK: Favorable prognosis for patients 12 to 18 months of age with stage 4 nonamplified MYCN neuroblastoma: a Children's Cancer Group Study. J Clin Oncol. 2005, 23: 6474-6480. 10.1200/JCO.2005.05.183.PubMedView ArticleGoogle Scholar
- London WB, Castleberry RP, Matthay KK, Look AT, Seeger RC, Shimada H, Thorner P, Brodeur G, Maris JM, Reynolds CP, Cohn SL: Evidence for an age cutoff greater than 365 days for neuroblastoma risk group stratification in the Children's Oncology Group. J Clin Oncol. 2005, 23: 6459-6465. 10.1200/JCO.2005.05.571.PubMedView ArticleGoogle Scholar
- Yagyu S, Gotoh T, Iehara T, Miyachi M, Katsumi Y, Tsubai-Shimizu S, Kikuchi K, Tamura S, Tsuchiya K, Imamura T, Misawa-Furihata A, Sugimoto T, Sawada T, Hosoi H: Circulating methylated-DCR2 gene in serum as an indicator of prognosis and therapeutic efficacy in patients with MYCN nonamplified neuroblastoma. Clin Cancer Res. 2008, 14: 7011-7019. 10.1158/1078-0432.CCR-08-1249.PubMedView ArticleGoogle Scholar
- Yang Q, Kiernan CM, Tian Y, Salwen HR, Chlenski A, Brumback BA, London WB, Cohn SL: Methylation of CASP8, DCR2, and HIN-1 in neuroblastoma is associated with poor outcome. Clin Cancer Res. 2007, 13: 3191-3197. 10.1158/1078-0432.CCR-06-2846.PubMedView ArticleGoogle Scholar
- Banelli B, Bonassi S, Casciano I, Mazzocco K, Di Vinci A, Scaruffi P, Brigati C, Allemanni G, Borza L, Tonini GP, Romani M: Outcome prediction and risk assessment by quantitative pyrosequencing methylation analysis of the SFN gene in advanced stage, high-risk, neuroblastic tumor patients. Int J Cancer. 2010, 126: 656-668. 10.1002/ijc.24768.PubMedView ArticleGoogle Scholar
- Kamimatsuse A, Matsuura K, Moriya S, Fukuba I, Yamaoka H, Fukuda E, Kamei N, Hiyama K, Sueda T, Hiyama E: Detection of CpG island hypermethylation of caspase-8 in neuroblastoma using an oligonucleotide array. Pediatr Blood Cancer. 2009, 52: 777-783. 10.1002/pbc.21977.PubMedView ArticleGoogle Scholar
- Yang Q, Zage P, Kagan D, Tian Y, Seshadri R, Salwen HR, Liu S, Chlenski A, Cohn SL: Association of epigenetic inactivation of RASSF1A with poor outcome in human neuroblastoma. Clin Cancer Res. 2004, 10: 8493-8500. 10.1158/1078-0432.CCR-04-1331.PubMedView ArticleGoogle Scholar
- Alaminos M, Davalos V, Cheung N-KV, Gerald WL, Esteller M: Clustering of gene hypermethylation associated with clinical risk groups in neuroblastoma. J Natl Cancer Inst. 2004, 96: 1208-1219. 10.1093/jnci/djh224.PubMedView ArticleGoogle Scholar
- Grau E, Martinez F, Orellana C, Canete A, Yañez Y, Oltra S, Noguera R, Hernandez M, Bermudez JD, Castel V: Epigenetic alterations in disseminated neuroblastoma tumour cells: influence of TMS1 gene hypermethylation in relapse risk in NB patients. J Cancer Res Clin Oncol. 2010, 136: 1415-1421. 10.1007/s00432-010-0796-9.PubMedView ArticleGoogle Scholar
- Petak I, Danam RP, Tillman DM, Vernes R, Howell SR, Berczi L, Kopper L, Brent TP, Houghton JA: Hypermethylation of the gene promoter and enhancer region can regulate Fas expression and sensitivity in colon carcinoma. Cell Death Differ. 2003, 10: 211-217. 10.1038/sj.cdd.4401132.PubMedView ArticleGoogle Scholar
- Chaopatchayakul P, Jearanaikoon P, Yuenyao P, Limpaiboon T: Aberrant DNA methylation of apoptotic signaling genes in patients responsive and nonresponsive to therapy for cervical carcinoma. Am J Obstet Gynecol. 2010, 202: 281.e1-9. 10.1016/j.ajog.2009.11.037.View ArticleGoogle Scholar
- Schaefer JT, Barthlen W, Schweizer P: Ceramide induces apoptosis in neuroblastoma cell cultures resistant to CD95 (Fas/APO-1)-mediated apoptosis. J Pediatr Surg. 2000, 35: 473-479. 10.1016/S0022-3468(00)90217-3.PubMedView ArticleGoogle Scholar
- Kisenge RR, Toyoda H, Kang J, Tanaka S, Yamamoto H, Azuma E, Komada Y: Expression of short-form caspase 8 correlates with decreased sensitivity to Fas-mediated apoptosis in neuroblastoma cells. Cancer Sci. 2003, 94: 598-605. 10.1111/j.1349-7006.2003.tb01489.x.PubMedView ArticleGoogle Scholar
- Takamizawa S, Okamoto S, Wen J, Bishop W, Kimura K, Sandler A: Overexpression of Fas-ligand by neuroblastoma: a novel mechanism of tumor-cell killing. J Pediatr Surg. 2000, 35: 375-379. 10.1016/S0022-3468(00)90044-7.PubMedView ArticleGoogle Scholar
- Bernassola F, Scheuerpflug C, Herr I, Krammer PH, Debatin KM, Melino G: Induction of apoptosis by IFNgamma in human neuroblastoma cell lines through the CD95/CD95L autocrine circuit. Cell Death Differ. 1999, 6: 652-660. 10.1038/sj.cdd.4400537.PubMedView ArticleGoogle Scholar
- Gross N, Balmas K, Beretta Brognara C, Tschopp J: Expression of Fas (APO-1/CD95) and Fas ligand (FasL) in human neuroblastoma. Med Pediatr Oncol. 2001, 36: 111-114. 10.1002/1096-911X(20010101)36:1<111::AID-MPO1026>3.0.CO;2-3.PubMedView ArticleGoogle Scholar
- Pedersen WA, Becker LE, Yeger H: Expression and distribution of peripherin protein in human neuroblastoma cell lines. Int J Cancer. 1993, 53: 463-470. 10.1002/ijc.2910530319.PubMedView ArticleGoogle Scholar
- Ambros IM, Hata J, Joshi VV, Roald B, Dehner LP, Tüchler H, Pötschger U, Shimada H: Morphologic features of neuroblastoma (Schwannian stroma-poor tumors) in clinically favorable and unfavorable groups. Cancer. 2002, 94: 1574-1583. 10.1002/cncr.10359.PubMedView ArticleGoogle Scholar
- 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, Ebus ME, Haneveld F, Lakeman A, Schild L, Molenaar P, Stroeken P, van Noesel MM, Ora I, Santo EE, Caron HN, Westerhout EM, Versteeg R: Sequencing of neuroblastoma identifies chromothripsis and defects in neuritogenesis genes. Nature. 2012, 483: 589-593. 10.1038/nature10910.PubMedView ArticleGoogle Scholar
- Chen Y, Stallings RL: Differential patterns of microRNA expression in neuroblastoma are correlated with prognosis, differentiation, and apoptosis. Cancer Res. 2007, 67: 976-983. 10.1158/0008-5472.CAN-06-3667.PubMedView ArticleGoogle Scholar
- Mestdagh P, Van Vlierberghe P, De Weer A, Muth D, Westermann F, Speleman F, Vandesompele J: A novel and universal method for microRNA RT-qPCR data normalization. Genome Biol. 2009, 10: R64-10.1186/gb-2009-10-6-r64.PubMedPubMed CentralView ArticleGoogle Scholar
- Schulte JH, Horn S, Otto T, Samans B, Heukamp LC, Eilers U-C, Krause M, Astrahantseff K, Klein-Hitpass L, Buettner R, Schramm A, Christiansen H, Eilers M, Eggert A, Berwanger B: MYCN regulates oncogenic MicroRNAs in neuroblastoma. Int J Cancer. 2008, 122: 699-704. 10.1002/ijc.23153.PubMedView ArticleGoogle Scholar
- Ragusa M, Majorana A, Banelli B, Barbagallo D, Statello L, Casciano I, Guglielmino MR, Duro LR, Scalia M, Magro G, Di Pietro C, Romani M, Purrello M: MIR152, MIR200B, and MIR338, human positional and functional neuroblastoma candidates, are involved in neuroblast differentiation and apoptosis. J Mol Med (Berl). 2010, 88: 1041-1053. 10.1007/s00109-010-0643-0.View ArticleGoogle Scholar
- Tan Y-C, Blumenfeld J, Rennert H: Autosomal dominant polycystic kidney disease: genetics, mutations and microRNAs. Biochim Biophys Acta. 2011, 1812: 1202-1212. 10.1016/j.bbadis.2011.03.002.PubMedView ArticleGoogle Scholar
- Hansen TB, Wiklund ED, Bramsen JB, Villadsen SB, Statham AL, Clark SJ, Kjems J: miRNA-dependent gene silencing involving Ago2-mediated cleavage of a circular antisense RNA. EMBO J. 2011, 30: 4414-4422. 10.1038/emboj.2011.359.PubMedPubMed CentralView ArticleGoogle Scholar
- Teitz T, Lahti JM, Kidd VJ: Aggressive childhood neuroblastomas do not express caspase-8: an important component of programmed cell death. J Mol Med (Berl). 2001, 79: 428-436. 10.1007/s001090100233.View ArticleGoogle Scholar
- Hoebeeck J, Michels E, Pattyn F, Combaret V, Vermeulen J, Yigit N, Hoyoux C, Laureys G, De Paepe A, Speleman F, Vandesompele J: Aberrant methylation of candidate tumor suppressor genes in neuroblastoma. Cancer Lett. 2009, 273: 336-346. 10.1016/j.canlet.2008.08.019.PubMedView ArticleGoogle Scholar
- Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, Hornik K, Hothorn T, Huber W, Iacus S, Irizarry R, Leisch F, Li C, Maechler M, Rossini AJ, Sawitzki G, Smith C, Smyth G, Tierney L, Yang JYH, Zhang J: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004, 5: R80-10.1186/gb-2004-5-10-r80.PubMedPubMed CentralView ArticleGoogle Scholar
- Morgan M, Anders S, Lawrence M, Aboyoun P, Pagès H, Gentleman R: ShortRead: a bioconductor package for input, quality assessment and exploration of high-throughput sequence data. Bioinformatics. 2009, 25: 2607-2608. 10.1093/bioinformatics/btp450.PubMedPubMed CentralView ArticleGoogle Scholar
- Lawrence M, Gentleman R, Carey V: rtracklayer: an R package for interfacing with genome browsers. Bioinformatics. 2009, 25: 1841-1842. 10.1093/bioinformatics/btp328.PubMedPubMed CentralView ArticleGoogle Scholar
- Hellemans J, Mortier G, De Paepe A, Speleman F, Vandesompele J: qBase relative quantification framework and software for management and automated analysis of real-time quantitative PCR data. Genome Biol. 2007, 8: R19-10.1186/gb-2007-8-2-r19.PubMedPubMed CentralView ArticleGoogle Scholar
- Lefever S, Vandesompele J, Speleman F, Pattyn F: RTPrimerDB: the portal for real-time PCR primers and probes. Nucleic Acids Res. 2009, 37: D942-945. 10.1093/nar/gkn777.PubMedPubMed CentralView ArticleGoogle Scholar
- Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, Speleman F: Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002, 3: RESEARCH0034Google Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.