Abstract
Quantitative total-body PET imaging of blood flow can be performed with freely diffusible flow radiotracers such as 15O-water and 11C-butanol, but their short half-lives necessitate close access to a cyclotron. Past efforts to measure blood flow with the widely available radiotracer 18F-fluorodeoxyglucose (FDG) were limited to tissues with high 18F-FDG extraction fraction. In this study, we developed an early-dynamic 18F-FDG PET method with high temporal resolution kinetic modeling to assess total-body blood flow based on deriving the vascular transit time of 18F-FDG and conducted a pilot comparison study against a 11C-butanol reference.
Methods The first two minutes of dynamic PET scans were reconstructed at high temporal resolution (60×1 s, 30×2 s) to resolve the rapid passage of the radiotracer through blood vessels. In contrast to existing methods that use blood-to-tissue transport rate (K1) as a surrogate of blood flow, our method directly estimates blood flow using a distributed kinetic model (adiabatic approximation to the tissue homogeneity model; AATH). To validate our 18F-FDG measurements of blood flow against a flow radiotracer, we analyzed total-body dynamic PET images of six human participants scanned with both 18F-FDG and 11C-butanol. An additional thirty-four total-body dynamic 18F-FDG PET scans of healthy participants were analyzed for comparison against literature blood flow ranges. Regional blood flow was estimated across the body and total-body parametric imaging of blood flow was conducted for visual assessment. AATH and standard compartment model fitting was compared by the Akaike Information Criterion at different temporal resolutions.
Results 18F-FDG blood flow was in quantitative agreement with flow measured from 11C-butanol across same-subject regional measurements (Pearson R=0.955, p<0.001; linear regression y=0.973x–0.012), which was visually corroborated by total-body blood flow parametric imaging. Our method resolved a wide range of blood flow values across the body in broad agreement with literature ranges (e.g., healthy cohort average: 0.51±0.12 ml/min/cm3 in the cerebral cortex and 2.03±0.64 ml/min/cm3 in the lungs, respectively). High temporal resolution (1 to 2 s) was critical to enabling AATH modeling over standard compartment modeling.
Conclusions Total-body blood flow imaging was feasible using early-dynamic 18F-FDG PET with high-temporal resolution kinetic modeling. Combined with standard 18F-FDG PET methods, this method may enable efficient single-tracer flow-metabolism imaging, with numerous research and clinical applications in oncology, cardiovascular disease, pain medicine, and neuroscience.
Introduction
Imaging blood flow has garnered considerable interest over the past 50 years as its dysfunction is characteristic in many diseases.1–3 PET imaging with a blood flow-specific radiotracer such as 11C-butanol or 15O-water is widely considered the gold standard for blood flow imaging.4–6 These flow radiotracers are freely diffusible across capillary membranes4–6 and accordingly the measured PET signal is closely proportional to blood flow. Blood flow can then be quantified by a simple one-tissue compartment model due to the complete or near-complete extraction of these freely diffusible flow radiotracers.4–6 Importantly, these flow radiotracers are highly extracted in tissue across the entire body and allow total-body imaging of blood flow.4,5,7
However, the short half-lives of the radioisotopes in flow radiotracers create practical challenges that hinder their broader accessibility. 15O-water has a half-life of 2.04 minutes, which necessitates an onsite cyclotron and a dose delivery system.8 11C-butanol comparably has a longer half-life of 20.40 minutes, but still requires nearby production, thus limiting access to urban or research PET centers. Other flow radiotracers, such as 82RbCl and 13N-ammonia, similarly have short radioisotope half-lives and high costs in addition to having non-linear uptake with flow.9 A blood flow imaging method using a widely available radiotracer such as 18F-fluorodeoxyglucose (FDG) may mitigate these challenges and open attractive opportunities for imaging of blood flow and glucose metabolism with a single-tracer dynamic scan.
Early dynamic 18F-FDG PET has been used to measure blood flow in selected tissue such as tumors,10 liver,11 and myocardium12 where 18F-FDG is moderately to highly extracted. The first 2-to-3-minute dynamic 18F-FDG PET signal is principally weighted towards the initial tissue delivery of the radiotracer13 and the higher regional extraction fraction makes the analysis amenable to simplified modeling like that of freely diffusible flow radiotracers. However, these approaches are not generally applicable to other regions like the brain, which has lower 18F-FDG extraction fraction.14,15
An intravenously injected tracer is delivered to local tissue vasculature at a rate equal to blood flow. Standard compartmental models neglect this transient process, but distributed kinetic models explicitly model the blood flow and transit time associated with the radiotracer traversing the blood vessels.16,17 Although described several decades ago, distributed models had limited application in PET due to the poor temporal resolution and statistical quality of time-activity curves measured with conventional PET scanners.18,19
Total-body PET has substantially greater sensitivity20–22 over conventional PET systems and allows high-temporal resolution dynamic imaging21,23 and kinetic modeling.13,24,25 This may revitalize opportunities to apply distributed kinetic models for blood flow estimation with 18F-FDG in various tissues. In this study, we developed an early dynamic 18F-FDG PET method for total-body blood flow imaging with high temporal resolution kinetic modeling and conducted a pilot validation against a 11C-butanol PET reference standard.
Methods
Total-body Dynamic PET
Two human cohorts were pooled in this study, each separately approved by the Institutional Review Board at the University of California, Davis. Written informed consent was obtained for all participants. All participants received total-body dynamic imaging on the uEXPLORER PET/CT system (United Imaging Healthcare, Shanghai, China) with the scan commencing immediately prior to bolus injection of the radiotracer.
The first cohort comprised six participants (4 female; mean age: 67±15 years) with chronic low-back myofascial pain who underwent total-body dynamic PET, receiving bolus injections of both 18F-FDG (98±9 MBq) and 11C-butanol PET (268±6 MBq) at two scanning sessions within 14 days (clinicaltrials.gov identifier: NCT05876858). The median interval between scans was 9 days (range: 0 to 14). Two participants were scanned on the same day with 11C-butanol scanning commencing first followed by at least a 3-hour interval before 18F-FDG PET to allow 11C to decay to negligible levels. The second cohort comprised of 34 healthy participants (21 female; mean age: 51±13 years) with no self-reported history of cancer or myocardial infarction in the past 5 years.26 Participants were scanned with total-body dynamic 18F-FDG PET (mean injected activity: 358±33 MBq, bolus injection) and was used for methodological development and validation against literature blood flow ranges. Two of the six participants from the first cohort and twenty of thirty-four participants in the second cohort self-identified as belonging to racial/ethnic minorities.26
For all dynamic scans, the first two minutes were reconstructed at high temporal resolution (HTR; 60×1 s, 30×2 s) using reconstruction software provided by the vendor. This involved a time-of-flight ordered subset expectation-maximum algorithm-based reconstruction without point spread function modeling and with 4 iterations, 20 subsets, and standard corrections for attenuation, scatter, randoms, dead time, and decay.22 We used a matrix size of 150×150×486 and an isotropic voxel size of 4 mm.
Tracer Kinetic Modeling of Blood Flow from Dynamic FDG Data
Existing methods to measure blood flow with 18F-FDG have been limited to selected tissue with high extraction fraction such that the blood-to-tissue transport rate K1 approximates blood flow directly10,11 or by non-linear calibration.12 K1 is defined as the product of blood flow (F) and extraction fraction (E): Equation (1) shows that K1 is a good approximation of blood flow only when E is close to 1.
18F-FDG K1 can be measured with early dynamic imaging and a standard one-tissue compartment (S1TC) model as the phosphorylation and dephosphorylation of 18F-FDG is not identifiable during the first few minutes of the dynamic scan.13,27 The impulse response function, where vb is the blood volume, K1 and k2the blood-to-tissue transport and clearance rates, respectively. Here, the value of vb at t = 0 reflects the compartmental assumption that radiotracer instantaneously and uniformly mixes in regional blood vessels.
In reality, the radiotracer requires a non-zero transit time to traverse the length of the blood vessels at a rate equal to blood flow. This can be explicitly modeled in distributed parameter models.16,17 Here, we used the adiabatic approximation to the tissue homogeneity (AATH) model,17 a distributed kinetic model with a closed-form time-domain solution that explicitly models blood flow and a mean vascular transit time. The impulse response function, R;S1TC(t) is (Figure 1):
where F is blood flow and Tc is the mean vascular transit time for the radiotracer to pass through the length of the blood vessels. The blood volume is accordingly the product of the volumetric blood flow rate and the average time required to traverse the vascular volume (vb = FTc).
The AATH impulse response function describes a finite-time vascular phase (0 ≤ t ≤Tc) during which the radiotracer traverses the blood vessels at a rate equal to the blood flow. After this mean vascular transit time (t > Tc), radiotracer exchanges between blood and tissue like a compartment model and thus the impulse response follows an exponential decay like the S1TC model. Accordingly, the impulse response of the AATH and S1TC mainly differ by the presence of a non-zero length vascular phase in the AATH model. We expect that the AATH and S1TC fittings may perform similarly at high extraction fractions as blood flow becomes tightly correlated with K1. For a general arterial input, Ca(t), the tissue time-activity curve, Q(t), can be derived as: where td is the time delay between radiotracer arrival at the measured arterial input location and local tissue vasculature. We used a basis function method for parameter estimation using parametric forms of each model as described previously28,29 and detailed in the Supplementary Materials. The AATH model was applied on both 18F-FDG and 11C-butanol.
Image Analysis
Total-body PET enabled non-invasive measurement of an image-derived input function for kinetic analysis. The ascending aorta was used for kinetic modeling of all tissue except the lungs for which a right ventricle input function was used.24,30,31 Early 18F-FDG kinetics were quantified by analyzing regional time-activity curves obtained from tissue segmentations in 10 regions of interest (Supplementary Materials).
Total-body parametric images of early kinetics were generated by voxel-wise kinetic modeling on 4-mm isotropic reconstructions. The dynamic images and parametric images were smoothed by the kernel method, which is analogous to nonlocal means denoising.32,33 Composite image priors were derived from multiple static PET images (FDG: 0-5, 5-20, 20-40, and 40-60 minutes; 11C-butanol: 0-1, 1-2, and 2-3 minutes) and we used 49 nearest neighbors within a 9×9×9 voxel neighbourhood as in our previous work.32,33
Evaluating Time-Activity Curve Fitting
We compared the quality of the AATH and S1TC model time-activity curve fits using the Akaike information criterion (AIC):34 where Q(t) and are the measured and fitted time-activity curves, respectively, M is the number of frames, tm is the midpoint time of the mth frame, and n is the number of model parameters. The AATH model comprised n 5 parameters (td, Tc, F, K1, k2) while the S1TC had n = 4 (td, vb, K1, k2). We computed the difference in AIC (AATH minus S1TC) for each region of interest. A lower AIC indicated better fitting after accounting for the number of model parameters and the residual fitting error. Practical identifiability analysis was also performed as in previous work27 to determine the reliability of AATH parameter estimates.
To evaluate the effect of temporal resolution on the suitability of the AATH model, we frame-averaged each measured regional time-activity curve in the healthy 18F-FDG PET cohort at 1, 2, 3, 5, and 10 s frame intervals. The resampled data was fitted with the AATH and S1TC models and AIC differences were compared for each region and frame interval.
Validating 18F-FDG Blood Flow Quantification
The mean and standard deviation of regional blood flow values estimated with the AATH model were computed for all participants. In participants with both 18F-FDG and 11C-butanol PET, we performed correlation and Bland-Altman analysis35 of regional blood flow estimates between radiotracers. Total-body blood flow parametric images were visually compared between radiotracers. For the healthy 18F-FDG PET cohort, we compared their average regional values against literature ranges (summarized in Supplementary Table 1) mainly derived from flow-tracer PET.
Results
Time-Activity Curve Fitting and Model Selection
An example high-temporal resolution 18F-FDG time-activity curve fitting in the cortical grey matter with the S1TC and AATH models is shown in Figure 1. The first-pass peak, which was accurately measured with high temporal resolution dynamic imaging, was better fitted with the AATH model compared to the S1TC model. Furthermore, the peak of the intravascular component (dashed red line) of the AATH fitted curve better aligned with the peak of the measured curve. The intravascular distribution of the S1TC-fitted curve was smaller than that of the AATH model fitting, and to compensate, the extravascular distribution of the S1TC-fitted curve grew larger than that of the AATH. In all regions of interest investigated, the AATH model was preferred on average over the S1TC model across 34 high-temporal resolution dynamic 18F-FDG scans of healthy participants (Figure 2a).
Effect of Temporal Resolution on Model Selection
Figure 2b illustrates the difference in AIC between the AATH and S1TC models at different temporal resolutions and tissue regions for our 34 healthy 18F-FDG cohort. By the AIC, the AATH model had improved fitting over the S1TC model at 1 to 2 s frame intervals, though the magnitude of AIC differences was less at 2 s and a few more individual cases preferred the S1TC model. At 3 s frame interval, the AATH and S1TC models were similarly preferred, but beyond 3 s, the S1TC model was clearly preferred by the AIC.
Validation of 18F-FDG PET Blood Flow Against 11C-Butanol PET
Correlation and Bland-Altman analysis between 11C-butanol and 18F-FDG blood flow across all 6 participants, each with 10 tissue regions, are shown in Figure 3. 18F-FDG blood flow estimated with our proposed method had strong quantitative agreement with the 11C-butanol reference measurement with a Pearson correlation coefficient of 0.955 (p<0.001) and a linear regression slope and intercept of 0.973 and –0.012, respectively. The mean difference (18F-FDG minus 11C-butanol) was –0.031 ml/min/cm3, indicating that our 18F-FDG blood flow measures marginally underestimated that of 11C-butanol on average. The Bland-Altman 95% limits of agreement were –0.445 to 0.383 ml/min/cm3 with the larger differences mainly driven by higher blood flow tissues. One participant had severe intra-frame respiratory motion during the 11C-butanol scan, which prevented accurate lung blood flow quantification and substantial overestimation (>1.0 ml/min/cm3) with our 18F-FDG method. Further analysis stratified by regions with similar blood flow are shown in Supplementary Figures 1 and 2. Standard 18F-FDG K1 did not strongly agree with 11C-butanol blood flow in general (Figure 3b).
Total-body Parametric Imaging of Blood Flow with 18F-FDG
Total-body parametric images of blood flow generated with 18F-FDG and 11C-butanol in the same participant are shown in Figure 4. Parametric images appeared similar both visually and in quantitative range across the body. A notable difference observed between the two blood flow maps was the absence of sagittal and transverse sinus in the 11C-butanol parametric image. This is likely due to the high extraction fraction of 11C-butanol in the brain resulting in its lower venous concentration. One participant had substantial differences in cerebral blood flow between 11C-butanol and 18F-FDG (Supplementary Figure 3), which may be from a combination of physiological and methodological factors.
Regional 18F-FDG Early Kinetics in Healthy Participants
The distribution of blood flow estimates with our 18F-FDG method in 34 healthy participants is plotted in Figure 5. On average, all tissues were within the expected range except the subcortical grey matter and lungs, which were slightly below and above the upper range of average blood flow values reported in literature (Figure 5 and Supplementary Table 1), respectively.4 The identifiability of regional blood flow estimates with our proposed method was overall excellent (absolute mean error < 5%, standard deviation < 15%) across tissue regions except the skeletal muscle (mean overestimation of 6.4%; Supplementary Table 2).
Regional 18F-FDG extraction fraction values in the healthy cohort are summarized in Table 1. 18F-FDG extraction fraction varied greatly between tissues across the body. Accordingly, S1TC 18F-FDG K1 was in general agreement with 11C-butanol blood flow only in tissues with high extraction fraction (Supplementary Figure 4).
Discussion
We developed an early dynamic 18F-FDG PET method for total-body blood flow imaging with high-temporal resolution kinetic modeling and validated against a 11C-butanol reference in a subset of participants scanned with both radiotracers. Conventional methods for 18F-FDG blood flow imaging have been limited to tissues with relatively high extraction fraction where blood-to-tissue transport rate K1 can approximate blood flow. Our proposed method instead uses a distributed kinetic model16,17 that explicitly accounts for the blood flow delivery rate of radiotracer to blood vessels, which was resolved with high-temporal resolution dynamic imaging. 18F-FDG blood flow estimates were in quantitative agreement with 11C-butanol in direct comparisons in the same subjects (Figure 3 and Figure 4). We further validated our method in 34 healthy participants showing regional blood flow values across the body were broadly within literature ranges (Figure 5). We report first data on 18F-FDG extraction fraction values across the body, which indeed varied between tissue types (4 to 95%; Table 1). To our knowledge, this is the first study to perform total-body blood flow imaging with 18F-FDG and compare against 11C-butanol flow tracer PET in the same subjects.
Our data indeed showed that high temporal resolution of 1 to 2 s was required for the AATH model to be preferred over the S1TC model based on the AIC metric (Figure 2). The temporal resolution may need to be closer to 1 s for tissue such as the lung where the right ventricle input function often has a very fast, sharp bolus. Total-body PET now allows the requisite temporal resolution for blood flow imaging across the body using the widely available 18F-FDG radiotracer. Our method is generally applicable across the body in contrast to conventional 18F-FDG blood flow estimation methods that require high extraction fraction. Directly using 18F-FDG K1 as a surrogate of blood flow does not generalize across all tissue types (Figure 3b).
Intravenously injected tracers are delivered to tissue vasculature by blood flow. Distributed models explicitly account for this process.16,18,19 Historically, the AATH distributed model used in this study has been employed for blood flow imaging using inert contrast agents and high temporal resolution dynamic CT or MRI.36,37 We now show that distributed modeling is applicable to a noninert metabolic radiotracer (18F-FDG) as well as a freely diffusible flow radiotracer (11C-butanol) provided the requisite temporal resolution is used. This suggests that our method may be generally applicable to a wide range of tracers, enabling single-tracer multiparametric imaging of biologically and physiologically meaningful parameters, such as flow-metabolic imaging1–3 with 18F-FDG or joint quantification of blood flow and amyloid burden with amyloid PET tracers.38
Our study did not examine liver and myocardial blood flow. Kinetic modeling of the liver is complicated by its additional portal vein input,39 which could not be accurately measured due to insufficient PET spatial resolution.39 Further development is required to enable hepatic blood flow measurements with our method; however, existing methods11,39 using 18F-FDG K1 as a surrogate of hepatic blood flow may be sufficient due to the high permeability of liver sinusoids.40 Our initial analysis of the myocardium suggested that spillover from the right and left ventricles were substantial at high temporal resolution dynamic imaging of 18F-FDG first pass, resulting in substantial blood flow overestimation. We will investigate methods for correcting spillover41 and cardiac motion42 in the future for better quantification of myocardial blood flow with the proposed method.
This study had limitations. First, the sample size of participants scanned with both 18F-FDG and 11C-butanol was small in this pilot study. Instead, the validity of our early dynamic 18F-FDG PET blood flow measurements was supported by comparisons of 34 additional healthy participants against literature blood flow ranges. Additional subjects scanned with both radiotracers will be analyzed in future studies. Second, participants were not recruited specifically for validation of 18F-FDG blood flow. One participant from the dual-tracer group was suspected to have physiological differences between 11C-butanol and 18F-FDG scans. Future studies will better account for physiological confounds by measuring pCO2, pO2, and heart rate among others. Lastly, we did not study patients with major diagnosed blood flow defects such as those with peripheral, carotid, or coronary artery disease among others. Further validation is required under these disease conditions.
Conclusion
This study presented the development of an early dynamic 18F-FDG PET method with high-temporal resolution kinetic modeling for total-body blood flow imaging. Utilizing the ubiquitous 18F-FDG radiotracer for blood flow imaging may mitigate the need for a costly and practically-challenging flow-tracer PET scan. In combination with standard 18F-FDG PET methods for glucose metabolic imaging, our proposed method may allow efficient single-tracer imaging of blood flow and metabolism, resulting in lower radiation exposure to the patient, shorter scan times, and less infrastructural requirements and cost. Our method may be generally applicable to other radiotracers, broadening the possibility of single-tracer multiparametric imaging of biologically and physiologically meaningful parameters from a single dynamic PET scan.
Data Availability
All data produced in the present study are available upon reasonable request to the authors.
Disclosure
The University of California, Davis has a research agreement and a revenue sharing agreement with United Imaging Healthcare. This work was supported in part by National Institutes of Health (NIH) grants R01 EB033435 and R61 AT012187. The image data of healthy participants were acquired under the support of NIH R01 CA206187 and P30 CA093373. The other authors declare no competing interests.
Key Points
Question
Can high-temporal resolution early dynamic 18F-fluorodeoxyglucose (FDG) PET kinetic analysis be used for total-body blood flow imaging?
Pertinent Findings
Blood flow estimates between 18F-FDG and gold standard 11C-butanol PET in the same participants showed good agreement across the body. Regional blood flow measurements with the proposed early dynamic 18F-FDG PET method in 34 healthy participants were within well-established reference ranges in tissues across the body.
Implications for Patient Care
Total-body blood flow imaging can be performed with the widely available 18F-FDG radiotracer, possibly mitigating the need for a dedicated flow radiotracer and expanding opportunities to efficiently study blood flow and glucose metabolism in combination with standard 18F-FDG metabolic imaging methods.
Supplementary Materials
Kinetic Parameter Estimation
Parametric forms of the standard one-tissue compartment (S1TC) model and adiabatic approximation to the tissue homogeneity (AATH) model time-activity curves can be derived by substituting Equations (2) and (3) into (4), respectively:
We interpreted the vvbb term and FF term as the intravascular distributions of the S1TC and AATH fitted time-activity curves, respectively, and the K1 term as the extravascular distribution.
We used a basis function method1,2 for all kinetic parameter estimation on time-activity curves of the dynamic scan’s first two minutes. For the AATH model, basis functions were computed by using grid searched values of ttdd from 0 to 50 s, TTcc from 3 to 50 s, and 100 logarithmically spaced values of k2 between 6×10-4 to 15 min-1. The remaining linear parameters (F, K1) were then estimated by a non-negative linear least squares algorithm.3 A similar procedure was followed for the S1TC model but without TTcc in the grid search and linearly estimating vb and (1 – vb) K1. For both radiotracers, we assumed that whole-blood tracer activity was equal to that in blood plasma over the first two minutes of the dynamic PET scan. 11C-butanol rapidly equilibrates uniformly between blood plasma and erythrocytes4 and for 18F-FDG, blood plasma is commonly approximated by the whole-blood image-derived arterial input function.
Tissue Segmentation
The lungs, renal cortex, spleen, and skeletal muscle (splenius capitis, psoas, thigh, calves), and bone marrow in the pelvis and lumbar vertebrae were manually delineated on 3D Slicer (Version 5.2)5 by referencing a combination of the total-body CT, dynamic PET, and 0-2 minute static PET images. For the brain, we used a deep learning-based 18F-FDG-PET/CT segmentation tool6 to delineate the 83 brain regions of the Hammersmith atlas,7 which were grouped to form masks of the cortical and subcortical grey matter, white matter, brainstem, and whole cerebellum. The grey and white matter in the cerebrum were distinguished by an Otsu threshold.8 In participants with both 18F-FDG and 11C-butanol PET, FDG brain masks were resampled to the 11C-butanol-PET brain space by co-registering9 the 0-2 minute static 18F-FDG-PET brain image to that of the 11C-butanol PET. Segmentations were visually inspected and manually adjusted as needed to avoid large vessels and organ boundaries where motion and spillover were more prevalent.
Acknowledgements
The authors gratefully acknowledge Dr. Benjamin A. Spencer and Dr. Yasser G. Abdelhafez for helpful discussions and technical support, as well as the technologists and staff at the EXPLORER Molecular Imaging Center, particularly Lynda E. Painting, for their assistance in patient consent and data acquisition.
Footnotes
Financial support: This work was supported in part by National Institutes of Health (NIH) grants R01 EB033435 and R61 AT012187. The image data of healthy participants were acquired under the support of NIH R01 CA206187 and P30 CA093373.