Optimal diffusion weighting for in vivo cardiac diffusion tensor imaging

Phase stabilization with motion compensated diffusion weighted imaging

PURPOSE

Diffusion encoding gradient waveforms can impart intra-voxel and inter-voxel dephasing owing to bulk motion, limiting achievable signal-to-noise and complicating multishot acquisitions. In this study, we characterize improvements in phase consistency via gradient moment nulling of diffusion encoding waveforms.

METHODS

Healthy volunteers received neuro (N = 10 $$ N=10 $$ ) and cardiac (N = 10 $$ N=10 $$ ) MRI. Three gradient moment nulling levels were evaluated: compensation for position (M 0 $$ {M}_0 $$ ), position + velocity (M 1 $$ {M}_1 $$ ), and position + velocity + acceleration (M 1 + M 2 $$ {M}_1+{M}_2 $$ ). Three experiments were completed: (Exp-1) Fixed Trigger Delay Neuro DWI; (Exp-2) Mixed Trigger Delay Neuro DWI; and (Exp-3) Fixed Trigger Delay Cardiac DWI. Significant differences (p < 0 . 05 $$ p<0.05 $$ ) of the temporal phase SD between repeated acquisitions and the spatial phase gradient across a given image were assessed.

RESULTS

M 0 $$ {M}_0 $$ moment nulling was a reference for all measures. In Exp-1, temporal phase SD forG z $$ {G}_z $$ diffusion encoding was significantly reduced withM 1 $$ {M}_1 $$ (35% of t-tests) andM 1 + M 2 $$ {M}_1+{M}_2 $$ (68% of t-tests). The spatial phase gradient was reduced in 23% of t-tests forM 1 $$ {M}_1 $$ and 2% of cases forM 1 + M 2 $$ {M}_1+{M}_2 $$ . In Exp-2, temporal phase SD significantly decreased withM 1 + M 2 $$ {M}_1+{M}_2 $$ gradient moment nulling only forG z $$ {G}_z $$ (83% of t-tests), but spatial phase gradient significantly decreased with onlyM 1 $$ {M}_1 $$ (50% of t-tests). In Exp-3,M 1 + M 2 $$ {M}_1+{M}_2 $$ gradient moment nulling significantly reduced temporal phase SD and spatial phase gradients (100% of t-tests), resulting in less signal attenuation and more accurate ADCs.

CONCLUSION

We characterized gradient moment nulling phase consistency for DWI. Using M1 for neuroimaging and M1 + M2 for cardiac imaging minimized temporal phase SDs and spatial phase gradients.

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.

Deep learning method with integrated invertible wavelet scattering for improving the quality ofin vivocardiac DTI

Objective.Respiratory motion, cardiac motion and inherently low signal-to-noise ratio (SNR) are major limitations ofin vivocardiac diffusion tensor imaging (DTI). We propose a novel enhancement method that uses unsupervised learning based invertible wavelet scattering (IWS) to improve the quality ofin vivocardiac DTI.Approach.Our method starts by extracting nearly transformation-invariant features from multiple cardiac diffusion-weighted (DW) image acquisitions using multi-scale wavelet scattering (WS). Then, the relationship between the WS coefficients and DW images is learned through a multi-scale encoder and a decoder network. Using the trained encoder, the deep features of WS coefficients of multiple DW image acquisitions are further extracted and then fused using an average rule. Finally, using the fused WS features and trained decoder, the enhanced DW images are derived.Main result.We evaluate the performance of the proposed method by comparing it with several methods on threein vivocardiac DTI datasets in terms of SNR, contrast to noise ratio (CNR), fractional anisotropy (FA), mean diffusivity (MD) and helix angle (HA). Comparing against the best comparison method, SNR/CNR of diastolic, gastric peristalsis influenced, and end-systolic DW images were improved by 1%/16%, 5%/6%, and 56%/30%, respectively. The approach also yielded consistent FA and MD values and more coherent helical fiber structures than the comparison methods used in this work.Significance.The ablation results verify that using the transformation-invariant and noise-robust wavelet scattering features enables us to effectively explore the useful information from the limited data, providing a potential mean to alleviate the dependence of the fusion results on the number of repeated acquisitions, which is beneficial for dealing with the issues of noise and residual motion simultaneously and therefore improving the quality ofinvivocardiac DTI. Code can be found inhttps://github.com/strawberry1996/WS-MCNN.

The effects of field strength on stimulated echo and motion-compensated spin-echo diffusion tensor cardiovascular magnetic resonance sequences

BACKGROUND

In-vivo diffusion tensor cardiovascular magnetic resonance (DT-CMR) is an emerging technique for microstructural tissue characterization in the myocardium. Most studies are performed at 3T, where higher signal-to-noise ratio (SNR) should benefit this signal-starved method. However, a few studies have suggested that DT-CMR is possible at 1.5T, where echo planar imaging artifacts may be less severe and 1.5T hardware is more widely available.

METHODS

We recruited 20 healthy volunteers and performed mid-ventricular short-axis DT-CMR at 1.5T and 3T. Acquisitions were performed at peak systole and end-diastole using both stimulated echo acquisition mode (STEAM) and motion-compensated spin-echo (MCSE) sequences at matched spatial resolutions. DT-CMR parameters were averaged over the left ventricle and compared between 1.5T and 3T sequences using both datasets with and without the blow reference data included.

RESULTS

Eleven (1.5T) and 12 (3T) diastolic MCSE acquisitions were rejected as the helix angle (HA) demonstrated <50% normal appearance circumferentially or the acquisition was abandoned due to poor image quality; a maximum of one acquisition was rejected for other datasets. Subjective HA map quality was significantly better at 3T than 1.5T for STEAM (p < 0.05), but not for MCSE and other DT-CMR quality measures were consistent with improvements in STEAM at 3T over 1.5T. When blow data were excluded, no significant differences in mean diffusivity were observed between field strengths, but fractional anisotropy was significantly higher at 1.5T than 3T for STEAM systole (p < 0.05). Absolute second eigenvector orientation (E2A, sheetlet angle) was significantly higher at 1.5T than 3T for MCSE systole and STEAM diastole, but significantly lower for STEAM systole (all p < 0.05). Transmural HA distribution was less steep at 1.5T than 3T for STEAM diastole data (p < 0.05). SNR was higher at 3T than 1.5T for all acquisitions (p < 0.05).

CONCLUSION

While 3T provides benefits in terms of SNR, both STEAM and MCSE can be performed at 1.5T. However, MCSE is unreliable in diastole at both field strengths and STEAM benefits from the improved SNR at 3T over 1.5T. Future clinical research studies may be able to leverage the wider availability of 1.5T CMR hardware where MCSE acquisitions are desirable.

Realistic numerical simulations of diffusion tensor cardiovascular magnetic resonance: The effects of perfusion and membrane permeability

PURPOSE

To study the sensitivity of diffusion tensor cardiovascular magnetic resonance (DT-CMR) to microvascular perfusion and changes in cell permeability.

METHODS

Monte Carlo (MC) random walk simulations in the myocardium have been performed to simulate self-diffusion of water molecules in histology-based media with varying extracellular volume fraction (ECV) and permeable membranes. The effect of microvascular perfusion on simulations of the DT-CMR signal has been incorporated by adding the contribution of particles traveling through an anisotropic capillary network to the diffusion signal. The simulations have been performed considering three pulse sequences with clinical gradient strengths: monopolar stimulated echo acquisition mode (STEAM), monopolar pulsed-gradient spin echo (PGSE), and second-order motion-compensated spin echo (MCSE).

RESULTS

Reducing ECV intensifies the diffusion restriction and incorporating membrane permeability reduces the anisotropy of the diffusion tensor. Widening the intercapillary velocity distribution results in increased measured diffusion along the cardiomyocytes long axis when the capillary networks are anisotropic. Perfusion amplifies the mean diffusivity for STEAM while the opposite is observed for short diffusion encoding time sequences (PGSE and MCSE).

CONCLUSION

The effect of perfusion on the measured diffusion tensor is reduced using an increased reference b-value. Our results pave the way for characterization of the response of DT-CMR to microstructural changes underlying cardiac pathology and highlight the higher sensitivity of STEAM to permeability and microvascular circulation due to its longer diffusion encoding time.

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.

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.

Comparison of interpolation methods of predominant cardiomyocyte orientation from in vivo and ex vivo cardiac diffusion tensor imaging data

Cardiac electrophysiology and cardiac mechanics both depend on the average cardiomyocyte long-axis orientation. In the realm of personalized medicine, knowledge of the patient-specific changes in cardiac microstructure plays a crucial role. Patient-specific computational modelling has emerged as a tool to better understand disease progression. In vivo cardiac diffusion tensor imaging (cDTI) is a vital tool to non-destructively measure the average cardiomyocyte long-axis orientation in the heart. However, cDTI suffers from long scan times, rendering volumetric, high-resolution acquisitions challenging. Consequently, interpolation techniques are needed to populate bio-mechanical models with patient-specific average cardiomyocyte long-axis orientations. In this work, we compare five interpolation techniques applied to in vivo and ex vivo porcine input data. We compare two tensor interpolation approaches, one rule-based approximation, and two data-driven, low-rank models. We demonstrate the advantage of tensor interpolation techniques, resulting in lower interpolation errors than do low-rank models and rule-based methods adapted to cDTI data. In an ex vivo comparison, we study the influence of three imaging parameters that can be traded off against acquisition time: in-plane resolution, signal to noise ratio, and number of acquired short-axis imaging slices.

Short-Term Repeatability of in Vivo Cardiac Intravoxel Incoherent Motion Tensor Imaging in Healthy Human Volunteers

BACKGROUND

Intravoxel incoherent motion (IVIM) tensor imaging is a promising technique for diagnosis and monitoring of cardiovascular diseases. Knowledge about measurement repeatability, however, remains limited.

PURPOSE

To evaluate short-term repeatability of IVIM tensor imaging in normal in vivo human hearts.

STUDY TYPE

Prospective.

POPULATION

Ten healthy subjects without history of heart diseases.

FIELD STRENGTH/SEQUENCE

Balanced steady-state free-precession cine sequence and single-shot spin-echo echo planar IVIM tensor imaging sequence (9 b-values, 0-400 seconds/mm² and six diffusion-encoding directions) at 3.0 T.

ASSESSMENT

Subjects were scanned twice with an interval of 15 minutes, leaving the scanner between studies. The signal-to-noise ratio (SNR) was evaluated in anterior, lateral, septal, and inferior segments of the left ventricle wall. Fractional anisotropy (FA), mean diffusivity (MD), mean fraction (MF), and helix angle (HA) in the four segments were independently measured by five radiologists.

STATISTICAL TESTS

IVIM tensor indexes were compared between observers using a one-way analysis of variance or between scans using a paired t-test (normal data) or a Wilcoxon rank-sum test (non-normal data). Interobserver agreement and test-retest repeatability were assessed using the intraclass correlation coefficient (ICC), within-subject coefficient of variation (WCV), and Bland-Altman limits of agreements.

RESULTS

SNR of inferior segment was significantly lower than the other three segments, and inferior segment was therefore excluded from repeatability analysis. Interobserver repeatability was excellent for all IVIM tensor indexes (ICC: 0.886-0.972; WCV: 0.62%-4.22%). Test-retest repeatability was excellent for MD of the self-diffusion tensor (D) and MF of the perfusion fraction tensor (f_p ) (ICC: 0.803-0.888; WCV: 1.42%-9.51%) and moderate for FA and MD of the pseudo-diffusion tensor (D^* ) (ICC: 0.487-0.532; WCV: 6.98%-10.89%). FA of D and f_p and HA of D presented good test-retest repeatability (ICC: 0.732-0.788; WCV: 3.28%-8.71%).

DATA CONCLUSION

The D and f_p indexes exhibited satisfactory repeatability, but further efforts were needed to improve repeatability of D^* indexes.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY: Stage 1.

Effect of b Value on Imaging Quality for Diffusion Tensor Imaging of the Spinal Cord at Ultrahigh Field Strength

Objective

To explore the optimal b value setting for diffusion tensor imaging of rats' spinal cord at ultrahigh field strength (7 T).

Methods

Spinal cord diffusion tensor imaging data were collected from 14 rats (5 healthy, 9 spinal cord injured) with a series of b values (200, 300, 400, 500, 600, 700, 800, 900, and 1000 s/mm²) under the condition that other scanning parameters were consistent. The image quality (including image signal-to-noise ratio and image distortion degree) and data quality (i.e., the stability and consistency of the DTI-derived parameters, referred to as data stability and data consistency) were quantitatively evaluated. The min-max normalization method was used to process the calculation results of the four indicators. Finally, the image and data quality under each b value were synthesized to determine the optimal b value.

Results

b = 200 s/mm² and b = 900 s/mm² ranked in the top two of the comprehensive evaluation, with the best image quality at b = 200 s/mm² and the best data quality at b = 900 s/mm².

Conclusion

Considering the shortcomings of the ability of low b values to reflect the microstructure, b = 900 s/mm² can be used as the optimal b value for 7 T spinal cord diffusion tensor scanning.

Cardiac multi-scale investigation of the right and left ventricle ex vivo: a review

The heart is a complex multi-scale system composed of components integrated at the subcellular, cellular, tissue and organ levels. The myocytes, the contractile elements of the heart, form a complex three-dimensional (3D) network which enables propagation of the electrical signal that triggers the contraction to efficiently pump blood towards the whole body. Cardiovascular diseases (CVDs), a major cause of mortality in developed countries, often lead to cardiovascular remodeling affecting cardiac structure and function at all scales, from myocytes and their surrounding collagen matrix to the 3D organization of the whole heart. As yet, there is no consensus as to how the myocytes are arranged and packed within their connective tissue matrix, nor how best to image them at multiple scales. Cardiovascular imaging is routinely used to investigate cardiac structure and function as well as for the evaluation of cardiac remodeling in CVDs. For a complete understanding of the relationship between structural remodeling and cardiac dysfunction in CVDs, multi-scale imaging approaches are necessary to achieve a detailed description of ventricular architecture along with cardiac function. In this context, ventricular architecture has been extensively studied using a wide variety of imaging techniques: ultrasound (US), optical coherence tomography (OCT), microscopy (confocal, episcopic, light sheet, polarized light), magnetic resonance imaging (MRI), micro-computed tomography (micro-CT) and, more recently, synchrotron X-ray phase contrast imaging (SR X-PCI). Each of these techniques have their own set of strengths and weaknesses, relating to sample size, preparation, resolution, 2D/3D capabilities, use of contrast agents and possibility of performing together with in vivo studies. Therefore, the combination of different imaging techniques to investigate the same sample, thus taking advantage of the strengths of each method, could help us to extract the maximum information about ventricular architecture and function. In this review, we provide an overview of available and emerging cardiovascular imaging techniques for assessing myocardial architecture ex vivo and discuss their utility in being able to quantify cardiac remodeling, in CVDs, from myocyte to whole organ.

Investigation of intravoxel incoherent motion tensor imaging for the characterization of the in vivo human heart

PURPOSE

To investigate intravoxel incoherent motion (IVIM) tensor imaging of the in vivo human heart and elucidate whether the estimation of IVIM tensors is affected by the complexity of pseudo-diffusion components in myocardium.

METHODS

The cardiac IVIM data of 10 healthy subjects were acquired using a diffusion weighted spin-echo echo-planar imaging sequence along 6 gradient directions with 10 b values (0~400 s/mm² ). The IVIM data of left ventricle myocardium were fitted to the IVIM tensor model. The complexity of myocardial pseudo-diffusion components was reduced through exclusion of low b values (0 and 5 s/mm² ) from the IVIM curve-fitting analysis. The fractional anisotropy, mean fraction/mean diffusivity, and Westin measurements of pseudo-diffusion tensors (f_p and D*) and self-diffusion tensor (D), as well as the angle between the main eigenvector of f_p (or D*) and that of D, were computed and compared before and after excluding low b values.

RESULTS

The fractional anisotropy values of f_p and D* without low b value participation were significantly higher (P < .001) than those with low b value participation, but an opposite trend was found for the mean fraction/diffusivity values. Besides, after removing low b values, the angle between the main eigenvector of f_p (or D*) and that of D became small, and both f_p and D* tensors presented significant decrease of spherical components and significant increase of linear components.

CONCLUSION

The presence of multiple pseudo-diffusion components in myocardium indeed influences the estimation of IVIM tensors. The IVIM tensor model needs to be further improved to account for the complexity of myocardial microcirculatory network and blood flow.

Diffusion Tensor Cardiovascular Magnetic Resonance in Cardiac Amyloidosis

Background Cardiac amyloidosis (CA) is a disease of interstitial myocardial infiltration, usually by light chains or transthyretin. We used diffusion tensor cardiovascular magnetic resonance (DT-CMR) to noninvasively assess the effects of amyloid infiltration on the cardiac microstructure. Methods DT-CMR was performed at diastole and systole in 20 CA, 11 hypertrophic cardiomyopathy, and 10 control subjects with calculation of mean diffusivity, fractional anisotropy, and sheetlet orientation (secondary eigenvector angle). Results Mean diffusivity was elevated and fractional anisotropy reduced in CA compared with both controls and hypertrophic cardiomyopathy (P<0.001). In CA, mean diffusivity was correlated with extracellular volume (r=0.68, P=0.004), and fractional anisotropy was inversely correlated with circumferential strain (r=-0.65, P=0.02). In CA, diastolic secondary eigenvector angle was elevated, and secondary eigenvector angle mobility was reduced compared with controls (both P<0.001). Diastolic secondary eigenvector angle was correlated with amyloid burden measured by extracellular volume in transthyretin, but not light chain amyloidosis. Conclusions DT-CMR can characterize the microstructural effects of amyloid infiltration and is a contrast-free method to identify the location and extent of the expanded disorganized myocardium. The diffusion biomarkers mean diffusivity and fractional anisotropy effectively discriminate CA from hypertrophic cardiomyopathy. DT-CMR demonstrated that failure of sheetlet relaxation in diastole correlated with extracellular volume in transthyretin, but not light chain amyloidosis. This indicates that different mechanisms may be responsible for impaired contractility in CA, with an amyloid burden effect in transthyretin, but an idiosyncratic effect in light chain amyloidosis. Consequently, DT-CMR offers a contrast-free tool to identify novel pathophysiology, improve diagnostics, and monitor disease through noninvasive microstructural assessment.

Motion-compensated b-tensor encoding for in vivo cardiac diffusion-weighted imaging

Motion is a major confound in diffusion-weighted imaging (DWI) in the body, and it is a common cause of image artefacts. The effects are particularly severe in cardiac applications, due to the nonrigid cyclical deformation of the myocardium. Spin echo-based DWI commonly employs gradient moment-nulling techniques to desensitise the acquisition to velocity and acceleration, ie, nulling gradient moments up to the 2nd order (M2-nulled). However, current M2-nulled DWI scans are limited to encode diffusion along a single direction at a time. We propose a method for designing b-tensors of arbitrary shapes, including planar, spherical, prolate and oblate tensors, while nulling gradient moments up to the 2nd order and beyond. The design strategy comprises initialising the diffusion encoding gradients in two encoding blocks about the refocusing pulse, followed by appropriate scaling and rotation, which further enables nulling undesired effects of concomitant gradients. Proof-of-concept assessment of in vivo mean diffusivity (MD) was performed using linear and spherical tensor encoding (LTE and STE, respectively) in the hearts of five healthy volunteers. The results of the M2-nulled STE showed that (a) the sequence was robust to cardiac motion, and (b) MD was higher than that acquired using standard M2-nulled LTE, where diffusion-weighting was applied in three orthogonal directions, which may be attributed to the presence of restricted diffusion and microscopic diffusion anisotropy. Provided adequate signal-to-noise ratio, STE could significantly shorten estimation of MD compared with the conventional LTE approach. Importantly, our theoretical analysis and the proposed gradient waveform design may be useful in microstructure imaging beyond diffusion tensor imaging where the effects of motion must be suppressed.

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.

Reproducibility of Neurite Orientation Dispersion and Density Imaging (NODDI) in rats at 9.4 Tesla

PURPOSE

Neurite Orientation Dispersion and Density Imaging (NODDI) is a diffusion MRI (dMRI) technique used to characterize tissue microstructure by compartmental modelling of neural water fractions. Intra-neurite, extra-neurite, and cerebral spinal fluid volume fractions are measured. The purpose of this study was to determine the reproducibility of NODDI in the rat brain at 9.4 Tesla.

METHODS

Eight data sets were successfully acquired on adult male Sprague Dawley rats. Each rat was scanned twice on a 9.4T Agilent MRI with a 7 ± 1 day separation between scans. A multi-shell diffusion protocol was implemented consisting of 108 total directions varied over two shells (b-values of 1000 s/mm2 and 2000 s/mm2). Three techniques were used to analyze the NODDI scalar maps: mean region of interest (ROI) analysis, whole brain voxel-wise analysis, and targeted ROI analyses (voxel-wise within a given ROI). The coefficient of variation (CV) was used to assess the reproducibility of NODDI and provide insight into necessary sample sizes and minimum detectable effect size.

RESULTS

CV maps for orientation dispersion index (ODI) and neurite density index (NDI) showed high reproducibility both between and within subjects. Furthermore, it was found that small biological changes (<5%) may be detected with feasible sample sizes (n < 6-10). In contrast, isotropic volume fraction (IsoVF) was found to have low reproducibility, requiring very large sample sizes (n > 50) for biological changes to be detected.

CONCLUSIONS

The ODI and NDI measured by NODDI in the rat brain at 9.4T are highly reproducible and may be sensitive to subtle changes in tissue microstructure.

Motion-Induced Signal Loss in In Vivo Cardiac Diffusion-Weighted Imaging

LEVEL OF EVIDENCE

5 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2020;51:319-320.

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.

High resolution in-vivo DT-CMR using an interleaved variable density spiral STEAM sequence

PURPOSE

Diffusion tensor cardiovascular magnetic resonance (DT-CMR) has a limited spatial resolution. The purpose of this study was to demonstrate high-resolution DT-CMR using a segmented variable density spiral sequence with correction for motion, off-resonance, and T2*-related blurring.

METHODS

A single-shot stimulated echo acquisition mode (STEAM) echo-planar-imaging (EPI) DT-CMR sequence at 2.8 × 2.8 × 8 mm³ and 1.8 × 1.8 × 8 mm³ was compared to a single-shot spiral at 2.8 × 2.8 × 8 mm³ and an interleaved spiral sequence at 1.8 × 1.8 × 8 mm³ resolution in 10 healthy volunteers at peak systole and diastasis. Motion-induced phase was corrected using the densely sampled central k-space data of the spirals. STEAM field maps and T2* measures were obtained using a pair of stimulated echoes each with a double spiral readout, the first used to correct the motion-induced phase of the second.

RESULTS

The high-resolution spiral sequence produced similar DT-CMR results and quality measures to the standard-resolution sequence in both cardiac phases. Residual differences in fractional anisotropy and helix angle gradient between the resolutions could be attributed to spatial resolution and/or signal-to-noise ratio. Data quality increased after both motion-induced phase correction and off-resonance correction, and sharpness increased after T2* correction. The high-resolution EPI sequence failed to provide sufficient data quality for DT-CMR reconstruction.

CONCLUSION

In this study, an in vivo DT-CMR acquisition at 1.8 × 1.8 mm² in-plane resolution was demonstrated using a segmented spiral STEAM sequence. Motion-induced phase and off-resonance corrections are essential for high-resolution spiral DT-CMR. Segmented variable density spiral STEAM was found to be the optimal method for acquiring high-resolution DT-CMR data.

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.

Increased diffusion and decreased perfusion in epiphyses of femoral heads detected via intravoxel incoherent motion after closed reduction in children with developmental dysplasia of the hip

Background Patients with developmental dysplasia of the hip (DDH) may have decreased blood supply to the femoral heads. Finding a non-invasive method to evaluate whether the femoral heads in patients with DDH are ischemic is paramount for orthopedic surgeons. Purpose To identify whether parameters reflecting perfusion and diffusion in intravoxel incoherent motion (IVIM) sequences can be used to assess ischemia in femoral heads in patients with DDH after closed reduction. Material and Methods Twenty-eight patients with DDH who had undergone closed reduction were enrolled. IVIM data were acquired using a 3-T magnetic resonance scanner, regions of interest were placed on the epiphyses; ADC_slow, ADC_fast, f, and ADC_fast×f were measured. A Mann-Whitney U test was performed to compare ADC_slow, ADC_fast, f, and ADC_fast×f between the lesion and control sides. Receiver operating characteristic curves were generated with respective cut-off values. The lesion sides were classified based on the International Hip Dysplasia Institute (IHDI) classification. ADC_slow, ADC_fast, f, and ADC_fast×f were compared among the groups. Results ADC_slow was higher and ADC_fast, f, and ADC_fast×f were lower on the lesion sides ( P = 0.000-0.002). The optimal cut-off value for ADC_fast×f, ADC_fast, ADC_slow, and f were 0.030, 0.626, 0.000251, and 0.636, respectively. Higher IHDI classification scores on the lesion side were associated with lower ADC_fast, f, and ADC_fast×f, and higher ADC_slow values. Conclusion IVIM is a promising method to investigate the perfusion and diffusion of epiphyses of femoral heads.

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.

Evaluation of the impact of strain correction on the orientation of cardiac diffusion tensors with in vivo and ex vivo porcine hearts

PURPOSE

To evaluate the importance of strain-correcting stimulated echo acquisition mode echo-planar imaging cardiac diffusion tensor imaging.

METHODS

Healthy pigs (n = 11) were successfully scanned with a 3D cine displacement-encoded imaging with stimulated echoes and a monopolar-stimulated echo-planar imaging diffusion tensor imaging sequence at 3 T during diastasis, peak systole, and strain sweet spots in a midventricular short-axis slice. The same diffusion tensor imaging sequence was repeated ex vivo after arresting the hearts in either a relaxed (KCl-induced) or contracted (BaCl₂ -induced) state. The displacement-encoded imaging with stimulated echoes data were used to strain-correct the in vivo cardiac diffusion tensor imaging in diastole and systole. The orientation of the primary (helix angles) and secondary (E2A) diffusion eigenvectors was compared with and without strain correction and to the strain-free ex vivo data.

RESULTS

Strain correction reduces systolic E2A significantly when compared without strain correction and ex vivo (median absolute E2A = 34.3° versus E2A = 57.1° (P = 0.01), E2A = 60.5° (P = 0.006), respectively). The systolic distribution of E2A without strain correction is closer to the contracted ex vivo distribution than with strain correction, root mean square deviation of 0.027 versus 0.038.

CONCLUSIONS

The current strain-correction model amplifies the contribution of microscopic strain to diffusion resulting in an overcorrection of E2A. Results show that a new model that considers cellular rearrangement is required. Magn Reson Med 79:2205-2215, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Assessment of Myocardial Microstructural Dynamics by In Vivo Diffusion Tensor Cardiac Magnetic Resonance

BACKGROUND

Cardiomyocytes are organized in microstructures termed sheetlets that reorientate during left ventricular thickening. Diffusion tensor cardiac magnetic resonance (DT-CMR) may enable noninvasive interrogation of in vivo cardiac microstructural dynamics. Dilated cardiomyopathy (DCM) is a condition of abnormal myocardium with unknown sheetlet function.

OBJECTIVES

This study sought to validate in vivo DT-CMR measures of cardiac microstructure against histology, characterize microstructural dynamics during left ventricular wall thickening, and apply the technique in hypertrophic cardiomyopathy (HCM) and DCM.

METHODS

In vivo DT-CMR was acquired throughout the cardiac cycle in healthy swine, followed by in situ and ex vivo DT-CMR, then validated against histology. In vivo DT-CMR was performed in 19 control subjects, 19 DCM, and 13 HCM patients.

RESULTS

In swine, a DT-CMR index of sheetlet reorientation (E2A) changed substantially (E2A mobility ~46°). E2A changes correlated with wall thickness changes (in vivo r² = 0.75; in situ r² = 0.89), were consistently observed under all experimental conditions, and accorded closely with histological analyses in both relaxed and contracted states. The potential contribution of cyclical strain effects to in vivo E2A was ~17%. In healthy human control subjects, E2A increased from diastole (18°) to systole (65°; p < 0.001; E2A mobility = 45°). HCM patients showed significantly greater E2A in diastole than control subjects did (48 ; p < 0.001) with impaired E2A mobility (23°; p < 0.001). In DCM, E2A was similar to control subjects in diastole, but systolic values were markedly lower (40° ; p < 0.001) with impaired E2A mobility (20°; p < 0.001).

CONCLUSIONS

Myocardial microstructure dynamics can be characterized by in vivo DT-CMR. Sheetlet function was abnormal in DCM with altered systolic conformation and reduced mobility, contrasting with HCM, which showed reduced mobility with altered diastolic conformation. These novel insights significantly improve understanding of contractile dysfunction at a level of noninvasive interrogation not previously available in humans. (J Am Coll Cardiol 2017;69:661–76) Published by Elsevier on behalf of the American College of Cardiology Foundation.

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

Supplemental Digital Content is available in the text.

Background—

The objective of this study is to assess the dynamic alterations of myocardial microstructure and strain between diastole and systole in patients with dilated cardiomyopathy relative to healthy controls using the magnetic resonance diffusion tensor imaging, myocardial tagging, and biomechanical modeling.

Methods and Results—

Dual heart-phase diffusion tensor imaging was successfully performed in 9 patients and 9 controls. Tagging data were acquired for the diffusion tensor strain correction and cardiac motion analysis. Mean diffusivity, fractional anisotropy, and myocyte aggregate orientations were compared between both cohorts. Cardiac function was assessed by left ventricular ejection fraction, torsion, and strain. Computational modeling was used to study the impact of cardiac shape on fiber reorientation and how fiber orientations affect strain. In patients with dilated cardiomyopathy, a more longitudinal orientation of diastolic myofiber aggregates was measured compared with controls. Although a significant steepening of helix angles (HAs) during contraction was found in the controls, consistent change in HAs during contraction was absent in patients. Left ventricular ejection fraction, cardiac torsion, and strain were significantly lower in the patients compared with controls. Computational modeling revealed that the dilated heart results in reduced HA changes compared with a normal heart. Reduced torsion was found to be exacerbated by steeper HAs.

Conclusions—

Diffusion tensor imaging revealed reduced reorientation of myofiber aggregates during cardiac contraction in patients with dilated cardiomyopathy relative to controls. Left ventricular remodeling seems to be an important factor in the changes to myocyte orientation. Steeper HAs are coupled with a worsening in strain and torsion. Overall, the findings provide new insights into the structural alterations in patients with dilated cardiomyopathy.

Diffusion MRI in the heart

Diffusion MRI provides unique information on the structure, organization, and integrity of the myocardium without the need for exogenous contrast agents. Diffusion MRI in the heart, however, has proven technically challenging because of the intrinsic non-rigid deformation during the cardiac cycle, displacement of the myocardium due to respiratory motion, signal inhomogeneity within the thorax, and short transverse relaxation times. Recently developed accelerated diffusion-weighted MR acquisition sequences combined with advanced post-processing techniques have improved the accuracy and efficiency of diffusion MRI in the myocardium. In this review, we describe the solutions and approaches that have been developed to enable diffusion MRI of the heart in vivo, including a dual-gated stimulated echo approach, a velocity- (M1 ) or an acceleration- (M2 ) compensated pulsed gradient spin echo approach, and the use of principal component analysis filtering. The structure of the myocardium and the application of these techniques in ischemic heart disease are also briefly reviewed. The advent of clinical MR systems with stronger gradients will likely facilitate the translation of cardiac diffusion MRI into clinical use. The addition of diffusion MRI to the well-established set of cardiovascular imaging techniques should lead to new and complementary approaches for the diagnosis and evaluation of patients with heart disease. © 2015 The Authors. NMR in Biomedicine published by John Wiley & Sons Ltd.

An intravoxel oriented flow model for diffusion-weighted imaging of the kidney

By combining intravoxel incoherent motion (IVIM) and diffusion tensor imaging (DTI) we introduce a new diffusion model called intravoxel oriented flow (IVOF) that accounts for anisotropy of diffusion and the flow-related signal. An IVOF model using a simplified apparent flow fraction tensor (IVOFf ) is applied to diffusion-weighted imaging of human kidneys. The kidneys of 13 healthy volunteers were examined on a 3 T scanner. Diffusion-weighted images were acquired with six b values between 0 and 800 s/mm(2) and 30 diffusion directions. Diffusivity and flow fraction were calculated for different diffusion models. The Akaike information criterion was used to compare the model fit of the proposed IVOFf model to IVIM and DTI. In the majority of voxels the proposed IVOFf model with a simplified apparent flow fraction tensor performs better than IVIM and DTI. Mean diffusivity is significantly higher in DTI compared with models that account for the flow-related signal. The fractional anisotropy of diffusion is significantly reduced when flow fraction is considered to be anisotropic. Anisotropy of the apparent flow fraction tensor is significantly higher in the renal medulla than in the cortex region. The IVOFf model describes diffusion-weighted data in the human kidney more accurately than IVIM or DTI. The apparent flow fraction in the kidney proved to be anisotropic.

Diffusion tensor imaging and histology of developing hearts

Diffusion tensor imaging (DTI) has emerged as a promising method for noninvasive quantification of myocardial microstructure. However, the origin and behavior of DTI measurements during myocardial normal development and remodeling remain poorly understood. In this work, conventional and bicompartmental DTI in addition to three-dimensional histological correlation were performed in a sheep model of myocardial development from third trimester to postnatal 5 months of age. Comparing the earliest time points in the third trimester with the postnatal 5 month group, the scalar transverse diffusivities preferentially increased in both left ventricle (LV) and right ventricle (RV): secondary eigenvalues D2 increased by 54% (LV) and 36% (RV), whereas tertiary eigenvalues D3 increased by 85% (LV) and 67% (RV). The longitudinal diffusivity D1 changes were small, which led to a decrease in fractional anisotropy by 41% (LV) and 33% (RV) in 5 month versus fetal hearts. Histological analysis suggested that myocardial development is associated with hyperplasia in the early stages of the third trimester followed by myocyte growth in the later stages up to 5 months of age (increased average myocyte width by 198%, myocyte length by 128%, and decreased nucleus density by 70% between preterm and postnatal 5 month hearts.) In a few histological samples (N = 6), correlations were observed between DTI longitudinal diffusivity and myocyte length (r = 0.86, P < 0.05), and transverse diffusivity and myocyte width (r = 0.96, P < 0.01). Linear regression analysis showed that transverse diffusivities are more affected by changes in myocyte size and nucleus density changes than longitudinal diffusivities, which is consistent with predictions of classical models of diffusion in porous media. Furthermore, primary and secondary DTI eigenvectors during development changed significantly. Collectively, the findings demonstrate a role for DTI to monitor and quantify myocardial development, and potentially cardiac disease. Copyright © 2016 John Wiley & Sons, Ltd.

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

Purpose To develop a clinically feasible whole-heart free-breathing diffusion-tensor (DT) magnetic resonance (MR) imaging approach with an imaging time of approximately 15 minutes to enable three-dimensional (3D) tractography. Materials and Methods The study was compliant with HIPAA and the institutional review board and required written consent from the participants. DT imaging was performed in seven healthy volunteers and three patients with pulmonary hypertension by using a stimulated echo sequence. Twelve contiguous short-axis sections and six four-chamber sections that covered the entire left ventricle were acquired by using simultaneous multisection (SMS) excitation with a blipped-controlled aliasing in parallel imaging readout. Rate 2 and rate 3 SMS excitation was defined as two and three times accelerated in the section axis, respectively. Breath-hold and free-breathing images with and without SMS acceleration were acquired. Diffusion-encoding directions were acquired sequentially, spatiotemporally registered, and retrospectively selected by using an entropy-based approach. Myofiber helix angle, mean diffusivity, fractional anisotropy, and 3D tractograms were analyzed by using paired t tests and analysis of variance. Results No significant differences (P > .63) were seen between breath-hold rate 3 SMS and free-breathing rate 2 SMS excitation in transmural myofiber helix angle, mean diffusivity (mean ± standard deviation, [0.89 ± 0.09] × 10^-3 mm²/sec vs [0.9 ± 0.09] × 10^-3 mm²/sec), or fractional anisotropy (0.43 ± 0.05 vs 0.42 ± 0.06). Three-dimensional tractograms of the left ventricle with no SMS and rate 2 and rate 3 SMS excitation were qualitatively similar. Conclusion Free-breathing DT imaging of the entire human heart can be performed in approximately 15 minutes without section gaps by using SMS excitation with a blipped-controlled aliasing in parallel imaging readout, followed by spatiotemporal registration and entropy-based retrospective image selection. This method may lead to clinical translation of whole-heart DT imaging, enabling broad application in patients with cardiac disease. ^© RSNA, 2016 Online supplemental material is available for this article.

The effects of noise in cardiac diffusion tensor imaging and the benefits of averaging complex data

There is growing interest in cardiac diffusion tensor imaging (cDTI), but, unlike other diffusion MRI applications, there has been little investigation of the effects of noise on the parameters typically derived. One method of mitigating noise floor effects when there are multiple image averages, as in cDTI, is to average the complex rather than the magnitude data, but the phase contains contributions from bulk motion, which must be removed first. The effects of noise on the mean diffusivity (MD), fractional anisotropy (FA), helical angle (HA) and absolute secondary eigenvector angle (E2A) were simulated with various diffusion weightings (b values). The effect of averaging complex versus magnitude images was investigated. In vivo cDTI was performed in 10 healthy subjects with b = 500, 1000, 1500 and 2000 s/mm(2). A technique for removing the motion-induced component of the image phase present in vivo was implemented by subtracting a low-resolution copy of the phase from the original images before averaging the complex images. MD, FA, E2A and the transmural gradient in HA were compared for un-averaged, magnitude- and complex-averaged reconstructions. Simulations demonstrated an over-estimation of FA and MD at low b values and an under-estimation at high b values. The transition is relatively signal-to-noise ratio (SNR) independent and occurs at a higher b value for FA (b = 1000-1250 s/mm(2)) than MD (b ≈ 250 s/mm(2)). E2A is under-estimated at low and high b values with a transition at b ≈ 1000 s/mm(2), whereas the bias in HA is comparatively small. The under-estimation of FA and MD at high b values is caused by noise floor effects, which can be mitigated by averaging the complex data. Understanding the parameters of interest and the effects of noise informs the selection of the optimal b values. When complex data are available, they should be used to maximise the benefit from the acquisition of multiple averages. The combination of complex data is also a valuable step towards segmented acquisitions.

Orientation dependence of microcirculation-induced diffusion signal in anisotropic tissues

PURPOSE

To seek a better understanding of the effect of organized capillary flow on the MR diffusion-weighted signal.

METHODS

A theoretical framework was proposed to describe the diffusion-weighted MR signal, which was then validated both numerically using a realistic model of capillary network and experimentally in an animal model of isolated perfused heart preparation with myocardial blood flow verified by means of direct arterial spin labeling measurements.

RESULTS

Microcirculation in organized tissues gave rise to an MR signal that could be described as a combination of the bi-exponential behavior of conventional intravoxel incoherent motion (IVIM) theory and diffusion tensor imaging (DTI) -like anisotropy of the vascular signal, with the flow-related pseudo diffusivity represented as the linear algebraic product between the encoding directional unit vector and an appropriate tensor entity. Very good agreement between theoretical predictions and both numerical and experimental observations were found.

CONCLUSION

These findings suggest that the DTI formalism of anisotropic spin motion can be incorporated into the classical IVIM theory to describe the MR signal arising from diffusion and microcirculation in organized tissues. Magn Reson Med 76:1252-1262, 2016. © 2015 Wiley Periodicals, Inc.

Spin echo versus stimulated echo diffusion tensor imaging of the in vivo human heart

Purpose

To compare signal‐to‐noise ratio (SNR) efficiency and diffusion tensor metrics of cardiac diffusion tensor mapping using acceleration‐compensated spin‐echo (SE) and stimulated echo acquisition mode (STEAM) imaging.

Methods

Diffusion weighted SE and STEAM sequences were implemented on a clinical 1.5 Tesla MR system. The SNR efficiency of SE and STEAM was measured (b = 50–450 s/mm²) in isotropic agar, anisotropic diffusion phantoms and the in vivo human heart. Diffusion tensor analysis was performed on mean diffusivity, fractional anisotropy, helix and transverse angles.

Results

In the isotropic phantom, the ratio of SNR efficiency for SE versus STEAM, SNR_t(SE/STEAM), was 2.84 ± 0.08 for all tested b‐values. In the anisotropic diffusion phantom the ratio decreased from 2.75 ± 0.05 to 2.20 ± 0.13 with increasing b‐value, similar to the in vivo decrease from 2.91 ± 0.43 to 2.30 ± 0.30. Diffusion tensor analysis revealed reduced deviation of helix angles from a linear transmural model and reduced transverse angle standard deviation for SE compared with STEAM. Mean diffusivity and fractional anisotropy were measured to be statistically different (P < 0.001) between SE and STEAM.

Conclusion

Cardiac DTI using motion‐compensated SE yields a 2.3–2.9× increase in SNR efficiency relative to STEAM and improved accuracy of tensor metrics. The SE method hence presents an attractive alternative to STEAM based approaches. Magn Reson Med 76:862–872, 2016. © 2015 The Authors. Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Heterogeneity of Fractional Anisotropy and Mean Diffusivity Measurements by In Vivo Diffusion Tensor Imaging in Normal Human Hearts

Background

Cardiac diffusion tensor imaging (cDTI) by cardiovascular magnetic resonance has the potential to assess microstructural changes through measures of fractional anisotropy (FA) and mean diffusivity (MD). However, normal variation in regional and transmural FA and MD is not well described.

Methods

Twenty normal subjects were scanned using an optimised cDTI sequence at 3T in systole. FA and MD were quantified in 3 transmural layers and 4 regional myocardial walls.

Results

FA was higher in the mesocardium (0.46 ±0.04) than the endocardium (0.40 ±0.04, p≤0.001) and epicardium (0.39 ±0.04, p≤0.001). On regional analysis, the FA in the septum was greater than the lateral wall (0.44 ±0.03 vs 0.40 ±0.05 p = 0.04). There was a transmural gradient in MD increasing towards the endocardium (epicardium 0.87 ±0.07 vs endocardium 0.91 ±0.08×10^-3 mm²/s, p = 0.04). With the lateral wall (0.87 ± 0.08×10^-3 mm²/s) as the reference, the MD was higher in the anterior wall (0.92 ±0.08×10^-3 mm²/s, p = 0.016) and septum (0.92 ±0.07×10^-3 mm²/s, p = 0.028). Transmurally the signal to noise ratio (SNR) was greatest in the mesocardium (14.5 ±2.5 vs endocardium 13.1 ±2.2, p<0.001; vs epicardium 12.0 ± 2.4, p<0.001) and regionally in the septum (16.0 ±3.4 vs lateral wall 11.5 ± 1.5, p<0.001). Transmural analysis suggested a relative reduction in the rate of change in helical angle (HA) within the mesocardium.

Conclusions

In vivo FA and MD measurements in normal human heart are heterogeneous, varying significantly transmurally and regionally. Contributors to this heterogeneity are many, complex and interactive, but include SNR, variations in cardiac microstructure, partial volume effects and strain. These data indicate that the potential clinical use of FA and MD would require measurement standardisation by myocardial region and layer, unless pathological changes substantially exceed the normal variation identified.

Comparison of diffusion tensor imaging by cardiovascular magnetic resonance and gadolinium enhanced 3D image intensity approaches to investigation of structural anisotropy in explanted rat hearts

BACKGROUND

Cardiovascular magnetic resonance (CMR) can through the two methods 3D FLASH and diffusion tensor imaging (DTI) give complementary information on the local orientations of cardiomyocytes and their laminar arrays.

METHODS

Eight explanted rat hearts were perfused with Gd-DTPA contrast agent and fixative and imaged in a 9.4T magnet by two types of acquisition: 3D fast low angle shot (FLASH) imaging, voxels 50 × 50 × 50 μm, and 3D spin echo DTI with monopolar diffusion gradients of 3.6 ms duration at 11.5 ms separation, voxels 200 × 200 × 200 μm. The sensitivity of each approach to imaging parameters was explored.

RESULTS

The FLASH data showed laminar alignments of voxels with high signal, in keeping with the presumed predominance of contrast in the interstices between sheetlets. It was analysed, using structure-tensor (ST) analysis, to determine the most (v1(ST)), intermediate (v2(ST)) and least (v3(ST)) extended orthogonal directions of signal continuity. The DTI data was analysed to determine the most (e1(DTI)), intermediate (e2(DTI)) and least (e3(DTI)) orthogonal eigenvectors of extent of diffusion. The correspondence between the FLASH and DTI methods was measured and appraised. The most extended direction of FLASH signal (v1(ST)) agreed well with that of diffusion (e1(DTI)) throughout the left ventricle (representative discrepancy in the septum of 13.3 ± 6.7°: median ± absolute deviation) and both were in keeping with the expected local orientations of the long-axis of cardiomyocytes. However, the orientation of the least directions of FLASH signal continuity (v3(ST)) and diffusion (e3(ST)) showed greater discrepancies of up to 27.9 ± 17.4°. Both FLASH (v3(ST)) and DTI (e3(DTI)) where compared to directly measured laminar arrays in the FLASH images. For FLASH the discrepancy between the structure-tensor calculated v3(ST) and the directly measured FLASH laminar array normal was of 9 ± 7° for the lateral wall and 7 ± 9° for the septum (median ± inter quartile range), and for DTI the discrepancy between the calculated v3(DTI) and the directly measured FLASH laminar array normal was 22 ± 14° and 61 ± 53.4°. DTI was relatively insensitive to the number of diffusion directions and to time up to 72 hours post fixation, but was moderately affected by b-value (which was scaled by modifying diffusion gradient pulse strength with fixed gradient pulse separation). Optimal DTI parameters were b = 1000 mm/s(2) and 12 diffusion directions. FLASH acquisitions were relatively insensitive to the image processing parameters explored.

CONCLUSIONS

We show that ST analysis of FLASH is a useful and accurate tool in the measurement of cardiac microstructure. While both FLASH and the DTI approaches appear promising for mapping of the alignments of myocytes throughout myocardium, marked discrepancies between the cross myocyte anisotropies deduced from each method call for consideration of their respective limitations.

In vivo cardiovascular magnetic resonance diffusion tensor imaging shows evidence of abnormal myocardial laminar orientations and mobility in hypertrophic cardiomyopathy

BACKGROUND

Cardiac diffusion tensor imaging (cDTI) measures the magnitudes and directions of intramyocardial water diffusion. Assuming the cross-myocyte components to be constrained by the laminar microstructures of myocardium, we hypothesized that cDTI at two cardiac phases might identify any abnormalities of laminar orientation and mobility in hypertrophic cardiomyopathy (HCM).

METHODS

We performed cDTI in vivo at 3 Tesla at end-systole and late diastole in 11 healthy controls and 11 patients with HCM, as well as late gadolinium enhancement (LGE) for detection of regional fibrosis.

RESULTS

Voxel-wise analysis of diffusion tensors relative to left ventricular coordinates showed expected transmural changes of myocardial helix-angle, with no significant differences between phases or between HCM and control groups. In controls, the angle of the second eigenvector of diffusion (E2A) relative to the local wall tangent plane was larger in systole than diastole, in accord with previously reported changes of laminar orientation. HCM hearts showed higher than normal global E2A in systole (63.9° vs 56.4° controls, p=0.026) and markedly raised E2A in diastole (46.8° vs 24.0° controls, p<0.001). In hypertrophic regions, E2A retained a high, systole-like angulation even in diastole, independent of LGE, while regions of normal wall thickness did not (LGE present 57.8°, p=0.0028, LGE absent 54.8°, p=0.0022 vs normal thickness 38.1°).

CONCLUSIONS

In healthy controls, the angles of cross-myocyte components of diffusion were consistent with previously reported transmural orientations of laminar microstructures and their changes with contraction. In HCM, especially in hypertrophic regions, they were consistent with hypercontraction in systole and failure of relaxation in diastole. Further investigation of this finding is required as previously postulated effects of strain might be a confounding factor.