Enhanced Extraction of Blood and Tissue Time-Activity Curves in Cardiac Mouse FDG PET Imaging by Means of Constrained Nonnegative Matrix Factorization

Otman Sarrhini; Pedro D'Orléans-Juste; Jacques A Rousseau; Jean-François Beaudoin; Roger Lecomte

doi:10.1155/2023/5366733

Enhanced Extraction of Blood and Tissue Time-Activity Curves in Cardiac Mouse FDG PET Imaging by Means of Constrained Nonnegative Matrix Factorization

Int J Biomed Imaging. 2023 Jun 15:2023:5366733. doi: 10.1155/2023/5366733. eCollection 2023.

Authors

Otman Sarrhini¹, Pedro D'Orléans-Juste², Jacques A Rousseau¹, Jean-François Beaudoin¹, Roger Lecomte^{1

3}

Affiliations

¹ Sherbrooke Molecular Imaging Center, Research Center of the Sherbrooke University Hospital (CRCHUS), Sherbrooke, QC, Canada.
² Department of Pharmacology-Physiology, Faculty of Medicine and Health Sciences, Université de Sherbrooke, Sherbrooke, QC, Canada.
³ Department of Nuclear Medicine and Radiobiology, Faculty of Medicine and Health Sciences, Université de Sherbrooke, Sherbrooke, QC, Canada.

Abstract

We propose an enhanced method to accurately retrieve time-activity curves (TACs) of blood and tissue from dynamic 2-deoxy-2-[¹⁸F]fluoro-D-glucose ([¹⁸F]FDG) positron emission tomography (PET) cardiac images of mice. The method is noninvasive and consists of using a constrained nonnegative matrix factorization algorithm (CNMF) applied to the matrix (A) containing the intensity values of the voxels of the left ventricle (LV) PET image. CNMF factorizes A into nonnegative matrices H and W, respectively, representing the physiological factors (blood and tissue) and their associated weights, by minimizing an extended cost function. We verified our method on 32 C57BL/6 mice, 14 of them with acute myocardial infarction (AMI). With CNMF, we could break down the mouse LV into myocardial and blood pool images. Their corresponding TACs were used in kinetic modeling to readily determine the [¹⁸F]FDG influx constant (K_i) required to compute the myocardial metabolic rate of glucose. The calculated K_i values using CNMF for the heart of control mice were in good agreement with those published in the literature. Significant differences in K_i values for the heart of control and AMI mice were found using CNMF. The values of the elements of W agreed well with the LV structural changes induced by ligation of the left coronary artery. CNMF was compared with the recently published method based on robust unmixing of dynamic sequences using regions of interest (RUDUR). A clear improvement of signal separation was observed with CNMF compared to the RUDUR method.