Multi-Scale hierarchical generation of PET parametric maps: application and testing on a [11C]DPN study

Voxel level quantification of [11C]CURB, a radioligand for Fatty Acid Amide Hydrolase, using high resolution positron emission tomography

[11C]CURB is a novel irreversible radioligand for imaging fatty acid amide hydrolase in the human brain. In the present work, we validate an algorithm for generating parametric map images of [11C]CURB acquired with a high resolution research tomograph (HRRT) scanner. This algorithm applies the basis function method on an irreversible two-tissue compartment model (k4 = 0) with arterial input function, i.e., BAFPIC. Monte Carlo simulations are employed to assess bias and variability of the binding macroparameters (Ki and λk3) as a function of the voxel noise level and the range of basis functions. The results show that for a [11C]CURB time activity curve with noise levels corresponding to a voxel of an image acquired with the HRRT and reconstructed with the filtered back projection algorithm, the implementation of BAFPIC requires the use of a constant vascular fraction of tissue (5%) and a cutoff for slow frequencies (0.06 min-1). With these settings, BAFPIC maintains the probabilistic distributions of the binding macroparameters with approximately Gaussian shape and minimizes the bias and variability for large physiological ranges of the rate constants of [11C]CURB. BAFPIC reduces the variability of Ki to a third of that given by Patlak plot, the standard graphical method for irreversible radioligands. Application to real data demonstrated an excellent correlation between region of interest and BAFPIC parametric data and agreed with the simulations results. Therefore, BAFPIC with a constant vascular fraction can be used to generate parametric maps of [11C]CURB images acquired with an HRRT provided that the limits of the basis functions are carefully selected.

A Variational Bayesian inference method for parametric imaging of PET data

In dynamic Positron Emission Tomography (PET) studies, compartmental models provide the richest information on the tracer kinetics of the tissue. Inverting such models at the voxel level is however quite challenging due to the low signal-to-noise ratio of the time activity curves. In this study, we propose the use of a Variational Bayesian (VB) approach to efficiently solve this issue and thus obtain robust quantitative parametric maps. VB was adapted to the non-uniform noise distribution of PET data. Moreover, we propose a novel hierarchical scheme to define the model parameter priors directly from the images in case such information are not available from the literature, as often happens with new PET tracers. VB was initially tested on synthetic data generated using compartmental models of increasing complexity, providing accurate (%bias<2%±2%, root mean square error<15%±5%) parameter estimates. When applied to real data on a paradigmatic set of PET tracers (L-[1-11C]leucine, [11C]WAY100635 and [18F]FDG), VB was able to generate reliable parametric maps even in presence of high noise in the data (unreliable estimates<11%±5%).

Bayesian Quantification of Contrast-Enhanced Ultrasound Images With Adaptive Inclusion of an Irreversible Component

Contrast Enhanced Ultrasound (CEUS) is a sensitive imaging technique to assess tissue vascularity and it can be particularly useful in early detection and grading of arthritis. In a recent study we have shown that a Gamma-variate can accurately quantify synovial perfusion and it is flexible enough to describe many heterogeneous patterns. However, in some cases the heterogeneity of the kinetics can be such that even the Gamma model does not properly describe the curve, with a high number of outliers. In this work we apply to CEUS data the single compartment recirculation model (SCR) which takes explicitly into account the trapping of the microbubbles contrast agent by adding to the single Gamma-variate model its integral. The SCR model, originally proposed for dynamic-susceptibility magnetic resonance imaging, is solved here at pixel level within a Bayesian framework using Variational Bayes (VB). We also include the automatic relevant determination (ARD) algorithm to automatically infer the model complexity (SCR vs. Gamma model) from the data. We demonstrate that the inclusion of trapping best describes the CEUS patterns in 50% of the pixels, with the other 50% best fitted by a single Gamma. Such results highlight the necessity of the use ARD, to automatically exclude the irreversible component where not supported by the data. VB with ARD returns precise estimates in the majority of the kinetics (88% of total percentage of pixels) in a limited computational time (on average, 3.6 min per subject). Moreover, the impact of the additional trapping component has been evaluated for the differentiation of rheumatoid and non-rheumatoid patients, by means of a support vector machine classifier with backward feature selection. The results show that the trapping parameter is always present in the selected feature set, and improves the classification.

Spectral Analysis of Dynamic PET Studies: A Review of 20 Years of Method Developments and Applications

In Positron Emission Tomography (PET), spectral analysis (SA) allows the quantification of dynamic data by relating the radioactivity measured by the scanner in time to the underlying physiological processes of the system under investigation. Among the different approaches for the quantification of PET data, SA is based on the linear solution of the Laplace transform inversion whereas the measured arterial and tissue time-activity curves of a radiotracer are used to calculate the input response function of the tissue. In the recent years SA has been used with a large number of PET tracers in brain and nonbrain applications, demonstrating that it is a very flexible and robust method for PET data analysis. Differently from the most common PET quantification approaches that adopt standard nonlinear estimation of compartmental models or some linear simplifications, SA can be applied without defining any specific model configuration and has demonstrated very good sensitivity to the underlying kinetics. This characteristic makes it useful as an investigative tool especially for the analysis of novel PET tracers. The purpose of this work is to offer an overview of SA, to discuss advantages and limitations of the methodology, and to inform about its applications in the PET field.

The predictive power of brain mRNA mappings for in vivo protein density: a positron emission tomography correlation study

Substantial efforts are being spent on postmortem mRNA transcription mapping on the assumption that in vivo protein distribution can be predicted from such data. We tested this assumption by comparing mRNA transcription maps from the Allen Human Brain Atlas with reference protein concentration maps acquired with positron emission tomography (PET) in two representative systems of neurotransmission (opioid and serotoninergic). We found a tight correlation between mRNA expression and specific binding with 5-HT1A receptors measured with PET, but for opioid receptors, the correlation was weak. The discrepancy can be explained by differences in expression regulation between the two systems: transcriptional mechanisms dominate the regulation in the serotoninergic system, whereas in the opioid system proteins are further modulated after transcription. We conclude that mRNA information can be exploited for systems where translational mechanisms predominantly regulate expression. Where posttranscriptional mechanisms are important, mRNA data have to be interpreted with caution. The methodology developed here can be used for probing assumptions about the relationship of mRNA and protein in multiple neurotransmission systems.