Longitudinal Evaluation of Fatty Acid Metabolism in Normal and Spontaneously Hypertensive Rat Hearts with Dynamic MicroSPECT Imaging

The goal of this project is to develop radionuclide molecular imaging technologies using a clinical pinhole SPECT/CT scanner to quantify changes in cardiac metabolism using the spontaneously hypertensive rat (SHR) as a model of hypertensive-related pathophysiology. This paper quantitatively compares fatty acid metabolism in hearts of SHR and Wistar-Kyoto normal rats as a function of age and thereby tracks physiological changes associated with the onset and progression of heart failure in the SHR model. The fatty acid analog, 123I-labeled BMIPP, was used in longitudinal metabolic pinhole SPECT imaging studies performed every seven months for 21 months. The uniqueness of this project is the development of techniques for estimating the blood input function from projection data acquired by a slowly rotating camera that is imaging fast circulation and the quantification of the kinetics of 123I-BMIPP by fitting compartmental models to the blood and tissue time-activity curves

Estimation of the parameter covariance matrix for aone-compartment cardiac perfusion model estimated from a dynamic sequencereconstructed using map iterative reconstruction algorithms

In dynamic cardiac SPECT estimates of kinetic parameters ofa one-compartment perfusion model are usually obtained in a two stepprocess: 1) first a MAP iterative algorithm, which properly models thePoisson statistics and the physics of the data acquisition, reconstructsa sequence of dynamic reconstructions, 2) then kinetic parameters areestimated from time activity curves generated from the dynamicreconstructions. This paper provides a method for calculating thecovariance matrix of the kinetic parameters, which are determined usingweighted least squares fitting that incorporates the estimated varianceand covariance of the dynamic reconstructions. For each transaxial slicesets of sequential tomographic projections are reconstructed into asequence of transaxial reconstructions usingfor each reconstruction inthe time sequence an iterative MAP reconstruction to calculate themaximum a priori reconstructed estimate. Time-activity curves for a sumof activity in a blood region inside the left ventricle and a sum in acardiac tissue region are generated. Also, curves for the variance of thetwo estimates of the sum and for the covariance between the two ROIestimates are generated as a function of time at convergence using anexpression obtained from the fixed-point solution of the statisticalerror of the reconstruction. A one-compartment model is fit to the tissueactivity curves assuming a noisy blood input function to give weightedleast squares estimates of blood volume fraction, wash-in and wash-outrate constants specifying the kinetics of 99mTc-teboroxime for theleftventricular myocardium. Numerical methods are used to calculate thesecond derivative of the chi-square criterion to obtain estimates of thecovariance matrix for the weighted least square parameter estimates. Eventhough the method requires one matrix inverse for each time interval oftomographic acquisition, efficient estimates of the tissue kineticparameters in a dynamic cardiac SPECT study can be obtained with presentday desk-top computers

Measuring Regional Changes in the Diastolic Deformation of the Left Ventricle of SHR Rats Using microPET Technology and Hyperelastic Warping

The objective of this research was to assess applicability of a technique known as hyperelastic warping for the measurement of local strains in the left ventricle (LV) directly from microPET image data sets. The technique uses differences in image intensities between template (reference) and target (loaded) image data sets to generate a body force that deforms a finite element (FE) representation of the template so that it registers with the target images. For validation, the template image was defined as the end-systolic microPET image data set from a Wistar Kyoto (WKY) rat. The target image was created by mapping the template image using the deformation results obtained from a FE model of diastolic filling. Regression analysis revealed highly significant correlations between the simulated forward FE solution and image derived warping predictions for fiber stretch (R2 = 0.96), circumferential strain (R2 = 0.96), radial strain (R2 = 0.93), and longitudinal strain (R2 = 0.76) (p<0.001for all cases). The technology was applied to microPET image data of two spontaneously hypertensive rats (SHR) and a WKY control. Regional analysis revealed that, the lateral freewall in the SHR subjects showed the greatest deformation compared with the other wall segments. This work indicates that warping can accurately predict the strain distributions during diastole from the analysis of microPET data sets

Spectral analysis for evaluation of myocardial tracers for medical imaging

Kinetic analysis of dynamic tracer data is performed with the goal of evaluating myocardial radiotracers for cardiac nuclear medicine imaging. Data from experiments utilizing the isolated rabbit heart model are acquired by sampling the venous blood after introduction of a tracer of interest and a reference tracer. We have taken the approach that the kinetics are properly characterized by an impulse response function which describes the difference between the reference molecule (which does not leave the vasculature) and the molecule of interest which is transported across the capillary boundary and is made available to the cell. Using this formalism we can model the appearance of the tracer of interest in the venous output of the heart as a convolution of the appearance of the reference tracer with the impulse response. In this work we parameterize the impulse response function as the sum of a large number of exponential functions whose predetermined decay constants form a spectrum, and each is required only to have a nonnegative coefficient. This approach, called spectral analysis, has the advantage that it allows conventional compartmental analysis without prior knowledge of the number of compartments which the physiology may require or which the data will support

Spectral analysis for evaluation of myocardial tracers for medical imaging

Kinetic analysis of dynamic tracer data is performed with the goal of evaluating myocardial radiotracers for cardiac nuclear medicine imaging. Data from experiments utilizing the isolated rabbit heart model are acquired by sampling the venous blood after introduction of a tracer of interest and a reference tracer. We have taken the approach that the kinetics are properly characterized by an impulse response function which describes the difference between the reference molecule (which does not leave the vasculature) and the molecule of interest which is transported across the capillary boundary and is made available to the cell. Using this formalism we can model the appearance of the tracer of interest in the venous output of the heart as a convolution of the appearance of the reference tracer with the impulse response. In this work we parameterize the impulse response function as the sum of a large number of exponential functions whose predetermined decay constants form a spectrum, and each is required only to have a nonnegative coefficient. This approach, called spectral analysis, has the advantage that it allows conventional compartmental analysis without prior knowledge of the number of compartments which the physiology may require or which the data will support

Flow heterogeneity following global no-flow ischemia in isolated rabbit heart

The purpose of this study was to evaluate flow heterogeneity and impaired reflow during reperfusion following 60 min global no-flow ischemia in the isolated rabbit heart. Radiolabeled microspheres were used to measure relative flow in small left ventricular (LV) segments in five ischemia + reperfused hearts and in five non-ischemic controls. Although variable in the post-ischemic hearts, flow heterogeneity was increased relative to pre-ischemia for the whole LV (0.92 plus or minus 0.41 vs. 0.37 plus or minus 0.07, P < 0.05) as well as the subendocardium (Endo) and subepicardium (Epi) considered separately (endo: 1.28 plus or minus 0.74 vs. 0.30 plus or minus 0.09; epi: 0.69 plus or minus 0.22 vs. 0.38 plus or minus 0.08; P < 0.05 for both comparisons) during early reperfusion. There were also segments with abnormally reduced reflow. The number of segments with abnormally reduced reflow increased as flow heterogeneity increased. Abnormally reduced reflow indicates that regional ischemia can persist despite restoration of normal global flow. In addition, the relationship between regional and global flow is altered and venous outflow is derived from regions with continued perfusion and not the whole LV. These observations emphasize the need to quantify regional reflow during reperfusion following sustained no-flow ischemia in the isolated rabbit heart

Scatter correction for positron emission mammography

In this paper we present a scatter correction method for a regularized list mode maximum likelihood reconstruction algorithm for the positron emission mammograph (PEM) that is being developed at our laboratory. The scatter events inside the object are modeled as additive Poisson random variables in the forward model of the reconstruction algorithm. The mean scatter sinogram is estimated using a Monte Carlo simulation program. With the assumption that the background activity is nearly uniform, the Monte Carlo scatter simulation only needs to run once for each PEM conguration. This saves computational time. The crystal scatters are modeled as a shift-invariant blurring in image domain because they are more localized. Thus, the useful information in the crystal scatters can be deconvolved in high-resolution reconstructions. The propagation of the noise from the estimated scatter sinograminto the reconstruction is analyzed the oretically. The results provide an easy way to calculate the required number of events in the Monte Carlo scatter simulation for a given noise level in the image. The analysis is also applicable to other scatter estimation methods, provided that the covariance of the estimated scatter sinogram is available

List mode reconstruction for PET with motion compensation: A simulation study

Motion artifacts can be a significant factor that limits the image quality in high-resolution PET. Surveillance systems have been developed to track the movements of the subject during a scan. Development of reconstruction algorithms that are able to compensate for the subject motion will increase the potential of PET. In this paper we present a list mode likelihood reconstruction algorithm with the ability of motion compensation. The subject moti is explicitly modeled in the likelihood function. The detections of each detector pair are modeled as a Poisson process with time vary ingrate function. The proposed method has several advantages over the existing methods. It uses all detected events and does not introduce any interpolation error. Computer simulations show that the proposed method can compensate simulated subject movements and that the reconstructed images have no visible motion artifacts