Evaluating the Accuracy of Susceptibility Maps Calculated from Single-echo versus Multi-echo Gradient-echo Acquisitions

For Susceptibility Mapping (SM), Laplacian-based methods (LBMs) can be used on single- or multi-echo gradient echo phase data. Previous studies have shown the advantage of using multi-echo versus single-echo data for noise reduction in susceptibility-weighted images and simulated data. Here, using simulated and acquired images, we compared the performance of two SM pipelines that used multi- or single-echo phase data and LBMs. We showed that the pipeline that fits the multi-echo data over time first and then applies LBMs gives more accurate local fields and $\chi$ maps than the pipelines that apply LBMs to single-echo phase data

Assessing white matter microstructural changes in idiopathic normal pressure hydrocephalus using voxel-based R2* relaxometry analysis

BackgroundR2* relaxometry and quantitative susceptibility mapping can be combined to distinguish between microstructural changes and iron deposition in white matter. Here, we aimed to explore microstructural changes in the white matter associated with clinical presentations such as cognitive impairment in patients with idiopathic normal-pressure hydrocephalus (iNPH) using R2* relaxometry analysis in combination with quantitative susceptibility mapping.MethodsWe evaluated 16 patients clinically diagnosed with possible or probable iNPH and 18 matched healthy controls (HC) who were chosen based on similarity in age and sex. R2* and quantitative susceptibility mapping were compared using voxel-wise and atlas-based one-way analysis of covariance (ANCOVA). Finally, partial correlation analyses were performed to assess the relationship between R2* and clinical presentations.ResultsR2* was lower in some white matter regions, including the bilateral superior longitudinal fascicle and sagittal stratum, in the iNPH group compared to the HC group. The voxel-based quantitative susceptibility mapping results did not differ between the groups. The atlas-based group comparisons yielded negative mean susceptibility values in almost all brain regions, indicating no clear paramagnetic iron deposition in the white matter of any subject. R2* and cognitive performance scores between the left superior longitudinal fasciculus (SLF) and right sagittal stratum (SS) were positively correlated. In addition to that, R2* and gait disturbance scores between left SS were negatively correlated.ConclusionOur analysis highlights the microstructural changes without iron deposition in the SLF and SS, and their association with cognitive impairment and gait disturbance in patients with iNPH

Magnetic susceptibility anisotropy of myocardium imaged by cardiovascular magnetic resonance reflects the anisotropy of myocardial filament α-helix polypeptide bonds.

BackgroundA key component of evaluating myocardial tissue function is the assessment of myofiber organization and structure. Studies suggest that striated muscle fibers are magnetically anisotropic, which, if measurable in the heart, may provide a tool to assess myocardial microstructure and function.MethodsTo determine whether this weak anisotropy is observable and spatially quantifiable with cardiovascular magnetic resonance (CMR), both gradient-echo and diffusion-weighted data were collected from intact mouse heart specimens at 9.4 Tesla. Susceptibility anisotropy was experimentally calculated using a voxelwise analysis of myocardial tissue susceptibility as a function of myofiber angle. A myocardial tissue simulation was developed to evaluate the role of the known diamagnetic anisotropy of the peptide bond in the observed susceptibility contrast.ResultsThe CMR data revealed that myocardial tissue fibers that were parallel and perpendicular to the magnetic field direction appeared relatively paramagnetic and diamagnetic, respectively. A linear relationship was found between the magnetic susceptibility of the myocardial tissue and the squared sine of the myofiber angle with respect to the field direction. The multi-filament model simulation yielded susceptibility anisotropy values that reflected those found in the experimental data, and were consistent that this anisotropy decreased as the echo time increased.ConclusionsThough other sources of susceptibility anisotropy in myocardium may exist, the arrangement of peptide bonds in the myofilaments is a significant, and likely the most dominant source of susceptibility anisotropy. This anisotropy can be further exploited to probe the integrity and organization of myofibers in both healthy and diseased heart tissue

Investigating the accuracy and precision of TE‐dependent versus multi‐echo QSM using Laplacian‐based methods at 3 T

Purpose: Multi‐echo gradient‐recalled echo acquisitions for QSM enable optimizing the SNR for several tissue types through multi‐echo (TE) combination or investigating temporal variations in the susceptibility (potentially reflecting tissue microstructure) by calculating one QSM image at each TE (TE‐dependent QSM). In contrast with multi‐echo QSM, applying Laplacian‐based methods (LBMs) for phase unwrapping and background field removal to single TEs could introduce nonlinear temporal variations (independent of tissue microstructure) into the measured susceptibility. Here, we aimed to compare the effect of LBMs on the QSM susceptibilities in TE‐dependent versus multi‐echo QSM. Methods: TE–dependent recalled echo data simulated in a numerical head phantom and gradient‐recalled echo images acquired at 3 T in 10 healthy volunteers. Several QSM pipelines were tested, including four distinct LBMs: sophisticated harmonic artifact reduction for phase data (SHARP), variable‐radius sophisticated harmonic artifact reduction for phase data (V‐SHARP), Laplacian boundary value background field removal (LBV), and one‐step total generalized variation (TGV). Results from distinct pipelines were compared using visual inspection, summary statistics of susceptibility in deep gray matter/white matter/venous regions of interest, and, in the healthy volunteers, regional susceptibility bias analysis and nonparametric tests. Results: Multi‐echo versus TE‐dependent QSM had higher regional accuracy, especially in high‐susceptibility regions and at shorter TEs. Everywhere except in the veins, a processing pipeline incorporating TGV provided the most temporally stable TE‐dependent QSM results with an accuracy similar to multi‐echo QSM. Conclusions: For TE‐dependent QSM, carefully choosing LBMs can minimize the introduction of LBM‐related nonlinear temporal susceptibility variations

Quantitative Susceptibility Mapping: Contrast Mechanisms and Clinical Applications.

Quantitative susceptibility mapping (QSM) is a recently developed MRI technique for quantifying the spatial distribution of magnetic susceptibility within biological tissues. It first uses the frequency shift in the MRI signal to map the magnetic field profile within the tissue. The resulting field map is then used to determine the spatial distribution of the underlying magnetic susceptibility by solving an inverse problem. The solution is achieved by deconvolving the field map with a dipole field, under the assumption that the magnetic field is a result of the superposition of the dipole fields generated by all voxels and that each voxel has its unique magnetic susceptibility. QSM provides improved contrast to noise ratio for certain tissues and structures compared to its magnitude counterpart. More importantly, magnetic susceptibility is a direct reflection of the molecular composition and cellular architecture of the tissue. Consequently, by quantifying magnetic susceptibility, QSM is becoming a quantitative imaging approach for characterizing normal and pathological tissue properties. This article reviews the mechanism generating susceptibility contrast within tissues and some associated applications

Fast quantitative susceptibility mapping with L1-regularization and automatic parameter selection

Purpose To enable fast reconstruction of quantitative susceptibility maps with total variation penalty and automatic regularization parameter selection. Methods ℓ[subscript 1]-Regularized susceptibility mapping is accelerated by variable splitting, which allows closed-form evaluation of each iteration of the algorithm by soft thresholding and fast Fourier transforms. This fast algorithm also renders automatic regularization parameter estimation practical. A weighting mask derived from the magnitude signal can be incorporated to allow edge-aware regularization. Results Compared with the nonlinear conjugate gradient (CG) solver, the proposed method is 20 times faster. A complete pipeline including Laplacian phase unwrapping, background phase removal with SHARP filtering, and ℓ[subscript 1]-regularized dipole inversion at 0.6 mm isotropic resolution is completed in 1.2 min using MATLAB on a standard workstation compared with 22 min using the CG solver. This fast reconstruction allows estimation of regularization parameters with the L-curve method in 13 min, which would have taken 4 h with the CG algorithm. The proposed method also permits magnitude-weighted regularization, which prevents smoothing across edges identified on the magnitude signal. This more complicated optimization problem is solved 5 times faster than the nonlinear CG approach. Utility of the proposed method is also demonstrated in functional blood oxygen level–dependent susceptibility mapping, where processing of the massive time series dataset would otherwise be prohibitive with the CG solver. Conclusion Online reconstruction of regularized susceptibility maps may become feasible with the proposed dipole inversion

Quantitative susceptibility mapping: Report from the 2016 reconstruction challenge

PURPOSE: The aim of the 2016 quantitative susceptibility mapping (QSM) reconstruction challenge was to test the ability of various QSM algorithms to recover the underlying susceptibility from phase data faithfully. METHODS: Gradient-echo images of a healthy volunteer acquired at 3T in a single orientation with 1.06 mm isotropic resolution. A reference susceptibility map was provided, which was computed using the susceptibility tensor imaging algorithm on data acquired at 12 head orientations. Susceptibility maps calculated from the single orientation data were compared against the reference susceptibility map. Deviations were quantified using the following metrics: root mean squared error (RMSE), structure similarity index (SSIM), high-frequency error norm (HFEN), and the error in selected white and gray matter regions. RESULTS: Twenty-seven submissions were evaluated. Most of the best scoring approaches estimated the spatial frequency content in the ill-conditioned domain of the dipole kernel using compressed sensing strategies. The top 10 maps in each category had similar error metrics but substantially different visual appearance. CONCLUSION: Because QSM algorithms were optimized to minimize error metrics, the resulting susceptibility maps suffered from over-smoothing and conspicuity loss in fine features such as vessels. As such, the challenge highlighted the need for better numerical image quality criteria

The effect of low resolution and coverage on the accuracy of susceptibility mapping

PURPOSE: Quantitative susceptibility mapping (QSM) has found increasing clinical applications. However, to reduce scan time, clinical acquisitions often use reduced resolution and coverage, particularly in the through-slice dimension. The effect of these factors on QSM has begun to be assessed using only balloon phantoms and downsampled brain images. Here, we investigate the effects (and their sources) of low resolution or coverage on QSM using both simulated and acquired images. METHODS: Brain images were acquired at 1 mm isotropic resolution and full brain coverage, and low resolution (up to 6 mm slice thickness) or coverage (down to 20 mm) in 5 healthy volunteers. Images at reduced resolution or coverage were also simulated in these volunteers and in a new, anthropomorphic, numerical phantom. Mean susceptibilities in 5 brain regions, including white matter, were investigated over varying resolution and coverage. RESULTS: The susceptibility map contrast decreased with increasing slice thickness and spacing, and with decreasing coverage below ~40 mm for 2 different QSM pipelines. Our simulations showed that calculated susceptibility values were erroneous at low resolution or very low coverage, because of insufficient sampling and overattenuation of the susceptibility-induced field perturbations. Susceptibility maps calculated from simulated and acquired images showed similar behavior. CONCLUSIONS: Both low resolution and low coverage lead to loss of contrast and errors in susceptibility maps. The widespread clinical practice of using low resolution and coverage does not provide accurate susceptibility maps. Simulations in images of healthy volunteers and in a new, anthropomorphic numerical phantom were able to accurately model low-resolution and low-coverage acquisitions

Investigating the effect of flow compensation and quantitative susceptibility mapping method on the accuracy of venous susceptibility measurement

Quantitative susceptibility mapping (QSM) is a promising non-invasive method for obtaining information relating to oxygen metabolism. However, the optimal acquisition sequence and QSM reconstruction method for reliable venous susceptibility measurements are unknown. Full flow compensation is generally recommended to correct for the influence of venous blood flow, although the effect of flow compensation on the accuracy of venous susceptibility values has not been systematically evaluated. In this study, we investigated the effect of different acquisition sequences, including different flow compensation schemes, and different QSM reconstruction methods on venous susceptibilities. Ten healthy subjects were scanned with five or six distinct QSM sequence designs using monopolar readout gradients and different flow compensation schemes. All data sets were processed using six different QSM pipelines and venous blood susceptibility was evaluated in whole-brain segmentations of the venous vasculature and single veins. The quality of vein segmentations and the accuracy of venous susceptibility values were analyzed and compared between all combinations of sequences and reconstruction methods. The influence of the QSM reconstruction method on average venous susceptibility values was found to be 2.7–11.6 times greater than the influence of the acquisition sequence, including flow compensation. The majority of the investigated QSM reconstruction methods tended to underestimate venous susceptibility values in the vein segmentations that were obtained. In summary, we found that multi-echo gradient-echo acquisition sequences without full flow compensation yielded venous susceptibility values comparable to sequences with full flow compensation. However, the QSM reconstruction method had a great influence on susceptibility values and thus needs to be selected carefully for accurate venous QSM