| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Molecular Biology and Genetics |
Pathway Identified in Childhood Astrocytomas by Angiogenesis Gene Profiling1 ,,2
Center for Cancer Research, Childrens Research Institute [S. K., K. M. P., C. L., T. J. M.], Research Center for Genetic Medicine, Childrens Research Institute [K. M. B., D. D., D. A. S.], and Department of Pathology [M. R. S.], Childrens National Medical Center, Washington, DC 20010; Department of Preventative Medicine, Division of Biostatistics, Vanderbilt University, Nashville, Tennessee 37232 [B. L.]; and Graduate Program in Genetics, George Washington University, Washington, DC 20052 [K. M. B.]
| ABSTRACT |
|---|
|
|
|---|
as well as high-level expression of FK506-binding protein (FKBP) 12 by HGA. Furthermore, 9 of 21 (43%) genes overexpressed by HGA were HIF/FKBP-associated genes. This group included the epidermal growth factor receptor (EGFR), which promotes HIF synthesis, as well as insulin-like growth factor-binding protein 2 (IGFBP2), a target gene of HIF activity. Differential protein expression of HIF-2
was validated in an independent group of 16 astrocytomas (P = 0.02). We conclude that the EGFR/FKBP12/HIF-2
pathway is important in childhood HGA and represents a potential new therapeutic target. | INTRODUCTION |
|---|
|
|
|---|
Studies demonstrating the crucial role of angiogenesis in cancer have been a major advance in our understanding of malignant tumor progression (3) . One of the key histopathological features that distinguish HGA from LGA is intense, increased angiogenesis. The invasiveness of HGA, another unique feature of this tumor in comparison with LGA, is associated with increased microvascular density and intratumoral hypoxia (4) . Thus, inhibitors of hypoxia-inducible angiogenesis factors could be important new therapeutic agents against HGA. In adults, the most commonly described regulators of brain tumor-derived angiogenesis are VEGF, platelet-derived growth factor, angiopoietin-2, and their respective receptors (5) . It is not known to what extent these same regulatory mechanisms exist in pediatric astrocytomas, especially because there are key differences in the common genetic abnormalities observed between adult and childhood astrocytomas that may greatly affect the development of angiogenesis. For instance, although p53 mutations appear to occur with equal frequency (2838%) in childhood and adult HGA, the frequency of these mutations in HGA of children younger than 4 years of age is significantly less (6) . Furthermore, the incidence of PTEN mutations in adult HGA is approximately 30%, whereas the same changes in childhood HGA are rare (7) . This is particularly important because studies have shown that wild-type p53 and PTEN expression inhibits angiogenesis induced by malignant astrocytoma cells (8 , 9) .
In this study we used oligonucleotide microarray profiling (12,626 probe sets) to identify angiogenesis-related genes predictive of childhood HGAs to gain a better understanding of the regulation of malignant angiogenesis in childhood HGA and to identify potential novel therapeutic targets.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Gene Expression Profiling.
Total RNA was isolated from the tissues and cell lines using TRIzol reagent as described previously (10)
. Briefly, a minimum of 1025 µg of total RNA was obtained from each tumor tissue sample and cell line. Ten µg of extracted RNA from each sample were synthesized into double-stranded cDNA with the SuperScript choice system (Life Technologies, Inc., Carlsbad, CA) using an oligo(dT) primer containing a T7 RNA polymerase promoter (Genset, San Diego, CA). The double-stranded cDNA was purified by phenol-chloroform extraction and in vitro transcribed using the ENZO Bioarray RNA transcript labeling kit (Affymetrix, Santa Clara, CA). The biotin-labeled cRNA was purified using the RNeasy kit (Qiagen, Valencia, CA) and then fragmented to approximately 200 bp by alkaline treatment [200 mM Tris-acetate (pH 8.2), 500 mM potassium acetate, and 150 mM magnesium acetate]. Each cRNA was hybridized to an Affymetrix U95Av2 array containing approximately 12,000 human transcripts and detected using a streptavidin-conjugated fluor and confocal laser scanning according to the manufacturers recommendations (Affymetrix). Data were extracted using Microarray Suite Version 4.0 (Affymetrix) and linearly scaled such that the average intensity across each chip equaled 800. Data from a chip were rejected if the ratio of the average intensity of the 3' to 5' ends of the human glyceraldehyde-3-phosphate dehydrogenase gene exceeded 3, or if the linear scaling factor exceeded 4. Genes not called "present" for any of the samples, as determined by the Affymetrix software, were excluded from further consideration. To avoid division by zero or negative numbers in the calculation of average fold differences between sample classes, values of 10 or less were set to 10. Subsequent statistical methods used the complete data set in which values less than 10 were not adjusted. All raw expression values can be downloaded from online supplementary data.2
Gene ontology annotations for genes represented on the U95Av2 chip were compiled from the Affymetrix NetAffx analysis website.6
RT-PCR Analysis.
PCR amplification was performed directly on the T7-dT primered double-stranded cDNA (originally generated for expression profiling) from three LGAs and three HGAs. Each cDNA was diluted 100-fold, and 5 µl of the diluted cDNA were used as template for each PCR. ß-Actin was used as internal control. Gene-specific primers from Invitrogen Life Technologies, Inc. [EGFR, 5'-tgcagcgatacagctcagaccc-3' (forward primer) and 5'-tttgggaacggactggtttatg-3' (reverse primer); IGFBP2, 5'-tgacaagcatggcctgtacaac-3' (forward primer) and 5'-ggatcagcttcccggtgttg-3' (reverse primer); HIF-2
, 5'-agcctccatctgccatcagtc-3' (forward primer) and 5'-cttgccatgcctgacaccttg-3' (reverse primer); and ß-actin, 5'-acgtgtacgcc-aacacagtg-3' (forward primer), 5'-gccatgccaatctcatctt-3' (reverse primer)] were used in a 50-µl reaction consisting of 1x PCR Gold Buffer, 2.25 mM MgCl2, 0.2 mM deoxynucleoside triphosphates, 2 units of TaqGold polymerase, and 12 pmol of each primer. Initial denaturation (94°C, 12 min) was followed by 30 cycles of denaturation (94°C, 30 s), annealing (57°C, 30 s), and extension (72°C, 45 s). This was followed by a final extension step (72°C, 5 min). Four µl of ß-actin product were combined with 8 µl of EGFR, IGFBP2, or HIF-2
product for each sample; electrophoresed on a 1.5% agarose gel stained with ethidium bromide; and visualized with UV light. Product bands were quantitated by volume densitometry using Quantity One software (Bio-Rad, Hercules, CA), and signals were normalized to actin.
Selection of Genes Involved in Angiogenesis.
The expression profile of 133 angiogenesis-related genes was specifically selected for analysis from the 9,198 total probe sets with detectable expression in the astrocytoma samples of the 12,626 total probe sets analyzed on the array. For genes represented on the array with more than one probe set, the probe set that was called "present" in the greatest number of samples by the Microarray Suite Version 4.0 software was chosen as the "best" representative for that gene. If the number of present calls was equal for more than one probe set, then the probe set with the highest average numerical expression value across all of the samples was chosen as the "best" representative. Designation of known genes associated with angiogenesis was based on extensive review of available published literature databases.
Hierarchical Clustering and PCA.
As a means to examine data structure, hierarchical clustering and PCA were performed on the complete microarray dataset (9,198 probe sets called present in at least one sample) as well as on the limited dataset of 133 preselected angiogenesis-specific genes. Hierarchical clustering dendrograms, as well as the graphical representation of the data for the statistically significant angiogenesis-related genes, were generated using downloaded Cluster and TreeView software7
(written by Michael Eisen). All clustering was performed using the uncentered-correlation and complete linkage clustering options. Data represented in clustergrams are normalized intensities. Two-dimensional PCA plots were generated using the PCA scores of the two components accounting for the greatest intersample variability as calculated by the Cleaver 1.0 website (Stanford University).8
Statistical Methods and Selection of Differentially Expressed Genes.
We used permutational Ps less than or equal to 0.05 to select our candidate gene list. Our permutation tests were performed to determine which genes showed significant differences in mean intensity between the classes. To determine a permutational P for each gene, we first calculated the mean and variances for the sample intensities within each class and calculated a t test statistic (
). The class assignments of the tumors were subsequently permuted, using sampling without replacement, 50,000 times. For each permutation we calculated a test statistic
*(i) for all i = 1, ..., 50,000 permutations. The permutational P is the probability of observing at least as large a value for the test statistic for any random permutation [
*(i)] as for our initial test statistic (
). This method has the benefit of not depending on distributional assumptions.
Protein Validation by Immunohistochemistry.
Sixteen formalin-fixed paraffin-embedded astrocytoma tissues (8 HGAs and 8 LGAs) as well as one normal control cortical brain were cut into 4-µm-thick sections. Standard immunohistochemistry was performed as described previously in detail (10)
. Antibodies for HIF-1
were obtained from US Biologicals (Swampscott, MA), and antibodies for HIF-2
were obtained from Novus Biologicals (Littleton, CO). Three authors independently reviewed the slides and graded the staining as absent (-), weak (+), moderate (++), and strong (+++). The staining was then correlated with the grades of tumors (WHO I-IV) and with RNA expression levels obtained by the microarray experiments. Significance of differential protein expression was evaluated by comparing the number of HGA and LGA specimens with positive (+, ++, and +++) and negative (-) staining using the one-tailed Fischers exact test.
| RESULTS |
|---|
|
|
|---|
, independently confirmed the microarray expression results (Fig. 2)
|
|
To explore the possibility that such a structure may exist in the smaller subset of 133 genes chosen for a role in angiogenesis, we performed hierarchical clustering and PCA on this subset of data (Fig. 1, C and E)
. As with the larger dataset, the two tumor classes clustered into separate dendrogram branches. PCA also distinguished the tumor classes using only a single component (principal component 1, explaining 24.2% of the variability). The most frequently occurring ontologies in the categories of biological process, cell component, and molecular function are summarized in the online supplementary data.
Immunohistochemical Staining for HIF Proteins.
We performed immunohistochemical protein validation for HIF-1
and HIF-2
in an independent set of 16 tumor tissues (8 HGAs and 8 LGAs) as well as a section of normal cortical brain. Representative images are shown in Fig.3
. Absent to minimal staining was observed for HIF-1
in normal brain (Fig. 3A)
as well as in LGA (Fig. 3B)
and HGA (Fig. 3C)
. Although there is some background staining, there is no definite cellular staining. In contrast, positive staining was observed for HIF-2
in astrocytoma specimens. Staining was observed primarily in HGA specimens (Fig. 3, F and G)
, often around the leading invasive tumor edges. Comparison of the number of positively (+, ++, and +++ combined) and negatively (-) staining slides for HIF-2
showed significant overexpression of HIF-2
protein in HGA in comparison with LGA (P = 0.02), as shown in Fig.3H
.
|
| DISCUSSION |
|---|
|
|
|---|
On the microarray, 9198 genes had above-threshold expression in at least one astrocytoma. Of these genes, 2401 (26%) were differentially expressed between HGA and LGA with a permutational P < 0.05. Thus, there is a marked difference in the expression patterns of the HGA and LGA specimens tested, and the largest proportion of differential genes is involved in signal transduction and gene transcription.
Using all probe sets (9198) expressed in the tumor samples, high- and low-grade tumors separated into different dendrogram branches by hierarchical clustering. Additionally, PCA using this same gene set yielded separation of the samples using only the one component explaining the largest proportion of variability in the dataset. The separation of the tumor classes using such unsupervised methods is not surprising, given the differences in expression of a large proportion of genes between malignant and benign astrocytomas. More striking were the considerable differences in the expression of the 133 angiogenesis-related genes that were chosen for examination a priori. Using only this gene list, hierarchical clustering and PCA separated the tumor classes into two distinct groups in a manner similar to that observed in the unsupervised methods using all genes. This indicates that the angiogenesis gene expression profile alone is an accurate predictor of pathological grade in astrocytoma. Not surprisingly, a large proportion of the angiogenesis-related genes (36%) showed marked expression differences between sample classes.
By focusing analysis on angiogenesis genes, we identified overexpression of HIF-2
in HGA for the first time in childhood astrocytoma. HIF-2
encodes a basic helix-loop-helix PAS domain protein transcription factor that regulates vascular maturation, remodeling, and stabilization of angiogenesis in response to hypoxia (12)
. Unlike the related transcription factor HIF-1
, which was found to be similarly expressed in HGA and LGA, HIF-2
was differentially overexpressed by HGA at the mRNA level (P = 0.003). This differential expression was confirmed to be significant at the protein level as well (P = 0.02). Expression of HIF-1
has been shown in adult HGA and oligodendroglioma (13)
. HIF-2
is strongly expressed in breast carcinoma and is associated with high vascular grade and poor survival (14)
.
The synthesis of HIF-2
is induced via a signal transduction pathway initiated by receptor tyrosine kinases, such as the EGFR. This leads to activation of phosphatidylinositol 3'-kinase, protein kinase B (AKT), and FRAP, also known as mammalian target for rapamycin (15)
. We show that EGFR mRNA is differentially up-regulated in HGA (P = 0.003). This is a consistent finding in adult HGA (16)
. Furthermore, we found multiple FKBPs, especially FKBP12, to be highly expressed in LGA, HGA, and astrocytoma cells. More importantly, the largest group of functionally related, differentially overexpressed genes is comprised of members of the EGFR/FKBP12/HIF-2
pathway, constituting nearly one-fourth of all of the up-regulated angiogenesis genes in HGA. Interestingly, the phosphatidylinositol 3'-kinase regulatory unit upstream of FKBP12, p85, is down-regulated in HGA, suggesting dysregulation of this pathway via EGFR overexpression. These data implicate the EGFR/FKBP12/HIF-2
pathway as an important promoter of angiogenesis in childhood HGA.
Expression of HIF-2
results in the induction of several target genes such as Glut-1, VEGF, and PGK1 (17)
. These genes are responsible for angiogenic and metabolic changes resulting in increased microvascular density and tumor growth. Glut-1 was not differentially expressed in the astrocytomas analyzed. However, PGK1 was highly expressed in HGA, LGA, and both cell lines and was differentially overexpressed in HGA (P = 0.01). PGK1 is secreted by tumor cells and reduces disulfide bonds in the serine proteinase plasmin to release the microvessel inhibitor angiostatin (18)
. Although VEGF was highly expressed in HGA, LGA, and astrocytoma cells, it was not differentially expressed. This is in contrast to that seen in adults (19)
. These data suggest that HIF-induced expression of PGK1, previously undescribed in astrocytoma, may be critical to the promotion of the malignant angiogenic phenotype in childhood HGA.
Other genes typically associated with angiogenesis in adult HGA, such as platelet-derived growth factor, angiopoietin-2, and their receptors, were also highly expressed yet not significantly different in expression levels between LGA and HGA. Interestingly, our study shows that 24% of all of the preselected angiogenesis-related genes that were differentially expressed between HGA and LGA belonged to the families of extracellular matrix and proteases such as collagen, matrix metalloproteinases, and urokinase.
Fig. 4
diagrams the pathway of EGFR/FKBP12/HIF-2
activation and highlights the large proportion of genes in this pathway that were found to be overexpressed by HGA in our study. Taken together, these findings could have direct therapeutic importance. For example, FK506 and rapamycin are structurally related immunosuppressants that block distinct stages in intracellular signaling pathways by binding to members of the FRAP/FKBP12 immunophilin family and down-regulating signals through S6 kinase and FRAP, affecting the expression level of mRNA important for progression from G1 to S phase (20)
. Because hypoxia enhances angiogenesis through FRAP/FKBP12 signaling, blocking this pathway may have potent antineoplastic effects. Thus, our findings raise the distinct possibility that rapamycin analogues and HIF-2 inhibitors could be useful as therapeutic agents in pediatric HGA.
|
| FOOTNOTES |
|---|
1 Supported by NIH Grant K08 CA92056-01A1, an Avery Research Scholarship from Childrens Research Institute, and the Childhood Brain Tumor Foundation [T. J. M. and K. M. P.] and by National Institute of Neurological Disorders and Stroke Grant 1U24NS43571-01 [to D. A. S. and K. M. B.]. ![]()
2 Supplementary data for this article are available at Cancer Research Online (http://cancerres.aacrjournals.org). ![]()
3 These authors contributed equally to this work. ![]()
4 To whom requests for reprints should be addressed, at Department of Hematology/Oncology, Center for Cancer Research, Childrens National Medical Center, 111 Michigan Avenue NW, Washington, DC 20010. Phone: (202) 884-2146; Fax: (202) 884-5685; E-mail: tmacdona{at}cnmc.org ![]()
5 The abbreviations used are: LGA, low-grade astrocytoma; HGA, high-grade astrocytoma; VEGF, vascular endothelial growth factor; HIF, hypoxia-inducible factor; PCA, principal components analysis; FKBP, FK506-binding protein; EGFR, epidermal growth factor receptor; FRAP, FK506-binding protein rapamycin-associated protein; PGK1, phosphoglycerate kinase 1; RT-PCR, reverse transcription-PCR. ![]()
8 http://classify.stanford.edu/. ![]()
Received 9/12/02. Accepted 2/14/03.
| REFERENCES |
|---|
|
|
|---|
in brain tumors: association with angiogenesis, invasion and progression. Cancer (Phila.), 88: 2606-2618, 2000.[Medline]
(HIF-2
) expression in tumor-infiltrative macrophages to tumor angiogenesis and the oxidative thymidine phosphorylase pathway in human breast cancer. Cancer Res., 62: 1326-1329, 2002.This article has been cited by other articles:
![]() |
H. K. Thorarinsdottir, M. Santi, R. McCarter, E. J. Rushing, R. Cornelison, A. Jales, and T. J. MacDonald Protein Expression of Platelet-Derived Growth Factor Receptor Correlates with Malignant Histology and PTEN with Survival in Childhood Gliomas Clin. Cancer Res., June 1, 2008; 14(11): 3386 - 3394. [Abstract] [Full Text] [PDF] |
||||
![]() |
V. Anelli, C. R. Gault, A. B. Cheng, and L. M. Obeid Sphingosine Kinase 1 Is Up-regulated during Hypoxia in U87MG Glioma Cells: ROLE OF HYPOXIA-INDUCIBLE FACTORS 1 AND 2 J. Biol. Chem., February 8, 2008; 283(6): 3365 - 3375. [Abstract] [Full Text] [PDF] |
||||
![]() |
S.-H. Kim, H. Nakagawa, A. Navaraj, Y. Naomoto, A. J.P. Klein-Szanto, A. K. Rustgi, and W. S. El-Deiry Tumorigenic Conversion of Primary Human Esophageal Epithelial Cells Using Oncogene Combinations in the Absence of Exogenous Ras Cancer Res., November 1, 2006; 66(21): 10415 - 10424. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. L. Covello, J. Kehler, H. Yu, J. D. Gordan, A. M. Arsham, C.-J. Hu, P. A. Labosky, M. C. Simon, and B. Keith HIF-2{alpha} regulates Oct-4: effects of hypoxia on stem cell function, embryonic development, and tumor growth. Genes & Dev., March 1, 2006; 20(5): 557 - 570. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Aghi, P. Gaviani, J. W. Henson, T. T. Batchelor, D. N. Louis, and F. G. Barker II Magnetic Resonance Imaging Characteristics Predict Epidermal Growth Factor Receptor Amplification Status in Glioblastoma Clin. Cancer Res., December 15, 2005; 11(24): 8600 - 8605. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Avellino, S. Romano, R. Parasole, R. Bisogni, A. Lamberti, V. Poggi, S. Venuta, and M. F. Romano Rapamycin stimulates apoptosis of childhood acute lymphoblastic leukemia cells Blood, August 15, 2005; 106(4): 1400 - 1406. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. W. Vogel, Z. Zhuang, J. Li, H. Okamoto, M. Furuta, Y.-S. Lee, W. Zeng, E. H. Oldfield, A. O. Vortmeyer, and R. J. Weil Proteins and Protein Pattern Differences between Glioma Cell Lines and Glioblastoma Multiforme Clin. Cancer Res., May 15, 2005; 11(10): 3624 - 3632. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. L. Covello, M. C. Simon, and B. Keith Targeted Replacement of Hypoxia-Inducible Factor-1{alpha} by a Hypoxia-Inducible Factor-2{alpha} Knock-in Allele Promotes Tumor Growth Cancer Res., March 15, 2005; 65(6): 2277 - 2286. [Abstract] [Full Text] [PDF] |
||||
![]() |
K.-K. Wong, Y.-M. Chang, Y. T.M. Tsang, L. Perlaky, J. Su, A. Adesina, D. L. Armstrong, M. Bhattacharjee, R. Dauser, S. M. Blaney, et al. Expression Analysis of Juvenile Pilocytic Astrocytomas by Oligonucleotide Microarray Reveals Two Potential Subgroups Cancer Res., January 1, 2005; 65(1): 76 - 84. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. R. Rhodes, J. Yu, K. Shanker, N. Deshpande, R. Varambally, D. Ghosh, T. Barrette, A. Pandey, and A. M. Chinnaiyan Large-scale meta-analysis of cancer microarray data identifies common transcriptional profiles of neoplastic transformation and progression PNAS, June 22, 2004; 101(25): 9309 - 9314. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. K. Lee, A. K. Hsu, J. Sajdak, J. Qin, and P. Pavlidis Coexpression Analysis of Human Genes Across Many Microarray Data Sets Genome Res., June 1, 2004; 14(6): 1085 - 1094. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. Hopfl, O. Ogunshola, and M. Gassmann HIFs and tumors--causes and consequences Am J Physiol Regulatory Integrative Comp Physiol, April 1, 2004; 286(4): R608 - R623. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Broniscer and A. Gajjar Supratentorial High-Grade Astrocytoma and Diffuse Brainstem Glioma: Two Challenges for the Pediatric Oncologist Oncologist, April 1, 2004; 9(2): 197 - 206. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Godard, G. Getz, M. Delorenzi, P. Farmer, H. Kobayashi, I. Desbaillets, M. Nozaki, A.-C. Diserens, M.-F. Hamou, P.-Y. Dietrich, et al. Classification of Human Astrocytic Gliomas on the Basis of Gene Expression: A Correlated Group of Genes with Angiogenic Activity Emerges As a Strong Predictor of Subtypes Cancer Res., October 15, 2003; 63(20): 6613 - 6625. [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 |