QSMxT: Robust masking and artifact reduction for quantitative susceptibility mapping

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.

Neurodesk: an accessible, flexible and portable data analysis environment for reproducible neuroimaging

Neuroimaging research requires purpose-built analysis software, which is challenging to install and may produce different results across computing environments. The community-oriented, open-source Neurodesk platform ( https://www.neurodesk.org/ ) harnesses a comprehensive and growing suite of neuroimaging software containers. Neurodesk includes a browser-accessible virtual desktop, command-line interface and computational notebook compatibility, allowing for accessible, flexible, portable and fully reproducible neuroimaging analysis on personal workstations, high-performance computers and the cloud.

Maximum spherical mean value filtering for whole-brain QSM

PURPOSE

To develop a tissue field-filtering algorithm, called maximum spherical mean value (mSMV), for reducing shadow artifacts in QSM of the brain without requiring brain-tissue erosion.

THEORY AND METHODS

Residual background field is a major source of shadow artifacts in QSM. The mSMV algorithm filters large field-magnitude values near the border, where the maximum value of the harmonic background field is located. The effectiveness of mSMV for artifact removal was evaluated by comparing existing QSM algorithms in numerical brain simulation as well as using in vivo human data acquired from 11 healthy volunteers and 93 patients.

RESULTS

Numerical simulation showed that mSMV reduces shadow artifacts and improves QSM accuracy. Better shadow reduction, as demonstrated by lower QSM variation in the gray matter and higher QSM image quality score, was also observed in healthy subjects and in patients with hemorrhages, stroke, and multiple sclerosis.

CONCLUSION

The mSMV algorithm allows QSM maps that are substantially equivalent to those obtained using SMV-filtered dipole inversion without eroding the volume of interest.

Improving quantitative susceptibility mapping for the identification of traumatic brain injury neurodegeneration at the individual level

BACKGROUND

Emerging evidence suggests that traumatic brain injury (TBI) is a major risk factor for developing neurodegenerative disease later in life. Quantitative susceptibility mapping (QSM) has been used by an increasing number of studies in investigations of pathophysiological changes in TBI. However, generating artefact-free quantitative susceptibility maps in brains with large focal lesions, as in the case of moderate-to-severe TBI (ms-TBI), is particularly challenging. To address this issue, we utilized a novel two-pass masking technique and reconstruction procedure (two-pass QSM) to generate quantitative susceptibility maps (QSMxT; Stewart et al., 2022, Magn Reson Med.) in combination with the recently developed virtual brain grafting (VBG) procedure for brain repair (Radwan et al., 2021, NeuroImage) to improve automated delineation of brain areas. We used QSMxT and VBG to generate personalised QSM profiles of individual patients with reference to a sample of healthy controls.

METHODS

Chronic ms-TBI patients (N = 8) and healthy controls (N = 12) underwent (multi-echo) GRE, and anatomical MRI (MPRAGE) on a 3T Siemens PRISMA scanner. We reconstructed the magnetic susceptibility maps using two-pass QSM from QSMxT. We then extracted values of magnetic susceptibility in grey matter (GM) regions (following brain repair via VBG) across the whole brain and determined if they deviate from a reference healthy control group [Z-score < -3.43 or > 3.43, relative to the control mean], with the aim of obtaining personalised QSM profiles.

RESULTS

Using two-pass QSM, we achieved susceptibility maps with a substantial increase in quality and reduction in artefacts irrespective of the presence of large focal lesions, compared to single-pass QSM. In addition, VBG minimised the loss of GM regions and exclusion of patients due to failures in the region delineation step. Our findings revealed deviations in magnetic susceptibility measures from the HC group that differed across individual TBI patients. These changes included both increases and decreases in magnetic susceptibility values in multiple GM regions across the brain.

CONCLUSIONS

We illustrate how to obtain magnetic susceptibility values at the individual level and to build personalised QSM profiles in ms-TBI patients. Our approach opens the door for QSM investigations in more severely injured patients. Such profiles are also critical to overcome the inherent heterogeneity of clinical populations, such as ms-TBI, and to characterize the underlying mechanisms of neurodegeneration at the individual level more precisely. Moreover, this new personalised QSM profiling could in the future assist clinicians in assessing recovery and formulating a neuroscience-guided integrative rehabilitation program tailored to individual TBI patients.

Improved dynamic distortion correction for fMRI using single-echo EPI and a readout-reversed first image (REFILL)

The boundaries between tissues with different magnetic susceptibilities generate inhomogeneities in the main magnetic field which change over time due to motion, respiration and system instabilities. The dynamically changing field can be measured from the phase of the fMRI data and corrected. However, methods for doing so need multi-echo data, time-consuming reference scans and/or involve error-prone processing steps, such as phase unwrapping, which are difficult to implement robustly on the MRI host. The improved dynamic distortion correction method we propose is based on the phase of the single-echo EPI data acquired for fMRI, phase offsets calculated from a triple-echo, bipolar reference scan of circa 3-10 s duration using a method which avoids the need for phase unwrapping and an additional correction derived from one EPI volume in which the readout direction is reversed. This Reverse-Encoded First Image and Low resoLution reference scan (REFILL) approach is shown to accurately measure B0 as it changes due to shim, motion and respiration, even with large dynamic changes to the field at 7 T, where it led to a > 20% increase in time-series signal to noise ratio compared to data corrected with the classic static approach. fMRI results from REFILL-corrected data were free of stimulus-correlated distortion artefacts seen when data were corrected with static field mapping. The method is insensitive to shim changes and eddy current differences between the reference scan and the fMRI time series, and employs calculation steps that are simple and robust, allowing most data processing to be performed in real time on the scanner image reconstruction computer. These improvements make it feasible to routinely perform dynamic distortion correction in fMRI.

Striatal Iron Deposition in Recreational MDMA (Ecstasy) Users

BACKGROUND

The common club drug MDMA (also known as ecstasy) enhances mood, sensory perception, energy, sociability, and euphoria. While MDMA has been shown to produce neurotoxicity in animal models, research on its potential neurotoxic effects in humans is inconclusive and has focused primarily on the serotonin system.

METHODS

We investigated 34 regular, largely pure MDMA users for signs of premature neurodegenerative processes in the form of increased iron load in comparison to a group of 36 age-, sex-, and education-matched MDMA-naïve control subjects. We used quantitative susceptibility mapping, a novel tool able to detect even small tissue (nonheme) iron accumulations. Cortical and relevant subcortical gray matter structures were grouped into 8 regions of interest and analyzed.

RESULTS

Significantly increased iron deposition in the striatum was evident in the MDMA user group. The effect survived correction for multiple comparisons and remained after controlling for relevant confounding factors, including age, smoking, and stimulant co-use. Although no significant linear relationship between measurements of the amounts of MDMA intake (hair analysis and self-reports) and quantitative susceptibility mapping values was observed, increased striatal iron deposition might nevertheless point to MDMA-induced neurotoxic processes. Additional factors (hyperthermia and simultaneous co-use of other substances) that possibly amplify neurotoxic effects of MDMA during the state of acute intoxication are discussed.

CONCLUSIONS

The demonstrated increased striatal iron accumulation may indicate that regular MDMA users have an increased risk potential for neurodegenerative diseases with progressing age.

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 quantitative susceptibility mapping (QSM) for clinical brain research. It is a consensus of the ISMRM 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 give rise to the need in the neuroimaging community for guidelines on implementation. This article describes relevant considerations and provides specific implementation recommendations for all steps in QSM data acquisition, processing, analysis, and presentation in scientific publications. We recommend that data be acquired using a monopolar 3D multi-echo GRE sequence, that phase images be saved and exported in DICOM format and unwrapped using an exact unwrapping approach. Multi-echo images should be combined before background removal, and a brain mask created using a brain extraction tool with the incorporation of phase-quality-based masking. Background fields should be removed within the brain mask 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 whole brain as a region of interest in the analysis, and QSM results should be reported with - as a minimum - the acquisition and processing specifications listed in the last section of the article. These recommendations should facilitate clinical QSM research and lead to increased harmonization in data acquisition, analysis, and reporting.

Neurodesk: An accessible, flexible, and portable data analysis environment for reproducible neuroimaging

Neuroimaging data analysis often requires purpose-built software, which can be challenging to install and may produce different results across computing environments. Beyond being a roadblock to neuroscientists, these issues of accessibility and portability can hamper the reproducibility of neuroimaging data analysis pipelines. Here, we introduce the Neurodesk platform, which harnesses software containers to support a comprehensive and growing suite of neuroimaging software (https://www.neurodesk.org/). Neurodesk includes a browser-accessible virtual desktop environment and a command line interface, mediating access to containerized neuroimaging software libraries on various computing platforms, including personal and high-performance computers, cloud computing and Jupyter Notebooks. This community-oriented, open-source platform enables a paradigm shift for neuroimaging data analysis, allowing for accessible, flexible, fully reproducible, and portable data analysis pipelines.

Multisite reproducibility of quantitative susceptibility mapping and effective transverse relaxation rate in deep gray matter at 3 T using locally optimized sequences in 24 traveling heads

Iron concentration in the human brain plays a crucial role in several neurodegenerative diseases and can be monitored noninvasively using quantitative susceptibility mapping (QSM) and effective transverse relaxation rate (R₂ *) mapping from multiecho T₂ *-weighted images. Large population studies enable better understanding of pathologies and can benefit from pooling multisite data. However, reproducibility may be compromised between sites and studies using different hardware and sequence protocols. This work investigates QSM and R₂ * reproducibility at 3 T using locally optimized sequences from three centers and two vendors, and investigates possible reduction of cross-site variability through postprocessing approaches. Twenty-four healthy subjects traveled between three sites and were scanned twice at each site. Scan-rescan measurements from seven deep gray matter regions were used for assessing within-site and cross-site reproducibility using intraclass correlation coefficient (ICC) and within-subject standard deviation (SDw) measures. In addition, multiple QSM and R₂ * postprocessing options were investigated with the aim to minimize cross-site sequence-related variations, including: mask generation approach, echo-timing selection, harmonizing spatial resolution, field map estimation, susceptibility inversion method, and linear field correction for magnitude images. The same-subject cross-site region of interest measurements for QSM and R₂ * were highly correlated (R²  ≥ 0.94) and reproducible (mean ICC of 0.89 and 0.82 for QSM and R₂ *, respectively). The mean cross-site SDw was 4.16 parts per billion (ppb) for QSM and 1.27 s^-1 for R₂ *. For within-site measurements of QSM and R₂ *, the mean ICC was 0.97 and 0.87 and mean SDw was 2.36 ppb and 0.97 s^-1 , respectively. The precision level is regionally dependent and is reduced in the frontal lobe, near brain edges, and in white matter regions. Cross-site QSM variability (mean SDw) was reduced up to 46% through postprocessing approaches, such as masking out less reliable regions, matching available echo timings and spatial resolution, avoiding the use of the nonconsistent magnitude contrast between scans in field estimation, and minimizing streaking artifacts.