Studying Dynamic Myofiber Aggregate Reorientation in Dilated Cardiomyopathy Using In Vivo Magnetic Resonance Diffusion Tensor Imaging

CMR to characterize myocardial structure and function in heart failure with preserved left ventricular ejection fraction

Despite remarkable progress in therapeutic drugs, morbidity, and mortality for heart failure (HF) remains high in developed countries. HF with preserved ejection fraction (HFpEF) now accounts for around half of all HF cases. It is a heterogeneous disease, with multiple aetiologies, and as such poses a significant diagnostic challenge. Cardiac magnetic resonance (CMR) has become a valuable non-invasive modality to assess cardiac morphology and function, but beyond that, the multi-parametric nature of CMR allows novel approaches to characterize haemodynamics and with magnetic resonance spectroscopy (MRS), the study of metabolism. Furthermore, exercise CMR, when combined with lung water imaging provides an in-depth understanding of the underlying pathophysiological and mechanistic processes in HFpEF. Thus, CMR provides a comprehensive phenotyping tool for HFpEF, which points towards a targeted and personalized therapy with improved diagnostics and prevention.

Cardiac diffusion-weighted and tensor imaging: a Society for Cardiovascular Magnetic Resonance (SCMR) special interest group consensus statement

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 3D microstructure of the heart. Free-contrast DWI is uniquely sensitized to the presence and displacement of water molecules within the myocardial tissue, including the intra-cellular, extra-cellular and intra-vascular 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) 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, post-processing 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.

Second-Order Motion-Compensated Echo-Planar Cardiac Diffusion-Weighted MRI: Usefulness of Compressed Sensitivity Encoding

BACKGROUND

Cardiac diffusion-weighted imaging (DWI) using second-order motion-compensated spin echo (M2C) can provide noninvasive in-vivo microstructural assessment, but limited by relatively low signal-to-noise ratio (SNR). Echo-planar imaging (EPI) with compressed sensitivity encoding (EPICS) could address these issues.

PURPOSE

To combine M2C DWI and EPCIS (M2C EPICS DWI), and compare image quality for M2C DWI.

STUDY TYPE

Prospective.

POPULATION

Ten ex-vivo hearts, 10 healthy volunteers (females, 5 [50%]; mean ± SD of age, 25 ± 4 years), and 12 patients with diseased hearts (female, 1 [8.3%]; mean ± SD of age, 44 ± 16 years; including coronary artery heart disease, congenital heart disease, dilated cardiomyopathy, amyloidosis, and myocarditis).

FIELD STRENGTH/SEQUENCE

3-T, M2C EPICS DWI, and M2C DWI.

ASSESSMENT

The apparent SNR (aSNR) and the rating scores were used to evaluate and compared image quality of all three groups. The aSNR was calculated usingaSNR = Mean intensity myocardium / Standard deviation myocardium $$ \mathrm{aSNR}={\mathrm{Mean}\ \mathrm{intensity}}_{\mathrm{myocardium}}/{\mathrm{Standard}\ \mathrm{deviation}}_{\mathrm{myocardium}} $$ , and the myocardium was segmented manually. Three observers independently rated subjective image quality using a 5-point Likert scale.

STATISTICAL TESTS

Bland-Altman analysis and paired t-tests. The threshold for statistical significance was set at P < 0.05.

RESULTS

In healthy volunteers, the aSNR with a b-value of 450 s/mm2 acquired by M2C EPICS DWI was significantly higher than M2C DWI at in-plane resolutions of 3.0 × 3.0, 2.5 × 2.5, and 2.0 × 2.0 mm2. In patients with diseased hearts, the aSNR ofM2C EPICS DWI was also significantly higher than that for M2C DWI (bias of M2C EPICS-M2C = 1.999, 95% limits of agreement, 0.362 to 3.636; mean ± SD, 7.80 ± 1.37 vs. 5.80 ± 0.81). The ADC values of M2C EPICS was significantly higher than M2C DWI in in-vivo hearts. Over 80% of the images with rating scores for M2C EPICS DWI were higher than M2C DWI in in-vivo hearts.

DATA CONCLUSION

Cardiac imaging by M2C EPICS DWI may demonstrate better overall image quality and higher aSNR than M2C DWI.

EVIDENCE LEVEL

2 TECHNICAL EFFICACY: Stage 1.

Development and validation of cardiac diffusion weighted magnetic resonance imaging for the diagnosis of myocardial injury in small animal models

Cardiac diffusion weighted-magnetic resonance imaging (DWI) has slowly developed due to its technical difficulties. However, this limitation could be overcome by advanced techniques, including a stimulated echo technique and a gradient moment nulling technique. This study aimed to develop and validate a high-order DWI sequence, using echo-planar imaging (EPI) and second-order motion-compensated (M012) diffusion gradient applied to cardiac imaging in small-sized animals with fast heart and respiratory rates, and to investigate the feasibility of cardiac DWI, diagnosing acute myocardial injury in isoproterenol-induced myocardial injury rat models. The M012 diffusion gradient sequence was designed for diffusion tensor imaging of the rat myocardium and validated in the polyvinylpyrrolidone phantom. Following sequence optimization, 23 rats with isoproterenol-induced acute myocardial injury and five healthy control rats underwent cardiac MRI, including cine imaging, T1 mapping, and DWI. Diffusion gradient was applied using a 9.4-T MRI scanner (Bruker, BioSpec 94/20, gradient amplitude = 440 mT/m, maximum slew rate = 3440 T/m/s) with double gating (electrocardiogram and respiratory gating). Troponin I was used as a serum biomarker for myocardial injury. Histopathologic examination of the heart was subsequently performed. The developed DWI sequence using EPI and M012 provided the interpretable images of rat hearts. The apparent diffusion coefficient (ADC) values were significantly higher in rats with acute myocardial injury than in the control group (1.847 ± 0.326 * 10-3 mm2/s vs. 1.578 ± 0.144 * 10-3 mm2/s, P < 0.001). Troponin I levels were increased in the blood samples of rats with acute myocardial injury (P < 0.001). Histopathologic examinations detected myocardial damage and subendocardial fibrosis in rats with acute myocardial injury. The newly developed DWI technique has the ability to detect myocardial injury in small animal models, representing high ADC values on the myocardium with isoproterenol-induced injury.

Advancing clinical translation of cardiac biomechanics models: a comprehensive review, applications and future pathways

Cardiac mechanics models are developed to represent a high level of detail, including refined anatomies, accurate cell mechanics models, and platforms to link microscale physiology to whole-organ function. However, cardiac biomechanics models still have limited clinical translation. In this review, we provide a picture of cardiac mechanics models, focusing on their clinical translation. We review the main experimental and clinical data used in cardiac models, as well as the steps followed in the literature to generate anatomical meshes ready for simulations. We describe the main models in active and passive mechanics and the different lumped parameter models to represent the circulatory system. Lastly, we provide a summary of the state-of-the-art in terms of ventricular, atrial, and four-chamber cardiac biomechanics models. We discuss the steps that may facilitate clinical translation of the biomechanics models we describe. A well-established software to simulate cardiac biomechanics is lacking, with all available platforms involving different levels of documentation, learning curves, accessibility, and cost. Furthermore, there is no regulatory framework that clearly outlines the verification and validation requirements a model has to satisfy in order to be reliably used in applications. Finally, better integration with increasingly rich clinical and/or experimental datasets as well as machine learning techniques to reduce computational costs might increase model reliability at feasible resources. Cardiac biomechanics models provide excellent opportunities to be integrated into clinical workflows, but more refinement and careful validation against clinical data are needed to improve their credibility. In addition, in each context of use, model complexity must be balanced with the associated high computational cost of running these models.

Cardiomyocyte orientation recovery at micrometer scale reveals long-axis fiber continuum in heart walls

Coordinated cardiomyocyte contraction drives the mammalian heart to beat and circulate blood. No consensus model of cardiomyocyte geometrical arrangement exists, due to the limited spatial resolution of whole heart imaging methods and the piecemeal nature of studies based on histological sections. By combining microscopy and computer vision, we produced the first-ever three-dimensional cardiomyocyte orientation reconstruction across mouse ventricular walls at the micrometer scale, representing a gain of three orders of magnitude in spatial resolution. We recovered a cardiomyocyte arrangement aligned to the long-axis direction of the outer ventricular walls. This cellular network lies in a thin shell and forms a continuum with longitudinally arranged cardiomyocytes in the inner walls, with a complex geometry at the apex. Our reconstruction methods can be applied at fine spatial scales to further understanding of heart wall electrical function and mechanics, and set the stage for the study of micron-scale fiber remodeling in heart disease.

Phenotyping heart failure by cardiac magnetic resonance imaging of cardiac macro- and microscopic structure: state of the art review

Heart failure demographics have evolved in past decades with the development of improved diagnostics, therapies, and prevention. Cardiac magnetic resonance (CMR) has developed in a similar timeframe to become the gold-standard non-invasive imaging modality for characterizing diseases causing heart failure. CMR techniques to assess cardiac morphology and function have progressed since their first use in the 1980s. Increasingly efficient acquisition protocols generate high spatial and temporal resolution images in less time. This has enabled new methods of characterizing cardiac systolic and diastolic function such as strain analysis, exercise real-time cine imaging and four-dimensional flow. A key strength of CMR is its ability to non-invasively interrogate the myocardial tissue composition. Gadolinium contrast agents revolutionized non-invasive cardiac imaging with the late gadolinium enhancement technique. Further advances enabled quantitative parametric mapping to increase sensitivity at detecting diffuse pathology. Novel methods such as diffusion tensor imaging and artificial intelligence-enhanced image generation are on the horizon. Magnetic resonance spectroscopy (MRS) provides a window into the molecular environment of the myocardium. Phosphorus (31P) spectroscopy can inform the status of cardiac energetics in health and disease. Proton (1H) spectroscopy complements this by measuring creatine and intramyocardial lipids. Hyperpolarized carbon (13C) spectroscopy is a novel method that could further our understanding of dynamic cardiac metabolism. CMR of other organs such as the lungs may add further depth into phenotypes of heart failure. The vast capabilities of CMR should be deployed and interpreted in context of current heart failure challenges.

A tomographic microscopy-compatible Langendorff system for the dynamic structural characterization of the cardiac cycle

Introduction

Cardiac architecture has been extensively investigated ex vivo using a broad spectrum of imaging techniques. Nevertheless, the heart is a dynamic system and the structural mechanisms governing the cardiac cycle can only be unveiled when investigating it as such.

Methods

This work presents the customization of an isolated, perfused heart system compatible with synchrotron-based X-ray phase contrast imaging (X-PCI).

Results

Thanks to the capabilities of the developed setup, it was possible to visualize a beating isolated, perfused rat heart for the very first time in 4D at an unprecedented 2.75 μm pixel size (10.6 μm spatial resolution), and 1 ms temporal resolution.

Discussion

The customized setup allows high-spatial resolution studies of heart architecture along the cardiac cycle and has thus the potential to serve as a tool for the characterization of the structural dynamics of the heart, including the effects of drugs and other substances able to modify the cardiac cycle.

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.

Characterization and correction of diffusion gradient-induced eddy currents in second-order motion-compensated echo-planar and spiral cardiac DTI

PURPOSE

Very high gradient amplitudes played out over extended time intervals as required for second-order motion-compensated cardiac DTI may violate the assumption of a linear time-invariant gradient system model. The aim of this work was to characterize diffusion gradient-related system nonlinearity and propose a correction approach for echo-planar and spiral spin-echo motion-compensated cardiac DTI.

METHODS

Diffusion gradient-induced eddy currents of 9 diffusion directions were characterized at b values of 150 s/mm² and 450 s/mm² for a 1.5 Tesla system and used to correct phantom, ex vivo, and in vivo motion-compensated cardiac DTI data acquired with echo-planar and spiral trajectories. Predicted trajectories were calculated using gradient impulse response function and diffusion gradient strength- and direction-dependent zeroth- and first-order eddy current responses. A reconstruction method was implemented using the predicted <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mi>k</mml:mi></mml:mrow> <mml:annotation>$$ k $$</mml:annotation></mml:semantics> </mml:math> -space trajectories to additionally include off-resonances and concomitant fields. Resulting images were compared to a reference reconstruction omitting diffusion gradient-induced eddy current correction.

RESULTS

Diffusion gradient-induced eddy currents exhibited nonlinear effects when scaling up the gradient amplitude and could not be described by a 3D basis alone. This indicates that a gradient impulse response function does not suffice to describe diffusion gradient-induced eddy currents. Zeroth- and first-order diffusion gradient-induced eddy current effects of up to -1.7 rad and -16 to +12 rad/m, respectively, were identified. Zeroth- and first-order diffusion gradient-induced eddy current correction yielded improved image quality upon image reconstruction.

CONCLUSION

The proposed approach offers correction of diffusion gradient-induced zeroth- and first-order eddy currents, reducing image distortions to promote improvements of second-order motion-compensated spin-echo cardiac DTI.

Myocardial mesostructure and mesofunction

The complex and highly organized structural arrangement of some five billion cardiomyocytes directs the coordinated electrical activity and mechanical contraction of the human heart. The characteristic transmural change in cardiomyocyte orientation underlies base-to-apex shortening, circumferential shortening, and left ventricular torsion during contraction. Individual cardiomyocytes shorten ∼15% and increase in diameter ∼8%. Remarkably, however, the left ventricular wall thickens by up to 30-40%. To accommodate this, the myocardium must undergo significant structural rearrangement during contraction. At the mesoscale, collections of cardiomyocytes are organized into sheetlets, and sheetlet shear is the fundamental mechanism of rearrangement that produces wall thickening. Herein, we review the histological and physiological studies of myocardial mesostructure that have established the sheetlet shear model of wall thickening. Recent developments in tissue clearing techniques allow for imaging of whole hearts at the cellular scale, whereas magnetic resonance imaging (MRI) and computed tomography (CT) can image the myocardium at the mesoscale (100 µm to 1 mm) to resolve cardiomyocyte orientation and organization. Through histology, cardiac diffusion tensor imaging (DTI), and other modalities, mesostructural sheetlets have been confirmed in both animal and human hearts. Recent in vivo cardiac DTI methods have measured reorientation of sheetlets during the cardiac cycle. We also examine the role of pathological cardiac remodeling on sheetlet organization and reorientation, and the impact this has on ventricular function and dysfunction. We also review the unresolved mesostructural questions and challenges that may direct future work in the field.

Evaluation of Right Ventricular Function and Myocardial Microstructure in Fetal Hypoplastic Left Heart Syndrome

Right ventricular (RV) function is one of the critical factors affecting the prognosis of fetuses with hypoplastic left heart syndrome (HLHS). Our study objectives included assessment of cardiac function and comprehensive measurement of cardiac microstructure. We retrospectively studied 42 fetuses diagnosed as HLHS by echocardiography. Myocardial deformation of the right ventricular wall was calculated automatically in offline software. Postmortem cardiac imaging for three control fetal hearts and four HLHS specimens was performed by a 9.4T DTI scanner. Myocardial deformation parameters of the RV (including strain, strain rate, and velocity) were significantly lower in HLHS fetuses (all p < 0.01). FA values increased (0.18 ± 0.01 vs. 0.21 ± 0.02; p < 0.01) in HLHS fetuses, but MD reduced (1.3 ± 0.15 vs. 0.88 ± 0.13; p < 0.001). The HLHS fetuses’ RV lateral base wall (−7.31 ± 51.91 vs. −6.85 ± 31.34; p = 0.25), middle wall (1.71 ± 50.92 vs. −9.38 ± 28.18; p < 0.001), and apical wall (−6.19 ± 46.61 vs. −11.16 ± 29.86, p < 0.001) had HA gradient ascent but HA gradient descent in the anteroseptal wall (p < 0.001) and inferoseptal wall (p < 0.001). RV basal lateral wall HA degrees were correlated with RVGLS (R2 = 0.97, p = 0.02). MD values were positively correlated with RVGLS (R2 = 0.93, p = 0.04). Our study found morphological and functional changes of the RV in HLHS fetuses, and cardiac function was related to the orientation patterns of myocardial fibers. It may provide insight into understanding the underlying mechanisms of impaired RV performance in HLHS.

A groupwise registration and tractography framework for cardiac myofiber architecture description by diffusion MRI: An application to the ventricular junctions

Background

Knowledge of the normal myocardial–myocyte orientation could theoretically allow the definition of relevant quantitative biomarkers in clinical routine to diagnose heart pathologies. A whole heart diffusion tensor template representative of the global myofiber organization over species is therefore crucial for comparisons across populations. In this study, we developed a groupwise registration and tractography framework to resolve the global myofiber arrangement of large mammalian sheep hearts. To demonstrate the potential application of the proposed method, a novel description of sub-regions in the intraventricular septum is presented.

Methods

Three explanted sheep (ovine) hearts (size ~12×8×6 cm3, heart weight ~ 150 g) were perfused with contrast agent and fixative and imaged in a 9.4T magnet. A group-wise registration of high-resolution anatomical and diffusion-weighted images were performed to generate anatomical and diffusion tensor templates. Diffusion tensor metrics (eigenvalues, eigenvectors, fractional anisotropy …) were computed to provide a quantitative and spatially-resolved analysis of cardiac microstructure. Then tractography was performed using deterministic and probabilistic algorithms and used for different purposes: i) Visualization of myofiber architecture, ii) Segmentation of sub-area depicting the same fiber organization, iii) Seeding and Tract Editing. Finally, dissection was performed to confirm the existence of macroscopic structures identified in the diffusion tensor template.

Results

The template creation takes advantage of high-resolution anatomical and diffusion-weighted images obtained at an isotropic resolution of 150 μm and 600 μm respectively, covering ventricles and atria and providing information on the normal myocardial architecture. The diffusion metric distributions from the template were found close to the one of the individual samples validating the registration procedure. Small new sub-regions exhibiting spatially sharp variations in fiber orientation close to the junctions of the septum and ventricles were identified. Each substructure was defined and represented using streamlines. The existence of a fiber-bundles in the posterior junction was validated by anatomical dissection. A complex structural organization of the anterior junction in comparison to the posterior junction was evidenced by the high-resolution acquisition.

Conclusions

A new framework combining cardiac template generation and tractography was applied on the whole sheep heart. The framework can be used for anatomical investigation, characterization of microstructure and visualization of myofiber orientation across samples. Finally, a novel description of the ventricular junction in large mammalian sheep hearts was proposed.

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.

Development of a cardiovascular magnetic resonance-compatible large animal isolated heart model for direct comparison of beating and arrested hearts

Cardiac motion results in image artefacts and quantification errors in many cardiovascular magnetic resonance (CMR) techniques, including microstructural assessment using diffusion tensor cardiovascular magnetic resonance (DT-CMR). Here, we develop a CMR-compatible isolated perfused porcine heart model that allows comparison of data obtained in beating and arrested states. Ten porcine hearts (8/10 for protocol optimisation) were harvested using a donor heart retrieval protocol and transported to the remote CMR facility. Langendorff perfusion in a 3D-printed chamber and perfusion circuit re-established contraction. Hearts were imaged using cine, parametric mapping and STEAM DT-CMR at cardiac phases with the minimum and maximum wall thickness. High potassium and lithium perfusates were then used to arrest the heart in a slack and contracted state, respectively. Imaging was repeated in both arrested states. After imaging, tissue was removed for subsequent histology in a location matched to the DT-CMR data using fiducial markers. Regular sustained contraction was successfully established in six out of 10 hearts, including the final five hearts. Imaging was performed in four hearts and one underwent the full protocol, including colocalised histology. The image quality was good and there was good agreement between DT-CMR data in equivalent beating and arrested states. Despite the use of autologous blood and dextran within the perfusate, T2 mapping results, DT-CMR measures and an increase in mass were consistent with development of myocardial oedema, resulting in failure to achieve a true diastolic-like state. A contiguous stack of 313 5-μm histological sections at and a 100-μm thick section showing cell morphology on 3D fluorescent confocal microscopy colocalised to DT-CMR data were obtained. A CMR-compatible isolated perfused beating heart setup for large animal hearts allows direct comparisons of beating and arrested heart data with subsequent colocalised histology, without the need for onsite preclinical facilities.

WA, Boyle J, Coll‐Font J, Dall'Armellina E, Ennis DB, Ferreira PF, Kalra P, Kolipaka A, Kozerke S, Lohr D, Mongeon F, Moulin K, Nguyen C, Nielles‐Vallespin S, Raterman B, Schreiber LM, Scott AD, Sosnovik DE, Stoeck CT, Tous C, Tunnicliffe EM, Weng AM, Croisille P, Viallon M, Schneider JE. Validation of cardiac diffusion tensor imaging sequences: A multicentre test-retest phantom study

Cardiac diffusion tensor imaging (DTI) is an emerging technique for the in vivo characterisation of myocardial microstructure, and there is a growing need for its validation and standardisation. We sought to establish the accuracy, precision, repeatability and reproducibility of state-of-the-art pulse sequences for cardiac DTI among 10 centres internationally. Phantoms comprising 0%-20% polyvinylpyrrolidone (PVP) were scanned with DTI using a product pulsed gradient spin echo (PGSE; N = 10 sites) sequence, and a custom motion-compensated spin echo (SE; N = 5) or stimulated echo acquisition mode (STEAM; N = 5) sequence suitable for cardiac DTI in vivo. A second identical scan was performed 1-9 days later, and the data were analysed centrally. The average mean diffusivities (MDs) in 0% PVP were (1.124, 1.130, 1.113) x 10^-3  mm² /s for PGSE, SE and STEAM, respectively, and accurate to within 1.5% of reference data from the literature. The coefficients of variation in MDs across sites were 2.6%, 3.1% and 2.1% for PGSE, SE and STEAM, respectively, and were similar to previous studies using only PGSE. Reproducibility in MD was excellent, with mean differences in PGSE, SE and STEAM of (0.3 ± 2.3, 0.24 ± 0.95, 0.52 ± 0.58) x 10^-5  mm² /s (mean ± 1.96 SD). We show that custom sequences for cardiac DTI provide accurate, precise, repeatable and reproducible measurements. Further work in anisotropic and/or deforming phantoms is warranted.

Impact of Intraventricular Septal Fiber Orientation on Cardiac Electromechanical Function

Cardiac fiber direction is an important factor determining the propagation of electrical activity, as well as the development of mechanical force. In this article, we imaged the ventricles of several species with special attention to the intraventricular septum to determine the functional consequences of septal fiber organization. First, we identified a dual-layer organization of the fiber orientation in the intraventricular septum of ex vivo sheep hearts using diffusion tensor imaging at high field MRI. To expand the scope of the results, we investigated the presence of a similar fiber organization in five mammalian species (rat, canine, pig, sheep, and human) and highlighted the continuity of the layer with the moderator band in large mammalian species. We implemented the measured septal fiber fields in three-dimensional electromechanical computer models to assess the impact of the fiber orientation. The downward fibers produced a diamond activation pattern superficially in the right ventricle. Electromechanically, there was very little change in pressure volume loops although the stress distribution was altered. In conclusion, we clarified that the right ventricular septum has a downwardly directed superficial layer in larger mammalian species, which can have modest effects on stress distribution.

NEW & NOTEWORTHY A dual-layer organization of the fiber orientation in the intraventricular septum was identified in ex vivo hearts of large mammals. The RV septum has a downwardly directed superficial layer that is continuous with the moderator band. Electrically, it produced a diamond activation pattern. Electromechanically, little change in pressure volume loops were noticed but stress distribution was altered. Fiber distribution derived from diffusion tensor imaging should be considered for an accurate strain and stress analysis.

Personalization of biomechanical simulations of the left ventricle by in-vivo cardiac DTI data: Impact of fiber interpolation methods

Simulations of cardiac electrophysiology and mechanics have been reported to be sensitive to the microstructural anisotropy of the myocardium. Consequently, a personalized representation of cardiac microstructure is a crucial component of accurate, personalized cardiac biomechanical models. In-vivo cardiac Diffusion Tensor Imaging (cDTI) is a non-invasive magnetic resonance imaging technique capable of probing the heart's microstructure. Being a rather novel technique, issues such as low resolution, signal-to noise ratio, and spatial coverage are currently limiting factors. We outline four interpolation techniques with varying degrees of data fidelity, different amounts of smoothing strength, and varying representation error to bridge the gap between the sparse in-vivo data and the model, requiring a 3D representation of microstructure across the myocardium. We provide a workflow to incorporate in-vivo myofiber orientation into a left ventricular model and demonstrate that personalized modelling based on fiber orientations from in-vivo cDTI data is feasible. The interpolation error is correlated with a trend in personalized parameters and simulated physiological parameters, strains, and ventricular twist. This trend in simulation results is consistent across material parameter settings and therefore corresponds to a bias introduced by the interpolation method. This study suggests that using a tensor interpolation approach to personalize microstructure with in-vivo cDTI data, reduces the fiber uncertainty and thereby the bias in the simulation results.

Manifold-based respiratory phase estimation enables motion and distortion correction of free-breathing cardiac diffusion tensor MRI

PURPOSE

For in vivo cardiac DTI, breathing motion and B₀ field inhomogeneities produce misalignment and geometric distortion in diffusion-weighted (DW) images acquired with conventional single-shot EPI. We propose using a dimensionality reduction method to retrospectively estimate the respiratory phase of DW images and facilitate both distortion correction (DisCo) and motion compensation.

METHODS

Free-breathing electrocardiogram-triggered whole left-ventricular cardiac DTI using a second-order motion-compensated spin echo EPI sequence and alternating directionality of phase encoding blips was performed on 11 healthy volunteers. The respiratory phase of each DW image was estimated after projecting the DW images into a 2D space with Laplacian eigenmaps. DisCo and motion compensation were applied to the respiratory sorted DW images. The results were compared against conventional breath-held T₂ half-Fourier single shot turbo spin echo. Cardiac DTI parameters including fractional anisotropy, mean diffusivity, and helix angle transmurality were compared with and without DisCo.

RESULTS

The left-ventricular geometries after DisCo and motion compensation resulted in significantly improved alignment of DW images with T₂ reference. DisCo reduced the distance between the left-ventricular contours by 13.2% ± 19.2%, P < .05 (2.0 ± 0.4 for DisCo and 2.4 ± 0.5 mm for uncorrected). DisCo DTI parameter maps yielded no significant differences (mean diffusivity: 1.55 ± 0.13 × 10^-3 mm² /s and 1.53 ± 0.13 × 10^-3 mm² /s, P = .09; fractional anisotropy: 0.375 ± 0.041 and 0.379 ± 0.045, P = .11; helix angle transmurality: 1.00% ± 0.10°/% and 0.99% ± 0.12°/%, P = .44), although the orientation of individual tensors differed.

CONCLUSION

Retrospective respiratory phase estimation with LE-based DisCo and motion compensation in free-breathing cardiac DTI resulting in significantly reduced geometric distortion and improved alignment within and across slices.

Biomechanical Modeling to Inform Pulmonary Valve Replacement in Tetralogy of Fallot Patients After Complete Repair

BACKGROUND

A biomechanical model of the heart can be used to incorporate multiple data sources (electrocardiography, imaging, invasive hemodynamics). The purpose of this study was to use this approach in a cohort of patients with tetralogy of Fallot after complete repair (rTOF) to assess comparative influences of residual right ventricular outflow tract obstruction (RVOTO) and pulmonary regurgitation on ventricular health.

METHODS

Twenty patients with rTOF who underwent percutaneous pulmonary valve replacement (PVR) and cardiovascular magnetic resonance imaging were included in this retrospective study. Biomechanical models specific to individual patient and physiology (before and after PVR) were created and used to estimate the RV myocardial contractility. The ability of models to capture post-PVR changes of right ventricular (RV) end-diastolic volume (EDV) and effective flow in the pulmonary artery (Q_eff) was also compared with expected values.

RESULTS

RV contractility before PVR (mean 66 ± 16 kPa, mean ± standard deviation) was increased in patients with rTOF compared with normal RV (38-48 kPa) (P < 0.05). The contractility decreased significantly in all patients after PVR (P < 0.05). Patients with predominantly RVOTO demonstrated greater reduction in contractility (median decrease 35%) after PVR than those with predominant pulmonary regurgitation (median decrease 11%). The model simulated post-PVR decreased EDV for the majority and suggested an increase of Q_eff-both in line with published data.

CONCLUSIONS

This study used a biomechanical model to synthesize multiple clinical inputs and give an insight into RV health. Individualized modeling allows us to predict the RV response to PVR. Initial data suggest that residual RVOTO imposes greater ventricular work than isolated pulmonary regurgitation.

Cardiovascular magnetic resonance imaging of functional and microstructural changes of the heart in a longitudinal pig model of acute to chronic myocardial infarction

BACKGROUND

We examined the dynamic response of the myocardium to infarction in a longitudinal porcine study using relaxometry, functional as well as diffusion cardiovascular magnetic resonance (CMR). We sought to compare non contrast CMR methods like relaxometry and in-vivo diffusion to contrast enhanced imaging and investigate the link of microstructural and functional changes in the acute and chronically infarcted heart.

METHODS

CMR was performed on five myocardial infarction pigs and four healthy controls. In the infarction group, measurements were obtained 2 weeks before 90 min occlusion of the left circumflex artery, 6 days after ischemia and at 5 as well as 9 weeks as chronic follow-up. The timing of measurements was replicated in the control cohort. Imaging consisted of functional cine imaging, 3D tagging, T2 mapping, native as well as gadolinium enhanced T1 mapping, cardiac diffusion tensor imaging, and late gadolinium enhancement imaging.

RESULTS

Native T1, extracellular volume (ECV) and mean diffusivity (MD) were significantly elevated in the infarcted region while fractional anisotropy (FA) was significantly reduced. During the transition from acute to chronic stages, native T1 presented minor changes (< 3%). ECV as well as MD increased from acute to the chronic stages compared to baseline: ECV: 125 ± 24% (day 6) 157 ± 24% (week 5) 146 ± 60% (week 9), MD: 17 ± 7% (day 6) 33 ± 14% (week 5) 29 ± 15% (week 9) and FA was further reduced: - 31 ± 10% (day 6) - 38 ± 8% (week 5) - 36 ± 14% (week 9). T2 as marker for myocardial edema was significantly increased in the ischemic area only during the acute stage (83 ± 3 ms infarction vs. 58 ± 2 ms control p < 0.001 and 61 ± 2 ms in the remote area p < 0.001). The analysis of functional imaging revealed reduced left ventricular ejection fraction, global longitudinal strain and torsion in the infarct group. At the same time the transmural helix angle (HA) gradient was steeper in the chronic follow-up and a correlation between longitudinal strain and transmural HA gradient was detected (r = 0.59 with p < 0.05). Comparing non-gadolinium enhanced data T2 mapping showed the largest relative change between infarct and remote during the acute stage (+ 33 ± 4% day 6, with p = 0.013 T2 vs. MD, p = 0.009 T2 vs. FA and p = 0.01 T2 vs. T1) while FA exhibited the largest relative change between infarct and remote during the chronic follow-up (+ 31 ± 2% week 5, with p = N.S. FA vs. MD, p = 0.03 FA vs. T2 and p = 0.003 FA vs. T1). Overall, diffusion parameters provided a higher contrast (> 23% for MD and > 27% for FA) during follow-up compared to relaxometry (T1 17-18%/T2 10-20%).

CONCLUSION

During chronic follow-up after myocardial infarction, cardiac diffusion tensor imaging provides a higher sensitivity for mapping microstructural alterations when compared to non-contrast enhanced relaxometry with the added benefit of providing directional tensor information to assess remodelling of myocyte aggregate orientations, which cannot be otherwise assessed.

Clinical Translation of Three-Dimensional Scar, Diffusion Tensor Imaging, Four-Dimensional Flow, and Quantitative Perfusion in Cardiac MRI: A Comprehensive Review

Cardiovascular magnetic resonance (CMR) imaging is a versatile tool that has established itself as the reference method for functional assessment and tissue characterisation. CMR helps to diagnose, monitor disease course and sub-phenotype disease states. Several emerging CMR methods have the potential to offer a personalised medicine approach to treatment. CMR tissue characterisation is used to assess myocardial oedema, inflammation or thrombus in various disease conditions. CMR derived scar maps have the potential to inform ablation therapy—both in atrial and ventricular arrhythmias. Quantitative CMR is pushing boundaries with motion corrections in tissue characterisation and first-pass perfusion. Advanced tissue characterisation by imaging the myocardial fibre orientation using diffusion tensor imaging (DTI), has also demonstrated novel insights in patients with cardiomyopathies. Enhanced flow assessment using four-dimensional flow (4D flow) CMR, where time is the fourth dimension, allows quantification of transvalvular flow to a high degree of accuracy for all four-valves within the same cardiac cycle. This review discusses these emerging methods and others in detail and gives the reader a foresight of how CMR will evolve into a powerful clinical tool in offering a precision medicine approach to treatment, diagnosis, and detection of disease.

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.

Cardiovascular magnetic resonance imaging: emerging techniques and applications

This review gives examples of emerging cardiovascular magnetic resonance (CMR) techniques and applications that have the potential to transition from research to clinical application in the near future. Four-dimensional flow CMR (4D-flow CMR) allows time-resolved three-directional, three-dimensional (3D) velocity-encoded phase-contrast imaging for 3D visualisation and quantification of valvular or intracavity flow. Acquisition times of under 10 min are achievable for a whole heart multidirectional data set and commercial software packages are now available for data analysis, making 4D-flow CMR feasible for inclusion in clinical imaging protocols. Diffusion tensor imaging (DTI) is based on the measurement of molecular water diffusion and uses contrasting behaviour in the presence and absence of boundaries to infer tissue structure. Cardiac DTI is capable of non-invasively phenotyping the 3D micro-architecture within a few minutes, facilitating transition of the method to clinical protocols. Hybrid positron emission tomography-magnetic resonance (PET-MR) provides quantitative PET measures of biological and pathological processes of the heart combined with anatomical, morphological and functional CMR imaging. Cardiac PET-MR offers opportunities in ischaemic, inflammatory and infiltrative heart disease.

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.

Diffusion tensor cardiovascular magnetic resonance in hypertrophic cardiomyopathy: a comparison of motion-compensated spin echo and stimulated echo techniques

Objectives

Diffusion tensor cardiovascular magnetic resonance (DT-CMR) interrogates myocardial microstructure. Two frequently used in vivo DT-CMR techniques are motion-compensated spin echo (M2-SE) and stimulated echo acquisition mode (STEAM). Whilst M2-SE is strain-insensitive and signal to noise ratio efficient, STEAM has a longer diffusion time and motion compensation is unnecessary. Here we compare STEAM and M2-SE DT-CMR in patients.

Materials and methods

Biphasic DT-CMR using STEAM and M2-SE, late gadolinium imaging and pre/post gadolinium T1-mapping were performed in a mid-ventricular short-axis slice, in ten hypertrophic cardiomyopathy (HCM) patients at 3 T.

Results

Adequate quality data were obtained from all STEAM, but only 7/10 (systole) and 4/10 (diastole) M2-SE acquisitions. Compared with STEAM, M2-SE yielded higher systolic mean diffusivity (MD) (p = 0.02) and lower fractional anisotropy (FA) (p = 0.02, systole). Compared with segments with neither hypertrophy nor late gadolinium, segments with both had lower systolic FA using M2-SE (p = 0.02) and trend toward higher MD (p = 0.1). The negative correlation between FA and extracellular volume fraction was stronger with STEAM than M2-SE (r² = 0.29, p < 0.001 STEAM vs. r² = 0.10, p = 0.003 M2-SE).

Discussion

In HCM, only STEAM reliably assesses biphasic myocardial microstructure. Higher MD and lower FA from M2-SE reflect the shorter diffusion times. Further work will relate DT-CMR parameters and microstructural changes in disease.

Electronic supplementary material

The online version of this article (10.1007/s10334-019-00799-3) contains supplementary material, which is available to authorized users.

Characterizing cardiac involvement in amyloidosis using cardiovascular magnetic resonance diffusion tensor imaging

BACKGROUND

In-vivo cardiovascular magnetic resonance (CMR) diffusion tensor imaging (DTI) allows imaging of alterations of cardiac fiber architecture in diseased hearts. Cardiac amyloidosis (CA) causes myocardial infiltration of misfolded proteins with unknown consequences for myocardial microstructure. This study applied CMR DTI in CA to assess microstructural alterations and their consequences for myocardial function compared to healthy controls.

METHODS

Ten patients with CA (8 AL, 2 ATTR) and ten healthy controls were studied using a diffusion-weighed second-order motion-compensated spin-echo sequence at 1.5 T. Additionally, left ventricular morphology, ejection fraction, strain and native T1 values were obtained in all subjects. In CA patients, T1 mapping was repeated after the administration of gadolinium for extracellular volume fraction (ECV) calculation. CMR DTI analysis was performed to yield the scalar diffusion metrics mean diffusivity (MD) and fractional anisotropy (FA) as well as the characteristics of myofiber orientation including helix, transverse and E2A sheet angle (HA, TA, E2A).

RESULTS

MD and FA were found to be significantly different between CA patients and healthy controls (MD 1.77 ± 0.17 10- 3 vs 1.41 ± 0.07 10- 3 mm2/s, p <  0.001; FA 0.25 ± 0.04 vs 0.35 ± 0.03, p <  0.001). MD demonstrated an excellent correlation with native T1 (r = 0.908, p <  0.001) while FA showed a significant correlation with ECV in the CA population (r = - 0.851, p <  0.002). HA exhibited a more circumferential orientation of myofibers in CA patients, in conjunction with a higher TA standard deviation and a higher absolute E2A sheet angle. The transmural HA slope was found to be strongly correlated with the global longitudinal strain (r = 0.921, p < 0.001).

CONCLUSION

CMR DTI reveals significant alterations of scalar diffusion metrics in CA patients versus healthy controls. Elevated MD and lower FA values indicate myocardial disarray with higher diffusion in CA that correlates well with native T1 and ECV measures. In CA patients, CMR DTI showed pronounced circumferential orientation of the myofibers, which may provide the rationale for the reduction of global longitudinal strain that occurs in amyloidosis patients. Accordingly, CMR DTI captures specific features of amyloid infiltration, which provides a deeper understanding of the microstructural consequences of CA.

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.

Myocardial Fiber Mapping of Rat Hearts, Using Apparent Backscatter, with Histologic Validation

Myocardial fiber architecture is a physiologically important regulator of ejection fraction, strain and pressure development. Apparent ultrasonic backscatter has been shown to be a useful method for recreating the myocardial fiber architecture in human-sized sheep hearts because of the dependence of its amplitude on the relative orientation of a myofiber to the angle of ultrasonic insonification. Thus, the anisotropy of the backscatter signal is linked to and provides information about the fiber orientation. In this study, we sought to determine whether apparent backscatter could be used to measure myofiber orientation in rodent hearts. Fixed adult-rat hearts were imaged intact, and both a transmural cylindrical core and transmural wedge of the left ventricular free wall were imaged. Cylindrical core samples confirmed that backscatter anisotropy could be measured in rat hearts. Ultrasound and histologic analysis of transmural myocardial wedge samples confirmed that the apparent backscatter could be reproducibly mapped to fiber orientation (angle of the fiber relative to the direction of insonification). These data provided a quantitative relationship between the apparent backscatter and fiber angle that was applied to whole-heart images. Myocardial fiber architecture was successfully measured in rat hearts. Quantifying myocardial fiber architecture, using apparent backscatter, provides a number of advantages, including its scalable use from rodents to man, its rapid low-cost acquisition and minimal contraindications. The method outlined in this study provides a method for investigators to begin detailed assessments of how the myocardial fiber architecture changes in preclinical disease models, which can be immediately translated into the clinic.

Characterization of Myocardial Microstructure and Function in an Experimental Model of Isolated Subendocardial Damage

Subendocardial damage is among the first cardiac manifestations of hypertension and is already present in asymptomatic disease states. Accordingly, markers of subendocardial impairment may facilitate early detection of cardiac damages and risk stratification under these conditions. This study aimed to investigate the impact of subendocardial damage on myocardial microstructure and function to elucidate early pathophysiologic processes and to identify corresponding diagnostic measures. Mice (n=38) were injected with isoproterenol to induce isolated subendocardial scarring or saline as corresponding control. Cardiac function and myocardial deformation were determined by high-frequency echocardiography. The cardiac stress response was assessed in a graded exercise test and during dobutamine stress echocardiography. Myocardial microstructure was studied ex vivo by 7 T diffusion tensor magnetic resonance imaging at a spatial resolution of 100×100×100 µm ³ . Results were correlated with histology and biomarker expression. Subendocardial fibrosis was accompanied by diastolic dysfunction, impaired longitudinal deformation (global peak longitudinal strain [LS]: -12.5±0.5% versus -15.6±0.5%; P<0.001) and elevated biomarker expression (ANP [atrial natriuretic peptide], Galectin-3, and ST2). Systolic function and cardiac stress response remained preserved. Diffusion tensor magnetic resonance imaging revealed a left-shift in helix angle towards lower values in isoproterenol-treated animals, which was mainly determined by subepicardial myofibers (mean helix angle: 2.2±0.8° versus 5.9±1.0°; P<0.01). Longitudinal strain and subepicardial helix angle were highly predictive for subendocardial fibrosis (sensitivity, 82%-92% and specificity, 89%-90%). The results indicate that circumscribed subendocardial damage alone can cause several hallmarks observed in cardiovascular high-risk patients. Microstructural remodeling under these conditions involves also remote regions, and corresponding changes in longitudinal strain and helix angle might serve as diagnostic markers.

Diffusion tensor cardiovascular magnetic resonance

Cardiac structure and function are complex and inter-related. Current in vivo techniques assess the heart on a macroscopic scale, but a novel technique called diffusion tensor cardiovascular magnetic resonance (DT-CMR) can now assess the cardiac microstructure non-invasively. It provides information on the helical arrangement of cardiomyocytes that drives torsion and offers dynamic assessment of the sheetlets (aggregated cardiomyocytes) that rotate through the cardiac cycle to facilitate wall thickening. Through diffusion biomarkers, the expansion and organisation of the underlying myocardium can be described. DT-CMR has already identified novel microstructural abnormalities in cardiomyopathy, and ischaemic and congenital heart disease. This new knowledge supports the potential of DT-CMR to improve diagnostics and prognostication in various cardiac diseases.

Postmortem cardiac magnetic resonance in sudden cardiac death

Postmortem imaging is increasingly used in forensic practice as good complementary tool to conventional autopsy investigations. Over the last decade, postmortem cardiac magnetic resonance (PMCMR) imaging was introduced in forensic investigations of natural deaths related to cardiovascular diseases, which represent the most common causes of death in developed countries. Postmortem CMR application has yielded interesting results in ischemic myocardium injury investigations and in visualizing other pathological findings in the heart. This review presents the actual state of postmortem imaging for cardiovascular pathologies in cases of sudden cardiac death (SCD), taking into consideration both the advantages and limitations of PMCMR application.

Novel insights into in-vivo diffusion tensor cardiovascular magnetic resonance using computational modeling and a histology-based virtual microstructure

Purpose

To develop histology‐informed simulations of diffusion tensor cardiovascular magnetic resonance (DT‐CMR) for typical in‐vivo pulse sequences and determine their sensitivity to changes in extra‐cellular space (ECS) and other microstructural parameters.

Methods

We synthesised the DT‐CMR signal from Monte Carlo random walk simulations. The virtual tissue was based on porcine histology. The cells were thickened and then shrunk to modify ECS. We also created idealised geometries using cuboids in regular arrangement, matching the extra‐cellular volume fraction (ECV) of 16–40%. The simulated voxel size was 2.8 × 2.8 × 8.0 mm³ for pulse sequences covering short and long diffusion times: Stejskal–Tanner pulsed‐gradient spin echo, second‐order motion‐compensated spin echo, and stimulated echo acquisition mode (STEAM), with clinically available gradient strengths.

Results

The primary diffusion tensor eigenvalue increases linearly with ECV at a similar rate for all simulated geometries. Mean diffusivity (MD) varies linearly, too, but is higher for the substrates with more uniformly distributed ECS. Fractional anisotropy (FA) for the histology‐based geometry is higher than the idealised geometry with low sensitivity to ECV, except for the long mixing time of the STEAM sequence. Varying the intra‐cellular diffusivity (D_IC) results in large changes of MD and FA. Varying extra‐cellular diffusivity or using stronger gradients has minor effects on FA. Uncertainties of the primary eigenvector orientation are reduced using STEAM.

Conclusions

We found that the distribution of ECS has a measurable impact on DT‐CMR parameters. The observed sensitivity of MD and FA to ECV and D_IC has potentially interesting applications for interpreting in‐vivo DT‐CMR parameters.

Left Ventricular Diastolic Myocardial Stiffness and End-Diastolic Myofibre Stress in Human Heart Failure Using Personalised Biomechanical Analysis

Understanding the aetiology of heart failure with preserved (HFpEF) and reduced (HFrEF) ejection fraction requires knowledge of biomechanical factors such as diastolic myocardial stiffness and stress. Cine CMR images and intra-ventricular pressure recordings were acquired in 8 HFrEF, 11 HFpEF and 5 control subjects. Diastolic myocardial stiffness was estimated using biomechanical models and found to be greater in HFrEF (6.4 ± 1.2 kPa) than HFpEF (2.7 ± 0.6 kPa, p < 0.05) and also greater than control (1.2 ± 0.4 kPa, p < 0.005). End-diastolic mid-ventricular myofibre stress derived from the personalised biomechanics model was higher in HFrEF (2.9 ± 0.3 kPa) than control (0.9 ± 0.3 kPa, p < 0.01). Chamber stiffness, measured from the slope of the diastolic pressure-volume relationship, is determined by the intrinsic tissue properties as well as the size and shape of the heart, and was unable to distinguish between any of the three groups (p > 0.05). Personalised biomechanical analysis may provide more specific information about myocardial mechanical behaviour than global chamber indices, which are confounded by variations in ventricular geometry.

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.

An in-vivo comparison of stimulated-echo and motion compensated spin-echo sequences for 3 T diffusion tensor cardiovascular magnetic resonance at multiple cardiac phases

BACKGROUND

Stimulated-echo (STEAM) and, more recently, motion-compensated spin-echo (M2-SE) techniques have been used for in-vivo diffusion tensor cardiovascular magnetic resonance (DT-CMR) assessment of cardiac microstructure. The two techniques differ in the length scales of diffusion interrogated, their signal-to-noise ratio efficiency and sensitivity to both motion and strain. Previous comparisons of the techniques have used high performance gradients at 1.5 T in a single cardiac phase. However, recent work using STEAM has demonstrated novel findings of microscopic dysfunction in cardiomyopathy patients, when DT-CMR was performed at multiple cardiac phases. We compare STEAM and M2-SE using a clinical 3 T scanner in three potentially clinically interesting cardiac phases.

METHODS

Breath hold mid-ventricular short-axis DT-CMR was performed in 15 subjects using M2-SE and STEAM at end-systole, systolic sweet-spot and diastasis. Success was defined by ≥50% of the myocardium demonstrating normal helix angles. From successful acquisitions DT-CMR results relating to tensor orientation, size and shape were compared between sequences and cardiac phases using non-parametric statistics. Strain information was obtained using cine spiral displacement encoding with stimulated echoes for comparison with DT-CMR results.

RESULTS

Acquisitions were successful in 98% of STEAM and 76% of M2-SE cases and visual helix angle (HA) map scores were higher for STEAM at the sweet-spot and diastasis. There were significant differences between sequences (p < 0.05) in mean diffusivity (MD), fractional anisotropy (FA), tensor mode, transmural HA gradient and absolute second eigenvector angle (E2A). Differences in E2A between systole and diastole correlated with peak radial strain for both sequences (p ≤ 0.01).

CONCLUSION

M2-SE and STEAM can be performed equally well at peak systole at 3 T using standard gradients, but at the sweet-spot and diastole STEAM is more reliable and image quality scores are higher. Differences in DT-CMR results are potentially due to differences in motion sensitivity and the longer diffusion time of STEAM, although the latter appears to be the dominant factor. The benefits of both sequences should be considered when planning future studies and sequence and cardiac phase specific normal ranges should be used for comparison.

Diffusion tensor cardiovascular magnetic resonance with a spiral trajectory: An in vivo comparison of echo planar and spiral stimulated echo sequences

PURPOSE

Diffusion tensor cardiovascular MR (DT-CMR) using stimulated echo acquisition mode (STEAM) with echo-planar-imaging (EPI) readouts is a low signal-to-noise-ratio (SNR) technique and therefore typically has a low spatial resolution. Spiral trajectories are more efficient than EPI, and could increase the SNR. The purpose of this study was to compare the performance of a novel STEAM spiral DT-CMR sequence with an equivalent established EPI technique.

METHODS

A STEAM DT-CMR sequence was implemented with a spiral readout and a reduced field of view. An in vivo comparison of DT-CMR parameters and data quality between EPI and spiral was performed in 11 healthy volunteers imaged in peak systole and diastasis at 3 T. The SNR was compared in a phantom and in vivo.

RESULTS

There was a greater than 49% increase in the SNR in vivo and in the phantom measurements (in vivo septum, systole: SNR_EPI  = 8.0 ± 2.2, SNR_spiral  = 12.0 ± 2.7; diastasis: SNR_EPI  = 8.1 ± 1.6, SNR_spiral  = 12.0 ± 3.7). There were no significant differences in helix angle gradient (HAG) (systole: HAG_EPI  = -0.79 ± 0.07 °/%; HAG_spiral  = -0.74 ± 0.16 °/%; P = 0.11; diastasis: HAG_EPI  = -0.63 ± 0.05 °/%; HAG_spiral  = -0.56 ± 0.14 °/%; P = 0.20), mean diffusivity (MD) in systole (MD_EPI  = 0.99 ± 0.06 × 10^-3 mm² /s, MD_spiral  = 1.00 ± 0.09 × 10^-3 mm² /s, P = 0.23) and secondary eigenvector angulation (E2A) (systole: E2A_EPI  = 61 ± 10 °; E2A_spiral  = 63 ± 10 °; P = 0.77; diastasis: E2A_EPI  = 18 ± 11 °; E2A_spiral  = 15 ± 8 °; P = 0.20) between the sequences. There was a small difference (≈ 20%) in fractional anisotropy (FA) (systole: FA_EPI  = 0.49 ± 0.03, FA_spiral  = 0.41 ± 0.04; P < 0.01; diastasis: FA_EPI  = 0.66 ± 0.05, FA_spiral  = 0.55 ± 0.03; P < 0.01) and mean diffusivity in diastasis (10%; MD_EPI  = 1.00 ± 0.12 × 10^-3 mm² /s, MD_spiral  = 1.10 ± 0.09 × 10^-3 mm² /s, P = 0.02).

CONCLUSION

This is the first study to demonstrate DT-CMR STEAM using a spiral trajectory. The SNR was increased by using a spiral rather than the more established EPI readout, and the DT-CMR parameters were largely similar between the two sequences. Magn Reson Med 80:648-654, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

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.