- Research article
- Open Access
CADM1 is a strong neuroblastoma candidate gene that maps within a 3.72 Mb critical region of loss on 11q23
BMC Cancer volume 8, Article number: 173 (2008)
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.
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.
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.
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.
Neuroblastoma (NB) is a rare but often highly aggressive tumour in children. Despite intensive gene copy number and mRNA expression studies thus far only two genes, namely MYCN and PHOX2B [2–4] have been found to be directly implicated in NB development. In order to provide clues for more effective therapies, insights into the molecular pathogenesis of this tumour are urgently needed. Analyses of recurrent patterns of somatically acquired DNA copy number alterations resulted in the delineation of three major genetic subgroups with predictive tumour behaviour (subtype 1, 2A and 2B), of which subtype 2A NB represents an aggressive subgroup of metastatic NB [5, 6] characterised by loss of 11q, gain of 17q and a normal MYCN copy number status .
Deletion of the long arm of chromosome 11 is found in 15–22% of sporadic NB [5–10] and has also been described in constitutional cases of NB [11, 12], suggesting the presence of one or more tumour suppressor gene(s) on chromosome 11. Functional evidence for a NB tumour suppressor gene at 11q came from microcell mediated chromosome transfer (MMCT) experiments, in which transfer of an intact chromosome 11 in a NB cell line with 11q loss resulted in a more differentiated phenotype . Array comparative genomic hybridisation (array CGH) analysis of these MMCT hybrids revealed 11q25 as a plausible location for a NB differentiation gene . Extensive microsatellite heterozygosity mapping studies however point at various critical regions of loss, located at 11q23.3  and within the chromosomal region 11q14-11q23 . Despite these mapping efforts, no genes with proven tumour suppressor activity in NB have been identified thus far.
In this study, we used an integrated approach combining high resolution copy number profiling and gene expression meta-analysis. This approach identified CADM1 as a strong candidate 11q tumour suppressor gene with prognostic power, which possibly exerts its effect through haplo-insufficiency.
Array CGH copy number profiling of NB patients and cell lines
A description of the primary NB tumour samples and NB cell lines as well as the array CGH procedure is given in Michels et al. . Twenty-five additional tumour cases were profiled for this study, including 4 stage 1, 3 stage 2, 6 stage 3, 9 stage 4 and 3 stage 4S tumours according to the International Neuroblastoma Staging System . The maximal size of the lost region is determined as the distance between the two normal clones flanking the lost clones. Detailed data for 75 tumours are published by Michels et al. . Detailed data for the 25 additionally profiled tumour cases are available in Additional File 1. Data for all tumours are also accessible from the webtool arrayCGHbase [18, 19]. Mapping data are based on Ensembl v44.
Meta-analysis of published NB gene expression datasets
Expression data were collected from seven independent gene expression studies in NB [20–26]. For each of the individual expression studies Cox regression analysis (enter method) was performed for each gene present in the defined SROs using the SPSS 15 software, allowing to pinpoint genes within our SROs having prognostic power. Next, in four studies for which genomic data were available [21, 22, 25, 26], tumours were also classified according to the three major NB subtypes: subtype 1 (characterised by mainly numerical aberrations and not 11q loss or MYCN amplification), 2A (characterised by 11q loss, without MYCN amplification) or 2B (characterised by MYCN amplification). Following this subclassification an independent Kruskal-Wallis analysis was performed for each of the four studies using the R software package with multiple testing correction. Genes were considered to be significant when a p-value < 0.5 was calculated (Additional file 2).
Following these analyses, an intersection was made between genes that were prognostically significant (Cox regression analysis) in more than three independent studies and genes which had subgroup discriminating power in more than two independent studies (Kruskal-Wallis analysis) to identify the strongest candidate genes (see results and discussion).
mRNA expression analysis for CADM1in cell lines untreated or treated with DAC
NB cells were plated at day 0 and treated 24 h later with 3 μM DAC during 3 days. The NB cell line panel consisted of CHP-134, CHP-901, CHP-902R, CLB-GA, GI-C-IN, IMR-32, LA-N-1, LA-N-2, N206, NBL-S, NGP, NLF, NMB, SH-SY5Y, SH-SY5Y TrkA, SJNB-1, SJNB-8, SJNB-12, SK-N-AS, SMS-KAN, SMS-KCNR and STA-NB-1.2. Cells were harvested and RNA was extracted for real-time quantitative PCR using the RNeasy Mini kit (Qiagen). Primers for CADM1 and four stably expressed reference genes (HPRT1, GAPDH, YWHAZ, HMBS) are available in the public RTPrimerDB database [27–29] (gene (RTPrimerDB-ID): HPRT1 (7), GAPDH (3), YWHAZ (7), HMBS (4) and CADM1 (3770)). DNAse treatment and real-time PCR analysis were performed as previously reported .
A cell line panel was composed including CHP-134, CHP-901, CLB-GA, GI-C-IN-1, LA-N-1, LA-N-5, NB-13, NB-19, N206, NBL-S, NGP, NLF, SHEP, SJNB-1, SJNB-6, SJNB-8, SJNB-10, SJNB-1.2, SK-N-AS, SKN-BE (1n), SK-N-BE (2c), SK-N-FI, SMS-KAN, SMS-KCNR, STA-NB-3, STA-NB-8, STA-NB-9, STA-NB-10, STA-NB-1.2, TR-14 and UHG-NP. Approximately 1 μg of DNA was modified with sodium bisulphite using the EZ DNA Methylation Kit (Zymo Research). The resulting modified DNA was used as a template in a PCR reaction. Briefly, PCR reaction was carried out in a 50 μl reaction containing 30 ng of bisulphite modified DNA, 1× Platinum Taq PCR reaction buffer (Invitrogen), 6 mM MgCl2, 200 μm of each dNTP, 1.25 U Platinum Taq polymerase (Invitrogen) and 300 nm of each primer. PCR cycling conditions were as follows: 4 min of initial denaturation at 93°C, followed by 40 cycles of 30 s denaturation at 93°C, 50 s of annealing at 57°C and 30 s of extension at 72°C, followed by a final extension step of 4 min at 72°C. Primers were developed using Bisearch [31, 32] and tested for specificity using our in house developed methBLAST software  and are available in the public methPrimerDB database (ID 243). The primers were designed to assess the methylation status of 20 CpG sites situated upstream of the transcription start 4 μl of PCR product was cloned into TOPO-TA vector and transformed into E. coli TOP10F' cells (Invitrogen) and 10 positive clones were randomly selected for direct sequencing using universal M13 primers. As a control, normal human genomic DNA and SssI methylase (New England Biolabs) treated DNA were processed respectively as a negative and a positive control for methylation after bisulphite modification and evaluated for overall bisulphite conversion using sequencing. A 100% conversion was reached (data not shown). Furthermore, we verified for preferential amplification of unmethylated or methylated alleles through sequencing of artificial mixture samples from normal DNA and different concentrations of methylated DNA. We could confirm no amplification bias (data not shown).
Results and Discussion
In an attempt to more accurately map 11q losses and to identify potential homozygous losses that were previously missed using low resolution analyses, we performed high resolution BAC array CGH copy number profiling on 100 primary NB tumours (16 stage 1, 12 stage 2, 21 stage 3, 40 stage 4 and 11 stage 4S tumours) and 29 cell lines using a 1 Mb BAC array supplemented with a chromosome 11q tiling path (position 63–72 Mb (11q13) and 84–134 Mb (11q14.1-qter)). Three patterns of chromosome 11 loss could be observed. Twenty-four (24%) tumours showed loss of the entire chromosome 11. These tumours consisted of predominantly low stage tumours with an overall pattern of numerical imbalances. In addition, also 3 cell lines (10%) showed whole chromosome 11 loss. Seventeen tumours (17%) and 16 cell lines (55%) showed partial 11q loss. Most deletions were large encompassing the distal end of 11q, in keeping with data from literature [6–10, 15, 34, 35]. In cell lines CLB-GA, SMS-KAN, and NB-5, small interstitial losses were detected (respectively minimal 42.90, 30.84, and 4.76 Mb in size). To our knowledge this is the first report of interstitial 11q loss in SMS-KAN and NB-5 cell lines. No homozygous 11q losses were observed.
Two common minimal regions of loss (shortest region of overlap, SRO) at 11q23 could be delineated (Figure 1). The telomeric SRO (SRO1, 3.54–3.72 Mb, 11q23.2-q23.3) is delineated by the breakpoints in N206 (centromeric border) and SMS-KAN (telomeric border) and contains 31 genes. More centromeric, a second SRO (SRO2, 1.65–1.79 Mb, 11q23.1-q23.2) is delineated by the breakpoints in STA-NB3 (centromeric border) and NB-5 (telomeric border) and contains 12 genes. Recent evidence indicates that also miRNAs can function as tumour suppressor genes. According to the miRBase, [36, 37] (August 2007, version 10.0) no miRNAs were however located within the defined SROs. The finding that the most distal SRO partially overlaps with the SRO defined by Guo et al.  and both SROs are positioned within the critical region defined by Maris et al.  validates our mapping efforts. Nevertheless, it is important to note that SRO mapping should always be interpreted with care as individual observations can dramatically reduce the size or location of the SRO. This is especially the case when including cell line data, in which alterations could be induced during culturing. In this respect, it should therefore be noted that NB cell lines have been proven to be very stable over time and that, although rare, interstitial 11q loss has also been described in primary NB tumours [6, 7, 15, 38], further underscoring the relevance of our defined SRO regions.
To prioritize potential NB 11q candidate genes within the two critical regions, we performed a meta-analysis on publicly available gene expression data of 692 tumours from seven independent studies [20–26]. An overview of the obtained results is given in Additional file 2. In a first analysis we determined which genes in our SROs have prognostic power using Cox regression analysis. In total, 24 genes were prognostically significant in one or more studies (see Additional file 2), of which four (CADM1, NCAM1, CD3E, ZNF259) were significant in three or more studies. In the four studies for which DNA copy number data were available, we determined in a second analysis which genes could discriminate between the various NB subgroups (subtype 1, 2A, 2B) (Kruskal-Wallis test). This multi-group analysis was specifically applied instead of a binary comparison between 11q normal versus 11q deleted samples in order to avoid combining tumours with a distinct biology and prognosis. Thirteen genes had discriminative power, of which five were significant in at least two studies (CADM1, NCAM1, BACE1, ZNF259, KIAA0999). Genes that were prognostic significant in more than three independent studies and had subgroup discriminating power in more than two independent studies were considered as strong candidate genes (CADM1, NCAM1 and ZNF259).
ZNF259 is a zinc finger protein involved in normal cell cycle progression . ZNF259 has been suggested as a modifier gene in spinal muscular atrophy , but no direct link with carcinogenesis has been established thus far. NCAM1 is a neural cell adhesion molecule with an important role in tumour migration and metastasis. In a rat model it was demonstrated that NCAM transfected glioma cells show a decreased invasion compared to low NCAM expressing control cells . Similarly, loss of NCAM in a transgenic mouse model of pancreatic β cell carcinogenesis caused the formation of lymph node metastasis, whereas NCAM overexpression prevented tumour metastasis . In several tumour types, including colon carcinoma , gastrointestinal neoplasia  and astrocytic tumours , downregulation of NCAM expression and an associated poor survival has been described. In chemoresistant NB cells an enhanced invasive capacity was observed due to downregulation of NCAM adhesion receptors . CADM1 encodes a cellular adhesion molecule with a role in synaptic formation of neural cells . This gene, also known as IGSF4 (immunoglobulin superfamily, member 4) or TSLC1 (tumour suppressor in lung cancer-1), is a bona fide tumour suppressor gene in non-small lung cancer and was identified by functional complementation of a frequently deleted region in nude mice . Reduced or lost expression has been described in a variety of human cancers and in vivo suppression of tumorigenesis after restoration of expression has been demonstrated in esophageal cancer, cervical cancer and prostate cancer [49–52].
To select the most promising candidate gene within the three remaining candidate genes, we integrated the gene expression profile of normal NB precursor cells (see Additional file 2). The transcriptome information of these normal fetal neuroblasts represents a unique resource for identifying genes contributing to NB development and the malignant phenotype . Remarkably CADM1 was the only gene in our candidate list that was differentially expressed in NB versus neuroblasts using the Rank Product algorithm. A significantly lower expression of CADM1 was observed when comparing unfavourable NB to neuroblasts . This is in agreement with the Kruskall-Wallis data showing CADM1 downregulation in unfavourable NB compared to favourable subtype 1 NB. Consequently, we consider CADM1 to be a compelling candidate NB tumour suppressor gene for further analysis.
In several cancer types, two-hit inactivation of CADM1 occurs by LOH and hypermethylation of the promoter region, whereas CADM1 mutations are rare. Indeed, it was shown that only two CADM1 sequence variants could be picked up in a panel of 31 NB cell lines , while in a study of 25 NB with proven CADM1 allelic loss no mutations were found in the remaining allele . Epigenetic modification has been described for several genes in NB [55, 56]. To study if the remaining allele is subject to methylation and resulting epigenetic silencing, we measured CADM1 mRNA expression in a panel of 22 NB cell lines before and after exposure to the demethylating agent 5-aza-2'-deoxycytidine (DAC). Four out of 13 cell lines with loss at the CADM1 locus showed a more than two fold upregulation of CADM1 mRNA expression after DAC treatment. This was also the case for two of the nine cell lines with normal CADM1 copy number. Representative examples of mRNA expression before and after treatment with DAC are given in Figure 2. As the maximal observed expression upregulation was 4.4 fold in cell line CHP-134, these data were inconclusive to predict the involvement of methylation associated CADM1 silencing in NB. Consequently, we performed bisulphite sequencing of twenty promoter region CpG dinucleotides in a panel of 31 NB cell lines (including the cell lines that were studied in the DAC experiment). Note, that the investigated CpG sites also cover the 93 bp consensus CpG region described by Kuramochi et al. . Methylation of all studied CpG sites was identified in 12.5, 20 and 11% of analysed clones in the cell lines NGP, SH-EP and SK-N-BE(2c). It should be noted that NGP was also present in the panel of cell lines treated with DAC, where a more than two fold upregulation was observed after DAC treatment. However, in all other cell lines only sporadic CpG methylation was detected at low frequency (< 30%). Based on this evidence and taken in account that cell lines typically show a higher extent of CpG island methylation compared to primary NB, epigenetic CADM1 inactivation probably does not play a major role in NB. The latter could explain the observed limited fold changes in the DAC re-expression data.
While CADM1 forms an excellent positional and functional candidate tumour suppressor gene for NB, we can not exclude a role for the 2 other identified candidates suppressor genes ZNF259 and NCAM1. Based on our own methylation and previously published mutation data, CADM1 is most likely not inactivated through a two-hit Knudson mechanism. However, novel concepts in tumour suppressor genetics such as haploinsufficiency currently challenge the two-hit paradigm and are gaining importance in cancer genetics. Hence, an intriguing possibility is sensitivity of CADM1 to gene dosage. In this respect, the observed reduced expression of CADM1 in unfavourable NB as compared to normal neuroblasts is interesting. CADM1 belongs to the family of immunoglobulin cell adhesion molecules and encodes a transmembrane glycoprotein involved in cell-cell interactions. As disruption or altering of cell-cell interactions is an important feature of tumour invasion and metastasis, it could therefore be hypothesised that one working CADM1 allele is insufficient to accomplish the normal gene function. This would make NB cells more aggressive and amenable for migration, invasion and metastasis formation. Furthermore, although the biological cascade of the CADM1 protein is not yet elucidated, it is interesting to note that EPB41L3 – an established binding partner of CADM1  – has been described as higher expressed in mass screening NB tumours versus aggressive NB  and favourable NB versus unfavourable NB . In 4/5 data sets, for which expression data on CADM1 and EPB41L3 were available, also a significant positive correlation between CADM1 and EPB41L3 mRNA expression was found (p < 0.05). Further research should therefore be performed to elucidate the involvement of CADM1 in NB.
Tumorigenesis is a multistep process characterised by a multitude of genetic and epigenetic alterations. In view of the limited success of classical genomic approaches towards the identification of NB tumour suppressor genes, we integrated high resolution array CGH with transcriptome analysis and identified CADM1 as a candidate 11q tumour suppressor gene with prognostic power. CADM1 does not seem to be inactivated through a two-hit mechanism in NB. However, due to the frequent involvement of CADM1 in human tumorigenesis and its presumed role in invasion, a hallmark of malignancy, further studies should be performed to confirm CADM1 involvement in NB and elucidate its mode of action.
Weiss WA, Godfrey T, Francisco C, Bishop JM: Genome-wide screen for allelic imbalance in a mouse model for neuroblastoma. Cancer Res. 2000, 60 (9): 2483-2487.
Mosse YP, Laudenslager M, Khazi D, Carlisle AJ, Winter CL, Rappaport E, Maris JM: Germline PHOX2B mutation in hereditary neuroblastoma. Am J Hum Genet. 2004, 75 (4): 727-730. 10.1086/424530.
Trochet D, Bourdeaut F, Janoueix-Lerosey I, Deville A, de Pontual L, Schleiermacher G, Coze C, Philip N, Frebourg T, Munnich A, Lyonnet S, Delattre O, Amiel J: Germline mutations of the paired-like homeobox 2B (PHOX2B) gene in neuroblastoma. Am J Hum Genet. 2004, 74 (4): 761-764. 10.1086/383253.
van Limpt V, Schramm A, van Lakeman A, Sluis P, Chan A, van Noesel M, Baas F, Caron H, Eggert A, Versteeg R: The Phox2B homeobox gene is mutated in sporadic neuroblastomas. Oncogene. 2004, 23 (57): 9280-9288.
Attiyeh EF, London WB, Mosse YP, Wang Q, Winter C, Khazi D, McGrady PW, Seeger RC, Look AT, Shimada H, Brodeur GM, Cohn SL, Matthay KK, Maris JM: Chromosome 1p and 11q deletions and outcome in neuroblastoma. N Engl J Med. 2005, 353 (21): 2243-2253. 10.1056/NEJMoa052399.
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 (10): 2280-2299. 10.1200/JCO.2005.06.104.
Guo C, White PS, Weiss MJ, Hogarty MD, Thompson PM, Stram DO, Gerbing R, Matthay KK, Seeger RC, Brodeur GM, Maris JM: Allelic deletion at 11q23 is common in MYCN single copy neuroblastomas. Oncogene. 1999, 18 (35): 4948-4957. 10.1038/sj.onc.1202887.
Plantaz D, Vandesompele J, Van Roy N, Lastowska M, Bown N, Combaret V, Favrot MC, Delattre O, Michon J, Benard J, Hartmann O, Nicholson JC, Ross FM, Brinkschmidt C, Laureys G, Caron H, Matthay KK, Feuerstein BG, Speleman F: Comparative genomic hybridization (CGH) analysis of stage 4 neuroblastoma reveals high frequency of 11q deletion in tumors lacking MYCN amplification. Int J Cancer. 2001, 91 (5): 680-686. 10.1002/1097-0215(200002)9999:9999<::AID-IJC1114>3.0.CO;2-R.
Vandesompele J, Speleman F, Van Roy N, Laureys G, Brinskchmidt C, Christiansen H, Lampert F, Lastowska M, Bown N, Pearson A, Nicholson JC, Ross F, Combaret V, Delattre O, Feuerstein BG, Plantaz D: Multicentre analysis of patterns of DNA gains and losses in 204 neuroblastoma tumors: how many genetic subgroups are there?. Med Pediatr Oncol. 2001, 36 (1): 5-10. 10.1002/1096-911X(20010101)36:1<5::AID-MPO1003>3.0.CO;2-E.
Vandesompele J, Van Roy N, Van Gele M, Laureys G, Ambros P, Heimann P, Devalck C, Schuuring E, Brock P, Otten J, Gyselinck J, De Paepe A, Speleman F: Genetic heterogeneity of neuroblastoma studied by comparative genomic hybridization. Genes Chromosomes Cancer. 1998, 23 (2): 141-152. 10.1002/(SICI)1098-2264(199810)23:2<141::AID-GCC7>3.0.CO;2-2.
Koiffmann CP, Gonzalez CH, Vianna-Morgante AM, Kim CA, Odone-Filho V, Wajntal A: Neuroblastoma in a boy with MCA/MR syndrome, deletion 11q, and duplication 12q. Am J Med Genet. 1995, 58 (1): 46-49. 10.1002/ajmg.1320580110.
Mosse Y, Greshock J, King A, Khazi D, Weber BL, Maris JM: Identification and high-resolution mapping of a constitutional 11q deletion in an infant with multifocal neuroblastoma. Lancet Oncol. 2003, 4 (12): 769-771. 10.1016/S1470-2045(03)01283-X.
Bader SA, Fasching C, Brodeur GM, Stanbridge EJ: Dissociation of suppression of tumorigenicity and differentiation in vitro effected by transfer of single human chromosomes into human neuroblastoma cells. Cell Growth Differ. 1991, 2 (5): 245-255.
De Preter K, Vandesompele J, Menten B, Carr P, Fiegler H, Edsjo A, Carter NP, Yigit N, Waelput W, Van Roy N, Bader S, Pahlman S, Speleman F: Positional and functional mapping of a neuroblastoma differentiation gene on chromosome 11. BMC Genomics. 2005, 6: 97-10.1186/1471-2164-6-97.
Maris JM, Guo C, White PS, Hogarty MD, Thompson PM, Stram DO, Gerbing R, Matthay KK, Seeger RC, Brodeur GM: Allelic deletion at chromosome bands 11q14-23 is common in neuroblastoma. Med Pediatr Oncol. 2001, 36 (1): 24-27. 10.1002/1096-911X(20010101)36:1<24::AID-MPO1007>3.0.CO;2-7.
Michels E, Vandesompele J, De Preter K, Hoebeeck J, Vermeulen J, Schramm A, Molenaar JJ, Menten B, Marques B, Stallings RL, Combaret V, Devalck C, De Paepe A, Versteeg R, Eggert A, Laureys G, Van Roy N, Speleman F: ArrayCGH-based classification of neuroblastoma into genomic subgroups. Genes Chromosomes Cancer. 2007, 46 (12): 1098-1108. 10.1002/gcc.20496.
Brodeur GM, Pritchard J, Berthold F, Carlsen NL, Castel V, Castelberry RP, De Bernardi B, Evans AE, Favrot M, Hedborg F, et al: Revisions of the international criteria for neuroblastoma diagnosis, staging, and response to treatment. J Clin Oncol. 1993, 11 (8): 1466-1477.
Menten B, Pattyn F, De Preter K, Robbrecht P, Michels E, Buysse K, Mortier G, De Paepe A, van Vooren S, Vermeesch J, Moreau Y, De Moor B, Vermeulen S, Speleman F, Vandesompele J: arrayCGHbase: an analysis platform for comparative genomic hybridization microarrays. BMC Bioinformatics. 2005, 6: 124-10.1186/1471-2105-6-124.
Berwanger B, Hartmann O, Bergmann E, Bernard S, Nielsen D, Krause M, Kartal A, Flynn D, Wiedemeyer R, Schwab M, Schafer H, Christiansen H, Eilers M: Loss of a FYN-regulated differentiation and growth arrest pathway in advanced stage neuroblastoma. Cancer Cell. 2002, 2 (5): 377-386. 10.1016/S1535-6108(02)00179-4.
De Preter K, Vandesompele J, Heimann P, Yigit N, Beckman S, Schramm A, Eggert A, Stallings RL, Benoit Y, Renard M, Paepe AD, Laureys G, Pahlman S, Speleman F: Human fetal neuroblast and neuroblastoma transcriptome analysis confirms neuroblast origin and highlights neuroblastoma candidate genes. Genome Biol. 2006, 7 (9): R84-10.1186/gb-2006-7-9-r84.
McArdle L, McDermott M, Purcell R, Grehan D, O'Meara A, Breatnach F, Catchpoole D, Culhane AC, Jeffery I, Gallagher WM, Stallings RL: Oligonucleotide microarray analysis of gene expression in neuroblastoma displaying loss of chromosome 11q. Carcinogenesis. 2004, 25 (9): 1599-1609. 10.1093/carcin/bgh173.
Oberthuer A, Berthold F, Warnat P, Hero B, Kahlert Y, Spitz R, Ernestus K, Konig R, Haas S, Eils R, Schwab M, Brors B, Westermann F, Fischer M: Customized oligonucleotide microarray gene expression-based classification of neuroblastoma patients outperforms current clinical risk stratification. J Clin Oncol. 2006, 24 (31): 5070-5078. 10.1200/JCO.2006.06.1879.
Ohira M, Oba S, Nakamura Y, Isogai E, Kaneko S, Nakagawa A, Hirata T, Kubo H, Goto T, Yamada S, Yoshida Y, Fuchioka M, Ishii S, Nakagawara A: Expression profiling using a tumor-specific cDNA microarray predicts the prognosis of intermediate risk neuroblastomas. Cancer Cell. 2005, 7 (4): 337-350. 10.1016/j.ccr.2005.03.019.
Schramm A, Schulte JH, Klein-Hitpass L, Havers W, Sieverts H, Berwanger B, Christiansen H, Warnat P, Brors B, Eils J, Eils R, Eggert A: Prediction of clinical outcome and biological characterization of neuroblastoma by expression profiling. Oncogene. 2005, 24 (53): 7902-7912. 10.1038/sj.onc.1208936.
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 NK, 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 (12): 6050-6062. 10.1158/0008-5472.CAN-05-4618.
Pattyn F, Robbrecht P, De Paepe A, Speleman F, Vandesompele J: RTPrimerDB: the real-time PCR primer and probe database, major update 2006. Nucleic Acids Res. 2006, 34 (Database issue): D684-8. 10.1093/nar/gkj155.
Pattyn F, Speleman F, De Paepe A, Vandesompele J: RTPrimerDB: the real-time PCR primer and probe database. Nucleic Acids Res. 2003, 31 (1): 122-123. 10.1093/nar/gkg011.
Vandesompele J, De Paepe A, Speleman F: Elimination of primer-dimer artifacts and genomic coamplification using a two-step SYBR green I real-time RT-PCR. Anal Biochem. 2002, 303 (1): 95-98. 10.1006/abio.2001.5564.
Aranyi T, Varadi A, Simon I, Tusnady GE: The BiSearch web server. BMC Bioinformatics. 2006, 7: 431-10.1186/1471-2105-7-431.
Tusnady GE, Simon I, Varadi A, Aranyi T: BiSearch: primer-design and search tool for PCR on bisulfite-treated genomes. Nucleic Acids Res. 2005, 33 (1): e9-10.1093/nar/gni012.
Pattyn F, Hoebeeck J, Robbrecht P, Michels E, De Paepe A, Bottu G, Coornaert D, Herzog R, Speleman F, Vandesompele J: methBLAST and methPrimerDB: web-tools for PCR based methylation analysis. BMC Bioinformatics. 2006, 7: 496-10.1186/1471-2105-7-496.
Breen CJ, O'Meara A, McDermott M, Mullarkey M, Stallings RL: Coordinate deletion of chromosome 3p and 11q in neuroblastoma detected by comparative genomic hybridization. Cancer Genet Cytogenet. 2000, 120 (1): 44-49. 10.1016/S0165-4608(99)00252-6.
Luttikhuis ME, Powell JE, Rees SA, Genus T, Chughtai S, Ramani P, Mann JR, McConville CM: Neuroblastomas with chromosome 11q loss and single copy MYCN comprise a biologically distinct group of tumours with adverse prognosis. Br J Cancer. 2001, 85 (4): 531-537. 10.1054/bjoc.2001.1960.
Griffiths-Jones S: The microRNA Registry. Nucleic Acids Res. 2004, 32 (Database issue): D109-11. 10.1093/nar/gkh023.
Griffiths-Jones S, Grocock RJ, van Dongen S, Bateman A, Enright AJ: miRBase: microRNA sequences, targets and gene nomenclature. Nucleic Acids Res. 2006, 34 (Database issue): D140-4. 10.1093/nar/gkj112.
Spitz R, Oberthuer A, Zapatka M, Brors B, Hero B, Ernestus K, Oestreich J, Fischer M, Simon T, Berthold F: Oligonucleotide array-based comparative genomic hybridization (aCGH) of 90 neuroblastomas reveals aberration patterns closely associated with relapse pattern and outcome. Genes Chromosomes Cancer. 2006, 45 (12): 1130-1142. 10.1002/gcc.20376.
Gangwani L: Deficiency of the zinc finger protein ZPR1 causes defects in transcription and cell cycle progression. J Biol Chem. 2006, 281 (52): 40330-40340. 10.1074/jbc.M608165200.
Doran B, Gherbesi N, Hendricks G, Flavell RA, Davis RJ, Gangwani L: Deficiency of the zinc finger protein ZPR1 causes neurodegeneration. Proc Natl Acad Sci U S A. 2006, 103 (19): 7471-7475. 10.1073/pnas.0602057103.
Owens GC, Orr EA, DeMasters BK, Muschel RJ, Berens ME, Kruse CA: Overexpression of a transmembrane isoform of neural cell adhesion molecule alters the invasiveness of rat CNS-1 glioma. Cancer Res. 1998, 58 (9): 2020-2028.
Perl AK, Dahl U, Wilgenbus P, Cremer H, Semb H, Christofori G: Reduced expression of neural cell adhesion molecule induces metastatic dissemination of pancreatic beta tumor cells. Nat Med. 1999, 5 (3): 286-291. 10.1038/6502.
Roesler J, Srivatsan E, Moatamed F, Peters J, Livingston EH: Tumor suppressor activity of neural cell adhesion molecule in colon carcinoma. Am J Surg. 1997, 174 (3): 251-257. 10.1016/S0002-9610(97)00142-6.
Fogar P, Basso D, Pasquali C, De Paoli M, Sperti C, Roveroni G, Pedrazzoli S, Plebani M: Neural cell adhesion molecule (N-CAM) in gastrointestinal neoplasias. Anticancer Res. 1997, 17 (2B): 1227-1230.
Sasaki H, Yoshida K, Ikeda E, Asou H, Inaba M, Otani M, Kawase T: Expression of the neural cell adhesion molecule in astrocytic tumors: an inverse correlation with malignancy. Cancer. 1998, 82 (10): 1921-1931. 10.1002/(SICI)1097-0142(19980515)82:10<1921::AID-CNCR16>3.0.CO;2-V.
Blaheta RA, Daher FH, Michaelis M, Hasenberg C, Weich EM, Jonas D, Kotchetkov R, Doerr HW, Cinatl J: Chemoresistance induces enhanced adhesion and transendothelial penetration of neuroblastoma cells by down-regulating NCAM surface expression. BMC Cancer. 2006, 6: 294-10.1186/1471-2407-6-294.
Biederer T, Sara Y, Mozhayeva M, Atasoy D, Liu X, Kavalali ET, Sudhof TC: SynCAM, a synaptic adhesion molecule that drives synapse assembly. Science. 2002, 297 (5586): 1525-1531. 10.1126/science.1072356.
Kuramochi M, Fukuhara H, Nobukuni T, Kanbe T, Maruyama T, Ghosh HP, Pletcher M, Isomura M, Onizuka M, Kitamura T, Sekiya T, Reeves RH, Murakami Y: TSLC1 is a tumor-suppressor gene in human non-small-cell lung cancer. Nat Genet. 2001, 27 (4): 427-430. 10.1038/86934.
Fukuhara H, Kuramochi M, Fukami T, Kasahara K, Furuhata M, Nobukuni T, Maruyama T, Isogai K, Sekiya T, Shuin T, Kitamura T, Reeves RH, Murakami Y: Promoter methylation of TSLC1 and tumor suppression by its gene product in human prostate cancer. Jpn J Cancer Res. 2002, 93 (6): 605-609.
Ito T, Shimada Y, Hashimoto Y, Kaganoi J, Kan T, Watanabe G, Murakami Y, Imamura M: Involvement of TSLC1 in progression of esophageal squamous cell carcinoma. Cancer Res. 2003, 63 (19): 6320-6326.
Murakami Y: Involvement of a cell adhesion molecule, TSLC1/IGSF4, in human oncogenesis. Cancer Sci. 2005, 96 (9): 543-552. 10.1111/j.1349-7006.2005.00089.x.
Steenbergen RD, Kramer D, Braakhuis BJ, Stern PL, Verheijen RH, Meijer CJ, Snijders PJ: TSLC1 gene silencing in cervical cancer cell lines and cervical neoplasia. J Natl Cancer Inst. 2004, 96 (4): 294-305.
Attiyeh EF, Mosse Y, Wang Q, Winter C, Deepa K, Pecor KM, Maris JM: Functional and sequence validation of candidate chromosome 11 neuroblastoma suppressor genes: Los Angeles, California.2006, , 45-
Gomyo H, Arai Y, Tanigami A, Murakami Y, Hattori M, Hosoda F, Arai K, Aikawa Y, Tsuda H, Hirohashi S, Asakawa S, Shimizu N, Soeda E, Sakaki Y, Ohki M: A 2-Mb sequence-ready contig map and a novel immunoglobulin superfamily gene IGSF4 in the LOH region of chromosome 11q23.2. Genomics. 1999, 62 (2): 139-146. 10.1006/geno.1999.6001.
Alaminos M, Davalos V, Cheung NK, Gerald WL, Esteller M: Clustering of gene hypermethylation associated with clinical risk groups in neuroblastoma. J Natl Cancer Inst. 2004, 96 (16): 1208-1219.
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 (36): 5619-5628. 10.1038/sj.onc.1208722.
Yageta M, Kuramochi M, Masuda M, Fukami T, Fukuhara H, Maruyama T, Shibuya M, Murakami Y: Direct association of TSLC1 and DAL-1, two distinct tumor suppressor proteins in lung cancer. Cancer Res. 2002, 62 (18): 5129-5133.
Krause A, Combaret V, Iacono I, Lacroix B, Compagnon C, Bergeron C, Valsesia-Wittmann S, Leissner P, Mougin B, Puisieux A: Genome-wide analysis of gene expression in neuroblastomas detected by mass screening. Cancer Lett. 2005, 225 (1): 111-120. 10.1016/j.canlet.2004.10.035.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2407/8/173/prepub
We thank Justine Nuytens for assistance with bisulphite sequencing. We thank Raymond Stallings, Miki Ohira, Rüdiger Spitz, John Maris and the Children's Oncology group for providing updated clinical patient information.
Jasmien Hoebeeck is a postdoctoral researcher supported by a grant of the Ghent University (BOF 01P07406). Jo Vandesompele and Katleen De Preter are postdoctoral researchers of the Fund for Scientific Research (FWO) Flanders.
This work was supported by the Kinderkankerfonds (a nonprofit childhood cancer foundation under Belgian law), The Fund for Scientific Research, Flanders (Krediet aan Navorsers 1.5.243.05), FWO-grant G.0185.04, BOF-grant 011F1200 and 011B4300, concerted research fund (GOA, nr. 12051203) and the Stichting tegen Kanker. This text presents research results of the Belgian program of Interuniversity Poles of Attraction initiated by the Belgian State, Prime Minister's Office, Science Policy Programming.
The authors declare that they have no competing interests.
EM carried out the genomic studies and drafted the manuscript. JH performed the methylation study. KDP carried out the meta-analysis of neuroblastoma mRNA expression data. GL, AS, AE and BB provided tumour material. JV and FS designed and coordinated the study, and edited the manuscript. All authors have reviewed the manuscript and FS, JV and ADP were the final editors of the manuscript.
Electronic supplementary material
Additional file 2: . Meta-analysis on publicly available neuroblastoma gene expression data sets for genes residing in the two delineated critical regions of loss. Genes that were statistically significant in one or more studies are represented and ordered from centromere to telomere. Black boxes represent a significant p-value (<0.05) in the respective study. Grey boxes represent absence of a gene specific probe on the expression array (nb: normal neuroblasts, NB: neuroblastoma, F: favorable and UF: unfavorable prognosis according to De Preter et al., 2006 ). (XLS 22 KB)
About this article
Cite this article
Michels, E., Hoebeeck, J., De Preter, K. et al. CADM1 is a strong neuroblastoma candidate gene that maps within a 3.72 Mb critical region of loss on 11q23. BMC Cancer 8, 173 (2008). https://doi.org/10.1186/1471-2407-8-173
- Spinal Muscular Atrophy
- MYCN Amplification
- Copy Number Profile
- Bisulphite Sequencing
- Microcell Mediate Chromosome Transfer