Search

PURPOSE OF REVIEW: The purpose of this review is to clarify the risks associated with MRI exams for patients with cardiac implantable electronic devices (CIEDs) and to provide information regarding the MRI examination protocol for patients with CIEDs. RECENT FINDINGS: Several prospective studies evaluated the feasibility of MRI exams for patients with CIEDs and reported no adverse events. These studies suggest that by following a specific MRI examination protocol and monitoring both CIED parameters and the patient’s symptoms, an MRI exam can be performed by appropriately trained personnel with an acceptable benefit-to-risk ratio. SUMMARY: Both MR unsafe and MR conditional CIEDs are commercially available, but there are no MR safe CIEDs. The potential risks faced by patients with CIEDs during an MRI exam are always present and warrant careful monitoring. Three magnetic fields in the MRI scanner interact with the device in ways that can damage the CIED or harm the patient. Due to safety concerns and out of an abundance of caution, the majority of MRI exams for patients with CIEDs are currently denied. However, when following a specific MRI exam protocol, these risks can be mitigated.

Retrospective gated cardiac computed tomography (CCT) images can provide high contrast and resolution images of the heart throughout the cardiac cycle. Feature tracking in retrospective CCT images using the temporal sparse free-form deformations (TSFFDs) registration method has previously been optimised for the left ventricle (LV). However, there is limited work on optimising nonrigid registration methods for feature tracking in the left atria (LA). This paper systematically optimises the sparsity weight (SW) and bending energy (BE) as two hyperparameters of the TSFFD method to track the LA endocardium from end-diastole (ED) to end-systole (ES) using 10-frame retrospective gated CCT images. The effect of two different control point (CP) grid resolutions was also investigated. TSFFD optimisation was achieved using the average surface distance (ASD), directed Hausdorff distance (DHD) and Dice score between the registered and ground truth surface meshes and segmentations at ES. For baseline comparison, the configuration optimised for LV feature tracking gave errors across the cohort of 0.826 ± 0.172 mm ASD, 5.882 ± 1.524 mm DHD, and 0.912 ± 0.033 Dice score. Optimising the SW and BE hyperparameters improved the TSFFD performance in tracking LA features, with case specific optimisations giving errors across the cohort of 0.750 ± 0.144 mm ASD, 5.096 ± 1.246 mm DHD, and 0.919 ± 0.029 Dice score. Increasing the CP resolution and optimising the SW and BE further improved tracking performance, with case specific optimisation errors of 0.372 ± 0.051 mm ASD, 2.739 ± 0.843 mm DHD and 0.949 ± 0.018 Dice score across the cohort. We therefore show LA feature tracking using TSFFDs is improved through a chamber-specific optimised configuration.

Despite positive initial outcomes emerging from preclinical and early clinical investigation of alginate hydrogel injection therapy as a treatment for heart failure, the lack of knowledge about the mechanism of action remains a major shortcoming that limits the efficacy of treatment design. To identify the mechanism of action, we examined previously unobtainable measurements of cardiac function from in vivo, ex vivo, and in silico states of clinically relevant heart failure (HF) in large animals. High-resolution ex vivo magnetic resonance imaging and histological data were used along with state-of-the-art subject-specific computational model simulations. Ex vivo data were incorporated in detailed geometric computational models for swine hearts in health (n=5), ischemic HF (n=5), and ischemic HF treated with alginate hydrogel injection therapy (n=5). Hydrogel injection therapy mitigated elongation of sarcomere lengths (1.68±0.10μm [treated] vs. 1.78±0.15μm [untreated], p99% accuracy and predicted small myofiber strain in the vicinity of the solidified hydrogel that was sustained for up to 13 mm away from the implant. These findings suggest that the solidified alginate hydrogel material acts as an LV mid-wall constraint that significantly reduces adverse LV remodeling compared to untreated HF controls without causing negative secondary outcomes to cardiac function. STATEMENT OF SIGNIFICANCE: Heart failure is considered a growing epidemic and hence an important health problem in the US and worldwide. Its high prevalence (5.8 million and 23 million, respectively) is expected to increase by 25% in the US alone by 2030. Heart failure is associated with high morbidity and mortality, has a 5-year mortality rate of 50%, and contributes considerably to the overall cost of health care ($53.1 billion in the US by 2030). Despite positive initial outcomes emerging from preclinical and early clinical investigation of alginate hydrogel injection therapy as a treatment for heart failure, the lack of knowledge concerning the mechanism of action remains a major shortcoming that limits the efficacy of treatment design. To understand the mechanism of action, we combined high-resolution ex vivo magnetic resonance imaging and histological data in swine with state-of-the-art subject-specific computational model simulations. The in silico models realistically simulated in vivo function with >99% accuracy and predicted small myofiber strain in the vicinity of the solidified hydrogel that was sustained for up to 13 mm away from the implant. These findings suggest that the solidified alginate hydrogel material acts as a left ventricular mid-wall constraint that significantly reduces adverse LV remodeling compared to untreated heart failure controls without causing negative secondary outcomes to cardiac function. Moreover, if the hydrogel can be delivered percutaneously rather than via the currently used open-chest procedure, this therapy may become routine for heart failure treatment. A minimally invasive procedure would be in the best interest of this patient population; i.e., one that cannot tolerate general anesthesia and surgery, and it would be significantly more cost-effective than surgery.

Aortic wall stiffening is a predictive marker for morbidity in hypertensive patients. Arterial pulse wave velocity (PWV) correlates with the level of stiffness and can be derived using non-invasive 4D-flow magnetic resonance imaging (MRI). The objectives of this study were twofold: to develop subject-specific thoracic aorta models embedded into an MRI-compatible flow circuit operating under controlled physiological conditions; and to evaluate how a range of aortic wall stiffness impacts 4D-flow-based quantification of hemodynamics, particularly PWV. Three aorta models were 3D-printed using a novel photopolymer material at two compliant and one nearly rigid stiffnesses and characterized via tensile testing. Luminal pressure and 4D-flow MRI data were acquired for each model and cross-sectional net flow, peak velocities, and PWV were measured. In addition, the confounding effect of temporal resolution on all metrics was evaluated. Stiffer models resulted in increased systolic pressures (112, 116, and 133 mmHg), variations in velocity patterns, and increased peak velocities, peak flow rate, and PWV (5.8–7.3 m/s). Lower temporal resolution (20 ms down to 62.5 ms per image frame) impacted estimates of peak velocity and PWV (7.31 down to 4.77 m/s). Using compliant aorta models is essential to produce realistic flow dynamics and conditions that recapitulated in vivo hemodynamics.

The analysis of quantitative hemodynamics and luminal pressure may add valuable information to aid treatment strategies and prognosis for aortic dissections. This work directly compared in vitro 4D-flow magnetic resonance imaging (MRI), catheter-based pressure measurements, and computational fluid dynamics that integrated fluid-structure interaction (CFD FSI). Experimental data was acquired with a compliant 3D-printed model of a type-B aortic dissection (TBAD) that was embedded into a flow circuit with tunable boundary conditions. In vitro flow and relative pressure information were used to tune the CFD FSI Windkessel boundary conditions. Results showed overall agreement of complex flow patterns, true to false lumen flow splits, and pressure distribution. This work demonstrates feasibility of a tunable experimental setup that integrates a patient-specific compliant model and provides a test bed for exploring critical imaging and modeling parameters that ultimately may improve the prognosis for patients with aortic dissections.

Background Cardiovascular disease is the leading cause of death in patients with Duchenne muscular dystrophy (DMD)—a fatal X-linked genetic disorder. Late gadolinium enhancement (LGE) imaging is the current gold standard for detecting myocardial tissue remodeling, but it is often a late finding. Current research aims to investigate cardiovascular magnetic resonance (CMR) biomarkers, including native (pre-contrast) T1 and extracellular volume (ECV) to evaluate the early on-set of microstructural remodeling and to grade disease severity. To date, native T1 measurements in DMD have been reported predominantly at 1.5T. This study uses 3T CMR: (1) to characterize global and regional myocardial pre-contrast T1 differences between healthy controls and LGE + and LGE− boys with DMD; and (2) to report global and regional myocardial post-contrast T1 values and myocardial ECV estimates in boys with DMD, and (3) to identify left ventricular (LV) T1-mapping biomarkers capable of distinguishing between healthy controls and boys with DMD and detecting LGE status in DMD. Methods Boys with DMD (N = 28, 13.2 ± 3.1 years) and healthy age-matched boys (N = 20, 13.4 ± 3.1 years) were prospectively enrolled and underwent a 3T CMR exam including standard functional imaging and T1 mapping using a modified Look-Locker inversion recovery (MOLLI) sequence. Pre-contrast T1 mapping was performed on all boys, but contrast was administered only to boys with DMD for post-contrast T1 and ECV mapping. Global and segmental myocardial regions of interest were contoured on mid LV T1 and ECV maps. ROI measurements were compared for pre-contrast myocardial T1 between boys with DMD and healthy controls, and for post-contrast myocardial T1 and ECV between LGE + and LGE− boys with DMD using a Wilcoxon rank-sum test. Results are reported as median and interquartile range (IQR). p-Values 1 mapping and LV function parameters in the tasks of distinguishing between healthy controls and boys with DMD, and detecting LGE status in DMD. The area under the curve is reported. Results Boys with DMD had significantly increased global native T1 [1332 (60) ms vs. 1289 (56) ms; p = 0.004] and increased within-slice standard deviation (SD) [100 (57) ms vs. 74 (27) ms; p = 0.001] compared to healthy controls. LGE− boys with DMD also demonstrated significantly increased lateral wall native T1 [1322 (68) ms vs. 1277 (58) ms; p = 0.001] compared to healthy controls. LGE + boys with DMD had decreased global myocardial post-contrast T1 [565 (113) ms vs 635 (126) ms; p = 0.04] and increased global myocardial ECV [32 (8) % vs. 28 (4) %; p = 0.02] compared to LGE− boys. In all classification tasks, T1-mapping biomarkers outperformed a conventional biomarker, LV ejection fraction. ECV was the best performing biomarker in the task of predicting LGE status (AUC = 0.95). Conclusions Boys with DMD exhibit elevated native T1 compared to healthy, sex- and age-matched controls, even in the absence of LGE. Post-contrast T1 and ECV estimates from 3T CMR are also reported here for pediatric patients with DMD for the first time and can distinguish between LGE + from LGE− boys. In all classification tasks, T1-mapping biomarkers outperform a conventional biomarker, LVEF.

Purpose Cardiomyocyte organization and performance underlie cardiac function, but the in vivo mobility of these cells during contraction and filling remains difficult to probe. Herein, a novel trigger delay (TD) scout sequence was used to acquire high in-plane resolution (1.6 mm) Spin-Echo (SE) cardiac diffusion tensor imaging (cDTI) at three distinct cardiac phases. The objective was to characterize cardiomyocyte organization and mobility throughout the cardiac cycle in healthy volunteers. Materials and methods Nine healthy volunteers were imaged with cDTI at three distinct cardiac phases (early systole, late systole, and diastasis). The sequence used a free-breathing Spin-Echo (SE) cDTI protocol (b-values = 350s/mm2, twelve diffusion encoding directions, eight repetitions) to acquire high-resolution images (1.6x1.6x8mm3) at 3T in ~7 minutes/cardiac phase. Helix Angle (HA), Helix Angle Range (HAR), E2 angle (E2A), Transverse Angle (TA), Mean Diffusivity (MD), diffusion tensor eigenvalues (λ1-2-3), and Fractional Anisotropy (FA) in the left ventricle (LV) were characterized. Results Images from the patient-specific TD scout sequence demonstrated that SE cDTI acquisition was possible at early systole, late systole, and diastasis in 78%, 100% and 67% of the cases, respectively. At the mid-ventricular level, mobility (reported as median [IQR]) was observed in HAR between early systole and late systole (76.9 [72.6, 80.5]° vs 96.6 [85.9, 100.3]°, p Conclusion We demonstrate that it is possible to probe cardiomyocyte mobility using multi-phase and high resolution cDTI. In healthy volunteers, aggregate cardiomyocytes re-orient themselves more longitudinally during contraction, while cardiomyocyte sheetlets tilt radially during wall thickening. These observations provide new insights into the three-dimensional mobility of myocardial microstructure during systolic contraction.

Author(s): Aliotta, Eric; Moulin, Kevin; Zhang, Zhaohuan; Ennis, Daniel B | Abstract: PURPOSE:To evaluate a technique for simultaneous quantitative T2 and apparent diffusion coefficient (ADC) mapping in the heart (T2 +ADC) using spin echo (SE) diffusion-weighted imaging (DWI). THEORY AND METHODS:T2 maps from T2 +ADC were compared with single-echo SE in phantoms and with T2 -prepared (T2 -prep) balanced steady-state free precession (bSSFP) in healthy volunteers. ADC maps from T2 +ADC were compared with conventional DWI in phantoms and in vivo. T2 +ADC was also demonstrated in a patient with acute myocardial infarction (MI). RESULTS:Phantom T2 values from T2 +ADC were closer to a single-echo SE reference than T2 -prep bSSFP (-2.3 ± 6.0% vs 22.2 ± 16.3%; P l 0.01), and ADC values were in excellent agreement with DWI (0.28 ± 0.4%). In volunteers, myocardial T2 values from T2 +ADC were significantly shorter than T2 -prep bSSFP (35.8 ± 3.1 vs 46.8 ± 3.8 ms; P l 0.01); myocardial ADC was not significantly (N.S.) different between T2 +ADC and conventional motion-compensated DWI (1.39 ± 0.18 vs 1.38 ± 0.18 mm2 /ms; P = N.S.). In the patient, T2 and ADC were both significantly elevated in the infarct compared with remote myocardium (T2 : 40.4 ± 7.6 vs 56.8 ± 22.0; P l 0.01; ADC: 1.47 ± 0.59 vs 1.65 ± 0.65 mm2 /ms; P l 0.01). CONCLUSION:T2 +ADC generated coregistered, free-breathing T2 and ADC maps in healthy volunteers and a patient with acute MI with no cost in accuracy, precision, or scan time compared with DWI. Magn Reson Med 79:654-662, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Purpose To design and evaluate eddy current–nulled convex optimized diffusion encoding (EN-CODE) gradient waveforms for efficient diffusion tensor imaging (DTI) that is free of eddy current–induced image distortions. Methods The EN-CODE framework was used to generate diffusion-encoding waveforms that are eddy current–compensated. The EN-CODE DTI waveform was compared with the existing eddy current–nulled twice refocused spin echo (TRSE) sequence as well as monopolar (MONO) and non–eddy current–compensated CODE in terms of echo time (TE) and image distortions. Comparisons were made in simulations, phantom experiments, and neuro imaging in 10 healthy volunteers. Results The EN-CODE sequence achieved eddy current compensation with a significantly shorter TE than TRSE (78 versus 96 ms) and a slightly shorter TE than MONO (78 versus 80 ms). Intravoxel signal variance was lower in phantoms with EN-CODE than with MONO (13.6 ± 11.6 versus 37.4 ± 25.8) and not different from TRSE (15.1 ± 11.6), indicating good robustness to eddy current–induced image distortions. Mean fractional anisotropy values in brain edges were also significantly lower with EN-CODE than with MONO (0.16 ± 0.01 versus 0.24 ± 0.02, P

BACKGROUND Cardiomyocytes are organized in microstructures termed sheetlets that reorientate during left ventricular thickening. Diffusion tensor cardiac magnetic resonance (DT-CMR) may enable noninvasive interrogation of in vivo cardiac microstructural dynamics. Dilated cardiomyopathy (DCM) is a condition of abnormal myocardium with unknown sheetlet function. OBJECTIVES This study sought to validate in vivo DT-CMR measures of cardiac microstructure against histology, characterize microstructural dynamics during left ventricular wall thickening, and apply the technique in hypertrophic cardiomyopathy (HCM) and DCM. METHODS In vivo DT-CMR was acquired throughout the cardiac cycle in healthy swine, followed by in situ and ex vivo DT-CMR, then validated against histology. In vivo DT-CMR was performed in 19 control subjects, 19 DCM, and 13 HCM patients. RESULTS In swine, a DT-CMR index of sheetlet reorientation (E2A) changed substantially (E2A mobility ~46°). E2A changes correlated with wall thickness changes (in vivo r2 = 0.75; in situ r2 = 0.89), were consistently observed under all experimental conditions, and accorded closely with histological analyses in both relaxed and contracted states. The potential contribution of cyclical strain effects to in vivo E2A was ~17%. In healthy human control subjects, E2A increased from diastole (18°) to systole (65°; p < 0.001; E2A mobility = 45°). HCM patients showed significantly greater E2A in diastole than control subjects did (48 ; p < 0.001) with impaired E2A mobility (23°; p < 0.001). In DCM, E2A was similar to control subjects in diastole, but systolic values were markedly lower (40° ; p < 0.001) with impaired E2A mobility (20°; p < 0.001). CONCLUSIONS Myocardial microstructure dynamics can be characterized by in vivo DT-CMR. Sheetlet function was abnormal in DCM with altered systolic conformation and reduced mobility, contrasting with HCM, which showed reduced mobility with altered diastolic conformation. These novel insights significantly improve understanding of contractile dysfunction at a level of noninvasive interrogation not previously available in humans. (J Am Coll Cardiol 2017;69:661–76) Published by Elsevier on behalf of the American College of Cardiology Foundation.

Purpose To characterize the effect of hepatic vessel flow using 4-dimensional (4D) flow magnetic resonance (MR) imaging and correlate their effect on microwave ablation volumes in an in vivo non-cirrhotic porcine liver model. Materials and Methods Microwave ablation antennas were placed under ultrasound guidance in each liver lobe of swine (n = 3 in each animal) for a total of 9 ablations. Pre- and post-ablation 4D flow MR imaging was acquired to quantify flow changes in the hepatic vasculature. Flow measurements, along with encompassed vessel size and vessel-antenna spacing, were then correlated with final ablation volume from segmented MR images. Results The linear regression model demonstrated that the preablation measurement of encompassed hepatic vein size (β = –0.80 ± 0.25, 95% confidence interval [CI] –1.15 to –0.22; P = .02) was significantly correlated to final ablation zone volume. The addition of hepatic vein flow rate found via 4D flow MRI (β = –0.83 ± 0.65, 95% CI –2.50 to 0.84; P = .26), and distance from antenna to hepatic vein (β = 0.26 ± 0.26, 95% CI –0.40 to 0.92; P = .36) improved the model accuracy but not significantly so (multivariate adjusted R2 = 0.70 vs univariate (vessel size) adjusted R2 = 0.63, P = .24). Conclusions Hepatic vein size in an encompassed ablation zone was found to be significantly correlated with final ablation zone volume. Although the univariate 4D flow MR imaging-acquired measurements alone were not found to be statistically significant, its addition to hepatic vein size improved the accuracy of the ablation volume regression model. Pre-ablation 4D flow MR imaging of the liver may assist in prospectively optimizing thermal ablation treatment.

Computational modeling of the heart requires accurately incorporating both gross anatomical detail and local microstructural information. Together, these provide the necessary data to build 3D meshes for simulation of cardiac mechanics and electrophysiology. Recent MRI advances make it possible to measure detailed heart motion in vivo, but in vivo microstructural imaging of the heart remains challenging. Consequently, the most detailed measurements of microstructural organization and microanatomical infarct details are obtained ex vivo. The objective of this work was to develop and evaluate a new method for restoring ex vivo ventricular geometry to match the in vivo configuration. This approach aids the integration of high-resolution ex vivo microstructural information with in vivo motion measurements. The method uses in vivo cine imaging to generate surface meshes, then creates a 3D printed left ventricular (LV) blood pool cast and a pericardial mold to restore the ex vivo cardiac geometry to a mid-diastasis reference configuration. The method was evaluated in healthy (N = 7) and infarcted (N = 3) swine. Dice similarity coefficients were calculated between in vivo and ex vivo images for the LV cavity (\(0.93\pm 0.01\)), right ventricle (RV) cavity (\(0.80\pm 0.05\)), and the myocardium (\(0.72\pm 0.04\)). The \(R^2\) coefficient between in vivo and ex vivo LV and RV cavity volumes were 0.95 and 0.91, respectively. These results suggest that this method adequately restores ex vivo geometry to match in vivo geometry. This approach permits a more precise incorporation of high-resolution ex vivo anatomical and microstructural data into computational models that use in vivo data for simulation of cardiac mechanics and electrophysiology.

Studies analyzing optimal voltage thresholds for scar detection with electroanatomic mapping frequently lack a gold standard for comparison.The purpose of this study was to use a porcine infarct model with ex vivo magnetic resonance imaging (MRI) integration to characterize the relationship between interelectrode spacing and bipolar voltage thresholds and examine the influence of 3-dimensional scar on unipolar voltages.Thirty-two combined endocardial-epicardial electroanatomic maps were created in 8 postinfarct porcine subjects (bipolar 2-mm, 5-mm, and 8-mm interelectrode spacing and unipolar) for comparison with ex vivo MRI. Two thresholds were compared: (1) 95% normal distribution and (2) best fit to MRI. Direct electrogram analysis was performed in regions across from MRI-defined scar and adjacent to scar border zone.A linear increase in optimal thresholds was observed with wider bipole spacing. The 95% thresholds for scar were lower than MRI-matched thresholds with moderate sensitivity for nontransmural scar (54% endo, 63% epi). Unipolar endocardial scar area exceeded MRI-defined scar, resulting in mismatched false scar in 5 of 8 (63%). Endocardial and epicardial unipolar voltages were lower than normal in regions adjacent and across from scar.Variations in interelectrode spacing necessitate tailored bipolar voltage thresholds to optimize scar detection. Statistical 95% thresholds appear to be conservative and not fully sensitive for the detection of scar defined by high-resolution ex vivo MRI. In the presence of endocardial scar, unipolar mapping to quantitatively characterize epicardial scar may be overly sensitive due to 3-dimensional spatial averaging.

Author(s): Aliotta, Eric; Wu, Holden H; Ennis, Daniel B | Abstract: PurposeTo evaluate convex optimized diffusion encoding (CODE) gradient waveforms for minimum echo time and bulk motion-compensated diffusion-weighted imaging (DWI).MethodsDiffusion-encoding gradient waveforms were designed for a range of b-values and spatial resolutions with and without motion compensation using the CODE framework. CODE, first moment (M1 ) nulled CODE-M1 , and first and second moment (M2 ) nulled CODE-M1 M2 were used to acquire neuro, liver, and cardiac ADC maps in healthy subjects (n=10) that were compared respectively to monopolar (MONO), BIPOLAR (M1 = 0), and motion-compensated (MOCO, M1 + M2 = 0) diffusion encoding.ResultsCODE significantly improved the SNR of neuro ADC maps compared with MONO (19.5 ± 2.5 versus 14.5 ± 1.9). CODE-M1 liver ADCs were significantly lower (1.3 ± 0.1 versus 1.8 ± 0.3 × 10-3 mm2 /s, ie, less motion corrupted) and more spatially uniform (6% versus 55% ROI difference) than MONO and had higher SNR than BIPOLAR (SNR = 14.9 ± 5.3 versus 8.0 ± 3.1). CODE-M1 M2 cardiac ADCs were significantly lower than MONO (1.9 ± 0.6 versus 3.8 ± 0.3 x10-3 mm2 /s) throughout the cardiac cycle and had higher SNR than MOCO at systole (9.1 ± 3.9 versus 7.0 ± 2.6) while reporting similar ADCs (1.5 ± 0.2 versus 1.4 ± 0.6 × 10-3 mm2 /s).ConclusionsCODE significantly improved SNR for ADC mapping in the brain, liver and heart, and significantly improved DWI bulk motion robustness in the liver and heart. Magn Reson Med 77:717-729, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Author(s): Li, Xinzhou; Perotti, Luigi E; Martinez, Jessica A; Duarte-Vogel, Sandra M; Ennis, Daniel B; Wu, Holden H | Abstract: PURPOSE:Real-time magnetic resonance imaging (MRI) is a promising alternative to X-ray fluoroscopy for guiding cardiovascular catheterization procedures. Major challenges, however, include the lack of guidewires that are compatible with the MRI environment, not susceptible to radiofrequency-induced heating, and reliably visualized. Preclinical evaluation of new guidewire designs has been conducted at 1.5T. Here we further evaluate the safety (device heating), device visualization, and procedural feasibility of 3T MRI-guided cardiovascular catheterization using a novel MRI-visible glass-fiber epoxy-based guidewire in phantoms and porcine models. METHODS:To evaluate device safety, guidewire tip heating (GTH) was measured in phantom experiments with different combinations of catheters and guidewires. In vivo cardiovascular catheterization procedures were performed in both healthy (N = 5) and infarcted (N = 5) porcine models under real-time 3T MRI guidance using a glass-fiber epoxy-based guidewire. The times for each procedural step were recorded separately. Guidewire visualization was assessed by measuring the dimensions of the guidewire-induced signal void and contrast-to-noise ratio (CNR) between the guidewire tip signal void and the blood signal in real-time gradient-echo MRI (specific absorption rate [SAR] = 0.04 W/kg). RESULTS:In the phantom experiments, GTH did not exceed 0.35°C when using the real-time gradient-echo sequence (SAR = 0.04 W/kg), demonstrating the safety of the glass-fiber epoxy-based guidewire at 3T. The catheter was successfully placed in the left ventricle (LV) under real-time MRI for all five healthy subjects and three out of five infarcted subjects. Signal void dimensions and CNR values showed consistent visualization of the glass-fiber epoxy-based guidewire in real-time MRI. The average time (minutes:seconds) for the catheterization procedure in all subjects was 4:32, although the procedure time varied depending on the subject's specific anatomy (standard deviation = 4:41). CONCLUSIONS:Real-time 3T MRI-guided cardiovascular catheterization using a new MRI-visible glass-fiber epoxy-based guidewire is feasible in terms of visualization and guidewire navigation, and safe in terms of radiofrequency-induced guidewire tip heating.

Introduction: Computational models of the heart increasingly require detailed microstructural information to capture the impact of tissue remodeling on cardiac electromechanics in, for example, hearts with myocardial infarctions. Myocardial infarctions are surrounded by the infarct border zone (BZ), which is a site of electromechanical property transition. Magnetic resonance imaging (MRI) is an emerging method for characterizing microstructural remodeling and focal myocardial infarcts and the BZ can be identified with late gadolinium enhanced (LGE) MRI. Microstructural remodeling within the BZ, however, remains poorly characterized by MRI due, in part, to the fact that LGE and DT-MRI are not always available for the same heart. Diffusion tensor MRI (DT-MRI) can evaluate microstructural remodeling by quantifying the DT apparent diffusion coefficient (ADC, increased with decreased cellularity), fractional anisotropy (FA, decreased with increased fibrosis), and tissue mode (decreased with increased fiber disarray). The purpose of this work was to use LGE MRI in post-infarct porcine hearts (N = 7) to segment remote, BZ, and infarcted myocardium, thereby providing a basis to quantify microstructural remodeling in the BZ and infarcted regions using co-registered DT-MRI. Methods: Chronic porcine infarcts were created by balloon occlusion of the LCx. 6-8 weeks post-infarction, MRI contrast was administered, and the heart was potassium arrested, excised, and imaged with LGE MRI (0.33 × 0.33 × 0.33 mm) and co-registered DT-MRI (1 × 1 × 3 mm). Myocardium was segmented as remote, BZ, or infarct by LGE signal intensity thresholds. DT invariants were used to evaluate microstructural remodeling by quantifying ADC, FA, and tissue mode. Results: The BZ significantly remodeled compared to both infarct and remote myocardium. BZ demonstrated a significant decrease in cellularity (increased ADC), significant decrease in tissue organization (decreased FA), and a significant increase in fiber disarray (decreased tissue mode) relative to remote myocardium (all p < 0.05). Microstructural remodeling in the infarct was similar, but significantly larger in magnitude (all p < 0.05). Conclusion: DT-MRI can identify regions of significant microstructural remodeling in the BZ that are distinct from both remote and infarcted myocardium.

Purpose The growth or rupture of cerebral aneurysms is thought to be a complex arterial wall response caused by an interaction of biological and hemodynamic factors. Current approaches that analyze the effects of flow stress on vascular cells are limited by their simplistic design, lacking patient-specific vascular anatomic features that contribute to disease. In this study, we compare wall shear stress (WSS) acquired using 4D Flow MRA from a 3D printed model with patient-specific aneurysm to gene expression of endothelial cells attached inside the same vascular model. Materials and methods Patient-specific vascular models were created from follow-up images in patients whose saccular aneurysms were confirmed to grow afterwards by CT angiography. The vascular models were connected to a flow-pump and a 4D Flow MRA was performed. WSS data was calculated from the acquisition and post-processed using a custom MATLAB environment. The same vascular models were coated with fibronectin and rotated in 3D with endothelium in culture incubator for cell lining. Viscosity-adjusted culture media was perfused using the same perfusion conditions with 4D Flow MRA. After perfusion, RNA was extracted from endothelial cells in the parent artery and the aneurysm and gene expression was examined using quantitative PCR (qPCR). In addition, the morphology of endothelial cells in specific regions of the model were observed with confocal microscopy. Results Endothelial cells in growing regions with low wall shear within the aneurysm were irregular in shape and size and associated with up-regulation of inflammatory genes compared with cells in the parent artery that exhibited typical spindle shape and aligned with the directionality of flow. Conclusions This research model enables a new approach to bridge the gap of biophysical flow phenomonology and the biological impact of complex flow patterns on endothelial cells using patient-derived imaging data. This approach is likely to be useful to determine flow-driven biological changes in vascular endothelium that contribute to aneurysm growth and rupture. Disclosures N. Kaneko: None. S. Tateshima: None. M. Loecher: None. J. Villablanca: None. D. Chen: None. Y. Komuro: None. W. Chang: None. D. Ennis: None. F. Vinuela: None. G. Duckwiler: None. J. Hinman: None.

Author(s): Ferreira, Pedro F; Nielles-Vallespin, Sonia; Scott, Andrew D; de Silva, Ranil; Kilner, Philip J; Ennis, Daniel B; Auger, Daniel A; Suever, Jonathan D; Zhong, Xiaodong; Spottiswoode, Bruce S; Pennell, Dudley J; Arai, Andrew E; Firmin, David N | Abstract: PurposeTo evaluate the importance of strain-correcting stimulated echo acquisition mode echo-planar imaging cardiac diffusion tensor imaging.MethodsHealthy pigs (n = 11) were successfully scanned with a 3D cine displacement-encoded imaging with stimulated echoes and a monopolar-stimulated echo-planar imaging diffusion tensor imaging sequence at 3 T during diastasis, peak systole, and strain sweet spots in a midventricular short-axis slice. The same diffusion tensor imaging sequence was repeated ex vivo after arresting the hearts in either a relaxed (KCl-induced) or contracted (BaCl2 -induced) state. The displacement-encoded imaging with stimulated echoes data were used to strain-correct the in vivo cardiac diffusion tensor imaging in diastole and systole. The orientation of the primary (helix angles) and secondary (E2A) diffusion eigenvectors was compared with and without strain correction and to the strain-free ex vivo data.ResultsStrain correction reduces systolic E2A significantly when compared without strain correction and ex vivo (median absolute E2A = 34.3° versus E2A = 57.1° (P = 0.01), E2A = 60.5° (P = 0.006), respectively). The systolic distribution of E2A without strain correction is closer to the contracted ex vivo distribution than with strain correction, root mean square deviation of 0.027 versus 0.038.ConclusionsThe current strain-correction model amplifies the contribution of microscopic strain to diffusion resulting in an overcorrection of E2A. Results show that a new model that considers cellular rearrangement is required. Magn Reson Med 79:2205-2215, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

In vivo visualization and quantification of surface tissue edema, or ‘tissue swelling’ following injury or as diagnostic criteria for the detection of disease remains a clinical challenge. The most common method for evaluating edema is by clinical examination, during which tissue swelling and turgor are routinely examined. This technique, however, is qualitative, difficult to standardize, and requires considerable clinical expertise. T2 — weighted (T2w) magnetic resonance imaging (MRI) has been used to visualize the hydration state of skin in vivo and provides comprehensive, layer-specific anatomical and physiological data. While clinical utility has been demonstrated in a range of applications, adoption has been limited due to instrument size, cost, and the technical challenges associated with surface tissue data interpretation. These factors provide a potential opportunity for THz medical imaging. In this talk we propose THz maps of tissue water content (TWC) as an adjunct to MRI. For applications where (1) surface tissue water content is correlated with disease/injury state and (2) technical challenges and/or adoption costs restrict the use of MRI, THz imaging can be leveraged to fill the diagnostic gap. Candidate applications are discussed including burn wound severity assessment, surgical flap viability prediction, and allergy testing.

Purpose: In this abstract, we discuss a biomechanical head and neck model that will be able to represent patient setup variations as well as physiologic changes and subsequently enable dose calculations on the deformed anatomy.Methods: We selected Multi Pose MRI as the imaging modality to aid in model development and validation. The MRI data allowed us to build a biomechanically predictive model that will enable accurate estimation of tumor position when seeded with CT data alone. The soft tissuecontrast and lack of ionizing radiation when using MRI enabled us to acquire extensive imaging datasets with a suitable variety of head pose variations. These poses were selected to encompass the clinical positioning variations so that the resulting model will accurately reflect internal organ motion and deformation. All images were acquired using an 8‐channel, 1.5T research MRI system in radiology. The imaging volume extended from about T3(upper thoracic vertebrae) to the top of the head, thereby covering the entire head and neck. Model components included: muscles, skeletal bones, lymph nodes, fat tissues, and organs such as salivary glands, tendons, andligaments. At first, one MRIimage dataset was selected as the reference image. The biometric properties (length, volume, mass, shape), hinge constraints of the bones, and the biomechanical properties of each of the anatomies were estimated using MRIs acquired at different head and neck poses. Results: The model's ability to represent different head and neck postures can be illustrated by observing the internal tissue deformations andthe model's ability to represent different postures. Conclusions: Results show that the biomechanical model was able to simulate different poses that may be exhibited during interfraction patient setup variations and intrafraction patient motion. Future work would focus on integrating dose calculations on the deforming model and validating the model deformations.

Author(s): Reyhan, Meral L; Wang, Zhe; Kim, Hyun J; Halnon, Nancy J; Finn, J Paul; Ennis, Daniel B | Abstract: PurposeCardiovascular magnetic resonance imaging exams can be performed during free-breathing. This may be especially important for boys with Duchenne muscular dystrophy (DMD) given their frequently limited breath-hold abilities. The impact of the respiratory compensation method on quantitative measurements of left ventricular (LV) rotational mechanics is incompletely understood. The purpose of this study was to evaluate differences in LV rotational mechanics acquired during breath-holding (BH), free-breathing with averaging (AVG), and free-breathing with respiratory bellows gating (BEL).MethodsLV short-axis tagged images from healthy subjects (N = 16) and DMD patients (N = 5) were acquired with BH, AVG, and BEL. LV twist and circumferential-longitudinal shear (CL-shear) angle were measured using the Fourier Analysis of STimulated echoes (FAST) method.ResultsPeak LV twist estimates using BEL were significantly lower compared with BH in both healthy subjects (10.2 ± 3.6 ° versus 12.9 ± 2.3 °, P = 0.003) and patients with DMD (8.6 ± 3.6 ° versus 10.5 ± 3.6 °, P = 0.004). AVG results were in between BEL and BH. No significant differences in CL-shear were detected between BEL and BH.ConclusionBreath-holding directly affects estimates of peak LV twist, but not CL-shear. Using a free-breathing strategy for the evaluation of cardiac function is important for intrasubject longitudinal studies, intersubject comparisons, and multicenter trials for patients with DMD. Magn Reson Med 77:864-869, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Author(s): Reyhan, Meral; Wang, Zhe; Li, Ming; Kim, Hyun J; Gupta, Himanshu; Lloyd, Steven G; Dell'Italia, Louis J; Denney, Thomas; Ennis, Daniel B | Abstract: PurposeTo evaluate the relationship between left ventricular (LV) twist, shear, and twist-per-volume and the severity of mitral regurgitation (MR). Primary MR is a valvular disorder that induces LV dysfunction. There exist several measures of LV rotational mechanics, but it remains unclear which measure of LV dysfunction best accords with the severity of MR. We hypothesized that LV systolic twist-per-volume slope would decrease with increasing severity of MR because of both decreases in rotational mechanics and increases in stroke volumes.Materials and methodsNormal subjects (n = 54), moderate MR patients (n = 29), and severe MR patients (n = 54) were studied. Magnetic resonance imaging (MRI) was performed on a 1.5T scanner and grid-tagged LV images were collected at the LV base and LV apex. Measures of LV rotational mechanics were derived from tagged images using Fourier Analysis of STimulated echoes (FAST).ResultsPeak systolic twist-per-volume slope was significantly different for all pairwise comparisons (P l 0.0001) and compared to normal subjects (-0.14 ± 0.05°/mL) was decreased in moderate MR (-0.12 ± 0.04°/mL) and further decreased in severe MR (-0.07 ± 0.03°/mL).ConclusionPeak systolic twist-per-volume slope significantly decreased with increasing severity of MR and is therefore a suitable quantitative imaging biomarker for LV dysfunction in patients with MR.

Purpose To assess the intra- and interscan reproducibility of LV twist using FAST. Assessing the reproducibility of the measurement of new MRI biomarkers is an important part of validation. Fourier Analysis of STimulated Echoes (FAST) is a new MRI tissue tagging method that has recently been shown to compare favorably with conventional estimates of left ventricular (LV) twist from cardiac tagged images, but with significantly reduced user interaction time. Materials and Methods Healthy volunteers (N = 10) were scanned twice using FAST over 1 week. On day 1, two measurements of LV twist were collected for intrascan comparisons. Measurements for LV twist were again collected on day 8 for interscan assessment. LV short-axis tagged images were acquired on a 3 Tesla (T) scanner to ensure detectability of tags during early and mid-diastole. Peak LV twist is reported as mean ± SD. Reproducibility was assessed using the concordance correlation coefficient (CCC) and the repeatability coefficient (RC) (95% confidence interval [CI] range). Results Mean peak twist measurements were 13.4 ± 4.3° (day 1, scan 1), 13.6 ± 3.7° (day 1, scan 2), and 13.0 ± 2.7° (day 8). Bland-Altman analysis resulted in intra- and interscan bias and 95% CI of −0.6° [−1.0°, 1.6°] and 1.4° (−1.0°, 3.0°), respectively. The Bland-Altman RC for peak LV twist was 2.6° and 4.0° for intra- and interscan, respectively. The CCC was 0.9 and 0.6 for peak LV twist for intra- and interscan, respectively. Conclusion FAST is a semi-automated method that provides a quick and quantitative assessment of LV systolic and diastolic twist that demonstrates high intrascan and moderate interscan reproducibility in preliminary studies. J. Magn. Reson. Imaging 2014;39:463–468. © 2013 Wiley Periodicals, Inc.

Purpose Cardiac cine balanced steady-state free precession (bSSFP) imaging uses a high flip angle (FA) to obtain high blood-myocardium signal-to-noise and contrast-to-noise ratios (CNR). Use of high FAs, however, results in substantially increased SAR. Our objective was to develop a variable FA bSSFP cardiac cine imaging technique with: (1) low SAR and blood-myocardium CNR similar to conventional constant FA bSSFP (CFA-bSSFP) or (2) increased blood-myocardium CNR compared to CFA-bSSFP with similar SAR. Methods Variable FA bSSFP cardiac cine imaging was achieved using an asynchronous k-space acquisition, which is asynchronous to the cardiac cycle (aVFA-bSSFP). Bloch simulations and phantom experiments were performed to compare the signal, resolution, and frequency response of the variable FA bSSFP and CFA-bSSFP schemes. Ten volunteers were imaged with different aVFA-bSSFP and asynchronous CFA-bSSFP schemes and compared to conventional segmented CFA-bSSFP. Results The SAR of aVFA-bSSFP is significantly decreased by 36% compared to asynchronous CFA-bSSFP (1.9 ± 0.2 vs. 3.0 ± 0.2 W/kg, P

Cardiovascular side effects secondary to cancer therapeutics are the leading causes of morbidity and mortality in cancer survivors. Anthracyclines (AC) serve as the backbone of many anti-cancer treatment strategies, but dose-dependent myocardial injury limits their use. Cumulative AC exposure can disrupt the dynamic equilibrium of the myocardial microarchitecture while repeated injury and repair leads to myocyte loss, interstitial myocardial fibrosis, and impaired contractility. Although children are assumed to have greater myocardial plasticity, AC exposure at a younger age portends worse prognosis. In older patients, there is lower overall survival once they develop cardiovascular disease. Because aberrations in the myocardial architecture predispose the heart to a decline in function, early detection with sensitive imaging tools is crucial and the implications for resource utilization are substantial. As a comprehensive imaging modality, cardiac magnetic resonance (CMR) imaging is able to go beyond quantification of ejection fraction and myocardial deformation to characterize adaptive microstructural and microvascular changes that are important to myocardial tissue health. Herein, we describe CMR as an established translational imaging tool that can be used clinically to characterize AC-associated myocardial remodeling.

Heart failure is a leading cause of death, yet its underlying electrophysiological (EP) mechanisms are not well understood. In this study, we use a multiscale approach to analyze a model of heart failure and connect its results to features of the electrocardiogram (ECG). The heart failure model is derived by modifying a previously validated electrophysiology model for a healthy rabbit heart. Specifically, in accordance with the heart failure literature, we modified the cell EP by changing both membrane currents and calcium handling. At the tissue level, we modeled the increased gap junction lateralization and lower conduction velocity due to downregulation of Connexin 43. At the biventricular level, we reduced the apex-to-base and transmural gradients of action potential duration (APD). The failing cell model was first validated by reproducing the longer action potential, slower and lower calcium transient, and earlier alternans characteristic of heart failure EP. Subsequently, we compared the electrical wave propagation in one dimensional cables of healthy and failing cells. The validated cell model was then used to simulate the EP of heart failure in an anatomically accurate biventricular rabbit model. As pacing cycle length decreases, both the normal and failing heart develop T-wave alternans, but only the failing heart shows QRS alternans (although moderate) at rapid pacing. Moreover, T-wave alternans is significantly more pronounced in the failing heart. At rapid pacing, APD maps show areas of conduction block in the failing heart. Finally, accelerated pacing initiated wave reentry and breakup in the failing heart. Further, the onset of VF was not observed with an upregulation of SERCA, a potential drug therapy, using the same protocol. The changes introduced at the cell and tissue level have increased the failing heart’s susceptibility to dynamic instabilities and arrhythmias under rapid pacing. However, the observed increase in arrhythmogenic potential is not due to a steepening of the restitution curve (not present in our model), but rather to a novel blocking mechanism., Author Summary Ventricular fibrillation (VF) is one of the leading causes of sudden death. During VF, the electrical wave of activation in the heart breaks up chaotically. Consequently, the heart is unable to contract synchronously and pump blood to the rest of the body. In our work we formulate and validate a model of heart failure (HF) that allows us to evaluate the arrhythmogenic potential of individual and combined electrophysiological changes. In diagnostic cardiology, the electrocardiogram (ECG) is one of the most commonly used tools for detecting abnormalities in the heart electrophysiology. One of our goals is to use our numerical model to link changes at the cellular and tissue level in a failing heart to a numerically computed ECG. This allows us to characterize the precursor to and the risk of VF. In order to understand the mechanisms underlying VF in HF, we design a test that simulates a HF patient performing physical exercise. We show that under fast heart rates with changes in pacing, HF patients are more prone to VF due to a new conduction blocking mechanism. In the long term, our mathematical model is suitable for investigating the effect of drug therapies in HF.

Catheter rupture during CT angiography has prompted policies prohibiting the use of electronic injectors with peripherally inserted central venous catheters (PICCs) not only for CT but also for MRI. Consequently, many institutions mandate hand injection for MR angiography, limiting precision of infusion rates and durations of delivery. To determine whether electronic injection of gadolinium-based contrast media through a range of small-caliber, single-lumen PICCs would be safe without risk of catheter rupture over the range of clinical protocols and determine whether programmed flow rates and volumes were realized when using PICCs for contrast delivery. Experiments were performed and recorded using the Medrad Spectris Solaris EP MR Injection System. PICC sizes, contrast media and flow rates were based on common institutional protocols. No catheters were damaged during any experiments. Mean difference between programmed and delivered volume was 0.07 ± 0.10 mL for all experiments. Reduced flow rates and prolonged injection durations were observed when the injector’s pressure-limiting algorithm was triggered, only in protocols outside the clinical range. PICCs commonly used in children can withstand in vitro power injection of gadolinium-based contrast media at protocols significantly above clinical levels.

Phase-contrast magnetic resonance imaging is subject to numerous sources of error, which decrease clinical confidence in the reported measures. This work outlines how stationary perivascular fat can impart a significant chemical shift induced phase-contrast magnetic resonance imaging measurement error using computational simulations, in vitro, and in vivo experiments. This chemical shift error does not subtract in phase difference processing, but can be minimized with proper parameter selection. The chemical shift induced phase errors largely depend on both the receiver bandwidth and the TE. Both theory and an in vivo comparison of the maximum difference in net forward flow between vessels with and without perivascular fat indicated that the effects of chemically shifted perivascular fat are minimized by the use of high bandwidth (814 Hz/px) and an in-phase TE (high BW-TE(IN)). In healthy volunteers (N = 10) high BW-TE(IN) significantly improves intrapatient net forward flow agreement compared with low bandwidth (401 Hz/px) and a mid-phase TE as indicated by significantly decreased measurement biases and limits of agreement for the ascending aorta (1.8 ± 0.5 mL vs. 6.4 ± 2.8 mL, P = 0.01), main pulmonary artery (2.0 ± 0.9 mL vs. 11.9 ± 5.8 mL, P = 0.04), the left pulmonary artery (1.3 ± 0.9 mL vs. 5.4 ± 2.5 mL, P = 0.003), and all vessels (1.7 ± 0.8 mL vs. 7.2 ± 4.4 mL, P = 0.001).

Reyhan et al. Journal of Cardiovascular Magnetic Resonance 2015, 17(Suppl 1):Q22 http://www.jcmr-online.com/content/17/S1/Q22 WALKING POSTER PRESENTATION Open Access The effect of free-breathing on left ventricular rotational mechanics in normal subjects and patients with duchenne muscular dystrophy Meral Reyhan 2,3 , Zhe Wang 2,1* , Hyun J Kim 2 , Nancy Halnon 4 , Paul J Finn 2 , Daniel B Ennis 2,1 From 18th Annual SCMR Scientific Sessions Nice, France. 4-7 February 2015 Background Patients with Duchenne Muscular Dystrophy (DMD) develop respiratory impairment and signs and symptoms of cardiac involvement at an early age. Quantitative mea- sures of cardiac function from MRI tagging can be used to stage disease progression and monitor the response to therapy. Traditional cardiac MRI exams, however, require repeated breath holding, which places a notable burden on patients with DMD and can impart intrathoracic pressure changes that alter cardiac function. The effect of free- breathing on quantitative measure of left ventricular (LV) rotational mechanics for DMD patients is incompletely understood. Objective: To evaluate differences in LV rota- tional mechanics acquired during breath holding (BH), free-breathing with averaging (AVG), and free-breathing with respiratory bellows gating (BEL). Methods 16 healthy volunteers (1 female, 27.8±3.9 years) and 5 DMD patients (all male, 11.3±3.7 years) were scanned on a 3.0T Siemens Scanner with BH, AVG and BEL. ORI-SPAMM (Reyhan. M et. al. JMRI 2014;39(2):339- 345) was used to acquire short-axis images at the base and apex of the heart. 280-330x280-330mm field-of- view, 6mm slice thickness, 192x192 acquisition matrix, 395 Hz/pixel receiver bandwidth, TE/TR =2.33-2.39/ 4.71-4.83 ms, 8 phase encode lines per segment, 12° imaging flip angle, 8 mm tag spacing, GRAPPA 2. LV Figure 1 Mean LV twist of (A) healthy volunteers and (B) patients with DMD for breath-held (BH, solid), free-breathing with averaging (AVG, dashed), and bellows gated (BEL, dash-dotted) methods. BH LV twist is generally higher compared to AVG or BEL, especially at its peak. Radiological Sciences, University of California Los Angeles, Los Angeles, CA, USA © 2015 Reyhan et al; licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http:// creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Purpose: To validate a novel method for the rapid and facile quantification of left ventricular (LV) twist from tagged magnetic resonance images and demonstrate the potential clinical utility in a series of 20 healthy volunteers. Materials and Methods: Cardiac magnetic resonance imaging (MRI) short-axis images were acquired with tissue tagging in 20 healthy subjects and six canines. The tagged images were processed using a novel Fourier Analysis of the STimulated echoes (FAST) method, which uses a series of Fourier-space operations to measure LV twist with limited user interaction. A subset of eight healthy subjects and the canine data were compared to results from previously validated “gold standard” software (FindTags). Interobserver and intraobserver coefficients of variation (CVINTER and CVINTRA), linear regression, and Bland–Altman analyses were used to assess agreement between observers and methods. Results: CVINTRA for peak systolic twist (2.9% and 2.6%) and CVINTER (4.3% and 4.2%) were all small. Linear regression analysis of the FAST and FindTags twist values indicated very good agreement in healthy subjects (R = 0.91) and in canines (R = 0.95). Bland–Altman comparison of the FAST and FindTags twist results indicated excellent agreement in healthy subjects (bias of −0.5°, 95% confidence intervals (−4.3°, 4.3°)) and canines (bias of 0.2°, 95% confidence intervals (−2.7°, 3.1°)). Peak systolic twist in healthy subjects averaged 10.5 ± 1.9° degrees. Conclusion: The FAST method for quantifying LV twist produces results that are not significantly different from the current “gold standard” in a fraction of the user interaction time and has demonstrated feasibility in human subjects. J. Magn. Reson. Imaging 2012;35:587-593. © 2011 Wiley Periodicals, Inc.

Background— Remodeling of the mitral annulus contributes to progression of mitral regurgitation (MR). In patients with moderate-to-severe MR, short-term treatment with β-blockers has been shown to increase left ventricular (LV) end-diastolic and end-systolic volume, and this could deleteriously increase mitral valve annular dimensions. The objective of this study was to quantify the effects of a short duration of β-blocker treatment on mitral annular dimensions and dynamics in patients with MR due to primary degenerative valve disease. Methods and Results— Twenty-five patients with moderate-to-severe degenerative MR and normal LV systolic function were studied in a double-blind crossover experiment using a β1-selective adrenergic blocker and placebo administered for 14±3 days. Cardiac MRI images were acquired after each treatment period to quantify mitral annular dimensions. At end diastole, there was no change in annular area (1659±331 versus 1632±299 mm 2 ; P P P P P P P P P Conclusions— Despite significant increases in LV end-diastolic and end-systolic volume, short-term β-blocker treatment of patients with moderate-to-severe MR reduced or preserved all mitral annular dimensions except SL at ES.

© 2015 IEEE. Cardiac electroanatomical mapping (EAM) to locate myocardial scar substrate during ventricular tachycardia (VT) ablation therapy is currently faced with a major limitation. The presence of epicardial fat leads to false-positive low-voltage maps that may interpreted as myocardial scar. Pre-procedural cardiac magnetic resonance (CMR) imaging can produce images of the heart in which both cardiac geometry and epicardial fat are particularly conspicuous. The goal of the present work is to build CMR based geometric models composed of a shell of the cardiac anatomy and epicardial fat. Herein, we first defined a black-blood imaging protocol that is feasible in patients with implanted devices. Then a semi-automated image segmentation method was developed and applied to extract ventricular surface anatomy and epicardial fat maps from which three-dimensional surfaces were built. The results suggest that deriving maps of epicardial fat from MRI data is feasible, accurate when compared to expert observers, and suitable for integration in a clinical catheter-based ablation procedure.

In mammalian hearts the capillaries are closely aligned with the muscle fibers. We report our observation of a main-field direc- tion-dependent contrast in MR spin-echo (SE) images of the heart in the presence of Ferumoxtran-10, an intravascular iron- oxide nanoparticle contrast agent (CA). We describe a novel MRI method for mapping the preferential orientation of capil- laries in the myocardial wall. The eigenvector corresponding to the minimum eigen value of the R2 relaxation rate tensor is consistent with the expected orientation of the capillary net- work. Preliminary results also demonstrate the feasibility of this method for in vivo application to rodent imaging. Magn Reson Med 55:725-730, 2006. Published 2006 Wiley-Liss, Inc. † The transmural helical structure of the myofibers in the heart and the consequent contraction pattern has been known for centuries (1-3), and so has the fact that the capillaries are closely aligned with the muscle fibers (3). We report our observations of a main magnetic field, B0, direction-dependent contrast in MR spin-echo (SE) images of the heart in the presence of intravascular iron-oxide nanoparticles. The main external field direction depen- dence makes it possible to calculate the tensor compo- nents of the direction-dependent transverse relaxation rate (R2 or 1/T2). The tertiary principal direction of the R2 tensor coincides with the capillary orientation within the myocardial wall, and the correlation between the two can be explained by the physics of proton spin relaxation around the mostly parallel capillary vessels. The nanoparticle contrast agent (CA) Ferumoxtran-10 is an aqueous solution that contains magnetite crystalline cores that are coated with dextran. They are too large to move across the pores of the endothelium of the capillaries and thus are confined to the intravascular space while they are slowly absorbed by the reticuloendothelial system of the liver. Because of their superparamagnetic nature they are fully magnetized in the field of a 1.5 Tesla MRI scanner (4), and create strong, heterogeneous, and localized fields around capillary vessels. In SE imaging, proton spins of tissue water undergo diffusive Brownian motion in these induced fields that accelerates T2 relaxation. In the heart the capillaries are morphologically nearly parallel tubes that are aligned with the muscle fibers, and the field pat- tern around such a capillary is highly dependent on its orientation relative to the external main field (5), as illus- trated in Fig. 1. The strongest local field gradients occur when the capillary is perpendicular to the main field, and as demonstrated by Monte-Carlo simulations of a capillary lattice in the heart (6), this condition should create the maximum T2 relaxation and the lowest signal level in SE images. We observed in both beating hearts and excised heart samples this variation of SE image intensity with orienta- tion, as well as heterogeneous image intensity in the myo- cardial wall. We ruled out other potential instrumentation factors that might cause such heterogeneity with control images acquired using short echo times (TEs) and thus minimal T2 relaxation. To characterize the directional de- pendence in the R2 relaxation rate, we derived the second- order R2 tensor from a set of images of an excised heart sample. A map of the tertiary eigenvector of the tensor, indicating the direction of minimum R2, was then plotted throughout the myocardial wall. The results can be fully explained by the anisotropic R2 relaxation in the presence of largely parallel capillaries, in the light of the well- known transmural helical patterns of the myofiber and capillary orientation of the heart (2) (i.e., the direction of minimum R2 relaxation coincides with the myofiber ori- entation in all locations of the heart). This observation provides a robust means of mapping the capillary orientation in vivo with simple SE imaging, and reducing susceptibility to physiological motion due to the lack of large motion-encoding gradients.

We describe a sequence of methods to produce a partial differential equation model of the electrical activation of the ventricles. In our framework, we incorporate the anatomy and cardiac microstructure obtained from magnetic resonance imaging and diffusion tensor imaging of a New Zealand White rabbit, the Purkinje structure and the Purkinje-muscle junctions, and an electrophysiologically accurate model of the ventricular myocytes and tissue, which includes transmural and apex-to-base gradients of action potential characteristics. We solve the electrophysiology governing equations using the finite element method and compute both a 6-lead precordial electrocardiogram (ECG) and the activation wavefronts over time. We are particularly concerned with the validation of the various methods used in our model and, in this regard, propose a series of validation criteria that we consider essential. These include producing a physiologically accurate ECG, a correct ventricular activation sequence, and the inducibility of ventricular fibrillation. Among other components, we conclude that a Purkinje geometry with a high density of Purkinje muscle junctions covering the right and left ventricular endocardial surfaces as well as transmural and apex-to-base gradients in action potential characteristics are necessary to produce ECGs and time activation plots that agree with physiological observations., Shankarjee Krishnamoorthi et. al

To test the hypothesis that alterations in electrical activation sequence contribute to depressed systolic function in the infarct border zone, we examined the anatomic correlation of abnormal electromechanics and infarct geometry in the canine post-myocardial infarction (MI) heart, using a high-resolution MR-based cardiac electromechanical mapping technique. Three to eight weeks after an MI was created in six dogs, a 247-electrode epicardial sock was placed over the ventricular epicardium under thoracotomy. MI location and geometry were evaluated with delayed hyperenhancement MRI. Three-dimensional systolic strains in epicardial and endocardial layers were measured in five short-axis slices with motion-tracking MRI (displacement encoding with stimulated echoes). Epicardial electrical activation was determined from sock recordings immediately before and after the MR scans. The electrodes and MR images were spatially registered to create a total of 160 nodes per heart that contained mechanical, transmural infarct extent, and electrical data. The average depth of the infarct was 55% (SD 11), and the infarct covered 28% (SD 6) of the left ventricular mass. Significantly delayed activation (>mean + 2SD) was observed within the infarct zone. The strain map showed abnormal mechanics, including abnormal stretch and loss of the transmural gradient of radial, circumferential, and longitudinal strains, in the region extending far beyond the infarct zone. We conclude that the border zone is characterized by abnormal mechanics directly coupled with normal electrical depolarization. This indicates that impaired function in the border zone is not contributed by electrical factors but results from mechanical interaction between ischemic and normal myocardium.

Diastolic and systolic left ventricular (LV) dysfunction often significantly contribute to disabling symptoms in familial hypertrophic cardiomyopathy (FHC). This study compares regional LV function (midwall circumferential strain) during systole and diastole in eight FHC patients and six normal volunteers (NVs) using MR tagging. A prospectively-gated fast gradient-echo sequence with an echo-train readout was modified to support complementary spatial modulation of magnetization (CSPAMM) tagging and full cardiac cycle data acquisition using the cardiac phase to order reconstruction (CAPTOR), thus providing tag persistence and data acquisition during the entire cardiac cycle. Total systolic strains in FHC patients were significantly reduced in septal and inferior regions (both P < 0.01). Early-diastolic strain rates were reduced in all regions of the FHC group (all P < 0.03). The combination of CSPAMM and CAPTOR allows regional indices of myocardial function to be quantified throughout the cardiac cycle. This technique reveals regional differences in systolic and diastolic impairment in FHC patients.

Background Cardiac diffusion tensor imaging (cDTI) measures the magnitudes and directions of intramyocardial water diffusion. Assuming the cross-myocyte components to be constrained by the laminar microstructures of myocardium, we hypothesized that cDTI at two cardiac phases might identify any abnormalities of laminar orientation and mobility in hypertrophic cardiomyopathy (HCM). Methods We performed cDTI in vivo at 3 Tesla at end-systole and late diastole in 11 healthy controls and 11 patients with HCM, as well as late gadolinium enhancement (LGE) for detection of regional fibrosis. Results Voxel-wise analysis of diffusion tensors relative to left ventricular coordinates showed expected transmural changes of myocardial helix-angle, with no significant differences between phases or between HCM and control groups. In controls, the angle of the second eigenvector of diffusion (E2A) relative to the local wall tangent plane was larger in systole than diastole, in accord with previously reported changes of laminar orientation. HCM hearts showed higher than normal global E2A in systole (63.9° vs 56.4° controls, p = 0.026) and markedly raised E2A in diastole (46.8° vs 24.0° controls, p

Terahertz (THz) detection has been proposed and applied to a variety of medical imaging applications in view of its unrivaled hydration profiling capabilities. Variations in tissue dielectric function have been demonstrated at THz frequencies to generate high contrast imagery of tissue, however, the source of image contrast remains to be verified using a modality with a comparable sensing scheme. To investigate the primary contrast mechanism, a pilot comparison study was performed in a burn wound rat model, widely known to create detectable gradients in tissue hydration through both injured and surrounding tissue. Parallel T2 weighted multi slice multi echo (T2w MSME) 7T Magnetic Resonance (MR) scans and THz surface reflectance maps were acquired of a full thickness skin burn in a rat model over a 5 hour time period. A comparison of uninjured and injured regions in the full thickness burn demonstrates a 3-fold increase in average T2 relaxation times and a 15% increase in average THz reflectivity, respectively. These results support the sensitivity and specificity of MRI for measuring in vivo burn tissue water content and the use of this modality to verify and understand the hydration sensing capabilities of THz imaging for acute assessments of the onset and evolution of diseases that affect the skin. A starting point for more sophisticated in vivo studies, this preliminary analysis may be used in the future to explore how and to what extent the release of unbound water affects imaging contrast in THz burn sensing.

Background Late gadolinium enhancement (LGE) magnetic resonance imaging (MRI) of ventricular scar has been shown to be accurate for detection and characterization of arrhythmia substrates. However, the majority of patients referred for ventricular tachycardia (VT) ablation have an implantable cardioverter-defibrillator (ICD), which obscures image integrity and the clinical utility of MRI. Objective The purpose of this study was to develop and validate a wideband LGE MRI technique for device artifact removal. Methods A novel wideband LGE MRI technique was developed to allow for improved scar evaluation on patients with ICDs. The wideband technique and the standard LGE MRI were tested on 18 patients with ICDs. VT ablation was performed in 13 of 18 patients with either endocardial and/or epicardial approach and the correlation between the scar identified on MRI and electroanatomic mapping (EAM) was analyzed. Results Hyperintensity artifact was present in 16 of 18 of patients using standard MRI, which was eliminated using the wideband LGE and allowed for MRI interpretation in 15 of 16 patients. All patients had ICD lead characteristics confirmed as unchanged post-MRI and had no adverse events. LGE scar was seen in 11 of 18 patients. Among the 15 patients in whom wideband LGE allowed visualization of myocardium, 10 had LGE scar and 5 had normal myocardium in the regions with image artifacts when using the standard LGE. The left ventricular scar size measurements using wideband MRI and EAM were correlated with R2= 0.83 and P =.00003. Conclusion Wideband LGE MRI improves the ability to visualize myocardium for clinical interpretation, which correlated well with EAM findings during VT ablation. © 2014 Heart Rhythm Society.

Accurate and early prediction of tissue viability is the most significant determinant of tissue flap survival in reconstructive surgery. Perturbation in tissue water content (TWC) is a generic component of the tissue response to such surgeries, and, therefore, may be an important diagnostic target for assessing the extent of flap viability in vivo. We have previously shown that reflective terahertz (THz) imaging, a non-ionizing technique, can generate spatially resolved maps of TWC in superficial soft tissues, such as cornea and wounds, on the order of minutes. Herein, we report the first in vivo pilot study to investigate the utility of reflective THz TWC imaging for early assessment of skin flap viability. We obtained longitudinal visible and reflective THz imagery comparing 3 bipedicled flaps (i.e. survival model) and 3 fully excised flaps (i.e. failure model) in the dorsal skin of rats over a postoperative period of 7 days. While visual differences between both models manifested 48 hr after surgery, statistically significant (p < 0.05, independent t-test) local differences in TWC contrast were evident in THz flap image sets as early as 24 hr. Excised flaps, histologically confirmed as necrotic, demonstrated a significant, yet localized, reduction in TWC in the flap region compared to non-traumatized skin. In contrast, bipedicled flaps, histologically verified as viable, displayed mostly uniform, unperturbed TWC across the flap tissue. These results indicate the practical potential of THz TWC sensing to accurately predict flap failure 24 hours earlier than clinical examination.

Prostate imaging commonly performed using transrectal ultrasound (TRUS), has been unable to image the prostate in real-time 3D. Transurethral ultrasound (TUUS), a new method of prostate imaging, may provide a solution. We developed a 3D reconstruction method for transurethral ultrasound using a radially-phased catheter transducer, and performed an in vitro validation and comparison with transrectal ultrasound (TRUS) and magnetic resonance imaging (MRI). The 3D reconstruction method was able to generate TUUS volumes in a near real-time fashion. Calculated volumes of a prostate phantom were similar in all three methods. Inserted fiducial markers were better visualized in TUUS. MR-TUUS fusion was demonstrated using current registration algorithms optimized for TRUS. Future improvements to the volume reconstruction and MR fusion accuracy may enable TUUS to become an aid in prostate cancer diagnosis.

Biomechanical models are employed for guiding the image integration of macro-scale images obtained from Magnetic Resonance Imaging (MRI) with micro-scale images obtained from Optical Coherence Microscopy (OCM) designed for use in image-guided clinical interventions for the head and neck region.

Impaired systolic ventricular function is common in patients diagnosed with heart failure (HF) or ischaemic heart disease. The diminished contractile performance with impaired contractility (systolic HF) can be induced by impaired filling function (diastolic HF) and the wall stress (both passive and active) may indicate the progression from diastolic HF to systolic HF. In order to better understand the distribution of active stress during ventricular contraction, a left ventricular (LV) finite element (FE) model incorporating LV fibre geometry and function was developed to parameterise a time-varying model of myocardial contraction by simulating LV mechanics. During systole, the isometric active stress monotonically increased to 95 kPa, and rapidly recovered during isovolumic relaxation. We also observed regional variations of the fibre length dependent contractile stress throughout the LV. The time-varying active stress curve thereby obtained enabled quantification of heart muscle performance. This type of integrative modelling enables the investigation of LV mechanics on an individualised basis.

The majority of patients with clinically diagnosed heart failure have normal systolic pump function and are commonly categorized as suffering from diastolic heart failure. The left ventricle (LV) remodels its structure and function to adapt to pathophysiological changes in geometry and loading conditions, which in turn can alter the passive ventricular mechanics. In order to better understand passive ventricular mechanics, a LV finite element (FE) model was customized to geometric data segmented from in vivo tagged magnetic resonance images (MRI) data and myofibre orientation derived from ex vivo diffusion tensor MRI (DTMRI) of a canine heart using nonlinear finite element fitting techniques. MRI tissue tagging enables quantitative evaluation of cardiac mechanical function with high spatial and temporal resolution, whilst the direction of maximum water diffusion in each voxel of a DTMRI directly corresponds to the local myocardial fibre orientation. Due to differences in myocardial geometry between in vivo and ex vivo imaging, myofibre orientations were mapped into the geometric FE model using host mesh fitting (a free form deformation technique). Pressure recordings, temporally synchronized to the tagging data, were used as the loading constraints to simulate the LV deformation during diastole. Simulation of diastolic LV mechanics allowed us to estimate the stiffness of the passive LV myocardium based on kinematic data obtained from tagged MRI. Integrated physiological modelling of this kind will allow more insight into mechanics of the LV on an individualized basis, thereby improving our understanding of the underlying structural basis of mechanical dysfunction under pathological conditions.

Patients suffering from dilated cardiomyopathy or myocardial infarction can develop left ventricular (LV) diastolic impairment. The LV remodels its structure and function to adapt to pathophysiological changes in geometry and loading conditions and this remodeling process can alter the passive ventricular mechanics. In order to better understand passive ventricular mechanics, a LV finite element model was developed to incorporate physiological and mechanical information derived from in vivo magnetic resonance imaging (MRI) tissue tagging, in vivo LV cavity pressure recording and ex vivo diffusion tensor MRI (DTMRI) of a canine heart. MRI tissue tagging enables quantitative evaluation of cardiac mechanical function with high spatial and temporal resolution, whilst the direction of maximum water diffusion (the primary eigenvector) in each voxel of a DTMRI directly correlates with the myocardial fibre orientation. This model was customized to the geometry of the canine LV during diastasis by fitting the segmented epicardial and endocardial surface data from tagged MRI using nonlinear finite element fitting techniques. Myofibre orientations, extracted from DTMRI of the same heart, were incorporated into this geometric model using a free form deformation methodology. Pressure recordings, temporally synchronized to the tissue tagging MRI data, were used to simulate the LV deformation during diastole. Simulation of the diastolic LV mechanics allowed us to estimate the stiffness of the passive LV myocardium based on kinematic data obtained from tagged MRI. This integrated physiological model will allow more insight into the regional passive diastolic mechanics of the LV on an individualized basis, thereby improving our understanding of the underlying structural basis of mechanical dysfunction in pathological conditions.

We measured leaflet displacements and used inverse finite-element analysis to define, for the first time, the material properties of mitral valve (MV) leaflets in vivo. Sixteen miniature radiopaque markers were sewn to the MV annulus, 16 to the anterior MV leaflet, and 1 on each papillary muscle tip in 17 sheep. Four-dimensional coordinates were obtained from biplane videofluoroscopic marker images (60 frames/s) during three complete cardiac cycles. A finite-element model of the anterior MV leaflet was developed using marker coordinates at the end of isovolumic relaxation (IVR; when the pressure difference across the valve is ∼0), as the minimum stress reference state. Leaflet displacements were simulated during IVR using measured left ventricular and atrial pressures. The leaflet shear modulus ( Gcirc-rad) and elastic moduli in both the commisure-commisure ( Ecirc) and radial ( Erad) directions were obtained using the method of feasible directions to minimize the difference between simulated and measured displacements. Group mean (±SD) values (17 animals, 3 heartbeats each, i.e., 51 cardiac cycles) were as follows: Gcirc-rad= 121 ± 22 N/mm2, Ecirc= 43 ± 18 N/mm2, and Erad= 11 ± 3 N/mm2( Ecirc> Erad, P < 0.01). These values, much greater than those previously reported from in vitro studies, may result from activated neurally controlled contractile tissue within the leaflet that is inactive in excised tissues. This could have important implications, not only to our understanding of mitral valve physiology in the beating heart but for providing additional information to aid the development of more durable tissue-engineered bioprosthetic valves.

The mitral annulus is a dynamic, saddle-shaped structure consisting of fibrous and muscular regions. Normal physiologic mechanisms of annular motion are incompletely understood, and more complete characterization is needed to provide rational basis for annuloplasty ring design and to enhance clinical outcomes.Seventeen sheep had radiopaque markers implanted; 16 around the annulus and 2 on middle anterior and posterior leaflet edges. Four-dimensional marker coordinates were acquired with biplanar videofluoroscopy at 60 Hz. Hinge angle was quantified between fibrous and muscular annular planes, with 0 degrees defined at end diastole, to characterize its contribution to alterations in mitral septal-lateral dimension and 2-dimensional total annular area throughout the cardiac cycle.During isovolumic contraction (pre-ejection), hinge angle abruptly increased, reaching maximum (steepest saddle shape, change 18 degrees +/- 13 degrees ) at peak left ventricular pressure. During ejection, hinge angle did not change; it then decreased during early filling (change 2 degrees +/- 2 degrees ). Septal-lateral dimension and total area paralleled hinge angle dynamics and leaflet distance (anterior to posterior marker). Pre-ejection septal-lateral reduction was 13% +/- 7% (3.3 +/- 1.5 mm) from 9% muscular dimension fall and 18 degrees +/- 13 degrees hinge angle increase.Pre-ejection increase in hinge angle contributes substantially to septal-lateral and total area reduction, facilitating leaflet coaptation. Semirigid annuloplasty rings or partial bands may preserve hinge motion, but possible recurrent annular dilatation could result in recurrent mitral regurgitation. Long-term clinical studies are required to determine who might benefit most from preserving intrinsic hinge motion without compromising repair durability.

The electrocardiogram (ECG) is one of the most significant outputs of a computational model of cardiac electrophysiology because it relates the numerical results to clinical data and is a universal tool for diagnosing heart diseases. One key features of the ECG is the T-wave, which is caused by longitudinal and transmural heterogeneity of the action potential duration (APD). Thus, in order to model a correct wave of repolarization, different cell properties resulting in different APDs must be assigned across the ventricular wall and longitudinally from apex to base. To achieve this requirement, a regional parametrization of the heart is necessary. We propose a robust approach to obtain the transmural and longitudinal segmentation in a general heart geometry without relying on ad hoc procedures. Our approach is based on auxiliary harmonic lifting analyses, already used in the literature to generate myocardial fiber orientations. Specifically, the solution of a sequence of Laplace boundary value problems allows parametrically controlled segmentation of both heart ventricles. The flexibility and simplicity of the proposed method is demonstrated through several representative examples, varying the locations and extents of the epicardial, midwall, and endocardial layers. Effects of the control parameters on the T-wave morphology are illustrated via computed ECGs.

Near-simultaneous measurements of electrical and mechanical activation over the entire ventricular surface are now possible using magnetic resonance imaging tagging and a multielectrode epicardial sock. This new electromechanical mapping technique is demonstrated in the ventricularly paced canine heart. A 128–electrode epicardial sock and pacing electrodes were placed on the hearts of four anesthetized dogs. In the magnetic resonance scanner, tagged cine images (8–15 ms/frame) and sock electrode recordings (1000 Hz) were acquired under right-ventricular pacing and temporally referenced to the pacing stimulus. Electrical recordings were obtained during intermittent breaks in image acquisition, so that both data sets represented the same physiologic state. Since the electrodes were not visible in the images, electrode recordings and cine images were spatially registered with Gd-DTPA markers attached to the sock. Circumferential strain was calculated at locations corresponding to electrodes. For each electrode location, electrical and mechanical activation times were calculated and relationships between the two activation patterns were demonstrated. This method holds promise for improving understanding of the relationships between the patterns of electrical activation and contraction in the heart. © 2003 Biomedical Engineering Society. PAC2003: 8761Pk, 8719Rr, 8719Hh