Dual-tracer PET using generalized factor analysis of dynamic sequences

Kinetic model-informed deep learning for multiplexed PET image separation

BACKGROUND

Multiplexed positron emission tomography (mPET) imaging can measure physiological and pathological information from different tracers simultaneously in a single scan. Separation of the multiplexed PET signals within a single PET scan is challenging due to the fact that each tracer gives rise to indistinguishable 511 keV photon pairs, and thus no unique energy information for differentiating the source of each photon pair.

METHODS

Recently, many applications of deep learning for mPET image separation have been concentrated on pure data-driven methods, e.g., training a neural network to separate mPET images into single-tracer dynamic/static images. These methods use over-parameterized networks with only a very weak inductive prior. In this work, we improve the inductive prior of the deep network by incorporating a general kinetic model based on spectral analysis. The model is incorporated, along with deep networks, into an unrolled image-space version of an iterative fully 4D PET reconstruction algorithm.

RESULTS

The performance of the proposed method was evaluated on a simulated brain image dataset for dual-tracer [18 F]FDG+[11 C]MET PET image separation. The results demonstrate that the proposed method can achieve separation performance comparable to that obtained with single-tracer imaging. In addition, the proposed method outperformed the model-based separation methods (the conventional voxel-wise multi-tracer compartment modeling method (v-MTCM) and the image-space dual-tracer version of the fully 4D PET image reconstruction algorithm (IS-F4D)), as well as a pure data-driven separation [using a convolutional encoder-decoder (CED)], with fewer training examples.

CONCLUSIONS

This work proposes a kinetic model-informed unrolled deep learning method for mPET image separation. In simulation studies, the method proved able to outperform both the conventional v-MTCM method and a pure data-driven CED with less training data.

Signal separation of simultaneous dual-tracer PET imaging based on global spatial information and channel attention

BACKGROUND

Simultaneous dual-tracer positron emission tomography (PET) imaging efficiently provides more complete information for disease diagnosis. The signal separation has long been a challenge of dual-tracer PET imaging. To predict the single-tracer images, we proposed a separation network based on global spatial information and channel attention, and connected it to FBP-Net to form the FBPnet-Sep model.

RESULTS

Experiments using simulated dynamic PET data were conducted to: (1) compare the proposed FBPnet-Sep model to Sep-FBPnet model and currently existing Multi-task CNN, (2) verify the effectiveness of modules incorporated in FBPnet-Sep model, (3) investigate the generalization of FBPnet-Sep model to low-dose data, and (4) investigate the application of FBPnet-Sep model to multiple tracer combinations with decay corrections. Compared to the Sep-FBPnet model and Multi-task CNN, the FBPnet-Sep model reconstructed single-tracer images with higher structural similarity, peak signal-to-noise ratio and lower mean squared error, and reconstructed time-activity curves with lower bias and variation in most regions. Excluding the Inception or channel attention module resulted in degraded image qualities. The FBPnet-Sep model showed acceptable performance when applied to low-dose data. Additionally, it could deal with multiple tracer combinations. The qualities of predicted images, as well as the accuracy of derived time-activity curves and macro-parameters were slightly improved by incorporating a decay correction module.

CONCLUSIONS

The proposed FBPnet-Sep model was considered a potential method for the reconstruction and signal separation of simultaneous dual-tracer PET imaging.

Direct reconstruction for simultaneous dual-tracer PET imaging based on multi-task learning

BACKGROUND

Simultaneous dual-tracer positron emission tomography (PET) imaging can observe two molecular targets in a single scan, which is conducive to disease diagnosis and tracking. Since the signals emitted by different tracers are the same, it is crucial to separate each single tracer from the mixed signals. The current study proposed a novel deep learning-based method to reconstruct single-tracer activity distributions from the dual-tracer sinogram.

METHODS

We proposed the Multi-task CNN, a three-dimensional convolutional neural network (CNN) based on a framework of multi-task learning. One common encoder extracted features from the dual-tracer dynamic sinogram, followed by two distinct and parallel decoders which reconstructed the single-tracer dynamic images of two tracers separately. The model was evaluated by mean squared error (MSE), multiscale structural similarity (MS-SSIM) index and peak signal-to-noise ratio (PSNR) on simulated data and real animal data, and compared to the filtered back-projection method based on deep learning (FBP-CNN).

RESULTS

In the simulation experiments, the Multi-task CNN reconstructed single-tracer images with lower MSE, higher MS-SSIM and PSNR than FBP-CNN, and was more robust to the changes in individual difference, tracer combination and scanning protocol. In the experiment of rats with an orthotopic xenograft glioma model, the Multi-task CNN reconstructions also showed higher qualities than FBP-CNN reconstructions.

CONCLUSIONS

The proposed Multi-task CNN could effectively reconstruct the dynamic activity images of two single tracers from the dual-tracer dynamic sinogram, which was potential in the direct reconstruction for real simultaneous dual-tracer PET imaging data in future.

Simultaneous radiomethylation of [11C]harmine and [11C]DASB and kinetic modeling approach for serotonergic brain imaging in the same individual

Simultaneous characterization of pathologies by multi-tracer positron emission tomography (PET) is among the most promising applications in nuclear medicine. Aim of this work was the simultaneous production of two PET-tracers in one module and test the relevance for human application. [11C]harmine and [11C]DASB were concurrently synthesized in a 'two-in-one-pot' reaction in quality for application. Dual-tracer protocol was simulated using 16 single PET scans in different orders of tracer application separated by different time intervals. Volume of distribution was calculated for single- and dual-tracer measurements using Logan's plot and arterial input function in 13 brain regions. The 'two-in-one-pot' reaction yielded equivalent amounts of both radiotracers with comparable molar activities. The simulations of the dual-tracer application were comparable to the single bolus injections in 13 brain regions, when [11C]harmine was applied first and [11C]DASB second, with an injection time interval of 45 min (rxy = 0.90). Our study shows the successful simultaneous dual-tracer production leading to decreased radiation burden and costs. The simulation of dual subject injection to quantify the monoamine oxidase-A and serotonin transporter distribution proved its high potential. Multi-tracer imaging may drive more sophisticated study designs and diminish the day-to-day differences in the same individual as well as increase PET scanner efficiency.

Technical note: A preliminary study of dual-tracer PET image reconstruction guided by FDG and/or MR kernels

PURPOSE

Clinically, single radiotracer positron emission tomography (PET) imaging is a commonly used examination method; however, since each radioactive tracer reflects the information of only one kind of cell, it easily causes false negatives or false positives in disease diagnosis. Therefore, reasonably combining two or more radiotracers is recommended to improve the accuracy of diagnosis and the sensitivity and specificity of the disease when conditions permit.

METHODS

This paper proposes incorporating ¹⁸ F-fluorodeoxyglucose (FDG) as a higher-quality PET image to guide the reconstruction of other lower-count ¹¹ C-methionine (MET) PET datasets to compensate for the lower image quality by a popular kernel algorithm. Specifically, the FDG prior is needed to extract kernel features, and these features were used to build a kernel matrix using a k-nearest-neighbor (kNN) search for MET image reconstruction. We created a 2-D brain phantom to validate the proposed method by simulating sinogram data containing Poisson random noise and quantitatively compared the performance of the proposed FDG-guided kernelized expectation maximization (KEM) method with the performance of Gaussian and non-local means (NLM) smoothed maximum likelihood expectation maximization (MLEM), MR-guided KEM, and multi-guided-S KEM algorithms. Mismatch experiments between FDG/MR and MET data were also carried out to investigate the outcomes of possible clinical situations.

RESULTS

In the simulation study, the proposed method outperformed the other algorithms by at least 3.11% in the signal-to-noise ratio (SNR) and 0.68% in the contrast recovery coefficient (CRC), and it reduced the mean absolute error (MAE) by 8.07%. Regarding the tumor in the reconstructed image, the proposed method contained more pathological information. Furthermore, the proposed method was still superior to the MR-guided KEM method in the mismatch experiments.

CONCLUSIONS

The proposed FDG-guided KEM algorithm can effectively utilize and compensate for the tissue metabolism information obtained from dual-tracer PET to maximize the advantages of PET imaging.

Quantitation of multiple injection dynamic PET scans: an investigation of the benefits of pooling data from separate scans when mapping kinetics

Multiple injection dynamic positron emission tomography (PET) scanning is used in the clinical management of certain groups of patients and in medical research. The analysis of these studies can be approached in two ways: (i) separate analysis of data from individual tracer injections, or (ii), concatenate/pool data from separate injections and carry out a combined analysis. The simplicity of separate analysis has some practical appeal but may not be statistically efficient. We use a linear model framework associated with a kinetic mapping scheme to develop a simplified theoretical understanding of separate and combined analysis. The theoretical framework is explored numerically using both 1D and 2D simulation models. These studies are motivated by the breast cancer flow-metabolism mismatch studies involving15O-water (H2O) and18F-Fluorodeoxyglucose (FDG) and repeat15O-H2O injections used in brain activation investigations. Numerical results are found to be substantially in line with the simple theoretical analysis: mean square error characteristics of alternative methods are well described by factors involving the local voxel-level resolution of the imaging data, the relative activities of the individual scans and the number of separate injections involved. While voxel-level resolution has dependence on scan dose, after adjustment for this effect, the impact of a combined analysis is understood in simple terms associated with the linear model used for kinetic mapping. This is true for both data reconstructed by direct filtered backprojection or iterative maximum likelihood. The proposed analysis has potential to be applied to the emerging long axial field-of-view PET scanners.

PET Parametric Imaging: Past, Present, and Future

Positron emission tomography (PET) is actively used in a diverse range of applications in oncology, cardiology, and neurology. The use of PET in the clinical setting focuses on static (single time frame) imaging at a specific time-point post radiotracer injection and is typically considered as semi-quantitative; e.g. standardized uptake value (SUV) measures. In contrast, dynamic PET imaging requires increased acquisition times but has the advantage that it measures the full spatiotemporal distribution of a radiotracer and, in combination with tracer kinetic modeling, enables the generation of multiparametric images that more directly quantify underlying biological parameters of interest, such as blood flow, glucose metabolism, and receptor binding. Parametric images have the potential for improved detection and for more accurate and earlier therapeutic response assessment. Parametric imaging with dynamic PET has witnessed extensive research in the past four decades. In this paper, we provide an overview of past and present activities and discuss emerging opportunities in the field of parametric imaging for the future.

Longitudinal mouse-PET imaging: a reliable method for estimating binding parameters without a reference region or blood sampling

Abstract

Longitudinal mouse PET imaging is becoming increasingly popular due to the large number of transgenic and disease models available but faces challenges. These challenges are related to the small size of the mouse brain and the limited spatial resolution of microPET scanners, along with the small blood volume making arterial blood sampling challenging and impossible for longitudinal studies. The ability to extract an input function directly from the image would be useful for quantification in longitudinal small animal studies where there is no true reference region available such as TSPO imaging.

Methods

Using dynamic, whole-body ¹⁸F-DPA-714 PET scans (60 min) in a mouse model of hippocampal sclerosis, we applied a factor analysis (FA) approach to extract an image-derived input function (IDIF). This mouse-specific IDIF was then used for 4D-resolution recovery and denoising (4D-RRD) that outputs a dynamic image with better spatial resolution and noise properties, and a map of the total volume of distribution (V_T) was obtained using a basis function approach in a total of 9 mice with 4 longitudinal PET scans each. We also calculated percent injected dose (%ID) with and without 4D-RRD. The V_T and %ID parameters were compared to quantified ex vivo autoradiography using regional correlations of the specific binding from autoradiography against V_T and %ID parameters.

Results

The peaks of the IDIFs were strongly correlated with the injected dose (Pearson R = 0.79). The regional correlations between the %ID estimates and autoradiography were R = 0.53 without 4D-RRD and 0.72 with 4D-RRD over all mice and scans. The regional correlations between the V_T estimates and autoradiography were R = 0.66 without 4D-RRD and 0.79 with application of 4D-RRD over all mice and scans.

Conclusion

We present a FA approach for IDIF extraction which is robust, reproducible and can be used in quantification methods for resolution recovery, denoising and parameter estimation. We demonstrated that the proposed quantification method yields parameter estimates closer to ex vivo measurements than semi-quantitative methods such as %ID and is immune to tracer binding in tissue unlike reference tissue methods. This approach allows for accurate quantification in longitudinal PET studies in mice while avoiding repeated blood sampling.

Electronic supplementary material

The online version of this article (10.1007/s00259-020-04755-5) contains supplementary material, which is available to authorized users.

Perspective on potential impact of HIV central nervous system latency on eradication

: Given the challenges of life-long adherence to suppressive HIV antiretroviral therapy (ART) and possibilities of comorbidities, such as HIV association neurocognitive disorder, HIV remission and eradication are desirable goals for people living with HIV. In some individuals, there is evidence that HIV persists and replicates in the CNS, impacting the success of HIV remission interventions. This article addresses the role of HIV CNS latency on HIV eradication, examines the effects of early ART, latency-modifying agents, antibody-based and T-cell enhancing therapies on the CNS as well as ART interruption in remission studies. We propose the integration of CNS monitoring into such studies in order to clarify the short-term and long-term neurological safety of experimental agents and treatment interruption, and to better characterize their effects on HIV CNS persistence.

Three-dimensional convolutional neural networks for simultaneous dual-tracer PET imaging

Dual-tracer positron emission tomography (PET) is a promising technique to measure the distribution of two tracers in the body by a single scan, which can improve the clinical accuracy of disease diagnosis and can also serve as a research tool for scientists. Most current research on dual-tracer PET reconstruction is based on mixed images pre-reconstructed by algorithms, which restricts the further improvement of the precision of reconstruction. In this study, we present a hybrid loss-guided deep learning based framework for dual-tracer PET imaging using sinogram data, which can achieve reconstruction by naturally unifying two functions: the reconstruction of the mixed images and the separation for individual tracers. Combined with volumetric dual-tracer images, we adopted a three-dimensional (3D) convolutional neural network (CNN) to learn full features, including spatial information and temporal information simultaneously. In addition, an auxiliary loss layer was introduced to guide the reconstruction of the dual tracers. We used Monte Carlo simulations with data augmentation to generate sufficient datasets for training and testing. The results were analyzed by the bias and variance both spatially (different regions of interest) and temporally (different frames). The analysis verified the feasibility of the 3D CNN framework for dual-tracer reconstruction. Furthermore, we compared the reconstruction results with a deep belief network (DBN), which is another deep learning based technique for the separation of dual-tracer images based on time-activity curves (TACs). The comparison results provide insights into the superior features and performance of the 3D CNN. Furthermore, we tested the [¹¹C]FMZ-[¹¹C]DTBZ images with three total-counts levels ([Formula: see text], [Formula: see text], [Formula: see text]), which indicate different noise ratios. The analysis results demonstrate that our method can successfully recover the respective distribution of lower total counts with nearly the same accuracy as that of the higher total counts in the total counts range we applied, which also also indicates the proposed 3D CNN framework is more robust to noise compared with DBN.

Dynamic whole-body PET imaging: principles, potentials and applications

PURPOSE

In this article, we discuss dynamic whole-body (DWB) positron emission tomography (PET) as an imaging tool with significant clinical potential, in relation to conventional standard uptake value (SUV) imaging.

BACKGROUND

DWB PET involves dynamic data acquisition over an extended axial range, capturing tracer kinetic information that is not available with conventional static acquisition protocols. The method can be performed within reasonable clinical imaging times, and enables generation of multiple types of PET images with complementary information in a single imaging session. Importantly, DWB PET can be used to produce multi-parametric images of (i) Patlak slope (influx rate) and (ii) intercept (referred to sometimes as "distribution volume"), while also providing (iii) a conventional 'SUV-equivalent' image for certain protocols.

RESULTS

We provide an overview of ongoing efforts (primarily focused on FDG PET) and discuss potential clinically relevant applications.

CONCLUSION

Overall, the framework of DWB imaging [applicable to both PET/CT(computed tomography) and PET/MRI (magnetic resonance imaging)] generates quantitative measures that may add significant value to conventional SUV image-derived measures, with limited pitfalls as we also discuss in this work.

Multi-Tracer Guided PET Image Reconstruction

Multi-tracer positron emission tomography (PET) has the potential to enhance PET imaging by providing complementary information from different physiological processes. However, one or more of the images may present high levels of noise. Guided image reconstruction methods transfer information from a guide image into the PET image reconstruction to encourage edge-preserving noise reduction. In this work we aim to reduce noise in poorer quality PET datasets via guidance from higher quality ones by using a weighted quadratic penalty approach. In particular, we applied this methodology to [18F]fluorodeoxyglucose (FDG) and [11C]methionine imaging of gliomas. 3D simulation studies showed that guiding the reconstruction of methionine datasets using pre-existing FDG images reduced reconstruction errors across the whole-brain (-8%) and within a tumour (-36%) compared to maximum likelihood expectation-maximisation (MLEM). Furthermore, guided reconstruction outperformed a comparable non-local means filter, indicating that regularising during reconstruction is preferable to post-reconstruction approaches. Hyperparameters selected from the 3D simulation study were applied to real data, where it was observed that the proposed FDG-guided methionine reconstruction allows for better edge preservation and noise reduction than standard MLEM. Overall, the results in this work demonstrate that transferring information between datasets in multi-tracer PET studies improves image quality and quantification performance.

[Kinetic cluster and α-divergence-based dynamic myocardial factorial analysis of positron-emission computed tomography images]

OBJECTIVE

We purpose a novel factor analysis method based on kinetic cluster and α-divergence measure for extracting the blood input function and the time-activity curve of the regional tissue from dynamic myocardial positron emission computed tomography(PET) images.

METHODS

Dynamic PET images were decomposed into initial factors and factor images by minimizing the α-divergence between the factor model and actual image data. The kinetic clustering as a priori constraint was then incorporated into the model to solve the nonuniqueness problem, and the tissue time-activity curves and the tissue space distributions with physiological significance were generated.

RESULTS

The model was applied to the 82RbPET myocardial perfusion simulation data and compared with the traditional model-based least squares measure and the minimal spatial overlap constraint. The experimental results showed that the proposed model performed better than the traditional model in terms of both accuracy and sensitivity.

CONCLUSION

This method can select the optimal measure by α value, and incorporate the prior information of the kinetic clustering of PET image pixels to obtain the accurate time-activity curves of the tissue, which has shown good performance in visual evaluation and quantitative evaluation.

Impact of contamination with long-lived radionuclides on PET kinetics modelling in multitracer studies

INTRODUCTION

An important issue in multitracer studies is the separation of signals from the different radiotracers. This is especially the case when an early tracer has a long physical half-life and kinetic modelling has to be performed, because the early tracer can confer a long-lived contaminating background not only to images but also to a measured input function derived from blood samples. In this study, we examined data from a sequential multitracer infection study involving In (t1/2=2.8 days), investigating the influence on gamma counting of blood samples and on the kinetic modelling of subsequent PET tracers. Blood sample counts were corrected by recounting the samples a few days later. A more optimal choice of energy window was also explored. The effect of correction versus noncorrection was investigated using a two-tissue kinetic model with irreversible uptake (K1, k2, k3).

RESULTS

K1 was least affected and k3 was most affected by the contamination, corresponding to the effect being relatively larger on the late part of the blood input function. A narrower energy window reduced the problem, but this will not be possible for all types of contaminating background.

CONCLUSION

Gamma counting of blood samples can lead to a contaminating background not observed in PET imaging and this background can affect kinetic modelling. If the contaminating tracer has a much longer half-life than the foreground tracer, then the problem can be solved by late recounting of the samples.

Positron Emission Tomography: Current Challenges and Opportunities for Technological Advances in Clinical and Preclinical Imaging Systems

Positron emission tomography (PET) imaging is based on detecting two time-coincident high-energy photons from the emission of a positron-emitting radioisotope. The physics of the emission, and the detection of the coincident photons, give PET imaging unique capabilities for both very high sensitivity and accurate estimation of the in vivo concentration of the radiotracer. PET imaging has been widely adopted as an important clinical modality for oncological, cardiovascular, and neurological applications. PET imaging has also become an important tool in preclinical studies, particularly for investigating murine models of disease and other small-animal models. However, there are several challenges to using PET imaging systems. These include the fundamental trade-offs between resolution and noise, the quantitative accuracy of the measurements, and integration with X-ray computed tomography and magnetic resonance imaging. In this article, we review how researchers and industry are addressing these challenges.

Clustering-initiated factor analysis application for tissue classification in dynamic brain positron emission tomography

The goal is to quantify the fraction of tissues that exhibit specific tracer binding in dynamic brain positron emission tomography (PET). It is achieved using a new method of dynamic image processing: clustering-initiated factor analysis (CIFA). Standard processing of such data relies on region of interest analysis and approximate models of the tracer kinetics and of tissue properties, which can degrade accuracy and reproducibility of the analysis. Clustering-initiated factor analysis allows accurate determination of the time-activity curves and spatial distributions for tissues that exhibit significant radiotracer concentration at any stage of the emission scan, including the arterial input function. We used this approach in the analysis of PET images obtained using (11)C-Pittsburgh Compound B in which specific binding reflects the presence of β-amyloid. The fraction of the specific binding tissues determined using our approach correlated with that computed using the Logan graphical analysis. We believe that CIFA can be an accurate and convenient tool for measuring specific binding tissue concentration and for analyzing tracer kinetics from dynamic images for a variety of PET tracers. As an illustration, we show that four-factor CIFA allows extraction of two blood curves and the corresponding distributions of arterial and venous blood from PET images even with a coarse temporal resolution.

EM reconstruction of dual isotope PET using staggered injections and prompt gamma positron emitters

PURPOSE

The aim of dual isotope positron emission tomography (DIPET) is to create two separate images of two coinjected PET radiotracers. DIPET shortens the duration of the study, reduces patient discomfort, and produces perfectly coregistered images compared to the case when two radiotracers would be imaged independently (sequential PET studies). Reconstruction of data from such simultaneous acquisition of two PET radiotracers is difficult because positron decay of any isotope creates only 511 keV photons; therefore, the isotopes cannot be differentiated based on the detected energy.

METHODS

Recently, the authors have proposed a DIPET technique that uses a combination of radiotracer A which is a pure positron emitter (such as(18)F or (11)C) and radiotracer B in which positron decay is accompanied by the emission of a high-energy (HE) prompt gamma (such as (38)K or (60)Cu). Events that are detected as triple coincidences of HE gammas with the corresponding two 511 keV photons allow the authors to identify the lines-of-response (LORs) of isotope B. These LORs are used to separate the two intertwined distributions, using a dedicated image reconstruction algorithm. In this work the authors propose a new version of the DIPET EM-based reconstruction algorithm that allows the authors to include an additional, independent estimate of radiotracer A distribution which may be obtained if radioisotopes are administered using a staggered injections method. In this work the method is tested on simple simulations of static PET acquisitions.

RESULTS

The authors' experiments performed using Monte-Carlo simulations with static acquisitions demonstrate that the combined method provides better results (crosstalk errors decrease by up to 50%) than the positron-gamma DIPET method or staggered injections alone.

CONCLUSIONS

The authors demonstrate that the authors' new EM algorithm which combines information from triple coincidences with prompt gammas and staggered injections improves the accuracy of DIPET reconstructions for static acquisitions so they reach almost the benchmark level calculated for perfectly separated tracers.