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

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.

Simultaneous non-contrast assessment of cardiac microstructure and perfusion in vivo in the human heart

BACKGROUND

Intravoxel incoherent motion (IVIM) imaging can provide information on cardiac microstructure and microvascular perfusion from a single examination. However, the spin echo-based approaches typically used for cardiac IVIM suffer from low sensitivity to changes in perfusion. The aim of this work was to develop a stimulated-echo (STEAM)-based method for IVIM and diffusion tensor cardiovascular magnetic resonance to simultaneously provide biomarkers of microstructure and perfusion in vivo in the human heart.

METHODS

Here we introduce a novel STEAM-IVIM sequence incorporating phase cycling to obtain true non-diffusion weighted images (b = 0 s/mm2). STEAM-IVIM imaging was performed at 20 b-values (0 to 1000 s/mm2) to enable accurate estimation of the IVIM parameters, and with six diffusion encoding directions to enable reconstruction of the diffusion tensor. 20 healthy subjects (8 female, median age 31 years) were imaged on a clinical 3T system with STEAM-IVIM. A simulation study was performed to investigate the optimal fitting algorithms for the IVIM parameters, which was subsequently used to create pixel-wise IVIM parameter maps for the in vivo acquisitions.

RESULTS

Good image quality across the myocardium was obtained for all b-values. Mean(±SD) IVIM parameter estimates were: diffusivity D = 0.83 ± 0.07 × 10-3 mm2/s, perfusion coefficient D* = 19.08 ± 6.48 × 10-3 mm2/s, perfusion fraction f = 19.72 ± 4.11%, and mean diffusion tensor parameters were: mean diffusivity = 0.88 ± 0.06 × 10-3 mm2/s, fractional anisotropy = 0.45 ± 0.04, absolute E2 angle = 55.29 ± 6.38º, helix angle gradient = -0.68 ± 0.18º/%.

CONCLUSION

Phase-cycled STEAM-IVIM enables fitting of cardiac diffusion tensor and perfusion parameters in healthy subjects and shows promise for the simultaneous detection of microstructural aberration and perfusion abnormalities in the presence of cardiac disease without the need for exogenous contrast agents.

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.

CMRsim-A python package for cardiovascular MR simulations incorporating complex motion and flow

PURPOSE

To present an open-source MR simulation framework that facilitates the incorporation of complex motion and flow for studying cardiovascular MR (CMR) acquisition and reconstruction.

METHODS

CMRsim is a Python package that allows simulation of CMR images using dynamic digital phantoms with complex motion as input. Two simulation paradigms are available, namely, numerical and analytical solutions to the Bloch equations, using a common motion representation. Competitive simulation speeds are achieved using TensorFlow for GPU acceleration. To demonstrate the capability of the package, one introductory and two advanced CMR simulation experiments are presented. The latter showcase phase-contrast imaging of turbulent flow downstream of a stenotic section and cardiac diffusion tensor imaging on a contracting left ventricle. Additionally, extensive documentation and example resources are provided.

RESULTS

The Bloch simulation with turbulent flow using approximately 1.5 million particles and a sequence duration of 710 ms for each of the seven different velocity encodings took a total of 29 min on a NVIDIA Titan RTX GPU. The results show characteristic phase contrast and magnitude modulation present in real data. The analytical simulation of cardiac diffusion tensor imaging with bulk-motion phase sensitivity took approximately 10 s per diffusion-weighted image, including preparation and loading steps. The results exhibit the expected alteration of diffusion metrics due to strain.

CONCLUSION

CMRsim is the first simulation framework that allows one to feasibly incorporate complex motion, including turbulent flow, to systematically study advanced CMR acquisition and reconstruction approaches. The open-source package features modularity and transparency, facilitating maintainability and extensibility in support of reproducible research.