Abstract
The investigation of cerebral metabolic rate of glucose (CMRGlu) at baseline and during specific tasks previously required separate scans with the drawback of high intrasubject variability. We aimed to validate a novel approach to assessing baseline glucose metabolism and task-specific changes in a single measurement with a constant infusion of 18F-FDG. Methods: Fifteen healthy subjects underwent two PET measurements with arterial blood sampling. As a reference, baseline CMRGlu was quantified from a 60-min scan after 18F-FDG bolus application using the Patlak plot (eyes closed). For the other scan, a constant radioligand infusion was applied for 95 min, during which the subjects opened their eyes at 10–20 min and 60–70 min and tapped their right thumb to their fingers at 35–45 min and 85–95 min. The constant-infusion scan was quantified in two steps. First, the general linear model was used to fit regional time–activity curves with regressors for baseline metabolism, task-specific changes for the eyes-open and finger-tapping conditions, and movement parameters. Second, the Patlak plot was used for quantification of CMRGlu. Multiplication of the baseline regressor by β-values from the general linear model yielded regionally specific time–activity curves for baseline metabolism. Further, task-specific changes in metabolism are directly proportional to changes in the slope of the time–activity curve and hence to changes in CMRGlu. Results: Baseline CMRGlu from the constant-infusion scan matched that from the bolus application (test–retest variability, 1.1% ± 24.7%), which was not the case for a previously suggested approach (variability, −39.9% ± 25.2%, P < 0.001). Task-specific CMRGlu increased in the primary visual and motor cortices for eyes open and finger tapping, respectively (P < 0.05, familywise error–corrected), with absolute changes of up to 2.1 μmol/100 g/min and 6.3% relative to baseline. For eyes open, a decreased CMRGlu was observed in default-mode regions (P < 0.05, familywise error–corrected). CMRGlu quantified with venous blood samples (n = 6) showed excellent agreement with results obtained from arterial samples (r > 0.99). Conclusion: Baseline glucose metabolism and task-specific changes can be quantified in a single measurement with constant infusion of 18F-FDG and venous blood sampling. The high sensitivity and regional specificity of the approach offer novel possibilities for functional and multimodal brain imaging.
The investigation of glucose metabolism in the human brain represents an important tool for the assessment of brain function and alterations thereof. These include tumor imaging (1), neurologic diseases such as epilepsy (2) and Alzheimer disease (3–5), and psychiatric disorders such as major depression (6). Hence, PET using the radiolabeled glucose analog 18F-FDG offers an excellent opportunity for imaging and quantification of glucose metabolism. 18F-FDG is considered an irreversible radioligand since it is trapped in the early metabolic pathway and dephosphorylation is energetically less favorable (k4 ≠ 0, described in Materials and Methods section and Discussion section) (7,8). After a bolus application, the uptake of the radioligand is rather stable after 30 min, enabling simplified methods such as SUVs and ratios to be used in clinical routine for the assessment of baseline metabolism.
Similarly, if a certain task is performed in the early uptake period, task-specific glucose metabolism can be quantified (8). Several applications have been proven successful in this context. These include changes in the visual cortex after visual stimulation (9–11), changes in amygdala metabolism when a subject is watching videos that evoke emotion (12), and changes during more complex paradigms such as passive problem solving (13) and car driving (14). However, the differentiation of task-dependent changes from baseline metabolism in these studies required two separate PET scans, introducing the issue of considerable intrasubject variability between repeated scans, commonly reaching 25% if the scans are several days apart (15,16). Therefore, the acquisition of baseline and task-specific metabolism in a single scan is highly desirable to increase both sensitivity and specificity.
One solution to the problem of repeated scans is use of a multiple-injection protocol to identify different physiologic states. In such a protocol, radioligand applications are temporally separated by 30–70 min in one continuous measurement (17–19). This, however, implies that imaging of each state and potential task performance consumes the same time, making it difficult to study task-related effects. Furthermore, the increase in scanning time and accompanied movement artifacts may limit the number of reasonably achievable states to two. In view of the sophisticated modeling requirements together with the increased radiation burden for subjects, the approach did not become a popular solution.
Recently, a novel technique with constant infusion of 18F-FDG was proposed (20). In this technique, the radioligand is constantly applied during the entire scan, enabling assessment of dynamic task-specific changes in glucose metabolism. The same concept was proposed earlier, although in a different context, using simulations. It was shown that a certain task may change radioligand uptake during the experiment, but constant infusion was not considered at that time (21). The advantage of constant radioligand infusion is, however, its circumvention of the issue that 18F-FDG uptake is most prominent early after a bolus injection and shows little change thereafter. Hence, a continuous supply of 18F-FDG provides free radioligand available to bind according to the metabolism at hand throughout the experiment.
Compared with well-established protocols using 15O-H2O or functional MRI, the assessment of neuronal activation with 18F-FDG may offer various benefits. The first two approaches quantify blood flow and the level of blood oxygenation, respectively, and are thus dependent on hemodynamics and neurovascular coupling. On the other hand, 18F-FDG is a surrogate for glucose metabolism. The different imaging modalities are therefore measuring different aspects of the energy budget, namely oxygen supply versus glucose utilization, with the assets and drawbacks discussed elsewhere (8). This complementary nature of the different approaches may be particularly interesting in multimodal imaging studies using hybrid PET/MR systems (22). Furthermore, 15O-H2O studies are technically more challenging and functional MRI requires MR compatibility; hence, 18F-FDG may represent a useful alternative.
Although the previous study provided the proof of principle, only 3 subjects were evaluated, and comparison with the commonly applied bolus application of 18F-FDG was missing (20). Therefore, we aimed to validate the technique of constant radioligand infusion for the quantification of glucose metabolism at baseline and during two different tasks in a single PET measurement.
MATERIALS AND METHODS
Subjects
Fifteen right-handed healthy subjects participated in this study (mean age ± SD, 24.5 ± 4.7 y; 8 women and 7 men). To rule out severe physical or mental disorders, all subjects were assessed by an experienced psychiatrist using the Structural Clinical Interview for DSM-IV. Additionally, they underwent routine medical examinations at the initial screening visit, which included blood laboratory testing, electrocardiography, and assessment of general physical and neurologic status (23). Female participants underwent a urine pregnancy test before each PET scan. The subjects reported no past or current substance abuse. All provided written informed consent after receiving a detailed explanation of the study protocol. The study was approved by the Ethics Committee of the Medical University of Vienna (ethics number 1916/2013), and the procedures were performed according to the Declaration of Helsinki.
Study Design
All subjects underwent two PET measurements in randomized order (6.1 ± 21.9 d apart). First, a 60-min measurement was performed using a bolus application in which the subjects wore an eye mask to keep their eyes closed and were instructed to stay awake. This scan served as a reference to validate the constant-infusion approach. The second measurement lasted for 95 min and was done with a constant radioligand infusion for the entire measurement (Volumed μVP7000; Arcomed). The infusion speed was set to 36 mL/h (20). During this scan, the subjects opened their eyes at 10–20 min and 60–70 min after injection start and repeatedly tapped their right thumb to their fingers at 35–45 min and 85–95 min after injection start, with the eyes closed and no finger movement at other times. None of the subjects reported falling asleep during the bolus or infusion scan.
PET and Blood Sampling
18F-FDG was freshly prepared every day according to a well-established method (24) using the E FASTlab platform (GE Healthcare). PET measurements were performed using an Advance (GE Healthcare) PET scanner as described previously (23). Briefly, head movement was minimized by placing each subject’s head in a cushioned polyurethane bowl with straps around the forehead. An attenuation correction scan was performed for 5 min with retractable 68Ge rod sources before radioligand injection. The dynamic PET scans were acquired in 3-dimensional mode with an intravenous (bolus or constant infusion) 3 MBq/kg injection of 18F-FDG. For the constant infusion, no initial bolus was applied and the infusion rate corresponded to 1.895 MBq/kg/h. PET images were reconstructed to consecutive frames of 1 min. The spatial resolution of reconstructed images was 4.36 mm full width half maximum at 1 cm from the center of the field of view (voxel size, 3.125 × 3.125 × 4.25 mm).
Blood Sampling
The glucose concentration in blood was measured before each PET scan. For the arterial input function, arterial blood samples were taken from the radial artery (25) during the PET scans. Automatic sampling was done for the first 3 min at a rate of 4 mL/min (Allogg). Manual samples were acquired at 5, 10, 20, 40, and 60 min for the bolus scan or at 10, 20, 35, 45, 60, 70, 85, and 95 min for the constant-infusion scan. For manual samples, whole-blood activity was measured in a γ-counter (Wizard2, 3″; Perkin Elmer), and plasma was separated by centrifugation and again measured. The arterial input function was constructed by combining whole-blood activities from automated and manual sampling, which were further multiplied by average plasma–to–whole-blood ratios. In addition, for 6 constant-infusion scans, manual venous blood samples were obtained from the cubital vein simultaneously with the arterial samples and were processed in the same way.
MRI
A T1-weighted structural MRI scan was acquired for all subjects (mMR; Siemens) and used for spatial normalization (magnetization-prepared rapid gradient-echo sequence; echo time/repetition time, 4.2/2,000 ms; voxel size, 1 × 1 × 1.1 mm).
Data Preprocessing and Regions of Interest
Data were processed in SPM12 with default parameters unless specified otherwise. PET images were corrected for head motion (quality, 1), and the mean image was coregistered to the structural MR images (normalized mutual information). The structural MR images were then spatially normalized to MNI (Montreal Neurologic Institute) space via segmentation (26), and the obtained transformation matrices were applied to the dynamic PET images. Spatially normalized PET images were smoothed with an 8 × 8 × 8 mm gaussian kernel. Regions of interest (ROI) were taken from probabilistic atlases provided in FSL (http://fsl.fmrib.ox.ac.uk/) and included frontal, temporal, parietal, occipital, cingulate, and somatosensory cortices; thalamus; caudate; putamen; amygdala-hippocampus complex; and cerebellar gray matter. Regional time–activity curves were obtained as averages across each region (Fig. 1).
Quantification with Bolus Application
The bolus measurement was quantified in PMOD, version 3.703 (http://pmod.com/), for all ROIs. Here, the influx constant,Eq. 1was estimated with the Gjedde–Patlak plot (27) asEq. 2with CT and CP being concentrations in tissue and plasma, respectively (Fig. 1). The plot of CT(T)/CP(T) versus integral/CP(T) becomes linear after time t*, that is, when the free compartment has reached its steady state. Hence, Ki can then be estimated as the slope. The arterial input function was fitted with a sum of 3 exponentials from the peak onward (25). The cerebral metabolic rate of glucose was further calculated asEq. 3with GluP and LC being the glucose concentration in plasma and the lumped constant, respectively. Here, the lumped constant was set to 0.89 (28,29), and although some studies may have used different values, this factor affects only the absolute values in cerebral metabolic rate of glucose (CMRGlu) but not Ki or task-specific changes in glucose metabolism. The actual lumped constant is thus of little relevance as long as it is consistently used within a study. For comparison, Ki and CMRGlu were also estimated using an unconstrained 2-tissue-compartment model with k4 ≠ 0 (3).
Quantification with Constant Infusion
Before data analysis, PET images from constant-infusion measurements were masked with SPM12’s gray matter prior converted into a binary code at 0.1 to exclude non–gray matter voxels from the analysis. Since 18F-FDG uptake occurs rather slowly with constant infusion and the paradigms lasted for 10 min each, sudden changes in time–activity curves were unlikely to be of physiologic origin. To reduce noise in PET signals (residual scatter, residual randoms, movement), a 12th-order finite impulse-response low-pass filter was applied with a cutoff frequency of 5 min (i.e., half the duration of the paradigm), which can be seen as a temporal smoothing.
Task-specific and baseline glucose metabolism from the constant-infusion measurement was quantified voxelwise in a 2-step procedure (Fig. 2) using Matlab R2011a (The MathWorks). The Discussion section presents the underlying assumptions. First, a general linear model (GLM) was used to fit each voxel’s time–activity curve, which was adopted from functional MRI analyses (30).Eq. 4The GLM design matrix included 4 predictors, yielding voxelwise β-estimates for each regressor. One baseline regressor (yielding βbase) was defined individually as follows. Using linear regression, the average time course across all gray matter voxels was modeled by a third-order polynomial, whereas task effects were modeled as nuisance variables. Hence, the baseline regressor estimated by the polynomial was free of task effects. One regressor was included for each task (βeye and βfinger), defined as a linear ramp function with slope of 1 during the task and slope of 0 otherwise (20). The task regressors were made orthogonal to the baseline regressor. Hence, the common variance is accounted for with the baseline regressor, whereas each task regressor models only the variance specific for each task. The fourth regressor was included to remove movement-related artifacts (βmove). A principal component analysis was computed on the 6 motion parameters (3 translation and 3 rotation), and the first principal component score was used as a regressor (supplemental data; supplemental materials are available at http://jnm.snmjournals.org). In the second step, the Patlak plot was used to estimate Ki. Here, βbase estimated from the GLM analysis was multiplied by the baseline regressor, which yielded regionally specific tissue time–activity curves for each voxel’s baseline metabolism (CT_base). As Ki equals the slope in Equation 2, baseline Ki_base is thus directly proportional to CT_baseEq. 5Further, task-specific changes in metabolism are proportional to changes in the slope of the time–activity curve (20), hence changing βbase by the value of βtask:Eq. 6from which it follows that task-specific metabolism is proportional toEq. 7The arterial input function for the Patlak plot was fitted with a sum of two exponentials asEq. 8with a and τ representing estimated parameters (Figs. 1 and 2). t* was fixed to 30 min for simplicity.
To evaluate the potential substitution of arterial blood samples by venous ones, baseline and task-specific Ki also was computed with venous blood samples using the same procedure.
For comparison, we also estimated baseline Ki and CMRGlu from the previously suggested approach (20). Basically, this included only the first part of the above-described method, namely the linear regression but without the Patlak plot analysis. Hence, Ki was previously defined as the slope of CT(T)/CP(T) versus time.
Statistics
The bias of baseline Ki between the Patlak plot and the 2-tissue-compartment model was calculated by regression analysis, that is, evaluating slope and intercept across all ROIs. For validation of the constant-infusion approach, the estimated baseline Ki and CMRGlu were compared with those from the bolus application. Test–retest values were calculated for each ROI asEq. 9Test–retest values were then compared with the original approach (20) by repeated-measures ANOVA with region and method as fixed factors and subject as random factor. Significant changes in task-related glucose metabolism were computed voxelwise in SPM12 using a 1-sample t test for each of the tasks (P < 0.05, familywise error [FWE]–corrected). To evaluate the substitution of arterial blood samples by venous ones, Ki and CMRGlu estimated with both blood samples were compared using linear regression implemented in Matlab across all ROIs plus 3 regions with task-relevant changes in metabolism.
RESULTS
For the bolus application, baseline Ki values were generally lower when estimated with the Patlak plot than with the 2-tissue-compartment model (slope, 0.73; intercept, 0.0004). For the constant-infusion approach, baseline metabolism was adequately fitted with the GLM (Fig. 3) and was consistent with values from the bolus application (test–retest Ki, 2.4% ± 16.0%; CMRGlu, 1.1% ± 24.7%; supplemental data). In contrast, the previous method (20) underestimated baseline metabolism (test–retest Ki, −39.0% ± 17.0%; CMRGlu, −39.9% ± 25.2%), which, across all ROIs, was significantly different from the approach presented here (F1,333 = 608.4 and 241.4, both P < 0.001, for Ki and CMRGlu, respectively).
Task-specific increases of up to 2.1 μmol/100 g/min in glucose metabolism were observed in the primary visual and primary motor cortices for the eyes-open condition and right-finger tapping, respectively, as well as in the cerebellum for both tasks (all P < 0.05, FWE-corrected voxel level) (Figs. 3 and 4). In these regions, percentage signal changes reached 6.1% ± 1.7% and 6.3% ± 2.8%, for the eyes-open condition and right-finger tapping, respectively (Table 1). Conversely, no relevant changes were identified in the motor cortex during the eyes-open condition or in the visual cortex during finger tapping, even without correction for multiple comparisons (−0.4%, Table 1). Interestingly, we also observed task-specific decreases in metabolism for the eyes-open condition, mostly in the frontal regions and the default-mode network (i.e., posterior cingulate and medial prefrontal cortices, Fig. 4; all P < 0.05, FWE-corrected voxel or cluster level).
Although differences between arterial and venous samples were observed (31,32), changes in ∫CP(t)dt (i.e., the factor determining Ki in the steady state) were highly correlated between both blood pools (r = 0.99; slope, 0.92; intercept, −0.19). Hence, quantification with venous blood samples showed excellent agreement with the results from arterial samples. This was true for baseline Ki and CMRGlu across all ROIs, as well as for increases in the primary visual cortex when the eyes were open, increases in the primary motor cortex during finger tapping, and decreases in the angular gyrus when the eyes were open (all r > 0.99; slope, 0.99…1.04; intercept, 0…1.31).
DISCUSSION
In this work, we have shown that baseline and task-specific changes in glucose metabolism can be robustly estimated in a single measurement with constant infusion of 18F-FDG. Task-related changes were regionally and temporally specific for the eyes-open condition and right-finger tapping, enabling separation of different metabolic demands.
Compared with the previously introduced approach (20), several extensions and validations were included here. The combination of the GLM with the subsequent Patlak plot analysis markedly decreased test–retest variability for baseline metabolism. The exact difference from the previous approach is, however, difficult to estimate because the previous study did not explicitly specify the baseline term. Our comparison with the bolus-application scans validates the use of a constant infusion to robustly estimate baseline metabolism. Proceeding from this validation, task-specific changes could be quantified because these are proportional to changes in the slope of the time–activity curve. Interestingly, the estimated percentage change in signal (∼6%) was lower than previously reported, possibly because of the simplicity of the tasks used here. For instance, light stimulation has been shown to produce 6%–12% increased metabolism in the visual cortex (9) and more than 20% during checkerboard stimulation (9,20). Again, the fact that our stimuli were rather subtle underlines the sensitivity of the approach, whereas a more pronounced effect can be expected from paradigms with stronger stimuli. On the other hand, the changes in glucose metabolism were higher than those usually observed (∼1%–2% (30)) in functional MRI experiments using the blood-oxygen level–dependent effect, even for 7 T imaging (33). Given the complementary nature of different imaging techniques (8,11,22,34,35), the presented approach may offer novel applications for simultaneous multimodal imaging with hybrid imaging systems. Finally, the excellent agreement that was found when venous and arterial blood samples were compared enables further simplification and thus greater accessibility of the method.
The approach also showed high regional specificity, as no metabolic changes in the visual cortex were observed during finger-tapping or in the motor cortex when eyes were opened. The increased glucose metabolism in the visual cortex is in line with various previous studies (9,11,35). In addition, increased metabolism in the cerebellum during the eyes-open condition and finger tapping may relate to the cerebellum’s involvement in attention (36), visual processing (37), and motor control (38). Furthermore, we observed decreases in glucose metabolism in regions of the frontal cortex and of the default-mode network. Such a deactivation during task performance was first described with PET blood flow and oxygen extraction studies (34,39) and in a marked number of functional MRI experiments thereafter (40). Interestingly, a recent study reported increased glucose metabolism when the eyes-open and eyes-closed conditions were compared (11). This disagreement with our finding is probably related to a difference in study design, namely recruitment of two different groups to evaluate two different states (11). Further considering an intrasubject variability of around 25% for repeated scans, we argue that the assessment of baseline and task-specific metabolism in a single measurement is essential for functional PET studies investigating the dynamic activation of distinct brain regions.
The method proposed here incorporates several assumptions and limitations. First, application of the Patlak plot does not take into account the fact that k4 may actually be higher than 0, resulting in a known underestimation of CMRGlu (here, ∼27%) (3,29). Although this is important to keep in mind when different studies are compared, the underestimation is a scaling factor and may thus be less problematic if applied consistently across subjects. Further, the arterial input function did not reach equilibrium after bolus application or during constant infusion. For the Patlak plot, this is not required, as changes in arterial input function need to be only small enough to allow the reversible part to reach a steady state (27,41). Here, the changes in arterial input function were 54% ± 40% lower, and linearity was reached more quickly for constant infusion than for bolus, supporting the applicability of the Patlak plot. Another assumption is that baseline metabolism can be modeled by a single baseline regressor, which is multiplied by regional β-estimates. Although the baseline regressor was cleaned of task effects and reasonable fits were obtained for all investigated ROIs, the assumption needs further validation in populations with regional alterations such as brain tumors. Next, we assume that the relationship between the neuronal response and glucose metabolism is linear time–invariant, similar to functional MRI (30). In this context, linearity implies that changes in the neuronal response will convert to changes in metabolism by the same factor—a reasonable assumption given their tight coupling (8). Conversely, changes in metabolism by two simultaneous events will sum together and add to the baseline metabolism. Time invariance refers to the temporal coupling between the stimulus and the 18F-FDG uptake. Although these assumptions are derived from functional MRI, PET imaging is less dependent on the prerequisite of neurovascular coupling, and 18F-FDG uptake was previously unrelated to pronounced changes in blood flow (20). Hence, linear time invariance seems to be a valid simplification for the proposed method. Finally, the approach assumes that both baseline and task-specific metabolism are at a steady state and that task-related glucose consumption is constant throughout the task. As discussed previously (20), this steady-state assumption is violated by almost all functional imaging methods but is a necessary approximation to enable quantitative description. On the other hand, if a task does not elicit a constant change in metabolism, regressors with shapes different from a linear ramp function may yield a better estimate of task-specific glucose consumption.
As this study is among the first of its kind, the approach could potentially benefit from further advances. For instance, the combination of a bolus application plus a constant infusion may yield a better equilibrium of plasma activity and increase the signal-to-noise ratio. Another improvement would be not only to start but also to end the paradigm with a baseline condition to increase confidence in the task-specific model fits. Finally, a more explicit stimulation (in this case, for example, explicit pacing of finger tapping (42)) may yield more robust changes in CMRGlu. It remains to be established how well the approach performs during more complex paradigms, especially when different conditions are being compared (e.g., activations and deactivations during recognition of emotions or objects (33)). The detection of task-relevant decreases in CMRGlu gives confidence that the approach is sensitive enough for such contrasts.
CONCLUSION
We have demonstrated that glucose metabolism at baseline and during several tasks can be quantified simultaneously using a constant 18F-FDG infusion and venous blood sampling. Estimating these different aspects of glucose metabolism in a single measurement decreases the intrasubject variability inherent in repeated scans. In view of its high sensitivity to subtle stimuli together with its regional specificity, the method offers novel applications in functional brain imaging.
DISCLOSURE
The costs of publication of this article were defrayed in part by the payment of page charges. Therefore, and solely to indicate this fact, this article is hereby marked “advertisement” in accordance with 18 USC section 1734. Gregor Gryglewski is the recipient of a DOC Fellowship from the Austrian Academy of Sciences at the Institute of Psychiatry and Psychotherapy, Medical University of Vienna. This research was supported by a grant from the Austrian Science Fund (FWF KLI 516) to Rupert Lanzenberger. This scientific project was performed with the support of the Medical Imaging Cluster of the Medical University of Vienna. No other potential conflict of interest relevant to this article was reported.
Acknowledgments
We are grateful to the technical team of the PET Centre, especially Harald Ibeschitz, Ingrid Leitinger, and Volker Weiss. We thank the physicians—Anna Höflich, Arkadiusz Komorowski, Marie Spies, Pia Baldinger, Christoph Kraus, and Dietmar Winkler—for medical support.
Footnotes
Published online Jul. 7, 2016.
- © 2016 by the Society of Nuclear Medicine and Molecular Imaging, Inc.
REFERENCES
- Received for publication March 23, 2016.
- Accepted for publication June 7, 2016.