Dynamic PET image reconstruction integrating temporal regularization associated with respiratory motion correction for applications in oncology

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.

Improved Clinical Workflow for Whole-Body Patlak Parametric Imaging Using Two Short Dynamic Acquisitions

Objective

We sought to explore the feasibility of shorter acquisition times using two short dynamic scans for a multiparametric PET study and the influence of quantitative performance in shortened dynamic PET.

Methods

Twenty-one patients underwent whole-body dynamic ¹⁸F-FDG PET/CT examinations on a PET/CT (Siemens Biograph Vision) with a total scan time of 75 min using continuous bed motion for Patlak multiparametric imaging. Two sets of Patlak multiparametric images were produced: the standard MR_FDG and DV_FDG images (MR_FDG-_std and DV_FDG-_std) and two short dynamic MR_FDG and DV_FDG images (MR_FDG-_tsd and DV_FDG-_tsd), which were generated by a 0–75 min post injection (p.i.) dynamic PET series and a 0–6 min + 60–75 min p.i. dynamic PET series, respectively. The maximum, mean, and peak values of the standard and two short dynamic multiparametric acquisitions were obtained and compared using Passing–Bablok regression and Bland–Altman analysis.

Results

High correlations were obtained between MR_FDG-_tsd and MR_FDG-_std, and between DV_FDG-_tsd and DV_FDG-_std for both normal organs and all lesions (0.962 ≦ Spearman’s rho ≦ 0.982, p < 0.0001). The maximum, mean, and peak values of the standard and two short dynamic multiparametric acquisitions were also in agreement. For normal organs, the Bland–Altman plot showed that the mean bias of MR_FDG-max, MR_FDG-mean, and MR_FDG-peak was -0.002 (95% CI: -0.032–0.027), -0.002 (95% CI: -0.026–0.023), and -0.002 (95% CI: -0.026–0.022), respectively. The mean bias of DV_FDG-max, DV_FDG-mean, and DV_FDG-peak was -3.3 (95% CI: -24.8–18.2), -1.4 (95% CI: -12.1–9.2), and -2.3 (95% CI: -15–10.4), respectively. For lesions, the Bland–Altman plot showed that the mean bias of MR_FDG-max, MR_FDG-mean, and MR_FDG-peak was -0.009 (95% CI: -0.056–0.038), -0.004 (95% CI: -0.039–0.031), and -0.004 (95% CI: -0.036–0.028), respectively. The mean bias of DV_FDG-max, DV_FDG-mean, and DV_FDG-peak was -8.4 (95% CI: -42.6–25.9), -4.8 (95% CI: -20.2–10.6), and -4.0 (95% CI: -23.7–15.6), respectively.

Conclusions

This study demonstrates the feasibility of using two short dynamic scans that include the first 0–6 min and 60–75 min scans p.i. for Patlak multiparametric images, which can increase patient throughout for parametric analysis.

Hepatic Positron Emission Tomography: Applications in Metabolism, Haemodynamics and Cancer

Evaluating in vivo the metabolic rates of the human liver has been a challenge due to its unique perfusion system. Positron emission tomography (PET) represents the current gold standard for assessing non-invasively tissue metabolic rates in vivo. Here, we review the existing literature on the assessment of hepatic metabolism, haemodynamics and cancer with PET. The tracer mainly used in metabolic studies has been [¹⁸F]2-fluoro-2-deoxy-D-glucose (¹⁸F-FDG). Its application not only enables the evaluation of hepatic glucose uptake in a variety of metabolic conditions and interventions, but based on the kinetics of ¹⁸F-FDG, endogenous glucose production can also be assessed. 14(R,S)-[¹⁸F]fluoro-6-thia-Heptadecanoic acid (¹⁸F-FTHA), ¹¹C-Palmitate and ¹¹C-Acetate have also been applied for the assessment of hepatic fatty acid uptake rates (¹⁸F-FTHA and ¹¹C-Palmitate) and blood flow and oxidation (¹¹C-Acetate). Oxygen-15 labelled water (¹⁵O-H₂O) has been used for the quantification of hepatic perfusion. ¹⁸F-FDG is also the most common tracer used for hepatic cancer diagnostics, whereas ¹¹C-Acetate has also shown some promising applications in imaging liver malignancies. The modelling approaches used to analyse PET data and also the challenges in utilizing PET in the assessment of hepatic metabolism are presented.

Short 2-[18F]Fluoro-2-Deoxy-D-Glucose PET Dynamic Acquisition Protocol to Evaluate the Influx Rate Constant by Regional Patlak Graphical Analysis in Patients With Non-Small-Cell Lung Cancer

Purpose: To test a short 2-[¹⁸F]Fluoro-2-deoxy-D-glucose (2-[¹⁸F]FDG) PET dynamic acquisition protocol to calculate K_i using regional Patlak graphical analysis in patients with non-small-cell lung cancer (NSCLC).

Methods: 24 patients with NSCLC who underwent standard dynamic 2-[¹⁸F]FDG acquisitions (60 min) were randomly divided into two groups. In group 1 (n = 10), a population-based image-derived input function (pIDIF) was built using a monoexponential trend (10–60 min), and a leave-one-out cross-validation (LOOCV) method was performed to validate the pIDIF model. In group 2 (n = 14), K_i was obtained by standard regional Patlak plot analysis using IDIF (0–60 min) and tissue response (10–60 min) curves from the volume of interests (VOIs) placed on descending thoracic aorta and tumor tissue, respectively. Moreover, with our method, the Patlak analysis was performed to obtain K_i,s using IDIF_Fitted curve obtained from PET counts (0–10 min) followed by monoexponential coefficients of pIDIF (10–60 min) and tissue response curve obtained from PET counts at 10 min and between 40 and 60 min, simulating two short dynamic acquisitions. Both IDIF and IDIF_Fitted curves were modeled to assume the value of 2-[¹⁸F]FDG plasma activity measured in the venous blood sampling performed at 45 min in each patient. Spearman's rank correlation, coefficient of determination, and Passing–Bablok regression were used for the comparison between K_i and K_i,s. Finally, K_i,s was obtained with our method in a separate group of patients (group 3, n = 8) that perform two short dynamic acquisitions.

Results: Population-based image-derived input function (10–60 min) was modeled with a monoexponential curve with the following fitted parameters obtained in group 1: a = 9.684, b = 16.410, and c = 0.068 min⁻¹. The LOOCV error was 0.4%. In patients of group 2, the mean values of K_i and K_i,s were 0.0442 ± 0.0302 and 0.33 ± 0.0298, respectively (R² = 0.9970). The Passing–Bablok regression for comparison between K_i and K_i,s showed a slope of 0.992 (95% CI: 0.94–1.06) and intercept value of −0.0003 (95% CI: −0.0033–0.0011).

Conclusions: Despite several practical limitations, like the need to position the patient twice and to perform two CT scans, our method contemplates two short 2-[¹⁸F]FDG dynamic acquisitions, a population-based input function model, and a late venous blood sample to obtain robust and personalized input function and tissue response curves and to provide reliable regional K_i estimation.

Quantitative PET in the 2020s: a roadmap

Positron emission tomography (PET) plays an increasingly important role in research and clinical applications, catalysed by remarkable technical advances and a growing appreciation of the need for reliable, sensitive biomarkers of human function in health and disease. Over the last 30 years, a large amount of the physics and engineering effort in PET has been motivated by the dominant clinical application during that period, oncology. This has led to important developments such as PET/CT, whole-body PET, 3D PET, accelerated statistical image reconstruction, and time-of-flight PET. Despite impressive improvements in image quality as a result of these advances, the emphasis on static, semi-quantitative 'hot spot' imaging for oncologic applications has meant that the capability of PET to quantify biologically relevant parameters based on tracer kinetics has not been fully exploited. More recent advances, such as PET/MR and total-body PET, have opened up the ability to address a vast range of new research questions, from which a future expansion of applications and radiotracers appears highly likely. Many of these new applications and tracers will, at least initially, require quantitative analyses that more fully exploit the exquisite sensitivity of PET and the tracer principle on which it is based. It is also expected that they will require more sophisticated quantitative analysis methods than those that are currently available. At the same time, artificial intelligence is revolutionizing data analysis and impacting the relationship between the statistical quality of the acquired data and the information we can extract from the data. In this roadmap, leaders of the key sub-disciplines of the field identify the challenges and opportunities to be addressed over the next ten years that will enable PET to realise its full quantitative potential, initially in research laboratories and, ultimately, in clinical practice.

FBP-Net for direct reconstruction of dynamic PET images

Dynamic positron emission tomography (PET) imaging can provide information about metabolic changes over time, used for kinetic analysis and auxiliary diagnosis. Existing deep learning-based reconstruction methods have too many trainable parameters and poor generalization, and require mass data to train the neural network. However, obtaining large amounts of medical data is expensive and time-consuming. To reduce the need for data and improve the generalization of network, we combined the filtered back-projection (FBP) algorithm with neural network, and proposed FBP-Net which could directly reconstruct PET images from sinograms instead of post-processing the rough reconstruction images obtained by traditional methods. The FBP-Net contained two parts: the FBP part and the denoiser part. The FBP part adaptively learned the frequency filter to realize the transformation from the detector domain to the image domain, and normalized the coarse reconstruction images obtained. The denoiser part merged the information of all time frames to improve the quality of dynamic PET reconstruction images, especially the early time frames. The proposed FBP-Net was performed on simulation and real dataset, and the results were compared with the state-of-art U-net and DeepPET. The results showed that FBP-Net did not tend to overfit the training set and had a stronger generalization.

Real-time control of respiratory motion: Beyond radiation therapy

Motion management in radiation oncology is an important aspect of modern treatment planning and delivery. Special attention has been paid to control respiratory motion in recent years. However, other medical procedures related to both diagnosis and treatment are likely to benefit from the explicit control of breathing motion. Quantitative imaging - including increasingly important tools in radiology and nuclear medicine - is among the fields where a rapid development of motion control is most likely, due to the need for quantification accuracy. Emerging treatment modalities like focussed-ultrasound tumor ablation are also likely to benefit from a significant evolution of motion control in the near future. In the present article an overview of available respiratory motion systems along with ongoing research in this area is provided. Furthermore, an attempt is made to envision some of the most expected developments in this field in the near future.

Partial-volume correction in dynamic PET-CT: effect on tumor kinetic parameter estimation and validation of simplified metrics

Background

Partial-volume effects generally result in an underestimation of tumor tracer uptake on PET-CT for small lesions, necessitating partial-volume correction (PVC) for accurate quantification. However, investigation of PVC in dynamic oncological PET studies to date is scarce. The aim of this study was to investigate PVC’s impact on tumor kinetic parameter estimation from dynamic PET-CT acquisitions and subsequent validation of simplified semi-quantitative metrics. Ten patients with EGFR-mutated non-small cell lung cancer underwent dynamic ¹⁸F-fluorothymidine PET-CT before, 7 days after, and 28 days after commencing treatment with a tyrosine kinase inhibitor. Parametric PVC was applied using iterative deconvolution without and with highly constrained backprojection (HYPR) denoising, respectively. Using an image-derived input function with venous parent plasma calibration, we estimated full kinetic parameters V_T, K₁, and k₃/k₄ (BP_ND) using a reversible two-tissue compartment model, and simplified metrics (SUV and tumor-to-blood ratio) at 50–60 min post-injection.

Results

PVC had a non-linear effect on measured activity concentrations per timeframe. PVC significantly changed each kinetic parameter, with a median increase in V_T of 11.8% (up to 25.1%) and 10.8% (up to 21.7%) without and with HYPR, respectively. Relative changes in kinetic parameter estimates vs. simplified metrics after applying PVC were poorly correlated (correlations 0.36–0.62; p < 0.01). PVC increased correlations between simplified metrics and V_T from 0.82 and 0.81 (p < 0.01) to 0.90 and 0.88 (p < 0.01) for SUV and TBR, respectively, albeit non-significantly. PVC also increased correlations between treatment-induced changes in simplified metrics vs. V_T at 7 (SUV) and 28 (SUV and TBR) days after treatment start non-significantly. Delineation on partial-volume corrected PET images resulted in a median decrease in metabolic tumor volume of 14.3% (IQR − 22.1 to − 7.5%), and increased the effect of PVC on kinetic parameter estimates.

Conclusion

PVC has a significant impact on tumor kinetic parameter estimation from dynamic PET-CT data, which differs from its effect on simplified metrics. However, it affected validation of these simplified metrics both as single measurements and as biomarkers of treatment response only to a small extent. Future dynamic PET studies should preferably incorporate PVC.

Trial registration

Dutch Trial Register, NTR3557.

Electronic supplementary material

The online version of this article (10.1186/s13550-019-0483-z) contains supplementary material, which is available to authorized users.