Voxelwise Principal Component Analysis of Dynamic [S-Methyl-(11)C]Methionine PET Data in Glioma Patients

SIMPLE SUMMARY: Recent works on dynamic amino acid positron emission tomography (PET) imaging of gliomas have highlighted characteristic behaviors of time-activity curves (TACs) extracted from the whole tumor w.r.t. the grade, genotype, and outcome. However, gliomas are known to be highly heterogeneous tumors. Here, we aim at highlighting similar dynamic behaviors at the voxel level within the tumor volume in [S-methyl-(11)C]methionine PET data of 33 glioma patients using principal component analysis (PCA). The PCA model was derived from TACs of 20 patients and subsequently applied to 13 other patients in whom our approach was shown to outperform classical pharmacokinetic modeling to this end. Our parameter-free approach provides additional parametric maps from dynamic methionine PET scans with little modification of the routine protocol and no arterial sampling. This early methodological work paves the way for various clinical studies on glioma heterogeneity with applications for treatment planning and response evaluation. ABSTRACT: Recent works have demonstrated the added value of dynamic amino acid positron emission tomography (PET) for glioma grading and genotyping, biopsy targeting, and recurrence diagnosis. However, most of these studies are based on hand-crafted qualitative or semi-quantitative features extracted from the mean time activity curve within predefined volumes. Voxelwise dynamic PET data analysis could instead provide a better insight into intra-tumor heterogeneity of gliomas. In this work, we investigate the ability of principal component analysis (PCA) to extract relevant quantitative features from a large number of motion-corrected [S-methyl-(11)C]methionine ([(11)C]MET) PET frames. We first demonstrate the robustness of our methodology to noise by means of numerical simulations. We then build a PCA model from dynamic [(11)C]MET acquisitions of 20 glioma patients. In a distinct cohort of 13 glioma patients, we compare the parametric maps derived from our PCA model to these provided by the classical one-compartment pharmacokinetic model (1TCM). We show that our PCA model outperforms the 1TCM to distinguish characteristic dynamic uptake behaviors within the tumor while being less computationally expensive and not requiring arterial sampling. Such methodology could be valuable to assess the tumor aggressiveness locally with applications for treatment planning and response evaluation. This work further supports the added value of dynamic over static [(11)C]MET PET in gliomas.