
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Imaging, Diagnosis, Prognosis |
Authors' Affiliations: Departments of 1 Tumour Genetics B030, 2 Molecular Genetics B060, 3 Theoretical Bioinformatics B080, and 4 Central Unit Biostatistics, Deutsches Krebsforschungszentrum, Heidelberg; 5 Department of Pediatric Oncology, University Children's Hospital, Cologne, Germany; and 6 Pediatric Oncology Branch, Oncogenomics Section, National Cancer Institute, Gaithersburg, Maryland
Requests for reprints: Kai-Oliver Henrich, Department of Tumour Genetics B030, Deutsches Krebsforschungszentrum, Im Neuenheimer Feld 280, 69120 Heidelberg, Germany. Phone: 49-6221-423220; Fax: 49-6221-423277; E-mail: k.henrich{at}dkfz.de.
| Abstract |
|---|
|
|
|---|
Experimental Design: Candidate genes localized within the deleted region were identified by sequence data analysis. Their expression was assessed in a cohort of 49 primary neuroblastomas using cDNA microarray analysis. Gene expression patterns associated with known prognostic markers and patient outcome were further evaluated by quantitative real-time reverse transcription-PCR in a cohort of 102 neuroblastomas.
Results: The commonly deleted region spans 261 kb and encompasses two genes, FLJ10737 and CAMTA1. We found no evidence for an association of FLJ10737 expression with established prognostic variables or outcome. In contrast, low CAMTA1 expression characterized tumors with 1p deletion, MYCN amplification, and advanced tumor stages 3 and 4. Moreover, low CAMTA1 expression was significantly associated with poor outcome (P < 0.001). In multivariate analysis of event-free survival, the prognostic information of low CAMTA1 expression was independent of 1p status, MYCN status, tumor stage, and age of the patient at diagnosis (hazard ratio, 3.52; 95% confidence interval, 1.21-10.28; P = 0.02).
Conclusions: Our data suggest that assessment of CAMTA1 expression may improve the prognostic models for neuroblastoma and that it will be important to define the biological function of CAMTA1 in this disease.
Deletions within 1p occur in
30% of all neuroblastomas, constitutional translocations and deletions have been reported as well (4, 5). Introduction of 1p chromosomal material into neuroblastoma cells resulted in reduced tumorigenicity (6), providing functional evidence for the role of 1p in neuroblastoma. Furthermore, deletion of 1p36 is associated with adverse prognostic markers such as amplified MYCN and has been identified as an independent predictor of decreased event-free survival (7). Thus, it is widely assumed that distal 1p harbors neuroblastoma-related genetic information.
In pursuit of defining the identity of the putative 1p neuroblastoma gene(s), substantial efforts have been undertaken to confirm and refine a 1p smallest region of overlapping deletions (SRO) in neuroblastoma. Early analyses (813) led to a preliminary composite SRO of about 5 Mb between the first retained distal marker D1S47 (9) and the first retained proximal marker D1S244 (12, 14). Subsequent studies have attempted to reduce the genomic complexity of the deleted region. A promising step in this direction was suggested by the identification of a homozygous deletion in a neuroblastoma cell line spanning
500 kb at D1S244 (15). However, thus far, extensive analysis of the six genes mapping to this region gave no solid evidence for any of them functioning as a neuroblastoma suppressor gene (1522).
Recent studies have made it possible to considerably narrow down a smallest region of overlapping heterozygous deletions in 1p36.3 mapping distal to the homozygous deletion at D1S244. The borders of a SRO that spans a region between D1S2660 (distal) and D1S214 (proximal) have been identified (2325). This region partially overlaps with a 754 kb SRO defined by D1S2731 (distal) and D1S2666 (proximal; ref. 26). Collectively, these results define a 261 kb composite SRO between D1S2731 and D1S214. We examined the expression of candidate genes mapping to this region in neuroblastomas to determine whether aberrant expression of these genes might play a role in neuroblastoma development. Potential prognostic information in these expression profiles was investigated by correlating them to clinicobiological variables and outcome in a large representative cohort of patients with neuroblastoma.
| Materials and Methods |
|---|
|
|
|---|
|
cDNA microarray analysis
Labeling of cDNA, fabrication of microarrays, hybridization, and image analysis were done as described previously (31). Briefly, mRNA was amplified with one round of a modified Eberwine RNA amplification procedure (32), and used for indirect fluorescent-labeling cDNA synthesis (33). Microarrays were fabricated printing a total of 42,578 cDNA clones (libraries purchased from Research Genetics, Huntsville, AL) on glass slides using a BioRobotics MicroGrid II spotter (Harvard Bioscience, Holliston, MA). Hybridization and washing of microarrays were done as described elsewhere (33). Images were acquired by an Agilent DNA Microarray scanner (Agilent) and analyzed using the Microarray Suite program (IPLab, Scanalytics, Fairfax, VA). Raw microarray data were normalized using the variance stabilization method (34), which has proven to be robust and reliable for normalization of microarray data (35). Briefly, values from each probe were transformed using the variance stabilizing asinh function, which is asymptotically equal to log(2x). For each chip, two parameters representing an offset and a scaling factor were estimated by the maximum likelihood method. Values were obtained after normalization on a logarithmic scale. Calculations were done with R (36) and Bioconductor (37) software using the package vsn (version 1.6.3). In order to check performance of the normalization and assess the quality of hybridizations, all data were subjected to quality assessment methods implemented in the Bioconductor package arrayMagic version 1.5.8 (38).
QPCR
Single-stranded cDNA was generated from total RNA of primary neuroblastomas using the Superscript II First-Strand Synthesis System (Invitrogen) according to the manufacturer's directions. QPCR was done on an ABI PRISM 7700 Sequence Detection System (Applied Biosystems, Foster City, CA) with SYBR Green chemistry using the standard curve method (user bulletin no. 2, ABI PRISM 7700 SDS). To prevent amplification from contaminating genomic DNA, primer sequences were selected allowing intron-spanning amplification. PCR reactions were run as duplicates for each sample and as triplicates for determination of standard curves. For normalization, the expression level of the target gene was divided by the geometric mean of expression levels of the housekeeping genes hypoxanthine phosphoribosyltransferase 1 (HPRT1) and succinate dehydrogenase complex, subunit A (SDHA; ref. 39). Primer sequences for CAMTA1 were: 5'-AGTGCAGAAAATGAAGAATGCG-3' (forward) and 5'-CAAAATTCTCCTGCTTGATTCG-3' (reverse).
Statistics
Analysis of cDNA microarray expression data. To test the association between candidate gene expression and established prognostic factors, the Wilcoxon rank sum test was used. Univariate survival analysis was done using Cox proportional hazards regression. For each candidate gene, the estimated hazard ratio and its 95% confidence interval (CI) were determined for an increase of expression levels from the lower to the upper quartile to obtain comparable results for all genes analyzed. A generalized likelihood ratio test was done to determine whether these changes in gene expression in relation to event-free survival and overall survival were statistically significant.
Analysis of QPCR expression data. To test the association between candidate gene expression and established prognostic factors, the Wilcoxon rank sum test was used. Correlation between the corresponding cDNA microarray and QPCR expression measurements was estimated by Pearson's product moment correlation coefficient. To identify a model describing the relationship between survival and candidate gene expression, the functional form of this relationship was tested by multivariable fractional polynomials (40) and maximally selected log-rank statistics (41). The resulting model was applied in further survival analyses. Multivariate Cox regression was used to investigate the prognostic power of candidate gene expression adjusting for established prognostic variables. To correct for overestimation of the hazard ratio estimate due to selection of the functional form of candidate gene expression variables, shrinkage of the parameter estimate was applied (42). Univariate survival analysis for established prognostic variables was done by means of the log-rank test. The Kaplan-Meier method was used to estimate event-free survival and overall survival curves.
Event-free survival was measured from date of diagnosis until occurrence of disease progression, relapse, or death due to neuroblastoma. Event-free survival times of patients who have not experienced an event until last follow-up were censored. Overall survival was measured from date of diagnosis until death due to neuroblastoma. Overall survival times of patients that were alive at the last follow-up were censored. An effect was considered as statistically significant if the P value of its corresponding test statistic was
5%. The statistical analyses were done using the software package R, version 2.0.1 (36).
| Results |
|---|
|
|
|---|
|
CAMTA1 and FLJ10737 expression as determined by microarray analysis. Expression patterns of CAMTA1 and FLJ10737 were derived from a cDNA microarray analysis7 of 49 neuroblastomas (Table 1) and correlated to clinicobiological data and outcome of the respective patients.
We found no significant association of FLJ10737 expression with MYCN status (Wilcoxon rank sum test, P = 0.61), stage (P = 0.23), or 1p status (P = 0.57). To determine whether FLJ10737 expression levels were related to differences in survival time distribution, hazard ratios were estimated by Cox proportional hazards regression comparing the upper to the lower quartile of the expression range. Differences in FLJ10737 expression were not statistically significant in relation to event-free survival (hazard ratio, lower quartile versus upper quartile: 1.25; 95% CI, 0.73-2.17; P = 0.41) or overall survival (hazard ratio, lower quartile versus upper quartile: 1.37; 95% CI, 0.61-3.03; P = 0.45).
CAMTA1 expression was significantly lower in MYCN-amplified tumors (P < 0.001), advanced stages (stages 3 and 4, P = 0.005), and tumors with 1p deletion (P < 0.001; Fig. 2A). Furthermore, low CAMTA1 expression was significantly associated with decreased event-free survival (hazard ratio, lower quartile versus upper quartile: 2.94; 95% CI, 1.64-5.26; P < 0.001) and decreased overall survival (hazard ratio, lower quartile versus upper quartile: 5.88; 95% CI, 2.13-16.67; P < 0.001).
|
The search for a model describing the relationship between survival and CAMTA1 expression using fractional polynomials and maximally selected log-rank statistics identified a cut point model to be most suitable: a CAMTA1 QPCR cutoff value of 4.44 separated two patient subgroups with significantly different outcome (maximally selected log-rank statistics, P < 0.001). Applying this CAMTA1 expression cutoff to Kaplan-Meier survival curve estimation revealed decreased survival probability and increased recurrence rate for patients with tumors expressing low levels of CAMTA1 (Fig. 3).
|
1 year (P = 0.002). Cox proportional hazards models were built based on 94 patients with complete data on these prognostic variables. Low CAMTA1 expression was selected as the only independent predictor for decreased event-free survival (Table 2). To further illustrate the prognostic independence of CAMTA1 expression from 1p deletion, multivariate survival analysis was done in the subgroup of 71 tumors with normal 1p status. Even within this cohort, CAMTA1 expression, together with tumor stage, emerged as an independent prognostic factor (Table 3).
|
|
| Discussion |
|---|
|
|
|---|
Expression analysis of 49 neuroblastomas by cDNA microarray gave no evidence for an association of FLJ10737 expression with known variables of poor outcome like MYCN amplification, advanced stages (3 and 4), or 1p deletion. Furthermore, we found no significant association of FLJ10737 expression with outcome. An absence of a correlation between FLJ10737 expression and 1p status might be explained by the existence of homeostatic mechanisms that compensate for copy number changes of this gene. Together, these results do not provide evidence for a contribution of deregulated FLJ10737 expression to neuroblastoma pathogenesis. However, sequence analysis will be necessary to clarify whether FLJ10737 mutations or sequence variants exist that exert influence on the disease.
Expression analysis of 49 neuroblastomas by cDNA microarray revealed significantly lower CAMTA1 expression in MYCN-amplified tumors, tumors of stages 3 and 4, and tumors with 1p deletion. The prognostic relevance of CAMTA1 expression indicated by the association with these variables of poor prognosis was confirmed by univariate analysis for both overall survival and event-free survival. The use of QPCR corroborated these results in a cohort of 102 patients.
However, these findings alone were not sufficient to indicate that CAMTA1 could add supplementary information to the existing risk stratification. Approximately 15% of 1p35-36 encoded genes are expressed at significantly lower levels in 1p-deleted compared with 1p-normal neuroblastomas (43). Low expression of such a gene, even if not involved in neuroblastoma development, is likely to be associated with poor outcome as it can act as a surrogate marker for 1p deletion [which itself might be a surrogate marker for inactivation or repression of other 1p gene(s)]. However, CAMTA1 expression contained prognostic information that was not covered by the variable "1p status." Multivariate survival analysis identified CAMTA1 expression as a predictor variable that was independent not only of 1p status but also of MYCN status, stage, and age at diagnosis. The prognostic independence of CAMTA1 expression from 1p deletion was further shown by multivariate survival analysis within the subgroup of 1p normal neuroblastomas.
All of the analyzed established prognostic variables (MYCN status, 1p status, stage, and age at diagnosis) were predictive of outcome in univariate survival analysis. This indicates that our group of patients was representative of an unselected neuroblastoma cohort and that the supplementary prognostic value of CAMTA1 expression reflects that in the general neuroblastoma population. Thus, the measurement of this variable should allow an additional biological stratification of tumors and help to assign patients to the appropriate therapy. The fact that the differential CAMTA1 expression with respect to survival probability could be detected by two independent methods, QPCR and cDNA microarray analysis, indicates the robustness of this prognostic marker. The prognostic value of CAMTA1 expression is further supported by a genome-wide microarray expression study that identified CAMTA1 as one of the 47 top-ranked differentially expressed transcripts separating progressive from regressive neuroblastoma phenotypes.7
CAMTA1 is a member of a recently characterized protein family designated as calmodulin-binding transcription activators (CAMTA; ref. 44). The physiologic roles of these proteins remain to be elucidated but it has been shown for plant CAMTAs that they (a) bind double-stranded DNA, (b) are predominantly localized in the nucleus, (c) bind calmodulin, and (d) could activate transcription. The highly conserved domain structure within this protein family suggests similar functions for CAMTA1 and its role as transcription activator could be confirmed in a yeast reporter system (44). These data indicate that CAMTA1 is a transcription factor which possibly responds to calcium signaling by direct binding of calmodulin. Such activities would well figure into pathways of tumorigenesis as many calcium-binding proteins are involved in cell cycle control and/or apoptotic cell death, and modulation of calmodulin levels has been recognized as a feature of cell transformation and of malignant cells (45, 46).
The CAMTA1 expression pattern found in neuroblastomas lets us speculate that down-regulation of CAMTA1 mediates a selective advantage in the course of neuroblastoma development. A variety of mechanisms could be responsible for low expression of CAMTA1 in unfavorable tumors.
A high proportion of this tumor entity displays 1p deletion. The significantly lower CAMTA1 levels in 1p-deleted tumors suggests a decrease of the gene's expression due to a reduction of its copy number. There is a growing body of literature on genes whose dosage reduction accounts for phenotypic changes that contribute to tumor development (47, 48). Whether reduced CAMTA1 gene dosage as a consequence of 1p deletion brings about such functional haploinsufficiency needs to be clarified. Low CAMTA1 expression in the absence of 1p deletion calls for additional regulation mechanisms in the respective tumors. Deregulation of upstream elements and transcription factors might account for differential CAMTA1 expression independent of 1p deletion. CAMTA1 levels were significantly lower in MYCN-amplified neuroblastomas and a direct or indirect regulation of CAMTA1 by MYCN might be a possible explanation for this observation. However, monitoring CAMTA1 levels in the MYCN-inducible neuroblastoma cell line Tet21N gave no evidence supporting this hypothesis (data not shown). In many tumor types, methylation of cytosine residues in gene-associated CpG islands is a common mechanism suppressing transcription of growth-regulating genes (49, 50). Whether transcription of CAMTA1 is suppressed by methylation remains to be investigated.
A potential role for CAMTA1 in tumor development is strongly supported by the findings of a recent loss of heterozygosity study on 205 gliomas: a 1p minimal deleted region was identified that spans a region of 150 kb between D1S2694 and D1S2666. Thus far, the only gene mapped to this region is CAMTA1 (51). Further support comes from a query of the microarray database Oncomine (52). A CAMTA1-specific search highlighted a study classifying subgroups of melanoma by genome-wide expression analysis (53). CAMTA1 expression was significantly reduced in melanoma of soft parts, a rare and aggressive melanoma subtype. This finding is intriguing considering that melanomas, as neuroblastomas, arise from neural crest cells, and suggests the need to investigate the prognostic relevance of CAMTA1 expression in this tumor.
Taken together, our data show that the expression of CAMTA1, a gene mapping to a 1p36 region commonly deleted in neuroblastoma, can represent a powerful prognostic variable that may complement the predictive value of established risk factors in neuroblastoma. We suggest that CAMTA1 reduction might contribute to the effect of 1p36 deletion on event-free survival (7). An understanding of the function and regulation of this gene and its product should give further insight into the biological features of neuroblastoma. Sequence analysis should address the question of whether CAMTA1 mutations or sequence variants exist that are associated with neuroblastoma.
| Acknowledgments |
|---|
| Footnotes |
|---|
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
Note: R. Wiedemeyer is presently at the Department of Medical Oncology, Dana-Farber Cancer Institute, Boston, MA.
7 Westermann F, Wei JS, Lutz W, et al. Tumor suppressor functions of the Rb-Skp2 pathway are disabled in high-risk but functional in low-risk neuroblastomas, manuscript in preparation. ![]()
Received 7/ 5/05; revised 9/13/05; accepted 10/11/05.
| References |
|---|
|
|
|---|
p36.2, but absence of mutations in primary tumors. Cytogenet Cell Genet 2000;89:626.[CrossRef][Medline]
gene located at 1p36.2 in neuroblastoma. Int J Oncol 2003;23:73744.[Medline]This article has been cited by other articles:
![]() |
F. Westermann, K.-O. Henrich, J. S. Wei, W. Lutz, M. Fischer, R. Konig, R. Wiedemeyer, V. Ehemann, B. Brors, K. Ernestus, et al. High Skp2 Expression Characterizes High-Risk Neuroblastomas Independent of MYCN Status Clin. Cancer Res., August 15, 2007; 13(16): 4695 - 4703. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. Gong, J. Han, K. Reddig, and H.-S. Li A Potential Dimerization Region of dCAMTA Is Critical for Termination of Fly Visual Response J. Biol. Chem., July 20, 2007; 282(29): 21253 - 21258. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. J. Huentelman, A. Papassotiropoulos, D. W. Craig, F. J. Hoerndli, J. V. Pearson, K.-D. Huynh, J. Corneveaux, J. Hanggi, C. R.A. Mondadori, A. Buchmann, et al. Calmodulin-binding transcription activator 1 (CAMTA1) alleles predispose human episodic memory performance Hum. Mol. Genet., June 15, 2007; 16(12): 1469 - 1477. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Schramm, J. Vandesompele, J. H. Schulte, S. Dreesmann, L. Kaderali, B. Brors, R. Eils, F. Speleman, and A. Eggert Translating Expression Profiling into a Clinically Feasible Test to Predict Neuroblastoma Outcome Clin. Cancer Res., March 1, 2007; 13(5): 1459 - 1465. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Asgharzadeh, R. Pique-Regi, R. Sposto, H. Wang, Y. Yang, H. Shimada, K. Matthay, J. Buckley, A. Ortega, and R. C. Seeger Prognostic Significance of Gene Expression Profiles of Metastatic Neuroblastomas Lacking MYCN Gene Amplification. J Natl Cancer Inst, September 6, 2006; 98(17): 1193 - 1203. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Cancer Research | Clinical Cancer Research |
| Cancer Epidemiology Biomarkers & Prevention | Molecular Cancer Therapeutics |
| Molecular Cancer Research | Cancer Prevention Research |
| Cancer Prevention Journals Portal | Cancer Reviews Online |
| Annual Meeting Education Book | Cell Growth & Differentiation |