Genome-wide promoter methylation analysis in neuroblastoma identifies prognostic methylation biomarkers.

BACKGROUND: 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. RESULTS: 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. CONCLUSIONS: 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.

Meta-analysis of neuroblastomas reveals a skewed ALK mutation spectrum in tumors with MYCN amplification.

PURPOSE: Activating mutations of the anaplastic lymphoma kinase (ALK) were recently described in neuroblastoma. We carried out a meta-analysis of 709 neuroblastoma tumors to determine their frequency and mutation spectrum in relation to genomic and clinical parameters, and studied the prognostic significance of ALK copy number and expression. EXPERIMENTAL DESIGN: The frequency and type of ALK mutations, copy number gain, and expression were analyzed in a new series of 254 neuroblastoma tumors. Data from 455 published cases were used for further in-depth analysis. RESULTS: ALK mutations were present in 6.9% of 709 investigated tumors, and mutations were found in similar frequencies in favorable [International Neuroblastoma Staging System (INSS) 1, 2, and 4S; 5.7%] and unfavorable (INSS 3 and 4; 7.5%) neuroblastomas (P = 0.087). Two hotspot mutations, at positions R1275 and F1174, were observed (49% and 34.7% of the mutated cases, respectively). Interestingly, the F1174 mutations occurred in a high proportion of MYCN-amplified cases (P = 0.001), and this combined occurrence was associated with a particular poor outcome, suggesting a positive cooperative effect between both aberrations. Furthermore, the F1174L mutant was characterized by a higher degree of autophosphorylation and a more potent transforming capacity as compared with the R1275Q mutant. Chromosome 2p gains, including the ALK locus (91.8%), were associated with a significantly increased ALK expression, which was also correlated with poor survival. CONCLUSIONS: ALK mutations occur in equal frequencies across all genomic subtypes, but F1174L mutants are observed in a higher frequency of MYCN-amplified tumors and show increased transforming capacity as compared with the R1275Q mutants.

methGraph: a genome visualization tool for PCR-based methylation assays.

Abnormalities in DNA methylation of CpG islands that play a role in gene regulation affect gene expression and hence play a role in disease, including cancer. Bisulfite-based DNA methylation analysis methods such as methylation-specific PCR (MSP) and bisulfite sequencing (BiSeq) are most commonly used to study gene-specific DNA methylation. Assessing specificity and visualizing the position of PCR primers in their genomic context is a laborious and tedious task, primarily due to the sequence changes induced during the bisulfite conversion. For this purpose, we developed methGraph, a web application for easy, fast and flexible visualization and accurate in silico quality evaluation of PCR-based methylation assays. The visualization process starts by submitting PCR primer sequences for specificity assessment and mapping on the genome using the BiSearch ePCR primer-search algorithm. The next step comprises the selection of relevant UCSC genome annotation tracks for display in the final graph. A custom track showing all individual CpG dinucleotides, representing their distribution in the CpG island is also provided. Finally, methGraph creates a BED file that is automatically uploaded to the UCSC genome browser, after which the resulting image files are extracted and made available for visualization and download. The generated high-quality figures can easily be customized and exported for use in publications or presentations. methGraph is available at http://mellfire.ugent.be/methgraph/.

The emerging molecular pathogenesis of neuroblastoma: implications for improved risk assessment and targeted therapy.

Neuroblastoma is one of the most common solid tumors of childhood, arising from immature sympathetic nervous system cells. The clinical course of patients with neuroblastoma is highly variable, ranging from spontaneous regression to widespread metastatic disease. Although the outcome for children with cancer has improved considerably during the past decades, the prognosis of children with aggressive neuroblastoma remains dismal. The clinical heterogeneity of neuroblastoma mirrors the biological and genetic heterogeneity of these tumors. Ploidy and MYCN amplification have been used as genetic markers for risk stratification and therapeutic decision making, and, more recently, gene expression profiling and genome-wide DNA copy number analysis have come into the picture as sensitive and specific tools for assessing prognosis. The applica tion of new genetic tools also led to the discovery of an important familial neuroblastoma cancer gene, ALK, which is mutated in approximately 8% of sporadic tumors, and genome-wide association studies have unveiled loci with risk alleles for neuroblastoma development. For some of the genomic regions that are deleted in some neuroblastomas, on 1p, 3p and 11q, candidate tumor suppressor genes have been identified. In addition, evidence has emerged for the contribution of epigenetic disturbances in neuroblastoma oncogenesis. As in other cancer entities, altered microRNA expression is also being recognized as an important player in neuroblastoma. The recent successes in unraveling the genetic basis of neuroblastoma are now opening opportunities for development of targeted therapies.

Chromosome 3p microsatellite allelotyping in neuroblastoma: a report on the technical hurdles.

Pinpointing critical regions of recurrent loss may help localize tumor suppressor genes. To determine the regions of loss on chromosome 3p in neuroblastoma, we performed loss of heterozygosity analysis using 16 microsatellite markers in a series of 65 primary tumors and 29 neuroblastoma cell lines. In this study, we report the results and discuss the technical hurdles that we encountered during data generation and interpretation that are of relevance for current studies or tests employing microsatellites. To provide functional support for the implication of 3p tumor suppressor genes in this childhood malignancy, we performed a microcell-mediated chromosome 3 transfer in neuroblastoma cells.

Aberrant methylation of candidate tumor suppressor genes in neuroblastoma.

CpG island hypermethylation has been recognized as an alternative mechanism for tumor suppressor gene inactivation. In this study, we performed methylation-specific PCR (MSP) to investigate the methylation status of 10 selected tumor suppressor genes in neuroblastoma. Seven of the investigated genes (CD44, RASSF1A, CASP8, PTEN, ZMYND10, CDH1, PRDM2) showed high frequencies (> or =30%) of methylation in 33 neuroblastoma cell lines. In 42 primary neuroblastoma tumors, the frequencies of methylation were 69%, CD44; 71%, RASSF1A; 56%, CASP8; 25%, PTEN; 15%, ZMYND10; 8%, CDH1; and 0%, PRDM2. Furthermore, CASP8 and CDH1 hypermethylation was significantly associated with poor event-free survival. Meta-analysis of 115 neuroblastoma tumors demonstrated a significant correlation between CASP8 methylation and MYCN amplification. In addition, there was a correlation between ZMYND10 methylation and MYCN amplification. The MSP data, together with optimized mRNA re-expression experiments (in terms of concentration and time of treatment and use of proper reference genes) further strengthen the notion that epigenetic alterations could play a significant role in NB oncogenesis. This study thus warrants the need for a global profiling of gene promoter hypermethylation to identify genome-wide aberrantly methylated genes in order to further understand neuroblastoma pathogenesis and to identify prognostic methylation markers.

CADM1 is a strong neuroblastoma candidate gene that maps within a 3.72 Mb critical region of loss on 11q23.

BACKGROUND: Recurrent loss of part of the long arm of chromosome 11 is a well established hallmark of a subtype of aggressive neuroblastomas. Despite intensive mapping efforts to localize the culprit 11q tumour suppressor gene, this search has been unsuccessful thus far as no sufficiently small critical region could be delineated for selection of candidate genes. METHODS: To refine the critical region of 11q loss, the chromosome 11 status of 100 primary neuroblastoma tumours and 29 cell lines was analyzed using a BAC array containing a chromosome 11 tiling path. For the genes mapping within our refined region of loss, meta-analysis on published neuroblastoma mRNA gene expression datasets was performed for candidate gene selection. The DNA methylation status of the resulting candidate gene was determined using re-expression experiments by treatment of neuroblastoma cells with the demethylating agent 5-aza-2'-deoxycytidine and bisulphite sequencing. RESULTS: Two small critical regions of loss within 11q23 at chromosomal band 11q23.1-q23.2 (1.79 Mb) and 11q23.2-q23.3 (3.72 Mb) were identified. In a first step towards further selection of candidate neuroblastoma tumour suppressor genes, we performed a meta-analysis on published expression profiles of 692 neuroblastoma tumours. Integration of the resulting candidate gene list with expression data of neuroblastoma progenitor cells pinpointed CADM1 as a compelling candidate gene. Meta-analysis indicated that CADM1 expression has prognostic significance and differential expression for the gene was noted in unfavourable neuroblastoma versus normal neuroblasts. Methylation analysis provided no evidence for a two-hit mechanism in 11q deleted cell lines. CONCLUSION: Our study puts CADM1 forward as a strong candidate neuroblastoma suppressor gene. Further functional studies are warranted to elucidate the role of CADM1 in neuroblastoma development and to investigate the possibility of CADM1 haploinsufficiency in neuroblastoma.

ArrayCGH-based classification of neuroblastoma into genomic subgroups.

High-resolution array comparative genomic hybridization (arrayCGH) profiling was performed on 75 primary tumors and 29 cell lines to gain further insight into the genetic heterogeneity of neuroblastoma and to refine genomic subclassification. Using a novel data-mining strategy, three major and two minor genomic subclasses were delineated. Eighty-three percent of tumors could be assigned to the three major genomic subclasses, corresponding to the three known clinically and biologically relevant subsets in neuroblastoma. The remaining subclasses represented (1) tumors with no/few copy number alterations or an atypical pattern of aberrations and (2) tumors with 11q13 amplification. Inspection of individual arrayCGH profiles showed that recurrent genomic imbalances were not exclusively associated with a specific subclass. Of particular notice were tumors with numerical imbalances typically observed in subtype 1 neuroblastoma, in association with genomic features of subtype 2A or 2B. A search for prognostically relevant genomic alterations disclosed 1q gain as a predictive marker for therapy failure within the group of subtype 2A and 2B tumors. In cell lines, a high incidence of 6q loss was observed, with a 3.87-5.32 Mb region of common loss within 6q25.1-6q25.2. Our study clearly illustrates the importance of genomic profiling in relation to tumor behavior in neuroblastoma. We propose that genome-wide assessment of copy number alterations should ideally be included in the genetic workup of neuroblastoma. Further multicentric studies on large tumor series are warranted in order to improve therapeutic stratification in conjunction with other features such as age at diagnosis, tumor stage, and gene expression signatures.

Real-time quantitative PCR as an alternative to Southern blot or fluorescence in situ hybridization for detection of gene copy number changes.

Changes in copy number of genes contribute to the pathogenesis of various genetic disorders and cancer. The status of a gene has not only diagnostic value but sometimes directs treatment stratification. Although, for many years, Southern blot and fluorescence in situ hybridization were the standard methods for the detection of deletion, duplication, or amplification of a gene, both methods have their own important limitations. Recently, realtime quantitative PCR has proven to be a good alternative for the detection of gene copy number changes. Its main advantages are the large dynamic range of accurate quantification, the absence of post-PCR manipulations, its high-throughput screening capacity and degree of automation, and the possibility to perform the assay on minimal amounts of sample DNA in just a few hours of time. In this chapter, we outline the procedure of how to develop an assay for the detection of gene copy number changes for your gene of interest. We illustrate the approach by describing a validated assay for the detection of germline VHL exon deletions and for determination of MYCN copy numbers in tumor samples.

High resolution tiling-path BAC array deletion mapping suggests commonly involved 3p21-p22 tumor suppressor genes in neuroblastoma and more frequent tumors.

The recurrent loss of 3p segments in neuroblastoma suggests the implication of 1 or more tumor suppressor genes but thus far few efforts have been made to pinpoint their detailed chromosomal position. To achieve this goal, array-based comparative genomic hybridization was performed on a panel of 23 neuroblastoma cell lines and 75 primary tumors using a tiling-path bacterial artificial chromosome array for chromosome 3p. A total of 45 chromosome 3 losses were detected, including whole chromosome losses, large terminal deletions and interstitial deletions. The latter, observed in cell lines as well as a number of distal deletions detected in primary tumors, allowed us to demarcate 3 minimal regions of loss of 3.6 Mb [3p21.31-p21.2, shortest regions of overlap (SRO)1], 1.4 Mb (3p22.3-3p22.2, SRO2) and 3.8 Mb (3p25.3-p25.1, SRO3) in size. The present data significantly extend previous findings and now firmly establish critical regions on 3p implicated in neuroblastoma. Interestingly, the 2 proximal regions coincide with previously defined SROs on 3p21.3 in more frequent tumors including lung and breast cancer. As such, similar tumor suppressor genes may play a critical role in development or progression of a variety of neoplasms, including neuroblastoma.

methBLAST and methPrimerDB: web-tools for PCR based methylation analysis.

BACKGROUND: DNA methylation plays an important role in development and tumorigenesis by epigenetic modification and silencing of critical genes. The development of PCR-based methylation assays on bisulphite modified DNA heralded a breakthrough in speed and sensitivity for gene methylation analysis. Despite this technological advancement, these approaches require a cumbersome gene by gene primer design and experimental validation. Bisulphite DNA modification results in sequence alterations (all unmethylated cytosines are converted into uracils) and a general sequence complexity reduction as cytosines become underrepresented. Consequently, standard BLAST sequence homology searches cannot be applied to search for specific methylation primers. RESULTS: To address this problem we developed methBLAST, a sequence similarity search program, based on the original BLAST algorithm but querying in silico bisulphite modified genome sequences to evaluate oligonucleotide sequence similarities. Apart from the primer specificity analysis tool, we have also developed a public database termed methPrimerDB for the storage and retrieval of validated PCR based methylation assays. The web interface allows free public access to perform methBLAST searches or database queries and to submit user based information. Database records can be searched by gene symbol, nucleotide sequence, analytical method used, Entrez Gene or methPrimerDB identifier, and submitter's name. Each record contains a link to Entrez Gene and PubMed to retrieve additional information on the gene, its genomic context and the article in which the methylation assay was described. To assure and maintain data integrity and accuracy, the database is linked to other reference databases. Currently, the database contains primer records for the most popular PCR-based methylation analysis methods to study human, mouse and rat epigenetic modifications. methPrimerDB and methBLAST are available at http://medgen.ugent.be/methprimerdb and http://medgen.ugent.be/methblast. CONCLUSION: We have developed two integrated and freely available web-tools for PCR based methylation analysis. methBLAST allows in silico assessment of primer specificity in PCR based methylation assays that can be stored in the methPrimerDB database, which provides a search portal for validated methylation assays.

The von Hippel-Lindau tumor suppressor gene expression level has prognostic value in neuroblastoma.

Deletions of the short arm of chromosome 3 are often observed in a specific subset of aggressive neuroblastomas (NBs) with loss of distal 11q and without MYCN amplification. The critical deleted region encompasses the locus of the von Hippel-Lindau gene (VHL, 3p25). Constitutional loss of function mutations in the VHL gene are responsible for the VHL syndrome, a dominantly inherited familial cancer syndrome predisposing to a variety of neoplasms, including pheochromocytoma. Pheochromocytomas are, like NB, derived from neural crest cells, but, unlike NB, consist of more mature chromaffin cells instead of immature neuroblasts. Further arguments for a putative role of VHL in NB are its function as oxygen sensitizer and the reported relation between hypoxia and dedifferentiation of NB cells, leading to a more aggressive phenotype. To test the possible involvement of VHL in NB, we did mRNA expression analysis and sought evidence for VHL gene inactivation. Although no evidence for a classic tumor suppressor role for VHL in NB could be obtained, a strong correlation was observed between reduced levels of VHL mRNA and low patient survival probability (p=0.013). Furthermore, VHL appears to have predictive power in NTRK1 (TRKA) positive tumor samples with presumed favorable prognosis, which makes it a potentially valuable marker for more accurate risk assessment in this subgroup of patients. The significance of the reduced VHL expression levels in relation to NB tumor biology remains unexplained, as functional analysis demonstrated no clear effect of the reduction in VHL mRNA expression on protein stability of its downstream target hypoxia-inducible factor alpha.

Rapid detection of VHL exon deletions using real-time quantitative PCR.

Various types of mutations exist that exert an effect on the normal function of a gene. Among these, exon/gene deletions often remain unnoticed in initial mutation screening. Until recently, no fast and efficient methods were available to detect this type of mutation. Molecular detection methods for gene copy number changes included Southern blot (SB) and fluorescence in situ hybridisation, both with their own intrinsic limitations. In this paper, we report the development and application of a fast, sensitive and high-resolution method for the detection of single exon or larger deletions in the VHL gene based on real-time quantitative PCR (Q-PCR). These deletions account for approximately one-fifth of all patients with the von Hippel-Lindau syndrome, a dominantly inherited highly penetrant familial cancer syndrome predisposing to specific malignancies including phaeochromocytomas and haemangioblastomas. Our VHL exon quantification strategy is based on SYBR Green I detection and normalisation using two reference genes with a normal copy number, that is, ZNF80 (3q13.31) and GPR15 (3q12.1). Choice of primer sequences and the use of two reference genes appears to be critical for accurate discrimination between 1 and 2 exon copies. In a blind Q-PCR study of 29 samples, all 14 deletions were detected, which is in perfect agreement with previously determined SB results. We propose Q-PCR as the method of choice for fast (within 3.5 h), accurate and sensitive (ng amount of input DNA) exon deletion screening in routine DNA diagnosis of VHL disease. Similar assays can be designed for deletion screening in other genetic disorders.

No evidence for involvement of SDHD in neuroblastoma pathogenesis.

BACKGROUND: Deletions in the long arm of chromosome 11 are observed in a subgroup of advanced stage neuroblastomas with poor outcome. The deleted region harbours the tumour suppressor gene SDHD that is frequently mutated in paraganglioma and pheochromocytoma, which are, like neuroblastoma, tumours originating from the neural crest. In this study, we sought for evidence for involvement of SDHD in neuroblastoma. METHODS: SDHD was investigated on the genome, transcriptome and proteome level using mutation screening, methylation specific PCR, real-time quantitative PCR based homozygous deletion screening and mRNA expression profiling, immunoblotting, functional protein analysis and ultrastructural imaging of the mitochondria. RESULTS: Analysis at the genomic level of 67 tumour samples and 37 cell lines revealed at least 2 bona-fide mutations in cell lines without allelic loss at 11q23: a 4bp-deletion causing skip of exon 3 resulting in a premature stop codon in cell line N206, and a Y93C mutation in cell line NMB located in a region affected by germline SDHD mutations causing hereditary paraganglioma. No evidence for hypermethylation of the SDHD promotor region was observed, nor could we detect homozygous deletions. Interestingly, SDHD mRNA expression was significantly reduced in SDHD mutated cell lines and cell lines with 11q allelic loss as compared to both cell lines without 11q allelic loss and normal foetal neuroblast cells. However, protein analyses and assessment of mitochondrial morphology presently do not provide clues as to the possible effect of reduced SDHD expression on the neuroblastoma tumour phenotype. CONCLUSIONS: Our study provides no indications for 2-hit involvement of SDHD in the pathogenesis of neuroblastoma. Also, although a haplo-insufficient mechanism for SDHD involvement in advanced stage neuroblastoma could be considered, the present data do not provide consistent evidence for this hypothesis.

Application of laser capture microdissection in genetic analysis of neuroblastoma and neuroblastoma precursor cells.

Recently developed quantitative and high-throughput technologies that allow automated and rapid screening of the whole genome, transcriptome and proteome have revolutionized the field of cancer genetics. At the same time, new challenges are met, e.g. the need for improved data analysis and standardization of tumor sample handling. Even if these issues are resolved, an 'old' problem in genetic tumor analysis remains, i.e. contamination of tumor samples by stromal and surrounding normal cells. To overcome this obstacle, laser capture microdissection (LCM) has been developed in order to procure the cells of interest from stained tissue sections with retention of morphology. In this review we describe the possible down-stream applications of LCM in the genetic analysis of neuroblastoma (NB). Special focus is given to MYCN copy number determination using real-time quantitative polymerase chain reaction (Q-PCR), analysis of 1p-, 3p- and 11q-deletions using loss of heterozygosity analysis and Q-PCR expression analysis of microdissected normal neuroblast cells and NB cells.