Publications

Current guidelines recommend an implantable cardioverter-defibrillator (ICD) according to the left ventricular ejection fraction (LVEF). However, they do not mandate volumetric LVEF assessment. We sought to determine whether volumetric LVEF measurement using cardiovascular magnetic resonance imaging (CMR-LVEF) is superior to conventional LVEF measurement using 2-dimensional transthoracic echocardiography (Echo-LVEF) for risk stratifying patients referred for primary prevention ICD. Patients who underwent primary prevention ICD implantation at our institution and had undergone preimplantation CMR-LVEF from November 2001 to February 2011 were identified. Volumetric CMR-LVEF was determined from cine short-axis data sets. CMR-LVEF and Echo-LVEF were extracted from the clinical reports. The end point was appropriate ICD discharge (shock and/or antitachycardia pacing). Of 48 patients, appropriate ICD discharge occurred in 9 (19%) within 29 ± 25 months (range 1 to 99, median 20). All patients met the Echo-LVEF criteria for ICD implantation; however 25% (95% confidence interval 13% to 37%) did not meet the CMR-LVEF criteria. None (0%) of these latter patients had received an appropriate ICD discharge. Using CMR-LVEF ≤30% as a threshold for ICD eligibility, 19 patients (40%) with a qualifying Echo-LVEF would not have been referred for ICD, and none (0%) received an ICD discharge.For primary prevention ICD implantation, volumetric CMR-LVEF might be superior to clinical Echo-LVEF for risk stratification and can identify a large minority of subjects in whom ICD implantation can be safely avoided. In conclusion, if confirmed by larger prospective series, volumetric methods such as CMR should be considered a superior "gatekeeper" for the identification of patients likely to benefit from primary prevention ICD implantation.

Electroanatomical voltage mapping (EAVM) is commonly performed prior to catheter ablation of scar-related ventricular tachycardia (VT) to locate the arrhythmic substrate and to guide the ablation procedure. EAVM is used to locate the position of the ablation catheter and to provide a 3-D reconstruction of left-ventricular anatomy and scar. However, EAVM measurements only represent the endocardial scar with no transmural or epicardial information. Furthermore, EAVM is a time-consuming procedure, with a high operator dependence and has low sampling density, i.e., spatial resolution. Late gadolinium enhancement (LGE) magnetic resonance imaging (MRI) allows noninvasive assessment of scar morphology that can depict 3-D scar architecture. Despite the potential use of LGE as a roadmap for VT ablation for identification of arrhythmogenic substrate, its utility has been very limited. To allow for identification of VT substrate, a correlation is needed between the substrates identified by EAVM as the gold standard and LGE-MRI scar characteristics. To do so, a system must be developed to fuse the datasets from these modalities. In this study, a registration pipeline for the fusion of LGE-MRI and EAVM data is presented. A novel surface registration algorithm is proposed, integrating the matching of global scar areas as an additional constraint in the registration process. A preparatory landmark registration is initially performed to expedite the convergence of the algorithm. Numerical simulations were performed to evaluate the accuracy of the registration in the presence of errors in identifying landmarks in EAVM or LGE-MRI datasets as well as additional errors due to respiratory or cardiac motion. Subsequently, the accuracy of the proposed fusion system was evaluated in a cohort of ten patients undergoing VT ablation where both EAVM and LGE-MRI data were available. Compared to landmark registration and surface registration, the presented method achieved significant improvement in registration error. The proposed data fusion system allows the fusion of EAVM and LGE-MRI data in VT ablation with registration errors less than 3.5  mm.

PURPOSE: To develop and evaluate a technique that uses the k-space lines rejected by prospective respiratory navigator (NAV) to improve the signal-to-noise ratio (SNR) without increasing the scan time. METHODS: In conventional image reconstruction, the motion-corrupted k-space lines rejected by the NAV are not used. In this study, a set of translational motion parameters for the NAV-rejected lines and a phase-corrected average for the k-space line are estimated jointly using a maximum-likelihood approach and the information from the corresponding accepted k-space lines. Left coronary artery images were acquired in 10 healthy adult subjects, and the proposed approach incorporating the NAV-rejected lines was compared with the conventional dataset with NAV-accepted lines only, as well as a simple average of all k-space lines, in terms of SNR, normalized vessel sharpness and qualitative image scores on a four-point scale (1 = poor, 4 = excellent). Late gadolinium enhancement images of the left atrium were also acquired in 21 patients with atrial fibrillation pre- or post-pulmonary vein isolation. Images reconstructed with the proposed, conventional, and simple averaging methods were compared in terms of SNR, and subjective image quality on a four-point scale. RESULTS: For coronary MRI, there was a significant improvement in SNR with the proposed technique, but no significant difference in normalized vessel sharpness or qualitative image scores were observed with respect to the conventional method. Simple averaging resulted in an SNR gain, but significant loss in vessel sharpness and image quality. For late gadolinium enhancement, there was a significant increase in SNR, but no significant differences were observed in subjective image quality scores between the proposed and conventional methods. There was an SNR gain, but image quality loss for simple averaging, when compared with the conventional technique. In both coronary MRI and late gadolinium enhancement, the SNR gain of the proposed method was not significantly different than the maximum theoretical SNR gain. CONCLUSION: The proposed technique improves SNR using the additional information from NAV-rejected k-space lines, while providing similar image quality to standard reconstruction using motion-free k-space data only, with no increase in scan time.

Respiratory motion compensation using diaphragmatic navigator gating with a 5 mm gating window is conventionally used for free-breathing cardiac MRI. Because of the narrow gating window, scan efficiency is low resulting in long scan times, especially for patients with irregular breathing patterns. In this work, a new retrospective motion compensation algorithm is presented to reduce the scan time for free-breathing cardiac MRI that increasing the gating window to 15 mm without compromising image quality. The proposed algorithm iteratively corrects for respiratory-induced cardiac motion by optimizing the sharpness of the heart. To evaluate this technique, two coronary MRI datasets with 1.3 mm(3) resolution were acquired from 11 healthy subjects (seven females, 25 ± 9 years); one using a navigator with a 5 mm gating window acquired in 12.0 ± 2.0 min and one with a 15 mm gating window acquired in 7.1 ± 1.0 min. The images acquired with a 15 mm gating window were corrected using the proposed algorithm and compared to the uncorrected images acquired with the 5 and 15 mm gating windows. The image quality score, sharpness, and length of the three major coronary arteries were equivalent between the corrected images and the images acquired with a 5 mm gating window (P-value > 0.05), while the scan time was reduced by a factor of 1.7.

Phase contrast (PC) cardiac MR is widely used for the clinical assessment of blood flow in cardiovascular disease. One of the challenges of PC cardiac MR is the long scan time which limits both spatial and temporal resolution. Compressed sensing reconstruction with accelerated PC acquisitions is a promising technique to increase the scan efficiency. In this study, we sought to use the sparsity of the complex difference of the two flow-encoded images as an additional constraint term to improve the compressed sensing reconstruction of the corresponding accelerated PC data acquisition. Using retrospectively under-sampled data, the proposed reconstruction technique was optimized and validated in vivo on 15 healthy subjects. Then, prospectively under-sampled data was acquired on 11 healthy subjects and reconstructed with the proposed technique. The results show that there is good agreement between the cardiac output measurements from the fully sampled data and the proposed compressed sensing reconstruction method using complex difference sparsity up to acceleration rate 5. In conclusion, we have developed and evaluated an improved reconstruction technique for accelerated PC cardiac MR that uses the sparsity of the complex difference of the two flow-encoded images.

PURPOSE: To investigate the use of radiofrequency (RF) gating in conjunction with a paramagnetic contrast agent to reduce the specific absorption rate (SAR) and increase the blood-myocardium contrast in balanced steady-state free precession (bSSFP) 3D cardiac cine. MATERIALS AND METHODS: RF gating was implemented by synchronizing the RF-excitation with an external respiratory sensor (bellows), which could additionally be used for respiratory gating. For reference, respiratory-gated 3D cine images were acquired without RF gating. Free-breathing 3D cine images were acquired in eight healthy subjects before and after contrast injection (Gd-BOPTA) and compared to breath-hold 2D cine. RESULTS: RF-gated 3D cine reduced the SAR by nearly 40% without introducing significant artifacts while providing left ventricle (LV) measurements similar to those obtained with 2D cine. The contrast-to-noise ratio (CNR) was significantly higher for 3D cine compared to 2D cine, both before and after contrast injection; however, no statistically significant CNR increase was observed for the postcontrast 3D cine compared to the precontrast acquisitions. CONCLUSION: Respiratory-triggered RF gating significantly reduces SAR in 3D cine acquisitions, which may enable a more widespread clinical use of 3D cine. Furthermore, CNR of 3D bSSFP cine is higher than of 2D and administration of Gd-BOPTA does not improve the CNR of 3D cine.

PURPOSE: To develop an improved chemical shift-based water-fat separation sequence using a water-selective inversion pulse for inversion recovery 3D contrast-enhanced cardiac magnetic resonance imaging (MRI). MATERIALS AND METHODS: In inversion recovery sequences the fat signal is substantially reduced due to the application of a nonselective inversion pulse. Therefore, for simultaneous visualization of water, fat, and myocardial enhancement in inversion recovery-based sequences such as late gadolinium enhancement imaging, two separate scans are used. To overcome this, the nonselective inversion pulse is replaced with a water-selective inversion pulse. Imaging was performed in phantoms, nine healthy subjects, and nine patients with suspected arrhythmogenic right ventricular cardiomyopathy plus one patient for tumor/mass imaging. In patients, images with conventional turbo-spin echo (TSE) with and without fat saturation were acquired prior to contrast injection for fat assessment. Subjective image scores (1 = poor, 4 = excellent) were used for image assessment. RESULTS: Phantom experiments showed a fat signal-to-noise ratio (SNR) increase between 1.7 to 5.9 times for inversion times of 150 and 300 msec, respectively. The water-selective inversion pulse retains the fat signal in contrast-enhanced cardiac MR, allowing improved visualization of fat in the water-fat separated images of healthy subjects with a score of 3.7 ± 0.6. Patient images acquired with the proposed sequence were scored higher when compared with a TSE sequence (3.5 ± 0.7 vs. 2.2 ± 0.5, P < 0.05). CONCLUSION: The water-selective inversion pulse retains the fat signal in inversion recovery-based contrast-enhanced cardiac MR, allowing simultaneous visualization of water and fat.

Accurate fusion of late gadolinium enhancement magnetic resonance imaging (MRI) and electro-anatomical voltage mapping (EAM) is required to evaluate the potential of MRI to identify the substrate of ventricular tachycardia. However, both datasets are not acquired at the same cardiac phase and EAM data is corrupted with respiratory motion limiting the accuracy of current rigid fusion techniques. Knowledge of cardiac and respiratory motion during EAM is thus required to enhance the fusion process. In this study, we propose a novel approach to characterize both cardiac and respiratory motion from EAM data using the temporal evolution of the 3D catheter location recorded from clinical EAM systems. Cardiac and respiratory motion components are extracted from the recorded catheter location using multi-band filters. Filters are calibrated for each EAM point using estimates of heart rate and respiratory rate. The method was first evaluated in numerical simulations using 3D models of cardiac and respiratory motions of the heart generated from real time MRI data acquired in 5 healthy subjects. An accuracy of 0.6-0.7 mm was found for both cardiac and respiratory motion estimates in numerical simulations. Cardiac and respiratory motions were then characterized in 27 patients who underwent LV mapping for treatment of ventricular tachycardia. Mean maximum amplitude of cardiac and respiratory motion was 10.2±2.7 mm (min = 5.5, max = 16.9) and 8.8±2.3 mm (min = 4.3, max = 14.8), respectively. 3D Cardiac and respiratory motions could be estimated from the recorded catheter location and the method does not rely on additional imaging modality such as X-ray fluoroscopy and can be used in conventional electrophysiology laboratory setting.

BACKGROUND: Despite the survival benefit of implantable-cardioverter-defibrillators (ICDs), the vast majority of patients receiving an ICD for primary prevention do not receive ICD therapy. We sought to assess the role of heterogeneous scar area (HSA) identified by late gadolinium enhancement cardiovascular magnetic resonance (LGE-CMR) in predicting appropriate ICD therapy for primary prevention of sudden cardiac death (SCD). METHODS: From September 2003 to March 2011, all patients who underwent primary prevention ICD implantation and had a pre-implantation LGE-CMR were identified. Scar size was determined using thresholds of 4 and 6 standard deviations (SD) above remote normal myocardium; HSA was defined using 3 different criteria; as the region between 2 SD and 4 SD (HSA2-4SD), between 2SD and 6SD (HSA2-6SD), and between 4SD and 6SD (HSA4-6SD). The end-point was appropriate ICD therapy. RESULTS: Out of 40 total patients followed for 25 ± 24 months, 7 had appropriate ICD therapy. Scar size measured by different thresholds was similar in ICD therapy and non-ICD therapy groups (P = NS for all). However, HSA2-4SD and HSA4-6SD were significantly larger in the ICD therapy group (P = 0.001 and P = 0.03, respectively). In multivariable model HSA2-4SD was the only significant independent predictor of ICD therapy (HR = 1.08, 95%CI: 1.00-1.16, P = 0.04). Kaplan-Meier analysis showed that patients with greater HSA2-4SD had a lower survival free of appropriate ICD therapy (P = 0.026). CONCLUSIONS: In primary prevention ICD implantation, LGE-CMR HSA identifies patients with appropriate ICD therapy. If confirmed in larger series, HSA can be used for risk stratification in primary prevention of SCD.