| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Regular Articles |
1 Division of Molecular Genetics (B060) and 2 Central Unit Biostatistics (C060), Deutsches Krebsforschungszentrum, Heidelberg, Germany, and Departments of 3 Neuropathology, 4 Neurooncological Surgery, NN Burdenko Neurosurgical Institute, Moscow, Russia
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Most pediatric medulloblastomas have abnormal karyotypes, including aneuploidy, unbalanced translocations, and numerous structural rearrangements in the form of double-minute, ring, and marker chromosomes (4
, 5)
. In 3050% of patients, deletion or rearrangement of part of chromosome 17 have been described previously (6, 7, 8)
. In most cases, 17p is lost, and head-to-head apposition of both 17q arms occurs, which had been previously considered as isochromosome 17q (9)
. Medulloblastomas with loss of 17p have breakpoints either proximally to the TP53 locus at 17p11.2 or distal to TP53 at 17p13.3. In
10% of medulloblastomas, gene amplifications were identified: the common amplified regions include the oncogenes MYC (8q24) and MYCN (2p24.1; Ref. 8
, 10
). Although the clinical significance of the majority of these abnormalities remains unclear, recent studies show that amplifications as well as high transcriptional levels of MYC are associated with an unfavorable survival outcome (10
, 11)
.
Recently, new insight in the prediction of medulloblastoma outcome was obtained by gene expression profiling experiments. The expression of genes characteristic for cerebellar differentiation or coding for extracellular matrix proteins was correlated with favorable prognosis, whereas the expression of genes related to cell cycle control, multidrug resistance, and ribosomal biogenesis was associated with poor prognosis (12) . Moreover, metastatic disease was associated with expression profiles, including platelet-derived growth factor receptor and genes that are involved in the RAS/MAPK signaling pathway (13) .
In the present study, we performed an expression profiling analysis of 35 medulloblastoma samples to identify potential candidate genes that predict for overall survival (OS) of the patients. Additionally, we applied the approach of tumor tissue microarray (TMA) experiments (14) , a modern high-throughput technology that allowed screening of 180 clinically well-defined medulloblastoma tumor specimens on a large scale. Nine promising candidate genes were analyzed for changes on the protein expression level by immunohistochemistry (IHC). In addition, fluorescence in situ hybridization (FISH) was performed for the detection of additional gene copies of four of them.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Clinical information on the patients included date of birth, date of operation, extent of resection determined by postoperative neuroradiological reports (gross total resection, no visible tumor; near-total resection, >90% of tumor removed; subtotal resection, 5190% of tumor removed), irradiation or chemotherapy, follow-up time, and outcome. Metastatic stage was studied by craniospinal neuroimaging and liquor cytopathology. Approval to link laboratory data to clinical data was obtained by the Institutional Review Board.
Table 1
shows the characteristics of patients with medulloblastoma included in the cDNA microarray (n = 35) and TMA experiments (n = 180). There was an overlap of 30 samples, which were included in both studies. All patients in both series received craniospinal irradiation after surgery. The postoperative irradiation included 36 Gy to the craniospinal axis and a boost of 5356 Gy to the posterior fossa. Although all patients enrolled in the cDNA microarray study were additionally treated with adjuvant chemotherapy, this was only applied to 112 patients (62.2%) included in the TMA study. Chemotherapy was routinely administered to all patients who were treated after 1995. The chemotherapy regimen comprises lomustine, cisplatin, and vincristine as published previously (2
, 16)
.
|
cDNA Microarray Experiments.
Samples from each primary tumor were taken at the time of initial operation, immediately frozen in liquid nitrogen, and then continuously stored at 80°C until usage for cDNA microarray experiments. To confirm the presence of viable cellular tumor (
90% neoplastic cells), cryosections of each medulloblastoma sample (
5-µm thick) were stained with H&E before RNA extraction. Total RNA was isolated according to a protocol that applies both Trizol reagent (Invitrogen, Karlsruhe, Germany) and RNeasy Midi spin columns (Qiagen, Hilden, Germany), followed by Oligotex resin (Qiagen) purification of poly-A RNA (17, 18, 19)
. The integrity and purity of the total RNA was analyzed on a 1% agarose gel. Total RNA quality and concentration were determined by absorption spectroscopy between 220 and 320 nm (1 A260 = 40 µg/ml single-stranded RNA).
For gene expression experiments, cDNA microarrays containing replicate spots of 4211 different gene specific fragments representing 2600 different genes with relevance to mitosis, cell cycle control, oncogenesis, or apoptosis were processed as described previously (17, 18, 19) . The tumor RNA was cohybridized with commercially available universal human reference RNA (Stratagene, La Jolla, CA), derived from 10 human cancer cell lines. Approximately 1 µg of medulloblastoma and 1 µg of reference mRNA were labeled with Cy3 and Cy5, respectively, using the Omniscript Reverse Transcriptase kit (Qiagen) and hybridized with 10 µg of C0t1 DNA, 30 µg of bovine liver tRNA, and 10 µg of oligo-dT nucleotides in an automated hybridization chamber (GeneTac; Genomic Solutions, Ann Arbor, MI). For all samples, we performed color switch experiments, where the tumor and reference DNA were labeled via Cy3- and Cy5-dUTP, respectively, and vice versa. Data sets for spots not recognized by the GenePix Pro 4.0 analysis software (Axon Instruments, Union City, CA) were excluded from additional considerations. Additionally, all remaining data sets were ranked according to spot homogeneity (as assayed by the ratio of median and mean fluorescence intensities), spot intensity, and the SD of log ratios for replicate spots. Those data points, ranked among the lower 20%, based on the criteria just described, were removed from the data set. For each hybridization, fluorescence ratios (Cy5/Cy3) were normalized by variance stabilization (20) . To combine experiments with switched dye labeling, the ratios of one experiment were inverted and averaged with the corresponding spots on the second array. The raw data of the microarray experiments are available online.5
Generation of TMAs.
As previously described (21)
, H&E-stained sections were made from each paraffin block to define representative tumor regions. From these regions, two tissue cylinders with a diameter of 0.6 mm were obtained and arrayed into a recipient block using the tissue chip microarrayer (Beecher Instruments, Silver Spring, MD). The recipient block was subsequently cut into sections of 5-µm on Histobond-silanized glass slides (Marienfeld, Lauda-Koenigshofen, Germany), which were treated with 3-aminopropylmethoxysilane (Fluka, Seelze, Germany) before cutting to support adhesion of the tissue samples.
IHC.
Paraffin sections were deparaffinized and rehydrated with graded ethanol. Antigen retrieval was performed by heating for 10 min in a microwave either in 10 mM sodium citrate or Tris-EDTA buffer. Endogenous peroxidase was inactivated by incubating the TMAs in 3% hydrogen peroxide. After blocking with normal serum, the sections were incubated with the primary antibodies at 4°C overnight. Antibodies against the following antigens were used: cMYC, clone 9-E11 (1:40): Novocastra (Newcastle upon Tyne, United Kingdom); cyclin B1, clone GNS1 (1:100), DAXX, clone H-7 (1:100), NRAS, clone F-155 (1:150), stathmin 1, clone FL149 (1:250), STK15, clone N20 (1:250): Santa Cruz Biotechnology (Santa Cruz, CA); cyclin D1, clone DSC6 (1:40), Ki-67, clone MIB-1 (1:100): DAKO (Glostrup, Denmark); and p53, clone BP53-11 (1:1000): Progene (Darra Queensland, Australia). For detection, TMAs were incubated for 30 min with biotinylated secondary antibody, and for 30 min with ABC Reagent (Vector Laboratories, Burlingame, CA). The slides were stained with 3,3'-diaminobenzidine (Linaris, Wertheim-Bettingen, Germany), counterstained with hematoxylin, and mounted. For negative controls, primary antibodies were substituted with commercially produced control reagents according to the manufacturers instructions.
An assessment of IHC was performed for each tissue core in a semiquantitative manner as a percentage of stained cells or nuclei/specimen independently by two different experimenters who were blinded to clinical data. Diffuse cytoplasmic and/or membranous immunostaining of DAXX, NRAS, stathmin 1, and STK15 as well as nuclear expression patterns of cyclin B1, cyclin D1, cMYC, p53, and Ki-67 of >10% of cells was considered as overexpression. Therefore, the analysis of our IHC experiments was only performed by a positive versus negative evaluation of stained cores, which guaranteed a high reproducibility of our results.
FISH.
Two-color FISH to the sections of the TMAs was performed using Spectrum Orange-labeled probes for CCND1, MYC [both Vysis (Downers Grove, IL)], NRAS (PAC RP51000E10), and STK15 [PAC RP5-1167H4, both (RZPD German Resource Center, Berlin, Germany)], and Spectrum Green-labeled centromeric probes of the chromosomes 1, 8, 11, and 20 (Vysis). Pretreatment of slides, hybridization, posthybridization processing, and signal detection were performed as described previously (22)
.
Medulloblastoma samples showing sufficient FISH efficiency (> 90% nuclei with signals) were evaluated, with 100 intact nuclei scored for the number of signals. Tumor signals were scored as gains when at least 10% of cells showed three or more signals of the oncogene probe. Three to seven signals/cell were considered as low level gains, whereas eight or more signals as well as tight signal clusters were counted as high level gains. The signals of the centromeric probes were used to control adequate hybridization and to exclude artifacts.
Statistical Analysis.
Primary end point of the study was OS. To identify single genes which expression levels were possibly predictive for OS, we applied genewise proportional hazards models. This approach takes into account that all patients have a different follow-up time. For each gene, the estimated hazard ratio and its 95% confidence interval 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 performed for each gene to determine whether these changes in gene expression in relation to OS were of statistical significance. We used a hazard ratio of >2.5 (i.e., 5:2) as cutoff level to identify a set of genes associated with unfavorable OS, whereas a hazard ratio of <0.4 (i.e., 2:5) was used to obtain genes predicting for favorable survival.
Survival probabilities were estimated by the method described by Kaplan and Meier. All genes analyzed by TMA that showed a statistically significant effect in the log-rank test (cyclin D1, stathmin 1, and STK15) were also tested in a multivariate proportional hazards regression model together with clinically defined variables such as age, gender, and metastatic stage. An effect was considered statistically significant if the P was
0.05. The statistical analyses were performed using the software package SPSS 11.5 (SPSS, Inc., Chicago, IL) as well as the software package R, version 1.7 (23)
.
| RESULTS |
|---|
|
|
|---|
|
|
|
|
|
|
|
|
|
|
| DISCUSSION |
|---|
|
|
|---|
In detail, we identified 54 and 27 genes with expression levels predicting for unfavorable and favorable prognosis in medulloblastoma, respectively. The genes predicting for poor outcome comprise those that are involved in cell cycle activation and mitosis, DNA metabolism, or code for viral oncogenes. Some of these genes (MYC, SOX4, and NRAS) have been previously reported to be involved in medulloblastoma tumorigenesis and associated with unfavorable outcome (10 , 11 , 13 , 25) . Concordantly, the potential role of platelet-derived growth factor receptor and downstream activation of the RAS/MAPK signaling pathway for metastasis in medulloblastoma was outlined by MacDonald et al. (13) . The list of unfavorable genes also includes MGMT, which codes for the DNA repair protein O6-methylguanine-DNA methyltransferase. Bacolod et al. (26) found that MGMT protein expression confers resistance to alkylating drugs in medulloblastoma-derived cell lines. The MGMT-mediated resistance to alkylating agents might be of importance because these drugs (e.g., lomustine) have been included in treatment protocols for medulloblastoma patients (2 , 16) .
Our set of genes, predicting for favorable prognosis, comprises TMSB4X, of which, expression is closely related to cerebellar differentiation (27) and PRDX1, which has been recently identified as candidate for tumor suppression (28) . Furthermore, we identified genes that have apoptosis-related functions (CASP10 and TRADD), are involved in inflammatory and immune response processes (IRF3 and IL11RA) or code for zinc finger proteins (ZNF9 and ZNF162).
Gains of chromosomes 7 or 17q were described in up to 60% of all medulloblastoma samples (5 , 8 , 24) . Interestingly, none of the genes associated with good prognosis were located on these chromosomal regions. In contrast, the set of genes predicting for adverse outcome included four genes located on chromosome 7 (RHEB2, YWHAG, PMS2, and PMSB1) and four genes on 17q (ARHGDIA, NSF, TOP2A, and TUBG1). This finding suggests that a cooperative selective growth advantage for medulloblastoma cases with alterations on chromosome 7 and 17q may exist, thus explaining the high frequency of aberrations involving both chromosomes along with the finding that high expression levels of genes located on both chromosomal regions correlated with an unfavorable survival.
In line with the results obtained by cDNA microarray experiments, a high expression of cyclin D1, STK15, and stathmin 1 on the protein level correlated with adverse prognosis in medulloblastoma. In contrast, no correlation between gene expression and IHC was found for cMYC, cyclin B1, DAXX, Ki-67, NRAS, and p53. To some extent, the discrepancies found for these genes on the protein and mRNA level might be related to different follow-up times of the patients included in the TMA and cDNA microarray studies (median follow-up time of 34 versus 24 months). On the other hand, the lack in correlation might be related to the short half-life time of some proteins. Formally, immunohistochemical heterogeneity of the tumor samples cannot be excluded, thus explaining the limitations of the TMA technology. Previous studies showed, however, that intratumor heterogeneity does not significantly affect the ability to detect clinicopathological correlations on TMAs, at least when two (as in our study) or more cores of each tumor are present on the TMAs (29) . To see an influence on OS, the total amount of the accumulated proteins probably needs to exceed a certain cellular cutoff level. Therefore, additional studies are needed to determine robust standardized cutoff values for different antibodies to define protein overexpression in medulloblastoma.
Recently, cyclin D1 was identified as an important mediator of Sonic hedgehog-induced proliferation in cerebellar granule cell precursors, which are a common target of transformation in medulloblastoma (30) . Because 2030% of sporadic medulloblastomas harbor mutation in patched and other elements of the Sonic hedgehog pathway (31) , cyclin D1 appears particularly intriguing with respect to medulloblastoma tumorigenesis. Both STK15 and stathmin 1 play an important role in the adequate segregation of chromosomes during mitosis. Mitotic defects might be a crucial event in medulloblastoma pathogenesis because aneuploidy is a common finding (5 , 24) . Stathmin 1 is a cytosolic phosphoprotein that is involved in cellular proliferation because its major function is to promote depolymerization of spindle microtubules. In Wilms tumors, stathmin 1 was found highly expressed in more advanced stages of disease (32) . In addition, high protein expression of stathmin 1 correlated with poor prognosis in breast cancer (33) . Several studies showed that stathmin 1 is a potential therapeutic target because a stathmin 1-transfected lung cancer cell line increased the sensitivity to Vinca alkaloids (34) , and antisense inhibition of stathmin 1 expression showed a synergistic apoptotic effect in combination with paclitaxel treatment (35) .
The STK15 gene is localized on chromosome 20q13, a region amplified in a variety of human malignancies, including medulloblastoma (5) . Additional gene copies of STK15 were described in primary tumors of the breast (36) , cervix (37) , bladder (38) , colon (39) , and brain (40) and reported to be associated with aneuploidy and tumor transformation. Because recent studies showed that overexpression of STK15 results in aberration of chromosomes and transformation of the cells, STK15 might contribute to an increasing karyotypic instability and tumor progression in medulloblastoma (36) .
When the prognostic significance of cyclin D1, STK15, and stathmin 1 protein expression was compared with clinically defined parameters, STK15 was identified as the only independent predictor for survival in medulloblastoma. Notably, STK15 expression was statistically more significant than the metastatic stage at diagnosis. This is in line with recent studies, showing that outcome predictors based on molecular markers are able to exceed well defined clinical parameters (3 , 11 , 12) .
In the current study, high level gains of MYC as well as low level gains of STK15 predicted for an unfavorable outcome of medulloblastoma patients. The frequency of MYC amplifications in the present study (6.7%) as well as their significance for survival is in line with previous studies (8 , 10 , 24) , whereas additional gene copies of STK15, which we found in approximately one-third of medulloblastoma samples have not been described previously. In addition, low level gains of STK15 were associated with a higher expression (1.7-fold) of this gene at the transcriptional level, suggesting a gene-dosis effect. Notably, increased gene copy numbers of STK15, as well as its high expression on the mRNA and protein level, correlated with an adverse outcome in medulloblastoma. Thus, STK15 may serve as a potent molecular target for additional medulloblastoma treatment because the inhibition of STK15 gene expression by antisense oligonucleotides resulted in the arrest of cell growth and increased apoptosis in cancer cell lines (41) .
In summary, the combination of gene expression profiling with TMA technology seems to be a powerful approach for the identification of potent candidate genes predicting for survival in medulloblastoma. Of the results highlighted by the various data analysis procedures, genes associated with cell cycle control and mitosis like cyclin D1, STK15, and stathmin 1 appear particularly intriguing with respect to medulloblastoma pathomechanism and might be promising targets to improve the cure rates in this disease.
| 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: K. Neben and A. Korshunov contributed equally to this work.
Requests for reprints: Peter Lichter, Deutsches Krebsforschungszentrum, Division of Molecular Genetics (B060), Im Neuenheimer Feld 280, D-69120 Heidelberg, Germany. Phone: 49-6221-42-4619; Fax: 49-6221-42-4639; E-mail: m.macleod{at}dkfz.de
5 Internet address: http://www.dkfz.de/kompl_genome/Other/medulloblastoma.zip ![]()
Received 12/18/03. Revised 2/ 6/04. Accepted 2/12/04.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
T. de Haas, N. Hasselt, D. Troost, H. Caron, M. Popovic, L. Zadravec-Zaletel, W. Grajkowska, M. Perek, M.-C. Osterheld, D. Ellison, et al. Molecular Risk Stratification of Medulloblastoma Patients Based on Immunohistochemical Analysis of MYC, LDHB, and CCNB1 Expression Clin. Cancer Res., July 1, 2008; 14(13): 4154 - 4160. [Abstract] [Full Text] [PDF] |
||||
![]() |
O. Gautschi, J. Heighway, P. C. Mack, P. R. Purnell, P. N. Lara Jr., and D. R. Gandara Aurora Kinases as Anticancer Drug Targets Clin. Cancer Res., March 15, 2008; 14(6): 1639 - 1648. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. Ogawa, K. Takenaka, H. Katakura, M. Adachi, Y. Otake, Y. Toda, H. Kotani, T. Manabe, H. Wada, and F. Tanaka Perimembrane Aurora-A Expression is a Significant Prognostic Factor in Correlation with Proliferative Activity in Non-Small-Cell Lung Cancer (NSCLC) Ann. Surg. Oncol., February 1, 2008; 15(2): 547 - 554. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. R. Pollack A Perspective on DNA Microarrays in Pathology Research and Practice Am. J. Pathol., August 1, 2007; 171(2): 375 - 385. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Fu, M. Bian, Q. Jiang, and C. Zhang Roles of Aurora Kinases in Mitosis and Tumorigenesis Mol. Cancer Res., January 1, 2007; 5(1): 1 - 10. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Reiter, P. Gais, U. Jutting, M. K. Steuer-Vogt, A. Pickhard, K. Bink, S. Rauser, S. Lassmann, H. Hofler, M. Werner, et al. Aurora Kinase A Messenger RNA Overexpression Is Correlated with Tumor Progression and Shortened Survival in Head and Neck Squamous Cell Carcinoma Clin. Cancer Res., September 1, 2006; 12(17): 5136 - 5141. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Okamoto, T. Utsunomiya, S. Wakiyama, M. Hashimoto, K. Fukuzawa, T. Ezaki, T. Hanai, H. Inoue, and M. Mori Specific Gene-Expression Profiles of Noncancerous Liver Tissue Predict the Risk for Multicentric Occurrence of Hepatocellular Carcinoma in Hepatitis C Virus-Positive Patients Ann. Surg. Oncol., July 1, 2006; 13(7): 947 - 954. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. Mendrzyk, A. Korshunov, A. Benner, G. Toedt, S. Pfister, B. Radlwimmer, and P. Lichter Identification of gains on 1q and epidermal growth factor receptor overexpression as independent prognostic markers in intracranial ependymoma. Clin. Cancer Res., April 1, 2006; 12(7): 2070 - 2079. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. Lopez-Rios, S. Chuai, R. Flores, S. Shimizu, T. Ohno, K. Wakahara, P. B. Illei, S. Hussain, L. Krug, M. F. Zakowski, et al. Global Gene Expression Profiling of Pleural Mesotheliomas: Overexpression of Aurora Kinases and P16/CDKN2A Deletion as Prognostic Factors and Critical Evaluation of Microarray-Based Prognostic Prediction. Cancer Res., March 15, 2006; 66(6): 2970 - 2979. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. Mendrzyk, B. Radlwimmer, S. Joos, F. Kokocinski, A. Benner, D. E. Stange, K. Neben, H. Fiegler, N. P. Carter, G. Reifenberger, et al. Genomic and Protein Expression Profiling Identifies CDK6 As Novel Independent Prognostic Marker in Medulloblastoma J. Clin. Oncol., December 1, 2005; 23(34): 8853 - 8862. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. M. Egan, P. A. Newcomb, C. B. Ambrosone, A. Trentham-Dietz, L. Titus-Ernstoff, J. M. Hampton, M. T. Kimura, and H. Nagase STK15 polymorphism and breast cancer risk in a population-based study Carcinogenesis, November 1, 2004; 25(11): 2149 - 2153. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||