Dataset Information

Decompose quantitative susceptibility mapping (QSM) to sub-voxel diamagnetic and paramagnetic components based on gradient-echo MRI data.

ABSTRACT:

Purpose

A method named DECOMPOSE-QSM is developed to decompose bulk susceptibility measured with QSM into sub-voxel paramagnetic and diamagnetic components based on a three-pool complex signal model.

Methods

Multi-echo gradient echo signal is modeled as a summation of three weighted exponentials corresponding to three types of susceptibility sources: reference susceptibility, diamagnetic and paramagnetic susceptibility relative to the reference. Paramagnetic component susceptibility (PCS) and diamagnetic component susceptibility (DCS) maps are constructed to represent the sub-voxel compartments by solving for linear and nonlinear parameters in the model.

Results

Numerical forward simulation and phantom validation confirmed the ability of DECOMPOSE-QSM to separate the mixture of paramagnetic and diamagnetic components. The PCS obtained from temperature-variant brainstem imaging follows the Curie's Law, which further validated the model and the solver. Initial in vivo investigation of human brain images showed the ability to extract sub-voxel PCS and DCS sources that produce visually enhanced contrast between brain structures comparing to threshold QSM.

SUBMITTER: Chen J

PROVIDER: S-EPMC8720043 | biostudies-literature | 2021 Nov

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Decompose quantitative susceptibility mapping (QSM) to sub-voxel diamagnetic and paramagnetic components based on gradient-echo MRI data.

Chen Jingjia J Gong Nan-Jie NJ Chaim Khallil Taverna KT Otaduy Maria Concepción García MCG Liu Chunlei C

NeuroImage 20210814

<h4>Purpose</h4>A method named DECOMPOSE-QSM is developed to decompose bulk susceptibility measured with QSM into sub-voxel paramagnetic and diamagnetic components based on a three-pool complex signal model.<h4>Methods</h4>Multi-echo gradient echo signal is modeled as a summation of three weighted exponentials corresponding to three types of susceptibility sources: reference susceptibility, diamagnetic and paramagnetic susceptibility relative to the reference. Paramagnetic component susceptibili ...[more]

PMID: 34403742

Similar Datasets

Project description:PurposeTo develop a method to reconstruct quantitative susceptibility mapping (QSM) from multi-echo, multi-flip angle data collected using strategically acquired gradient echo (STAGE) imaging.MethodsThe proposed QSM reconstruction algorithm, referred to as "structurally constrained Susceptibility Weighted Imaging and Mapping" scSWIM, performs an ℓ 1 and ℓ 2 regularization-based reconstruction in a single step. The unique contrast of the T1 weighted enhanced (T1WE) image derived from STAGE imaging was used to extract reliable geometry constraints to protect the basal ganglia from over-smoothing. The multi-echo multi-flip angle data were used for improving the contrast-to-noise ratio in QSM through a weighted averaging scheme. The measured susceptibility values from scSWIM for both simulated and in vivo data were compared to the: original susceptibility model (for simulated data only), the multi orientation COSMOS (for in vivo data only), truncated k-space division (TKD), iterative susceptibility weighted imaging and mapping (iSWIM), and morphology enabled dipole inversion (MEDI) algorithms. Goodness of fit was quantified by measuring the root mean squared error (RMSE) and structural similarity index (SSIM). Additionally, scSWIM was assessed in ten healthy subjects.ResultsThe unique contrast and tissue boundaries from T1WE and iSWIM enable the accurate definition of edges of high susceptibility regions. For the simulated brain model without the addition of microbleeds and calcium, the RMSE was best at 5.21ppb for scSWIM and 8.74ppb for MEDI thanks to the reduced streaking artifacts. However, by adding the microbleeds and calcium, MEDI's performance dropped to 47.53ppb while scSWIM performance remained the same. The SSIM was highest for scSWIM (0.90) and then MEDI (0.80). The deviation from the expected susceptibility in deep gray matter structures for simulated data relative to the model (and for the in vivo data relative to COSMOS) as measured by the slope was lowest for scSWIM + 1%(-1%); MEDI + 2%(-11%) and then iSWIM -5%(-10%). Finally, scSWIM measurements in the basal ganglia of healthy subjects were in agreement with literature.ConclusionThis study shows that using a data fidelity term and structural constraints results in reduced noise and streaking artifacts while preserving structural details. Furthermore, the use of STAGE imaging with multi-echo and multi-flip data helps to improve the signal-to-noise ratio in QSM data and yields less artifacts.

Project description:Background Myocardial oxygenation imaging could help determine the presence of microvascular dysfunction associated with increased cardiovascular risk. However, it is challenging to depict the potentially small oxygenation alterations with current noninvasive cardiac MRI blood oxygen level-dependent (BOLD) techniques. Purpose To demonstrate the cardiac application of a gradient-echo spin-echo (GESE) echo-planar imaging sequence for dynamic and quantitative heartbeat-to-heartbeat BOLD MRI and evaluate the sequence in populations both healthy and with hypertension in combination with a breath hold-induced CO2 intervention. Materials and Methods GESE echo-planar imaging sequence was performed in 18 healthy participants and in eight prospectively recruited participants with hypertension on a 3.0-T MRI system. T2 and T2* maps were calculated per heartbeat with a four-parameter fitting technique. Septal regions of interests were used to determine T2 and T2* values per heartbeat and examined over the course of a breath hold to determine BOLD changes. T2 and T2* changes of healthy participants and participants with hypertension were compared by using a nonparametric Mann-Whitney test. Results GESE echo-planar imaging approach gave spatially stable T2 and T2* maps per heartbeat for healthy participants and participants with hypertension, with mean T2 values of 43 msec ± 5 (standard deviation) and 46 msec ± 9, respectively, and mean T2* values of 28 msec ± 5 and 22 msec ± 5, respectively. The healthy participants exhibited increasing T2 and T2* values over the course of a breath hold with a mean positive slope of 0.2 msec per heartbeat ± 0.1 for T2 and 0.2 msec per heartbeat ± 0.1 for T2*, whereas for participants with hypertension these dynamic T2 and T2* values had a mean negative slope of -0.2 msec per heartbeat ± 0.2 for T2 and -0.1 msec per heartbeat ± 0.2 for T2*. The difference in these mean slopes between healthy participants and participants with hypertension was significant for both T2 (P < .001) and T2* (P < .001). Conclusion Gradient-echo spin-echo echo-planar imaging sequence provided quantitative T2 and T2* maps per heartbeat and enabled dynamic heartbeat-to-heartbeat blood oxygen level-dependent (BOLD)-response imaging by analyzing changes in T2 and T2* over the time of a breath-hold intervention. This approach could identify differences in the BOLD response between healthy participants and participants with hypertension. © RSNA, 2020 Online supplemental material is available for this article. See also the editorial by Friedrich in this issue.

Dataset Information

Decompose quantitative susceptibility mapping (QSM) to sub-voxel diamagnetic and paramagnetic components based on gradient-echo MRI data.

Purpose

Methods

Results

Publications

Decompose quantitative susceptibility mapping (QSM) to sub-voxel diamagnetic and paramagnetic components based on gradient-echo MRI data.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets