|
|
||||||||
Clinical Investigations |
1 Division of Oncological Diagnostics and Therapy, German Cancer Research Center, Heidelberg, Germany
2 Department of Orthopedic Clinics, University of Heidelberg, Heidelberg, Germany
3 Department of Nuclear Medicine, University of Zürich, Zürich, Switzerland
| ABSTRACT |
|---|
|
|
|---|
Key Words: 18F-FDG PET bone lesion fractal dimension kinetics
| INTRODUCTION |
|---|
|
|
|---|
The purpose of this study was to assess whether dynamic quantitative 18F-FDG PET studies with noninvasive measurement of the input function, a dedicated data evaluation using the classical 2-tissue-compartment model, and a noncompartment approach can help to improve the differential diagnosis of malignant and benign bone lesions. In particular, we tried to identify those pharmacokinetic parameters of 18F-FDG that may be helpful for an accurate differentiation of malignant and benign bone lesions. Furthermore, we compared the semiquantitative analysis based on a single static measurement with the analysis of the full kinetic 18F-FDG data. In general, a diagnostic procedure should provide a high posterior probability of true-positive results at a low level of prior probability of disease. Because of the dependency of sensitivity and specificity on the prevalence of disease, we applied the Bayesian statistics to the data to assess the gain in information provided by PET.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Data Acquisition
Dynamic PET studies were performed after intravenous application of 300370 MBq 18F-FDG for 60 min (10 frames of 1 min, 5 frames of 2 min, and 8 frames of 5 min). 18F-FDG was prepared according to the method of Toorongian et al. (5). All patients fasted for at least 4 h before PET. The blood glucose level was measured immediately before 18F-FDG application and was within the normal range in all patients.
A dedicated PET system (ECAT EXACT HR+; Siemens, Erlangen, Germany) based on the block detector technology with a craniocaudal field of view of 15.3 cm, operated in septa-extended mode (2-dimensional mode), was used for patient studies. The system allows the simultaneous acquisition of 63 transversal slices with a theoretic slice thickness of 2.4 mm. The system consists of 4 rings, and each of the rings has 72 bismuth germanate detector blocks. A single block detector is divided into an 8 x 8 matrix. The crystal size of a single detector element is 4.39 x 4.05 x 30 mm. Transmission scans for a total of 10 min were obtained with 3 rotating germanium pin sources before the first radionuclide application for the attenuation correction of the acquired emission tomographic images.
The PET data were transferred by file transfer protocol to a subnet server. A webinterface was used to start and distribute the reconstruction tasks on several computer systems, where the reconstruction programs were running in a semiparallel mode (6,7). All PET images were attenuation corrected and reconstructed using an iterative reconstruction algorithm (weighted least-square method, ordered subsets, 4 subsets, 6 iterations, 256 x 256 image matrix) running on Pentium platforms (Pentium III [600 MHz, double processor, 512-MB random access memory]) and Windows NT (Version 4.0 SP6; Microsoft, Redmond, WA).
Data Analysis
Evaluation of the dynamic PET data was performed using the software package PMod, provided by cooperation with the University of Zürich (Zürich, Switzerland) (8,9). Timeactivity curves were created using volumes of interest (VOIs). 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 not redrawn. In general, a detailed quantitative evaluation of tracer kinetics requires the use of compartment modeling. Patlak analysis and a 2-tissue-compartment model are standard methodologies for the quantification of dynamic 18F-FDG studies (10,11). For the basic analysis, we used the semiquantitative approach based on the calculation of a distribution value, for which the term "standardized uptake value" (SUV) was introduced by Strauss and Conti (1): SUV = tissue concentration (MBq/g)/(injected dose [MBq]/body weight [g]). The 55- to 60-min uptake value served for quantification of the 18F-FDG SUV data.
One problem in patient studies is the 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). We performed compartment analysis to gain more information about the tracer distribution. For the input function the mean value of the VOI data obtained from an arterial vessel was used. A vessel VOI consisted of at least 10 ROIs on sequential PET images. The input data were then fitted using a 3-exponential function to reduce noise. In patients with an abdominal or a thoracic lesion, the descending aorta was used for this purpose because the spillover from other organs is low and the descending aorta extends from the upper chest to the lower abdomen. The input function is generally a problem, particularly if the lesions are located in the extremities, where the vessel diameter is relatively small. For lesions located in the extremities (e.g., legs), we used a VOI consisting of at least 10 ROIs over the femoral artery or another arterial vessel, which was well delineated in the field of view. Eight of 37 tumors and 10 of 46 benign lesions were excluded from the final evaluation because of problems with the input function. Low counting rates and a relatively low increase of the 18F-FDG uptake in the vessel VOI in the early phase (up to 3 min after injection) were some of the problems. The recovery coefficient is 0.85 for a diameter of 8 mm and for the system described above. Partial-volume correction was used for small vessels (diameter < 8 mm) on the basis of phantom measurements of the recovery function. The diameter of the vessels was assessed by contrast-enhanced CT or MR images that were available for the PET study. The 18F-FDG influx (Ki) was calculated using the transport rates from the 2-tissue-compartment model according to the following formula: Ki = ([K1 x k3]/[k2 + k3]). The metabolic rate of glucose according to Patlak and Blasberg (10) was not calculated because of the high variation of
. The transport constant K1 and the rate constants k2, k3, and k4 were calculated using a 2-tissue-compartment model based on a method implemented in the PMod software taking into account the vascular fraction (VB) within a VOI. Details of the applied compartment models are described elsewhere (9).
In addition to the compartment analysis, we used a noncompartment model based on the fractal dimension (13). As shown by other investigators, the fractal dimension is a parameter for the heterogeneity. It was recently shown that the fractal dimension is an appropriate procedure to describe the heterogeneity of blood flow in animal models (14). A Java-based module was implemented in the PMod software to calculate the fractal dimension for the timeactivity data (15). The program is based on the box-counting method (13). The values of the fractal dimension vary from 0 to 2 and are a parameter for a deterministic or more chaotic distribution of the tracer activity. For this purpose, no input function is needed.
Statistical evaluation of the data was performed using the Statistica software package (Version 6.0; StatSoft Co., Hamburg, Germany) on a personal computer (Pentium III [600 MHz, double processor, 512-MB random access memory]) running with Windows NT (Version 4.0 SP6). Descriptive statistics and boxwhiskers plots were used for the analysis of the data. The Student t test was applied to all evaluated parameters (SUV, VB, K1, k2, k3, k4, Ki, fractal dimension) to find out which parameters are significant for the differentiation of malignant and benign lesions. Differences were considered significant for P < 0.05. Discriminant analysis (DA) was used to determine the diagnostic accuracy of an 18F-FDG study using all evaluated parameters (SUV, VB, K1, k2, k3, k4, Ki, fractal dimension) with regard to the final histologic diagnosis. The Bayesian theorem was applied to the results to evaluate the information provided by PET with respect to different levels of the prevalence of disease (16). Details about the discriminant function and the Bayesian analysis are described in the Appendix.
| RESULTS |
|---|
|
|
|---|
|
|
|
|
|
|
|
| DISCUSSION |
|---|
|
|
|---|
Dehdashti et al. (18) studied 20 patients with intraosseous lesions and reported that the SUV was helpful in the diagnosis of benign and malignant intraosseous lesions. With a 2.0 cutoff value for the SUV, the sensitivity was 93% and the specificity was 80% in 15 patients with malignant lesions and 5 patients with benign bone lesions. However, benign lesions such as a fibrous dysplasia and infections showed an SUV in the malignant range (false-positive), whereas 3 bone metastases from prostate, breast, and bladder carcinoma did not show enhanced 18F-FDG uptake (false-negative). Schulte et al. (3) reported on 202 histologically verified bone lesions and found a sensitivity of 93% and a specificity of 66.7% using a cutoff level of 3.0 for the tumor-to-background ratio. Low-grade sarcomas (grade I) and plasmocytomas revealed false-negative results, whereas various benign lesions, such as aneurysmatic bone cysts, fibromas, and parathyroid osteopathy, showed an enhanced 18F-FDG uptake. The data show that sensitivity and specificity depend on the variety of histologic subtypes included in a study and on the cutoff level defined by the investigators. Using our results with a cutoff of 1.0 SUV, the sensitivity was 81.58% and the specificity was 54.17%. In contrast, for a cutoff of 2.0 SUV, the sensitivity was 68.42% and specificity was 85.42%. These data show that the heuristic selection of a cutoff level is critical with respect to sensitivity and specificity. Therefore, statistical methods such as the DA should find preferential use for the data analysis.
A 2-tissue-compartment model is a generally accepted method for an accurate, detailed kinetic analysis of the 18F-FDG metabolism. To limit the burden for the patient, we chose the retrieval of the input function from the image data. Ohtake et al. (12) showed that the image-based data obtained from a vessel VOI consisting of at least 7 consequent ROIs correlate well with those obtained by arterial and venous blood sampling. We used VOIs instead of ROIs to maximize the information retrieved from the images. The high resolution, the small pixel size used for reconstruction, and the use of VOIs limit the partial-volume effects in our study (recovery coefficient, 0.85 for lesions of >8 mm in diameter). For small vessels with a diameter of <8 mm, a partial-volume correction was performed on the basis of phantom measurements of the recovery function.
However, no data have been provided about the impact of a quantitative procedure on diagnostic accuracy. The role of quantitative dynamic 18F-FDG studies with arterial blood sampling and calculation of metabolic rates was studied by Kole et al. (4). The authors reported the dynamic data, including arterial blood sampling for the input function and the lack of correlation between the metabolic rate of 18F-FDG and the aggressiveness of the neoplasms in 19 malignant and 7 benign bone lesions (4). Furthermore, the authors concluded that it was not possible to differentiate between malignant and benign bone tumors using the metabolic rate of glucose consumption. However, the authors confined the evaluation to the metabolic rate and did not include the kinetic parameters of a full 2-tissue-compartment model.
The use of several kinetic parameters obtained from the dynamic 18F-FDG data provides more information about 18F-FDG pharmacokinetics than the SUV of a single acquisition. The transport constant K1 is a parameter for the transport capacity of 18F-FDG, and the rate constant k3 is associated with the phosphorylation rate of the radiopharmaceutical. The blood volume in a tumor tissue is a parameter that modulates the uptake of the tracer. Therefore, the use of the VB of 18F-FDG is another parameter that can improve diagnostic accuracy. In addition to compartment analysis, the fractal dimension may help to quantify heterogeneity. In general, tumors showed a higher fractal dimensionas shown in the parametric image of Figure 1than benign processes. The increased fractal dimension is indicative of a more chaotic distribution of 18F-FDG. Furthermore, the example of Figure 1 reveals that the area of chaotic metabolism of 18F-FDG is larger than the area of enhanced 18F-FDG uptake expressed in SUV.
We used the t test for the basic analysis of the data and to determine the most statistically significant kinetic parameters for the differentiation of benign and malignant lesions. Interestingly, SUV, VB, Ki, and fractal dimension were significant parameters for the differentiation of these 2 groups. Some investigators have used the fractal dimension as a parameter for the assessment of spatial heterogeneity. Kleen et al. (14) used the fractal dimension as a scale-independent factor to measure spatial heterogeneity of blood flow. We used the fractal dimension to characterize the kinetics of 18F-FDG in all lesions. The mean value for the fractal dimension in malignant lesions is significantly higher than that in benign lesions. According to the data, the 18F-FDG turnover in benign, metabolically inactive lesions is more deterministic than that in tumors and in inflammations. The advantage of the fractal dimension is the lack of an input function and the good reproducibility of the values. DA, using only the fractal dimension data as input variable, revealed an improvement of sensitivity and accuracy in comparison to the analysis with SUV.
Analysis of the pharmacokinetic data revealed some interesting aspects about the bone lesions. To our knowledge, there are no reports about transport rates (K1k4) in bone lesions. To analyze the predictive value of 18F-FDG kinetics for the differential diagnosis (benign vs. malignant), we applied the DA to the data and compared the predicted classification with the histologically observed classification in each patient (Table 3). The results revealed a sensitivity of 75.86%, a specificity of 97.22%, and an accuracy of 87.69% when all transport constants, VB, Ki, fractal dimension, and SUV were used as input variables (nonlinear analysis). The data indicate that a negative PET study is reliable for the exclusion of a bone tumor. The coefficients of the nonlinear discriminant function are presented in the Appendix and can also be used for a prospective classification of a lesion into benign or malignant.
In addition to the impact of the evaluation parameters on the diagnostic accuracy of 18F-FDG PET, the prevalence of disease is another aspect that must be considered. The Bayesian analysis was applied to the data of this study to analyze sensitivity and specificity of the PET results with respect to different levels of the prevalence of disease. Theoretically, a diagnostic procedure should provide a high rate of true-positive results at a low level of prevalence of disease (prior probability of disease). Analysis of our data generally showed a high posterior probability of true-positive results at relatively low levels of prior probability for almost all parameters. Using only the 55- to 60-min SUV for evaluation, a posterior probability true-positive value of 0.8 was achieved at a 0.39 level of prior probability (Fig. 4A). However, the posterior probability of a false-negative value was not acceptable in comparison with the results obtained with the other kinetic parameters. In contrast, the fractal dimension was superior to the SUV concerning the posterior probability of false-negative results (Fig. 4A). However, the best results were obtained using the combination of all kinetic parameters with a posterior probability true-positive value of 0.8 at a prior probability of 0.235 (Fig. 4A). According to our data, the full kinetic analysis of an 18F-FDG study is the most accurate method for the differentiation of malignant and benign lesions at a low level of prior probability of disease. The Bayesian analysis showed a clear dependency of the PET results on the prevalence of disease (prior probability) and on the kinetic parameters used for the quantification of the PET images. This explains the highly variable specificity reported in the literature for 18F-FDG PET studies of bone lesions.
| CONCLUSION |
|---|
|
|
|---|
| APPENDIX |
|---|
|
|
|---|
![]() |
![]() |
Those variables with the largest standardized regression coefficients are the ones that contribute most to the prediction of group membership.
To discriminate malignant from benign lesions using all kinetic parameters, the following equations were applied:
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
The coefficients given in these 2 equations were calculated for the data included in this study. We classify the case as belonging to the group for which it has the highest classification score (www.statsoft.com/textbook/stdiscan.html).
Bayesian Analysis
The Bayesian analysis takes into account the different prevalence of disease and its influence on diagnostic accuracy. The Bayesian theorem was used to assess the performance of the different quantification parameters regarding different levels of disease prevalence.
The equation used to calculate the posterior probability of true-positive results (P1) for each quantification parameter was:
![]() |
The equation used to calculate the posterior probability of false-negative results (P2) for each quantification parameter was:
![]() |
The prevalence of disease P(D+) was varied from 0 to 1 in Equations P1 and P2: a=true-positive; b=false-negative; c=false-positive; d=true-negative.
The gain in information was calculated using the formula:
![]() |
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
For correspondence or reprints contact: Antonia Dimitrakopoulou-Strauss, MD, Department of Oncological Diagnostics and Therapy, E0105, Medical PET Group, Biological Imaging, German Cancer Research Center, Im Neuenheimer Feld 280, Heidelberg, D-69120 Germany.
E-mail: a.dimitrakopoulou-strauss{at}dkfz.de
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
S. Zhao, Y. Kuge, M. Kohanawa, T. Takahashi, Y. Zhao, M. Yi, K. Kanegae, K.-i. Seki, and N. Tamaki Usefulness of 11C-Methionine for Differentiating Tumors from Granulomas in Experimental Rat Models: A Comparison with 18F-FDG and 18F-FLT J. Nucl. Med., January 1, 2008; 49(1): 135 - 141. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. W. Hetts, S. D. Hilchey, R. Wilson, and B. Franc Case 110: Nonossifying Fibroma Radiology, April 1, 2007; 243(1): 288 - 292. [Full Text] [PDF] |
||||
![]() |
R. K. Heck Jr., T. D. Peabody, and M. A. Simon Staging of Primary Malignancies of Bone CA Cancer J Clin, November 1, 2006; 56(6): 366 - 375. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. Even-Sapir Imaging of Malignant Bone Involvement by Morphologic, Scintigraphic, and Hybrid Modalities J. Nucl. Med., August 1, 2005; 46(8): 1356 - 1367. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Zhao, Y. Kuge, T. Mochizuki, T. Takahashi, K. Nakada, M. Sato, T. Takei, and N. Tamaki Biologic Correlates of Intratumoral Heterogeneity in 18F-FDG Distribution with Regional Expression of Glucose Transporters and Hexokinase-II in Experimental Tumor J. Nucl. Med., April 1, 2005; 46(4): 675 - 682. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. G. Strauss, A. Dimitrakopoulou-Strauss, D. Koczan, L. Bernd, U. Haberkorn, V. Ewerbeck, and H.-J. Thiesen 18F-FDG Kinetics and Gene Expression in Giant Cell Tumors J. Nucl. Med., September 1, 2004; 45(9): 1528 - 1535. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. G. Strauss, A. Dimitrakopoulou-Strauss, and U. Haberkorn Shortened PET Data Acquisition Protocol for the Quantification of 18F-FDG Kinetics J. Nucl. Med., December 1, 2003; 44(12): 1933 - 1939. [Abstract] [Full Text] [PDF] |
||||
![]() |
S.-Y. Ma, L.-C. See, C.-H. Lai, H.-H. Chou, C.-S. Tsai, K.-K. Ng, S. Hsueh, W.-J. Lin, J.-T. Chen, and T.-C. Yen Delayed 18F-FDG PET for Detection of Paraaortic Lymph Node Metastases in Cervical Cancer Patients J. Nucl. Med., November 1, 2003; 44(11): 1775 - 1783. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| JOURNAL OF NUCLEAR MEDICINE TECHNOLOGY | THE JOURNAL OF NUCLEAR MEDICINE |