Purpose

Myocardial perfusion imaging (MPI) with single photon emission computed tomography (SPECT) is routinely used for stress testing in nuclear medicine. Recently, our group extended its potential going from 3D visual qualitative image analysis to 4D spatiotemporal reconstruction of dynamically acquired data to capture the time variation of the radiotracer concentration and the estimated myocardial blood flow (MBF) and coronary flow reserve (CFR). However, the quality of reconstructed image is compromised due to cardiac deformation and respiration. The work presented here develops an algorithm that reconstructs the dynamic sequence of separate respiratory and cardiac phases and evaluates the algorithm with data simulated with a Monte Carlo simulation for the continuous image acquisition and processing with a slowly rotating SPECT camera.

Methods

A clinically realistic Monte Carlo (MC) simulation is developed using the 4D Extended Cardiac Torso (XCAT) digital phantom with respiratory and cardiac motion to model continuous data acquisition of dynamic cardiac SPECT with slowly rotating gamma cameras by incorporating deformation and displacement of the myocardium due to cardiac and respiratory motion. We extended our previously developed 4D maximum-likelihood expectation-maximization (MLEM) reconstruction algorithm for a data set binned from a continuous list mode (LM) simulation with cardiac and respiratory information. Our spatiotemporal image reconstruction uses splines to explicitly model the temporal change of the tracer for each cardiac and respiratory gate that delineates the myocardial spatial position as the tracer washes in and out. Unlike in a fully list-mode data acquisition and reconstruction the accumulated photons are binned over a specific but very short time interval corresponding to each cardiac and respiratory gate. Reconstruction results are presented showing the dynamics of the tracer in the myocardium as it continuously deforms. These results are then compared with the conventional 4D spatiotemporal reconstruction method that models only the temporal changes of the tracer activity. Mean Stabilized Activity (MSA), signal to noise ratio (SNR) and Bias for the myocardium activities for three different target-to-background ratios (TBRs) are evaluated. Dynamic quantitative indices such as wash-in (K₁) and wash-out (k₂) rates at each gate were also estimated.

Results

The MSA and SNR are higher with higher TBRs while biases were improved with higher TBRs to less than 10%. The correlation between exhalation-inhalation sequence with the ground truth during respiratory cycle was excellent. Our reconstruction method showed better resolved myocardial walls during diastole to systole as compared to the ungated 4D image. Estimated values of K₁ and k₂ were also consistent with the ground truth.

Conclusion

The continuous image acquisition for dynamic scan using conventional two-head gamma cameras can provide valuable information for MPI. Our study demonstrated the viability of using a continuous image acquisition method on a widely used clinical two-head SPECT system. Our reconstruction method showed better resolved myocardial walls during diastole to systole as compared to the ungated 4D image. Precise implementation of reconstruction algorithms, better segmentation techniques by generating images of different tissue types and background activity would improve the feasibility of the method in real clinical environment.

Myocardial perfusion imaging (MPI) using slow rotating large field of view cameras requires spatiotemporal reconstruction of dynamically acquired data to capture the time variation of the radiotracer concentration. In vivo, MPI contains additional degrees of freedom involving unavoidable motion of the heart due to quasiperiodic beating and the effects of respiration, which can severely degrade the quality of the images. This work develops a technique for a single photon emission computed tomography (SPECT) that reconstructs the distribution of the radiotracer concentration in the myocardium using a tensor product of different sets of basis functions that approximately describe the spatiotemporal variation of the radiotracer concentration and the motion of the heart. In this study the temporal B-spline basis functions are chosen to reflect the dynamics of the radiotracer, while the intrinsic deformation and the extrinsic motion of the heart are described by a product of a discrete set of Gaussian basis functions. Reconstruction results are presented showing the dynamics of the tracer in the myocardium as it deforms due to cardiac beating, and is displaced due to respiratory motion. These results are compared with the conventional 4D-spatiotemporal reconstruction method that models only the temporal changes of the tracer activity. The higher dimensional reconstruction method proposed here improves bias, yet the signal-to-noise ratio (SNR) decreases slightly due to redistribution of the counts over the cardiac-respiratory gates. Additionally, there is a trade-off between the number of gates and the number of projections per gate to achieve high contrast images.

Background

Risk stratification and early detection of cardiac allograft vasculopathy (CAV) is essential in orthotopic heart transplantation (OHT) patients. This study assesses the changes in myocardial blood flow (MBF) noninvasively in OHT patients using quantitative cardiac PET with regadenoson.

Methods

Twelve patients (Group 1) (8 males, 4 females, mean age 55 ± 7 years) with no history of post OHT myocardial ischemia were enrolled 5.4 ± 2.0 years after OHT. Fifteen patients (Group 2) (9 males, 6 females, mean age 71 ± 9 years) with intermediate pretest probability but not documented evidence for coronary artery disease (CAD) were also included to serve as control. Global and regional MBFs were assessed using dynamic ¹³N-NH₃ PET at rest and during regadenoson-induced hyperemia. The coronary flow reserve (CFR) was also calculated as the ratio of hyperemic to resting MBF.

Results

Mean regadenoson-induced rate-pressure products were similar in both groups, while there was an increase in resting rate-pressure product in Group 1 patients. Both mean and median values of resting MBF were higher in Group 1 than Group 2 patients (1.33 ± 0.31 and 1.01 ± 0.21 mL/min/g for Groups 1 and 2, respectively, P < .001), while mean hyperemic MBF values were similar in both Groups (2.68 ± 0.84 and 2.64 ± 0.94 mL/min/g, P = NS) but median hyperemic MBF values were lower in Group 1 than Group 2 patients (2.0 vs. 2.60 mL/min/g, P = .018). Both mean and median CFR values demonstrated a significant reduction for Group 1 compared to Group 2 patients (2.07 ± 0.74 vs 2.63 ± 0.48, P = .025).

Conclusions

This study suggests that the MBF in OHT patients may be abnormal at resting state with diminished CFR. This hints that the epicardial and microvascular coronary subsystem may be exacerbated after OHT leading to the gradual progression of CAV.

Background

Risk stratification and early detection of cardiac allograft vasculopathy (CAV) is essential in orthotopic heart transplantation (OHT) patients. This study assesses the changes in myocardial blood flow (MBF) noninvasively in OHT patients using quantitative cardiac PET with regadenoson.

Methods

Twelve patients (Group 1) (8 males, 4 females, mean age 55 ± 7 years) with no history of post OHT myocardial ischemia were enrolled 5.4 ± 2.0 years after OHT. Fifteen patients (Group 2) (9 males, 6 females, mean age 71 ± 9 years) with intermediate pretest probability but not documented evidence for coronary artery disease (CAD) were also included to serve as control. Global and regional MBFs were assessed using dynamic ¹³N-NH₃ PET at rest and during regadenoson-induced hyperemia. The coronary flow reserve (CFR) was also calculated as the ratio of hyperemic to resting MBF.

Results

Mean regadenoson-induced rate-pressure products were similar in both groups, while there was an increase in resting rate-pressure product in Group 1 patients. Both mean and median values of resting MBF were higher in Group 1 than Group 2 patients (1.33 ± 0.31 and 1.01 ± 0.21 mL/min/g for Groups 1 and 2, respectively, P < .001), while mean hyperemic MBF values were similar in both Groups (2.68 ± 0.84 and 2.64 ± 0.94 mL/min/g, P = NS) but median hyperemic MBF values were lower in Group 1 than Group 2 patients (2.0 vs. 2.60 mL/min/g, P = .018). Both mean and median CFR values demonstrated a significant reduction for Group 1 compared to Group 2 patients (2.07 ± 0.74 vs 2.63 ± 0.48, P = .025).

Conclusions

This study suggests that the MBF in OHT patients may be abnormal at resting state with diminished CFR. This hints that the epicardial and microvascular coronary subsystem may be exacerbated after OHT leading to the gradual progression of CAV.

Background

SPECT myocardial perfusion imaging (MPI) is a clinical mainstay that is typically performed with static imaging protocols and visually or semi-quantitatively assessed for perfusion defects based upon the relative intensity of myocardial regions. Dynamic cardiac SPECT presents a new imaging technique based on time-varying information of radiotracer distribution, which permits the evaluation of regional myocardial blood flow (MBF) and coronary flow reserve (CFR). In this work, a preliminary feasibility study was conducted in a small patient sample designed to implement a unique combined static-dynamic single-dose one-day visit imaging protocol to compare quantitative dynamic SPECT with static conventional SPECT for improving the diagnosis of coronary artery disease (CAD).

Methods

Fifteen patients (11 males, four females, mean age 71 ± 9 years) were enrolled for a combined dynamic and static SPECT (Infinia Hawkeye 4, GE Healthcare) imaging protocol with a single dose of ^99mTc-tetrofosmin administered at rest and a single dose administered at stress in a one-day visit. Out of 15 patients, eleven had selective coronary angiography (SCA), 8 within 6 months and the rest within 24 months of SPECT imaging, without intervening symptoms or interventions. The extent and severity of perfusion defects in each myocardial region was graded visually. Dynamically acquired data were also used to estimate the MBF and CFR. Both visually graded images and estimated CFR were tested against SCA as a reference to evaluate the validity of the methods.

Results

Overall, conventional static SPECT was normal in ten patients and abnormal in five patients, dynamic SPECT was normal in 12 patients and abnormal in three patients, and CFR from dynamic SPECT was normal in nine patients and abnormal in six patients. Among those 11 patients with SCA, conventional SPECT was normal in 5, 3 with documented CAD on SCA with an overall accuracy of 64%, sensitivity of 40% and specificity of 83%. Dynamic SPECT image analysis also produced a similar accuracy, sensitivity, and specificity. CFR was normal in 6, each with CAD on SCA with an overall accuracy of 91%, sensitivity of 80%, and specificity of 100%. The mean CFR was significantly lower for SCA detected abnormal than for normal patients (3.86±1.06 vs 1.94±0. 0.67, P < 0.001).

Conclusions

The visually assessed image findings in static and dynamic SPECT are subjective, and may not reflect direct physiologic measures of coronary lesion based on SCA. The CFR measured with dynamic SPECT is fully objective, with better sensitivity and specificity, available only with the data generated from the dynamic SPECT method.

Background

SPECT myocardial perfusion imaging (MPI) is a clinical mainstay that is typically performed with static imaging protocols and visually or semi-quantitatively assessed for perfusion defects based upon the relative intensity of myocardial regions. Dynamic cardiac SPECT presents a new imaging technique based on time-varying information of radiotracer distribution, which permits the evaluation of regional myocardial blood flow (MBF) and coronary flow reserve (CFR). In this work, a preliminary feasibility study was conducted in a small patient sample designed to implement a unique combined static-dynamic single-dose one-day visit imaging protocol to compare quantitative dynamic SPECT with static conventional SPECT for improving the diagnosis of coronary artery disease (CAD).

Methods

Fifteen patients (11 males, four females, mean age 71 ± 9 years) were enrolled for a combined dynamic and static SPECT (Infinia Hawkeye 4, GE Healthcare) imaging protocol with a single dose of ^99mTc-tetrofosmin administered at rest and a single dose administered at stress in a one-day visit. Out of 15 patients, eleven had selective coronary angiography (SCA), 8 within 6 months and the rest within 24 months of SPECT imaging, without intervening symptoms or interventions. The extent and severity of perfusion defects in each myocardial region was graded visually. Dynamically acquired data were also used to estimate the MBF and CFR. Both visually graded images and estimated CFR were tested against SCA as a reference to evaluate the validity of the methods.

Results

Overall, conventional static SPECT was normal in ten patients and abnormal in five patients, dynamic SPECT was normal in 12 patients and abnormal in three patients, and CFR from dynamic SPECT was normal in nine patients and abnormal in six patients. Among those 11 patients with SCA, conventional SPECT was normal in 5, 3 with documented CAD on SCA with an overall accuracy of 64%, sensitivity of 40% and specificity of 83%. Dynamic SPECT image analysis also produced a similar accuracy, sensitivity, and specificity. CFR was normal in 6, each with CAD on SCA with an overall accuracy of 91%, sensitivity of 80%, and specificity of 100%. The mean CFR was significantly lower for SCA detected abnormal than for normal patients (3.86±1.06 vs 1.94±0. 0.67, P < 0.001).

Conclusions

The visually assessed image findings in static and dynamic SPECT are subjective, and may not reflect direct physiologic measures of coronary lesion based on SCA. The CFR measured with dynamic SPECT is fully objective, with better sensitivity and specificity, available only with the data generated from the dynamic SPECT method.

Our approach differs from the usual global measure of cardiac efficiency by using PET/MRI to measure efficiency of small pieces of cardiac tissue whose limiting size is equal to the spatial resolution of the PET scanner. We initiated a dynamic cardiac PET study immediately prior to the injection of 15.1 mCi of ¹¹C-acetate acquiring data for 25 minutes while simultaneously acquiring MRI cine data. 1) A 3D finite element (FE) biomechanical model of the imaged heart was constructed by utilizing nonrigid deformable image registration to alter the Dassault Systèmes FE Living Heart Model (LHM) to fit the geometry in the cardiac MRI cine data. The patient specific FE cardiac model with estimates of stress, strain, and work was transformed into PET/MRI format. 2) A 1-tissue compartment model was used to calculate wash-in (K₁) and the linear portion of the decay in the PET ¹¹C-acetate time activity curve (TAC) was used to calculate the wash-out k₂(mono) rate constant. K₁ was used to calculate blood flow and k₂(mono) was used to calculate myocardial volume oxygen consumption ( MVO₂ ). 3) Estimates of stress and strain were used to calculate Myocardial Equivalent Minute Work ( MEMW ) and Cardiac Efficiency = MEMW/MVO₂ was then calculated for 17 tissue segments of the left ventricle. The global MBF was 0.96 ± 0.15 ml/min/gm and MVO₂ ranged from 8 to 17 ml/100gm/min. Six central slices of the MRI cine data provided a range of MEMW of 0.1 to 0.4 joules/gm/min and a range of Cardiac Efficiency of 6 to 18%.

Long-distance intracellular axonal transport is predominantly microtubule-based, and its impairment is linked to neurodegeneration. In this study, we present theoretical arguments that suggest that near the axon boundaries (walls), the effective viscosity can become large enough to impede cargo transport in small (but not large) caliber axons. Our theoretical analysis suggests that this opposition to motion increases rapidly as the cargo approaches the wall. We find that having parallel microtubules close enough together to enable a cargo to simultaneously engage motors on more than one microtubule dramatically enhances motor activity, and thus minimizes the effects of any opposition to transport. Even if microtubules are randomly placed in axons, we find that the higher density of microtubules found in small-caliber axons increases the probability of having parallel microtubules close enough that they can be used simultaneously by motors on a cargo. The boundary effect is not a factor in transport in large-caliber axons where the microtubule density is lower.