|
|
||||||||
Basic Science Investigations |
1 Clinical Cooperation Unit Nuclear Medicine, German Cancer Research Center, Heidelberg, Germany
2 Institute of Immunology, University Rostock, Rostock, Germany
3 Department I, University Orthopedic Clinic, Heidelberg, Germany
| ABSTRACT |
|---|
|
|
|---|
Key Words: giant cell tumor 18F-FDG quantification gene expression
| INTRODUCTION |
|---|
|
|
|---|
The differential diagnosis includes osteosarcoma and metastatic carcinoma with giant cell tumor as well as benign diseases, such as nonossifying fibroma and benign fibrous histiocytoma. Giant cell tumors of the spinal region are especially difficult to differentiate from other lesions. Tubbs et al. emphasized that the metastases noted in their study also had no distinguishing radiologic features (1).
Some studies have been performed to assess the role of gene expression in giant cell tumors. Tian et al. examined the expression of ADAM12, a member of the disintegrin and metalloprotease family, in 20 giant cell tumors and found that ADAM12 is involved in the fusion process for mononuclear stromal cells in giant cell tumors (2). Roux et al. evaluated the receptor activator of nuclear factor
B and noted that the expression of the gene may have a role in the pathogenesis of giant cell tumors (3).
The roles of vascular endothelial growth factor (VEGF) and matrix metalloproteinase 9 in the extent of bone destruction and the likelihood of local tumor recurrence were evaluated by Kumta et al. for patients with bone tumors, including 14 patients with giant cell tumors (4). The authors emphasized that recurrent tumors showed higher expression of VEGF and matrix metalloproteinase 9. One limitation of some published studies is that the evaluations were usually confined to a few genes. This problem can be solved by the use of gene chip technology, which allows the assessment of a large number of genes simultaneously. Dimitrakopoulou-Strauss et al. reported preliminary results obtained with gene chip analysis for tumor patients and demonstrated a correlation of cyclin A and VEGF-C with 18F-FDG uptake (5). Furthermore, a correlation among glucose transporters, oncogenes, and genes related to multidrug resistance was observed (6).
The purpose of this study was to assess 18F-FDG kinetics in giant cell tumors and to compare 18F-FDG kinetic data with data from a gene chip analysis of tumor specimens obtained by surgery shortly after a PET study. The major aims were to obtain quantitative data about 18F-FDG kinetics in these tumors and to compare 18F-FDG kinetics with gene expression.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Dynamic PET studies were performed for 60 min after the intravenous application of 300370 MBq of 18F-FDG by use of a 23-frame protocol (10 frames of 1 min, 5 frames of 2 min, and 8 frames of 5 min). A dedicated PET system (ECAT EXACT HR+; Siemens Co.) with an axial field of view of 15.3 cm, operated with septa extended (2-dimensional mode), was used for the PET examinations. The system allows the simultaneous acquisition of 63 transverse slices with a theoretical slice thickness of 2.4 mm. Transmission scans were obtained for a total of 10 min before radiotracer application for attenuation correction of the acquired emission tomographic data. All PET images were attenuation corrected, and an image matrix of 256 x 256 pixels was used for iterative image reconstruction. The reconstructed cross-sectional images were converted to standardized uptake value (SUV) images by use of the formula SUV = tissue concentration (Bq/g)/(injected dose [Bq]/body weight [g]) (7).
Dynamic PET data were evaluated by use of the software package PMod (provided through cooperation with the University of Zurich) (8,9). Visual analysis was performed by evaluating the metabolically active areas on transaxial, coronal, and sagittal images. Timeactivity curves were created by use of volumes of interest (VOIs). Generally, a VOI consists of several regions of interest (ROIs) over the target area. Irregular ROIs were drawn manually. To compensate for possible patient motion during the acquisition time, the original ROIs were visually repositioned but were not redrawn.
Quantitative evaluation of tracer kinetics requires compartment modeling. A 2-tissue compartment model is the generally accepted standard model for the quantification of dynamic 18F-FDG examinations (10,11). One problem with patient studies is accurate measurement of the input function, which theoretically requires arterial blood sampling. However, the input function can be retrieved from the image data with good accuracy (12). For the input function, the mean value of the VOI data obtained from a large arterial vessel was used. A vessel VOI consisted of at least 7 ROIs in sequential PET images. The recovery coefficient was 0.85 for a diameter of 8 mm and for the system described above. For lesions located in the arms or in the spine, a VOI over the aorta was used to obtain input data. For lesions located in the legs, partial volume correction was applied when a small vessel had to be used to retrieve the input data. Noise in the input curve has an effect on parameter estimates. Therefore, we used a preprocessing tool that is available in the PMod software and that allows a fit of the input curve by summing up to 3 decaying exponential functions to reduce noise.
The constants k1k4 (18F-FDG transport into cells, 18F-FDG transport out of cells, phosphorylation of intracellular 18F-FDG, and dephosphorylation of intracellular 18F-FDG, respectively) were calculated by use of a 2-compartment model implemented in the PMod software, taking into account the vascular fraction (VB) in a VOI as an additional variable. Details about the applied compartment model are described by Burger and Buck (8). One advantage of the PMod software is the graphical interface, which allows interactive configuration of the kinetic model by the user as well as the application of some preprocessing steps, for example, setting up initial values and limits for the fit parameters. Visual evaluation of each plot was performed to check the quality of the fit. Each model curve was compared with the corresponding timeactivity curve, and the total
2 difference was used as the cost function, where the criterion was to minimize the summed squares (
2) of the differences between the measured curve and the model curve (
2 was usually less than 1). This means that the squared residual values (measured value minus estimated value) are multiplied by weights. In theory, the weight should be related to the SE of a measurement. The distribution at each individual point is taken to be gaussian, with an SD to be specified. The residual covariance was dependent on the kinetic parameter and typically was less than 10% for k1. The model parameters were accepted when k1k4 were less than 1 and the VB values exceeded 0. The unit for rate constants k1k4 is 1 min1, whereas VB reflects the fraction of blood within the evaluated volume.
Besides compartment analysis, we applied a noncompartment model to the kinetic data to obtain the fractal dimension. The software module applies the box plot model to the kinetic data and calculates the fractal dimension of the timeactivity curve. The fractal dimension values vary from 0 to 2 and are parameters for a deterministic or more chaotic distribution of tracer activity over time. No input function is needed for the fractal dimension model. Besides the kinetic analysis, the SUV in the last frame, 5660 min after tracer injection, was used for statistical evaluation.
Surgery was performed shortly after the PET examination, and tumor specimens were obtained from 5 patients. The surgeons were aware of the PET results and removed a specimen from the target area already evaluated with PET. The tumor tissue was transported in liquid nitrogen, and total RNA was extracted for further processing by gene chip analysis. RNA extraction was done by use of an RNeasy kit (Qiagen, Inc.). Thereafter, the quality of isolated RNA was evaluated photometrically by use of the ratio of the measurement at 280 nm to the measurement at 260 nm and on an agarose gel. We used a U95A gene chip (Affymetrix Inc.), which provides quantitative information about more than 12,000 gene sequences. RNA concentrations were measured spectrophotometrically at 260 nm. RNA probes were labeled according to the suppliers instructions (Affymetrix) and hybridized as described by Teague et al. (13). Specifically, first-strand synthesis was performed by use of a T7.(dT)24 primer and SuperScript II/ribonuclease H reverse transcriptase (Invitrogen GmbH) with a total RNA sample of 10 µg. Second-strand synthesis was performed by use of Escherichia coli DNA polymerase I, E. coli DNA ligase, and ribonuclease H. Fragment end polishing was performed by use of T4 polymerase. An in vitro transcription reaction was used to incorporate biotin-11-cytidine triphosphate and biotin-16-uridine triphosphate into the complementary RNA probe (bioarray high-yield RNA transcript labeling kit; Enzo), resulting in linear amplification of the double-stranded complementary DNA template. The fragmented complementary RNA was hybridized overnight (45°C) to HG U95Av2 arrays (Affymetrix). Washes were performed by use of GeneChip Fluidics Station 400 (Affymetrix) according to the manufacturers protocol. Staining was done with R-phycoerythrin-streptavidin (Molecular Probes; MoBiTec) and was followed by an antibody amplification procedure with a biotinylated antistreptavidin antibody (Vector Laboratories) and goat IgG (Sigma-Aldrich Chemie GmbH). Scanning (GeneArray Scanner 2500; Hewlett-Packard GmbH) was performed at a 3-µm resolution with 488-nm excitation and 570-nm emission wavelengths. Intensity values were scaled to a target intensity of 500 by use of MicroArray Suite 5.0 (Affymetrix) software. The scan data analysis was done by use of MicroDB 3.0 and Data Mining Tool 3.0 (Affymetrix); additionally, GeneMaths (Applied Maths BVBA) was used for cluster analysis later on.
Gene chip expression data were normalized to the expression of the ß2-microglobulin gene (Affymetrix code 34644_at; Homo sapiens messenger RNA for ß2-microglobulin) by use of the following formula: relative expression value = 1,000 x (expression value for a gene/expression value for the ß2-microglobulin gene). Analyses of both PET data and gene chip data were performed by use of dedicated software programmed by our group by using Delphi 7 Studio (Borland Software Corp.). The program provides a correlation analysis of all gene expression data and corresponding PET data (SUV, VB, k1k4, and fractal dimension).
Statistical evaluation of PET data was performed by using a Statistica software package (version 6.0; StatSoft Co.) on a personal computer (Pentium IV [Intel Corp.] 2.4-GHz processor and 1 GB of random-access memory) running with Windows XP Professional SP1 (Microsoft Co.). The program package was used to calculate the basic statistics for kinetic data. Furthermore, the program package was also used to perform regression analysis for PET data and individual gene expression data.
| RESULTS |
|---|
|
|
|---|
|
|
|
|
|
|
|
| DISCUSSION |
|---|
|
|
|---|
In addition to radiographic examinations, PET with 18F-FDG may be helpful for gaining diagnostic information by assessing tumor metabolism. However, only limited information is available about the use of 18F-FDG for giant cell tumors. Kern et al. was one of the first groups to describe the uptake of 18F-FDG in a cystic giant cell tumor of the femur (20). The authors evaluated 5 patients with bone and soft-tissue tumors and noted for giant cell tumors a glucose utilization rate comparable to that of neurofibromas and lower than that of grade I liposarcomas. Dimitrakopoulou-Strauss et al. evaluated patients with soft-tissue sarcomas and reported a median SUV of 0.7, with a range of 0.32.2, for 10 patients with grade I liposarcomas (21). We noted a median SUV of 4.6 for the 19 giant cell tumors that we studied; this SUV is significantly higher than that found by those authors for grade I liposarcomas. Interestingly, the median value for k1 was 1.8 times higher in our study than the k1 reported for soft-tissue sarcomas; this result indicates significantly higher 18F-FDG transport in giant cell tumors than in sarcomas. In a study of 37 malignant bone tumors, Wu et al. reported a median SUV of 3.1 (range, 0.44612.290), which is more comparable to our data for giant cell tumors (22). Aoki et al. evaluated the SUV for 52 primary benign and malignant bone tumors, including 5 giant cell tumors; they reported a mean SUV of 4.6, which is in excellent agreement with our finding of a mean SUV of 4.8 (23). However, those authors performed only SUV measurements, and no kinetic data were provided. According to our experience with the assessment of 18F-FDG kinetics in malignant tumors, mainly VB and k1 contribute to 18F-FDG uptake in giant cell tumors.
Several studies have been performed to evaluate gene coexpression. Ito et al. evaluated the coexpression of GLUT-1 and matrix metalloproteinase 2 (MMP-2) in 8 human cancer cell lines and 80 human cancer specimens (24). They concluded that MMP-2 expression and cell invasiveness are associated with GLUT-1 expression in human cancer cell lines. Interestingly, only 45 human tumor specimens (56%) showed the coexpression of GLUT-1 and MMP-2, while 19 expressed neither GLUT-1 nor MMP-2. However, with regard to an expected coexpression fraction of only 56%, we would not assume a high correlation for GLUT-1 and MMP-2. Indeed, in a review of our gene chip data, we noted a correlation of only 0.36 for GLUT-1 and MMP-2; this value was not significant. However, we noted a high, significant correlation of 0.93 between MMP-2 expression and the vascular fraction (VB) calculated from the PET kinetic data. This result suggests a link between the VB and MMP-2 activity.
Angiogenesis and glucose transport as well as the role of hypoxia have been studied by several authors. Pedersen et al. assessed 18F-FDG uptake as well as the expression of glucose transporters and VEGF in 2 small-cell lung cancer cell lines (25). They noted a significant upregulation of VEGF-A and GLUT-1 mRNAs as well as GLUT-1 and GLUT-3 proteins during hypoxia. 18F-FDG uptake changed according to the GLUT-1 protein level. Burke et al. evaluated the change in gene expression during hypoxia in human macrophages (26). They used macrophages and a complementary DNA array and found that the mRNA levels were upregulated for several genes, including those for GLUT-1, VEGF-A, matrix metalloproteinase 7, neuromedin B receptor, and the DNA-binding protein inhibitor Id2 (26). The data suggest that 1 parameter, such as hypoxia in tumors, may have an impact on a group of genes, resulting in different effects on the biologic properties of a tumor lesion.
Overall, the literature results provided evidence for a linkage of glucose metabolism, glucose transporter expression, angiogenesis-related gene expression, and hypoxia. The results that we obtained for the 5 giant cell tumors that we studied demonstrated that 18F-FDG uptake, as measured by SUV, was correlated with VEGF-A expression. Furthermore, the fractal dimension for the 18F-FDG kinetic data also was associated with VEGF-A expression, and a very high correlation coefficient was noted for these parameters (r = 0.9707). The correlation of VEGF-A and GLUT-1 was 0.83, which was significant at a P value of <0.1. However, our results showed that the vascular fraction (VB) was correlated with both VEGF-A (r = 0.86; P < 0.1 [significant]) and VEGF-B (r = 0.87; P < 0.1 [significant]) but not with VEGF-C. These results are in agreement with the functions of VEGF-A and VEGF-B, which are related to the formation of new vessels, whereas VEGF-C is mainly related to the lymphatic spread of tumors. The data suggest that VEGF-A and VEGF-B have a major impact on 18F-FDG kinetics and that gene expression levels may be classified on the basis of SUV, VB, and fractal dimension, obtained by quantitative analysis of 18F-FDG kinetics, as predictor variables.
Interestingly, sFLT-1, which binds to VEGF-A and acts as a VEGF-A antagonist, was correlated with k3 (Fig. 4A), possibly because sFLT-1 and VEGF-A may show simultaneous enhancement of gene expression in tumor cells. This finding was reported by Inoue et al. for human hematopoietic cell lines (27). They noted that 16 of 17 cell lines showing enhanced expression of sFLT-1 also showed enhanced expression of VEGF-A. Therefore, the overall effect on tumor angiogenesis is dependent on the impact of all of the parameters. Interestingly, k3 was also associated with HIF-1. HIF-1 expression is usually enhanced in cells with hypoxia and is one of the reasons for increased expression of VEGF-A. HIF-1 is known to be linked to the expression of glucose transporters. Our data also show that k3 is dependent on parameters linked to angiogenesis.
Prognostic parameters are generally important for improving therapy management in patients with tumor lesions. Takahashi et al. evaluated VEGF-A expression in gastrointestinal stromal tumors of the stomach and showed that it was an indicator of a poor prognosis (28). Several studies have been performed with giant cell tumors to detect parameters relevant for prognosis. Zheng et al. evaluated giant cell tumors and noted that 3 major isoforms were present in the tumors (29). We noted the highest median values for the expression of VEGF-A splice variant 189 and VEGF-A, whereas VEGF-B expression was 5-fold lower and VEGF-C expression was only 6% of VEGF-A expression. Zheng et al. were able to show that VEGF-A expression was correlated with clinical stage (29). Therefore, the overexpression of VEGF-A is likely to be associated with a more advanced tumor stage. On the basis of the previously discussed correlation of VEGF-A and VEGF-B with 18F-FDG kinetic data, the quantitative evaluation of PET 18F-FDG kinetics can help to detect tumors with a higher likelihood of an advanced stage.
Microvessel counting is used as an indicator of a poor prognosis for several tumor types (30). However, the number of microvessels alone may not be directly associated with VEGF-A expression or a poor prognosis. Marion-Audibert et al. assessed microvascular density in endocrine tumors of the pancreas and noted a significant difference in the survival curves, with an association of shorter survival and lower vessel density (31). Furthermore, no close correlation was observed between vascular density and VEGF-A expression for this tumor type. Those authors concluded that on the basis of their data, other factors may be involved in the regulation of tumor-associated angiogenesis. Sulh et al. evaluated the proliferative index and the vascular density for 7 nonrecurrent and 13 recurrent giant cell tumors (32). They noted no significant difference in the mean values for the groups and concluded that neither parameter is helpful for predicting tumor recurrence. However, the data from Sulh et al. showed that the SDs for both groups were high; therefore, the lack of significance may have been related to the high variability of the parameters in both groups.
Besides VEGF-A, VEGF-B also binds to the receptor VEGFR1, also known as Flt-1 (33). Silvestre et al. evaluated VEGF-B and emphasized that VEGF-B promotes angiogenesis in association with the activation of Akt and eNOS-related pathways (34). We noted a significant correlation of the vascular fraction, as measured by VB, with VEGF-B expression (Fig. 3A). Although the median expression of VEGF-B (relative expression value, 44.5) was lower than that of VEGF-A (relative expression value, 220.5) for the 5 giant cell tumors that we studied, the median expression of VEGF-B splice variant 186 (relative expression value, 339.9) exceeded even that of VEGF-A. The data indicate a major role of both VEGF-A and VEGF-B in angiogenesis for these tumors.
| CONCLUSION |
|---|
|
|
|---|
| FOOTNOTES |
|---|
For correspondence or reprints contact: Ludwig G. Strauss, MD, Medical PET Group-Biological Imaging (E0601), Clinical Cooperation Unit Nuclear Medicine, German Cancer Research Center, Im Neuenheimer Feld 280, D-69120 Heidelberg, Germany.
E-mail: lgs{at}ads-lgs.com
| REFERENCES |
|---|
|
|
|---|
Related articles in JNM:
This article has been cited by other articles:
![]() |
L. G. Strauss, D. Koczan, S. Klippel, L. Pan, C. Cheng, S. Willis, U. Haberkorn, and A. Dimitrakopoulou-Strauss Impact of Angiogenesis-Related Gene Expression on the Tracer Kinetics of 18F-FDG in Colorectal Tumors J. Nucl. Med., August 1, 2008; 49(8): 1238 - 1244. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. J. Beer, S. Lorenzen, S. Metz, K. Herrmann, P. Watzlowik, H.-J. Wester, C. Peschel, F. Lordick, and M. Schwaiger Comparison of Integrin {alpha}v 3 Expression and Glucose Metabolism in Primary and Metastatic Lesions in Cancer Patients: A PET Study Using 18F-Galacto-RGD and 18F-FDG J. Nucl. Med., January 1, 2008; 49(1): 22 - 29. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Dimitrakopoulou-Strauss, P. Hohenberger, U. Haberkorn, H. R. Macke, M. Eisenhut, and L. G. Strauss 68Ga-Labeled Bombesin Studies in Patients with Gastrointestinal Stromal Tumors: Comparison with 18F-FDG J. Nucl. Med., August 1, 2007; 48(8): 1245 - 1250. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Maeda, U. Tateishi, T. Terauchi, C. Hamashima, N. Moriyama, Y. Arai, E. E. Kim, and K. Sugimura Unsuspected Bone and Soft Tissue Lesions Identified at Cancer Screening using Positron Emission Tomography Jpn. J. Clin. Oncol., March 1, 2007; 37(3): 207 - 215. [Abstract] [Full Text] [PDF] |
||||
![]() |
A.M. McKinney, P. Reichert, J. Short, T. Dhurairaj, K. SantaCruz, Z. McKinney, and S. Kieffer Metachronous, Multicentric Giant Cell Tumor of the Sphenoid Bone with Histologic, CT, MR Imaging, and Positron-Emission Tomography/CT Correlation AJNR Am. J. Neuroradiol., November 1, 2006; 27(10): 2199 - 2201. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| JOURNAL OF NUCLEAR MEDICINE TECHNOLOGY | THE JOURNAL OF NUCLEAR MEDICINE |