Advertisement
Journal of Clinical Oncology  
Search for:
Limit by:
  Browse by Subject or Issue
Home Search or Browse JCO My JCO Subscriptions Customer Service Site Map

Journal of Clinical Oncology, Vol 26, No 18 (June 20), 2008: pp. 3015-3024
© 2008 American Society of Clinical Oncology.
DOI: 10.1200/JCO.2007.15.7164

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Purchase Article
Right arrow View Shopping Cart
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a colleague
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Save to my personal folders
Right arrow Download to citation manager
Right arrowRights & Permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Murat, A.
Right arrow Articles by Hegi, M. E.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Murat, A.
Right arrow Articles by Hegi, M. E.
Social Bookmarking
 Add to CiteULike   Add to Complore   Add to Connotea   Add to Del.icio.us   Add to Digg   Add to Facebook   Add to Reddit   Add to Technorati   Add to Twitter  
What's this?

Stem Cell–Related "Self-Renewal" Signature and High Epidermal Growth Factor Receptor Expression Associated With Resistance to Concomitant Chemoradiotherapy in Glioblastoma

Anastasia Murat, Eugenia Migliavacca, Thierry Gorlia, Wanyu L. Lambiv, Tal Shay, Marie-France Hamou, Nicolas de Tribolet, Luca Regli, Wolfgang Wick, Mathilde C.M. Kouwenhoven, Johannes A. Hainfellner, Frank L. Heppner, Pierre-Yves Dietrich, Yitzhak Zimmer, J. Gregory Cairncross, Robert-Charles Janzer, Eytan Domany, Mauro Delorenzi, Roger Stupp, Monika E. Hegi

From the Laboratory of Tumor Biology and Genetics, Centre Universitaire Romand de Neurochirurgie; Division of Neuropathology; Multidisciplinary Oncology Center, Centre Hospitalier Universitaire Vaudois and University of Lausanne, Lausanne; the National Center of Competence in Research Molecular Oncology at the Swiss Institute of Experimental Cancer Research; the Swiss Institute of Bioinformatics, Epalinges; the Centre Universitaire Romand de Neurochirurgie, the Service of Oncology, University Hospital Geneva, Geneva; the Department of Clinical Research and Radiation Oncology, Inselspital and University of Berne, Berne; the Institute of Neuropathology, University Hospital Zurich, Zurich, Switzerland; Data Center, European Organisation for Research and Treatment of Cancer (EORTC), Brussels, Belgium; Department of Physics of Complex Systems, Weizmann Institute of Science, Rehovot, Israel; the Department of General Neurology and Hertie Institute for Clinical Brain Research, University of Tübingen, Tübingen; Department of Neurooncology, University Clinic Heidelberg, Heidelberg; Institute of Neuropathology, Charité –Universitätsmedizin Berlin, Berlin, Germany; Department of Neurology, Erasmus, Rotterdam, the Netherlands; Institute of Neurology, Medical University of Vienna, Vienna, Austria; and the Department of Clinical Neurosciences, University of Calgary, Calgary, Canada

Corresponding author: Monika E. Hegi, PhD, Laboratory of Tumor Biology and Genetics, Centre Universitaire Romand de Neurochirurgie, Centre Hospitalier Universitaire Vaudois and University of Lausanne, 46, rue du Bugnon, Lausanne 1011/Switzerland; e-mail: Monika.Hegi{at}chuv.ch


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 PATIENTS AND METHODS
 RESULTS
 DISCUSSION
 AUTHORS' DISCLOSURES OF...
 AUTHOR CONTRIBUTIONS
 Appendix
 Glossary Terms
 REFERENCES
 
Purpose Glioblastomas are notorious for resistance to therapy, which has been attributed to DNA-repair proficiency, a multitude of deregulated molecular pathways, and, more recently, to the particular biologic behavior of tumor stem-like cells. Here, we aimed to identify molecular profiles specific for treatment resistance to the current standard of care of concomitant chemoradiotherapy with the alkylating agent temozolomide.

Patients and Methods Gene expression profiles of 80 glioblastomas were interrogated for associations with resistance to therapy. Patients were treated within clinical trials testing the addition of concomitant and adjuvant temozolomide to radiotherapy.

Results An expression signature dominated by HOX genes, which comprises Prominin-1 (CD133), emerged as a predictor for poor survival in patients treated with concomitant chemoradiotherapy (n = 42; hazard ratio = 2.69; 95% CI, 1.38 to 5.26; P = .004). This association could be validated in an independent data set. Provocatively, the HOX cluster was reminiscent of a "self-renewal" signature (P = .008; Gene Set Enrichment Analysis) recently characterized in a mouse leukemia model. The HOX signature and EGFR expression were independent prognostic factors in multivariate analysis, adjusted for the O-6-methylguanine-DNA methyltransferase (MGMT) methylation status, a known predictive factor for benefit from temozolomide, and age. Better outcome was associated with gene clusters characterizing features of tumor-host interaction including tumor vascularization and cell adhesion, and innate immune response.

Conclusion This study provides first clinical evidence for the implication of a "glioma stem cell" or "self-renewal" phenotype in treatment resistance of glioblastoma. Biologic mechanisms identified here to be relevant for resistance will guide future targeted therapies and respective marker development for individualized treatment and patient selection.


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 PATIENTS AND METHODS
 RESULTS
 DISCUSSION
 AUTHORS' DISCLOSURES OF...
 AUTHOR CONTRIBUTIONS
 Appendix
 Glossary Terms
 REFERENCES
 
In glioblastoma, introduction of combined chemoradiotherapy of concomitant and adjuvant temozolomide (TMZ) and radiotherapy (TMZ/RT->TMZ) has allowed significant prolongation of survival,1,2 in particular in patients with an epigenetically silenced O-6-methylguanine-DNA methyltransferase (MGMT) DNA repair gene.3,4 However, outcome remains unsatisfactory, and ongoing clinical trials explore modulation of MGMT or the addition of targeted agents.1,5 Recognizing molecular tumor signatures of underlying biologic processes associated with resistance in patients treated with this new standard therapy will allow the identification of potential targets for improvement of therapy.

Recent concepts for cancer development suggest that a minority population of cancer stem-like cells may determine the biologic behavior of tumors, including response to therapy. Failure to cure cancer has been attributed to the fact that therapies are aimed at the tumor bulk without significantly harming tumor stem-like cells,6 supported by experimental evidence in a respective mouse model showing that this glioblastoma subpopulation of cells is more resistant to radiotherapy.7 Facilitated by markers differentiating stem-cells and progenitors of the different lineages, the origin of leukemic stem cells has been traced back to hematopoietic stem cells as well as progenitor populations that have acquired self-renewal properties.8,9 In contrast, the origin and concept of glioma stem-like cells remains to be fully elucidated. CD133 positivity has been postulated to be a glioma stem-cell marker, considering that this subpopulation of glioma-derived cells seems to have a higher potential to generate and maintain tumors in vivo.7,10

With the goal of identifying molecular mechanisms of treatment resistance, we report here on the analysis of glioblastoma-derived gene expression profiles from patients treated within two prospective clinical trials.1,2 We have identified several biologic processes associated with resistance or responsiveness to combined chemoradiotherapy that provide important information guiding novel treatment strategies and aiming at individualized therapy. Intriguingly, an expression signature associated with resistance shows high similarity with a stem cell–related self-renewal signature.9 Our data provide the first clinical evidence to our knowledge for the involvement of a tumor stem-cell phenotype in the escape of glioblastoma from chemoradiotherapy.


    PATIENTS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 PATIENTS AND METHODS
 RESULTS
 DISCUSSION
 AUTHORS' DISCLOSURES OF...
 AUTHOR CONTRIBUTIONS
 Appendix
 Glossary Terms
 REFERENCES
 
Tumor Samples and Patient Characteristics
Gene expression profiles were established from 80 frozen glioblastoma samples obtained from 76 patients, comprising 70 tumors from initial surgery and 10 samples resected at recurrence, and from four non-neoplastic brain tissue samples. All patients were treated within a phase II or a randomized phase III trial1,2 and provided written informed consent for molecular studies of their tumor. The protocol was approved by the ethics committee at each center. Sixty-eight patients with complete molecular and clinical information were included in survival analysis, with a median age of 51 years (range, 26 to 70 years) at enrollment. Thereof, 42 received TMZ/RT->TMZ, and 26 were randomly assigned to RT only. Second-line therapy frequently involved alkylating agents including TMZ. Patient characteristics are summarized in Table A1(online only). The validation set comprised 76 independent patients of the European Organisation for Research and Treatment of Cancer (EORTC)/National Cancer Institute of Canada (NCIC) study,1 39 randomly assigned to TMZ/RT->TMZ, and 37 to RT (median age, 54 years; range, 25-69 years), whose glioblastomas were available on a tissue microarray (TMA). There was no difference in survival compared with the general trial population, neither in the test population nor the validation set (P > .2).

Gene Expression Profiling
Total RNA was extracted from frozen tumor sections (10 to 15 mg; Qiagen RNeasy-Lipid Tissue Kit; Hilden, Germany). The first section served as reference for diagnosis and tumor content (> 70%). Probes were prepared with the Enzo BioArray-High Yield Kit (Enzo Life Sciences, Farmingdale, NY) for double amplification and were hybridized to Affymetrix HG-133Plus2.0 GeneChips (Affymetrix, Santa Clara, CA). The microarray data is deposited in the Gene Expression Omnibus (GEO) database at http://www.ncbi.nlm.nih.gov/geo/ (accession-number GSE7696). Quantitative reverse transcription polymerase-chain reaction (qRT-PCR) was performed on the ABI Prism7900 with SYBR Green (Applied Biosystems, Foster City, CA). Primers are listed in Table A2 (online only). Results were normalized to the expression of the EIF2C3, DNAJA4, and B2M genes that exhibit little variation in this data set.

Data Analysis and Statistical Methods
Analyses were carried out in R, a free software environment available at http://www.r-project.org/, SAS (V9.1.3; SAS Institute, Cary, NC), or Coupled Two Way Clustering (CTWC),11-13 available at http://ctwc.weizmann.ac.il. The expression intensities for all probe sets from Affymetrix CEL-files were estimated using robust multiarray average with probe-level quantile normalization followed by median polish summarization as implemented in the BioConductor software (http://www.bioconductor.org/).

The expression matrix of 84 samples and 3,860 most variable probe sets (standard deviation, > 0.75) was input into CTWC using default parameters and two levels of clustering. CTWC analysis can be viewed at: http://bcf.isb-sib.ch/projects/cancer/glio/. Probe sets comprising stable gene clusters emerging from CTWC served as input for supervised analyses.

Details on the Cox and partial least square models, and published data sets used9,14-16 are described in the Appendix and Table A3 (online only). The Benjamini-Hochberg procedure was applied for multiple testing correction (false-discovery rates).17

TMA and Immunohistochemistry
The TMA was constructed with glioblastoma from the EORTC/NCIC-trial.1 HOXA10 (Santa Cruz Biotechnology, Santa Cruz, CA; sc-17159; dilution, 1:200) expression was determined by immunohistochemistry (citrate buffer: pH 6.0, 15 minutes in pressure cooker) and scored without knowledge of clinical information on a scale of 0 to 6 (0, no expression; 1, weak nuclear expression in < 20% cells; 2, strong in < 20% cells; 3, weak in 20% to 50% cells; 4, strong in 20% to 50% cells; 5, weak in > 50% cells; 6, strong in > 50% cells). Dichotomization for survival analysis was no to low score (1, 2) v intermediate to high expression (scores 2-6). Frozen sections of neurospheres were fixed with acetone and stained for HOXA10 and CD133 (Santa Cruz; sc-17159, dilution 1:100; Miltenyi AC133-2, 293C3, dilution 1:50).

Glioblastoma-Derived Neurospheres
Fresh glioblastoma tissue was dissociated in presence of papain and DNase I basically as described previously.18 Cells were cultured under stem-cell conditions to form spheres using Dulbecco's modified Eagle medium/F12 medium containing B27 supplement and 20 ng/mL of both epidermal growth factor and fibroblast growth factor 2. Neurospheres from six glioblastomas propagated between 4 weeks and 14 months were used for immunostaining.


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 PATIENTS AND METHODS
 RESULTS
 DISCUSSION
 AUTHORS' DISCLOSURES OF...
 AUTHOR CONTRIBUTIONS
 Appendix
 Glossary Terms
 REFERENCES
 
Gene Expression Signatures Associated With Tumor Resistance
There was no obvious association of patient characteristics or survival with genetic subtypes evident from the sample dendrogram S1(G1), clustering all samples (S1) and all genes (G1) passing a variation filter (Fig 1). The methylation status of the MGMT gene promoter appears to be randomly distributed. All stable gene clusters emerging from this analysis are listed in Table 1, named by the predominant biologic function suggested by the genes they comprise, whereas their inter-relationship is visualized in Fig 2A. Similar gene clusters are obtained when using only the subset of glioblastoma clustered in S4 [G1(S4); Appendix Table A4, online only]. Cluster S4 comprises most glioblastoma (69 of 80), but not the nontumoral tissues that form their own stable cluster, S3 (Fig 1). Similar gene clusters are present in other glioblastoma data sets as visualized in Fig 2B for data sets published by the groups of Nelson and Aldape, respectively.14,15All 18 nonoverlapping, stable gene clusters [G1(S1)] were interrogated for association with survival using Cox proportional hazards, adjusted for age (> 50 years) and MGMT methylation status.3,4 Seven gene clusters were most influential for explaining survival in patients assigned to TMZ/RT->TMZ (Table 1; gene lists, Appendix Tables A5-A10, online only). A respective partial least square model yielded comparable results (Appendix Fig A1, online only). However, the MGMT methylation status was yet the most influential predictor of survival (Table 2; Appendix Fig A1, online only). For this study, we focused on the two most significant clusters characterized by HOX (homeobox) gene (G28/G98) and EGFR (epidermal growth factor receptor) gene expression (G25), respectively, and here we briefly comment on the biologic and clinical implications denoted by the other relevant clusters.


Figure 1
View larger version (36K):
[in this window]
[in a new window]
[PowerPoint Slide for Teaching]
 
Fig 1. Sample dendrogram S1(G1) by Coupled Two-Way Clustering (CTWC). Sample dendrogram S1(G1) emerging from clustering all 84 samples S1 with all genes G1. Stable sample clusters S2 to S6 emerge. Clinical information: age at diagnosis (green, < 50 years; brown, > 50 years), overall survival in months (green, < 9; red, 9-18; pink > 18); MGMT methylation status (gray, methylated; black, unmethylated; white, unknown); sex (red, female); and sample Id# (sample identification number, nontumoral brain tissue, yellow; recurrent glioblastoma, black; color code tumor pairs from same patient, black only recurrent glioblastoma).

 

View this table:
[in this window]
[in a new window]

 
Table 1. Cox Analysis Using Stable Gene Clusters Identified in G1(S1) in the Cohort of Treated Patients (n = 42)

 

Figure 2
View larger version (102K):
[in this window]
[in a new window]
[PowerPoint Slide for Teaching]
 
Fig 2. (A): The distance matrix of all filtered genes G1 (Euclidian distance) heatmap; blue indicates short distance (more similar), red indicates large distance. The distance ranges from 0 (corresponding to Pearson correlation = 1; ie, full similarity of the two expression profiles) to 2 (corresponding to anticorrelated profiles, Pearson = –1). Zero correlation indicated by orange. Annotation is based on gene dendrogram G1(S1), not shown. The distance matrix visualizes relationships between clusters as indicated by lines and arrows (eg, G7 is in the center of a supercluster, reflecting close relationship with G9, G12, G14 and G2). In addition G7 is related to G29, but not G23 or G21. (B) Respective distance matrices for two published data sets with 48 and 54 glioblastomas14,15 comprising the common probe sets (3,649) and ordered according to the matrix in Figure 2A. All three data sets appear similar, particularly visible for the largest clusters annotated in Figure 2A.

 

View this table:
[in this window]
[in a new window]

 
Table 2. Overall Cox Proportional Hazards Models for the Cohort of 42 Patients Treated With Concomitant and Adjuvant Temozolomide and Radiotherapy

 
Self-Renewal Signature Associated With Resistance to Chemoradiotherapy in Glioblastoma
Increased expression of cluster G28, dominated by HOX genes and comprising the cell-cycle checkpoint gene GADD45G,19 was found to be associated with worse outcome (P = .004; hazard ratio [HR] = 2.69; 95% CI, 1.38 to 5.26; Table 1). HOX genes are essential in axis determination during embryonic development and are known to be involved in cancer including glioblastomas.20,21 The interaction term between this HOX cluster and chemoradiotherapy was significant (P = .001) in the Cox model when evaluating all 68 patients from the two treatment arms, implying that high expression of HOX genes may be predictive for resistance to TMZ/RT->TMZ therapy.

The HOX gene clusters G28 and G98 emerging from clustering all genes (G1) either with all samples (S1) or only with glioblastoma clustered in S4 are almost identical (Pearson correlation, 0.98; 19 of 21 probe-sets; Appendix Table A11, online only). Intriguingly, G98 in addition comprises Prominin 1, encoding the putative glioma stem–cell marker CD133 (Fig 3A), suggesting that in a subpopulation of glioblastoma, concerted upregulation of HOX genes might be associated with a tumor stem-like cell phenotype. In accordance, we find high HOXA10 protein expression in glioblastoma-derived neurospheres, cultured under stem-cell conditions, as displayed in Fig 3B together with CD133 expression. To include information on Prominin 1, we show results on G98 for all following analyses. Fig 3C visualizes the association of short survival with enhanced expression of G98.


Figure 3
View larger version (36K):
[in this window]
[in a new window]
[PowerPoint Slide for Teaching]
 
Fig 3. Relation of glioblastoma derived HOX gene signature with survival. (A) Stable cluster G98 was found by clustering all genes G1 with samples clustered in S4 (69 glioblastomas). Red indicates relative overexpression and blue relative underexpression. (B) Glioblastoma-derived neurospheres (40x magnification) cultured under stem-cell conditions exhibit strong nuclear staining for HOXA10 and express CD133 as determined by immunohistochemistry. Representative neurospheres from two glioblastoma are shown. (C) Kaplan-Meier survival estimates of the 42 patients treated with concomitant and adjuvant temozolomide (TMZ) and radiotherapy (TMZ/RT->TMZ) separated into high (High HOX) versus low (Low HOX) expressors of G98 genes (dichotomized according to Coupled Two-Way Clustering sample dendrogram). The P value of the log-rank test is shown for the two groups, stratified by the MGMT methylation status. M-MGMT, methylated MGMT; U-MGMT, unmethylated MGMT.

 
Validation in Independent Data Sets
The association of the HOX signature with resistance to treatment was subsequently validated in a sample set of the trial not available for initial discovery, arrayed on TMA. HOXA10 was evaluated by immunohistochemistry (Fig 4) as a representative of the correlated set of HOX genes. Strong nuclear expression was often observed in patches of tumor cells situated in the vicinity of blood vessels. High HOXA10 expression was associated with worse outcome in patients randomly assigned to TMZ/RT->TMZ therapy (n = 39; HR = 2.57; 95% CI, 1.21 to 5.47; P = .014; Fig 5). The patients randomly assigned to RT only did not show such a relationship, suggesting that high HOXA10 expression may be predictive for resistance to a synergistic effect of concomitant chemoradiotherapy, in concordance with the significant interaction term between treatment and expression of G28 or G98 (Appendix Table A12, online only). Similar HOX clusters can be identified in the Nelson and Aldape glioblastoma data sets.14,15 Correlating G98 gene expression with outcome in these data sets totaling 102 glioblastoma revealed a trend for worse outcome (P = .09; HR = 1.29; 95% CI, 0.97 to 1.72; Appendix Fig A2, online only). Of note, in contrast to our data, these patients were treated before the TMZ/RT->TMZ regimen was established. In accordance with better survival, anaplastic glioma (WHO grade 3) profiled in these publications revealed significantly lower expression of G98 genes compared with glioblastoma (WHO grade 4, P < .001 for the Aldape data set14; P = .002 for the Nelson data set15; Wilcoxon rank sum test with continuity correction; Appendix Fig A3, online only). However, within grade 3 glioma of the two data sets increased expression of G98 genes was associated with worse outcome (n = 44; P = .007; HR = 3.35; 95% CI, 1.39 to 8.08; Appendix Fig A2).


Figure 4
View larger version (82K):
[in this window]
[in a new window]
[PowerPoint Slide for Teaching]
 
Fig 4. HOXA10 expression determined by immunohistochemistry on tumor microarray. Diameter of tissues in upper panel is 0.6 mm. (A) High nuclear expression; (B) focal high nuclear expression; and (C) no expression.

 

Figure 5
View larger version (11K):
[in this window]
[in a new window]
[PowerPoint Slide for Teaching]
 
Fig 5. HOXA10 expression and overall survival in independent patient cohort. Kaplan-Meier survival estimates of (A) 39 patients randomly assigned to concomitant and adjuvant temozolomide (TMZ) and radiotherapy (TMZ/RT->TMZ) and (B) 37 patients randomly assigned to RT, separated into high versus low expressors of HOXA10 as determined by immunohistochemistry on tumor microarray.

 
HOX Signature Reminiscent of Self-Renewal
The core of the HOX cluster was found to be part of the top 20 genes of a self-renewal signature identified by Krivtsov et al9 in murine MLL-AF9-induced leukemic stem cells derived from committed progenitors. Our G98-derived signature was significantly enriched in genes discriminating self-renewal versus non–self-renewal in this expression data set according to Gene Set Enrichment Analysis (GSEA; P = .008; Appendix Fig A4 and Table A13, online only) and the Wilcoxon two-sample test (G98; P < .001).22 The relevance of our signature was further demonstrated in a human data set16 in which G98, similar to the original murine self-renewal signature,9 was able to significantly differentiate MLL-rearranged–acute myeloid leukemia from acute myeloid leukemia (GSEA, P < .001; Wilcoxon two-sample test P = .01; Appendix Table A14, online only).

Interestingly, a significant, although low, correlation between the mean DNA copy number of the two BAC (bacterial artificial chromosome) clones (GS1-213H12 and CTB-23D20) bordering the HOXA gene locus on chromosome 7 and the mean expression of the HOXA genes was observed in a set of 60 glioblastoma (Pearson correlation coefficient r = 0.27; P = .03; manuscript in preparation). These flanking BACs were more amplified than their neighbors (Appendix Fig A4, online only). Hence, the herein-proposed HOX-dominated self-renewal signature for glioblastoma may in part be acquired by increased gene dosage.

High EGFR Expression Is Associated With Tumor Resistance
Two gene clusters representing amplification-mediated overexpression of proto-oncogenes were associated with tumor resistance: G13 (P = .02; HR = 1.1; 95% CI, 1.0 to 1.2) characterized by coordinated upregulation of contiguous genes on chromosome 12q13-15, comprising the proto-oncogenes CDK4 and MDM2, and G25, dominated by EGFR probe sets (P = .002; HR = 2.8; 95% CI, 1.4 to 5.4; Table 1; Fig 6). The array-derived EGFR measurement (probe set 201983_s_at) was confirmed by qRT-PCR (Pearson correlation, 0.89). Expression of the constitutively activated EGFRvIII mutant (18 of 70; 26%), as measured by qRT-PCR, did not further influence outcome prediction (P = .94).


Figure 6
View larger version (17K):
[in this window]
[in a new window]
[PowerPoint Slide for Teaching]
 
Fig 6. Association of EGFR expression with survival. Kaplan-Meier curves of the 42 patients in the concomitant and adjuvant temozolomide (TMZ) and radiotherapy (TMZ/RT->TMZ) cohort divided between low and high expression of EGFR (probe set 201983_s_at, dichotomized according to a Gaussian mixture model). The P value from log-rank test is shown for two groups defined by EGFR expression and stratified according to the MGMT methylation status. M-MGMT, methylated MGMT; U-MGMT, unmethylated MGMT.

 
Cluster G18, associated with tumor resistance (P = .03; HR = 1.94; 95% CI, 1.07 to 3.51), displayed some correlation with EGFR expression (r = 0.57; Fig 2A) in particular with Aquaporin 4 (AQP4). AQP4 has been associated with brain tumor related edema.23 Aquaporins require activation of mitogen-activated protein kinase signaling that may be mediated by EGFR activation.24 Another family member, AQP1, has been linked with tumor angiogenesis and cell migration25 and was associated with worse outcome in this study (P = .003; HR = 2.44, 95% CI, 1.36 to 4.04) and in the two external data sets (combined P = .009; HR = 1.51; 95% CI, 1.11 to 2.06).

Blood Vessels Markers Associated With Better Outcome
G7 is characterized by genes associated with endothelial cells, basement membranes, signaling pathways of vascular development and angiogenesis, and tumor-derived endothelial markers26,27 (Table A5). Hence, G7 may differentiate tumors according to their angiogenic pattern that may be indicative for drug perfusion, and therefore show association with benefit from treatment (Table 1). G7 is in the center of a "super cluster" (Fig 2A) constituted of several stable gene clusters with biologically related features, such as hypoxia-regulated genes G9, and the myeloid progenitor/adhesion cluster G14, which is also correlated with better outcome (Table 1). Beside genes related to cell adhesion, mesenchymal stem cells (PRR16), and homeobox genes (MEOX1, MKX), G14 includes aldehyde dehydrogenase28 and bone morphogenetic protein 5 gene (BMP5)29 both associated with differentiation of stem cells. The cluster comprises positive (MEOX1) and negative regulators (SOSTDC1) of BMPs, which recently have been shown to inhibit tumorigenic potential of human brain tumor stem cells by promoting their differentiation.29 This cluster may reflect the perivascular microenvironment proposed recently to serve as niche for brain tumor stem cells.30

Innate Immune Response Associated With Better Survival
Alongside markers of innate immunity and macrophages (CD11b), G24 comprises numerous cell-surface receptors known as markers for M2-polarized macrophages,31 such as CD163 (Appendix Table A9, online only). M2 polarization mediates tolerance and downregulates inflammation, alleviating immune surveillance.31,32 A wide range of CD163-positive cells can be observed in glioblastoma (Appendix Fig A5, online only). The cluster also contains probes encoding major histocompatibility complex class II surface molecules, but lacks expression of costimulatory molecules critical for T-cell activation.33 This immune signature may be relevant for strategies of tumor vaccination.

HOX Signature and EGFR Expression Are Independent Prognostic Factors
Multivariate analysis suggests that the HOX signature and EGFR expression, respectively, are independent prognostic factors for poor outcome in TMZ/RT->TMZ–treated glioblastoma patients, explaining 67% of the survival outcome variations together with the MGMT status and age (Table 2).


    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 PATIENTS AND METHODS
 RESULTS
 DISCUSSION
 AUTHORS' DISCLOSURES OF...
 AUTHOR CONTRIBUTIONS
 Appendix
 Glossary Terms
 REFERENCES
 
Gene expression signatures relevant for treatment resistance to TMZ/RT->TMZ have been identified in a prospectively treated population of glioblastoma patients. Although epigenetic inactivation of the MGMT gene promoter remained the most prominent predictive factor, expression signatures allowed identification of patient subgroups that may benefit from specific additional therapies targeting particular mechanisms of resistance.

As an independent predictive factor of resistance, we have identified a HOX-dominated gene cluster, evocative of a self-renewal signature recently described for leukemia.9 Strong HOXA10 expression of glioblastoma-derived neurospheres is in line with a role of HOX genes in the glioma stem-like cell compartment. These findings provide the first clinical evidence to our knowledge for the relevance of a stem-like cell phenotype in treatment resistance of glioblastoma.

In leukemia, expression of translocation-related fusion proteins led to MLL-mediated chromatin remodeling associated with re-expression of HOX genes.9,34,35 In glioblastoma, however, such fusion proteins have not been described, and no indications from the present data set link MLL expression with the self-renewal signature. We provide evidence that the HOX-dominated gene signature emerges with malignant progression to glioblastoma, and may be acquired in some glioblastomas by low-level amplification, the latter supporting the notion that gliomas may also arise from progenitors, in agreement with mouse models.36 The identification of GADD45G as part of the HOX signature may provide further evidence for an enhanced DNA repair potential that has recently been associated with radiation resistance of glioma stem cells.7

These molecular and clinical data underscore the importance of the self-renewal phenotype that could be explored as a potential treatment target.37 First efforts blunting glioma stem cell–related self-renewal properties of tumors suggest that strategies forcing differentiation (eg, mediated by cytokines such as BMP4) may be promising.29

Targeting resistance to TMZ/RT->TMZ associated with overexpression of the EGFR gene is of particular clinical interest because this alteration affects a large proportion of patients, and inhibitors of several key targets in downstream signaling pathways are already available or are in advanced clinical development.38 Similarly, inhibitors for CDK4 and MDM2 are under investigation, and both proto-oncogenes showed amplification-mediated overexpression, identified here as associated with worse outcome.

Surprisingly, good prognosis was associated with increased expression of a signature for tumor endothelium markers. This signature may predict improved cytotoxic activity by means of better perfusion of the tumor with the chemotherapy agent TMZ. This is in accordance with the concept suggesting that antiangiogenic agents may temporarily lead to "normalization" of aberrant tumor vasculature, resulting in more efficient delivery of drugs and oxygen to the tumor.39,40 In a recent trial, addition of the antiangiogenic integrin-inhibitor cilengitide appears to confer increased antitumor activity in conjunction with TMZ/RT->TMZ in patients with a methylated MGMT gene promoter.5

Another interesting insight of our study suggests infiltration of M2-polarized macrophages into the tumors. The altered capacity of these glioma-infiltrating macrophages to induce effective antitumor T-cell response may obstruct therapeutic strategies aimed at boosting adaptive immunity against the tumor. M2 polarization is driven by tumor- and T-cell–derived cytokines,31 consistent with the well-known expression of the immunosuppressive cytokines transforming growth factor-β and interleukin-10 in malignant glioma.41 Thus, for effective immunotherapy/vaccination, full resection of the tumor may be required to remove the microenvironment conferring immunosuppression and tolerance.

The molecular signatures identified in this study associated with outcome underline the need for development of multimodality treatments targeting not only the tumor cells, but including strategies aimed at the glioma stem-like cell compartment, and interfering with tumor-host interaction that provides the specialized microenvironment relevant for the maintenance of tumor stem-like cells (the stem-cell niche), angiogenesis, and immune response. The hypotheses generated in this study need to be tested in prospective trials. Our findings may guide a rational choice of agents, targets, trial design, and appropriate patient selection, incorporating biomarkers defining mechanisms of response and resistance.


    AUTHORS’ DISCLOSURES OF POTENTIAL CONFLICTS OF INTEREST
 TOP
 ABSTRACT
 INTRODUCTION
 PATIENTS AND METHODS
 RESULTS
 DISCUSSION
 AUTHORS' DISCLOSURES OF...
 AUTHOR CONTRIBUTIONS
 Appendix
 Glossary Terms
 REFERENCES
 
Although all authors completed the disclosure declaration, the following author(s) indicated a financial or other interest that is relevant to the subject matter under consideration in this article. Certain relationships marked with a "Urdquo; are those for which no compensation was received; those relationships marked with a "C" were compensated. For a detailed description of the disclosure categories, or for more information about ASCO's conflict of interest policy, please refer to the Author Disclosure Declaration and the Disclosures of Potential Conflicts of Interest section in Information for Contributors.

Employment or Leadership Position: None Consultant or Advisory Role: Roger Stupp, Schering-Plough (C), Merck Serono (C); Monika E. Hegi, Oncomethylome Sciences (C) Stock Ownership: None Honoraria: Wolfgang Wick, Schering-Plough, Essex-Pharma; J. Gregory Cairncross, Schering-Plough; Roger Stupp, Schering-Plough, Merck Serono; Monika E. Hegi, Schering-Plough Research Funding: Monika E. Hegi, Schering-Plough, Oncomethylome Sciences Expert Testimony: None Other Remuneration: None


    AUTHOR CONTRIBUTIONS
 TOP
 ABSTRACT
 INTRODUCTION
 PATIENTS AND METHODS
 RESULTS
 DISCUSSION
 AUTHORS' DISCLOSURES OF...
 AUTHOR CONTRIBUTIONS
 Appendix
 Glossary Terms
 REFERENCES
 
Conception and design: Anastasia Murat, Eugenia Migliavacca, Monika E. Hegi

Financial support: Nicolas de Tribolet, Monika E. Hegi

Administrative support: Roger Stupp, Monika E. Hegi

Provision of study materials or patients: Nicolas de Tribolet, Luca Regli, Wolfgang Wick, Mathilde C.M. Kouwenhoven, Johannes A. Hainfellner, Frank L. Heppner, Pierre-Yves Dietrich, Yitzhak Zimmer, Gregory Cairncross, Roger Stupp

Collection and assembly of data: Anastasia Murat, Eugenia Migliavacca, Marie-France Hamou, Robert-Charles Janzer, Roger Stupp, Monika E. Hegi

Data analysis and interpretation: Anastasia Murat, Eugenia Migliavacca, Thierry Gorlia, Wanyu L. Lambiv, Tal Shay, Marie-France Hamou, Eytan Domany, Mauro Delorenzi, Monika E. Hegi

Manuscript writing: Anastasia Murat, Eugenia Migliavacca, Thierry Gorlia, Eytan Domany, Mauro Delorenzi, Roger Stupp, Monika E. Hegi

Final approval of manuscript: Anastasia Murat, Eugenia Migliavacca, Thierry Gorlia, Wanyu L. Lambiv, Tal Shay, Marie-France Hamou, Nicolas de Tribolet, Luca Regli, Wolfgang Wick, Mathilde C.M. Kouwenhoven, Johannes A. Hainfellner, Frank L. Heppner, Yitzhak Zimmer, Gregory Cairncross, Robert-Charles Janzer, Eytan Domany, Mauro Delorenzi, Roger Stupp, Monika E. Hegi


    Appendix
 TOP
 ABSTRACT
 INTRODUCTION
 PATIENTS AND METHODS
 RESULTS
 DISCUSSION
 AUTHORS' DISCLOSURES OF...
 AUTHOR CONTRIBUTIONS
 Appendix
 Glossary Terms
 REFERENCES
 
Go


Figure 7
View larger version (15K):
[in this window]
[in a new window]
[PowerPoint Slide for Teaching]
 
Fig A1. Partial least square (PLS) model for prediction of overall survival. Cluster X weights were obtained by averaging the X weights of their constituent genes. The clusters are numbered according to the legend and the description in Table 1. Age and MGMT methylation status were used as covariates. Clusters of genes with X weights that are nearest to the PLS factors, represented by the axes (PLS factor 1 and 2), and the farthest from the center of the plane contribute most to the PLS regression. The first two factors shown explain 66% of the survival outcome variations. Clusters of genes grouped in the upper and right side of the plane have a positive association with shorter survival (ie, higher hazards), while those in the lower and left side are positively associated with longer survival (ie, lower hazards).

 
Go


Figure 8
View larger version (9K):
[in this window]
[in a new window]
[PowerPoint Slide for Teaching]
 
Fig A2. Expression of G98 and association with survival in external data sets. High expression of the self-renewal cluster G98 is significantly associated with worse outcome in two independent data sets comprising a total of 146 malignant glioma (P = .007, hazard ratio = 1.46, 95% CI, 1.11 to 1.92; model adjusted for tumor grade WHO grade 3 and 4; age, and stratified for the data set, Nelson [Freije WA, Castro-Vargas FE, Fang Z, et al: Cancer Res 64:6503-6510, 2004]; Aldape [Phillips HS, Kharbanda S, Chen R, et al: Cancer Cell 9:157-173, 2006]). The forest plots visualize the relationship between expression of G98 and outcome in the two data sets separated for tumor grade. The hazard ratio for glioblastoma (WHO grade 4) is 1.29 (95% CI, 0.97 to 1.72; P = .09; left panel). A hazard ratio of 3.35 (95% CI, 1.39 to 8.08; P = .007) is observed for the subset of anaplastic glioma (WHO grade 3; right panel).

 
Go


Figure 9
View larger version (13K):
[in this window]
[in a new window]
[PowerPoint Slide for Teaching]
 
Fig A3. Box plot for expression of G98 in grade 3 and grade 4 glioma of external data sets. G98 expression significantly differentiates anaplastic glioma (WHO grade 3) from glioblastoma (WHO grade 4). Glioblastomas exhibit significantly higher expression of G98 in both external data sets (Wilcoxon rank sum test with continuity correction, Nelson [Freije WA, Castro-Vargas FE, Fang Z, et al: Cancer Res 64:6503-6510, 2004] P < .001; Aldape [Phillips HS, Kharbanda S, Chen R, et al: Cancer Cell 9:157-173, 2006] P = .002).

 
Go


Figure 10
View larger version (32K):
[in this window]
[in a new window]
[PowerPoint Slide for Teaching]
 
Fig A4. Array-based comparative genomic hybridization data for the chromosomal region around the HOXA gene cluster on chromosome 7. Every row is a marker, ordered by chromosomal location. Every column is a sample (n = 60), ordered by the copy number of BAC GS1-213H12 (bacterial artificial chromosome) using SPIN software (Tsafrir D, Tsafrir I, Ein-Dor L, et al: Bioinformatics 21:2301-2308, 2005). Color code represents the log2 (tumor/normal) ratio of the copy number, red corresponding to amplification. The HOXA gene cluster resides between markers GS1-213H12 and CTB-23D20 that appear to be amplified more frequently and more strongly than the neighboring markers.

 
Go


Figure 11
View larger version (88K):
[in this window]
[in a new window]
[PowerPoint Slide for Teaching]
 
Fig A5. CD163 expression in glioblastoma. Immunohistochemistry for CD163 (Novocastra, Newcastle upon Tyne, United Kingdom; NCL-163; dilution 1:800) on representative paraffin-embedded glioblastoma. Diameter of tissues in upper panel is 0.6 mm. CD163 is a marker for M2-polarized macrophages.

 
Go


View this table:
[in this window]
[in a new window]

 
Table A1. Summary of Patient Information

 
Go


View this table:
[in this window]
[in a new window]

 
Table A2. Primer Sequences Used in Quantitative Reverse Transcription Polymerase-Chain Reaction

 
Go


View this table:
[in this window]
[in a new window]

 
Table A3. Cox Survival Models Used for Individual Data Sets

 
Go


View this table:
[in this window]
[in a new window]

 
Table A4. Cox Analysis Using G1(S4) Clusters in the Temozolomide- and Radiotherapy-Treated Patient Cohort (n = 42)

 
Go


View this table:
[in this window]
[in a new window]

 
Table A5. Cluster G07

 
Go


View this table:
[in this window]
[in a new window]

 
Table A6. Cluster G13

 
Go


View this table:
[in this window]
[in a new window]

 
Table A7. Cluster G14

 
Go


View this table:
[in this window]
[in a new window]

 
Table A8. Cluster G18

 
Go


View this table:
[in this window]
[in a new window]

 
Table A9. Cluster G24

 
Go


View this table:
[in this window]
[in a new window]

 
Table A10. Cluster G25

 
Go


View this table:
[in this window]
[in a new window]

 
Table A11. G28 Versus G98 Probe Sets

 
Go


View this table:
[in this window]
[in a new window]

 
Table A12. Estimated Hazard Ratios for Treatment Effect in MGMT-Methylated Patients at Different Expression Levels of G98: Multivariate Cox Proportional Hazards Model (n = 68)

 
Go


View this table:
[in this window]
[in a new window]

 
Table A13. Mouse Probe Sets Used for GSEA in the Krivtsov Data Set: G98 Mouse Homologues

 
Go


View this table:
[in this window]
[in a new window]

 
Table A14. Human Probe Sets Used for GSEA in the Ross Data Set: G98 U133A Corresponding Probe Sets to U133 Plus 2.0

 
Partial Least Square Regression
Partial least square (PLS) regression is a technique that combines features from principal component analysis and multiple linear regression (Wold H: Partial Least Squares. New York, NY Wiley, 1985). It is particularly useful to predict an outcome from a large set of highly correlated predictors. In this study, the PLS procedure in SAS was used to define combinations of the genes expressions (ie, factors) that attempt to explain the genes expression variability and the survival outcome at the same time. For each patient, the martingale residuals obtained from a Cox regression with a constant as the only predictor were used as outcome variable to be explained (Wold H: Partial Least Squares. New York, NY Wiley, 1985). This approach was derived from the works of Therneau et al (Biometrika 77:147-160, 1990) and Leblanc et al (Biometrics 48:411-425, 1992) with the classification and regression tree algorithm and applied to PLS regression. This allowed us to account for the effect of censoring on survival estimates.

The optimal number of PLS factors is generally obtained by cross validation, but in this small sample (n = 42 or 68), this method was not conclusive (ie, no optimal number of factors could be obtained). Therefore, we empirically chose to fit a PLS model with two factors to minimize overfitting and obtain interpretable results. The two factors explained 66% of the survival outcome variations. The contribution of each gene was assessed by their coefficients in the linear regression and their X weights plotted in a factorial plane. Arbitrarily, clusters were given X weights by averaging the X weights of their constituent genes.

Description of External Data Sets
Two external gene expression data sets (GeneChip U133A&B, Affymetrix) of malignant glioma were used for validation of survival factors, referred to as "Nelson" (Freije WA, Castro-Vargas FE, Fang Z, et al: Cancer Res 64:6503-6510, 2004) and "Aldape" data sets (Phillips HS, Kharbanda S, Chen R, et al: Cancer Cell 9:157-173, 2006), respectively. Only cases with diagnosis of glioblastoma, astrocytoma grade III, and mixed oligoastrocytoma grade III were included. The Nelson data set utilized here comprised 17 grade III and 54 grade IV gliomas. The median age was 43 years at diagnosis for all 71 patients, or 51 years for glioblastoma patients, respectively. The Aldape data set consisted of 75 patients, 27 grade III, 48 grade IV with a median age at diagnosis of 47 years for all cases, or 49 years for glioblastoma patients, respectively. This data set was enriched for long-term survivors according to the authors. Adjustments in the Cox models for these two external glioma validation sets comprised age as a dichotomous variable (> 50 years) and tumor grade.

In addition, we used a mouse expression data set published by Krivtsov et al (Krivtsov AV, Twomey D, Feng Z, et al: Nature 442:818-822, 2006) and a human leukemia data set by Ross et al (Ross ME, Mahfouz R, Onciu M, et al: Blood 104:3679-3687, 2004). For both human malignant glioma validation sets, only the chips that passed quality criteria based on the BioConductor library affyPLM (Gentleman R, Carey V, Huber W, et al: Bioinformatics and Computational Biology Solutions Using R and Bioconductor, Springer, 2005) were included. The external data sets were each separately normalized using robust multiarray average, when Affymetrix CEL files were available (Aldape, Nelson, and Ross). For the Krivtsov mouse expression data set (Krivtsov AV, Twomey D, Feng Z, et al: Nature 442:818-822, 2006), the Affymetrix Microarray Suite version 5.0 processed data, as submitted in the Gene Expression Omnibus database, were used. For the human-murine comparison, we mapped the human gene sets to homolog murine gene sets using the annotation provided by Affymetrix on http//www.affymetrix.com (accessed November 14, 2006; Table A13). For GSEA analysis, the data sets were filtered to retain only the most variant probe sets (standard deviation, > 0.5).

Supervised Analysis
For each covariate, the hazard ratio, its 95% CI, and the associated Wald P were examined. The Benjamini-Hochberg procedure (Benjamini Y, Hochberg Y: J R Stat Soc B 57:289-300, 1995) was applied as multiple testing correction method to convert P values into estimated false-discovery rates. The HRs for a gene (or the mean of a cluster) were standardized by referring them to a variation corresponding to an interquartile range, to obtain an interpretable value, for microarray as well as for real-time PCR data. Therefore, the interpretation of an estimated HR of 1.5 is that the hazard rate increases by 50% for an increase of 1 interquartile range of the log-scale gene expression, independently of the expression level at which it is calculated.


    Glossary Terms
 TOP
 ABSTRACT
 INTRODUCTION
 PATIENTS AND METHODS
 RESULTS
 DISCUSSION
 AUTHORS' DISCLOSURES OF...
 AUTHOR CONTRIBUTIONS
 Appendix
 Glossary Terms
 REFERENCES
 

HOX (homeobox) genes:
The proteins encoded by these genes, the homeodomain proteins, play important roles in the developmental processes.

MGMT:
The DNA repair enzyme, O6-methylguanine DNA methyltransferase, which confers resistance to alkylating agents. Thus, cells are protected from the toxicity of alkylating agents, which frequently target the O6 position of guanine in DNA.

Clustering:
Organization of data consisting of many variables (multivariate data) into classes with similar patterns. Hierarchical clustering creates a dendrogram based on pairwise similarities in gene expression within a set of samples. Samples within a cluster are more similar to one another than to samples outside the cluster. The vertical length of branches in the tree represents the extent of similarity between the samples. Thus, shorter the branch length, the fewer the differences between the samples.

Tissue microarray:
Used to analyze the expression of genes of interest simultaneously in multiple tissue samples, tissue microarrays consist of hundreds of individual tissue samples placed on slides ranging from 2 to 3 mm in diameter. Using conventional histochemical and molecular detection techniques, tissue microarrays are powerful tools to evaluate the expression of genes of interest in tissue samples. In cancer research, tissue microarrays are used to analyze the frequency of a molecular alteration in different tumor type, to evaluate prognostic markers, and to test potential diagnostic markers.

EGFR (epidermal growth factor receptor):
Also known as HER-1, EGFR belongs to a family of receptors (HER-2, HER-3, HER-4 are other members of the family) and binds to the EGF, TGF-{alpha}, and other related proteins, leading to the generation of proliferative and survival signals within the cell. It also belongs to the larger family of tyrosine kinase receptors and is generally overexpressed in several solid tumors of epithelial origin.

CD133:
CD133, encoded by the Prominin-1 gene, is used as a marker for stem cells of normal tissues such as neural and hematopoietic stem cells, but also as a marker for tumor stem-like cells of distinct origin such as brain tumors and breast and colon cancer.

Cancer stem-like cells:
A cancer cell that has the potential to transfer disease or to form tumors after transplantation. Cancer stem-like cells have the potential to self-renew (see self-renewal), forming additional tumorigenic cancer cells of similar phenotype, and to give rise to phenotypically diverse cancer cells with more limited potential.

Self-renewal:
Process by which a progenitor gives rise to daughter progenitors of equivalent developmental potential (ie, multipotent stem cells self-renew by dividing to generate one or two multipotent daughter cells).

Tumor-host interaction:
The evolving cross talk between different cell types within the tumor and its host environment. Cancer cells can subvert the normal tissue architecture and reprogram the activity of stromal cells to their advantage, such as promoting the formation of new blood vessels or the recruitment of inflammatory cells, which are key events promoting cancer cell survival, invasion, and metastasis.


    ACKNOWLEDGMENTS
 
We thank A.-C. Diserens, P. Descombes, and D. Chollet for excellent technical support; D.R. Macdonald, A. Guha, and A. Merlo for providing frozen tumor tissues; and the clinical colleagues and nurses involved in the studies at the participating hospitals.


    NOTES
 
Supported by the Swiss National Science Foundation (Grant No. 3100AO-108266/1; M.E.H.); Jacqueline Seroussi Memorial Foundation for Cancer Research (M.E.H.), Translational Research Fund of the EORTC (M.E.H.); the Nélia and Amadeo Barletta Foundation (M.E.H., R.S.); National Center of Competence in Research (NCCR) Molecular Oncology at the Swiss Institute for Experimental Cancer Research (ISREC) (M.E.H., M.D.); Ridgefield Foundation (E.D.); EC FP6 (E.D.); and the VITAL-IT project of the Swiss Institute of Bioinformatics (E.M., M.D.).

A.M. and E.M. contributed equally to this article.

Terms in blue are defined in the glossary, found at the end of this article and online at www.jco.org.

Authors’ disclosures of potential conflicts of interest and author contributions are found at the end of this article.


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 PATIENTS AND METHODS
 RESULTS
 DISCUSSION
 AUTHORS' DISCLOSURES OF...
 AUTHOR CONTRIBUTIONS
 Appendix
 Glossary Terms
 REFERENCES
 
1. Stupp R, Mason WP, van den Bent MJ, et al: Radiotherapy plus concomitant and adjuvant temozolomide for glioblastoma. N Engl J Med 352:987-996, 2005[Abstract/Free Full Text]

2. Stupp R, Dietrich P, Ostermann Kraljevic S, et al: Promising survival for patients with newly diagnosed glioblastoma multiforme treated with concomitant radiation plus temozolomide followed by adjuvant temozolomide. J Clin Oncol 20:1375-1382, 2002[Abstract/Free Full Text]

3. Hegi ME, Diserens AC, Godard S, et al: Clinical trial substantiates the predictive value of O-6-methylguanine-DNA methyltransferase promoter methylation in glioblastoma patients treated with temozolomide. Clin Cancer Res 10:1871-1874, 2004[Abstract/Free Full Text]

4. Hegi ME, Diserens AC, Gorlia T, et al: MGMT gene silencing and benefit from temozolomide in glioblastoma. N Engl J Med 352:997-1003, 2005[Abstract/Free Full Text]

5. Stupp R, Hegi ME, Gilbert MR, et al: Chemoradiotherapy in malignant glioma: Standard of care and future directions. J Clin Oncol 25:4127-4136, 2007[Abstract/Free Full Text]

6. Reya T, Morrison SJ, Clarke MF, et al: Stem cells, cancer, and cancer stem cells. Nature 414:105-111, 2001[CrossRef][Medline]

7. Bao S, Wu Q, McLendon RE, et al: Glioma stem cells promote radioresistance by preferential activation of the DNA damage response. Nature 444:756-760, 2006[CrossRef][Medline]

8. Bonnet D, Dick JE: Human acute myeloid leukemia is organized as a hierarchy that originates from a primitive hematopoietic cell. Nat Med 3:730-737, 1997[CrossRef][Medline]

9. Krivtsov AV, Twomey D, Feng Z, et al: Transformation from committed progenitor to leukaemia stem cell initiated by MLL-AF9. Nature 442:818-822, 2006[CrossRef][Medline]

10. Singh SK, Hawkins C, Clarke ID, et al: Identification of human brain tumour initiating cells. Nature 432:396-401, 2004[CrossRef][Medline]

11. Getz G, Domany E: Coupled two-way clustering server. Bioinformatics 19:1153-1154, 2003[Abstract/Free Full Text]

12. Getz G, Levine E, Domany E: Coupled two-way clustering analysis of gene microarray data. Proc Natl Acad Sci U S A 97:12079-12084, 2000[Abstract/Free Full Text]

13. Blatt M, Wiseman S, Domany E: Superparamagnetic clustering of data. Phys Rev Lett 76:3251-3254, 1996[CrossRef][Medline]

14. Freije WA, Castro-Vargas FE, Fang Z, et al: Gene expression profiling of gliomas strongly predicts survival. Cancer Res 64:6503-6510, 2004[Abstract/Free Full Text]

15. Phillips HS, Kharbanda S, Chen R, et al: Molecular subclasses of high-grade glioma predict prognosis, delineate a pattern of disease progression, and resemble stages in neurogenesis. Cancer Cell 9:157-173, 2006[CrossRef][Medline]

16. Ross ME, Mahfouz R, Onciu M, et al: Gene expression profiling of pediatric acute myelogenous leukemia. Blood 104:3679-3687, 2004[Abstract/Free Full Text]

17. Benjamini Y, Hochberg Y: Controlling the false discovery rate: A practical and powerful approach to multiple testing. J R Stat Soc B 57:289-300, 1995

18. Clement V, Sanchez P, de Tribolet N, et al: HEDGEHOG-GLI1 signaling regulates human glioma growth, cancer stem cell self-renewal, and tumorigenicity. Curr Biol 17:165-172, 2007[CrossRef][Medline]

19. Vairapandi M, Balliet AG, Hoffman B, et al: GADD45b and GADD45g are cdc2/cyclinB1 kinase inhibitors with a role in S and G2/M cell cycle checkpoints induced by genotoxic stress. J Cell Physiol 192:327-338, 2002[CrossRef][Medline]

20. Abate-Shen C: Deregulated homeobox gene expression in cancer: Cause or consequence? Nat Rev Cancer 2:777-785, 2002[CrossRef][Medline]

21. Abdel-Fattah R, Xiao A, Bomgardner D, et al: Differential expression of HOX genes in neoplastic and non-neoplastic human astrocytes. J Pathol 209:15-24, 2006[CrossRef][Medline]

22. Subramanian A, Tamayo P, Mootha VK, et al: Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A 102:15545-15550, 2005[Abstract/Free Full Text]

23. Manley GT, Fujimura M, Ma T, et al: Aquaporin-4 deletion in mice reduces brain edema after acute water intoxication and ischemic stroke. Nat Med 6:159-163, 2000[CrossRef][Medline]

24. Herrlich A, Leitch V, King LS: Role of proneuregulin 1 cleavage and human epidermal growth factor receptor activation in hypertonic aquaporin induction. Proc Natl Acad Sci U S A 101:15799-15804, 2004[Abstract/Free Full Text]

25. Saadoun S, Papadopoulos MC, Hara-Chikuma M, et al: Impairment of angiogenesis and cell migration by targeted aquaporin-1 gene disruption. Nature 434:786-792, 2005[CrossRef][Medline]

26. St Croix B, Rago C, Velculescu V, et al: Genes expressed in human tumor endothelium. Science 289:1197-1202, 2000[Abstract/Free Full Text]

27. Madden SL, Cook BP, Nacht M, et al: Vascular gene expression in nonneoplastic and malignant brain. Am J Pathol 165:601-608, 2004[Abstract/Free Full Text]

28. Chute JP, Muramoto GG, Whitesides J, et al: Inhibition of aldehyde dehydrogenase and retinoid signaling induces the expansion of human hematopoietic stem cells. Proc Natl Acad Sci U S A 103:11707-11712, 2006[Abstract/Free Full Text]

29. Piccirillo SG, Reynolds BA, Zanetti N, et al: Bone morphogenetic proteins inhibit the tumorigenic potential of human brain tumour-initiating cells. Nature 444:761-765, 2006[CrossRef][Medline]

30. Calabrese C, Poppleton H, Kocak M, et al: A perivascular niche for brain tumor stem cells. Cancer Cell 11:69-82, 2007[CrossRef][Medline]

31. Mantovani A, Sozzani S, Locati M, et al: Macrophage polarization: Tumor-associated macrophages as a paradigm for polarized M2 mononuclear phagocytes. Trends Immunol 23:549-555, 2002[CrossRef][Medline]

32. Van Ginderachter JA, Movahedi K, Hassanzadeh Ghassabeh G, et al: Classical and alternative activation of mononuclear phagocytes: Picking the best of both worlds for tumor promotion. Immunobiology 211:487-501, 2006[CrossRef][Medline]

33. Hussain SF, Yang D, Suki D, et al: The role of human glioma-infiltrating microglia/macrophages in mediating antitumor immune responses. Neuro-oncol 8:261-279, 2006[Abstract/Free Full Text]

34. Dorrance AM, Liu S, Yuan W, et al: Mll partial tandem duplication induces aberrant Hox expression in vivo via specific epigenetic alterations. J Clin Invest 116:2707-2716, 2006[CrossRef][Medline]

35. Krivtsov AV, Armstrong SA: MLL translocations, histone modifications and leukaemia stem-cell development. Nat Rev Cancer 7:823-833, 2007[CrossRef][Medline]

36. Bachoo RM, Maher EA, Ligon KL, et al: Epidermal growth factor receptor and Ink4a/Arf: Convergent mechanisms governing terminal differentiation and transformation along the neural stem cell to astrocyte axis. Cancer Cell 1:269-277, 2002[CrossRef][Medline]

37. Stupp R, Hegi ME: Targeting brain-tumor stem cells. Nat Biotechnol 25:193-194, 2007[CrossRef][Medline]

38. Sathornsumetee S, Rich JN: New treatment strategies for malignant gliomas. Expert Rev Anticancer Ther 6:1087-1104, 2006[CrossRef][Medline]

39. Jain RK: Normalization of tumor vasculature: An emerging concept in antiangiogenic therapy. Science 307:58-62, 2005[Abstract/Free Full Text]

40. Batchelor TT, Sorensen AG, di Tomaso E, et al: AZD2171, a pan-VEGF receptor tyrosine kinase inhibitor, normalizes tumor vasculature and alleviates edema in glioblastoma patients. Cancer Cell 11:83-95, 2007[CrossRef][Medline]

41. Kjellman C, Olofsson SP, Hansson O, et al: Expression of TGF-beta isoforms, TGF-beta receptors, and SMAD molecules at different stages of human glioma. Int J Cancer 89:251-258, 2000[CrossRef][Medline]

Submitted December 11, 2007; accepted March 11, 2008.


Add to CiteULike CiteULike   Add to Complore Complore   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us   Add to Digg Digg   Add to Facebook Facebook   Add to Reddit Reddit   Add to Technorati Technorati   Add to Twitter Twitter    What's this?


This article has been cited by other articles:


Home page
Clin. Cancer Res.Home page
J. Felsberg, M. Rapp, S. Loeser, R. Fimmers, W. Stummer, M. Goeppert, H.-J. Steiger, B. Friedensdorf, G. Reifenberger, and M. C. Sabel
Prognostic Significance of Molecular Markers and Extent of Resection in Primary Glioblastoma Patients
Clin. Cancer Res., November 1, 2009; 15(21): 6683 - 6693.
[Abstract] [Full Text] [PDF]


Home page
Neuro OncologyHome page
H. Colman, L. Zhang, E. P. Sulman, J. M. McDonald, N. L. Shooshtari, A. Rivera, S. Popoff, C. L. Nutt, D. N. Louis, J. G. Cairncross, et al.
A multigene predictor of outcome in glioblastoma
Neuro Oncology, October 20, 2009; (2009) nop007v1.
[Abstract] [Full Text] [PDF]


Home page
Ann OncolHome page
B. Neyns, J. Sadones, E. Joosens, F. Bouttens, L. Verbeke, J.-F. Baurain, L. D'Hondt, T. Strauven, C. Chaskis, P. In't Veld, et al.
Stratified phase II trial of cetuximab in patients with recurrent high-grade glioma
Ann. Onc., September 1, 2009; 20(9): 1596 - 1603.
[Abstract] [Full Text] [PDF]


Home page
Mol Cancer ResHome page
G. Frosina
DNA Repair and Resistance of Gliomas to Chemotherapy and Radiotherapy
Mol. Cancer Res., July 1, 2009; 7(7): 989 - 999.
[Abstract] [Full Text] [PDF]


Home page
Acta Biochim Biophys SinHome page
W. Zhuang, Z. Qin, and Z. Liang
The role of autophagy in sensitizing malignant glioma cells to radiation therapy
Acta Biochim Biophys Sin, May 1, 2009; 41(5): 341 - 351.
[Abstract] [Full Text] [PDF]


Home page
Mol Cancer ResHome page
M. Ropolo, A. Daga, F. Griffero, M. Foresta, G. Casartelli, A. Zunino, A. Poggi, E. Cappelli, G. Zona, R. Spaziante, et al.
Comparative Analysis of DNA Repair in Stem and Nonstem Glioma Cell Cultures
Mol. Cancer Res., March 1, 2009; 7(3): 383 - 392.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Purchase Article
Right arrow View Shopping Cart
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a colleague
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Save to my personal folders
Right arrow Download to citation manager
Right arrowRights & Permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Murat, A.
Right arrow Articles by Hegi, M. E.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Murat, A.
Right arrow Articles by Hegi, M. E.
Social Bookmarking
 Add to CiteULike   Add to Complore   Add to Connotea   Add to Del.icio.us   Add to Digg   Add to Facebook   Add to Reddit   Add to Technorati   Add to Twitter  
What's this?

About
JCO
 Editorial
Roster
 Advertising
Information
 Librarians &
Institutions
 Rights &
Permissions
 PDA Services

Copyright © 2008 by the American Society of Clinical Oncology, Online ISSN: 1527-7755. Print ISSN: 0732-183X
Terms and Conditions of Use
  HighWire Press HighWire Press™ assists in the publication of JCO Online