Neuroblastoma is remarkable for its clinical heterogeneity and is characterized by genomic alterations that are strongly correlated with tumor behavior. The specific genes that influence neuroblastoma biology and are targeted by genomic alterations remain largely unknown. We quantified mRNA expression in a highly annotated series of 101 prospectively collected diagnostic neuroblastoma primary tumors using an oligonucleotide-based microarray. Genomic copy number status at the prognostically relevant loci 1p36, 2p24 (MYCN), 11q23, and 17q23 was determined by PCR and was aberrant in 26, 20, 40, and 38 cases, respectively. In addition, 72 diagnostic neuroblastoma primary tumors assayed in a different laboratory were used as an independent validation set. Unsupervised hierarchical clustering showed that gene expression was highly correlated with genomic alterations and clinical markers of tumor behavior. The vast majority of samples with MYCN amplification and 1p36 loss of heterozygosity (LOH) clustered together on a terminal node of the sample dendrogram, whereas the majority of samples with 11q deletion clustered separately and both of these were largely distinct from the copy number neutral group of tumors. Genes involved in neurodevelopment were broadly overrepresented in the more benign tumors, whereas genes involved in RNA processing and cellular proliferation were highly represented in the most malignant cases. By combining transcriptomic and genomic data, we showed that LOH at 1p and 11q was associated with significantly decreased expression of 122 (61%) and 88 (27%) of the genes mapping to 1p35-36 and all of 11q, respectively, suggesting that multiple genes may be targeted by LOH events. A total of 71 of the 1p35-36 genes were also differentially expressed in the independent validation data set, providing a prioritized list of candidate neuroblastoma suppressor genes. Taken together, these data are consistent with the hypotheses that the neuroblastoma transcriptome is a sensitive marker of underlying tumor biology and that chromosomal deletion events in this cancer likely target multiple genes through alteration in mRNA dosage. Lead positional candidates for neuroblastoma suppressor genes can be inferred from these data, but the potential multiplicity of transcripts involved has significant implications for ongoing gene discovery strategies. (Cancer Res 2006; 66(12): 6050-62)

BackgroundNeuroblastomas are characterized by hemizygous 1p deletions, suggesting that a tumor suppressor gene resides in this region. We previously mapped the smallest region of consistent deletion to a 2-Mb region of 1p36.31 that encodes 23 genes. Based on mutation analysis, expression pattern, and putative function, we identified CHD5 as the best tumor suppressor gene candidate.MethodsWe determined the methylation status of the CHD5 gene promoter in NLF and IMR5 (with 1p deletion) and SK-N-SH and SK-N-FI neuroblastoma cell lines using methylation-specific sequencing and measured CHD5 mRNA expression by reverse transcription polymerase chain reaction in cells treated with or without 5-aza-2-deoxycytidine, an inhibitor of DNA methylation. We transfected the cells with CHD5 and antisense (AS) CHD5 DNA to assess the effect of CHD5 overexpression and suppression, respectively, on colony formation in soft agar and growth of xenograft tumors in athymic mice. We also analyzed the association of CDH5 expression with outcomes of 99 neuroblastoma patients. Statistical tests were two-sided.ResultsCHD5 expression was very low or absent in neuroblastoma cell lines. The CHD5 promoter was highly methylated in NLF and IMR5 lines, and CHD5 expression increased after treatment with 5-aza-2-deoxycytidine. Clonogenicity and tumor growth were abrogated in NLF and IMR5 cells overexpressing CHD5 compared with antisense CHD5 (clonogenicity: mean no. of colonies per plate, NLF-CHD5, 43 colonies, 95% confidence interval [CI] = 35 to 51 colonies, vs NLF-CHD5-AS, 74 colonies, 95% CI = 62 to 86 colonies, P < .001; IMR5-CHD5, 11 colonies, 95% CI = 2 to 20 colonies, vs IMR5-CHD5-AS, 39 colonies, 95% CI = 17 to 60 colonies, P = .01; tumor growth, n = 10 mice per group: mean tumor size at 5 weeks, NLF-CHD5, 0.36 cm3, 95% CI = 0.17 to 0.44 cm3, vs NLF-CHD5-AS, 1.65 cm3, 95% CI = 0.83 to 2.46 cm3, P = .002; IMR5-CHD5, 0.28 cm3, 95% CI = 0.18 to 0.38 cm3, vs IMR5-CHD5-AS, 1.15 cm3, 95% CI = 0.43 to 1.87 cm3; P = .01). High CHD5 expression was strongly associated with favorable event-free and overall survival (P < .001), even after correction for MYCN amplification and 1p deletion (P = .027).ConclusionsCHD5 is the strongest candidate tumor suppressor gene that is deleted from 1p36.31 in neuroblastomas, and inactivation of the second allele may occur by an epigenetic mechanism.

Substantial genomic and functional evidence from primary tumors and cell lines indicates that a consistent region of distal chromosome 1p is deleted in a sizable proportion of human neuroblastomas, suggesting that this region contains one or more tumor suppressor genes. To determine systematically and precisely the location and extent of 1p deletion in neuroblastomas, we performed allelic loss studies of 737 primary neuroblastomas and genotype analysis of 46 neuroblastoma cell lines. Together, the results defined a single region within 1p36.3 that was consistently deleted in 25% of tumors and 87% of cell lines. Two neuroblastoma patients had constitutional deletions of distal 1p36 that overlapped the tumor-defined region. The tumor-and constitutionally-derived deletions together defined a smallest region of consistent deletion (SRD) between D1S2795 and D1S253. The 1p36.3 SRD was deleted in all but one of the 184 tumors with 1p deletion. Physical mapping and DNA sequencing determined that the SRD minimally spans an estimated 729 kb. Genomic content and sequence analysis of the SRD identified 15 characterized, nine uncharacterized, and six predicted genes in the region. The RNA expression profiles of 21 of the genes were investigated in a variety of normal tissues. The SHREW1 and KCNAB2 genes both had tissue-restricted expression patterns, including expression in the nervous system. In addition, a novel gene (CHD5) with strong homology to proteins involved in chromatin remodeling was expressed mainly in neural tissues. Together, these results suggest that one or more genes involved in neuroblastoma tumorigenesis or tumor progression are likely contained within this region.

Neuroblastomas are characterized by 1p deletions, suggesting that a tumor suppressor gene (TSG) resides in this region. We have mapped the smallest region of deletion (SRD) to a 2 Mb region of 1p36.31 using microsatellite and single nucleotide polymorphisms. We have identified 23 genes in this region, and we have analysed these genes for mutations and RNA expression patterns to identify candidate TSGs. We sequenced the coding exons of these genes in 30 neuroblastoma cell lines. Although rare mutations were found in 10 of the 23 genes, none showed a pattern of genetic change consistent with homozygous inactivation. We examined the expression of these 23 genes in 20 neuroblastoma cell lines, and most showed readily detectable expression, and no correlation with 1p deletion. However, 7 genes showed uniformly low expression in the lines, and 2 genes (CHD5, RNF207) had virtually absent expression, consistent with the expected pattern for a TSG. Our mutation and expression analysis in neuroblastoma cell lines, combined with expression analysis in normal tissues, putative function and prior implication in neuroblastoma pathogenesis, suggests that the most promising TSG deleted from the 1p36 SRD is CHD5, but TNFRSF25, CAMTA1 and AJAP1 are also viable candidates.