|
|
||||||||
BASIC SCIENCE INVESTIGATIONS |
PET-Center and Department of Medicine V, Aarhus University Hospital, and Institute for Experimental Clinical Research, Aarhus University, Aarhus, Denmark; and Department of Mathematics, University of Queensland, Brisbane, Australia
| ABSTRACT |
|---|
|
|
|---|
Key Words: liver glucose metabolism FDG kinetics compartment model PET liver blood flow
| INTRODUCTION |
|---|
|
|
|---|
This study was designed to compare the dual-input function with the arterial-input function for kinetic parameter estimation using 3-O-[11C]methylglucose (MG) and FDG. The errors made with the arterial-input function were examined, and the importance of using the dual-input function was shown. The hypothesis tested was that the use of the dual-input function, instead of the conventional arterial-input function, would improve model fitting and provide physiologic parameter estimates.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Animal Preparation
Each pig was anesthetized by intramuscular injection of 10 mg/kg (milligram per kilogram of body weight) ketamine and 0.5 mg/kg midazolam followed by intravenous injections of 5 mg/kg and 0.25 mg/kg, respectively. The anesthesia was maintained using intravenous infusions of 48 mg/min/kg ketamine and 0.40.8 mg/min/kg midazolam, and the pig was ventilated with air containing 40% oxygen using a Servo 900 respirator (Siemens-Elema, Solna, Sweden). Each pig was placed on its back, and a catheter was placed in the caval vein through the right femoral vein for infusions. A catheter was placed into the aorta through the right femoral artery for arterial blood sampling corresponding to the hepatic artery (HA) and for blood pressure measurements. A catheter was placed through a 10-cm-long abdominal incision into the portal vein (PV) through the splenic vein for PV blood sampling. The position of the PV catheter was controlled by fluoroscopy before each experiment. Finally, two sonographic transit-time flow-meter probes (CardioMed; In Vivo Aps, Oslo, Norway) were placed around the HA and PV through the abdominal incision. The pig was then allowed to recover for 1 h before administration of the first tracer. The oxygen saturation and pH in the arterial blood samples were measured every hour and were adjusted toward >98% and 7.45, respectively, by changing the amount of air delivered from the respirator. The arterial blood glucose level was measured three times during each scan and was in the range of 5.06.7 mmol/L. The pigs were covered and placed on a thermostatically controlled heat blanket set to keep the body temperature between 38.5°C and 39.5°C.
PET Examination and Blood Sampling
All six pigs were scanned after 500 MBq [15O]carbon monoxide (CO) inhalation, 500 MBq MG 15-s intravenous injection, and 300 MBq FDG 15-s intravenous injection. Between scans, we waited six times the radioactive half-lives to allow the previous tracer to decay. PET was performed using an ECAT EXACT HR-47 camera (CTI, Knoxville, TN/Siemens Medical Systems, Inc., Hoffman Estates, IL). The PET camera was calibrated using a phantom containing a 68Ge/68Ga solution with a known activity concentration. The pig was positioned with the liver inside the 15-cm field of view. A 15-min transmission scan with external rod sources was performed before the first emission scan and used for photon attenuation correction. After the CO inhalation, eight 20-s frames were recorded with corresponding arterial blood samples. For the two glucose scans, dynamic data acquisition was started at the injection time and included 38 frames: 18 x 10 s, 4 x 30 s, 5 x 1 min, 6 x 5 min, and 5 x 10 min, for a total scanning time of 90 min. The PET data were reconstructed with filtered backprojection using a Hanning filter with a cutoff at 0.3 Nyquist frequency. The resulting three-dimensional images contained 128 x 128 x 47 voxels with a size of 2.4 x 2.4 x 3.1 mm and a central spatial resolution of 6.7 mm full width at half maximum. In the pig, MG and FDG do not enter red blood cells (S. Keiding, unpublished data), in agreement with previous findings for glucose and galactose (5). The activity concentrations were measured in the blood, and the clearance of the tracer from the blood to tissue therefore refers to the blood flow. HA and PV blood samples of 1 mL each were collected manually for 88 min starting at the injection time in the following intervals: 18 x 5 s, 3 x 10 s, 2 x 30 s, 1 x 2 min, 1 x 3 min, 1 x 4 min, 1 x 6 min, and 7 x 10 min. Radioactivity concentrations were measured using a well counter (Packard Instrument Co., Meriden, CT), which was calibrated using a 68Ge/68Ga solution with a known activity concentration. All measurements in the blood samples were decay corrected to the start of the tracer administration. Blood flows in HA and PV were measured continuously throughout the experiment. After finishing the measurements, the liver was removed from the anesthetized animal and weighted.
Image Analysis
Regions of interest were defined using transaxial slices of the mean image consisting of the sum of all frames 1290 min after injection. This procedure made it easier to distinguish the liver from the surrounding tissue in comparison with using only a single frame from the study. The regions of interest were drawn to contain liver tissue without including extrahepatic tissue and summed to a volume of interest of about 200 mL. Liver tissue timeactivity curves (TACs) were generated from these volumes of interest.
Modeling
CO.
CO is transported by the red blood cells, which cannot cross the endothelial cells. During the steady state, the hepatic blood volume, Vb, can be calculated using the liver tissue activity concentration, M, and the blood input activity concentration, Cin:
![]() | (Eq. 1) |
MG.
The bloodtissue exchange of the glucose analogs takes place in the hepatic sinusoids. Blood is separated from the hepatocytes by the endothelial cells permeable for plasma and an extracellular volume including the space of Disse (6). MG is a nonmetabolizable glucose analog that enters and leaves the hepatocytes by the hepatic GLUT-2 transport protein but is not phosphorylated (7). The dynamic data obtained after the MG injection were analyzed using a one-tissue compartment model accounting for MG in an extended blood volume and in an intracellular space. The extended blood volume parameter accounts for MG in the extracellular volume, including the hepatic blood volume and space of Disse. The model has three free parameters: the clearance into the cell, K1; the reverse rate constant, k2; and the extended blood volume, V0. Denoting tissue activity concentration by M(t) and blood input activity concentration by Cin(t), this model predicts:
![]() | (Eq. 2) |
denotes a convolution integral. The external PET measurement of M(t) is direct, whereas the input to the liver, Cin(t), is inferred from blood samples, being either the dual-input or the arterial-input function.
FDG.
FDG enters the hepatocytes by the hepatic GLUT-2 transport protein. In the hepatocytes, FDG is phosphorylated through hexokinase to FDG-6PO4 (7), which can be dephosphorylated by glucose-6-phosphatase (8). Thus, FDG can be transported out of the cells again. Any further metabolism through the glycolytic series is assumed not to occur within the experimental period. The FDG data were analyzed using a two-tissue compartmental model (9,10), accounting for FDG in the extended blood volume and in the intracellular space, and a pool representing FDG-6PO4. In addition to the parameters in the MG model (K1, k2, and V0), the FDG model includes rate constants for the phosphorylation of FDG, k3, and the dephosphorylation of FDG-6PO4, k4. The two-tissue compartmental model then predicts:
![]() |
![]() |
![]() |
![]() | (Eq. 4) |
GjeddePatlak Representation.
In the liver MG and FDG experiments discussed in this article, two short-time processes occur. After a 15-s injection, the dynamic phase of the input lasts about 1 min (Fig. 1). The equilibration between the extracellular and intracellular pools of the blood-borne precursor has the relaxation time 1/k2, also about 1 min. For a nonmetabolized tracer, such as MG, or an irreversibly metabolized tracer, such as FDG, in most extrahepatic tissues, no other time scales are defined. Reversible metabolism introduces another, usually much longer, time scale of equilibration between the precursor and the metabolite pool. A typical example is the characteristic time given by the rate constant for dephosphorylation of FDG in the liver, 1/k4
100 min, which is longer than the duration of our experiment.
|
![]() | (Eq. 5) |
(t) is a nonlinear, monotonically increasing function of time. The ordinate M(t)/Cin(t) changes rapidly on the short time scales and then slowly, or not at all, on the long time scale. A particularly clear-cut application of the GjeddePatlak plot is to a precursor exchanged between the blood and cells with parameters K1 and k2 and then metabolized irreversibly in the cells with a rate constant of k3. In that case, beyond the short time scales, the plot has the form:
![]() | (Eq. 6) |
![]() | (Eq. 7) |
Parameter Estimation and Statistical Criteria.
Model parameters were estimated by minimizing the weighted residual sum of squares (WRSS), and PET data were weighted in proportion to the frame durations. Plots of the weighted residuals against time were examined for systematic errors. Identifiability of the model parameters was examined using sensitivity functions:
![]() | (Eq. 8) |
The best model fitted to the data is not necessarily the model producing the smallest WRSS, because adding more parameters generally decreases the WRSS. Therefore, we identified the statistically favorable model based on two criteria that included penalty functions proportional to the number of parameters in the model (13). The Akaike (14) and Schwarz (15) criteria (Eqs. 9 and 10, respectively) were used to measure the quality of a fit:
![]() | (Eq. 9) |
![]() | (Eq. 10) |
Results are given as the mean ± SEM. Pairwise comparisons were performed using t tests; P < 0.05 was considered statistically significant.
| RESULTS |
|---|
|
|
|---|
TimeActivity Curves
The dual-input concentrations were calculated as the flow-weighted input concentrations from HA and PV according to:
![]() | (Eq. 11) |
Typical PET TACs for tissue activity concentrations after MG and FDG injections are shown in Figure 2. Higher activity concentrations at late time points show that FDG is retained longer in the liver than is MG. This result reflects the transport of FDG into and out of the intracellular FDG-6PO4 pool.
|
|
|
1 min, so data from the first few minutes were not included in the analysis. This early phase also includes the dynamic phase, for which the inputs are very different (Fig. 1). Parameter estimates and the statistical criteria are presented in Table 2. Examples of GjeddePatlak plots are shown in Figures 4 and 5. The quantities K and V pertaining to the asymptote are seen to be unaffected by the errors arising from the use of the arterial input. This result is a consequence of the differences between arterial and dual inputs being confined to the short time scales.
|
|
|
| DISCUSSION |
|---|
|
|
|---|
The conventional arterial input leads to absurdly low estimates for the extracellular volume in the liver, which we refuted by separate measurements with CO, and leads to systematic underestimation of rate constants for rapid bloodtissue exchange. In general, kinetic parameters determined using compartmental models and arterial input are underestimated, as shown in Table 1. The FDG rate constants obtained using arterial input correspond to values obtained from FDG liver studies in fasting humans, for which corrected TACs from the left ventricular cavity of the heart were used as the input function (4). In the same study, on the basis of AIC and SC, it was found that a two-tissue compartment model without vascular volume, but with a time-delayed input function, was the best kinetic model to describe the hepatic FDG PET data. In view of our results, these findings can be explained by their use of a single input function. Using our data and the dual-input function, we found AIC of 89 ± 8 and SC of 97 ± 8 including vascular volume, and AIC of 134 ± 9 and SC of 140 ± 12 excluding vascular volume. Thus, the statistical criteria clearly favor the more physiologic model including a vascular volume. Because our blood sampling sites were close to the liver, no time delay was needed.
In our experiments, when integrating the arterial input and the dual input over the entire sampling period, we found that the difference between the integrals was always <1%. This result is important, because any loss of tracer in the intestines would make the use of arterial concentrations clearly wrong for both compartmental and GjeddePatlak analysis. Our GjeddePatlak analysis showed a zero slope when used to analyze MG data. This result was expected, because MG has a completely reversible metabolism. In the more complicated cases, such as FDG metabolism with reversible phosphorylation, the steady state may not be reached within the time scale of these experiments, and the GjeddePatlak plot may not have an asymptotic straight line. Its interpretation therefore involves problems that do not concern us here, except inasmuch as the division of the plot into short and long time scales is preserved, and the parameter estimates are unaffected by the use of arterial samples. For MG and FDG, the GjeddePatlak analysis will produce the same parameter estimates for both arterial input and dual input, because the analysis is confined to time scales for which the reversible pools have reached a steady state, and by that time, the arterial-input and dual-input functions are similar. The GjeddePatlak analysis, which is often used for evaluating liver tumors (1,3), can therefore be used for the quantification of the forward metabolic clearance using only arterial samples.
| CONCLUSION |
|---|
|
|
|---|
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
For correspondence or reprints contact: Ole L. Munk, MSc, PET-Center, Aarhus University Hospital, Nørrebrogade 44, DK 8000 Aarhus, Denmark.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
M. Sorensen, O. L. Munk, F. V. Mortensen, A. K. Olsen, D. Bender, L. Bass, and S. Keiding Hepatic uptake and metabolism of galactose can be quantified in vivo by 2-[18F]fluoro-2-deoxygalactose positron emission tomography Am J Physiol Gastrointest Liver Physiol, July 1, 2008; 295(1): G27 - G36. [Abstract] [Full Text] [PDF] |
||||
![]() |
K Murase, S Miyazaki, and X Yang An efficient method for calculating kinetic parameters in a dual-input single-compartment model Br. J. Radiol., May 1, 2007; 80(953): 371 - 375. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. Iozzo, A. Gastaldelli, M. J. Jarvisalo, J. Kiss, R. Borra, E. Buzzigoli, A. Viljanen, G. Naum, T. Viljanen, V. Oikonen, et al. 18F-FDG Assessment of Glucose Disposal and Production Rates During Fasting and Insulin Stimulation: A Validation Study J. Nucl. Med., June 1, 2006; 47(6): 1016 - 1022. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Bertoldo, J. Price, C. Mathis, S. Mason, D. Holt, C. Kelley, C. Cobelli, and D. E. Kelley Quantitative Assessment of Glucose Transport in Human Skeletal Muscle: Dynamic Positron Emission Tomography Imaging of [O-Methyl-11C]3-O-Methyl-D-Glucose J. Clin. Endocrinol. Metab., March 1, 2005; 90(3): 1752 - 1759. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. Iozzo, A. K. Turpeinen, T. Takala, V. Oikonen, J. Bergman, T. Gronroos, E. Ferrannini, P. Nuutila, and J. Knuuti Defective Liver Disposal of Free Fatty Acids in Patients with Impaired Glucose Tolerance J. Clin. Endocrinol. Metab., July 1, 2004; 89(7): 3496 - 3502. [Abstract] [Full Text] [PDF] |
||||
![]() |
O. L. Munk, L. Bass, H. Feng, and S. Keiding Determination of Regional Flow by Use of Intravascular PET Tracers: Microvascular Theory and Experimental Validation for Pig Livers J. Nucl. Med., November 1, 2003; 44(11): 1862 - 1870. [Abstract] [Full Text] [PDF] |
||||
![]() |
O. L. Munk, S. Keiding, and L. Bass Impulse-response function of splanchnic circulation with model-independent constraints: theory and experimental validation Am J Physiol Gastrointest Liver Physiol, October 1, 2003; 285(4): G671 - G680. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. Iozzo, K. Hallsten, V. Oikonen, K. A. Virtanen, R. Parkkola, J. Kemppainen, O. Solin, F. Lonnqvist, E. Ferrannini, J. Knuuti, et al. Effects of Metformin and Rosiglitazone Monotherapy on Insulin-Mediated Hepatic Glucose Uptake and Their Relation to Visceral Fat in Type 2 Diabetes Diabetes Care, July 1, 2003; 26(7): 2069 - 2074. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. Iozzo, K. Hallsten, V. Oikonen, K. A. Virtanen, J. Kemppainen, O. Solin, E. Ferrannini, J. Knuuti, and P. Nuutila Insulin-Mediated Hepatic Glucose Uptake Is Impaired in Type 2 Diabetes: Evidence for a Relationship with Glycemic Control J. Clin. Endocrinol. Metab., May 1, 2003; 88(5): 2055 - 2060. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. Iozzo, F. Geisler, V. Oikonen, M. Maki, T. Takala, O. Solin, E. Ferrannini, J. Knuuti, and P. Nuutila Insulin Stimulates Liver Glucose Uptake in Humans: An 18F-FDG PET Study J. Nucl. Med., May 1, 2003; 44(5): 682 - 689. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Kaarstad, D. Bender, L. Bentzen, O. L. Munk, and S. Keiding Metabolic Fate of 18F-FDG in Mice Bearing Either SCCVII Squamous Cell Carcinoma or C3H Mammary Carcinoma J. Nucl. Med., July 1, 2002; 43(7): 940 - 947. [Abstract] [Full Text] [PDF] |
||||
![]() |
O. L. Munk, S. Keiding, G. Brix, M. E. Bellemann, S. I. Ziegler, and U. Haberkorn Quantification of 18F-FDG Uptake in the Liver Using Dynamic PET J. Nucl. Med., March 1, 2002; 43(3): 439 - 441. [Full Text] [PDF] |
||||
![]() |
D. Bender, O. L. Munk, H.-Q. Feng, and S. Keiding Metabolites of 18F-FDG and 3-O-11C-Methylglucose in Pig Liver J. Nucl. Med., November 1, 2001; 42(11): 1673 - 1678. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| JOURNAL OF NUCLEAR MEDICINE TECHNOLOGY | THE JOURNAL OF NUCLEAR MEDICINE |