Direct estimation of kinetic parametric images for dynamic PET.

Dynamic positron emission tomography (PET) can monitor spatiotemporal distribution of radiotracer in vivo. The spatiotemporal information can be used to estimate parametric images of radiotracer kinetics that are of physiological and biochemical interests. Direct estimation of parametric images from raw projection data allows accurate noise modeling and has been shown to offer better image quality than conventional indirect methods, which reconstruct a sequence of PET images first and then perform tracer kinetic modeling pixel-by-pixel. Direct reconstruction of parametric images has gained increasing interests with the advances in computing hardware. Many direct reconstruction algorithms have been developed for different kinetic models. In this paper we review the recent progress in the development of direct reconstruction algorithms for parametric image estimation. Algorithms for linear and nonlinear kinetic models are described and their properties are discussed

Current and Future Use of Long Axial Field-of-View Positron Emission Tomography/Computed Tomography Scanners in Clinical Oncology

The latest technical development in the field of positron emission tomography/computed tomography (PET/CT) imaging has been the extension of the PET axial field-of-view. As a result of the increased number of detectors, the long axial field-of-view (LAFOV) PET systems are not only characterized by a larger anatomical coverage but also by a substantially improved sensitivity, compared with conventional short axial field-of-view PET systems. In clinical practice, this innovation has led to the following optimization: (1) improved overall image quality, (2) decreased duration of PET examinations, (3) decreased amount of radioactivity administered to the patient, or (4) a combination of any of the above. In this review, novel applications of LAFOV PET in oncology are highlighted and future directions are discussed.</p

Direct inference of Patlak parametric images in whole-body PET/CT imaging using convolutional neural networks

Purpose: This study proposed and investigated the feasibility of estimating Patlak-derived influx rate constant (Ki) from standardized uptake value (SUV) and/or dynamic PET image series. Methods: Whole-body 18F-FDG dynamic PET images of 19 subjects consisting of 13 frames or passes were employed for training a residual deep learning model with SUV and/or dynamic series as input and Ki-Patlak (slope) images as output. The training and evaluation were performed using a nine-fold cross-validation scheme. Owing to the availability of SUV images acquired 60 min post-injection (20 min total acquisition time), the data sets used for the training of the models were split into two groups: “With SUV” and “Without SUV.” For “With SUV” group, the model was first trained using only SUV images and then the passes (starting from pass 13, the last pass, to pass 9) were added to the training of the model (one pass each time). For this group, 6 models were developed with input data consisting of SUV, SUV plus pass 13, SUV plus passes 13 and 12, SUV plus passes 13 to 11, SUV plus passes 13 to 10, and SUV plus passes 13 to 9. For the “Without SUV” group, the same trend was followed, but without using the SUV images (5 models were developed with input data of passes 13 to 9). For model performance evaluation, the mean absolute error (MAE), mean error (ME), mean relative absolute error (MRAE%), relative error (RE%), mean squared error (MSE), root mean squared error (RMSE), peak signal-to-noise ratio (PSNR), and structural similarity index (SSIM) were calculated between the predicted Ki-Patlak images by the two groups and the reference Ki-Patlak images generated through Patlak analysis using the whole acquired data sets. For specific evaluation of the method, regions of interest (ROIs) were drawn on representative organs, including the lung, liver, brain, and heart and around the identified malignant lesions. Results: The MRAE%, RE%, PSNR, and SSIM indices across all patients were estimated as 7.45 ± 0.94%, 4.54 ± 2.93%, 46.89 ± 2.93, and 1.00 ± 6.7 × 10−7, respectively, for models predicted using SUV plus passes 13 to 9 as input. The predicted parameters using passes 13 to 11 as input exhibited almost similar results compared to the predicted models using SUV plus passes 13 to 9 as input. Yet, the bias was continuously reduced by adding passes until pass 11, after which the magnitude of error reduction was negligible. Hence, the predicted model with SUV plus passes 13 to 9 had the lowest quantification bias. Lesions invisible in one or both of SUV and Ki-Patlak images appeared similarly through visual inspection in the predicted images with tolerable bias. Conclusion: This study concluded the feasibility of direct deep learning-based approach to estimate Ki-Patlak parametric maps without requiring the input function and with a fewer number of passes. This would lead to shorter acquisition times for WB dynamic imaging with acceptable bias and comparable lesion detectability performance.</p

Feasibility of using abbreviated scan protocols with population-based input functions for accurate kinetic modeling of [18F]-FDG datasets from a long axial FOV PET scanner.

BACKGROUND Accurate kinetic modeling of 18F-fluorodeoxyglucose ([18F]-FDG) positron emission tomography (PET) data requires accurate knowledge of the available tracer concentration in the plasma during the scan time, known as the arterial input function (AIF). The gold standard method to derive the AIF requires collection of serial arterial blood samples, but the introduction of long axial field of view (LAFOV) PET systems enables the use of non-invasive image-derived input functions (IDIFs) from large blood pools such as the aorta without any need for bed movement. However, such protocols require a prolonged dynamic PET acquisition, which is impractical in a busy clinical setting. Population-based input functions (PBIFs) have previously shown potential in accurate Patlak analysis of [18F]-FDG datasets and can enable the use of shortened dynamic imaging protocols. Here, we exploit the high sensitivity and temporal resolution of a LAFOV PET system and explore the use of PBIF with abbreviated protocols in [18F]-FDG total body kinetic modeling. METHODS Dynamic PET data were acquired in 24 oncological subjects for 65 min following the administration of [18F]-FDG. IDIFs were extracted from the descending thoracic aorta, and a PBIF was generated from 16 datasets. Five different scaled PBIFs (sPBIFs) were generated by scaling the PBIF with the AUC of IDIF curve tails using various portions of image data (35-65, 40-65, 45-65, 50-65, and 55-65 min post-injection). The sPBIFs were compared with the IDIFs using the AUCs and Patlak Ki estimates in tumor lesions and cerebral gray matter. Patlak plot start time (t*) was also varied to evaluate the performance of shorter acquisitions on the accuracy of Patlak Ki estimates. Patlak Ki estimates with IDIF and t* = 35 min were used as reference, and mean bias and precision (standard deviation of bias) were calculated to assess the relative performance of different sPBIFs. A comparison of parametric images generated using IDIF and sPBIFs was also performed. RESULTS There was no statistically significant difference between AUCs of the IDIF and sPBIFs (Wilcoxon test: P > 0.05). Excellent agreement was shown between Patlak Ki estimates obtained using sPBIF and IDIF. Using the sPBIF55-65 with the Patlak model, 20 min of PET data (i.e., 45 to 65 min post-injection) achieved  0.99 and peak signal-to-noise ratio > 55 dB. CONCLUSION We demonstrate the feasibility of performing accurate [18F]-FDG Patlak analysis using sPBIFs with only 20 min of PET data from a LAFOV PET scanner

Methodological considerations in quantification of oncological FDG PET studies

Contains fulltext : 87741.pdf (publisher's version ) (Closed access) Contains fulltext : 87741-1.pdf (postprint version ) (Open Access)PURPOSE: This review aims to provide insight into the factors that influence quantification of glucose metabolism by FDG PET images in oncology as well as their influence on repeated measures studies (i.e. treatment response assessment), offering improved understanding both for clinical practice and research. METHODS: Structural PubMed searches have been performed for the many factors affecting quantification of glucose metabolism by FDG PET. Review articles and references lists have been used to supplement the search findings. RESULTS: Biological factors such as fasting blood glucose level, FDG uptake period, FDG distribution and clearance, patient motion (breathing) and patient discomfort (stress) all influence quantification. Acquisition parameters should be adjusted to maximize the signal to noise ratio without exposing the patient to a higher than strictly necessary radiation dose. This is especially challenging in pharmacokinetic analysis, where the temporal resolution is of significant importance. The literature is reviewed on the influence of attenuation correction on parameters for glucose metabolism, the effect of motion, metal artefacts and contrast agents on quantification of CT attenuation-corrected images. Reconstruction settings (analytical versus iterative reconstruction, post-reconstruction filtering and image matrix size) all potentially influence quantification due to artefacts, noise levels and lesion size dependency. Many region of interest definitions are available, but increased complexity does not necessarily result in improved performance. Different methods for the quantification of the tissue of interest can introduce systematic and random inaccuracy. CONCLUSIONS: This review provides an up-to-date overview of the many factors that influence quantification of glucose metabolism by FDG PET.01 juli 201

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-[18F]Fluoro-2-deoxy-D-glucose (2-[18F]FDG) PET dynamic acquisition protocol to calculate Ki using regional Patlak graphical analysis in patients with non-small-cell lung cancer (NSCLC). Methods: 24 patients with NSCLC who underwent standard dynamic 2-[18F]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), Ki 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 Ki,s using IDIFFitted 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 IDIFFitted curves were modeled to assume the value of 2-[18F]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 Ki and Ki,s. Finally, Ki,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−1. The LOOCV error was 0.4%. In patients of group 2, the mean values of Ki and Ki,s were 0.0442 ± 0.0302 and 0.33 ± 0.0298, respectively (R2 = 0.9970). The Passing–Bablok regression for comparison between Ki and Ki,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-[18F]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 Ki estimation

Quantitative PET image reconstruction employing nested expectation-maximization deconvolution for motion compensation

Bulk body motion may randomly occur during PET acquisitions introducing blurring, attenuation-emission mismatches and, in dynamic PET, discontinuities in the measured time activity curves between consecutive frames. Meanwhile, dynamic PET scans are longer, thus increasing the probability of bulk motion. In this study, we propose a streamlined 3D PET motion-compensated image reconstruction (3D-MCIR) framework, capable of robustly deconvolving intra-frame motion from a static or dynamic 3D sinogram. The presented 3D-MCIR methods need not partition the data into multiple gates, such as 4D MCIR algorithms, or access list-mode (LM) data, such as LM MCIR methods, both associated with increased computation or memory resources. The proposed algorithms can support compensation for any periodic and non-periodic motion, such as cardio-respiratory or bulk motion, the latter including rolling, twisting or drifting. Inspired from the widely adopted point-spread function (PSF) deconvolution 3D PET reconstruction techniques, here we introduce an image-based 3D generalized motion deconvolution method within the standard 3D maximum-likelihood expectation-maximization (ML-EM) reconstruction framework. In particular, we initially integrate a motion blurring kernel, accounting for every tracked motion within a frame, as an additional MLEM modeling component in the image space (integrated 3D-MCIR). Subsequently, we replaced the integrated model component with a nested iterative Richardson-Lucy (RL) image-based deconvolution method to accelerate the MLEM algorithm convergence rate (RL-3D-MCIR). The final method was evaluated with realistic simulations of whole-body dynamic PET data employing the XCAT phantom and real human bulk motion profiles, the latter estimated from volunteer dynamic MRI scans. In addition, metabolic uptake rate Ki parametric images were generated with the standard Patlak method. Our results demonstrate significant improvement in contrast-to-noise ratio (CNR) and noise-bias performance in both dynamic and parametric images. The proposed nested RL-3D-MCIR method is implemented on the Software for Tomographic Image Reconstruction (STIR) open-source platform and is scheduled for public release

FDG-PET Quantification of Lung Inflammation with Image-Derived Blood Input Function in Mice

Dynamic FDG-PET imaging was used to study inflammation in lungs of mice following administration of a virulent strain of Klebsiella (K.) pneumoniae. Net whole-lung FDG influx constant (Ki) was determined in a compartment model using an image-derived blood input function. Methods. K. pneumoniae (~3 x 105 CFU) was intratracheally administered to six mice with 6 other mice serving as controls. Dynamic FDG-PET and X-Ray CT scans were acquired 24 hr after K. pneumoniae administration. The experimental lung time activity curves were fitted to a 3-compartment FDG model to obtain Ki. Following imaging, lungs were excised and immunohistochemistry analysis was done to assess the relative presence of neutrophils and macrophages. Results. Mean Ki for control and K. pneumoniae infected mice were (5.1 ± 1.2) ×10−3 versus (11.4 ± 2.0) ×10−3 min−1, respectively, revealing a 2.24 fold significant increase (P = 0.0003) in the rate of FDG uptake in the infected lung. Immunohistochemistry revealed that cellular lung infiltrate was almost exclusively neutrophils. Parametric Ki maps by Patlak analysis revealed heterogeneous inflammatory foci within infected lungs. Conclusion. The kinetics of FDG uptake in the lungs of mice can be noninvasively quantified by PET with a 3-compartment model approach based on an image-derived input function

Radiogenomic analysis of primary breast cancer reveals [18F]-fluorodeoxglucose dynamic flux-constants are positively associated with immune pathways and outperform static uptake measures in associating with glucose metabolism

Background: PET imaging of 18F-fluorodeoxygucose (FDG) is used widely for tumour staging and assessment of treatment response, but the biology associated with FDG uptake is still not fully elucidated. We therefore carried out gene set enrichment analyses (GSEA) of RNA sequencing data to find KEGG pathways associated with FDG uptake in primary breast cancers. Methods: Pre-treatment data were analysed from a window-of-opportunity study in which 30 patients underwent static and dynamic FDG-PET and tumour biopsy. Kinetic models were fitted to dynamic images, and GSEA was performed for enrichment scores reflecting Pearson and Spearman coefficients of correlations between gene expression and imaging. Results: A total of 38 pathways were associated with kinetic model flux-constants or static measures of FDG uptake, all positively. The associated pathways included glycolysis/gluconeogenesis (‘GLYC-GLUC’) which mediates FDG uptake and was associated with model flux-constants but not with static uptake measures, and 28 pathways related to immune-response or inflammation. More pathways, 32, were associated with the flux-constant K of the simple Patlak model than with any other imaging index. Numbers of pathways categorised as being associated with individual micro-parameters of the kinetic models were substantially fewer than numbers associated with flux-constants, and lay around levels expected by chance. Conclusions: In pre-treatment images GLYC-GLUC was associated with FDG kinetic flux-constants including Patlak K, but not with static uptake measures. Immune-related pathways were associated with flux-constants and static uptake. Patlak K was associated with more pathways than were the flux-constants of more complex kinetic models. On the basis of these results Patlak analysis of dynamic FDG-PET scans is advantageous, compared to other kinetic analyses or static imaging, in studies seeking to infer tumour-to-tumour differences in biology from differences in imaging. Trial registration NCT01266486, December 24th 2010