Global-two-stage filtering of clinical PET parametric maps: Application to [11C]-(R)-PK11195

INTRODUCTION In Positron Emission Tomography (PET) quantification of physiological parameters at the voxel level may result in unreliable estimates due to the high noise of voxel time activity curves. Global-Two-Stage (GTS), an estimation technique belonging to the group of "population approaches", can be used to tackle this problem. GTS was previously tested on simulated PET data and yielded substantial improvements when compared to standard estimation approaches such as Weighted NonLinear Least Squares (WNLLS) and Basis Function Method (BFM). In this work GTS performance is assessed in a clinical context using the neuroinflammation marker [(11)C]-(R)-PK11195 applied to a cohort of Huntington's disease (HD) patients with and without symptoms. MATERIALS AND METHODS Parametric maps of binding potential (BP(ND)) of 12 normal controls (NC), 9 symptomatic and 9 presymptomatic HD patients were generated by applying a modified reference tissue model that accounts for tracer vascular activity in both reference and target tissues (SRTMV). GTS was then applied to SRTMV maps and its performance compared with that of SRTMV. Three smoothed versions of SRTMV, obtained by filtering the original SRTMV maps with Gaussian filters of 3 mm, 5 mm and 7 mm Full Width Half Maximum (FWHM), were also included in the comparison. Since striatal degeneration is the hallmark of HD, sensitivity was assessed for all methods by computing the mean of z-scores in caudate, putamen and globus pallidus in the voxel-by-voxel statistical comparison of BP(ND) between HD and NC. RESULTS Application of GTS to parametric maps brought a substantial qualitative improvement to SRTMV maps to the extent that anatomical structures often became visible. In addition, most parameter estimates that were outside the physiological range with SRTMV were corrected by GTS. GTS yielded a 2.3-fold increase in sensitivity with respect to SRTMV for the symptomatic cohort (mean of striatal z-scores of 0.76 for SRTMV and 1.79 for GTS) and an even more substantial increase for the presymptomatic cohort (mean of striatal z-scores of 0.34 for SRTMV and 0.96 for GTS). The sensitivity of GTS was similar to the one obtained with a filter of 7 mm FWHM applied to the initial SRTMV maps but GTS images were not characterized by the notable loss of resolution typical of smoothed maps. GTS, additionally, does not require to change/define settings according to the tracer and level of noise, whereas the choice of the FWHM value of the Gaussian filter normally employed in the smoothing procedure is typically arbitrary. CONCLUSIONS GTS is a powerful and robust tool for improving the quality of parametric maps in PET. The method is particularly appealing in that it can be applied to any tracer and estimation method, provided that initial estimates of the parameter vector and of its covariance are available. Although the benefits of GTS are far from being exhaustively assessed, the significant improvements obtained both on real and simulated data suggest that it could become an important tool for dynamic PET in the future.

[1]  T. Spinks,et al.  Physical characteristics of the ECAT EXACT3D positron tomograph. , 2000, Physics in medicine and biology.

[2]  Vincent J. Cunningham,et al.  Parametric Imaging of Ligand-Receptor Binding in PET Using a Simplified Reference Region Model , 1997, NeuroImage.

[3]  Alessandra Bertoldo,et al.  PET Parametric Imaging Improved by Global-Two-Stage Method , 2009, Annals of Biomedical Engineering.

[4]  Yun Zhou,et al.  Linear ridge regression with spatial constraint for generation of parametric images in dynamic positron emission tomography studies , 2001 .

[5]  Alexander Hammers,et al.  A systematic comparison of kinetic modelling methods generating parametric maps for [11C]-(R)-PK11195 , 2007, NeuroImage.

[6]  Marvin Bergsneider,et al.  Improved Parametric Image Generation Using Spatial-Temporal Analysis of Dynamic PET Studies , 2002, NeuroImage.

[7]  Alessandra Bertoldo,et al.  "Population" approach improves parameter estimation of kinetic models from dynamic PET data , 2004, IEEE Transactions on Medical Imaging.

[8]  R. Banati,et al.  Visualising microglial activation in vivo , 2002, Glia.

[9]  F. Turkheimer,et al.  Reference and target region modeling of [11C]-(R)-PK11195 brain studies. , 2007, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[10]  P G Bhide,et al.  Early and Progressive Accumulation of Reactive Microglia in the Huntington Disease Brain , 2001, Journal of neuropathology and experimental neurology.

[11]  R. P. Maguire,et al.  Consensus Nomenclature for in vivo Imaging of Reversibly Binding Radioligands , 2007, Journal of cerebral blood flow and metabolism : official journal of the International Society of Cerebral Blood Flow and Metabolism.

[12]  Hervé Boutin,et al.  Nuclear imaging of neuroinflammation: a comprehensive review of [11C]PK11195 challengers , 2008, European Journal of Nuclear Medicine and Molecular Imaging.

[13]  G. McLachlan,et al.  The EM algorithm and extensions , 1996 .

[14]  D J Brooks,et al.  Microglial activation correlates with severity in Huntington disease , 2006, Neurology.

[15]  S. Kety The theory and applications of the exchange of inert gas at the lungs and tissues. , 1951, Pharmacological reviews.

[16]  M. Reivich,et al.  THE [14C]DEOXYGLUCOSE METHOD FOR THE MEASUREMENT OF LOCAL CEREBRAL GLUCOSE UTILIZATION: THEORY, PROCEDURE, AND NORMAL VALUES IN THE CONSCIOUS AND ANESTHETIZED ALBINO RAT 1 , 1977, Journal of neurochemistry.

[17]  Akihiko Uchiyama,et al.  MAP-based kinetic analysis for voxel-by-voxel compartment model estimation: Detailed imaging of the cerebral glucose metabolism using FDG , 2006, NeuroImage.

[18]  Alexander Gerhard,et al.  Microglial activation in presymptomatic Huntington's disease gene carriers. , 2005, Brain : a journal of neurology.

[19]  F E Turkheimer,et al.  Multiresolution Analysis of Emission Tomography Images in the Wavelet Domain , 1999, Journal of cerebral blood flow and metabolism : official journal of the International Society of Cerebral Blood Flow and Metabolism.

[20]  Alexander Hammers,et al.  Three‐dimensional maximum probability atlas of the human brain, with particular reference to the temporal lobe , 2003, Human brain mapping.

[21]  F E Turkheimer,et al.  Statistical Modeling of Positron Emission Tomography Images in Wavelet Space , 2000, Journal of cerebral blood flow and metabolism : official journal of the International Society of Cerebral Blood Flow and Metabolism.

[22]  Alessandra Bertoldo,et al.  Novel Reference Region Model Reveals Increased Microglial and Reduced Vascular Binding of 11C-(R)-PK11195 in Patients with Alzheimer's Disease , 2008, Journal of Nuclear Medicine.

[23]  D. Newport,et al.  A Single Scatter Simulation Technique for Scatter Correction in 3D PET , 1996 .

[24]  Jeffrey A Fessler,et al.  Improving PET Receptor Binding Estimates from Logan Plots Using Principal Component Analysis , 2008, Journal of cerebral blood flow and metabolism : official journal of the International Society of Cerebral Blood Flow and Metabolism.

[25]  T. Hebert,et al.  A generalized EM algorithm for 3-D Bayesian reconstruction from Poisson data using Gibbs priors. , 1989, IEEE transactions on medical imaging.

[26]  Dimitris Visvikis,et al.  PET Image Denoising Using a Synergistic Multiresolution Analysis of Structural (MRI/CT) and Functional Datasets , 2008, Journal of Nuclear Medicine.

[27]  J F Boisvieux,et al.  Alternative approaches to estimation of population pharmacokinetic parameters: comparison with the nonlinear mixed-effect model. , 1984, Drug metabolism reviews.

[28]  Jeffrey A. Fessler Penalized weighted least-squares image reconstruction for positron emission tomography , 1994, IEEE Trans. Medical Imaging.

[29]  Ken D. Sauer,et al.  A local update strategy for iterative reconstruction from projections , 1993, IEEE Trans. Signal Process..

[30]  P. R. Bevington,et al.  Data Reduction and Error Analysis for the Physical Sciences , 1969 .

[31]  Marie Davidian,et al.  Nonlinear Models for Repeated Measurement Data , 1995 .

[32]  M. Mintun,et al.  A quantitative model for the in vivo assessment of drug binding sites with positron emission tomography , 1984, Annals of neurology.

[33]  D. Rubin,et al.  Maximum likelihood from incomplete data via the EM - algorithm plus discussions on the paper , 1977 .