Multi-modality imaging for assessment of the microcirculation in peripheral artery disease: Bench to clinical practice.

Peripheral artery disease (PAD) is a highly prevalent disorder with a high risk of mortality and amputation despite the introduction of novel medical and procedural treatments. Microvascular disease (MVD) is common among patients with PAD, and despite the established role as a predictor of amputations and mortality, MVD is not routinely assessed as part of current standard practice. Recent pre-clinical and clinical perfusion and molecular imaging studies have confirmed the important role of MVD in the pathogenesis and outcomes of PAD. The recent advancements in the imaging of the peripheral microcirculation could lead to a better understanding of the pathophysiology of PAD, and result in improved risk stratification, and our evaluation of response to therapies. In this review, we will discuss the current understanding of the anatomy and physiology of peripheral microcirculation, and the role of imaging for assessment of perfusion in PAD, and the latest advancements in molecular imaging. By highlighting the latest advancements in multi-modality imaging of the peripheral microcirculation, we aim to underscore the most promising imaging approaches and highlight potential research opportunities, with the goal of translating these approaches for improved and personalized management of PAD in the future.

Time and event-specific deep learning for personalized risk assessment after cardiac perfusion imaging.

Standard clinical interpretation of myocardial perfusion imaging (MPI) has proven prognostic value for predicting major adverse cardiovascular events (MACE). However, personalizing predictions to a specific event type and time interval is more challenging. We demonstrate an explainable deep learning model that predicts the time-specific risk separately for all-cause death, acute coronary syndrome (ACS), and revascularization directly from MPI and 15 clinical features. We train and test the model internally using 10-fold hold-out cross-validation (n = 20,418) and externally validate it in three separate sites (n = 13,988) with MACE follow-ups for a median of 3.1 years (interquartile range [IQR]: 1.6, 3.6). We evaluate the model using the cumulative dynamic area under receiver operating curve (cAUC). The best model performance in the external cohort is observed for short-term prediction - in the first six months after the scan, mean cAUC for ACS and all-cause death reaches 0.76 (95% confidence interval [CI]: 0.75, 0.77) and 0.78 (95% CI: 0.78, 0.79), respectively. The model outperforms conventional perfusion abnormality measures at all time points for the prediction of death in both internal and external validations, with improvement increasing gradually over time. Individualized patient explanations are visualized using waterfall plots, which highlight the contribution degree and direction for each feature. This approach allows the derivation of individual event probability as a function of time as well as patient- and event-specific risk explanations that may help draw attention to modifiable risk factors. Such a method could help present post-scan risk assessments to the patient and foster shared decision-making.

Direct Risk Assessment From Myocardial Perfusion Imaging Using Explainable Deep Learning.

BACKGROUND: Myocardial perfusion imaging (MPI) is frequently used to provide risk stratification, but methods to improve the accuracy of these predictions are needed. OBJECTIVES: The authors developed an explainable deep learning (DL) model (HARD MACE [major adverse cardiac events]-DL) for the prediction of death or nonfatal myocardial infarction (MI) and validated its performance in large internal and external testing groups. METHODS: Patients undergoing single-photon emission computed tomography MPI were included, with 20,401 patients in the training and internal testing group (5 sites) and 9,019 in the external testing group (2 different sites). HARD MACE-DL uses myocardial perfusion, motion, thickening, and phase polar maps combined with age, sex, and cardiac volumes. The primary outcome was all-cause mortality or nonfatal MI. Prognostic accuracy was evaluated using area under the receiver-operating characteristic curve (AUC). RESULTS: During internal testing, patients with normal perfusion and elevated HARD MACE-DL risk were at higher risk than patients with abnormal perfusion and low HARD MACE-DL risk (annualized event rate, 2.9% vs 1.2%; P < 0.001). Patients in the highest quartile of HARD MACE-DL score had an annual rate of death or MI (4.8%) 10-fold higher than patients in the lowest quartile (0.48% per year). In external testing, the AUC for HARD MACE-DL (0.73; 95% CI: 0.71-0.75) was higher than a logistic regression model (AUC: 0.70), stress total perfusion deficit (TPD) (AUC: 0.65), and ischemic TPD (AUC: 0.63; all P < 0.01). Calibration, a measure of how well predicted risk matches actual risk, was excellent in both groups (Brier score, 0.079 for internal and 0.070 for external). CONCLUSIONS: The DL model predicts death or MI directly from MPI, by estimating patient-level risk with good calibration and improved accuracy compared with traditional quantitative approaches. The model incorporates mechanisms to explain to the physician which image regions contribute to the adverse event prediction.

Novel Cardiac Computed Tomography Methods for the Assessment of Anthracycline Induced Cardiotoxicity.

Anthracyclines are among the most frequently utilized anti-cancer therapies; however, their use is frequently associated with off-target cardiotoxic effects. Cardiac computed tomography (CCT) is a validated and rapidly evolving technology for the evaluation of cardiac structures, coronary anatomy and plaque, cardiac function and preprocedural planning. However, with emerging new techniques, CCT is rapidly evolving to offer information beyond the evaluation of cardiac structure and epicardial coronary arteries to provide details on myocardial deformation, extracellular volume, and coronary vasoreactivity. The potential for molecular imaging in CCT is also growing. In the current manuscript we review these emerging computed tomography techniques and their potential role in the evaluation of anthracycline-induced cardiotoxicity.

Computed Tomographic Angiography Assessment of Epicardial Coronary Vasoreactivity for Early Detection of Doxorubicin-Induced Cardiotoxicity.

BACKGROUND: The vascular endothelium is a novel target for the detection, management, and prevention of doxorubicin (DOX)-induced cardiotoxicity. OBJECTIVES: The study aimed to: 1) develop a methodology by computed tomography angiography (CTA) to evaluate stress-induced changes in epicardial coronary diameter; and 2) apply this to a chronic canine model of DOX-induced cardiotoxicity to assess vascular toxicity. METHODS: To develop and validate quantitative methods, sequential retrospectively gated coronary CTAs were performed in 16 canines. Coronary diameters were measured at prespecified distances during rest, adenosine (ADE) (280 µg/kg/min), rest 30 min post-ADE, and dobutamine (DOB) (5 µg/kg/min). A subgroup of 8 canines received weekly intravenous DOX (1 mg/kg) for 12 to 15 weeks, followed by rest-stress CTA at cumulative doses of ∼4-mg/kg (3 to 5 mg/kg), ∼8-mg/kg (7 to 9 mg/kg), and ∼12-mg/kg (12 to 15 mg/kg) of DOX. Echocardiograms were performed at these timepoints to assess left ventricular ejection fraction and global longitudinal strain. RESULTS: Under normal conditions, epicardial coronary arteries reproducibly dilated in response to ADE (left anterior descending coronary artery [LAD]: 12 ± 2%, left circumflex coronary artery [LCx]: 13 ± 2%, right coronary artery [RCA]: 14 ± 2%) and DOB (LAD: 17 ± 3%, LCx: 18 ± 2%, RCA: 15 ± 3%). With DOX, ADE vasodilator responses were impaired after ∼4-mg/kg (LAD: -3 ± 1%, LCx: 0 ± 2%, RCA: -5 ± 2%) and ∼8-mg/kg (LAD: -3 ± 1%, LCx: 0 ± 1%, RCA: -2 ± 2%). The DOB dilation response was preserved at ∼4-mg/kg of DOX (LAD: 18 ± 4%, LCx: 11 ± 3%, RCA: 11 ± 2%) but tended to decrease at ∼8-mg/kg of DOX (LAD: 4 ± 2%, LCx: 8 ± 3%, RCA: 3 ± 2%). A significant left ventricular ejection fraction reduction was observed only at 12 to 15 mg/kg DOX (baseline: 63 ± 2%, 12-mg/kg: 45 ± 3%). Global longitudinal strain was abnormal at ∼4-mg/kg of DOX (p = 0.011). CONCLUSIONS: CTA can reliably assess epicardial coronary diameter in response to pharmacological stressors, providing a noninvasive functional index of coronary vasoreactivity. Impaired epicardial vasodilation occurs early in DOX-induced cardiotoxicity.

Serial Assessment of Coronary Flow Reserve by Rubidium-82 Positron Emission Tomography Predicts Mortality in Heart Transplant Recipients.

OBJECTIVES: This study aimed to evaluate the long-term prognostic value of serial assessment of coronary flow reserve (CFR) by rubidium Rb 82 (82Rb) positron emission tomography (PET) in heart transplantation (HT) patients. BACKGROUND: Cardiac allograft vasculopathy is a major determinant of late mortality in HT recipients. The long-term prognostic value of serial CFR quantification by PET imaging in HT patients is unknown. METHODS: A total of 89 patients with history of HT (71% men, 7.0 ± 5.7 years post-HT, age 57 ± 11 years) scheduled for dynamic rest and stress (dipyridamole) 82Rb PET between March 1, 2008 and July 31, 2009 (PET-1) were prospectively enrolled in a single-center study. PET myocardial perfusion studies were reprocessed using U.S. Food and Drug Administration-approved software (Corridor 4DM, version 2017) for calculation of CFR. Follow-up PET (PET-2) imaging was performed in 69 patients at 1.9 ± 0.3 years following PET-1. Patients were categorized based on CFR values considering CFR ≤1.5 as low and CFR >1.5 as high CFR. RESULTS: Forty deaths occurred during the median follow-up time of 8.6 years. Low CFR at PET-1 was associated with a 2.77-fold increase in all-cause mortality (95% confidence interval [CI]: 1.34 to 5.74; p = 0.004). CFR decreased over time in patients with follow-up imaging (PET-1: 2.11 ± 0.74 vs. PET-2: 1.81 ± 0.61; p = 0.003). Twenty-five patients were reclassified based on PET-1 and PET-2 (high to low CFR: n = 18, low to high CFR: n = 7). Overall survival was similar in patients reclassified from high to low as patients with low to low CFR, whereas patients reclassified from low to high had similar survival as patients with high to high CFR. In multivariate Cox regression of patients with PET-2, higher baseline CFR (hazard ratio [HR] for a 0.73 unit (one SD) increase: 0.36, 95% CI: 0.16 to 0.82) and reduction in CFR from PET-1 to PET-2 (HR for a 0.79 unit (one SD) decrease: 1.50 to 7.84) were independent predictors of all-cause mortality. CONCLUSIONS: Serial assessment of CFR by 82Rb PET independently predicts long-term mortality in HT patients.

Quantitative Assessment of Coronary Microvascular Function: Dynamic Single-Photon Emission Computed Tomography, Positron Emission Tomography, Ultrasound, Computed Tomography, and Magnetic Resonance Imaging.

A healthy, functional microcirculation in combination with nonobstructed epicardial coronary arteries is the prerequisite of normal myocardial perfusion. Quantitative assessment in myocardial perfusion and determination of absolute myocardial blood flow can be achieved noninvasively using dynamic imaging with multiple imaging modalities. Extensive evidence supports the clinical value of noninvasively assessing indices of coronary flow for diagnosing coronary microvascular dysfunction; in certain diseases, the degree of coronary microvascular impairment carries important prognostic relevance. Although, currently positron emission tomography is the most commonly used tool for the quantification of myocardial blood flow, other modalities, including single-photon emission computed tomography, computed tomography, magnetic resonance imaging, and myocardial contrast echocardiography, have emerged as techniques with great promise for determination of coronary microvascular dysfunction. The following review will describe basic concepts of coronary and microvascular physiology, review available modalities for dynamic imaging for quantitative assessment of coronary perfusion and myocardial blood flow, and discuss their application in distinct forms of coronary microvascular dysfunction.

T1-refBlochi: high resolution 3D post-contrast T1 myocardial mapping based on a single 3D late gadolinium enhancement volume, Bloch equations, and a reference T1.

BACKGROUND: High resolution 3D T1 mapping is important for assessment of diffuse myocardial fibrosis in left atrium or other thin-walled structures. In this work, we investigated a fast single-TI 3D high resolution T1 mapping method that directly transforms a 3D late gadolinium enhancement (LGE) volume to a 3D T1 map. METHODS: The proposed method, T1-refBlochi, is based on Bloch equation modeling of the LGE signal, a single-point calibration, and assumptions that proton density and T2* are relatively uniform in the heart. Several sources of error of this method were analyzed mathematically and with simulations. Imaging was performed in phantoms, eight swine and five patients, comparing T1-refBlochi to a standard spin-echo T1 mapping, 3D multi-TI T1 mapping, and 2D ShMOLLI, respectively. RESULTS: The method has a good accuracy and adequate precision, even considering various sources of error. In phantoms, over a range of protocols, heart-rates and T1 s, the bias ±1SD was -3 ms ± 9 ms. The porcine studies showed excellent agreement between T1-refBlochi and the multi-TI method (bias ±1SD = -6 ± 22 ms). The proton density and T2* weightings yielded ratios for scar/blood of 0.94 ± 0.01 and for myocardium/blood of 1.03 ± 0.02 in the eight swine, confirming that sufficient uniformity of proton density and T2* weightings exists among heterogeneous tissues of the heart. In the patients, the mean T1 bias ±1SD in myocardium and blood between T1-refBlochi and ShMOLLI was -9 ms ± 21 ms. CONCLUSION: T1-refBlochi provides a fast single-TI high resolution 3D T1 map of the heart with good accuracy and adequate precision.

PET-CMR in heart failure - synergistic or redundant imaging?

Imaging in heart failure (HF) provides data for diagnosis, prognosis and disease monitoring. Both MRI and nuclear imaging techniques have been successfully used for this purpose in HF. Positron Emission Tomography-Cardiac Magnetic Resonance (PET-CMR) is an example of a new multimodality diagnostic imaging technique with potential applications in HF. The threshold for adopting a new diagnostic tool to clinical practice must necessarily be high, lest they exacerbate costs without improving care. New modalities must demonstrate clinical superiority, or at least equivalence, combined with another important advantage, such as lower cost or improved patient safety. The purpose of this review is to outline the current status of multimodality PET-CMR with regard to HF applications, and determine whether the clinical utility of this new technology justifies the cost.

High-Sensitivity and High-Resolution SPECT/CT Systems Provide Substantial Dose Reduction Without Compromising Quantitative Precision for Assessment of Myocardial Perfusion and Function.

UNLABELLED: There is increasing concern about radiation exposure from myocardial perfusion SPECT (MPS). γ-cameras with solid-state cadmium-zinc-telluride (CZT) detectors have better count sensitivity and spatial resolution than conventional sodium iodine detectors, allowing for significant reductions in radiotracer dose or acquisition time. This study aimed to demonstrate the capability of a hybrid CZT SPECT/64-slice CT system for dose reduction and to determine the maximal reduction possible without compromising image quality or the quantification precision of clinical MPS. METHODS: The imaging data of patients with normal myocardial perfusion and 30 patients with mid-sized to large perfusion defects who had undergone stress (99m)Tc-tetrofosmin MPS were postprocessed. Low-dose (361 ± 60 MBq) and high-dose (725 ± 142 MBq) (99m)Tc-tetrofosmin scans were included, with 6-min and 4-min scanning times, respectively. List-mode SPECT data were reconstructed with CT-based attenuation correction and with full as well as 50% and 75% reductions in acquisition time to simulate the corresponding relative dose reductions. The reconstructed SPECT images were analyzed to calculate global MPS defect size and regional defect size for 3 coronary artery territories-left anterior descending, left circumflex, and right-as well as left ventricular (LV) volume and ejection fraction. RESULTS: For patients with normal MPS results, there were no differences in defect size, LV volume, or ejection fraction, regardless of whether 50% or 75% reduction was used. For patients with abnormal MPS results, at a 50% reduction there was a significant difference in global defect size but not in regional defect size in the left anterior descending, left circumflex, and right coronary artery territories, whereas at a 75% reduction the difference was statistically significant in all territories, including the difference in global defect size. Nonetheless, differences in the defect size were minimal. The LV end-diastolic and end-systolic volumes and LV ejection fraction were not significantly different, regardless of whether 50% or 75% dose reduction was used. CONCLUSION: Ultra-low-dose (<190 MBq) MPS even with short imaging times (<6 min) is feasible using a hybrid CZT SPECT/CT camera without compromising image quality or significantly altering quantification of myocardial perfusion or LV function. We demonstrated that an additional 50% reduction in the current low-dose recommendations from the American Society of Nuclear Cardiology guidelines for (99m)Tc-labeled MPS is highly feasible while retaining short imaging protocols.

Scatter and crosstalk corrections for (99m)Tc/(123)I dual-radionuclide imaging using a CZT SPECT system with pinhole collimators.

PURPOSE: The energy spectrum for a cadmium zinc telluride (CZT) detector has a low energy tail due to incomplete charge collection and intercrystal scattering. Due to these solid-state detector effects, scatter would be overestimated if the conventional triple-energy window (TEW) method is used for scatter and crosstalk corrections in CZT-based imaging systems. The objective of this work is to develop a scatter and crosstalk correction method for (99m)Tc/(123)I dual-radionuclide imaging for a CZT-based dedicated cardiac SPECT system with pinhole collimators (GE Discovery NM 530c/570c). METHODS: A tailing model was developed to account for the low energy tail effects of the CZT detector. The parameters of the model were obtained using (99m)Tc and (123)I point source measurements. A scatter model was defined to characterize the relationship between down-scatter and self-scatter projections. The parameters for this model were obtained from Monte Carlo simulation using SIMIND. The tailing and scatter models were further incorporated into a projection count model, and the primary and self-scatter projections of each radionuclide were determined with a maximum likelihood expectation maximization (MLEM) iterative estimation approach. The extracted scatter and crosstalk projections were then incorporated into MLEM image reconstruction as an additive term in forward projection to obtain scatter- and crosstalk-corrected images. The proposed method was validated using Monte Carlo simulation, line source experiment, anthropomorphic torso phantom studies, and patient studies. The performance of the proposed method was also compared to that obtained with the conventional TEW method. RESULTS: Monte Carlo simulations and line source experiment demonstrated that the TEW method overestimated scatter while their proposed method provided more accurate scatter estimation by considering the low energy tail effect. In the phantom study, improved defect contrasts were observed with both correction methods compared to no correction, especially for the images of (99m)Tc in dual-radionuclide imaging where there is heavy contamination from (123)I. In this case, the nontransmural defect contrast was improved from 0.39 to 0.47 with the TEW method and to 0.51 with their proposed method and the transmural defect contrast was improved from 0.62 to 0.74 with the TEW method and to 0.73 with their proposed method. In the patient study, the proposed method provided higher myocardium-to-blood pool contrast than that of the TEW method. Similar to the phantom experiment, the improvement was the most substantial for the images of (99m)Tc in dual-radionuclide imaging. In this case, the myocardium-to-blood pool ratio was improved from 7.0 to 38.3 with the TEW method and to 63.6 with their proposed method. Compared to the TEW method, the proposed method also provided higher count levels in the reconstructed images in both phantom and patient studies, indicating reduced overestimation of scatter. Using the proposed method, consistent reconstruction results were obtained for both single-radionuclide data with scatter correction and dual-radionuclide data with scatter and crosstalk corrections, in both phantom and human studies. CONCLUSIONS: The authors demonstrate that the TEW method leads to overestimation in scatter and crosstalk for the CZT-based imaging system while the proposed scatter and crosstalk correction method can provide more accurate self-scatter and down-scatter estimations for quantitative single-radionuclide and dual-radionuclide imaging.

Multimodality imaging approach for serial assessment of regional changes in lower extremity arteriogenesis and tissue perfusion in a porcine model of peripheral arterial disease.

BACKGROUND: A standard quantitative imaging approach to evaluate peripheral arterial disease does not exist. Quantitative tools for evaluating arteriogenesis in vivo are not readily available, and the feasibility of monitoring serial regional changes in lower extremity perfusion has not been examined. METHODS AND RESULTS: Serial changes in lower extremity arteriogenesis and muscle perfusion were evaluated after femoral artery occlusion in a porcine model using single photon emission tomography (SPECT)/CT imaging with postmortem validation of in vivo findings using gamma counting, postmortem imaging, and histological analysis. Hybrid 201Tl SPECT/CT imaging was performed in pigs (n=8) at baseline, immediately postocclusion, and at 1 and 4 weeks postocclusion. CT imaging was used to identify muscle regions of interest in the ischemic and nonischemic hindlimbs for quantification of regional changes in CT-defined arteriogenesis and quantification of 201Tl perfusion. Four weeks postocclusion, postmortem tissue 201Tl activity was measured by gamma counting, and immunohistochemistry was performed to assess capillary density. Relative 201Tl retention (ischemic/nonischemic) was reduced immediately postocclusion in distal and proximal muscles and remained lower in calf and gluteus muscles 4 weeks later. Analysis of CT angiography revealed collateralization at 4 weeks within proximal muscles (P<0.05). SPECT perfusion correlated with tissue gamma counting at 4 weeks (P=0.01). Increased capillary density was seen within the ischemic calf at 4 weeks (P=0.004). CONCLUSIONS: 201Tl SPECT/CT imaging permits serial, regional quantification of arteriogenesis and resting tissue perfusion after limb ischemia. This approach may be effective for detection of disease and monitoring therapy in peripheral arterial disease.

Assessment of left ventricular 2D flow pathlines during early diastole using spatial modulation of magnetization with polarity alternating velocity encoding: a study in normal volunteers and canine animals with myocardial infarction.

A high-temporal resolution 2D flow pathline analysis method to study early diastolic filling is presented. Filling patterns in normal volunteers (n = 8) and canine animals [baseline (n = 1) and infarcted (n = 6)] are studied. Data are acquired using spatial modulation of magnetization with polarity alternating velocity encoding, which permits simultaneous quantification of 1D blood velocities (using phase contrast encoding) and myocardial strain (using spatial modulation of magnetization tagging and harmonic phase analysis) at high-temporal resolution of 14 ms within a single breath hold. Virtual emitter particles, released from the mitral valve plane every time frame during rapid filling, are tracked to depict the 2D pathlines on the imaged plane. The pathline regional distribution is compared with myocardial longitudinal strains and to regional pressure gradients. Quantitative analysis of net kinetic energy of pathlines is finally performed. Our results demonstrate a linear correlation (r(2) = 0.85) between pathline spatial distribution and myocardial strain. Peak net kinetic energy of 0.06 ± 0.01 mJ in normal volunteers, 0.043 mJ in baseline dog, 0.143 ± 0.03 mJ in infarcted dogs with nominal flow dysfunction, and 0.016 ± 0.007 mJ in infarcted dogs with severe flow dysfunction is observed. In conclusion, 2D pathline analysis provides a direct regional assessment of early diastolic filling patterns and is sensitive to abnormalities in early diastolic filling.

Assessment of early diastolic strain-velocity temporal relationships using spatial modulation of magnetization with polarity alternating velocity encoding (SPAMM-PAV).

A novel MR imaging technique, spatial modulation of magnetization with polarity alternating velocity encoding (SPAMM-PAV), is presented to simultaneously examine the left ventricular early diastolic temporal relationships between myocardial deformation and intra-cavity hemodynamics with a high temporal resolution of 14 ms. This approach is initially evaluated in a dynamic flow and tissue mimicking phantom. A comparison of regional longitudinal strains and intra-cavity pressure differences (integration of computed in-plane pressure gradients within a selected region) in relation to mitral valve inflow velocities is performed in eight normal volunteers. Our results demonstrate that apical regions have higher strain rates (0.145 ± 0.005 %/ms) during the acceleration period of rapid filling compared to mid-ventricular (0.114 ± 0.007 %/ms) and basal regions (0.088 ± 0.009 %/ms), and apical strain curves plateau at peak mitral inflow velocity. This pattern is reversed during the deceleration period, when the strain-rates in the basal regions are the highest (0.027 ± 0.003 %/ms) due to ongoing basal stretching. A positive base-to-apex gradient in peak pressure difference is observed during acceleration, followed by a negative base-to-apex gradient during deceleration. These studies shed insight into the regional volumetric and pressure difference changes in the left ventricle during early diastolic filling.

Vascular Tree Reconstruction by Minimizing A Physiological Functional Cost.

The reconstruction of complete vascular trees from medical images has many important applications. Although vessel detection has been extensively investigated, little work has been done on how connect the results to reconstruct the full trees. In this paper, we propose a novel theoretical framework for automatic vessel connection, where the automation is achieved by leveraging constraints from the physiological properties of the vascular trees. In particular, a physiological functional cost for the whole vascular tree is derived and an efficient algorithm is developed to minimize it. The method is generic and can be applied to different vessel detection/segmentation results, e.g. the classic rigid detection method as adopted in this paper. We demonstrate the effectiveness of this method on both 2D and 3D data.