
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Imaging, Diagnosis, Prognosis |
Authors' Affiliations: Departments of 1 Human Genetics, 2 Radiological Sciences, 3 Biostatistics, 4 Neurology, 5 Molecular and Medical Pharmacology, 6 Neurological Surgery, and 7 Pharmacology and Laboratory Medicine, University of California at Los Angeles and David Geffen School of Medicine at the University of California at Los Angeles, Los Angeles, California
Requests for reprints: Whitney B. Pope, Department of Radiological Sciences, University of California at Los Angeles, 10833 LeConte Avenue, BL-153/CHS, Los Angeles, CA 90095-1271. Phone: 310-794-7923; Fax: 310-206-5958; E-mail: wpope{at}mednet.ucla.edu.
| Abstract |
|---|
|
|
|---|
Experimental Design: Using DNA microarray analysis, we examined VEGF and related gene expression in 71 newly diagnosed malignant gliomas and analyzed the relationship to edema and survival.
Results and Conclusions: VEGF expression was predictive of survival in tumors with little or no edema [Cox proportional hazard model, 6.88; 95% confidence interval (95% CI), 2.61-18.1; P < 0.0001], but not in tumors with extensive edema. The expression of several proangiogenic genes, including adrenomedullin (correlation coefficient, 0.80), hypoxia-inducible factor-1A (0.51), and angiopoietin-2 (0.44), was correlated with VEGF expression (all with P < 0.0001), whereas that of several antiangiogenic genes was inversely correlated. The expression of six genes was increased greater than 3-fold in edematous versus nonedematous tumors in the absence of increased VEGF expression. The most increased, neuronal pentraxin 2 (NPTX2, 7-fold change), was predictive of survival in tumors with the highest levels of edema, in contrast to VEGF (hazard ratio, 2.73; 95% CI, 1.49-5.02; P = 0.049). NPTX2 was tightly correlated with expression of the water channel aquaporin-3 (0.74, P < 0.0001). These results suggest that there are both VEGF-dependent and VEGF-independent pathways of edema production in gliomas and may explain why edema is not reduced in some patients following anti-VEGF treatment.
Several imaging features of GBM, including edema, negatively correlate with survival (4). Edema is also associated with shortened survival in grade III gliomas. The vascular endothelial growth factor (VEGF) is a key mediator of tumor angiogenesis and edema. Multiple VEGF isoforms are generated by alternative splicing, although the first five exons are conserved across all variants (5). Through the activation of receptor tyrosine kinases, VEGF impacts endothelial cell permeability, activation, survival, proliferation, invasion, and migration, all of which play a significant role in tumor progression. VEGF is also a potent mediator of vascular permeability, being 50 times more effective than histamine (6). Malignant gliomas express both VEGF and its receptors (7), and glioblastoma cell lines have been shown to secrete VEGF (8). Anti-VEGF therapy has been shown to improve survival in patients with metastatic colorectal cancer, breast cancer, and lung cancer (6).
Microarray analysis offers the ability to assess gene expression and to correlate these data with imaging features, including edema. Although many glioma patients show dramatic reduction in tumor and edema following anti-VEGF therapy, other patients show no response (9, 10). Therefore, we sought to use microarray analysis to investigate the possibility of VEGF-independent pathways linked to edema and survival in patients with malignant glioma.
| Materials and Methods |
|---|
|
|
|---|
Microarray data. Using the TRIzol reagent (Invitrogen Life Technologies), total RNA was extracted from the tumor samples, and processed using an RNeasy column (Qiagen). cDNA and cRNA were generated using standard protocols (11). All samples were processed, scanned, and quality checked as previously described (12).
For analysis of gene expression measures, affymetrix data were normalized using the justRMA method provided by the Bioconductor group (13). Most samples were processed on the HG-U133A arrays. Recent samples were processed on the HG-U133 2.0 arrays. These two array pools were normalized as above. It was observed that the average brightness of these two pools was similar but not identical. Thus, the smaller pool was uniformly scaled to match the average brightness of the median brightness array from the larger group. Because the HG-U133A platform is a subset of the HG-U133 2.0 platform, only data from probe sets that were shared between the two platforms were used.
Probe sets designed to detect VEGF expression were assessed for their ability to represent VEGF expression. All four representative probe sets behaved similarly for all tests. The probe set with the greatest sensitivity as evidenced by the largest coefficient of variation was selected to represent the results presented here.
To confirm VEGF expression levels, 25 tumor samples spanning the range of VEGF expression values were analyzed with real-time PCR. The ratio of VEGF to actin was calculated, and the Pearson correlation coefficient was determined.
Statistical analysis. The Kaplan-Meier method was used to estimate the survival distributions (14). Log-rank tests were used to test the difference between stratified survival groups. To assess which covariates affect survival, we used multivariable Cox proportional hazard models (15). Hazard ratios correspond to risk of death, and thus, an increased hazard ratio implies a worse prognosis. The proportional hazard assumption was tested using scaled Schoenfeld residuals (16). For each covariate, the relative hazard rate and the associated P value were examined. For all analyses, a P value of <0.05 was accepted as significant. Statistical analyses were carried out with the freely available software packages R.8 The relationship between edema and VEGF expression levels was cross-validated with the use of a Kruskal-Wallis test. Univariate differences in covariates were tested across categorical groupings by using the Kruskal-Wallis test (17). Distributions of covariates across categorical groupings were visualized with box plots.
| Results |
|---|
|
|
|---|
Expression of VEGF in grade III and IV gliomas. VEGF expression was highly variable in GBM. The range of VEGF expression was 78 to 5,584 (n = 52), with a mean of 1,288 and SD of 1,206. For grade III gliomas (n = 19), the mean VEGF expression was 313 with a SD of 541. Thus, VEGF expression was
4.1 times higher in GBM than grade III gliomas (P < 0.0001), in agreement with a previous report (18). For all malignant glioma, the median expression of VEGF was 508. Of the tumors with greater than this level of VEGF expression (n = 35), 33 were GBM.
Correlation of VEGF expression with edema. Figure 1 illustrates grading of edema, with examples of edema grades 0, 1, and 2 categorized per Materials and Methods. For grade III tumors, 79% were edema grade 0, 16% were edema grade 1, and 5% were edema grade 2. For GBM, those percentages were 23%, 23%, and 54%, respectively. For all malignant gliomas, VEGF expression was correlated with edema (Kruskal-Wallis test, P = 0.0011). When subdivided into edema grades, the mean expression of VEGF in grade 0 (n = 27) was 479 (SD, 582; range, 64-1,873), mean expression of VEGF in grade 1 (n = 15) was 1,373 (SD, 1182; range, 84-3,605), and mean expression of VEGF in grade 2 (n = 29) was 1,359 (SD, 1,353; range 94-5,584). Note that some tumors had high levels of edema, without elevated VEGF expression. There was a large range of VEGF expression within the different edema subgroupings, particularly for tumors with grade 2 edema (Fig. 2 ).
|
|
|
|
Expression of several antiangiogenic genes was anticorrelated with VEGF expression: HIF-1
subunit inhibitor (HIF-1AN; –0.54, P < 0.0001), brain-specific angiogenesis inhibitor 1 (BAI1; –0.45, P < 0.0001), phosphatase and tensin homologue (PTEN; –0.36, P < 0.0001), BAI3 (–0.34, P < 0.003), somatostatin (–0.31, P < 0.008), and BAI2 (–0.25, P = 0.034).
Survival analysis showed that for tumors with edema <2, HIG-2, angiopoetin-2, IGFBP-5, IGFBP-7, neuropilin, and RBBP-8 were all significant predictors of survival, whereas HIF-1A was not. HIF-1A was significant only for the edema grade 1 group. Kaplan-Meier analysis segregated by edema grades shows the similarities between the survival curves of several of the proangiogenic genes and that of VEGF (Fig. 3).
Genes associated with edema in low VEGF tumors. Some tumors with low VEGF expression were found to have large edema volumes, suggesting VEGF-independent edema pathways in malignant gliomas. To generate a list of potential genes involved in producing edema without VEGF involvement, we analyzed the subset of tumors with low VEGF (<400), comparing expression levels in tumors with edema versus tumors without edema (see Table 1 ). Neuronal pentraxin 2 (NPTX2) was identified as the gene in which expression was most increased (7-fold) in the edematous versus nonedematous tumors. NPTX2 was positively correlated with edema in all malignant gliomas (0.24, P = 0.04; Kruskall-Wallis P = 0.0008; Fig. 2) and within the GBM subset (Kruskall-Wallis P = 0.0009; Fig. 2). NPTX2 was not significantly correlated with VEGF expression. To assess the effects of NPTX2 expression on survival, we again used the Kaplan-Meier method, sorting tumors by edema grade. Unlike VEGF, NPTX2 was shown to be a negative predictor of survival in the highest edema grade (edema grade 2, Fig. 3). NPTX2 was also significant for survival in all gliomas, as well as the GBM subset. The hazard ratio for NPTX2 using the Cox model was 2.73 (P = 0.049; 95% CI, 1.49-5.02). Both NPTX2 (P = 0.024) and edema (P < 0.0001) were independent predictors of survival when used together in a Cox multivariable proportional hazard model. These results suggest that NPTX2 may mediate edema independent of VEGF expression. The other genes with greater than 3-fold increase in the edematous tumors were the tissue inhibitor of metalloproteinase 1, matrix Gla protein, CD163 antigen, cyclic AMP-regulated phosphoprotein 21 kDa, and neurofilament 3. The expression of these genes did not show a significant correlation with survival in the high-edema tumors.
|
Gene expression and survival in GBM. As stated above, we found a significant relationship between survival and VEGF expression in low-edema tumors, whereas high-edema tumors showed a relationship between survival and both NPTX2 and aquaporin 3 expression. To determine if these relationships were affected by tumor grade, we also analyzed these data for the subset of GBM patients (n = 52; Fig. 4). Kaplan-Meier analysis showed a relationship between survival and VEGF in GBM patients with little or no edema (P = 0.04). GBM patients with grade 2 edema also showed a significant relationship between survival and both NPTX2 (P = 0.014) and aquaporin 3 (P = 0.026) expression as shown for all malignant glioma patients. Lastly, in a Cox multivariable proportional hazard model using VEGF and NPTX2, both were independent predictors of survival (VEGF, P = 0.011; NPTX2, P = 0.016). For all gliomas, these values are P < 0.0001 for VEGF and P = 0.004 for NPTX2.
Relationship between steroid dose, edema, and gene expression. Corticosteroids are used to diminish edema and mass effect in brain tumor patients. Steroid dose was determined for patients at the time of the MRI used for edema grading by reviewing the medical records. Of the 71 patients, 34 were not on steroids, and 34 were on steroids (22 on decadron, doses from 4 to 60 mg/day, 10 on decadron with unknown dose, 1 on prednisone, and 1 on dexamethasone), and no data on steroids were available for 3 patients. Of the patients on steroids, 16 had edema grade 2 and, thus, did not have edema levels diminished by steroid administration. There was no significant difference in steroid use by edema grade for the patients with known steroid dose (n = 56; Kruskal-Wallis test, P = 0.37). There was no significant relationship between steroid use and VEGF (P = 0.16), NPTX2 (P = 0.22), or aquaporin 3 (P = 0.91) expression. We also analyzed Kaplan-Meier survival data in patients who were not on steroids (n = 34). For edema grade 2, the relationship between survival and NPTX2 maintained significance (P = 0.0049), although aquaporin 3 did not (P = 0.15), most likely because of the reduced statistical power. When edema grade 2 patients were added to the no-steroid group (i.e., the subset of patients in which steroids were not given or had no effect on edema grade), aquaporin 3 expression maintained a significant relationship to survival (P = 0.018). For VEGF, the relationship with survival for the edema <2 groups showed a low P value (P = 0.055) and showed a similar survival curve for all gliomas with edema <2. Therefore, we conclude that differences in steroid treatment do not account for the relationship between VEGF, NPTX2, and aquaporin 3 expression and survival.
| Discussion |
|---|
|
|
|---|
In the current study, we used microarray data to show a large range of VEGF expression in malignant gliomas. We found that VEGF expression was correlated with edema. However, some tumors with a large amount of peritumoral edema had low levels of VEGF expression. This is consistent with our hypothesis that some gliomas cause edema through VEGF-independent pathways, potentially accounting for failure of these tumors to respond to anti-VEGF therapy.
In addition to being a potent stimulator of vascular permeability, VEGF is thought to play a pivotal role in angiogenesis in malignant gliomas (19, 22). We found that VEGF expression levels were predictive of survival independent of edema in tumors with little or no edema. This was true both for malignant gliomas and for only GBM. This supports the conclusion that whereas part of the deleterious effect of VEGF expression may be due to the promotion of edema and mass effect, other properties of VEGF, presumably its role in angiogenesis, may also contribute to shortened survival.
To investigate the mechanism underlying VEGF-induced angiogenesis, we assessed the correlation between VEGF expression and several putative pro- and antiangiogenic genes. Of these, adrenomedullin was the most tightly correlated with VEGF expression. This was not unexpected because both adrenomedullin and VEGF are thought to be induced by HIF (19), and HIF was also correlated with VEGF expression. Hypoxia seems crucial in the transformation of glial tumors into more aggressive phenotypes (reviewed in ref. 19). We found that another gene product associated with hypoxia, HIG-2, was also well correlated with VEGF expression. Others have shown a link between HIG-2 and renal cell carcinoma (23). No reports implicating HIG-2 in glioma progression have been published.
Several antiangiogenesis genes were inversely correlated with VEGF expression, including HIF-1AN, BAI1-3, PTEN, and somatostatin. Interestingly, PTEN has been shown to inhibit adrenomedullin expression in glioma cell lines (24) and to increase VEGF expression in a mouse model for asthma (25). PTEN is also associated with chromosome 10q deletion and a poor prognosis (26). We analyzed the relationship between edema, survival, and the expression of these angiogenesis-related genes as we had done for VEGF. Several were significant predictors of survival, with many showing the same pattern as VEGF, that is, associated with survival in tumors with little or no edema. Correlation between VEGF expression and the expression of these angiogenesis-related genes, as well as their similar survival curves, suggests a functional relationship that merits further investigation. These findings also support the idea that angiogenesis is a result of the balance between pro- and antiangiogenic gene products and identifies a subset of genes that may impact survival and be more closely linked to VEGF-mediated angiogenesis in human gliomas.
Given our hypothesis that some malignant gliomas have VEGF-independent pathways leading to edema, we looked for gene expression enriched in edematous versus nonedematous tumors with low VEGF expression. This led to the identification of NPTX2, which was increased 7-fold in the edematous tumor group. This was the greatest percentage increase of all genes analyzed. NPTX2 is normally expressed in the central nervous system and is a member of a family of proteins related to C reactive protein and other acute-phase inflammatory mediators (27, 28) and, thus, is a good candidate for mediation of tumoral edema. NPTX2 was found to be correlated with edema in all gliomas and in only GBM patients. NPTX2 was not correlated with VEGF expression, supporting a VEGF-independent role for this gene product. Most importantly, increased NPTX2 was associated with poorer survival in tumors with the highest levels of edema. This is the reverse of VEGF expression, which was predictive of survival only for tumors with little or no edema. It will be of interest to determine if NPTX2 is elevated in tumors that do not respond to anti-VEGF therapy.
We also searched for genes that correlate with NPTX2 expression, looking for candidates that could be involved in edema production. This led to the identification of aquaporin 3, which was the gene with expression levels that most highly correlated (0.74, P < 0.0001) with those of NPTX2. Aquaporin 3 is a member of a group of transmembrane proteins that act as water and solute channels and, thus, could potentially be involved in edema regulation (reviewed in ref. 29). Others have suggested a relationship between the expression of other aquaporin isoforms and edema in traumatic brain injury (30). Additionally, aquaporin-1 has been implicated in increased metastatic potential in a mouse tumor model (31) and increased in astrocytomas (32). More investigation will be required to understand the possible role of aquaporin 3 in brain edema in tumor patients, as well as its functional relationship to NPTX2. It is notable, however, that aquaporin 3 is induced by epidermal growth factor in fibroblast culture (33), and amplification of the epidermal growth factor receptor has been implicated in glioma progression (34, 35).
In conclusion, although VEGF expression is correlated with edema, it is an independent predictor of survival, presumably due to its angiogenic properties. Of the many putative pro- and antiangiogenic genes, we have found several that seem closely related to VEGF expression, suggesting these are genes are involved in VEGF-mediated angiogenesis. Because VEGF expression has a greater impact on survival in tumors with little or no edema, our data support the hypothesis that VEGF may be crucial in the "angiogenic switch" that occurs when tumors degenerate into more aggressive, edema-producing phenotypes. How these genetic factors interplay with prognostically important clinical data, including age, treatment, and Karnofsky performance scale, remains to be determined, but this is a focus of ongoing investigation. Our data also suggest that VEGF-independent edema pathways may be important in limiting survival. In this exploratory study, with its limitation of small sample size, we found evidence that NPTX2 and aquaporin 3 are two gene products that merit further investigation as potential mediators of continued high-edema states.
| Footnotes |
|---|
Note: M.R.J. Carlson and W.B. Pope contributed equally to this work.
Received 11/21/06; revised 2/ 7/07; accepted 2/ 7/07.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
A. S. Chi, A. G. Sorensen, R. K. Jain, and T. T. Batchelor Angiogenesis as a Therapeutic Target in Malignant Gliomas Oncologist, June 1, 2009; 14(6): 621 - 636. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. L. Nghiemphu, W. Liu, Y. Lee, T. Than, C. Graham, A. Lai, R. M. Green, W. B. Pope, L. M. Liau, P. S. Mischel, et al. Bevacizumab and chemotherapy for recurrent glioblastoma: A single-institution experience Neurology, April 7, 2009; 72(14): 1217 - 1222. [Abstract] [Full Text] [PDF] |
||||
![]() |
W. B. Pope, J. H. Chen, J. Dong, M. R. J. Carlson, A. Perlina, T. F. Cloughesy, L. M. Liau, P. S. Mischel, P. Nghiemphu, A. Lai, et al. Relationship between Gene Expression and Enhancement in Glioblastoma Multiforme: Exploratory DNA Microarray Analysis Radiology, October 1, 2008; 249(1): 268 - 277. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Diehn, C. Nardini, D. S. Wang, S. McGovern, M. Jayaraman, Y. Liang, K. Aldape, S. Cha, and M. D. Kuo From the Cover: Identification of noninvasive imaging surrogates for brain tumor gene-expression modules PNAS, April 1, 2008; 105(13): 5213 - 5218. [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 | Meeting Abstracts Online |