Evaluation of non-Gaussian diffusion in cardiac MRI

Cardiac q-space trajectory imaging by motion-compensated tensor-valued diffusion encoding in human heart in vivo

PURPOSE

Tensor-valued diffusion encoding can probe more specific features of tissue microstructure than what is available by conventional diffusion weighting. In this work, we investigate the technical feasibility of tensor-valued diffusion encoding at high b-values with q-space trajectory imaging (QTI) analysis, in the human heart in vivo.

METHODS

Ten healthy volunteers were scanned on a 3T scanner. We designed time-optimal gradient waveforms for tensor-valued diffusion encoding (linear and planar) with second-order motion compensation. Data were analyzed with QTI. Normal values and repeatability were investigated for the mean diffusivity (MD), fractional anisotropy (FA), microscopic FA (μFA), isotropic, anisotropic and total mean kurtosis (MKi, MKa, and MKt), and orientation coherence (Cc ). A phantom, consisting of two fiber blocks at adjustable angles, was used to evaluate sensitivity of parameters to orientation dispersion and diffusion time.

RESULTS

QTI data in the left ventricular myocardium were MD = 1.62 ± 0.07 μm2 /ms, FA = 0.31 ± 0.03, μFA = 0.43 ± 0.07, MKa = 0.20 ± 0.07, MKi = 0.13 ± 0.03, MKt = 0.33 ± 0.09, and Cc  = 0.56 ± 0.22 (mean ± SD across subjects). Phantom experiments showed that FA depends on orientation dispersion, whereas μFA was insensitive to this effect.

CONCLUSION

We demonstrated the first tensor-valued diffusion encoding and QTI analysis in the heart in vivo, along with first measurements of myocardial μFA, MKi, MKa, and Cc . The methodology is technically feasible and provides promising novel biomarkers for myocardial tissue characterization.

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.

Assessing Myocardial Microstructure With Biophysical Models of Diffusion MRI

Biophysical models are a promising means for interpreting diffusion weighted magnetic resonance imaging (DW-MRI) data, as they can provide estimates of physiologically relevant parameters of microstructure including cell size, volume fraction, or dispersion. However, their application in cardiac microstructure mapping (CMM) has been limited. This study proposes seven new two-compartment models with combination of restricted cylinder models and a diffusion tensor to represent intra- and extracellular spaces, respectively. Three extended versions of the cylinder model are studied here: cylinder with elliptical cross section (ECS), cylinder with Gamma distributed radii (GDR), and cylinder with Bingham distributed axes (BDA). The proposed models were applied to data in two fixed mouse hearts, acquired with multiple diffusion times, q-shells and diffusion encoding directions. The cylinderGDR-pancake model provided the best performance in terms of root mean squared error (RMSE) reducing it by 25% compared to diffusion tensor imaging (DTI). The cylinderBDA-pancake model represented anatomical findings closest as it also allows for modelling dispersion. High-resolution 3D synchrotron X-ray imaging (SRI) data from the same specimen was utilized to evaluate the biophysical models. A novel tensor-based registration method is proposed to align SRI structure tensors to the MR diffusion tensors. The consistency between SRI and DW-MRI parameters demonstrates the potential of compartment models in assessing physiologically relevant parameters.

MRI phantom for tissue simulation with respect to diffusion coefficient and kurtosis - Validation with injection of liposomal theranostics

This study presents gelatine-based and agar-based phantoms with an addition of glycerol, safflower oil, silicone oil and cellulose microcrystalline with a potential to cover the entire range of tissue diffusion coefficients and kurtosis values. Forty types of phantoms were prepared and examined for NMR relaxation times T₁ and T₂ and diffusional metrics D, K and ADC. Wide ranges of values of D (0.0003-0.0031 mm²s^-1), K (0.00-7.24) and ADC (0.0002-0.0031 mm²s^-1) were observed. Two of the phantoms closely mimic muscle and cortical gray matter with respect to water diffusion parameters. Although many of the presented phantoms display both D and K values within the range of human tissues, they match different tissues with respect to D and K. The imaging results for the gray matter simulating phantom injected with the liposomal solution, bear a resemblance to the particle size effect described in the literature. The phantoms presented in this work are simple in preparation and affordable tissue-simulating materials to be used primarily in development of diffusion kurtosis-based MRI methods and possibly in a preliminary assessment of MRI contrast agents. Further adjustments of the chemical compositions could potentially lead to development of new types of phantoms mimicking diffusional properties of more kinds of soft tissues.

Diffusion-weighted imaging in hypertrophic cardiomyopathy: association with high T2-weighted signal intensity in addition to late gadolinium enhancement

Diffusion-weighted imaging (DWI) has been confirmed to be associated with late gadolinium enhancement (LGE) in hypertrophic cardiomyopathy (HCM). In this context, we aimed to study whether DWI could reflect the active tissue injury and edema information of HCM which were usually indicated by T2 weighted images. Forty HCM patients were examined using a 3.0 T magnetic resonance scanner. Cine, T2-weighted short tau inversion recovery (T2-STIR), DWI and LGE sequences were acquired. T1 mapping was also included to quantify the focal and diffuse fibrosis. Cardiac troponin I (cTnI) was tested to assess the recently myocardial injury. Student's t-test, Mann-Whitney U test, One-way analysis, Kruskal-Wallis analysis, the Spearman correlation analysis, and multivariable regression were used in this study. The apparent diffusion coefficient (ADC) was significantly elevated in the cTnI positive group (P = 0.01) and correlated with LGE (ρ = 0.312, P = 0.049) and HighT2 extent (ρ = 0.443, P = 0.004) in the global level. In the segmental analysis, the ADC significantly differentiated the segments with and without HighT2/LGE presence (P = 0.00). The average ADC values were higher in segments with HighT2 and LGE coexistence than in those with only LGE presence (P < 0.05). Multivariable regression indicated that segmental HighT2 and LGE were both contributing factors to the ADC values. In this study of HCM, we confirmed that ADC as a molecular diffusion parameter reflects the replacement fibrosis of myocardium. Moreover, it also reveals edema extent and its association with serum cTnI.

Regional variations in ex-vivo diffusion tensor anisotropy are associated with cardiomyocyte remodeling in rats after left ventricular pressure overload

BACKGROUND

Pressure overload left ventricular (LV) hypertrophy is characterized by increased cardiomyocyte width and ventricle wall thickness, however the regional variation of this remodeling is unclear. Cardiovascular magnetic resonance (CMR) diffusion tensor imaging (DTI) may provide a non-invasive, comprehensive, and geometrically accurate method to detect regional differences in structural remodeling in hypertrophy. We hypothesized that DTI parameters, such as fractional and planar anisotropy, would reflect myocyte remodeling due to pressure overload in a regionally-dependent manner.

METHODS

We investigated the regional distributions of myocyte remodeling in rats with or without transverse aortic constriction (TAC) via direct measurement of myocyte dimensions with confocal imaging of thick tissue sections, and correlated myocyte cross-sectional area and other geometric features with parameters of diffusivity from ex-vivo DTI in the same regions of the same hearts.

RESULTS

We observed regional differences in several parameters from DTI between TAC hearts and SHAM controls. Consistent with previous studies, helix angles from DTI correlated strongly with those measured directly from histological sections (p < 0.001, R2 = 0.71). There was a transmural gradient in myocyte cross-sectional area in SHAM hearts that was diminished in the TAC group. We also found several regions of significantly altered DTI parameters in TAC LV compared to SHAM, especially in myocyte sheet angle dispersion and planar anisotropy. Among others, these parameters correlated significantly with directly measured myocyte aspect ratios.

CONCLUSIONS

These results show that structural remodeling in pressure overload LV hypertrophy is regionally heterogeneous, especially transmurally, with a greater degree of remodeling in the sub-endocardium compared to the sub-epicardium. Additionally, several parameters derived from DTI correlated significantly with measurements of myocyte geometry from direct measurement in histological sections. We suggest that DTI may provide a non-invasive, comprehensive method to detect regional structural myocyte LV remodeling during disease.

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

Purpose

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

Methods

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

Results

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

Conclusion

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

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.

Diffusion kurtosis imaging provides quantitative assessment of the microstructure changes of disc degeneration: an in vivo experimental study

OBJECTIVE

Our aim was to assess the microstructural changes of intervertebral disc degeneration induced by annulus needle puncture in rats by diffusion kurtosis imaging (DKI).

METHODS

Eighteen rats (36 discs) were punctured percutaneously at the intervertebral disc between C6/7, C7/8 (C-coccygeal vertebrae) with a 21-gauge needle. The rats were divided into six groups according to the time after the puncture: 3 h, 48 h, 3 days, 7 days, 10 days and 14 days. There were six discs in three rats in the control group. The rats' tail was imaged at 3T MRI with T₂-weighted and diffusion-weighted and diffusion kurtosis imaging (DWI)/DKI sequences. The discs were categorized using a five-grade degeneration system based on the T₂ images. The height of the discs and the parameters in DWI/DKI were measured and compared between the different time points. The histological images were also obtained from the discs.

RESULTS

The histological study revealed that the discs in the rat of the punctured groups were degenerated. The T2 grades of different groups presented an increasing trend from 7 to 10 days after puncture (R² = 0.9424, P < 0.001), while the DWI/DKI parameters changes were consistent with the histological changes at the different time points and showed significant differences between the different groups (P < 0.05).

CONCLUSIONS

DKI provides quantitative assessment of the microstructure changes of disc degeneration, and it is a non-invasive method. The DKI multi-parameter analysis is sensitive to discs changes caused by puncture. These slides can be retrieved under Electronic Supplementary Material.