|
|
||||||||
Basic Science Investigation |
Mallinckrodt Institute of Radiology, Washington University School of Medicine, St. Louis, Missouri
Correspondence: For correspondence or reprints contact: Michael J. Welch, PhD, Mallinckrodt Institute of Radiology, Washington University School of Medicine, Campus Box 8225, 510 S. Kingshighway Blvd., St. Louis, MO 63110. E-mail: welchm{at}wustl.edu
| ABSTRACT |
|---|
|
|
|---|
Key Words: rodent microPET blood input function factor analysis
| INTRODUCTION |
|---|
|
|
|---|
Investigating changes in organ function, including myocardial function (e.g., blood flow, glucose metabolism, and free fatty acid metabolism), with the appropriate kinetic models requires accurate measurement of the blood input function and myocardial timeactivity curve. Even with current detector technology, which provides improved spatial resolution, cardiovascular research with small-animal PET still is limited by the inherent spatial resolution of the system (partial-volume effects and spillover), as well as the rapid movement of the beating heart. The size of the rodent heart is small compared with the resolution of the system, making it difficult to extract artifact-free blood-pool input functions and myocardial timeactivity curves from region of interests (ROIs) drawn on dynamic small-animal PET images.
A number of methods for measuring tracer activity concentrations in blood have been examined for application in kinetic models. Microblood sampling techniques are used routinely by our group (8). These techniques provide accurate blood input function measurements, but when multiple tracer studies are performed in the same animal, this technique is limited by the number of blood samples that can be withdrawn. Another method is the use of a ß-probe, which lies inside (9) or alongside (10,11) a blood vessel. The placement of an arterialvenous shunt also was attempted but was complicated by both the counts originating from the tissue (12) and the complex microsurgical procedure. A noninvasive approach with a theoretic basis (predefined shape for the blood input function) and computer simulation for extracting the blood input function from mice imaged with a small-animal PET scanner has been described (13). This method was validated with 18F-FDG kinetics in normal mice, and the theoretic basis of this method is applicable only to a limited number of animal models. The conventional noninvasive method currently used in a clinical PET setting requires the drawing of ROIs. However, when applied to dynamic small-animal PET images, this method is limited by large partial-volume effects and spillover from the rodent myocardium. This deficiency can be overcome by the use of factor analysis (FA) for the noninvasive extraction of the myocardial timeactivity curve and the blood input function.
In the 1980s, FA was introduced in nuclear medicine (1418). FA describes the covariance relationships among multiple variables in terms of a few underlying, but unobservable, random quantities called factors. All variables within a particular group are highly correlated. Variables in different groups have relatively small correlations. This method has been used to decompose dynamic images into component images and to determine their timeactivity curves by analysis of the variance in the data (14,15). It has been reported that the FA method successfully extracts the input functions and myocardial timeactivity curves from dynamic PET images of canines (16,17), humans (16), monkeys (18), and rats (19). In the rat study (19), the blood input functions were extracted with the FA method and were compared with those obtained from dynamic images by the ROI method. In the present study, the FA method was applied to dynamic small-animal PET images of rats and mice, and the timeactivity curves generated were compared with those obtained from traditional blood samples.
| MATERIALS AND METHODS |
|---|
|
|
|---|
5 µL at the center of the field of view on the basis of the 10-ns timing window and is less than
14 µL within a 1-cm radius. A complete performance evaluation of the scanner was done by Knoess et al. (20).
Animal Studies
To investigate the accuracy of the FA method in determining the blood input function and myocardial timeactivity curve, 2 groups of rodents were studied. The rodents had free access to standard rat chow and water. In group 1, 6 SpragueDawley rats (weight [mean ± SD], 305 ± 91 g) were studied with 18F-FDG (n = 3) and 1-11C-acetate (n = 3). In group 2, 6 BALB/c mice (weight, 20.8 ± 1.4 g) were studied with 18F-FDG (n = 3) and 1-11C-acetate (n = 3). The radiotracers (18F-FDG and 1-11C-acetate) were produced in our laboratory with a commercially available module (CTI Molecular Imaging).
Animal Preparations and Small-Animal PET Imaging
The use of animals was approved by the Animal Studies Committee at Washington University School of Medicine. Approximately 30 min before undergoing PET, rodents were anesthetized with inhalant anesthesia (1.5% isoflurane). Catheters were placed in the jugular vein for radiotracer administration and in the carotid artery for blood withdrawal during the scan. Rodents were placed in a supine position in the small-animal PET scanner.
All dynamic small-animal PET images were obtained with a microPET R4 housed in a temperature-controlled suite. After an intravenous bolus injection (100150 µL) of 18F-FDG (26.8 ± 4.5 MBq for rats and 24.1 ± 3.8 MBq for mice), list mode data were acquired for 60 min. The dynamic small-animal PET image data were sorted into 3-dimensional (3D) sinograms as follows: 24 frames of 5 s and 58 frames of 60 s. The frame durations were chosen to allow small intervals at early time points to show the rapid changes in the blood-pool input function and sampling at larger intervals at later time points to show the uptake characteristics of the myocardium. After an intravenous bolus injection of 1-11C-acetate (32.7 ± 4.8 MBq for rats and 20.5 ± 5.5 MBq for mice) and because of the rapid kinetics of the tracer in the heart, list mode data were acquired for 20 min. Dynamic sorting of the data into 3D sinograms was performed as follows: 24 frames of 5 s, 12 frames of 15 s, 10 frames of 30 s, and 10 frames of 60 s. The frame durations were chosen on the basis of the dynamic imaging protocols used for 1-11C-acetate in humans and large mammals studied at our institute.
The images were reconstructed by Fourier rebinning of the 3D sinograms followed by 2-dimensional filtered backprojection with a ramp filter (cutoff at Nyquist frequency). To obtain sufficient pixel density for analysis, the original images were enlarged with zoom factors of 3 and 4.225 for the rat and mouse studies, respectively, resulting in pixel sizes of 0.28 x 0.28 mm for the rat studies and 0.20 x 0.20 mm for the mouse studies. The slice thickness was 1.2 mm. At the time of the study, the only attenuation correction method available for this system was based on a 68Ge point source and led to excessive noise in the transmission scans. Therefore, attenuation and scatter corrections were not applied. In general, attenuation correction factors are constant over time and should not change the shape of the timeactivity curves. The timeactivity curves were adjusted for their dynamic range on the basis of the maximum and minimum of the dynamic images and rescaled at the end of the FA procedure, as will be described in the next section. These steps effectively accounted for the scale introduced by attenuation and partially compensated for the bias introduced by scatter. Although the system is now equipped with appropriate attenuation and scatter corrections and the images it provides are more quantitative, we do not believe the absence of attenuation and scatter corrections significantly changed the results of the present study.
To compare blood-pool input functions obtained with the FA method to blood activity concentrations measured directly, arterial blood samples (mean ± SD volumes of blood were 6.1 ± 0.4 µL for rats and 5.9 ± 0.3 µL for mice) were withdrawn from the carotid artery catheter into Wiretrol II precision micropipette tubes (Drummond Scientific Co.). The dead volume of the catheter and sampling line was 50100 µL, depending on the length of the catheter line and the size of the rodents used. During blood sampling, the catheter was not flushed, but a void sample was taken to ensure true blood sampling at each time point. For both types of animals and both radiotracers, blood samples were rapidly drawn for the scan duration (810). The micropipette tubes were preweighed and postweighed, and counting was done with a Beckman 3000
-counter for true blood activity concentration determinations.
FA Algorithm
FA assumes that the measured images are constructed from the superposition of a limited number of structures (or factors). For the present application, these structures correspond to the right ventricle, left ventricle, and myocardium. Each structure has a different timeactivity curve, which is dependent on the kinetics of the tracer used. In the present study, we used the methodology of Di Paola et al. (15), which was extended by several groups (1619). The FA method creates physiologic timeactivity curves and their corresponding images in 4 steps. The flow chart in Figure 1 shows the processing steps used for FA.
|
Second, after image processing, a principal components analysis was performed. The covariance matrix, eigenvectors, and eigenvalues (loadings) of the normalized time-dependent pixels were calculated for the 3 structures.
The third step was the determination of new factor axes that corresponded to the best representation without negative components for the 3 structures describing the dynamic data. These axes were calculated from the set of loading values for which the triangle subtended by the combination of any 3 loading values within the loading value space led to the maximum area. Thus, triangles were formed by any 3 loading values chosen out of 1,024 loading values for the mouse studies or 1,600 loading values for the rat studies. The normalized timeactivity curves and the corresponding factor images were generated with these new factor axes. The final timeactivity curves then were normalized by linear interpolation of the maximum and minimum values in the original images.
Finally, a constant given by the difference between the last data point for the FA blood input function and the last blood sample was added to left ventricle and right ventricle timeactivity curves.
Statistical Analysis
Correlations and percentage errors for input functions extracted from blood samples and the FA method were calculated with linear least squares regression analysis.
| RESULTS |
|---|
|
|
|---|
28%). Correlation functions with a slope of 1 are represented by dashed lines on Figure 5 for comparison.
|
|
|
|
| DISCUSSION |
|---|
|
|
|---|
For several decades, rodent blood sampling techniques have been used with a variety of tracers to quantify glucose metabolism and blood flow (2426) in rats and in mice (27,28). A comparison of glucose metabolic rates in rats and humans demonstrated a good correlation (28), showing that rapid arterial blood sampling in rodents produces an accurate input function. Huang et al. (13) described a theoretic approach to blood flow measurements, combined with limited blood sampling, in individual animals over a large population. This statistical analysisbased method discriminated inaccurate input functions from real input functions on the basis of a comparison of simulated blood sampling and image-derived input functions, each with its own associated errors. Bentourkia et al. (19) showed that the blood input function in rats could be extracted from the left ventricle with FA, but this method was not validated against blood sampling.
In the rat studies in this report, we could extract the 3 factor images (right ventricle, left ventricle, and myocardium) and generate timeactivity curves for each factor image using the FA method. The FA input function agreed well with that obtained from blood sampling in the 18F-FDG and 1-11C-acetate rat studies. This finding shows that although PET images of rats are subject to partial-volume effects and cardiac motion blurring, the FA method overcomes these potential limitations and is able to produce a true blood input function from either right ventricle or left ventricle factors.
The FA input function (extracted from the right ventricle factor) agreed well with that obtained from direct blood sampling in the 18F-FDG and 1-11C-acetate mouse studies. However, the left ventricle timeactivity curve had a shape similar to the summation of the right ventricle blood input function and myocardial timeactivity curve, indicating that the left ventricle timeactivity curve was contaminated by the myocardial timeactivity curve. As a result, the factor image of the left ventricle showed contributions from the factor images of both the blood pool and the myocardium. Because of the very small size of the mouse heart compared with the resolution of the system, it was extremely difficult to distinguish suborgan structures and to extract their timeactivity curves from dynamic PET images with any image-based methods, including FA. The application of FA in mice therefore showed limitations when an imaging system with insufficient resolution was used. It is possible that the next generation of scanners will allow for more accurate FA in mice.
As mentioned in the Results section, the slope of the linear regression lines was consistently greater than unity. In general, the image-derived timeactivity curves yielded higher peak maxima than the well counter blood sample data, suggesting an underestimation of the peaks even when rapid blood sampling techniques were used and leading to the slopes observed in Figure 5. Similar observations were made by Wu et al. (18), who compared the input function derived from an ROI in the left ventricle of small primates with the input function derived from blood sampling and well counting.
As anticipated, the shapes of FA myocardial timeactivity curves obtained from mouse and rat studies were similar. Residual spillover from the blood pool still existed and could be seen as a peak at early time points on the FA timeactivity curves. There is no current standard for extracting myocardial timeactivity curves from rodent data for comparison with myocardial timeactivity curves generated from FA. In both radiotracer studies, myocardial timeactivity curves were not corrected for residual partial-volume effects and spillover, and the FA blood-pool input function was not corrected for metabolite activity. Therefore, the measured ratio of myocardium to blood-pool activity concentrations appeared to be smaller than the expected value. Nevertheless, for extracting the blood input function from dynamic images, the FA method has an obvious advantage over the ROI method. FA allows a clear separation of the different compartments for an imaging system characterized by a given spatial resolution. FA and image-based extraction of blood input functions and myocardial timeactivity curves of rodents may be improved by the use of cardiac gated dynamic images (29) or a higher-resolution small-animal PET scanner (21).
Such advances may further improve the performance of the FA method and lead to the replacement of multiple blood sampling in rodent studies.
| CONCLUSION |
|---|
|
|
|---|
This technique has been used to quantify myocardial blood flow in rats (30), in which the measurement has been validated by comparison with the gold standard radioactive microsphere method. In addition, it has been used to quantify myocardial metabolism in the Zucker rat model of diabetes (31). In these 2 instances, this technique has been used to quantify important cardiac metabolic parameters.
| ACKNOWLEDGMENTS |
|---|
| References |
|---|
|
|
|---|
Related articles in JNM:
This article has been cited by other articles:
![]() |
S. J. Kim, J. S. Lee, K. C. Im, S.-Y. Kim, S.-A. Park, S. J. Lee, S. J. Oh, D. S. Lee, and D. H. Moon Kinetic Modeling of 3'-Deoxy-3'-18F-Fluorothymidine for Quantitative Cell Proliferation Imaging in Subcutaneous Tumor Models in Mice J. Nucl. Med., December 1, 2008; 49(12): 2057 - 2066. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. L. Franc, P. D. Acton, C. Mari, and B. H. Hasegawa Small-Animal SPECT and SPECT/CT: Important Tools for Preclinical Investigation J. Nucl. Med., October 1, 2008; 49(10): 1651 - 1663. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. I. Shoghi, R. J. Gropler, T. Sharp, P. Herrero, N. Fettig, Y. Su, M. S. Mitra, A. Kovacs, B. N. Finck, and M. J. Welch Time Course of Alterations in Myocardial Glucose Utilization in the Zucker Diabetic Fatty Rat with Correlation to Gene Expression of Glucose Transporters: A Small-Animal PET Investigation J. Nucl. Med., August 1, 2008; 49(8): 1320 - 1327. [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] |
||||
![]() |
G. Z. Ferl, X. Zhang, H.-M. Wu, and S.-C. Huang Estimation of the 18F-FDG Input Function in Mice by Use of Dynamic Small-Animal PET and Minimal Blood Sample Data J. Nucl. Med., December 1, 2007; 48(12): 2037 - 2045. [Abstract] [Full Text] [PDF] |
||||
![]() |
H.-M. Wu, G. Sui, C.-C. Lee, M. L. Prins, W. Ladno, H.-D. Lin, A. S. Yu, M. E. Phelps, and S.-C. Huang In Vivo Quantitation of Glucose Metabolism in Mice Using Small-Animal PET and a Microfluidic Device J. Nucl. Med., May 1, 2007; 48(5): 837 - 845. [Abstract] [Full Text] [PDF] |
||||
![]() |
W. K. Schiffer, M. M. Mirrione, and S. L. Dewey Optimizing Experimental Protocols for Quantitative Behavioral Imaging with 18F-FDG in Rodents J. Nucl. Med., February 1, 2007; 48(2): 277 - 287. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. J. Welch, J. S. Lewis, J. Kim, T. L. Sharp, C. S. Dence, R. J. Gropler, and P. Herrero Assessment of Myocardial Metabolism in Diabetic Rats Using Small-Animal PET: A Feasibility Study J. Nucl. Med., April 1, 2006; 47(4): 689 - 697. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. Herrero, J. Kim, T. L. Sharp, J. A. Engelbach, J. S. Lewis, R. J. Gropler, and M. J. Welch Assessment of Myocardial Blood Flow Using 15O-Water and 1-11C-Acetate in Rats with Small-Animal PET J. Nucl. Med., March 1, 2006; 47(3): 477 - 485. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| JOURNAL OF NUCLEAR MEDICINE TECHNOLOGY | THE JOURNAL OF NUCLEAR MEDICINE |