Image-derived and arterial blood sampled input functions for quantitative PET imaging of the angiotensin II subtype 1 receptor in the kidney

Automated extraction of the arterial input function from brain images for parametric PET studies

BACKGROUND

Accurate measurement of the arterial input function (AIF) is crucial for parametric PET studies, but the AIF is commonly derived from invasive arterial blood sampling. It is possible to use an image-derived input function (IDIF) obtained by imaging a large blood pool, but IDIF measurement in PET brain studies performed on standard field of view scanners is challenging due to lack of a large blood pool in the field-of-view. Here we describe a novel automated approach to estimate the AIF from brain images.

RESULTS

Total body 18F-FDG PET data from 12 subjects were split into a model adjustment group (n = 6) and a validation group (n = 6). We developed an AIF estimation framework using wavelet-based methods and unsupervised machine learning to distinguish arterial and venous activity curves, compared to the IDIF from the descending aorta. All of the automatically extracted AIFs in the validation group had similar shape to the IDIF derived from the descending aorta IDIF. The average area under the curve error and normalised root mean square error across validation data were - 1.59 ± 2.93% and 0.17 ± 0.07.

CONCLUSIONS

Our automated AIF framework accurately estimates the AIF from brain images. It reduces operator-dependence, and could facilitate the clinical adoption of parametric PET.

A short 18F-FDG imaging window triple injection neuroimaging protocol for parametric mapping in PET

BACKGROUND

In parametric PET, kinetic parameters are extracted from dynamic PET images. It is not commonly used in clinical practice because of long scan times and the requirement for an arterial input function (AIF). To address these limitations, we designed an 18F-fluorodeoxyglucose (18F-FDG) triple injection dynamic PET protocol for brain imaging with a standard field of view PET scanner using a 24-min imaging window and an input function modeled using measurements from a region of interest placed over the left ventricle.

METHODS

To test the protocol in 6 healthy participants, we examined the quality of voxel-based maps of kinetic parameters in the brain generated using the two-tissue compartment model and compared estimated parameter values with previously published values. We also utilized data from a 36-min validation imaging window to compare (1) the modeled AIF against the input function measured in the validation window; and (2) the net influx rate ([Formula: see text]) computed using parameter estimates from the short imaging window against the net influx rate obtained using Patlak analysis in the validation window.

RESULTS

Compared to the AIF measured in the validation window, the input function estimated from the short imaging window achieved a mean area under the curve error of 9%. The voxel-wise Pearson's correlation between [Formula: see text] estimates from the short imaging window and the validation imaging window exceeded 0.95.

CONCLUSION

The proposed 24-min triple injection protocol enables parametric 18F-FDG neuroimaging with noninvasive estimation of the AIF from cardiac images using a standard field of view PET scanner.

Medicinal (Radio) Chemistry: Building Radiopharmaceuticals for the Future

Radiopharmaceuticals are increasingly playing a leading role in diagnosing, monitoring, and treating disease. In comparison with conventional pharmaceuticals, the development of radiopharmaceuticals does follow the principles of medicinal chemistry in the context of imaging-altered physiological processes. The design of a novel radiopharmaceutical has several steps similar to conventional drug discovery and some particularity. In the present work, we revisited the insights of medicinal chemistry in the current radiopharmaceutical development giving examples in oncology, neurology, and cardiology. In this regard, we overviewed the literature on radiopharmaceutical development to study overexpressed targets such as prostate-specific membrane antigen and fibroblast activation protein in cancer; β-amyloid plaques and tau protein in brain disorders; and angiotensin II type 1 receptor in cardiac disease. The work addresses concepts in the field of radiopharmacy with a special focus on the potential use of radiopharmaceuticals for nuclear imaging and theranostics.

An update on the use of image-derived input functions for human PET studies: new hopes or old illusions?

BACKGROUND

The need for arterial blood data in quantitative PET research limits the wider usability of this imaging method in clinical research settings. Image-derived input function (IDIF) approaches have been proposed as a cost-effective and non-invasive alternative to gold-standard arterial sampling. However, this approach comes with its own limitations-partial volume effects and radiometabolite correction among the most important-and varying rates of success, and the use of IDIF for brain PET has been particularly troublesome.

MAIN BODY

This paper summarizes the limitations of IDIF methods for quantitative PET imaging and discusses some of the advances that may make IDIF extraction more reliable. The introduction of automated pipelines (both commercial and open-source) for clinical PET scanners is discussed as a way to improve the reliability of IDIF approaches and their utility for quantitative purposes. Survey data gathered from the PET community are then presented to understand whether the field's opinion of the usefulness and validity of IDIF is improving. Finally, as the introduction of next-generation PET scanners with long axial fields of view, ultra-high sensitivity, and improved spatial and temporal resolution, has also brought IDIF methods back into the spotlight, a discussion of the possibilities offered by these state-of-the-art scanners-inclusion of large vessels, less partial volume in small vessels, better description of the full IDIF kinetics, whole-body modeling of radiometabolite production-is included, providing a pathway for future use of IDIF.

CONCLUSION

Improvements in PET scanner technology and software for automated IDIF extraction may allow to solve some of the major limitations associated with IDIF, such as partial volume effects and poor temporal sampling, with the exciting potential for accurate estimation of single kinetic rates. Nevertheless, until individualized radiometabolite correction can be performed effectively, IDIF approaches remain confined at best to a few tracers.

Quantitation of dynamic total-body PET imaging: recent developments and future perspectives

BACKGROUND

Positron emission tomography (PET) scanning is an important diagnostic imaging technique used in disease diagnosis, therapy planning, treatment monitoring, and medical research. The standardized uptake value (SUV) obtained at a single time frame has been widely employed in clinical practice. Well beyond this simple static measure, more detailed metabolic information can be recovered from dynamic PET scans, followed by the recovery of arterial input function and application of appropriate tracer kinetic models. Many efforts have been devoted to the development of quantitative techniques over the last couple of decades.

CHALLENGES

The advent of new-generation total-body PET scanners characterized by ultra-high sensitivity and long axial field of view, i.e., uEXPLORER (United Imaging Healthcare), PennPET Explorer (University of Pennsylvania), and Biograph Vision Quadra (Siemens Healthineers), further stimulates valuable inspiration to derive kinetics for multiple organs simultaneously. But some emerging issues also need to be addressed, e.g., the large-scale data size and organ-specific physiology. The direct implementation of classical methods for total-body PET imaging without proper validation may lead to less accurate results.

CONCLUSIONS

In this contribution, the published dynamic total-body PET datasets are outlined, and several challenges/opportunities for quantitation of such types of studies are presented. An overview of the basic equation, calculation of input function (based on blood sampling, image, population or mathematical model), and kinetic analysis encompassing parametric (compartmental model, graphical plot and spectral analysis) and non-parametric (B-spline and piece-wise basis elements) approaches is provided. The discussion mainly focuses on the feasibilities, recent developments, and future perspectives of these methodologies for a diverse-tissue environment.

Total-body parametric imaging using the patlak model: Feasibility of reduced scan time

PURPOSE

This study explored the feasibility of reducing the scan time of Patlak parametric imaging on the uEXPLORER.

METHODS

A total of 65 patients (27 females and 38 males, age 56.1±10.4) were recruited in this study. ¹⁸ F-FDG was injected and its dose was adjusted by body weight (4.07 MBq / kg).Total-body dynamic scanning was performed on the uEXPLORER Total-Body PET/CT scanner with a total scan time of 60 minutes from the injection. The image derived input function (IDIF) was obtained from the aortic arch. The voxelwise Patlak analysis was applied to generate the K_i images designated as G_IDIF with different acquisition times (20-60, 30-60, 40-60, and 44-60 min). The population-based input function (PBIF) was constructed from the mean value of the IDIF from the population, and K_i images designated as G_PBIF were generated using the PBIF. Non-localmeans (NLM) denoising was applied to the generated images to get two extra groups of (NLM-designated) images: G_IDIF+NLM and G_PBIF+NLM . Two radiologists evaluated the overall image quality, noise, and lesion detectability of the K_i images from different groups. The 20-60 min scans in G_IDIF were selected as the gold standard for each patient. We determined that image quality is at sufficient level if all the lesions can be recognized and meet the clinical criteria. K_i values in muscle and lesion were compared across different groups to evaluate the quantitative accuracy.

RESULTS

The overall image quality, image noise, and lesion conspicuity were significantly better in long time series than short time series in all 4 groups (all p<0.001). The K_i images in the G_IDIF and G_PBIF groups generated from 30-min scans showed diagnostic value equivalent to the 40-min scans of G_IDIF . While the image quality of the 16-min scans was poor, all lesions could still be detected. No significant difference was found between K_i values estimated with G_IDIF and G_PBIF in muscle and lesion regions (all p>0.5). After applying the NLM filter, The coefficient of variation could be reduced on the order of [1%, 15%, 19%,37%] and [110%, 125%, 94%, 69%] with four acquisition time schemes for lesion and muscle. The reduction percentage did not have a substantial difference in IDIF and PBIF group. The K_i images in the G_IDIF+NLM and G_PBIF+NLM groups generated from the 20-min acquisitions showed acceptable quality. All lesions could be found on the NLM processed images of the 16-min scans. No significant difference was found between K_i values produced with G_IDIF+NLM and G_PBIF+NLM in muscle and lesion regions(all p>0.7).

CONCLUSIONS

The K_i images generated by the PBIF-based Patlak model using a 20-min dynamic scan with the NLM filter achieved a similar diagnostic efficiency to images with G_IDIF from 40-min dynamic data, and there is no significant difference between K_i images generated using IDIF or PBIF (p>0.5). This article is protected by copyright. All rights reserved.

Whole-body Parametric Imaging of FDG PET using uEXPLORER with Reduced Scan Time

Parametric imaging of the net influx rate (K_i) in ¹⁸F-FDG PET has been shown to provide improved quantification and specificity for cancer detection compared with SUV imaging. Current methods of generating parametric images usually require a long dynamic scanning time. With the recently developed uEXPLORER scanner, a dramatic increase in sensitivity has reduced the noise in dynamic imaging, making it more robust to use a nonlinear estimation method and flexible protocols. In this work, we explored 2 new possible protocols besides the standard 60-min one for the possibility of reducing scanning time for K_i imaging. Methods: The gold standard protocol (protocol 1) was conventional dynamic scanning with a 60-min scanning time. The first proposed protocol (protocol 2) included 2 scanning periods: 0–4 min and 54–60 min after injection. The second proposed protocol (protocol 3) consisted of a single scanning period from 50 to 60 min after injection, with a second injection applied at 56 min. The 2 new protocols were simulated from the 60-min standard scans. A hybrid input function combining the population-based input function and the image-derived input function (IDIF) was used. The results were also compared with the IDIF acquired from protocol 1. A previously developed maximum-likelihood approach was used to estimate the K_i images. In total, 7 cancer patients imaged using the uEXPLORER scanner were enrolled in this study. Lesions were identified from the patient data, and the lesion K_i values were compared among the different protocols. Results: The acquired hybrid input function was comparable in shape to the IDIF for each patient. The average difference in area under the curve was about 3%, suggesting good quantitative accuracy. The visual difference between the K_i images generated using IDIF and those generated using the hybrid input function was also minimal. The acquired K_i images using different protocols were visually comparable. The average K_i difference in the lesions was 2.8% ± 2.1% for protocol 2 and 1% ± 2.2% for protocol 3. Conclusion: The results suggest that it is possible to acquire K_i images using the nonlinear estimation approach with a much-reduced scanning time. Between the 2 new protocols, the protocol with dual injection shows the greatest promise in terms of practicality.

Minimally invasive quantification of cerebral P2X7R occupancy using dynamic [18F]JNJ-64413739 PET and MRA-driven image derived input function

[¹⁸F]JNJ-64413739 has been evaluated as PET-ligand for in vivo quantification of purinergic receptor subtype 7 receptor (P2X7R) using Logan graphical analysis with a metabolite-corrected arterial plasma input function. In the context of a P2X7R PET dose occupancy study, we evaluated a minimally invasive approach by limiting arterial sampling to baseline conditions. Meanwhile, post dose distribution volumes (V_T) under blocking conditions were estimated by combining baseline blood to plasma ratios and metabolite fractions with an MR angiography driven image derived input function (IDIF). Regional postdose V_T,IDIF values were compared with corresponding V_T,AIF estimates using a arterial input function (AIF), in terms of absolute values, test–retest reliability and receptor occupancy. Compared to an invasive AIF approach, postdose V_T,IDIF values and corresponding receptor occupancies showed only limited bias (Bland–Altman analysis: 0.06 ± 0.27 and 3.1% ± 6.4%) while demonstrating a high correlation (Spearman ρ = 0.78 and ρ = 0.98 respectively). In terms of test–retest reliability, regional intraclass correlation coefficients were 0.98 ± 0.02 for V_T,IDIF compared to 0.97 ± 0.01 for V_T,AIF. These results confirmed that a postdose IDIF, guided by MR angiography and using baseline blood and metabolite data, can be considered for accurate [¹⁸F]JNJ-64413739 PET quantification in a repeated PET study design, thus avoiding multiple invasive arterial sampling and increasing dosing flexibility.

Simplified protocol for whole-body Patlak parametric imaging with 18 F-FDG PET/CT: Feasibility and error analysis

PURPOSE

Parametric imaging using the Patlak model has been shown to provide improved lesion detectability and specificity. The Patlak model requires both tissue time-activity curves (TACs) after equilibrium and knowledge of the input function from the start of injection. Therefore, the conventional dynamic scanning protocol typically starts from the radiotracer injection all the way to equilibrium. In this paper, we propose the use of hybrid population-based and model-based input function estimation and evaluate its use for whole-body Patlak analysis, in order to reduce the total scan time and simplify clinical Patlak parametric imaging protocols. Possible quantitative errors caused by the simplified scanning protocol were also analyzed both theoretically and with the use of clinical data.

MATERIALS AND METHODS

Clinical data from 24 patients referred for tumor staging were included in this study. The patients underwent a whole-body dynamic PET study, 20 min after FDG injection (0.13 mCi/kg). The proposed whole-body scanning protocol includes 6 passes with 4-5 bed positions, depending on the size of the patient, with 2 min for each bed position. An input function from the literature was selected as the shape of the population-based input function. The descending aorta from the corresponding CT image was segmented and applied on the reconstructed dynamic PET images to acquire an image-based input function, which was later fitted using an exponential model. Due to the late scan time, only the later portion of the input function was available, which was used to scale the population-based input function. The hybrid input function was used to derive the whole-body Patlak images. Assuming a given error in the population-based input function, its influence on the final Patlak images were also derived theoretically and verified using the clinical data sets. Finally, the image quality of the reconstructed Patlak slope image was evaluated by an experienced radiologist in four different aspects: image artifacts, image noise, lesion sharpness, and lesion detectability.

RESULTS

It was found that errors in the population-based input function only affect the absolute scale of the Patlak slope image. The induced error is proportional to the percentage area-under-curve (AUC) error in the input function. These findings were also confirmed by numerical analysis. The predicted global scale was in good agreement with results from both image-based Patlak and direct Patlak approach. The fractions of the AUC from the early portion population-based input function were also found to be around 18% of the total AUC of the input function, further limiting the propagation of quantitation error from population-based input function to the final Patlak slope image. The reconstructed Patlak images were also found by the radiologist to provide excellent confidence in lesion detection tasks.

CONCLUSIONS

We have proposed a simplified whole-body scanning protocol that utilizes both population-based input function and model-based input function. The error from the population-based function was found to only affect the global scale and the overall quantitative impact can be predicted using our proposed formulas.

Total-Body Quantitative Parametric Imaging of Early Kinetics of 18F-FDG

Parametric imaging has been shown to provide better quantitation physiologically than SUV imaging in PET. With the increased sensitivity from a recently developed total-body PET scanner, whole-body scans with higher temporal resolution become possible for dynamic analysis and parametric imaging. In this paper, we focus on deriving the parameter k ₁ using compartmental modeling and on developing a method to acquire whole-body ¹⁸F-FDG PET parametric images using only the first 90 s of the postinjection scan data with the total-body PET system. Methods: Dynamic projections were acquired with a time interval of 1 s for the first 30 s and a time interval of 2 s for the following minute. Image-derived input functions were acquired from the reconstructed dynamic sequences in the ascending aorta. A 1-tissue-compartment model with 4 parameters (k ₁, k ₂, blood fraction, and delay time) was used. A maximum-likelihood-based estimation method was developed with the 1-tissue-compartment model solution. The accuracy of the acquired parameters was compared with the ones estimated using a 2-tissue-compartment irreversible model with 1-h-long data. Results: All 4 parametric images were successfully calculated using data from 2 volunteers. By comparing the time-activity curves acquired from the volumes of interest, we showed that the parameters estimated using our method were able to predict the time-activity curves of the early dynamics of ¹⁸F-FDG in different organs. The delay-time effects for different organs were also clearly visible in the reconstructed delay-time image with delay variations of as large as 40 s. The estimated parameters using both 90-s data and 1-h data agreed well for k ₁ and blood fraction, whereas a large difference in k ₂ was found between the 90-s and 1-h data, suggesting k ₂ cannot be reliably estimated from the 90-s scan. Conclusion: We have shown that with total-body PET and the increased sensitivity, it is possible to estimate parametric images based on the very early dynamics after ¹⁸F-FDG injection. The estimated k ₁ might potentially be used clinically as an indicator for identifying abnormalities.

Novel 18F-Labeled PET Imaging Agent FV45 Targeting the Renin-Angiotensin System

Renin-angiotensin system (RAS) plays an important role in the regulation of blood pressure and hormonal balance. Using positron emission tomography (PET) technology, it is possible to monitor the physiological and pathological distribution of angiotensin II type 1 receptors (AT₁), which reflects the functionality of RAS. A new ¹⁸F-labeled PET tracer derived from the clinically used AT₁ antagonist valsartan showing the least possible chemical alteration from the valsartan structure has been designed and synthesized with several strategies, which can be applied for the syntheses of further derivatives. Radioligand binding study showed that the cold reference FV45 (K _i 14.6 nM) has almost equivalent binding affinity as its lead valsartan (K _i 11.8 nM) and angiotensin II (K _i 1.7 nM). Successful radiolabeling of FV45 in a one-pot radiofluorination followed by the deprotection procedure with 21.8 ± 8.5% radiochemical yield and >99% radiochemical purity (n = 5) enabled a distribution study in rats and opened a path to straightforward large-scale production. A fast and clear kidney uptake could be observed, and this renal uptake could be selectively blocked by pretreatment with AT₁-selective antagonist valsartan. Overall, as the first ¹⁸F-labeled PET tracer based on a derivation from clinically used drug valsartan with almost identical chemical structure, [¹⁸F]FV45 will be a new tool for assessing the RAS function by visualizing AT₁ receptor distributions and providing further information regarding cardiovascular system malfunction as well as possible applications in inflammation research and cancer diagnosis.