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

Recommended implementation of quantitative susceptibility mapping for clinical research in the brain: A consensus of the ISMRM electro-magnetic tissue properties study group

This article provides recommendations for implementing QSM for clinical brain research. It is a consensus of the International Society of Magnetic Resonance in Medicine, Electro-Magnetic Tissue Properties Study Group. While QSM technical development continues to advance rapidly, the current QSM methods have been demonstrated to be repeatable and reproducible for generating quantitative tissue magnetic susceptibility maps in the brain. However, the many QSM approaches available have generated a need in the neuroimaging community for guidelines on implementation. This article outlines considerations and implementation recommendations for QSM data acquisition, processing, analysis, and publication. We recommend that data be acquired using a monopolar 3D multi-echo gradient echo (GRE) sequence and that phase images be saved and exported in Digital Imaging and Communications in Medicine (DICOM) format and unwrapped using an exact unwrapping approach. Multi-echo images should be combined before background field removal, and a brain mask created using a brain extraction tool with the incorporation of phase-quality-based masking. Background fields within the brain mask should be removed using a technique based on SHARP or PDF, and the optimization approach to dipole inversion should be employed with a sparsity-based regularization. Susceptibility values should be measured relative to a specified reference, including the common reference region of the whole brain as a region of interest in the analysis. The minimum acquisition and processing details required when reporting QSM results are also provided. These recommendations should facilitate clinical QSM research and promote harmonized data acquisition, analysis, and reporting.

Super-resolution QSM in little or no additional time for imaging (NATIve) using 2D EPI imaging in 3 orthogonal planes

Quantitative Susceptibility Mapping has the potential to provide additional insights into neurological diseases but is typically based on a quite long (5-10 min) 3D gradient-echo scan which is highly sensitive to motion. We propose an ultra-fast acquisition based on three orthogonal (sagittal, coronal and axial) 2D simultaneous multi-slice EPI scans with 1 mm in-plane resolution and 3 mm thick slices. Images in each orientation are corrected for susceptibility-related distortions and co-registered with an iterative non-linear Minimum Deformation Averaging (Volgenmodel) approach to generate a high SNR, super-resolution data set with an isotropic resolution of close to 1 mm. The net acquisition time is 3 times the volume acquisition time of EPI or about 12 s, but the three volumes could also replace "dummy scans" in fMRI, making it feasible to acquire QSM in little or No Additional Time for Imaging (NATIve). NATIve QSM values agreed well with reference 3D GRE QSM in the basal ganglia in healthy subjects. In patients with multiple sclerosis, there was also a good agreement between the susceptibility values within lesions and control ROIs and all lesions which could be seen on 3D GRE QSMs could also be visualized on NATIve QSMs. The approach is faster than conventional 3D GRE by a factor of 25-50 and faster than 3D EPI by a factor of 3-5. As a 2D technique, NATIve QSM was shown to be much more robust to motion than the 3D GRE and 3D EPI, opening up the possibility of studying neurological diseases involving iron accumulation and demyelination in patients who find it difficult to lie still for long enough to acquire QSM data with conventional methods.

Modeling venous bias in resting state functional MRI metrics

Resting-state (rs) functional magnetic resonance imaging (fMRI) is used to detect low-frequency fluctuations in the blood oxygen-level dependent (BOLD) signal across brain regions. Correlations between temporal BOLD signal fluctuations are commonly used to infer functional connectivity. However, because BOLD is based on the dilution of deoxyhemoglobin, it is sensitive to veins of all sizes, and its amplitude is biased by draining veins. These biases affect local BOLD signal location and amplitude, and may also influence BOLD-derived connectivity measures, but the magnitude of this venous bias and its relation to vein size and proximity is unknown. Here, veins were identified using high-resolution quantitative susceptibility maps and utilized in a biophysical model to investigate systematic venous biases on common local rsfMRI-derived measures. Specifically, we studied the impact of vein diameter and distance to veins on the amplitude of low-frequency fluctuations (ALFF), fractional ALFF (fALFF), Hurst exponent (HE), regional homogeneity (ReHo), and eigenvector centrality values in the grey matter. Values were higher across all distances in smaller veins, and decreased with increasing vein diameter. Additionally, rsfMRI values associated with larger veins decrease with increasing distance from the veins. ALFF and ReHo were the most biased by veins, while HE and fALFF exhibited the smallest bias. Across all metrics, the amplitude of the bias was limited in voxel-wise data, confirming that venous structure is not the dominant source of contrast in these rsfMRI metrics. Finally, the models presented can be used to correct this venous bias in rsfMRI metrics.

Assessing robustness of quantitative susceptibility-based MRI radiomic features in patients with multiple sclerosis

Multiple Sclerosis (MS) is an autoimmune demyelinating disease characterised by changes in iron and myelin content. These biomarkers are detectable by Quantitative Susceptibility Mapping (QSM), an advanced Magnetic Resonance Imaging technique detecting magnetic properties. When analysed with radiomic techniques that exploit its intrinsic quantitative nature, QSM may furnish biomarkers to facilitate early diagnosis of MS and timely assessment of progression. In this work, we explore the robustness of QSM radiomic features by varying the number of grey levels (GLs) and echo times (TEs), in a sample of healthy controls and patients with MS. We analysed the white matter in total and within six clinically relevant tracts, including the cortico-spinal tract and the optic radiation. After optimising the number of GLs (n = 64), at least 65% of features were robust for each Volume of Interest (VOI), with no difference (p > .05) between left and right hemispheres. Different outcomes in feature robustness among the VOIs depend on their characteristics, such as volume and variance of susceptibility values. This study validated the processing pipeline for robustness analysis and established the reliability of QSM-based radiomics features against GLs and TEs. Our results provide important insights for future radiomics studies using QSM in clinical applications.

Investigating the effect of oblique image acquisition on the accuracy of QSM and a robust tilt correction method

PURPOSE

Quantitative susceptibility mapping (QSM) is used increasingly for clinical research where oblique image acquisition is commonplace, but its effects on QSM accuracy are not well understood.

THEORY AND METHODS

The QSM processing pipeline involves defining the unit magnetic dipole kernel, which requires knowledge of the direction of the main magnetic fieldB ^ 0 $$ {\hat{\boldsymbol{B}}}_{\mathbf{0}} $$ with respect to the acquired image volume axes. The direction ofB ^ 0 $$ {\hat{\boldsymbol{B}}}_{\mathbf{0}} $$ is dependent on the axis and angle of rotation in oblique acquisition. Using both a numerical brain phantom and in vivo acquisitions in 5 healthy volunteers, we analyzed the effects of oblique acquisition on magnetic susceptibility maps. We compared three tilt-correction schemes at each step in the QSM pipeline: phase unwrapping, background field removal and susceptibility calculation, using the RMS error and QSM-tuned structural similarity index.

RESULTS

Rotation of wrapped phase images gave severe artifacts. Background field removal with projection onto dipole fields gave the most accurate susceptibilities when the field map was first rotated into alignment withB ^ 0 $$ {\hat{\boldsymbol{B}}}_{\mathbf{0}} $$ . Laplacian boundary value and variable-kernel sophisticated harmonic artifact reduction for phase data background field removal methods gave accurate results without tilt correction. For susceptibility calculation, thresholded k-space division, iterative Tikhonov regularization, and weighted linear total variation regularization, all performed most accurately when local field maps were rotated into alignment withB ^ 0 $$ {\hat{\boldsymbol{B}}}_{\mathbf{0}} $$ before susceptibility calculation.

CONCLUSION

For accurate QSM, oblique acquisition must be taken into account. Rotation of images into alignment withB ^ 0 $$ {\hat{\boldsymbol{B}}}_{\mathbf{0}} $$ should be carried out after phase unwrapping and before background-field removal. We provide open-source tilt-correction code to incorporate easily into existing pipelines: https://github.com/o-snow/QSM_TiltCorrection.git.

Implications of quantitative susceptibility mapping at 7 Tesla MRI for microbleeds detection in cerebral small vessel disease

Background

Cerebral microbleeds (MBs) are a hallmark of cerebral small vessel disease (CSVD) and can be found on T2*-weighted sequences on MRI. Quantitative susceptibility mapping (QSM) is a postprocessing method that also enables MBs identification and furthermore allows to differentiate them from calcifications.

Aims

We explored the implications of using QSM at submillimeter resolution for MBs detection in CSVD.

Methods

Both 3 and 7 Tesla (T) MRI were performed in elderly participants without MBs and patients with CSVD. MBs were quantified on T2*-weighted imaging and QSM. Differences in the number of MBs were assessed, and subjects were classified in CSVD subgroups or controls both on 3T T2*-weighted imaging and 7T QSM.

Results

48 participants [mean age (SD) 70.9 (8.8) years, 48% females] were included: 31 were healthy controls, 6 probable cerebral amyloid angiopathy (CAA), 9 mixed CSVD, and 2 were hypertensive arteriopathy [HA] patients. After accounting for the higher number of MBs detected at 7T QSM (Median = Mdn; Mdn7T−QSM = 2.5; Mdn3T−T2 = 0; z = 4.90; p < 0.001) and false positive MBs (6.1% calcifications), most healthy controls (80.6%) demonstrated at least one MB and more MBs were discovered in the CSVD group.

Conclusions

Our observations suggest that QSM at submillimeter resolution improves the detection of MBs in the elderly human brain. A higher prevalence of MBs than so far known in healthy elderly was revealed.

Multi-echo dipole inversion for magnetic susceptibility mapping

PURPOSE

Reconstructing tissue magnetic susceptibility (QSM) from MRI phase data involves solving multiple consecutive ill-posed inverse problems such as phase unwrapping, background field removal, and field-to-source inversion. Multi-echo acquisitions present an additional challenge, as the magnetization field is typically computed from the multiple phase data prior to reconstructing the susceptibility map. Processing the multiple phase data introduces errors during the field estimation, violating assumptions of the subsequent inverse problems, manifesting as streaking artifacts in the susceptibility map. To address this challenge, we propose a multi-echo field-to-source forward model that forgoes the field estimation step. Moreover, we propose a fully general underestimation correction step to recover susceptibility sources that were regularized away during the field-to-source inversion.

METHODS

The multi-echo forward model and correction step were validated on the QSM Challenge 2.0 datasets and compared to the standard single field-to-source model in in vivo human brains using different types of deconvolution algorithms.

RESULTS

On the QSM Challenge 2.0 datasets the multi-echo forward model and correction step attain state-of-the-art results on all metrics by a wide margin. Experiments in in vivo brains show that the multi-echo model is in agreement with the single field-to-source model and that the proposed forward model and correction step can be used with any available dipole inversion method.

CONCLUSION

A multi-echo field-to-source forward model forgoes the need to fit multi-echo phase data and achieves state-of-the-art results on the QSM Challenge 2.0 data. Underestimated low-frequency susceptibility distributions can be partially recovered using a correction step.

Multi-echo quantitative susceptibility mapping: how to combine echoes for accuracy and precision at 3 Tesla

Purpose

To compare different multi‐echo combination methods for MRI QSM. Given the current lack of consensus, we aimed to elucidate how to optimally combine multi‐echo gradient‐recalled echo signal phase information, either before or after applying Laplacian‐base methods (LBMs) for phase unwrapping or background field removal.

Methods

Multi‐echo gradient‐recalled echo data were simulated in a numerical head phantom, and multi‐echo gradient‐recalled echo images were acquired at 3 Tesla in 10 healthy volunteers. To enable image‐based estimation of gradient‐recalled echo signal noise, 5 volunteers were scanned twice in the same session without repositioning. Five QSM processing pipelines were designed: 1 applied nonlinear phase fitting over TEs before LBMs; 2 applied LBMs to the TE‐dependent phase and then combined multiple TEs via either TE‐weighted or SNR‐weighted averaging; and 2 calculated TE‐dependent susceptibility maps via either multi‐step or single‐step QSM and then combined multiple TEs via magnitude‐weighted averaging. Results from different pipelines were compared using visual inspection; summary statistics of susceptibility in deep gray matter, white matter, and venous regions; phase noise maps (error propagation theory); and, in the healthy volunteers, regional fixed bias analysis (Bland–Altman) and regional differences between the means (nonparametric tests).

Results

Nonlinearly fitting the multi‐echo phase over TEs before applying LBMs provided the highest regional accuracy of χ and the lowest phase noise propagation compared to averaging the LBM‐processed TE‐dependent phase. This result was especially pertinent in high‐susceptibility venous regions.

Conclusion

For multi‐echo QSM, we recommend combining the signal phase by nonlinear fitting before applying LBMs.

Click here for author‐reader discussions

Diagnostic accuracy of quantitative susceptibility mapping in multiple system atrophy: The impact of echo time and the potential of histogram analysis

•

We performed histogram analysis on χ maps at different TEs on MSA patients and HC.

•

We found altered χ distribution in Pu, SN, GP, CN for MSAp and in SN, DN for MSAc.

•

QSM diagnostic accuracy is TE-dependent and is enhanced at short TEs.

•

Short TEs capture rapidly-decaying contributions of high χ sources.

•

Histogram features detect χ spatial heterogeneities improving diagnostic accuracy.

The non-invasive quantification of iron stores via Quantitative Susceptibility Mapping (QSM) could play an important role in the diagnosis and the differential diagnosis of atypical Parkinsonisms. However, the susceptibility (χ) values measured via QSM depend on echo time (TE). This effect relates to the microstructural organization within the voxel, whose composition can be altered by the disease. Moreover, pathological iron deposition in a brain area may not be spatially uniform, and conventional Region of Interest (ROI)-based analysis may fail in detecting alterations. Therefore, in this work we evaluated the impact of echo time on the diagnostic accuracy of QSM on a population of patients with Multiple System Atrophy (MSA) of either Parkinsonian (MSAp) or cerebellar (MSAc) phenotypes. In addition, we tested the potential of histogram analysis to improve QSM classification accuracy.

We enrolled 32 patients (19 MSAp and 13 MSAc) and 16 healthy controls, who underwent a 7T MRI session including a gradient-recalled multi-echo sequence for χ mapping. Nine histogram features were extracted from the χ maps computed for each TE in atlas-based ROIs covering deep brain nuclei, and compared among groups.

Alterations of susceptibility distribution were found in the Putamen, Substantia Nigra, Globus Pallidus and Caudate Nucleus for MSAp and in the Substantia Nigra and Dentate Nucleus for MSAc. Increased iron deposition was observed in a larger number of ROIs for the two shortest TEs and the standard deviation, the 75^th and the 90^th percentile were the most informative features yielding excellent diagnostic accuracy with area under the ROC curve > 0.9.

In conclusion, short TEs may enhance QSM diagnostic performances, as they can capture variations in rapidly-decaying contributions of high χ sources. The analysis of histogram features allowed to reveal fine heterogeneities in the spatial distribution of susceptibility alteration, otherwise undetected by a simple evaluation of ROI χ mean values.

Complementing canonical fMRI with functional Quantitative Susceptibility Mapping (fQSM) in modern neuroimaging research

Functional Quantitative Susceptibility Mapping (fQSM) allows for the quantitative measurement of time-varying magnetic susceptibility across cortical and subcortical brain structures with a potentially higher spatial specificity than conventional fMRI. While the usefulness of fQSM with General Linear Model and "On/Off" paradigms has been assessed, little is known about the potential applications and limitations of this technique in more sophisticated experimental paradigms and analyses, such as those currently used in modern neuroimaging. To thoroughly characterize fQSM activations, here we used 7T MRI, tonotopic mapping, as well as univariate (i.e., GLM and population Receptive Field) and multivariate (Representational Similarity Analysis; RSA) analyses. Although fQSM detected less tone-responsive voxels than fMRI, they were more consistently localized in gray matter. Also, the majority of active gray matter voxels exhibited negative fQSM response, signaling the expected oxyhemoglobin increase, whereas positive fQSM activations were mainly in white matter. Though fMRI- and fQSM-based tonotopic maps were overall comparable, the representation of frequency tunings in tone-sensitive regions was significantly more balanced for fQSM. Lastly, RSA revealed that frequency information from the auditory cortex could be successfully retrieved by using either methods. Overall, fQSM produces complementary results to conventional fMRI, as it captures small-scale variations in the activation pattern which inform multivariate measures. Although positive fQSM responses deserve further investigation, they do not impair the interpretation of contrasts of interest. The quantitative nature of fQSM, its spatial specificity and the possibility to simultaneously acquire canonical fMRI support the use of this technique for longitudinal and multicentric studies and pre-surgical mapping.

QSMxT: Robust masking and artifact reduction for quantitative susceptibility mapping

PURPOSE

Quantitative susceptibility mapping (QSM) estimates the spatial distribution of tissue magnetic susceptibilities from the phase of a gradient-echo signal. QSM algorithms require a signal mask to delineate regions with reliable phase for subsequent susceptibility estimation. Existing masking techniques used in QSM have limitations that introduce artifacts, exclude anatomical detail, and rely on parameter tuning and anatomical priors that narrow their application. Here, a robust masking and reconstruction procedure is presented to overcome these limitations and enable automated QSM processing. Moreover, this method is integrated within an open-source software framework: QSMxT.

METHODS

A robust masking technique that automatically separates reliable from less reliable phase regions was developed and combined with a two-pass reconstruction procedure that operates on the separated sources before combination, extracting more information and suppressing streaking artifacts.

RESULTS

Compared with standard masking and reconstruction procedures, the two-pass inversion reduces streaking artifacts caused by unreliable phase and high dynamic ranges of susceptibility sources. It is also robust across a range of acquisitions at 3 T in volunteers and phantoms, at 7 T in tumor patients, and in an in silico head phantom, with significant artifact and error reductions, greater anatomical detail, and minimal parameter tuning.

CONCLUSION

The two-pass masking and reconstruction procedure separates reliable from less reliable phase regions, enabling a more accurate QSM reconstruction that mitigates artifacts, operates without anatomical priors, and requires minimal parameter tuning. The technique and its integration within QSMxT makes QSM processing more accessible and robust to streaking artifacts.

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.

QSM reconstruction challenge 2.0: A realistic in silico head phantom for MRI data simulation and evaluation of susceptibility mapping procedures

PURPOSE

To create a realistic in silico head phantom for the second QSM reconstruction challenge and for future evaluations of processing algorithms for QSM.

METHODS

We created a digital whole-head tissue property phantom by segmenting and postprocessing high-resolution (0.64 mm isotropic), multiparametric MRI data acquired at 7 T from a healthy volunteer. We simulated the steady-state magnetization at 7 T using a Bloch simulator and mimicked a Cartesian sampling scheme through Fourier-based processing. Computer code for generating the phantom and performing the MR simulation was designed to facilitate flexible modifications of the phantom in the future, such as the inclusion of pathologies as well as the simulation of a wide range of acquisition protocols. Specifically, the following parameters and effects were implemented: TR and TE, voxel size, background fields, and RF phase biases. Diffusion-weighted imaging phantom data are provided, allowing future investigations of tissue-microstructure effects in phase and QSM algorithms.

RESULTS

The brain part of the phantom featured realistic morphology with spatial variations in relaxation and susceptibility values similar to the in vivo setting. We demonstrated some of the phantom's properties, including the possibility of generating phase data with nonlinear evolution over TE due to partial-volume effects or complex distributions of frequency shifts within the voxel.

CONCLUSION

The presented phantom and computer programs are publicly available and may serve as a ground truth in future assessments of the faithfulness of quantitative susceptibility reconstruction algorithms.

RESUMEN: A flexible class of multi-parameter qMRI protocols

PURPOSE

To introduce a class of fast 3D quantitative MRI (qMRI) schemes (RESUME^N, for N=1,…,4) that allow for a thorough characterization of microstructural properties of brain tissues.

METHODS

An arbitrary multi-echo GRE acquisition optimized for quantitative susceptibility mapping (QSM) is complemented with an appropriate low flip-angle GRE sequence drawn from four possible choices. The acquired signals are processed to analytically derive the longitudinal relaxation (R₁) and free induction decay (R₂^∗) rates, as well as the proton density (PD) and QSM. A comprehensive modeling of the excitation and B₁^- profiles and of the RF-spoiling is included in the acquisition and processing pipeline.

RESULTS

The RESUME^N maps appear homogeneous throughout the field-of-view and exhibit comparable values and high SNR across the considered range of N values.

CONCLUSIONS

The introduced schemes represent a class of robust and flexible strategies to derive a thorough and fast qMRI study, suitable for a whole-brain acquisition with isotropic voxel resolution of 700 μm in less than 15 min.

Quantitative MRI susceptibility mapping reveals cortical signatures of changes in iron, calcium and zinc in malformations of cortical development in children with drug-resistant epilepsy

OBJECTIVE

Malformations of cortical development (MCD), including focal cortical dysplasia (FCD), are the most common cause of drug-resistant focal epilepsy in children. Histopathological lesion characterisation demonstrates abnormal cell types and lamination, alterations in myelin (typically co-localised with iron), and sometimes calcification. Quantitative susceptibility mapping (QSM) is an emerging MRI technique that measures tissue magnetic susceptibility (χ) reflecting it's mineral composition. We used QSM to investigate abnormal tissue composition in a group of children with focal epilepsy with comparison to effective transverse relaxation rate (R2*) and Synchrotron radiation X-ray fluorescence (SRXRF) elemental maps. Our primary hypothesis was that reductions in χ would be found in FCD lesions, resulting from alterations in their iron and calcium content. We also evaluated deep grey matter nuclei for changes in χ with age.

METHODS

QSM and R2* maps were calculated for 40 paediatric patients with suspected MCD (18 histologically confirmed) and 17 age-matched controls. Patients' sub-groups were defined based on concordant electro-clinical or histopathology data. Quantitative investigation of QSM and R2* was performed within lesions, using a surface-based approach with comparison to homologous regions, and within deep brain regions using a voxel-based approach with regional values modelled with age and epilepsy as covariates. Synchrotron radiation X-ray fluorescence (SRXRF) was performed on brain tissue resected from 4 patients to map changes in iron, calcium and zinc and relate them to MRI parameters.

RESULTS

Compared to fluid-attenuated inversion recovery (FLAIR) or T1-weighted imaging, QSM improved lesion conspicuity in 5% of patients. In patients with well-localised lesions, quantitative profiling demonstrated decreased χ, but not R2*, across cortical depth with respect to the homologous regions. Contra-lateral homologous regions additionally exhibited increased χ at 2-3 mm cortical depth that was absent in lesions. The iron decrease measured by the SRXRF in FCDIIb lesions was in agreement with myelin reduction observed by Luxol Fast Blue histochemical staining. SRXRF analysis in two FCDIIb tissue samples showed increased zinc and calcium in one patient, and decreased iron in the brain region exhibiting low χ and high R2* in both patients. QSM revealed expected age-related changes in the striatum nuclei, substantia nigra, sub-thalamic and red nucleus.

CONCLUSION

QSM non-invasively revealed cortical/sub-cortical tissue alterations in MCD lesions and in particular that χ changes in FCDIIb lesions were consistent with reduced iron, co-localised with low myelin and increased calcium and zinc content. These findings suggest that measurements of cortical χ could be used to characterise tissue properties non-invasively in epilepsy lesions.