Diffusion Tractography of the Entire Left Ventricle by Using Free-breathing Accelerated Simultaneous Multisection Imaging

Cardiac diffusion-weighted and tensor imaging: A consensus statement from the special interest group of the Society for Cardiovascular Magnetic Resonance

Thanks to recent developments in cardiovascular magnetic resonance (CMR), cardiac diffusion-weighted magnetic resonance is fast emerging in a range of clinical applications. Cardiac diffusion-weighted imaging (cDWI) and diffusion tensor imaging (cDTI) now enable investigators and clinicians to assess and quantify the tridimensional microstructure of the heart. Free-contrast DWI is uniquely sensitized to the presence and displacement of water molecules within the myocardial tissue, including the intracellular, extracellular, and intravascular spaces. CMR can determine changes in microstructure by quantifying: a) mean diffusivity (MD)-measuring the magnitude of diffusion; b) fractional anisotropy (FA)-specifying the directionality of diffusion; c) helix angle (HA) and transverse angle (TA)-indicating the orientation of the cardiomyocytes; d) absolute sheetlet angle (E2A) and E2A mobility-measuring the alignment and systolic-diastolic mobility of the sheetlets, respectively. This document provides recommendations for both clinical and research cDWI and cDTI, based on published evidence when available and expert consensus when not. It introduces the cardiac microstructure focusing on the cardiomyocytes and their role in cardiac physiology and pathophysiology. It highlights methods, observations, and recommendations in terminology, acquisition schemes, postprocessing pipelines, data analysis, and interpretation of the different biomarkers. Despite the ongoing challenges discussed in the document and the need for ongoing technical improvements, it is clear that cDTI is indeed feasible, can be accurately and reproducibly performed and, most importantly, can provide unique insights into myocardial pathophysiology.

Soft robotic artificial left ventricle simulator capable of reproducing myocardial biomechanics

The heart's intricate myocardial architecture has been called the Gordian knot of anatomy, an impossible tangle of intricate muscle fibers. This complexity dictates equally complex cardiac motions that are difficult to mimic in physical systems. If these motions could be generated by a robotic system, then cardiac device testing, cardiovascular disease studies, and surgical procedure training could reduce their reliance on animal models, saving time, costs, and lives. This work introduces a bioinspired soft robotic left ventricle simulator capable of reproducing the minutiae of cardiac motion while providing physiological pressures. This device uses thin-filament artificial muscles to mimic the multilayered myocardial architecture. To demonstrate the device's ability to follow the cardiac motions observed in the literature, we used canine myocardial strain data as input signals that were subsequently applied to each artificial myocardial layer. The device's ability to reproduce physiological volume and pressure under healthy and heart failure conditions, as well as effective simulation of a cardiac support device, were experimentally demonstrated in a left-sided mock circulation loop. This work also has the potential to deliver faithful simulated cardiac motion for preclinical device and surgical procedure testing, with the potential to simulate patient-specific myocardial architecture and motion.

Diffusion Tensor Phenomapping of the Healthy and Pressure-Overloaded Human Heart

Current techniques to image the microstructure of the heart with diffusion tensor MRI (DTI) are highly under-resolved. We present a technique to improve the spatial resolution of cardiac DTI by almost 10-fold and leverage this to measure local gradients in cardiomyocyte alignment or helix angle (HA). We further introduce a phenomapping approach based on voxel-wise hierarchical clustering of these gradients to identify distinct microstructural microenvironments in the heart. Initial development was performed in healthy volunteers (n=8). Thereader, subjects with severe but well-compensated aortic stenosis (AS, n=10) were compared to age-matched controls (CTL, n=10). Radial HA gradient was significantly reduced in AS (8.0±0.8°/mm vs. 10.2±1.8°/mm, p=0.001) but the other HA gradients did not change significantly. Four distinct microstructural clusters could be idenJfied in both the CTL and AS subjects and did not differ significantly in their properties or distribution. Despite marked hypertrophy, our data suggest that the myocardium in well-compensated AS can maintain its microstructural coherence. The described phenomapping approach can be used to characterize microstructural plasticity and perturbation in any organ system and disease.

Right Ventricular Architectural Remodeling and Functional Adaptation in Pulmonary Hypertension

BACKGROUND

Global indices of right ventricle (RV) function provide limited insights into mechanisms underlying RV remodeling in pulmonary hypertension (PH). While RV myocardial architectural remodeling has been observed in PH, its effect on RV adaptation is poorly understood.

METHODS

Hemodynamic assessments were performed in 2 rodent models of PH. RV free wall myoarchitecture was quantified using generalized Q-space imaging and tractography analyses. Computational models were developed to predict RV wall strains. Data from animal studies were analyzed to determine the correlations between hemodynamic measurements, RV strains, and structural measures.

RESULTS

In contrast to the PH rats with severe RV maladaptation, PH rats with mild RV maladaptation showed a decrease in helical range of fiber orientation in the RV free wall (139º versus 97º; P=0.029), preserved global circumferential strain, and exhibited less reduction in right ventricular-pulmonary arterial coupling (0.029 versus 0.017 mm/mm Hg; P=0.037). Helical range correlated positively with coupling (P=0.036) and stroke volume index (P<0.01). Coupling correlated with global circumferential strain (P<0.01) and global radial strain (P<0.01) but not global longitudinal strain.

CONCLUSIONS

Data analysis suggests that adaptive RV architectural remodeling could improve RV function in PH. Our findings suggest the need to assess RV architecture within routine screenings of PH patients to improve our understanding of its prognostic and therapeutic significance in PH.

An Overview of Deep Learning Methods for Left Ventricle Segmentation

Cardiac health diseases are one of the key causes of death around the globe. The number of heart patients has considerably increased during the pandemic. Therefore, it is crucial to assess and analyze the medical and cardiac images. Deep learning architectures, specifically convolutional neural networks have profoundly become the primary choice for the assessment of cardiac medical images. The left ventricle is a vital part of the cardiovascular system where the boundary and size perform a significant role in the evaluation of cardiac function. Due to automatic segmentation and good promising results, the left ventricle segmentation using deep learning has attracted a lot of attention. This article presents a critical review of deep learning methods used for the left ventricle segmentation from frequently used imaging modalities including magnetic resonance images, ultrasound, and computer tomography. This study also demonstrates the details of the network architecture, software, and hardware used for training along with publicly available cardiac image datasets and self-prepared dataset details incorporated. The summary of the evaluation matrices with results used by different researchers is also presented in this study. Finally, all this information is summarized and comprehended in order to assist the readers to understand the motivation and methodology of various deep learning models, as well as exploring potential solutions to future challenges in LV segmentation.

Accelerating Cardiac Diffusion Tensor Imaging With a U-Net Based Model: Toward Single Breath-Hold

BACKGROUND

In vivo cardiac diffusion tensor imaging (cDTI) characterizes myocardial microstructure. Despite its potential clinical impact, considerable technical challenges exist due to the inherent low signal-to-noise ratio.

PURPOSE

To reduce scan time toward one breath-hold by reconstructing diffusion tensors for in vivo cDTI with a fitting-free deep learning approach.

STUDY TYPE

Retrospective.

POPULATION

A total of 197 healthy controls, 547 cardiac patients.

FIELD STRENGTH/SEQUENCE

A 3 T, diffusion-weighted stimulated echo acquisition mode single-shot echo-planar imaging sequence.

ASSESSMENT

A U-Net was trained to reconstruct the diffusion tensor elements of the reference results from reduced datasets that could be acquired in 5, 3 or 1 breath-hold(s) (BH) per slice. Fractional anisotropy (FA), mean diffusivity (MD), helix angle (HA), and sheetlet angle (E2A) were calculated and compared to the same measures when using a conventional linear-least-square (LLS) tensor fit with the same reduced datasets. A conventional LLS tensor fit with all available data (12 ± 2.0 [mean ± sd] breath-holds) was used as the reference baseline.

STATISTICAL TESTS

Wilcoxon signed rank/rank sum and Kruskal-Wallis tests. Statistical significance threshold was set at P = 0.05. Intersubject measures are quoted as median [interquartile range].

RESULTS

For global mean or median results, both the LLS and U-Net methods with reduced datasets present a bias for some of the results. For both LLS and U-Net, there is a small but significant difference from the reference results except for LLS: MD 5BH (P = 0.38) and MD 3BH (P = 0.09). When considering direct pixel-wise errors the U-Net model outperformed significantly the LLS tensor fit for reduced datasets that can be acquired in three or just one breath-hold for all parameters.

DATA CONCLUSION

Diffusion tensor prediction with a trained U-Net is a promising approach to minimize the number of breath-holds needed in clinical cDTI studies.

EVIDENCE LEVEL

4 TECHNICAL EFFICACY: Stage 1.

The relationship between myocardial microstructure and strain in chronic infarction using cardiovascular magnetic resonance diffusion tensor imaging and feature tracking

BACKGROUND

Cardiac diffusion tensor imaging (cDTI) using cardiovascular magnetic resonance (CMR) is a novel technique for the non-invasive assessment of myocardial microstructure. Previous studies have shown myocardial infarction to result in loss of sheetlet angularity, derived by reduced secondary eigenvector (E2A) and reduction in subendocardial cardiomyocytes, evidenced by loss of myocytes with right-handed orientation (RHM) on helix angle (HA) maps. Myocardial strain assessed using feature tracking-CMR (FT-CMR) is a sensitive marker of sub-clinical myocardial dysfunction. We sought to explore the relationship between these two techniques (strain and cDTI) in patients at 3 months following ST-elevation MI (STEMI).

METHODS

32 patients (F = 28, 60 ± 10 years) underwent 3T CMR three months after STEMI (mean interval 105 ± 17 days) with second order motion compensated (M2), free-breathing spin echo cDTI, cine gradient echo and late gadolinium enhancement (LGE) imaging. HA maps divided into left-handed HA (LHM, - 90 < HA < - 30), circumferential HA (CM, - 30° < HA < 30°), and right-handed HA (RHM, 30° < HA < 90°) were reported as relative proportions. Global and segmental analysis was undertaken.

RESULTS

Mean left ventricular ejection fraction (LVEF) was 44 ± 10% with a mean infarct size of 18 ± 12 g and a mean infarct segment LGE enhancement of 66 ± 21%. Mean global radial strain was 19 ± 6, mean global circumferential strain was - 13 ± - 3 and mean global longitudinal strain was - 10 ± - 3. Global and segmental radial strain correlated significantly with E2A in infarcted segments (p = 0.002, p = 0.011). Both global and segmental longitudinal strain correlated with RHM of infarcted segments on HA maps (p < 0.001, p = 0.003). Mean Diffusivity (MD) correlated significantly with the global infarct size (p < 0.008). When patients were categorised according to LVEF (reduced, mid-range and preserved), all cDTI parameters differed significantly between the three groups.

CONCLUSION

Change in sheetlet orientation assessed using E2A from cDTI correlates with impaired radial strain. Segments with fewer subendocardial cardiomyocytes, evidenced by a lower proportion of myocytes with right-handed orientation on HA maps, show impaired longitudinal strain. Infarct segment enhancement correlates significantly with E2A and RHM. Our data has demonstrated a link between myocardial microstructure and contractility following myocardial infarction, suggesting a potential role for CMR cDTI to clinically relevant functional impact.

Exercise-induced CITED4 expression is necessary for regional remodeling of cardiac microstructural tissue helicity

Both exercise-induced molecular mechanisms and physiological cardiac remodeling have been previously studied on a whole heart level. However, the regional microstructural tissue effects of these molecular mechanisms in the heart have yet to be spatially linked and further elucidated. We show in exercised mice that the expression of CITED4, a transcriptional co-regulator necessary for cardioprotection, is regionally heterogenous in the heart with preferential significant increases in the lateral wall compared with sedentary mice. Concordantly in this same region, the heart's local microstructural tissue helicity is also selectively increased in exercised mice. Quantification of CITED4 expression and microstructural tissue helicity reveals a significant correlation across both sedentary and exercise mouse cohorts. Furthermore, genetic deletion of CITED4 in the heart prohibits regional exercise-induced microstructural helicity remodeling. Taken together, CITED4 expression is necessary for exercise-induced regional remodeling of the heart's microstructural helicity revealing how a key molecular regulator of cardiac remodeling manifests into downstream local tissue-level changes.

In Vivo Super-Resolution Cardiac Diffusion Tensor MRI: A Feasibility Study

A super-resolution (SR) technique is proposed for imaging myocardial fiber architecture with cardiac magnetic resonance. Images were acquired with a motion-compensated cardiac diffusion tensor imaging (cDTI) sequence. The heart left ventricle was covered with three stacks of thick slices, in short axis, horizontal and vertical long axes orientations, respectively. The three low-resolution stacks (2 × 2 × 8 mm³) were combined into an isotropic volume (2 × 2 × 2 mm³) by a super-resolution reconstruction. For in vivo measurements, each slice was acquired during a breath-hold period. Bulk motion was corrected by optimizing a similarity metric between intensity profiles from all intersecting slices in the dataset. The benefit of the proposed approach was evaluated using a numerical heart phantom, a physical helicoidal phantom with artificial fibers, and six healthy subjects. The SR technique showed improved results compared to the native scans, in terms of image quality and cDTI metrics. In particular, the myocardial helix angle (HA) was more accurately estimated in the physical phantom (HA = 41.5° ± 1.1°, with the ground truth being 42.0°). In vivo, it resulted in a sharper rate of change of HA across the myocardial wall (−0.993°/% ± 0.007°/% against −0.873°/% ± 0.010°/%).

Accelerated free-breathing 3D whole-heart magnetic resonance angiography with a radial phyllotaxis trajectory, compressed sensing, and curvelet transform

PURPOSE

To develop and validate an accelerated free-breathing 3D whole-heart magnetic resonance angiography (MRA) technique using a radial k-space trajectory with compressed sensing and curvelet transform.

METHOD

A 3D radial phyllotaxis trajectory was implemented to traverse the centerline of k-space immediately before the segmented whole-heart MRA data acquisition at each cardiac cycle. The k-space centerlines were used to correct the respiratory-induced heart motion in the acquired MRA data. The corrected MRA data were then reconstructed by a novel compressed sensing algorithm using curvelets as the sparsifying domain. The proposed 3D whole-heart MRA technique (radial CS curvelet) was then prospectively validated against compressed sensing with a conventional wavelet transform (radial CS wavelet) and a standard Cartesian acquisition in terms of scan time and border sharpness.

RESULTS

Fifteen patients (females 10, median age 34-year-old) underwent 3D whole-heart MRA imaging using a standard Cartesian trajectory and our proposed radial phyllotaxis trajectory. Scan time for radial phyllotaxis was significantly shorter than Cartesian (4.88 ± 0.86 min. vs. 6.84 ± 1.79 min., P-value = 0.004). Radial CS curvelet border sharpness was slightly lower than Cartesian and, for the majority of vessels, was significantly better than radial CS wavelet (P-value < 0.050).

CONCLUSION

The proposed technique of 3D whole-heart MRA acquisition with a radial CS curvelet has a shorter scan time and slightly lower vessel sharpness compared to the Cartesian acquisition with radial profile ordering, and has slightly better sharpness than radial CS wavelet. Future work on this technique includes additional clinical trials and extending this technique to 3D cine imaging.

Optimized 64-channel array configurations for accelerated simultaneous multislice acquisitions in 3T cardiac MRI

PURPOSE

Three 64-channel cardiac coils with different detector array configurations were designed and constructed to evaluate acceleration capabilities in simultaneous multislice (SMS) imaging for 3T cardiac MRI.

METHODS

Three 64-channel coil array configurations obtained from a simulation-guided design approach were constructed and systematically evaluated regarding their encoding capabilities for accelerated SMS cardiac acquisitions at 3T. Array configuration A_Uni-sized consists of uniformly distributed equally sized loops in an overlapped arrangement, B_Gapped uses a gapped array design with symmetrically distributed equally sized loops, and C_Dense has non-uniform loop density and size, where smaller elements were centered over the heart and larger elements were placed surrounding the target region. To isolate the anatomic variation from differences in the coil configurations, all three array coils were built with identical semi-adjustable housing segments. The arrays' performance was compared using bench-level measurements and imaging performance tests, including signal-to-noise ratio (SNR) maps, array element noise correlation, and SMS acceleration capabilities. Additionally, all cardiac array coils were evaluated on a healthy volunteer.

RESULTS

The array configuration C_Dense with the non-uniformly distributed loop density showed the best overall cardiac imaging performance in both SNR and SMS encoding power, when compared to the other constructed arrays. The diffusion weighted cardiac acquisitions on a healthy volunteer support the favorable accelerated SNR performance of this array configuration.

CONCLUSION

Our results indicate that optimized highly parallel cardiac arrays, such as the 64-channel coil with a non-uniform loop size and density improve highly accelerated SMS cardiac MRI in comparison to symmetrically distributed loop array designs.

Free-breathing diffusion tensor MRI of the whole left ventricle using second-order motion compensation and multitasking respiratory motion correction

PURPOSE

We aimed to develop a novel free-breathing cardiac diffusion tensor MRI (DT-MRI) approach, M2-MT-MOCO, capable of whole left ventricular coverage that leverages second-order motion compensation (M2) diffusion encoding and multitasking (MT) framework to efficiently correct for respiratory motion (MOCO).

METHODS

Imaging was performed in 16 healthy volunteers and 3 heart failure patients with symptomatic dyspnea. The healthy volunteers were scanned to compare the accuracy of interleaved multislice coverage of the entire left ventricle with a single-slice acquisition and the accuracy of the free-breathing conventional MOCO and MT-MOCO approaches with reference breath-hold DT-MRI. Mean diffusivity (MD), fractional anisotropy (FA), helix angle transmurality (HAT), and intrascan repeatability were quantified and compared.

RESULTS

In all subjects, free-breathing M2-MT-MOCO DT-MRI yielded DWI of the entire left ventricle without bulk motion-induced signal loss. No significant differences were seen in the global values of MD, FA, and HAT in the multislice and single-slice acquisitions. Furthermore, global quantification of MD, FA, and HAT were also not significantly different between the MT-MOCO and breath-hold, whereas conventional MOCO yielded significant differences in MD, FA, and HAT with MT-MOCO and FA with breath-hold. In heart failure patients, M2-MT-MOCO DT-MRI was feasible yielding higher MD, lower FA, and lower HAT compared with healthy volunteers. Substantial agreement was found between repeated scans across all subjects for MT-MOCO.

CONCLUSION

M2-MT-MOCO enables free-breathing DT-MRI of the entire left ventricle in 10 min, while preserving quantification of myocardial microstructure compared to breath-held and single-slice acquisitions and is feasible in heart failure patients.

Myofiber strain in healthy humans using DENSE and cDTI

PURPOSE

Myofiber strain, E_ff , is a mechanistically relevant metric of cardiac cell shortening and is expected to be spatially uniform in healthy populations, making it a prime candidate for the evaluation of local cardiomyocyte contractility. In this study, a new, efficient pipeline was proposed to combine microstructural cDTI and functional DENSE data in order to estimate E_ff in vivo.

METHODS

Thirty healthy volunteers were scanned with three long-axis (LA) and three short-axis (SA) DENSE slices using 2D displacement encoding and one SA slice of cDTI. The total acquisition time was 11 minutes ± 3 minutes across volunteers. The pipeline first generates 3D SA displacements from all DENSE slices which are then combined with cDTI data to generate a cine of myofiber orientations and compute E_ff . The precision of the post-processing pipeline was assessed using a computational phantom study. Transmural myofiber strain was compared to circumferential strain, E_cc , in healthy volunteers using a Wilcoxon sign rank test.

RESULTS

In vivo, computed E_ff was found uniform transmurally compared to E_cc (-0.14[-0.15, -0.12] vs -0.18 [-0.20, -0.16], P < .001, -0.14 [-0.16, -0.12] vs -0.16 [-0.17, -0.13], P < .001 and -0.14 [-0.16, -0.12] vs E_{cc_C} = -0.14 [-0.15, -0.11], P = .002, E_{ff_C} vs E_{cc_C} in the endo, mid, and epi layers, respectively).

CONCLUSION

We demonstrate that it is possible to measure in vivo myofiber strain in a healthy human population in 10 minutes per subject. Myofiber strain was observed to be spatially uniform in healthy volunteers making it a potential biomarker for the evaluation of local cardiomyocyte contractility in assessing cardiovascular dysfunction.

Probing cardiomyocyte mobility with multi-phase cardiac diffusion tensor MRI

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<0.001). E2A also changed significantly between early systole, late systole, and diastasis (27.7 [20.8, 35.1]° vs 45.2 [42.1, 49]° vs 20.7 [16.6, 26.4]°, p<0.001).

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.

Magnetic Resonance-Based Characterization of Myocardial Architecture

Advances in technology have made it possible to image the microstructure of the heart with diffusion-weighted magnetic resonance. The technique provides unique insights into the cellular architecture of the myocardium and how this is perturbed in a range of disease contexts. In this review, the physical basis of diffusion MRI and the challenges of implementing it in the beating heart are discussed. Cutting edge acquisition and analysis techniques, as well as the results of initial clinical studies, are reported.

Longitudinal assessment of tissue properties and cardiac diffusion metrics of the ex vivo porcine heart at 7 T: Impact of continuous tissue fixation using formalin

In this study we aimed to assess the effects of continuous formalin fixation on diffusion and relaxation metrics of the ex vivo porcine heart at 7 T. Magnetic resonance imaging was performed on eight piglet hearts using a 7 T whole body system. Hearts were measured fresh within 3 hours of cardiac arrest followed by immersion in 10% neutral buffered formalin. T₂^* and T₂ were assessed using a gradient multi-echo and multi-echo spin echo sequence, respectively. A spin echo and a custom stimulated echo sequence were employed to assess diffusion time-dependent changes in metrics of cardiac diffusion tensor imaging. SNR was determined for b = 0 images. Scans were performed for 5 mm thick apical, midcavity and basal slices (in-plane resolution: 1 mm) and repeated 7, 15, 50, 100 and 200 days postfixation. Eigenvalues of the apparent diffusion coefficient (ADC) and fractional anisotropy (FA) decreased significantly (P < 0.05) following fixation. Relative to fresh hearts, FA values 7 and 200 days postfixation were 90% and 80%, while respective relative ADC values at those fixation stages were 78% and 92%. Statistical helix and sheetlet angle distributions as well as respective mean and median values showed no systematic influence of continuous formalin fixation. Similar to changes in the ADC, values for T₂ , T₂^* and SNR dropped initially postfixation. Respective relative values compared with fresh hearts at day 7 were 64%, 79% and 68%, whereas continuous fixation restored T₂ , T₂^* and SNR leading to relative values of 74%, 100%, and 81% at day 200, respectively. Relaxation parameters and diffusion metrics are significantly altered by continuous formalin fixation. The preservation of microstructure metrics following prolonged fixation is a key finding that may enable future studies of ventricular remodeling in cardiac pathologies.

Identification of Myocardial Disarray in Patients With Hypertrophic Cardiomyopathy and Ventricular Arrhythmias

Background

Myocardial disarray is a likely focus for fatal arrhythmia in hypertrophic cardiomyopathy (HCM). This microstructural abnormality can be inferred by mapping the preferential diffusion of water along cardiac muscle fibers using diffusion tensor cardiac magnetic resonance (DT-CMR) imaging. Fractional anisotropy (FA) quantifies directionality of diffusion in 3 dimensions. The authors hypothesized that FA would be reduced in HCM due to disarray and fibrosis that may represent the anatomic substrate for ventricular arrhythmia.

Objectives

This study sought to assess FA as a noninvasive in vivo biomarker of HCM myoarchitecture and its association with ventricular arrhythmia.

Methods

A total of 50 HCM patients (47 ± 15 years of age, 77% male) and 30 healthy control subjects (46 ± 16 years of age, 70% male) underwent DT-CMR in diastole, cine, late gadolinium enhancement (LGE), and extracellular volume (ECV) imaging at 3-T.

Results

Diastolic FA was reduced in HCM compared with control subjects (0.49 ± 0.05 vs. 0.52 ± 0.03; p = 0.0005). Control subjects had a mid-wall ring of high FA. In HCM, this ring was disrupted by reduced FA, consistent with published histology demonstrating that disarray and fibrosis invade circumferentially aligned mid-wall myocytes. LGE and ECV were significant predictors of FA, in line with fibrosis contributing to low FA. Yet FA adjusted for LGE and ECV remained reduced in HCM (p = 0.028). FA in the hypertrophied segment was reduced in HCM patients with ventricular arrhythmia compared to patients without (n = 15; 0.41 ± 0.03 vs. 0.46 ± 0.06; p = 0.007). A decrease in FA of 0.05 increased odds of ventricular arrhythmia by 2.5 (95% confidence interval: 1.2 to 5.3; p = 0.015) in HCM and remained significant even after correcting for LGE, ECV, and wall thickness (p = 0.036).

Conclusions

DT-CMR assessment of left ventricular myoarchitecture matched patterns reported previously on histology. Low diastolic FA in HCM was associated with ventricular arrhythmia and is likely to represent disarray after accounting for fibrosis. The authors propose that diastolic FA could be the first in vivo marker of disarray in HCM and a potential independent risk factor.

Improved compressed sensing and super-resolution of cardiac diffusion MRI with structure-guided total variation

Purpose

Structure‐guided total variation is a recently introduced prior that allows reconstruction of images using knowledge of the location and orientation of edges in a reference image. In this work, we demonstrate the advantages of a variant of structure‐guided total variation known as directional total variation (DTV), over traditional total variation (TV), in the context of compressed‐sensing reconstruction and super‐resolution.

Methods

We compared TV and DTV in retrospectively undersampled ex vivo diffusion tensor imaging and diffusion spectrum imaging data from healthy, sham, and hypertrophic rat hearts.

Results

In compressed sensing at an undersampling factor of 8, the RMS error of mean diffusivity and fractional anisotropy relative to the fully sampled ground truth were 44% and 20% lower in DTV compared with TV. In super‐resolution, these values were 29% and 14%, respectively. Similarly, we observed improvements in helix angle, transverse angle, sheetlet elevation, and sheetlet azimuth. The RMS error of the diffusion kurtosis in the undersampled data relative to the ground truth was uniformly lower (22% on average) with DTV compared to TV.

Conclusion

Acquiring one fully sampled non‐diffusion‐weighted image and 10 diffusion‐weighted images at 8× undersampling would result in an 80% net reduction in data needed. We demonstrate efficacy of the DTV algorithm over TV in reducing data sampling requirements, which can be translated into higher apparent resolution and potentially shorter scan times. This method would be equally applicable in diffusion MRI applications outside the heart.

Motion and eddy current-induced signal dephasing in in vivo cardiac DTI

PURPOSE

To address motion in cardiac DWI, stimulated-echo acquisition mode (STEAM) and second-order motion-compensated spin-echo (SE) sequences have been proposed. Despite applying motion-compensation strategies, residual motion can cause misleading signal attenuation. The purpose of this study is to estimate the motion-induced error in both sequences by analysis of image phase.

METHODS

Diffusion-weighted motion-compensated SE sequences and STEAM imaging was applied in vivo with diffusion encoding along 3 orthogonal directions. A b-value range of 100 to 600 s/mm² and trigger delays of 25%, 50%, and 75% of end systole and middiastole were used. Eddy-current contributions were obtained from phantom measurements. After computation of motion-induced phase maps, the amount of signal dephasing was computed from phase gradients, and the resulting errors in diffusion tensor parameters were calculated.

RESULTS

Motion-induced dephasing from the STEAM sequence showed less dependency on the b-value and no dependency on the heart phase, whereas SE imaging performed best at 75% end systole followed by 50% end systole and middiastole. For a typical experimental setting, errors of 3.3%/3.0% mean diffusivity, 4.9%/4.8% fractional anisotropy, 2.9º/3.2º helix angulation, 0.8º/0.7º transverse angulation, and 9.9º/10.0º sheet angulation (SE/STEAM) were calculated.

CONCLUSION

Image phase contains valuable information regarding uncompensated motion and eddy currents in cardiac DTI. Although the trigger delay window for SE is narrower compared with the STEAM-based approach, imaging in both systole and diastole is feasible and both sequences perform similarly if the trigger delays are selected carefully with SE.

Cardiac Diffusion: Technique and Practical Applications

The 3D microarchitecture of the cardiac muscle underlies the mechanical and electrical properties of the heart. Cardiomyocytes are arranged helically through the depth of the wall, and their shortening leads to macroscopic torsion, twist, and shortening during cardiac contraction. Furthermore, cardiomyocytes are organized in sheetlets separated by shear layers, which reorientate, slip, and shear during macroscopic left ventricle (LV) wall thickening. Cardiac diffusion provides a means for noninvasive interrogation of the 3D microarchitecture of the myocardium. The fundamental principle of MR diffusion is that an MRI signal is attenuated by the self-diffusion of water in the presence of large diffusion-encoding gradients. Since water molecules are constrained by the boundaries in biological tissue (cell membranes, collagen layers, etc.), depicting their diffusion behavior elucidates the shape of the myocardial microarchitecture they are embedded in. Cardiac diffusion therefore provides a noninvasive means to understand not only the dynamic changes in cardiac microstructure of healthy myocardium during cardiac contraction but also the pathophysiological changes in the presence of disease. This unique and innovative technology offers tremendous potential to enable improved clinical diagnosis through novel microstructural and functional assessment. in vivo cardiac diffusion methods are immediately translatable to patients, opening new avenues for diagnostic investigation and treatment evaluation in a range of clinically important cardiac pathologies. This review article describes the 3D microstructure of the LV, explains in vivo and ex vivo cardiac MR diffusion acquisition and postprocessing techniques, as well as clinical applications to date. Level of Evidence: 1 Technical Efficacy: Stage 3 J. Magn. Reson. Imaging 2019. J. Magn. Reson. Imaging 2020;52:348-368.

Spin echo based cardiac diffusion imaging at 7T: An ex vivo study of the porcine heart at 7T and 3T

Purpose of this work was to assess feasibility of cardiac diffusion tensor imaging (cDTI) at 7 T in a set of healthy, unfixed, porcine hearts using various parallel imaging acceleration factors and to compare SNR and derived cDTI metrics to a reference measured at 3 T. Magnetic resonance imaging was performed on 7T and 3T whole body systems using a spin echo diffusion encoding sequence with echo planar imaging readout. Five reference (b = 0 s/mm2) images and 30 diffusion directions (b = 700 s/mm2) were acquired at both 7 T and 3 T using a GRAPPA acceleration factor R = 1. Scans at 7 T were repeated using R = 2, R = 3, and R = 4. SNR evaluation was based on 30 reference (b = 0 s/mm2) images of 30 slices of the left ventricle and cardiac DTI metrics were compared within AHA segmentation. The number of hearts scanned at 7 T and 3 T was n = 11. No statistically significant differences were found for evaluated helix angle, secondary eigenvector angle, fractional anisotropy and apparent diffusion coefficient at the different field strengths, given sufficiently high SNR and geometrically undistorted images. R≥3 was needed to reduce susceptibility induced geometric distortions to an acceptable amount. On average SNR in myocardium of the left ventricle was increased from 29±3 to 44±6 in the reference image (b = 0 s/mm2) when switching from 3 T to 7 T. Our study demonstrates that high resolution, ex vivo cDTI is feasible at 7 T using commercial hardware.

An Overview on Image Registration Techniques for Cardiac Diagnosis and Treatment

Image registration has been used for a wide variety of tasks within cardiovascular imaging. This study aims to provide an overview of the existing image registration methods to assist researchers and impart valuable resource for studying the existing methods or developing new methods and evaluation strategies for cardiac image registration. For the cardiac diagnosis and treatment strategy, image registration and fusion can provide complementary information to the physician by using the integrated image from these two modalities. This review also contains a description of various imaging techniques to provide an appreciation of the problems associated with implementing image registration, particularly for cardiac pathology intervention and treatments.

Myocardial Scar Delineation Using Diffusion Tensor Magnetic Resonance Tractography

Background

Late gadolinium enhancement (LGE) is the current standard for myocardial scar delineation. In this study, we introduce the tractographic propagation angle (PA), a metric of myofiber curvature (degrees/unit distance) derived from diffusion tensor imaging (DTI), and compare its use to LGE and invasive scar assessment by endocardial voltage mapping.

Methods and Results

DTI was performed on 7 healthy human volunteers, 5 patients with myocardial infarction, 6 normal mice, and 7 mice with myocardial infarction. LGE to delineate the infarct and border zones was performed with a 2‐dimensional inversion recovery gradient‐echo sequence. Ex vivo DTI was performed on 5 normal human and 5 normal sheep hearts. Endocardial electroanatomic mapping and subsequent ex vivo DTI was performed on 5 infarcted sheep hearts. PA in the normal human hearts varied smoothly and was generally <4. The mean PA in the infarct zone was significantly elevated (10.34±1.02 versus 4.05±0.45, P<0.05). Regions with a PA ≤4 consistently had a bipolar voltage ≥1.5 mV, whereas those with PA values between 4 and 10 had voltages between 0.5 and 1.5 mV. A PA threshold >4 was the most accurate DTI‐derived measure of infarct size and demonstrated the greatest correlation with LGE (r=0.95).

Conclusions

We found a strong correlation between infarct size by PA and LGE in both mice and humans. There was also an inverse relationship between PA values and endocardial voltage. The use of PA may enable myocardial scar delineation and characterization of arrhythmogenic substrate without the need for exogenous contrast agents.

Changes in overall ventricular myocardial architecture in the setting of a porcine animal model of right ventricular dilation

BACKGROUND

Chronic pulmonary regurgitation often leads to myocardial dysfunction and heart failure. It is not fully known why secondary hypertrophy cannot fully protect against the increase in wall stress brought about by the increased end-diastolic volume in ventricular dilation. It has been assumed that mural architecture is not deranged in this situation, but we hypothesised that there might be a change in the pattern of orientation of the aggregations of cardiomyocytes, which would contribute to contractile impairment.

METHODS

We created pulmonary valvular regurgitation by open chest, surgical suturing of its leaflets in seven piglets, performing sham operations in seven control animals. Using cardiovascular magnetic resonance imaging after 12 weeks of recovery, we demonstrated significantly increased right ventricular volumes in the test group. After sacrifice, diffusion tensor imaging of their hearts permitted measurement of the orientation of the cardiomyocytes.

RESULTS

The helical angles in the right ventricle approached a more circumferential orientation in the setting of right ventricular RV dilation (p = 0.007), with an increased proportion of surface-parallel cardiomyocytes. In contrast, this proportion decreased in the left ventricle. Also in the left ventricle a higher proportion of E3 angles with a value around zero was found, and conversely a lower proportion of angles was found with a numerical higher value. In the dilated right ventricle the proportion of E3 angles around -90° is increased, while the proportion around 90° is decreased.

CONCLUSION

Contrary to traditional views, there is a change in the orientation of both the left ventricular and right ventricular cardiomyocytes subsequent to right ventricular dilation. This will change their direction of contraction and hinder the achievement of normalisation of cardiomyocytic strain, affecting overall contractility. We suggest that the aetiology of the cardiac failure induced by right vetricular dilation may be partly explained by morphological changes in the myocardium itself.

When should we use contrast material in cardiac MRI?

At present, most of the cardiac magnetic resonance imaging (MRI) examinations rely on contrast-enhanced protocols, but noncontrast alternatives are emerging. Late gadolinium enhancement (LGE) imaging for the detection of myocardial scar can be considered the main cause for the embedding of cardiac MRI into the clinical routine. The novel noncontrast technique of native T₁ mapping shows promise for tissue characterization in ischemic and nonischemic cardiomyopathy and may provide additional information over conventional LGE imaging. Technical issues, including measurements variability, still need to be resolved to facilitate a wide clinical application. Ischemia detection can be performed with contrast-based stress perfusion and contrast-free stress wall motion imaging. For coronary magnetic resonance angiography (MRA), protocols with and without contrast material have been developed. Research on coronary atherosclerotic plaque characterization has introduced new applications of contrast material. For MRA of the aorta, which traditionally relied on contrast administration, several noncontrast protocols have become available. This review provides an overview of when to use contrast material in cardiac and cardiac-related vascular MRI, summarizes the major imaging building blocks, and describes the diagnostic value of the available contrast-enhanced and noncontrast techniques. Contrast material in cardiac MRI should be used for LGE imaging for tissue characterization in ischemic or nonischemic cardiomyopathy and may be used for stress perfusion imaging for the detection of ischemia. In cardiac-related vascular MRI, use of contrast material should be avoided, unless high-quality angiography is required that cannot be obtained with noncontrast protocols.

LEVEL OF EVIDENCE

5 Technical Efficacy: Stage 3 J. Magn. Reson. Imaging 2017;46:1551-1572.