A cluster of genes located in 1p36 are down-regulated in neuroblastomas with poor prognosis, but not due to CpG island methylation
Molecular Cancervolume 4, Article number: 10 (2005)
A common feature of neuroblastoma tumours are partial deletions of the short arm of chromosome 1 (1p-deletions). This is indicative of a neuroblastoma tumour suppressor gene being located in the region. Several groups including our have been studying candidate neuroblastoma genes in the region, but no gene/genes have yet been found that fulfil the criteria for being a neuroblastoma tumour suppressor. Since frequent mutations have not been detected, we have now analyzed the expression and promoter CpG island methylation status of the genes UBE4B, KIF1B, PGD, APITD1, DFFA and PEX14 in the 1p36.22 region in order to find an explanation for a possible down-regulation of this region.
The current study shows that gene transcripts in high stage neuroblastoma tumours are significantly down-regulated compared to those in low stage tumours in the 1p36.22 region. CpG island methylation does not seem to be the mechanism of down-regulation for most of the genes tested, since no methylation was detected in the fragments analyzed. One exception is the CpG island of APITD1. Methylation of this gene is also seen in blood from control individuals and is therefore not believed to participate in tumour development.
The genes UBE4B, KIF1B, PGD, APITD1, DFFA and PEX14 are down-regulated in high stage NB tumours, a feature that can not be explained by CpG island methylation.
Neuroblastoma (NB) is the most common paediatric solid tumour, responsible for 15% of cancer deaths of childhood. It is a tumour of the postganglionic sympathetic nervous system that develops from immature or dedifferentiated neural crest-derived cells . The distal part of chromosome 1p shows loss of heterozygosity (LOH) in 20–40% of NB tumours and has therefore been alleged to contain one or more tumour suppressor genes. We and others have previously analyzed the chromosomal region 1p36.2-3 [2–12] and we have recently focused on the gene region involving the genes: UBE4B-KIF1B-PGD-APITD1-CORT-DFFA-PEX14. These genes have been analyzed for mutations and a few have been found in rare tumours [13–17].
In search of tumour suppressor genes, the focus has in the last years moved towards epigenetics and methylation of promoter regions. Methylation of cytosines in CpG-dinucleotides is a common modification in mammalian genomes. Methylated cytosines are more susceptible to deamination, which have lead to an erosion of the number of CpG-dinucleotides. The vast majority of CpGs resides within repetitive elements and is methylated. There are also stretches of DNA rich in CpG that are gene associated, i.e. CpG islands, which are normally unmethylated . Methylation is generally associated with repression of transcription. Gene regulation by methylation includes tissue-specific regulation during development and processes as X-chromosome inactivation and genomic imprinting, reviewed by Herman and co-workers. . Cancer is associated with a genome-wide hypomethylation and a more gene-specific hypermethylation. Hypermethylation of CpG islands has been shown to be a common mechanism for the inactivation of tumour suppressor genes and is found in a wide range of tumour types [20–23]. According to the Knudson two-hit hypothesis, two successive mutations are required to inactivate a tumour suppressor gene and turn a normal cell into a malignant one. Inactivation could be due to deletions or mutations. Epigenetic events, such as hypermethylation of promoter-associated CpG islands have come in focus during the last decade as a route to inactivation .
The most common way to analyze methylation status is based on bisulphite modification of DNA . In this reaction, unmethylated cytosines are deaminated to uracil, while methylated cytosines remain unconverted. The region of interest can be amplified using non-discriminating primers, amplifying both methylated and unmethylated templates in one reaction, or with methylation-specific PCR (MSP) in which methylated and unmethylated templates are amplified in separate reactions .
Some genes have been analyzed in NB tumours focusing on methylation status. For example, CASP8, on 2q33, was one of the first genes found to be methylated in neuroblastoma, with a frequency of about 40 % in primary NB tumours . On 3p21.3, RASSF1A and BLU have been shown to be frequently methylated [27–29]. In a study by van Noesel and co-workers, 34 genes in 12 different chromosomal regions were analyzed in neuroblastoma cell lines . A total of six genes that were methylated in at least three of the 22 cell lines were found. CASP8, as already known, was one of these genes, also FLIP at 2q33 was methylated as well as four genes in the chromosome region 8p21; DR4, DR5, DCR1 and DCR2. Genes on 1p36 were also included in this study, none of which were found to be methylated. No genes in the NB tumour suppressor candidate region on 1p36.22 were included. Alaminos and co-workers have in a study screened 45 genes in NB cell lines . 34 genes were found to be methylated in at least one of the ten cell lines. The analysis also demonstrated that the percentage of methylated genes was higher in cell lines with MYCN-amplification than in those without. UBE4B, located in 1p36.22, was included in this study and found to be unmethylated.
The genes in the 1p36 consensus deleted region are all, except for CORT, associated with a CpG island in their respective promoter regions. The fact that this region is deleted in a subset of NB tumours and that frequent mutations are not found, makes it a possible candidate region for hypermethylation of genes. In the current study, we wanted to analyze the expression pattern of these genes in the 1p36.22 region using TaqMan technology and the methylation status of promoter regions associated with CpG islands, in a panel of NB cell lines and primary NB tumours.
GUSB was selected as an endogenous control for real-time PCR quantification, and further used as an internal reference for normalization. Real time-PCR studies of all the genes in the region on cDNA-samples showed a significant reduction of mRNA levels in high stage NB tumours compared to low stage tumour mRNA levels, ranging from 49% for PEX14 to 79% for APITD1 (Table 4) .
The average expression in NB cell lines compared to primary tumours varied between the different genes in the region. Cell lines showed a higher expression of APITD1 and DFFA compared to primary tumours. Lower expression in cell lines was seen in KIF1B and equal levels of cell lines and stage 2 primary tumours in UBE4B, PGD and PEX14 (Table 4; Fig. 1). The 1p-deleted cell line SK-N-AS generally showed lower expression of the genes in the region compared to the other NB cell lines tested (Fig. 1).
The analyzed fragments of the genes UBE4B, KIF1B, PGD, DFFA and PEX14 were not methylated in the panel of NB primary tumours and cell lines. For APITD1, the region -393 to -222 relative to initiation codon showed methylation. This region was therefore amplified and cloned with a T/A cloning kit in order to give a value of the percentage of clones being methylated (Fig. 2). The region contains 16 CpG sites and these sites show various degree of methylation, in different clones as well as in different samples. Methylation of cytosine in front of guanine is concentrated to the beginning of this "methylated region"; the three first cytosine bases showing the highest degree of methylated clones (Fig 3). No methylation was seen up-stream of -393 in the fragments analyzed. At some CpG sites in the region, methylation was never detected. Only one sample had no clones with methylation for any of the sites; the stage 2 tumour 20S9. Blood from control individuals also showed various degree of methylation at these sites. The cell lines did not show methylation down-stream of CpG -212 in the analyzed fragments, clones from tumours were not methylated down-stream of CpG -311. Healthy controls showed a smaller region of methylation, no methylated clones were detected down-stream of CpG -280.
Deletion of parts of chromosome 1p is a common feature of neuroblastoma tumours. We have previously narrowed down the shortest region of overlap of deletions (SRO) to 25 cM in our tumour material [6, 7]. By including a germ cell tumour with 1p deletion we could confine the SRO to 5 cM . This SRO was further confirmed when a cell line with a homozygous deletion of 500 kb was found within this region . We have previously screened the 500 kb region for mutations and only rare have been detected. The aim of this study was to investigate the expression of the genes in the region and to correlate this with the methylation status. TaqMan expression data showed a clear down-regulation in gene expression of genes in the region tested when comparing high stage NB tumours (stage 3 and 4) and low stage tumours (stage 2). This down-regulation was seen regardless if the high stage tumour was 1p-deleted or not (Fig. 1). This finding further supports that this region might be involved in NB development and progression. We can speculate about some central mechanism for the expression of the genes in the entire region that is affected in tumours with poor prognosis, either due to the deletion of one allele and/or by some other mechanism. An attractive explanation for this would be methylation of CpG islands that can lead to inactivation of the respective gene. This does not seem to be the case though, since none of the genes, APITD1 excepted, showed methylation. Furthermore, DNA from healthy controls is also methylated in the same region of the APITD1 gene, making methylation a less likely explanation for the difference in APITD1 gene expression between advanced and favourable primary NB tumours.
CpG methylation is not the only possible explanation for the down-regulation in gene expression shown in high stage tumours. Analysis of chromatin modifications could be the next step in a further analysis of this region, in search of a mechanism that could lead to down-regulation of the genes. Deficiencies in transcription factors and upstream elements could also account for the decrease in gene transcripts. Another possibility is that the variations are due to differences in the treatment of the tumours before and after surgery. Radiation and chemotherapy could probably have a vast affect on the expression patterns of a substantial amount of genes in the tumour cell. It is however unlikely that all the genes in the analyzed region would be affected by the treatment alone, rather this might explain the differences in the level of down-regulation between different genes.
The difference in expression between NB cell lines and primary NB could be explained by the difference in cell environment between cells in culture and in vivo. In primary tumours, there are always some non-NB cells in the RNA preparations that might affect the gene expression results. On the other hand, genetic events happen to the cells when they get immortalized in cell culture, some genes tend to be down-regulated and others up-regulated. We can also speculate that some of these genes, for example APITD1, might even be so crucial that only cells with expression of this gene are capable of surviving in culture.
In conclusion, in our aim to find a mechanism that could inactivate the genes in the candidate gene region we analyzed the expression status and the methylation profile of six genes. The techniques used were TaqMan real-time-RT-PCR technology and bisulphite DNA sequencing. Of the genes analyzed all were down-regulated in high staged NB tumours as compared to low stage tumours. Promoter methylation was not detected in the genes analyzed, except for the CpG island of APITD1. This methylation does not seem to be tumour specific, since methylation was also detected in healthy blood controls. Hence, the six genes UBE4B, KIF1B, PGD, APITD1, DFFA and PEX14 are down-regulated in high stage NB tumours, a feature that can not be explained by methylation, rather by a mechanism still remaining to be discovered.
Cell lines and patients
A panel of 10 tumours from primary NBs (4 Stage 2, all from patients with no evidence of decease at last follow-up, and 6 Stage 3 or 4 tumours, all from patients with adverse outcome) and 7 NB cell lines (IMR-32, SK-N-AS, SK-N-BE (2), SK-N-SH, SK-N-DZ, SK-N-F1 and SH-SY-5Y) were analyzed (Table 1).
Total RNA was extracted from frozen (-70°C) NB tumour tissue using RNeasy RNA extraction kit (Qiagen, Hilden, Germany). 2.4 μg total-RNA of each sample was reversed transcribed to cDNA using Superscript II (Amersham, Buckinghamshire, UK) and random hexamer primers, all according to supplier's protocol. All cDNAs were quality tested by amplification of the housekeeping genes UNPH and GAPDH.
Real time PCR -Endogenous control
To select the most appropriate endogenous control for the real-time PCR quantification analysis, we tested eight different primary NB samples of different stages for their expression levels of ten commonly used housekeeping genes with TaqMan Human Endogenous Control Plate, (Applied Biosystems, Foster City, CA). Analysis was performed according to supplier's protocol. GUSB (β-glucuronidase) and B2M (β2-microglobulin) showed least variations in ΔCT levels, and were expressed at constant levels in all samples regardless of NB-stage. GUSB was selected, and further used as an internal reference for normalization in the real-time PCR quantification analysis (Abel et al., submitted).
Real time PCR- TaqMan
TaqMan primers and probes were derived from the commercially available "TaqMan® Assays-on-Demand™ Gene Expression Products" (URL:http://myscience.appliedbiosystems.com). Real-time PCR was performed in 384-well plates using ABI PRISM® 7900HT Sequence Detection System (Applied Biosystems). Amplification reactions (10 μl) were carried out in duplicate with 0.1 μl template cDNA according to manufacturers protocol (Applied Biosystems). In each assay, a standard curve with six cDNA dilutions was recorded and two non-template controls were included.
Quantification was performed by the standard-curve method. The mean CT-value for duplicates were calculated, and the gene concentration (or gene copy numbers) of test samples was interpolated based on standard curves. All samples were normalized by dividing the concentration of the test gene with the concentration of the housekeeping gene β-glucuronidase (GUSB) in the same cDNA sample.
The logarithms of the expression levels were compared with Student's two-sided t-test on each group of tumours; low stage and high stage tumours.
DNA was phenol extracted with the use of phase lock gel (Eppendorf AG, Hamburg, Germany) according to standard procedure and was, with some minor changes, modified according to previously published papers [24, 33]. Briefly, 1 μg of genomic DNA was digested with restriction endonucleases that cut close but outside the region of interest. The DNA was then denaturated in 0.3 M freshly prepared NaOH at 40°C for 15 minutes. Sodium metabisulphite (Sigma-Aldrich CO, St Louis, MO) and urea, at a final concentration of 1.73 M and 5.36 M respectively, were added in order to sulphonate the unmethylated cytosines, along with hydroquinone (0.5 mM). Conversion was carried out at 55°C for 16 hours, with a temperature rise to 95°C for 30 seconds every third hour. DNA was purified with Wizard DNA clean up system (Promega Corporation, Madison, WI) according to the manufactures instructions and desulphonated in 0.3 M NaOH at 37°C for 15 minutes. DNA was then precipitated in ethanol, resuspended in distilled H2O and stored at -20°C.
Promoter analysis and DNA amplification
The putative promoter regions of the genes were predicted using Genomatix Promoter Inspector software (URL: http://www.genomatix.de) and CpG islands with MethPrimer software (URL: http://www.urogene.org/methprimer/; Table 2) . These regions, or parts of them, were amplified with one primer pair, or if needed, with semi-nested primers (Table 3). The methylation status was analyzed using bisulphite sequencing. Conditions for PCR amplification were 1× PCR Gold Buffer (Applied Biosystems), 0.5 mM dNTPs, 2.0–3.0 mM MgCl2, 0.4 μM of forward and reverse primers respectively and 1 unit of AmpliTaq Gold, in a total volume of 50 μl. Reactions were denatured at 95°C for 10 min followed by 5 cycles of 95°C for 1 min, 49–55°C for 2 min, 72°C for 3 min and 30 cycles of 95°C for 30 sec, 49–55°C for 2 min, 72°C for 1 min 30 sec and ending with 10 min extension at 72°C. The PCR products were immediately sequenced or cloned into a sequencing vector using the TOPO T/A cloning kit (Life Technologies Invitrogen, Carlsbad, CA), where after 10–30 clones were picked. PCR products were purified with ExoSAP-IT™ (USB Corporation, Cleveland, Ohio) and sequencing was carried out using forward or reverse primer with ABI Prism BigDye™ cycle sequencing Ready Reaction Kit (Applied Biosystems). The samples were analyzed in an ABI 3100 or an ABI 3730 Genetic Analyzer.
Maris JM, Matthay KK: Molecular biology of neuroblastoma. J Clin Oncol. 1999, 17: 2264-2279.
Caron H, Peter M, van Sluis P, Speleman F, de Kraker J, Laureys G, Michon J, Brugieres L, Voute PA, Westerveld A: Evidence for two tumour suppressor loci on chromosomal bands 1p35-36 involved in neuroblastoma: one probably imprinted, another associated with N-myc amplification. Hum Mol Genet. 1995, 4: 535-539.
Caron H, Spieker N, Godfried M, Veenstra M, van Sluis P, de Kraker J, Voute P, Versteeg R: Chromosome bands 1p35-36 contain two distinct neuroblastoma tumor suppressor loci, one of which is imprinted. Genes Chromosomes Cancer. 2001, 30: 168-174.
Maris JM, White PS, Beltinger CP, Sulman EP, Castleberry RP, Shuster JJ, Look AT, Brodeur GM: Significance of chromosome 1p loss of heterozygosity in neuroblastoma. Cancer Res. 1995, 55: 4664-4669.
Maris JM, Guo C, Blake D, White PS, Hogarty MD, Thompson PM, Rajalingam V, Gerbing R, Stram DO, Matthay KK, Seeger RC, Brodeur GM: Comprehensive analysis of chromosome 1p deletions in neuroblastoma. Med Pediatr Oncol. 2001, 36: 32-36.
Martinsson T, Sjöberg RM, Hedborg F, Kogner P: Deletion of chromosome 1p loci and microsatellite instability in neuroblastomas analyzed with short-tandem repeat polymorphisms. Cancer Res. 1995, 55: 5681-5686.
Martinsson T, Sjöberg RM, Hallstensson K, Nordling M, Hedborg F, Kogner P: Delimitation of a critical tumour suppressor region at distal 1p in neuroblastoma tumours. Eur J Cancer. 1997, 33: 1997-2001.
White PS, Thompson PM, Seifried BA, Sulman EP, Jensen SJ, Guo C, Maris JM, Hogarty MD, Allen C, Biegel JA, Matise TC, Gregory SG, Reynolds CP, Brodeur GM: Detailed molecular analysis of 1p36 in neuroblastoma. Med Pediatr Oncol. 2001, 36: 37-41.
White PS, Maris JM, Beltinger C, Sulman E, Marshall HN, Fujimori M, Kaufman BA, Biegel JA, Allen C, Hilliard C: A region of consistent deletion in neuroblastoma maps within human chromosome 1p36.2-36.3. Proc Natl Acad Sci U S A. 1995, 92: 5520-5524.
Bauer A, Savelyeva L, Claas A, Praml C, Berthold F, Schwab M: Smallest region of overlapping deletion in 1p36 in human neuroblastoma: a 1 Mbp cosmid and PAC contig. Genes Chromosomes Cancer. 2001, 31: 228-239.
Ejeskär K, Sjöberg RM, Abel F, Kogner P, Ambros PF, Martinsson T: Fine mapping of a tumour suppressor candidate gene region in 1p36.2-3, commonly deleted in neuroblastomas and germ cell tumours. Med Pediatr Oncol. 2001, 36: 61-66.
Spieker N, Beitsma M, Van Sluis P, Chan A, Caron H, Versteeg R: Three chromosomal rearrangements in neuroblastoma cluster within a 300-kb region on 1p36.1. Genes Chromosomes Cancer. 2001, 31: 172-181.
Krona C, Ejeskär K, Abel F, Kogner P, Bjelke J, Björk E, Sjöberg RM, Martinsson T: Screening for gene mutations in a 500 kb neuroblastoma tumor suppressor candidate region in chromosome 1p; mutation and stage-specific expression in UBE4B/UFD2. Oncogene. 2003, 22: 2343-2351.
Ejeskär K, Abel F, Sjöberg RM, Bäckström J, Kogner P, Martinsson T: Fine mapping of the human preprocortistatin gene (CORT) to neuroblastoma consensus deletion region 1p36.3-->p36.2, but absence of mutations in primary tumors. Cytogenet Cell Genet. 2000, 89: 62-66.
Abel F, Sjöberg RM, Ejeskär K, Krona C, Martinsson T: Analyses of apoptotic regulators CASP9 and DFFA at 1P36.2, reveal rare allele variants in human neuroblastoma tumours. Br J Cancer. 2002, 86: 596-604.
Krona C, Ejeskär K, Carén H, Abel F, Sjöberg RM, Martinsson T: A novel 1p36.2 located gene, APITD1, with tumour-suppressive properties and a putative p53-binding domain, shows low expression in neuroblastoma tumours. Br J Cancer. 2004, 91: 1119-1130.
Abel F, Sjöberg RM, Krona C, Nilsson S, Martinsson T: Mutations in the N-terminal domain of DFF45 in a primary germ cell tumor and in neuroblastoma tumors. Int J Oncol. 2004, 25: 1297-1302.
Plass C: Cancer epigenomics. Hum Mol Genet. 2002, 11: 2479-2488.
Herman JG, Baylin SB: Gene silencing in cancer in association with promoter hypermethylation. N Engl J Med. 2003, 349: 2042-2054.
Costello JF, Plass C: Methylation matters. J Med Genet. 2001, 38: 285-303.
Jones PA, Laird PW: Cancer epigenetics comes of age. Nat Genet. 1999, 21: 163-167.
Tycko B: Epigenetic gene silencing in cancer. J Clin Invest. 2000, 105: 401-407.
Esteller M: CpG island hypermethylation and tumor suppressor genes: a booming present, a brighter future. Oncogene. 2002, 21: 5427-5440.
Clark SJ, Harrison J, Paul CL, Frommer M: High sensitivity mapping of methylated cytosines. Nucleic Acids Res. 1994, 22: 2990-2997.
Herman JG, Graff JR, Myohanen S, Nelkin BD, Baylin SB: Methylation-specific PCR: a novel PCR assay for methylation status of CpG islands. Proc Natl Acad Sci U S A. 1996, 93: 9821-9826.
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.
Harada K, Toyooka S, Maitra A, Maruyama R, Toyooka KO, Timmons CF, Tomlinson GE, Mastrangelo D, Hay RJ, Minna JD, Gazdar AF: Aberrant promoter methylation and silencing of the RASSF1A gene in pediatric tumors and cell lines. Oncogene. 2002, 21: 4345-4349.
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.
Agathanggelou A, Dallol A, Zochbauer-Muller S, Morrissey C, Honorio S, Hesson L, Martinsson T, Fong KM, Kuo MJ, Yuen PW, Maher ER, Minna JD, Latif F: Epigenetic inactivation of the candidate 3p21.3 suppressor gene BLU in human cancers. Oncogene. 2003, 22: 1580-1588.
van Noesel MM, van Bezouw S, Voute PA, Herman JG, Pieters R, Versteeg R: Clustering of hypermethylated genes in neuroblastoma. Genes Chromosomes Cancer. 2003, 38: 226-233.
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: 1208-1219.
Ohira M, Kageyama H, Mihara M, Furuta S, Machida T, Shishikura T, Takayasu H, Islam A, Nakamura Y, Takahashi M, Tomioka N, Sakiyama S, Kaneko Y, Toyoda A, Hattori M, Sakaki Y, Ohki M, Horii A, Soeda E, Inazawa J, Seki N, Kuma H, Nozawa I, Nakagawara A: Identification and characterization of a 500-kb homozygously deleted region at 1p36.2-p36.3 in a neuroblastoma cell line. Oncogene. 2000, 19: 4302-4307.
Paulin R, Grigg GW, Davey MW, Piper AA: Urea improves efficiency of bisulphite-mediated sequencing of 5'-methylcytosine in genomic DNA. Nucleic Acids Res. 1998, 26: 5009-5010.
Li LC, Dahiya R: MethPrimer: designing primers for methylation PCRs. Bioinformatics. 2002, 18: 1427-1431.
This work was supported by grants from the Swedish Cancer Society, the Children's Cancer Foundation, the King Gustav V Jubilee Clinic Cancer Research Foundation, the Assar Gabrielsson Foundation, the Wilhelm and Martina Lundgren Research Foundation and the Sahlgrenska University Foundation. We would like to thank the Swegene Gothenburg Genomics resource unit for access to the ABI 7900HT Sequence Detection System and the ABI 3730 Genetic Analyzer.
HC participated in the design of the study, carried out the methylation analysis and drafted the manuscript. SF, LH and FL participated in parts of the methylation study. KE and CK carried out the TaqMan runs. KE also contributed to drafting the manuscript. RMS participated in the cloning. TM coordinated the study. All authors read and approved the final manuscript.