Decompose quantitative susceptibility mapping (QSM) to sub-voxel diamagnetic and paramagnetic components based on gradient-echo MRI data.
Ontology highlight
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 |
REPOSITORIES: biostudies-literature
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.
| S-EPMC7645168 | biostudies-literature
Project description:Recent studies have shown that there is a direct link between the orientation of the nerve fibers in white matter (WM) and the contrast observed in magnitude and phase images acquired using gradient echo MRI. Understanding the origin of this link is of great interest because it could offer access to a new diagnostic tool for investigating tissue microstructure. Since it has been suggested that myelin is the dominant source of this contrast, creating an accurate model for characterizing the effect of the myelin sheath on the evolution of the NMR signal is an essential step toward fully understanding WM contrast. In this study, we show by comparison of the results of simulations and experiments carried out on human subjects at 7T, that the magnitude and phase of signals acquired from WM in vivo can be accurately characterized by (i) modeling the myelin sheath as a hollow cylinder composed of material having an anisotropic magnetic susceptibility that is described by a tensor with a radially oriented principal axis, and (ii) adopting a two-pool model in which the water in the sheath has a reduced T(2) relaxation time and spin density relative to its surroundings, and also undergoes exchange. The accuracy and intrinsic simplicity of the hollow cylinder model provides a versatile framework for future exploitation of the effect of WM microstructure on gradient echo contrast in clinical MRI.
| S-EPMC3494918 | biostudies-literature
Project description:In MRI, the main magnetic field polarizes the electron cloud of a molecule, generating a chemical shift for observer protons within the molecule and a magnetic susceptibility inhomogeneity field for observer protons outside the molecule. The number of water protons surrounding a molecule for detecting its magnetic susceptibility is vastly greater than the number of protons within the molecule for detecting its chemical shift. However, the study of tissue magnetic susceptibility has been hindered by poor molecular specificities of hitherto used methods based on MRI signal phase and T2* contrast, which depend convolutedly on surrounding susceptibility sources. Deconvolution of the MRI signal phase can determine tissue susceptibility but is challenged by the lack of MRI signal in the background and by the zeroes in the dipole kernel. Recently, physically meaningful regularizations, including the Bayesian approach, have been developed to enable accurate quantitative susceptibility mapping (QSM) for studying iron distribution, metabolic oxygen consumption, blood degradation, calcification, demyelination, and other pathophysiological susceptibility changes, as well as contrast agent biodistribution in MRI. This paper attempts to summarize the basic physical concepts and essential algorithmic steps in QSM, to describe clinical and technical issues under active development, and to provide references, codes, and testing data for readers interested in QSM.
| S-EPMC4297605 | biostudies-literature
Project description:PURPOSE:Quantitative susceptibility mapping (QSM) MRI is a tool that can characterize changes in susceptibility, an intrinsic property which is associated with compositional changes in the tissue. Current QSM estimation of cortical bone is challenging because conventional clinical MRI cannot acquire signal in cortical bone. This study aimed to implement Cones 3D ultrashort echo time MRI (UTE-MRI) for ex vivo QSM measurements in human tibial cortical bone, investigating the correlations of QSM with volumetric intracortical bone mineral density (BMD). MATERIALS AND METHODS:Nine tibial midshaft cortical bone specimens (25 mm long specimens cut at the mid-point of tibial shaft, 67 ± 20 years old, 5 women and 4 men) were scanned on a clinical 3 T MRI scanner for QSM measurement. The specimens were also scanned on a high-resolution micro-computed tomography (μCT) scanner for volumetric BMD estimation. QSM and μCT results were compared at approximately nine regions of interest (ROIs) per specimen. RESULTS:Average 3D UTE-MRI QSM showed significantly strong correlation with volumetric BMD (R = -0.82, P < 0.01) and bone porosity (R = 0.72, P < 0.01). Combining all data points together (77 ROIs), QSM showed significant moderate to strong correlation with volumetric BMD after correction for interdependencies in specimens (R = -0.70, P < 0.01). The corrections were required because the data points were not independent in each specimen. Similarly, the correlation between QSM and porosity was significant (R = 0.68, P < 0.01). CONCLUSIONS:These results suggest that the Cones 3D UTE-MRI QSM technique can potentially serve as a novel and accurate tool to assess intracortical bone mineral density whilst avoiding ionizing radiation.
| S-EPMC6689249 | biostudies-literature
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.
| S-EPMC7053244 | biostudies-literature
Project description:In this paper we demonstrate that gradient echo MRI with hyperpolarized (3)He gas can be used for simultaneously extracting in vivo information about lung ventilation properties, alveolar geometrical parameters, and blood vessel network structure. This new approach is based on multi-gradient-echo experimental measurements of hyperpolarized (3)He gas MRI signal from human lungs and a proposed theoretical model of this signal. Based on computer simulations of (3)He atoms diffusing in the acinar airway tree in the presence of an inhomogeneous magnetic field induced by the susceptibility differences between lung tissue (alveolar septa, blood vessels) and lung airspaces, we derive analytical expressions relating the time-dependent MR signal to the geometrical parameters of acinar airways and the blood vessel network. Data obtained on eight healthy volunteers are in good agreement with literature values. This information is complementary to the information obtained by means of the in vivo lung morphometry technique with hyperpolarized 3He diffusion MRI previously developed by our group, and opens new opportunities to study lung microstructure in health and disease.
| S-EPMC4232999 | biostudies-other
Project description:PURPOSE:To evaluate the echo dependence of 3D ultrashort echo time (TE) quantitative susceptibility mapping (3D UTE-QSM) and effective transverse relaxation rate ( R2*) measurement in the setting of high concentrations of iron oxide nanoparticles. METHODS:A phantom study with iron concentrations ranging from 2 to 22?mM was performed using a 3D UTE Cones sequence. Simultaneous QSM processing with morphology-enabled dipole inversion (MEDI) and R2* single exponential fitting was conducted offline with the acquired 3D UTE data. The dependence of UTE-QSM and R2* on echo spacing (?TE) and first TE (TE1 ) was investigated. RESULTS:A linear relationship was observed between UTE-QSM measurement and iron concentration up to 22?mM only, with the minimal TE1 of 0.032 ms and ?TE of less than 0.1 ms. A linear relationship was observed between R2* and iron concentration up to 22?mM only when TE1 was less than 0.132 ms and ?TE was less than 1.2 ms. UTE-QSM with MEDI processing showed strong dependence on ?TE and TE1 , especially at high iron concentrations. CONCLUSION:UTE-QSM is more sensitive than R2* measurement to TE selection. Both an ultrashort TE1 and a small ?TE are needed to achieve accurate QSM for high iron concentrations. Magn Reson Med 79:2315-2322, 2018. © 2018 International Society for Magnetic Resonance in Medicine.
| S-EPMC5811338 | biostudies-literature
Project description:PurposeTo compare cortical gray matter oxygen extraction fraction (OEF) estimated from 2 MRI methods: (1) the quantitative susceptibility mapping (QSM) plus quantitative blood oxygen level dependent imaging (qBOLD) (QSM+qBOLD or QQ), and (2) the dual-gas calibrated-BOLD (DGCB) in healthy subjects; and to investigate the validity of iso-cerebral metabolic rate of oxygen consumption assumption during hypercapnia using QQ.MethodsIn 10 healthy subjects, 3 tesla MRI including a multi-echo gradient echo sequence at baseline and hypercapnia for QQ, as well as an EPI dual-echo pseudo-continuous arterial spin labeling for DGCB, were performed under a hypercapnic and a hyperoxic condition. OEFs from QQ and DGCB were compared using region of interest analysis and paired t test. For QQ, cerebral metabolic rate of oxygen consumption = cerebral blood flow*OEF*arterial oxygen content was generated for both baseline and hypercapnia, which were compared.ResultsAverage OEF in cortical gray matter across 10 subjects from QQ versus DGCB was 35.5 ± 6.7% versus 38.0 ± 9.1% (P = .49) at baseline and 20.7 ± 4.4% versus 28.4 ± 7.6% (P = .02) in hypercapnia: OEF in cortical gray matter was significantly reduced as measured in QQ (P < .01) and in DGCB (P < .01). Cerebral metabolic rate of oxygen consumption (in μmol O2 /min/100 g) was 168.2 ± 54.1 at baseline from DGCB and was 153.1 ± 33.8 at baseline and 126.4 ± 34.2 (P < .01) in hypercapnia from QQ.ConclusionThe differences in OEF obtained from QQ and DGCB are small and nonsignificant at baseline but are statistically significant during hypercapnia. In addition, QQ shows a cerebral metabolic rate of oxygen consumption decrease (17.4%) during hypercapnia.
| S-EPMC7722021 | biostudies-literature
Project description:This work introduces a novel algorithm for deconvolution of the BOLD signal in multi-echo fMRI data: Multi-echo Sparse Paradigm Free Mapping (ME-SPFM). Assuming a linear dependence of the BOLD percent signal change on the echo time (TE) and using sparsity-promoting regularized least squares estimation, ME-SPFM yields voxelwise time-varying estimates of the changes in the apparent transverse relaxation (ΔR2⁎) without prior knowledge of the timings of individual BOLD events. Our results in multi-echo fMRI data collected during a multi-task event-related paradigm at 3 Tesla demonstrate that the maps of R2⁎ changes obtained with ME-SPFM at the times of the stimulus trials show high spatial and temporal concordance with the activation maps and BOLD signals obtained with standard model-based analysis. This method yields estimates of ΔR2⁎ having physiologically plausible values. Owing to its ability to blindly detect events, ME-SPFM also enables us to map ΔR2⁎ associated with spontaneous, transient BOLD responses occurring between trials. This framework is a step towards deciphering the dynamic nature of brain activity in naturalistic paradigms, resting-state or experimental paradigms with unknown timing of the BOLD events.
| S-EPMC6819276 | biostudies-literature
Project description:Quantitative susceptibility mapping (QSM) is increasingly used to measure variation in tissue composition both in the brain and in other areas of the body in a range of disease pathologies. Although QSM measurements were originally believed to be independent of the echo time (TE) used in the gradient-recalled echo (GRE) acquisition from which they are derived; recent literature (Sood et al., 2016) has shown that these measurements can be highly TE-dependent in a number of brain regions. In this work we systematically investigate possible causes of this effect through analysis of apparent frequency and QSM measurements derived from data acquired at multiple TEs in vivo in healthy brain regions and in cerebral microbleeds (CMBs); QSM data acquired in a gadolinium-doped phantom; and in QSM data derived from idealized simulated phase data. Apparent frequency measurements in the optic radiations (OR) and central corpus callosum (CC) were compared to those predicted by a 3-pool white matter model, however the model failed to fully explain contrasting frequency profiles measured in the OR and CC. Our results show that TE-dependent QSM measurements can be caused by a failure of phase unwrapping algorithms in and around strong susceptibility sources such as CMBs; however, in healthy brain regions this behavior appears to result from intrinsic non-linear phase evolution in the MR signal. From these results we conclude that care must be taken when deriving frequency and QSM measurements in strong susceptibility sources due to the inherent limitations in phase unwrapping; and that while signal compartmentalization due to tissue microstructure and content is a plausible cause of TE-dependent frequency and QSM measurements in healthy brain regions, better sampling of the MR signal and more complex models of tissue are needed to fully exploit this relationship.
| S-EPMC5368016 | biostudies-literature