|
|
|||||||||
Clinical Investigations |
Department of Cardiology and PET Center, Institute for Cardiovascular Research-VU, VU Medical Center, Amsterdam, The Netherlands
| ABSTRACT |
|---|
|
|
|---|
Key Words: PET input function FDG myocardial metabolic rate of glucose
| INTRODUCTION |
|---|
|
|
|---|
Quantitative analysis requires an input function, for which online arterial blood sampling (AS) serves as the gold standard. Whenever possible, however, AS should be avoided because of rare potential complications associated with such an invasive procedure (1012). In addition, frequent blood sampling contributes to the radiation burden of personnel. Alternative methods used for characterization of the input function are sampling of arterialized venous blood (13) and blood timeactivity curves derived from PET images themselves, using regions of interest (ROIs) on vascular structures such as the left ventricle (LV) (14,15), left atrium (LA) (16), or aorta (1719). Image-derived input functions (IDIF) are attractive because of their noninvasive nature and because the scanning protocol can be simplified. Several reports have supported the use of IDIF in clinical practice, provided that appropriate partial-volume and spillover corrections are applied (14,17). It should be noted that, since those reports, spatial resolution of PET scanners has improved substantially.
To date, no study has determined which vascular structure is best suited for defining ROIs in cardiac FDG PET studies. Recently, a study comparing the aortic arch, LV, and LA indicated that, in the majority of cases, the aortic arch was the best vascular structure for defining IDIF for oncologic applications (18). Although this finding was supported by independent quality control measures, no direct comparison with an externally measured arterial input function was provided. Moreover, patients were fasted for more than 6 h, significantly reducing potential spillover effects compared with cardiac studies.
This study was conducted to evaluate the accuracy of IDIF from the LV, LA, ascending aorta (AA), and descending aorta (DA) for the calculation of myocardial MRGlu (during a hyperinsulinemic, euglycemic clamp). As determined by online continuously withdrawn arterial blood, the arterial input function was used as the gold standard.
| Materials and Methods |
|---|
|
|
|---|
|
Preparation of FDG
FDG was synthesized by the BV Cyclotron VU (Amsterdam, The Netherlands) according to a modified method previously described by Hamacher et al. (21) using nuclear interface equipment. Radiochemical purity was >97%.
Scanning Protocol
All scans were performed in 2-dimensional mode, using an ECAT EXACT HR+ scanner (Siemens/CTI, Knoxville, TN). This scanner acquires 63 planes of data over an axial field of view of 15 cm.
Using rotating 68Ge line sources, a 10-min transmission scan was performed after a short transmission scan for patient positioning. Subsequently, 370 MBq FDG dissolved in 5 mL saline were injected intravenously (followed by a 47-mL saline flush at a rate of 2 mL/s). The dynamic scan consisted of 39 frames with variable frame length (6 x 5, 6 x 10, 3 x 20, 5 x 30, 5 x 60, 8 x 150, 6 x 300 s) for a total time of 60 min. Arterial 18F concentrations were monitored continuously using a fully automated blood-sampling device (Veenstra Instruments, Joure, The Netherlands) (22). Withdrawal rate was 5 mL/min during the first 10 min and 2.5 mL/min thereafter. At set times (5, 10, 20, 30, and 50 min after injection), continuous withdrawal was interrupted for collection of manual samples, after which the arterial line was flushed with heparinized saline. These samples were measured in a well counter, cross-calibrated against the PET scanner, and used for online calibration of the blood sampler and for determination of plasmatowhole-blood 18F ratios. In addition, these samples were used for measurement of plasma glucose levels (hexokinase method, Hitachi 747; Boeringer Mannheim, Mannheim, Germany). All dynamic scan data were corrected for physical decay of 18F and for dead time, scatter, randoms, and photon attenuation. The images were reconstructed using filtered backprojection with a Hanning filter at the Nyquist frequency. This resulted in a transaxial spatial resolution of approximately 7-mm full width at half maximum.
Image Analysis
ROIs were defined on several planes over the LV, LA, AA, and DA on early transaxial images (2550 s after injection of FDG).
To avoid spillover from the left atrium, only planes in which the LA was not visible were used for the AA and DA. The distal DA was not used because its decreasing size might have resulted in considerable partial-volume effects. ROIs over the LV were copied to the last frame to check for possible inclusion of the myocardial wall. For uniformity, all blood-pool ROIs had a circular shape and were placed at the center of each structure (although the LV is not circular in a transaxial slice, the use of noncircular ROIs did provide similar results). Timeactivity curves were generated by projecting the ROIs onto the complete dynamic dataset. All timeactivity curves of the individual blood-pool ROIs were compared with each other in each study. No attempt was made to select only input functions that were close to the arterial samples, because the purpose of this study was to evaluate the use of IDIF in lieu of AS. Therefore, only obvious outlying factors, such as a disproportionalrise of the tail of the curve, were rejected. The individual ROIs of each blood pool were grouped, creating a volume of interest. Then, to avoid bias by statistical noise, each volume of interest was automatically copied twice and the copies were placed directly adjacent to the original. To avoid a large shift from the center of the structure, ROIs were partly overlapping and the positioning of each volume of interest was verified. In this way, 3 individual volumes of interest were created, all identical with respect to plane selection. Their timeactivity curves were averaged and the average curve was considered to represent the optimal IDIF of that particular blood pool.
To determine the effect of ROI size on partial-volume effects and statistical noise, this procedure was performed for 3 different ROI sizes (diameters, 5, 10, and 15 mm). First, volumes of interest were created with 15-mm ROIs. ROI size was then reduced to create volumes of interest with smaller ROIs, but at the very same position and identical with respect to plane selection. Thus, 12 IDIF were generated for each study (3 different ROI sizes for 4 vascular structures).
The last 3 frames were summed (4560 min after injection) and resliced, resulting in a single-frame short-axis image, which was used to define the myocardial tissue ROIs. These ROIs were defined over several planes and were grouped to create a single, whole-myocardium ROI. This ROI was projected onto the complete dynamic resliced dataset, generating a whole-myocardium timeactivity curve.
Input functions derived from the same blood pool might differ significantly, depending on how ROIs were defined. As a consequence, a variation in MRGlu might be observed. To evaluate possible interobserver variability in MRGlu, a second investigator defined ROIs over the LV, LA, AA, and DA independently, the only requirement being the use of ROIs 15 mm in diameter. Both investigators used the same tissue ROIs.
Data Analysis
To correct for delay, the peaks of all image-derived vascular timeactivity curves in each patient were averaged and the difference in arrival time of the peak of this averaged curve and that of the sampler curve was used to shift the sampler curve.
Both the sampler and the image-derived vascular timeactivity curves were transformed to plasma timeactivity curves using the arterial plasmatowhole-blood 18F ratios.
To determine MRGlu for the whole-myocardium ROI, both compartmental and Patlak graphical analyses (23) were used. For nonlinear regression, the standard 2-tissue compartment model (13), with and without setting k4 to zero, was used. In both fits, a spillover blood volume parameter was incorporated to account for intravascular activity. For model selection, both the Akaike information criterion (24) and the Schwarz criterion (25) were used. Based on these selection criteria, the best model was used for comparison of compartmental and Patlak analyses using the online blood sampler data as input function. The Patlak analysis was performed over the period from 10 to 60 min after injection.
The ratio of MRGlu obtained with the IDIF and MRGlu obtained with AS was calculated for all patients with the Patlak graphical analysis (23). Limits of agreement were assessed by means of the analysis described by Bland and Altman (26). This analysis was also used to assess the agreement of the kinetic analyses for determining MRGlu and to evaluate the interobserver variability introduced with IDIF.
| Results |
|---|
|
|
|---|
TimeActivity Curves
The highest peak activity was observed in the LA, followed by the LV, AA, and DA, respectively. All image-derived timeactivity curves showed an earlier and higher peak than that of the online blood-sampling device (Fig. 1). LA and DA timeactivity curves from small ROIs showed a slightly higher peak than timeactivity curves from larger ROIs (6% and 9%, respectively), whereas no significant differences were observed in LV and AA timeactivity curves.
|
|
|
In 1 patient, the parameter estimation from nonlinear regression analysis had to be rejected because of clear abnormal results. With Patlak analysis (23), no deficient fits were observed. Good correlation was observed between the results of compartmental and Patlak analysis (y = 1.01x; R2 = 0.99). Spread of MRGlu values (without the single outlier) are shown in Figure 3 (mean difference = -0.01, 2 SD = 0.01 µmol/mL/min). The range of SEs for the individual MRGlu as obtained from the compartmental and Patlak analyses was 0.001.52 µmol/mL/min and 0.000.02 µmol/mL/min, respectively. In 8 of the 18 studies, the SE was much higher with compartmental analysis; in the remaining 10 studies, it was comparable with that obtained with Patlak analysis.
|
|
The individual ratios of all patients, using the large ROI (diameter, 15 mm), are shown in Figure 4. For LV and DA, the result of 1 study was regarded as an obvious outlier (ratios of 1.18 and 0.66, respectively) and therefore not used for further calculations. The highest variation was observed when the DA was used for characterization of the input function.
|
|
|
| Discussion |
|---|
|
|
|---|
In all image-derived vascular timeactivity curves, earlier and higher peak activity of FDG was observed than in the arterial plasma timeactivity curve, probably because of delay and dispersion of tracer in the arm and blood-sampling device. After shifting the arterial curve for delay, however, the difference of the integrated values for the first 5 min after injection of FDG was, on average, only a few percentage points. When Patlak graphical analysis (23) is used for calculation of MRGlu, the impact of such a difference in peak activity will be small.
In contrast, overestimation of FDG activity at later time points, seen in timeactivity curves derived from the LV and LA, will have a major effect on the calculation of MRGlu in both Patlak (23) and compartmental analysis. These effects of spillover from myocardial tissue into the LV blood pool had already been recognized in early FDG PET studies (14,27,34). In an attempt to minimize partial- volume and spillover effects, small ROIs are frequently used. The present results indicate, however, that a smaller ROI will only modestly reduce spillover effects, whereas statistical noise in the timeactivity curves will increase. Timeactivity curves from the LV were especially sensitive to statistical noise, probably because of the limited number of planes that could be used for defining ROIs. Timeactivity curves derived from the AA suffered less from statistical noise and corresponded most consistently with the arterial plasma timeactivity curve (Fig. 2).
In this study, results of Patlak graphical analysis (23) showed good agreement with compartmental analysis, in agreement with previous reports (14,18,28). Consequently, because Patlak graphical analysis is less sensitive to noise, it was used for the comparison of different input functions in determining MRGlu.
For each patient in this study, ratios of MRGlu obtained with IDIF and MRGlu obtained with the arterial plasma timeactivity curve as input function were calculated. The mean ratios for AA and DA were close to 1, indicating that these results were comparable with the results obtained with arterial sampling (Table 3). In contrast, LA and LV resulted in significant underestimation of myocardial MRGlu. A similar underestimation of MRGlu was observed when short-axis ROIs were used instead of transaxial LV ROIs. MRGlu values for small ROIs were only slightly higher than those obtained with larger ROIs. Apparently, there was only a minor reduction in spillover effects. This is partly because of the smaller overlap between adjacent circles for the small ROIs compared with the large ROIs; the difference between the effective ROIs used in the analysis was somewhat smaller then suggested by the difference in radii. In addition, partial-volume effects for LV ROIs will be affected by cardiac and respiratory movement. As a result of higher noise levels, a larger variation was observed in the results for the smaller ROIs (Table 3). Therefore, large ROIs (diameter, approximately 15 mm) appeared to be more useful for characterization of the input function. The IDIF that had the highest agreement with the online arterial blood sampler appeared to be the AA. Although the DA was a good alternative, there was a larger variability in the results (Table 3; Fig. 4).
Another important advantage of AA and DA was the lower interobserver variation of MRGlu than with LA and LV, probably because of the smaller size of the aortic structures, forcing the investigators to define almost identical ROIs. Interobserver variability could be improved using short-axis ROIs instead of transaxial LV ROIs but still was significantly higher than with AA.
This study indicates that without correction for partial- volume and spillover effects, the AA is the best structure for defining IDIF. Several methods for partial-volume and spillover correction have been developed (14,15,17,2729,35) and their use might have led to better results for LV and LA. Assessment of the various correction methods was, however, beyond the scope of this study. Although, in theory, partial-volume and spillover correction methods will allow the use of larger ROIs, it should be noted that any correction will also lead to some degree of statistical degradation of signals (i.e., subtraction of counts). Therefore, the purpose of this study was to determine which vascular structure would require the smallest correction.
| Conclusion |
|---|
|
|
|---|
| FOOTNOTES |
|---|
For correspondence or reprints contact: Arno P. van der Weerdt, MD, Department of Cardiology, Room 6N120, VU Medical Center, Postbus 7057, 1007 MB Amsterdam, The Netherlands.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
D. Vriens, L.-F. de Geus-Oei, W. J.G. Oyen, and E. P. Visser A Curve-Fitting Approach to Estimate the Arterial Plasma Input Function for the Assessment of Glucose Metabolic Rate and Response to Treatment J. Nucl. Med., December 1, 2009; 50(12): 1933 - 1939. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. Zanotti-Fregonara, R. Maroy, C. Comtat, S. Jan, V. Gaura, A. Bar-Hen, M.-J. Ribeiro, and R. Trebossen Comparison of 3 Methods of Automated Internal Carotid Segmentation in Human Brain PET Studies: Application to the Estimation of Arterial Input Function J. Nucl. Med., March 1, 2009; 50(3): 461 - 467. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. J. de Langen, M. Lubberink, R. Boellaard, M. D. Spreeuwenberg, E. F. Smit, O. S. Hoekstra, and A. A. Lammertsma Reproducibility of Tumor Perfusion Measurements Using 15O-Labeled Water and PET J. Nucl. Med., November 1, 2008; 49(11): 1763 - 1768. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. Tomasi, P. Edison, A. Bertoldo, F. Roncaroli, P. Singh, A. Gerhard, C. Cobelli, D. J. Brooks, and F. E. Turkheimer Novel Reference Region Model Reveals Increased Microglial and Reduced Vascular Binding of 11C-(R)-PK11195 in Patients with Alzheimer's Disease J. Nucl. Med., August 1, 2008; 49(8): 1249 - 1256. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y.-H. D. Fang and R. F. Muzic Jr. Spillover and Partial-Volume Correction for Image-Derived Input Functions for Small-Animal 18F-FDG PET Studies J. Nucl. Med., April 1, 2008; 49(4): 606 - 614. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Schroeder, M. F. Vidal Melo, G. Musch, R. S. Harris, J. G. Venegas, and T. Winkler Image-Derived Input Function for Assessment of 18F-FDG Uptake by the Inflamed Lung J. Nucl. Med., November 1, 2007; 48(11): 1889 - 1896. [Abstract] [Full Text] [PDF] |
||||
![]() |
B.-K. Teo, Y. Seo, S. L. Bacharach, J. A. Carrasquillo, S. K. Libutti, H. Shukla, B. H. Hasegawa, R. A. Hawkins, and B. L. Franc Partial-Volume Correction in PET: Validation of an Iterative Postreconstruction Method with Phantom and Patient Data J. Nucl. Med., May 1, 2007; 48(5): 802 - 810. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. M. Sondergaard, M. Bottcher, M. Marie Madsen, O. Schmitz, S. B. Hansen, T. T. Nielsen, and H. E. Botker Impact of Type 2 Diabetes on Myocardial Insulin Sensitivity to Glucose Uptake and Perfusion in Patients with Coronary Artery Disease J. Clin. Endocrinol. Metab., December 1, 2006; 91(12): 4854 - 4861. [Abstract] [Full Text] [PDF] |
||||
![]() |
L.-F. de Geus-Oei, E. P. Visser, P. F.M. Krabbe, B. A. van Hoorn, E. B. Koenders, A. T. Willemsen, J. Pruim, F. H.M. Corstens, and W. J.G. Oyen Comparison of Image-Derived and Arterial Input Functions for Estimating the Rate of Glucose Metabolism in Therapy-Monitoring 18F-FDG PET Studies J. Nucl. Med., June 1, 2006; 47(6): 945 - 949. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. J. Hoekstra, S. G. Stroobants, E. F. Smit, J. Vansteenkiste, H. van Tinteren, P. E. Postmus, R. P. Golding, B. Biesma, F. J.H.M. Schramel, N. van Zandwijk, et al. Prognostic Relevance of Response Evaluation Using [18F]-2-Fluoro-2-Deoxy-D-Glucose Positron Emission Tomography in Patients With Locally Advanced Non-Small-Cell Lung Cancer J. Clin. Oncol., November 20, 2005; 23(33): 8362 - 8370. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Lubberink, R. Boellaard, A. P. van der Weerdt, F. C. Visser, and A. A. Lammertsma Quantitative Comparison of Analytic and Iterative Reconstruction Methods in 2- and 3-Dimensional Dynamic Cardiac 18F-FDG PET J. Nucl. Med., December 1, 2004; 45(12): 2008 - 2015. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. Kostakoglu and S. J. Goldsmith 18F-FDG PET Evaluation of the Response to Therapy for Lymphoma and for Breast, Lung, and Colorectal Carcinoma J. Nucl. Med., February 1, 2003; 44(2): 224 - 239. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. J. Hoekstra, O. S. Hoekstra, S. G. Stroobants, J. Vansteenkiste, J. Nuyts, E. F. Smit, M. Boers, J. W.R. Twisk, and A. A. Lammertsma Methods to Monitor Response to Chemotherapy in Non-Small Cell Lung Cancer with 18F-FDG PET J. Nucl. Med., October 1, 2002; 43(10): 1304 - 1309. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | RSS | TABLE OF CONTENTS |
| JOURNAL OF NUCLEAR MEDICINE TECHNOLOGY | THE JOURNAL OF NUCLEAR MEDICINE |