Quantitative susceptibility mapping (QSM): Decoding MRI data for a tissue magnetic biomarker

A subject-specific unsupervised deep learning method for quantitative susceptibility mapping using implicit neural representation

Quantitative susceptibility mapping (QSM) is an MRI-based technique that estimates the underlying tissue magnetic susceptibility based on phase signal. Deep learning (DL)-based methods have shown promise in handling the challenging ill-posed inverse problem for QSM reconstruction. However, they require extensive paired training data that are typically unavailable and suffer from generalization problems. Recent model-incorporated DL approaches also overlook the non-local effect of the tissue phase in applying the source-to-field forward model due to patch-based training constraint, resulting in a discrepancy between the prediction and measurement and subsequently suboptimal QSM reconstruction. This study proposes an unsupervised and subject-specific DL method for QSM reconstruction based on implicit neural representation (INR), referred to as INR-QSM. INR has emerged as a powerful framework for learning a high-quality continuous representation of the signal (image) by exploiting its internal information without training labels. In INR-QSM, the desired susceptibility map is represented as a continuous function of the spatial coordinates, parameterized by a fully-connected neural network. The weights are learned by minimizing a loss function that includes a data fidelity term incorporated by the physical model and regularization terms. Additionally, a novel phase compensation strategy is proposed for the first time to account for the non-local effect of tissue phase in data consistency calculation to make the physical model more accurate. Our experiments show that INR-QSM outperforms traditional established QSM reconstruction methods and the compared unsupervised DL method both qualitatively and quantitatively, and is competitive against supervised DL methods under data perturbations.

Quantitative MRI methods for the assessment of structure, composition, and function of musculoskeletal tissues in basic research and preclinical applications

Osteoarthritis (OA) is a disabling chronic disease involving the gradual degradation of joint structures causing pain and dysfunction. Magnetic resonance imaging (MRI) has been widely used as a non-invasive tool for assessing OA-related changes. While anatomical MRI is limited to the morphological assessment of the joint structures, quantitative MRI (qMRI) allows for the measurement of biophysical properties of the tissues at the molecular level. Quantitative MRI techniques have been employed to characterize tissues' structural integrity, biochemical content, and mechanical properties. Their applications extend to studying degenerative alterations, early OA detection, and evaluating therapeutic intervention. This article is a review of qMRI techniques for musculoskeletal tissue evaluation, with a particular emphasis on articular cartilage. The goal is to describe the underlying mechanism and primary limitations of the qMRI parameters, their association with the tissue physiological properties and their potential in detecting tissue degeneration leading to the development of OA with a primary focus on basic and preclinical research studies. Additionally, the review highlights some clinical applications of qMRI, discussing the role of texture-based radiomics and machine learning in advancing OA research.

SIPAS: A comprehensive susceptibility imaging process and analysis studio

Quantitative susceptibility mapping (QSM) is a rising MRI-based technology and quite a few QSM-related algorithms have been proposed to reconstruct maps of tissue susceptibility distribution from phase images. In this paper, we develop a comprehensive susceptibility imaging process and analysis studio (SIPAS) that can accomplish reliable QSM processing and offer a standardized evaluation system. Specifically, SIPAS integrates multiple methods for each step, enabling users to select algorithm combinations according to data conditions, and QSM maps could be evaluated by two aspects, including image quality indicators within all voxels and region-of-interest (ROI) analysis. Through a sophisticated design of user-friendly interfaces, the results of each procedure are able to be exhibited in axial, coronal, and sagittal views in real-time, meanwhile ROIs can be displayed in 3D rendering visualization. The accuracy and compatibility of SIPAS are demonstrated by experiments on multiple in vivo human brain datasets acquired from 3T, 5T, and 7T MRI scanners of different manufacturers. We also validate the QSM maps obtained by various algorithm combinations in SIPAS, among which the combination of iRSHARP and SFCR achieves the best results on its evaluation system. SIPAS is a comprehensive, sophisticated, and reliable toolkit that may prompt the QSM application in scientific research and clinical practice.

Unveiling metabolic patterns in dementia: Insights from high-resolution quantitative blood-oxygenation-level-dependent MRI

BACKGROUND

Oxygen extraction fraction (OEF) and deoxyhemoglobin (DoHb) levels reflect variations in cerebral oxygen metabolism in demented patients.

PURPOSE

Delineating the metabolic profiles evident throughout different phases of dementia necessitates an integrated analysis of OEF and DoHb levels. This is enabled by leveraging high-resolution quantitative blood oxygenation level dependent (qBOLD) analysis of magnitude images obtained from a multi-echo gradient-echo MRI (mGRE) scan performed on a 3.0 Tesla scanner.

METHODS

Achieving superior spatial resolution in qBOLD necessitates the utilization of an mGRE scan with only four echoes, which in turn limits the number of measurements compared to the parameters within the qBOLD model. Consequently, it becomes imperative to discard non-essential parameters to facilitate further analysis. This process entails transforming the qBOLD model into a format suitable for fitting the log-magnitude difference (L-MDif) profiles of the four echo magnitudes present in each brain voxel. In order to bolster spatial specificity, the log-difference qBOLD model undergoes refinement into a representative form, termed as r-qBOLD, particularly when applied to class-averaged L-MDif signals derived through k-means clustering of L-MDif signals from all brain voxels into a predetermined number of clusters. The agreement between parameters estimated using r-qBOLD for different cluster sizes is validated using Bland-Altman analysis, and the model's goodness-of-fit is evaluated using aχ 2 ${\chi ^2}$ -test. Retrospective MRI data of Alzheimer's disease (AD), mild cognitive impairment (MCI), and non-demented patients without neuropathological disorders, pacemakers, other implants, or psychiatric disorders, who completed a minimum of three visits prior to MRI enrolment, are utilized for the study.

RESULTS

Utilizing a cohort comprising 30 demented patients aged 65-83 years in stages 4-6 representing mild, moderate, and severe stages according to the clinical dementia rating (CDR), matched with an age-matched non-demented control group of 18 individuals, we conducted joint observations of OEF and DoHb levels estimated using r-qBOLD. The observations elucidate metabolic signatures in dementia based on OEF and DoHb levels in each voxel. Our principal findings highlight the significance of spatial patterns of metabolic profiles (metabolic patterns) within two distinct regimes: OEF levels exceeding the normal range (S1-regime), and OEF levels below the normal range (S2-regime). The S1-regime, accompanied by low DoHb levels, predominantly manifests in fronto-parietal and perivascular regions with increase in dementia severity. Conversely, the S2-regime, accompanied by low DoHb levels, is observed in medial temporal (MTL) regions. Other regions with abnormal metabolic patterns included the orbitofrontal cortex (OFC), medial-orbital prefrontal cortex (MOPFC), hypothalamus, ventro-medial prefrontal cortex (VMPFC), and retrosplenial cortex (RSP). Dysfunction in the OFC and MOPFC indicated cognitive and emotional impairment, while hypothalamic involvement potentially indicated preclinical dementia. Reduced metabolic activity in the RSP suggested early-stage AD related functional abnormalities.

CONCLUSIONS

Integrated analysis of OEF and DoHb levels using r-qBOLD reveals distinct metabolic signatures across dementia phases, highlighting regions susceptible to neuronal loss, vascular involvement, and preclinical indicators.

Quantitative susceptibility mapping in the brain reflects spatial expression of genes involved in iron homeostasis and myelination

Quantitative susceptibility mapping (QSM) is an MRI modality used to non-invasively measure iron content in the brain. Iron exhibits a specific anatomically varying pattern of accumulation in the brain across individuals. The highest regions of accumulation are the deep grey nuclei, where iron is stored in paramagnetic molecule ferritin. This form of iron is considered to be what largely contributes to the signal measured by QSM in the deep grey nuclei. It is also known that QSM is affected by diamagnetic myelin contents. Here, we investigate spatial gene expression of iron and myelin related genes, as measured by the Allen Human Brain Atlas, in relation to QSM images of age-matched subjects. We performed multiple linear regressions between gene expression and the average QSM signal within 34 distinct deep grey nuclei regions. Our results show a positive correlation (p < .05, corrected) between expression of ferritin and the QSM signal in deep grey nuclei regions. We repeated the analysis for other genes that encode proteins thought to be involved in the transport and storage of iron in the brain, as well as myelination. In addition to ferritin, our findings demonstrate a positive correlation (p < .05, corrected) between the expression of ferroportin, transferrin, divalent metal transporter 1, several gene markers of myelinating oligodendrocytes, and the QSM signal in deep grey nuclei regions. Our results suggest that the QSM signal reflects both the storage and active transport of iron in the deep grey nuclei regions of the brain.

Whole brain multiparametric mapping in two minutes using a dual-flip-angle stack-of-stars blipped multi-gradient-echo acquisition

A new MRI technique is presented for three-dimensional fast simultaneous whole brain mapping of myelin water fraction (MWF), T1, proton density (PD), R2*, magnetic susceptibility (QSM), and B1 transmit field (B1+). Phantom and human (N = 9) datasets were acquired using a dual-flip-angle blipped multi-gradient-echo (DFA-mGRE) sequence with a stack-of-stars (SOS) trajectory. Images were reconstructed using a subspace-based algorithm with a locally low-rank constraint. A novel joint-sparsity-constrained multicomponent T2*-T1 spectrum estimation (JMSE) algorithm is proposed to correct for the T1 saturation effect and B1+/B1- inhomogeneities in the quantification of MWF. A tissue-prior-based B1+ estimation algorithm was adapted for B1 correction in the mapping of T1 and PD. In the phantom study, measurements obtained at an acceleration factor (R) of 12 using prospectively under-sampled SOS showed good consistency (R2 > 0.997) with Cartesian reference for R2*/T1app/M0app. In the in vivo study, results of retrospectively under-sampled SOS with R = 6, 12, 18, showed good quality (structure similarity index measure > 0.95) compared with those of fully-sampled SOS. Besides, results of prospectively under-sampled SOS with R = 12 showed good consistency (intraclass correlation coefficient > 0.91) with Cartesian reference for T1/PD/B1+/MWF/QSM/R2*, and good reproducibility (coefficient of variation < 7.0 %) in the test-retest analysis for T1/PD/B1+/MWF/R2*. This study has demonstrated the feasibility of simultaneous whole brain multiparametric mapping with a two-minute scan using the DFA-mGRE SOS sequence, which may overcome a major obstacle for neurological applications of multiparametric MRI.

Multiparametric MRI dataset for susceptibility-based radiomic feature extraction and analysis

Multiple sclerosis (MS) is a progressive demyelinating disease impacting the central nervous system. Conventional Magnetic Resonance Imaging (MRI) techniques (e.g., T2w images) help diagnose MS, although they sometimes reveal non-specific lesions. Quantitative MRI techniques are capable of quantifying imaging biomarkers in vivo, offering the potential to identify specific signs related to pre-clinical inflammation. Among those techniques, Quantitative Susceptibility Mapping (QSM) is particularly useful for studying processes that influence the magnetic properties of brain tissue, such as alterations in myelin concentration. Because of its intrinsic quantitative nature, it is particularly well-suited to be analyzed through radiomics, including techniques that extract a high number of complex and multi-dimensional features from radiological images. The dataset presented in this work provides information about normal-appearing white matter (NAWM) in a cohort of MS patients and healthy controls. It includes QSM-based radiomic features from NAWM and its tracts, and MR sequences necessary to implement the pipeline: T1w, T2w, QSM, DWI. The workflow is outlined in this article, along with an application showing feature reliability assessment.

Improved test-retest reliability of R 2 * $$ {\mathrm{R}}_2^{\ast } $$ and susceptibility quantification using multishot multi-echo 3D EPI

PURPOSE

This study aimed to evaluate the potential of 3D EPI for improving the reliability ofT 2 * $$ {\mathrm{T}}_2^{\ast } $$ -weighted data and quantification ofR 2 * $$ {\mathrm{R}}_2^{\ast } $$ decay rate and susceptibility (χ) compared with conventional gradient-echo (GRE)-based acquisition.

METHODS

Eight healthy subjects in a wide age range were recruited. Each subject received repeated scans for both GRE and EPI acquisitions with an isotropic 1 mm resolution at 3 T. Maps ofR 2 * $$ {\mathrm{R}}_2^{\ast } $$ and χ were quantified, and their interscan differences were used to evaluate the test-retest reliability. Interprotocol differences ofR 2 * $$ {\mathrm{R}}_2^{\ast } $$ and χ between GRE and EPI were also measured voxel by voxel and in selected regions of interest to test the consistency between the two acquisition methods.

RESULTS

The quantifications ofR 2 * $$ {\mathrm{R}}_2^{\ast } $$ and χ using EPI protocols showed increased test-retest reliability with higher EPI factors up to 5 as performed in the experiment and were consistent with those based on GRE.

CONCLUSION

The result suggests that multishot multi-echo 3D EPI can be a useful alternative acquisition method forT 2 * $$ {\mathrm{T}}_2^{\ast } $$ -weighted MRI and quantification ofR 2 * $$ {\mathrm{R}}_2^{\ast } $$ and χ with reduced scan time, improved test-retest reliability, and similar accuracy compared with commonly used 3D GRE.

Phase-sensitive deep reconstruction method for rapid multiparametric MR fingerprinting and quantitative susceptibility mapping in the brain

INTRODUCTION

This study explores the potential of Magnetic Resonance Fingerprinting (MRF) with a novel Phase-Sensitivity Deep Reconstruction Network (PS-DRONE) for simultaneous quantification of T1, T2, Proton Density, B1+, phase and quantitative susceptibility mapping (QSM).

METHODS

Data were acquired at 3 T in vitro and in vivo using an optimized EPI-based MRF sequence. Phantom experiments were conducted using a standardized phantom for T1 and T2 maps and a custom-made agar-based gadolinium phantom for B1 and QSM maps. In vivo experiments included five healthy volunteers and one patient diagnosed with brain metastasis. PSDRONE maps were compared to reference maps obtained through standard imaging sequences.

RESULTS

Total scan time was 2 min for 32 slices and a resolution of [1 mm, 1 mm, 4.5 mm]. The reconstruction of T1, T2, Proton Density, B1+ and phase maps were reconstructed within 1 s. In the phantoms, PS-DRONE analysis presented accurate and strongly correlated T1 and T2 maps (r = 0.99) compared to the reference maps. B1 maps from PS-DRONE showed slightly higher values, though still correlated (r = 0.6) with the reference. QSM values showed a small bias but were strongly correlated (r = 0.99) with reference data. In the in vivo analysis, PS-DRONE-derived T1 and T2 values for gray and white matter matched reference values in healthy volunteers. PS-DRONE B1 and QSM maps showed strong correlations with reference values.

CONCLUSION

The PS-DRONE network enables concurrent acquisition of T1, T2, PD, B1+, phase and QSM maps, within 2 min of acquisition time and 1 s of reconstruction time.

Reduced oxygen extraction fraction in deep cerebral veins associated with cognitive impairment in multiple sclerosis

Studying the relationship between cerebral oxygen utilization and cognitive impairment is essential to understanding neuronal functional changes in the disease progression of multiple sclerosis (MS). This study explores the potential of using venous susceptibility in internal cerebral veins (ICVs) as an imaging biomarker for cognitive impairment in relapsing-remitting MS (RRMS) patients. Quantitative susceptibility mapping derived from fully flow-compensated MRI phase data was employed to directly measure venous blood oxygen saturation levels (SvO2) in the ICVs. Results revealed a significant reduction in the susceptibility of ICVs (212.4 ± 30.8 ppb vs 239.4 ± 25.9 ppb) and a significant increase of SvO2 (74.5 ± 1.89% vs 72.4 ± 2.23%) in patients with RRMS compared with age- and sex-matched healthy controls. Both the susceptibility of ICVs (r = 0.508, p = 0.031) and the SvO2 (r = -0.498, p = 0.036) exhibited a moderate correlation with cognitive decline in these patients assessed by the Paced Auditory Serial Addition Test, while no significant correlation was observed with clinical disability measured by the Expanded Disability Status Scale. The findings suggest that venous susceptibility in ICVs has the potential to serve as a specific indicator of oxygen metabolism and cognitive function in RRMS. .

Plug-and-Play latent feature editing for orientation-adaptive quantitative susceptibility mapping neural networks

Quantitative susceptibility mapping (QSM) is a post-processing technique for deriving tissue magnetic susceptibility distribution from MRI phase measurements. Deep learning (DL) algorithms hold great potential for solving the ill-posed QSM reconstruction problem. However, a significant challenge facing current DL-QSM approaches is their limited adaptability to magnetic dipole field orientation variations during training and testing. In this work, we propose a novel Orientation-Adaptive Latent Feature Editing (OA-LFE) module to learn the encoding of acquisition orientation vectors and seamlessly integrate them into the latent features of deep networks. Importantly, it can be directly Plug-and-Play (PnP) into various existing DL-QSM architectures, enabling reconstructions of QSM from arbitrary magnetic dipole orientations. Its effectiveness is demonstrated by combining the OA-LFE module into our previously proposed phase-to-susceptibility single-step instant QSM (iQSM) network, which was initially tailored for pure-axial acquisitions. The proposed OA-LFE-empowered iQSM, which we refer to as iQSM+, is trained in a simulated-supervised manner on a specially-designed simulation brain dataset. Comprehensive experiments are conducted on simulated and in vivo human brain datasets, encompassing subjects ranging from healthy individuals to those with pathological conditions. These experiments involve various MRI platforms (3T and 7T) and aim to compare our proposed iQSM+ against several established QSM reconstruction frameworks, including the original iQSM. The iQSM+ yields QSM images with significantly improved accuracies and mitigates artifacts, surpassing other state-of-the-art DL-QSM algorithms. The PnP OA-LFE module's versatility was further demonstrated by its successful application to xQSM, a distinct DL-QSM network for dipole inversion. In conclusion, this work introduces a new DL paradigm, allowing researchers to develop innovative QSM methods without requiring a complete overhaul of their existing architectures.

Quantitative susceptibility mapping through model-based deep image prior (MoDIP)

The data-driven approach of supervised learning methods has limited applicability in solving dipole inversion in Quantitative Susceptibility Mapping (QSM) with varying scan parameters across different objects. To address this generalization issue in supervised QSM methods, we propose a novel training-free model-based unsupervised method called MoDIP (Model-based Deep Image Prior). MoDIP comprises a small, untrained network and a Data Fidelity Optimization (DFO) module. The network converges to an interim state, acting as an implicit prior for image regularization, while the optimization process enforces the physical model of QSM dipole inversion. Experimental results demonstrate MoDIP's excellent generalizability in solving QSM dipole inversion across different scan parameters. It exhibits robustness against pathological brain QSM, achieving over 32 % accuracy improvement than supervised deep learning methods. It is also 33 % more computationally efficient and runs 4 times faster than conventional DIP-based approaches, enabling 3D high-resolution image reconstruction in under 4.5 min.

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.

DF-QSM: Data Fidelity based Hybrid Approach for Improved Quantitative Susceptibility Mapping of the Brain

Quantitative Susceptibility Mapping (QSM) is an advanced magnetic resonance imaging (MRI) technique to quantify the magnetic susceptibility of the tissue under investigation. Deep learning methods have shown promising results in deconvolving the susceptibility distribution from the measured local field obtained from the MR phase. Although existing deep learning based QSM methods can produce high-quality reconstruction, they are highly biased toward training data distribution with less scope for generalizability. This work proposes a hybrid two-step reconstruction approach to improve deep learning based QSM reconstruction. The susceptibility map prediction obtained from the deep learning methods has been refined in the framework developed in this work to ensure consistency with the measured local field. The developed method was validated on existing deep learning and model-based deep learning methods for susceptibility mapping of the brain. The developed method resulted in improved reconstruction for MRI volumes obtained with different acquisition settings, including deep learning models trained on constrained (limited) data settings.

So You Want to Image Myelin Using MRI: Magnetic Susceptibility Source Separation for Myelin Imaging

In MRI, researchers have long endeavored to effectively visualize myelin distribution in the brain, a pursuit with significant implications for both scientific research and clinical applications. Over time, various methods such as myelin water imaging, magnetization transfer imaging, and relaxometric imaging have been developed, each carrying distinct advantages and limitations. Recently, an innovative technique named as magnetic susceptibility source separation has emerged, introducing a novel surrogate biomarker for myelin in the form of a diamagnetic susceptibility map. This paper comprehensively reviews this cutting-edge method, providing the fundamental concepts of magnetic susceptibility, susceptibility imaging, and the validation of the diamagnetic susceptibility map as a myelin biomarker that indirectly measures myelin content. Additionally, the paper explores essential aspects of data acquisition and processing, offering practical insights for readers. A comparison with established myelin imaging methods is also presented, and both current and prospective clinical and scientific applications are discussed to provide a holistic understanding of the technique. This work aims to serve as a foundational resource for newcomers entering this dynamic and rapidly expanding field.

Distinguishing microgliosis and tau deposition in the mouse brain using paramagnetic and diamagnetic susceptibility source separation

Tauopathies, including Alzheimer's disease (AD), are neurodegenerative disorders characterized by hyperphosphorylated tau protein aggregates in the brain. In addition to protein aggregates, microglia-mediated inflammation and iron dyshomeostasis are other pathological features observed in AD and other tauopathies. It is known that these alterations at the subcellular level occur much before the onset of macroscopic tissue atrophy or cognitive deficits. The ability to detect these microstructural changes with MRI therefore has substantive importance for improved characterization of disease pathogenesis. In this study, we demonstrate that quantitative susceptibility mapping (QSM) with paramagnetic and diamagnetic susceptibility source separation has the potential to distinguish neuropathological alterations in a transgenic mouse model of tauopathy. 3D multi-echo gradient echo data were acquired from fixed brains of PS19 (Tau) transgenic mice and age-matched wild-type (WT) mice (n = 5 each) at 11.7 T. The multi-echo data were fit to a 3-pool complex signal model to derive maps of paramagnetic component susceptibility (PCS) and diamagnetic component susceptibility (DCS). Group-averaged signal fraction and composite susceptibility maps showed significant region-specific differences between the WT and Tau mouse brains. Significant bilateral increases in PCS and |DCS| were observed in specific hippocampal and cortical sub-regions of the Tau mice relative to WT controls. Comparison with immunohistological staining for microglia (Iba1) and phosphorylated-tau (AT8) further indicated that the PCS and DCS differences corresponded to regional microgliosis and tau deposition in the PS19 mouse brains, respectively. The results demonstrate that quantitative susceptibility source separation may provide sensitive imaging markers to detect distinct pathological alterations in tauopathies.

Evaluation of brain iron deposition in different cerebral arteries of acute ischaemic stroke patients using quantitative susceptibility mapping

AIM

To investigate differences in iron deposition between infarct and normal cerebral arterial regions in acute ischaemic stroke (AIS) patients using quantitative susceptibility mapping (QSM).

MATERIALS AND METHODS

Forty healthy controls and 40 AIS patients were recruited, and their QSM images were obtained. There were seven regions of interest (ROIs) in AIS patients, including the infarct regions of responsible arteries (R1), the non-infarct regions of responsible arteries (R2), the contralateral symmetrical sites of lesions (R3), and the non-responsible cerebral arterial regions (R4, R5, R6, R7). For the healthy controls, the cerebral arterial regions corresponding to the AIS patient group were selected as ROIs. The differences in corresponding ROI susceptibilities between AIS patients and healthy controls and the differences in susceptibilities between infarcted and non-infarct regions in AIS patients were compared.

RESULTS

The susceptibilities of infarct regions in AIS patients were significantly higher than those in healthy controls (p<0.0001). There was no significant difference in non-infarct regions between the two groups (p>0.05). The susceptibility of the infarct regions in AIS patients was significantly higher than those of the non-infarct region of responsible artery and non-responsible cerebral arterial regions (p<0.01).

CONCLUSIONS

Abnormal iron deposition detected by QSM in the infarct regions of AIS patients may not affect iron levels in the non-infarct regions of responsible arteries and normal cerebral arteries, which may open the door for potential new diagnostic and treatment strategies.

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.

Simultaneous time-of-flight MR angiography and quantitative susceptibility mapping with key time-of-flight features

A technique for combined time-of-flight (TOF) MR angiography (MRA) and quantitative susceptibility mapping (QSM) was developed with key features of standard three-dimensional (3D) TOF acquisitions, including multiple overlapping thin slab acquisition (MOTSA), ramped RF excitation, and venous saturation. The developed triple-echo 3D TOF-QSM sequence enabled TOF-MRA, susceptibility-weighted imaging (SWI), QSM, and R2* mapping. The effects of ramped RF, resolution, flip angle, venous saturation, and MOTSA were studied on QSM. Six volunteers were scanned at 3 T with the developed sequence, conventional TOF-MRA, and conventional SWI. Quantitative comparison of susceptibility values on QSM and normalized arterial and venous vessel-to-background contrasts on TOF and SWI were performed. The ramped RF excitation created an inherent phase variation in the raw phase. A generic correction factor was computed to remove the phase variation to obtain QSM without artifacts from the TOF-QSM sequence. No statistically significant difference was observed between the developed and standard QSM sequence for susceptibility values. However, maintaining standard TOF features led to compromises in signal-to-noise ratio for QSM and SWI, arising from the use of MOTSA rather than one large 3D slab, higher TOF spatial resolution, increased TOF background suppression due to larger flip angles, and reduced venous signal from venous saturation. In terms of vessel contrast, veins showed higher normalized contrast on SWI derived from TOF-QSM than the standard SWI sequence. While fast flowing arteries had reduced contrast compared with standard TOF-MRA, no statistical difference was observed for slow flowing arteries. Arterial contrast differences largely arise from the longer TR used in TOF-QSM over standard TOF-MRA to accommodate additional later echoes for SWI. In conclusion, although the sequence has a longer TR and slightly lower arterial contrast, provided an adequate correction is made for ramped RF excitation effects on phase, QSM may be performed from a multiecho sequence that includes all key TOF features, thus enabling simultaneous TOF-MRA, SWI, QSM, and R2* map computation.

Exploiting gradient-echo frequency evolution: Probing white matter microstructure and extracting bulk susceptibility-induced frequency for quantitative susceptibility mapping

PURPOSE

This work is to investigate the microstructure-induced frequency shift in white matter (WM) with crossing fibers and to separate the microstructure-related frequency shift from the bulk susceptibility-induced frequency shift by model fitting the gradient-echo (GRE) frequency evolution for potentially more accurate quantitative susceptibility mapping (QSM).

METHODS

A hollow-cylinder fiber model (HCFM) with two fiber populations was developed to investigate GRE frequency evolutions in WM voxels with microstructural orientation dispersion. The simulated and experimentally measured TE-dependent local frequency shift was then fitted to a simplified frequency evolution model to obtain a microstructure-related frequency difference parameter (∆ f $$ \Delta f $$ ) and a TE-independent bulk susceptibility-induced frequency shift (C f $$ {C}_f $$ ). The obtainedC f $$ {C}_f $$ was then used for QSM reconstruction. Reconstruction performances were evaluated using a numerical head phantom and in vivo data and then compared to other multi-echo combination methods.

RESULTS

GRE frequency evolutions and∆ f $$ \Delta f $$ -based tissue parameters in both parallel and crossing fibers determined from our simulations were comparable to those observed in vivo. The TE-dependent frequency fitting method outperformed other multi-echo combination methods in estimatingC f $$ {C}_f $$ in simulations. The fitted∆ f $$ \Delta f $$ ,C f $$ {C}_f $$ , and QSM could be improved further by navigator-based B0 fluctuation correction.

CONCLUSION

A HCFM with two fiber populations can be used to characterize microstructure-induced frequency shifts in WM regions with crossing fibers. HCFM-based TE-dependent frequency fitting provides tissue contrast related to microstructure (∆ f $$ \Delta f $$ ) and in addition may help improve the quantification accuracy ofC f $$ {C}_f $$ and the corresponding QSM.

Neuroimaging of Parkinson's disease by quantitative susceptibility mapping

Parkinson's disease (PD) is a common neurodegenerative disease, and apart from a few rare genetic causes, its pathogenesis remains largely unclear. Recent scientific interest has been captured by the involvement of iron biochemistry and the disruption of iron homeostasis, particularly within the brain regions specifically affected in PD. The advent of Quantitative Susceptibility Mapping (QSM) has enabled non-invasive quantification of brain iron in vivo by MRI, which has contributed to the understanding of iron-associated pathogenesis and has the potential for the development of iron-based biomarkers in PD. This review elucidates the biochemical underpinnings of brain iron accumulation, details advancements in iron-sensitive MRI technologies, and discusses the role of QSM as a biomarker of iron deposition in PD. Despite considerable progress, several challenges impede its clinical application after a decade of QSM studies. The initiation of multi-site research is warranted for developing robust, interpretable, and disease-specific biomarkers for monitoring PD disease progression.

Assessment for Carotid Atherosclerotic Plaque Using Vessel Wall Magnetic Resonance Imaging: A Multireader ROC Study to Determine Optimal Sequence for Detecting Vessel Wall Calcification

INTRODUCTION

We aimed to compare conventional vessel wall MR imaging techniques and quantitative susceptibility mapping (QSM) to determine the optimal sequence for detecting carotid artery calcification.

METHODS

Twenty-two patients who underwent carotid vessel wall MR imaging and neck CT were enrolled. Four slices of 6-mm sections from the bilateral internal carotid bifurcation were subdivided into 4 segments according to clock position (0-3, 3-6, 6-9, and 9-12) and assessed for calcification. Two blinded radiologists independently reviewed a total of 704 segments and scored the likelihood of calcification using a 5-point scale on spin-echo imaging, FLASH, and QSM. The observer performance for detecting calcification was evaluated by a multireader, multiple-case receiver operating characteristic study. Weighted κ statistics were calculated to assess interobserver agreement.

RESULTS

QSM had a mean area under the receiver operating characteristic curve of 0.85, which was significantly higher than that of any other sequence (p < 0.01) and showed substantial interreader agreement (κ = 0.68). A segment with a score of 3-5 was defined as positive, and a segment with a score of 1-2 was defined as negative; the sensitivity and specificity of QSM were 0.75 and 0.87, respectively.

CONCLUSION

QSM was the most reliable MR sequence for the detection of plaque calcification.

Quantitative susceptibility mapping for detection of kidney stones, hemorrhage differentiation, and cyst classification in ADPKD

BACKGROUND AND PURPOSE

The objective is to demonstrate feasibility of quantitative susceptibility mapping (QSM) in autosomal dominant polycystic kidney disease (ADPKD) patients and to compare imaging findings with traditional T1/T2w magnetic resonance imaging (MRI).

METHODS

Thirty-three consecutive patients (11 male, 22 female) diagnosed with ADPKD were initially selected. QSM images were reconstructed from the multiecho gradient echo data and compared to co-registered T2w, T1w, and CT images. Complex cysts were identified and classified into distinct subclasses based on their imaging features. Prevalence of each subclass was estimated.

RESULTS

QSM visualized two renal calcifications measuring 9 and 10 mm and three pelvic phleboliths measuring 2 mm but missed 24 calcifications measuring 1 mm or less and 1 larger calcification at the edge of the field of view. A total of 121 complex T1 hyperintense/T2 hypointense renal cysts were detected. 52 (43%) Cysts appeared hyperintense on QSM consistent with hemorrhage; 60 (49%) cysts were isointense with respect to simple cysts and normal kidney parenchyma, while the remaining 9 (7%) were hypointense. The presentation of the latter two complex cyst subtypes is likely indicative of proteinaceous composition without hemorrhage.

CONCLUSION

Our results indicate that QSM of ADPKD kidneys is possible and uniquely suited to detect large renal calculi without ionizing radiation and able to identify properties of complex cysts unattainable with traditional approaches.

Quantitative susceptibility mapping in multiple sclerosis: A systematic review and meta-analysis

BACKGROUND

Quantitative susceptibility mapping (QSM) is a quantitative measure based on magnetic resonance imaging sensitive to iron and myelin content. This makes QSM a promising non-invasive tool for multiple sclerosis (MS) in research and clinical practice.

OBJECTIVE

We performed a systematic review and meta-analysis on the use of QSM in MS.

METHODS

Our review was prospectively registered on PROSPERO (CRD42022309563). We searched five databases for studies published between inception and 30th April 2023. We identified 83 English peer-reviewed studies that applied QSM images on MS cohorts. Fifty-five included studies had at least one of the following outcome measures: deep grey matter QSM values in MS, either compared to healthy controls (HC) (k = 13) or correlated with the score on the Expanded Disability Status Scale (EDSS) (k = 7), QSM lesion characteristics (k = 22) and their clinical correlates (k = 17), longitudinal correlates (k = 11), histological correlates (k = 7), or correlates with other imaging techniques (k = 12). Two meta-analyses on deep grey matter (DGM) susceptibility data were performed, while the remaining findings could only be analyzed descriptively.

RESULTS

After outlier removal, meta-analyses demonstrated a significant increase in the basal ganglia susceptibility (QSM values) in MS compared to HC, caudate (k = 9, standardized mean difference (SDM) = 0.54, 95 % CI = 0.39-0.70, I2 = 46 %), putamen (k = 9, SDM = 0.38, 95 % CI = 0.19-0.57, I2 = 59 %), and globus pallidus (k = 9, SDM = 0.48, 95 % CI = 0.28-0.67, I2 = 60 %), whereas thalamic QSM values exhibited a significant reduction (k = 12, SDM = -0.39, 95 % CI = -0.66--0.12, I2 = 84 %); these susceptibility differences in MS were independent of age. Further, putamen QSM values positively correlated with EDSS (k = 4, r = 0.36, 95 % CI = 0.16-0.53, I2 = 0 %). Regarding rim lesions, four out of seven studies, representing 73 % of all patients, reported rim lesions to be associated with more severe disability. Moreover, lesion development from initial detection to the inactive stage is paralleled by increasing, plateauing (after about two years), and gradually decreasing QSM values, respectively. Only one longitudinal study provided clinical outcome measures and found no association. Histological data suggest iron content to be the primary source of QSM values in DGM and at the edges of rim lesions; further, when also considering data from myelin water imaging, the decrease of myelin is likely to drive the increase of QSM values within WM lesions.

CONCLUSIONS

We could provide meta-analytic evidence for DGM susceptibility changes in MS compared to HC; basal ganglia susceptibility is increased and, in the putamen, associated with disability, while thalamic susceptibility is decreased. Beyond these findings, further investigations are necessary to establish the role of QSM in MS for research or even clinical routine.

Quantifying cerebral microbleeds using quantitative susceptibility mapping from magnetization-prepared rapid gradient-echo

T1-weighted magnetization-prepared rapid gradient-echo (MPRAGE) is commonly included in brain studies for structural imaging using magnitude images; however, its phase images can provide an opportunity to assess microbleed burden using quantitative susceptibility mapping (QSM). This potential application for MPRAGE-based QSM was evaluated using in vivo and simulated measurements. Possible factors affecting image quality were also explored. Detection sensitivity was evaluated against standard multiecho gradient echo (MEGE) QSM using 3-T in vivo data of 15 subjects with a combined total of 108 confirmed microbleeds. The two methods were compared based on the microbleed size and susceptibility measurements. In addition, simulations explored the detection sensitivity of MPRAGE-QSM at different representative magnetic field strengths and echo times using microbleeds of different size, susceptibility, and location. Results showed that in vivo microbleeds appeared to be smaller (× 0.54) and of higher mean susceptibility (× 1.9) on MPRAGE-QSM than on MEGE-QSM, but total susceptibility estimates were in closer agreement (slope: 0.97, r2 : 0.94), and detection sensitivity was comparable. In simulations, QSM at 1.5 T had a low contrast-to-noise ratio that obscured the detection of many microbleeds. Signal-to-noise ratio (SNR) levels at 3 T and above resulted in better contrast and increased detection. The detection rates for microbleeds of minimum one-voxel diameter and 0.4-ppm susceptibility were 0.55, 0.80, and 0.88 at SNR levels of 1.5, 3, and 7 T, respectively. Size and total susceptibility estimates were more consistent than mean susceptibility estimates, which showed size-dependent underestimation. MPRAGE-QSM provides an opportunity to detect and quantify the size and susceptibility of microbleeds of at least one-voxel diameter at B0  of 3 T or higher with no additional time cost, when standard T2 *-weighted images are not available or have inadequate spatial resolution. The total susceptibility measure is more robust against sequence variations and might allow combining data from different protocols.

Comparison of three magnetic resonance imaging measures of brain iron in healthy and cocaine use disorder participants

Several magnetic resonance imaging (MRI) measures for quantifying endogenous nonheme brain iron have been proposed. These correspond to distinct physical properties with varying sensitivities and specificities to iron. Moreover, they may depend not only on tissue iron concentration, but also on the intravoxel spatial pattern of iron deposition, which is complex in many brain regions. Here, the three MRI brain iron measures of R 2 * , magnetic field correlation (MFC), and magnetic susceptibility are compared in several deep gray matter regions for both healthy participants (HPs) and individuals with cocaine use disorder (CUD). Their concordance is assessed from their correlations with each other and their relative dependencies on age. In addition, associations between the iron measures and microstructure in adjacent white matter regions are investigated by calculating their correlations with diffusion MRI measures from the internal capsule, and associations with cognition are determined by using results from a battery of standardized tests relevant to CUD. It is found that all three iron measures are strongly correlated with each other for the considered gray matter regions, but with correlation coefficients substantially less than one indicating important differences. The age dependencies of all three measures are qualitatively similar in most regions, except for the red nucleus, where the susceptibility has a significantly stronger correlation with age than R 2 * . Weak to moderate correlations are seen for the iron measures with several of the diffusion and cognitive measures, with the strongest correlations being obtained for R 2 * . The iron measures differ little between the HP and CUD groups, although susceptibility is significantly lower in the red nucleus for the CUD group. For the comparisons made, the iron measures behave similarly in most respects, but with notable quantitative differences. It is suggested that these differences may be, in part, attributable to a higher sensitivity to the spatial pattern of iron deposition for R 2 * and MFC than for susceptibility. This is supported most strongly by a sharp contrast between the values of the iron measures in the globus pallidus relative to those in the red nucleus. The observed correlations of the iron measures with diffusion and cognitive scores point to possible connections between gray matter iron, white matter microstructure, and cognition.

Harmonic field extension for QSM with reduced spatial coverage using physics-informed generative adversarial network

Quantitative susceptibility mapping (QSM) is frequently employed in investigating brain iron related to brain development and diseases within deep gray matter (DGM). Nonetheless, the acquisition of whole-brain QSM data is time-intensive. An alternative approach, focusing the QSM specifically on areas of interest such as the DGM by reducing the field-of-view (FOV), can significantly decrease scan times. However, severe susceptibility value underestimations have been reported during QSM reconstruction with a limited FOV, largely attributable to artifacts from incorrect background field removal in the boundary region. This presents a considerable barrier to the clinical use of QSM with small spatial coverages using conventional methods alone. To mitigate the propagation of these errors, we proposed a harmonic field extension method based on a physics-informed generative adversarial network. Both quantitative and qualitative results demonstrate that our method outperforms conventional methods and delivers results comparable to those obtained with full FOV. Furthermore, we demonstrate the versatility of our method by applying it to data acquired prospectively with limited FOV and to data from patients with Parkinson's disease. The method has shown significant improvements in local field results, with QSM outcomes. In a clear illustration of its feasibility and effectiveness in real clinical environments, our proposed method addresses the prevalent issue of susceptibility underestimation in QSM with small spatial coverage.

Iron Overload in Brain: Transport Mismatches, Microbleeding Events, and How Nanochelating Therapies May Counteract Their Effects

Iron overload in many brain regions is a common feature of aging and most neurodegenerative diseases. In this review, the causes, mechanisms, mathematical models, and possible therapies are summarized. Indeed, physiological and pathological conditions can be investigated using compartmental models mimicking iron trafficking across the blood-brain barrier and the Cerebrospinal Fluid-Brain exchange membranes located in the choroid plexus. In silico models can investigate the alteration of iron homeostasis and simulate iron concentration in the brain environment, as well as the effects of intracerebral iron chelation, determining potential doses and timing to recover the physiological state. Novel formulations of non-toxic nanovectors with chelating capacity are already tested in organotypic brain models and could be available to move from in silico to in vivo experiments.

A spatially adaptive regularization based three-dimensional reconstruction network for quantitative susceptibility mapping

Objective.Quantitative susceptibility mapping (QSM) is a new imaging technique for non-invasive characterization of the composition and microstructure ofin vivotissues, and it can be reconstructed from local field measurements by solving an ill-posed inverse problem. Even for deep learning networks, it is not an easy task to establish an accurate quantitative mapping between two physical quantities of different units, i.e. field shift in Hz and susceptibility value in ppm for QSM.Approach. In this paper, we propose a spatially adaptive regularization based three-dimensional reconstruction network SAQSM. A spatially adaptive module is specially designed and a set of them at different resolutions are inserted into the network decoder, playing a role of cross-modality based regularization constraint. Therefore, the exact information of both field and magnitude data is exploited to adjust the scale and shift of feature maps, and thus any information loss or deviation occurred in previous layers could be effectively corrected. The network encoding has a dynamic perceptual initialization, which enables the network to overcome receptive field intervals and also strengthens its ability to detect features of various sizes.Main results. Experimental results on the brain data of healthy volunteers, clinical hemorrhage and simulated phantom with calcification demonstrate that SAQSM can achieve more accurate reconstruction with less susceptibility artifacts, while perform well on the stability and generalization even for severe lesion areas.Significance. This proposed framework may provide a valuable paradigm to quantitative mapping or multimodal reconstruction.

Model resolution-based deconvolution for improved quantitative susceptibility mapping

Quantitative susceptibility mapping (QSM) utilizes the relationship between the measured local field and the unknown susceptibility map to perform dipole deconvolution. The aim of this work is to introduce and systematically evaluate the model resolution-based deconvolution for improved estimation of the susceptibility map obtained using the thresholded k-space division (TKD). A two-step approach has been proposed, wherein the first step involves the TKD susceptibility map computation and the second step involves the correction of this susceptibility map using the model-resolution matrix. The TKD-estimated susceptibility map can be expressed as the weighted average of the true susceptibility map, where the weights are determined by the rows of the model-resolution matrix, and hence a deconvolution of the TKD susceptibility map using the model-resolution matrix yields a better approximation to the true susceptibility map. The model resolution-based deconvolution is realized using closed-form, iterative, and sparsity-regularized implementations. The proposed approach was compared with L2 regularization, TKD, rescaled TKD in superfast dipole inversion, the modulated closed-form method, and iterative dipole inversion, as well as sparsity-regularized dipole inversion. It was observed that the proposed approach showed a substantial reduction in the streaking artifacts across 94 test volumes considered in this study. The proposed approach also showed better error reduction and edge preservation compared with other approaches. The proposed model resolution-based deconvolution compensates for the truncation of zero coefficients in the dipole kernel at the magic angle and hence provides a closer approximation to the true susceptibility map compared with other direct methods.

Direct comparison between 18F-Flortaucipir tau PET and quantitative susceptibility mapping in progressive supranuclear palsy

PURPOSE

The pattern of flortaucipir tau PET uptake is topographically similar to the pattern of magnetic susceptibility in progressive supranuclear palsy (PSP); both with increased signal in subcortical structures such as the basal ganglia and midbrain, suggesting that they may be closely related. However, their relationship remains unknown since no studies have directly compared these two modalities in the same PSP cohort. We hypothesized that some flortaucipir uptake in PSP is associated with magnetic susceptibility, and hence iron deposition. The aim of this study was to evaluate the regional relationship between flortaucipir uptake and magnetic susceptibility and to examine the effects of susceptibility on flortaucipir uptake in PSP.

METHODS

Fifty PSP patients and 67 cognitively normal controls were prospectively recruited and underwent three Tesla MRI and flortaucipir tau PET scans. Quantitative susceptibility maps were reconstructed from multi-echo gradient-echo MRI images. Region of interest (ROI) analysis was performed to obtain flortaucipir and susceptibility values in the subcortical regions. Relationships between flortaucipir and susceptibility signals were evaluated using partial correlation analysis in the subcortical ROIs and voxel-based analysis in the whole brain. The effects of susceptibility on flortaucipir uptake were examined by using the framework of mediation analysis.

RESULTS

Both flortaucipir and susceptibility were greater in PSP compared to controls in the putamen, pallidum, subthalamic nucleus, red nucleus, and cerebellar dentate (p<0.05). The ROI-based and voxel-based analyses showed that these two signals were positively correlated in these five regions (r = 0.36-0.59, p<0.05). Mediation analysis showed that greater flortaucipir uptake was partially explained by susceptibility in the putamen, pallidum, subthalamic nucleus, and red nucleus, and fully explained in the cerebellar dentate.

CONCLUSIONS

These results suggest that some of the flortaucipir uptake in subcortical regions in PSP is related to iron deposition. These findings will contribute to our understanding of the mechanisms underlying flortaucipir tau PET findings in PSP and other neurodegenerative diseases.

Reduced Oxygen Extraction Fraction in Deep Cerebral Veins Associated with Cognitive Impairment in Multiple Sclerosis

Studying the relationship between cerebral oxygen utilization and cognitive impairment is essential to understanding neuronal functional changes in the disease progression of multiple sclerosis (MS). This study explores the potential of using venous susceptibility in internal cerebral veins (ICVs) as an imaging biomarker for cognitive impairment in relapsing-remitting MS (RRMS) patients. Quantitative susceptibility mapping derived from fully flow-compensated MRI phase data was employed to directly measure venous blood oxygen saturation levels (SvO2) in the ICVs. Results revealed a significant reduction in the susceptibility of ICVs (212.4 ± 30.8 ppb vs 239.4 ± 25.9 ppb) and a significant increase of SvO2 (74.5 ± 1.89 % vs 72.4 ± 2.23 %) in patients with RRMS compared with age- and sex-matched healthy controls. Both the susceptibility of ICVs (r = 0.646, p = 0.004) and the SvO2 (r = -0.603, p = 0.008) exhibited a strong correlation with cognitive decline in these patients assessed by the Paced Auditory Serial Addition Test, while no significant correlation was observed with clinical disability measured by the Expanded Disability Status Scale. The findings suggest that venous susceptibility in ICVs has the potential to serve as a specific indicator of oxygen metabolism and cognitive function in RRMS.

mcLARO: Multi-contrast learned acquisition and reconstruction optimization for simultaneous quantitative multi-parametric mapping

PURPOSE

To develop a method for rapid sub-millimeter T1 , T2 ,T 2 * $$ {\mathrm{T}}_2^{\ast } $$ , and QSM mapping in a single scan using multi-contrast learned acquisition and reconstruction optimization (mcLARO).

METHODS

A pulse sequence was developed by interleaving inversion recovery and T2 magnetization preparations and single-echo and multi-echo gradient echo acquisitions, which sensitized k-space data to T1 , T2 ,T 2 * $$ {\mathrm{T}}_2^{\ast } $$ , and magnetic susceptibility. The proposed mcLARO optimized both the multi-contrast k-space under-sampling pattern and image reconstruction based on image feature fusion using a deep learning framework. The proposed mcLARO method withR = 8 $$ R=8 $$ under-sampling was validated in a retrospective ablation study and compared with other deep learning reconstruction methods, including MoDL and Wave-MoDL, using fully sampled data as reference. Various under-sampling ratios in mcLARO were investigated. mcLARO was also evaluated in a prospective study using separately acquired conventionally sampled quantitative maps as reference standard.

RESULTS

The retrospective ablation study showed improved image sharpness of mcLARO compared to the baseline network without the multi-contrast sampling pattern optimization or image feature fusion module. The under-sampling ratio experiment showed that higher under-sampling ratios resulted in blurrier images and lower precision of quantitative values. The prospective study showed that small or negligible bias and narrow 95% limits of agreement on regional T1 , T2 ,T 2 * $$ {\mathrm{T}}_2^{\ast } $$ , and QSM values by mcLARO (5:39 mins) compared to reference scans (40:03 mins in total). mcLARO outperformed MoDL and Wave-MoDL in terms of image sharpness, noise suppression, and artifact removal.

CONCLUSION

mcLARO enabled fast sub-millimeter T1 , T2 ,T 2 * $$ {\mathrm{T}}_2^{\ast } $$ , and QSM mapping in a single scan.

A latent code based multi-variable modulation network for susceptibility mapping

Quantitative susceptibility mapping (QSM) is a technique for obtaining quantitative information on tissue susceptibility and has shown promising potential for clinical applications, in which the magnetic susceptibility is calculated by solving an ill-posed inverse problem. Recently, deep learning-based methods are proposed to address this issue, but the diversity of data distribution was not well considered, and thus the model generalization is limited in clinical applications. In this paper, we propose a Latent Code based Multi-Variable modulation network for QSM reconstruction (LCMnet). Particularly, a specific modulation module is exploited to incorporate three variables, i.e., field map, magnitude image, and initial susceptibility. The latent code in the modulated convolution is learned from feature maps of the field data using the encoder-decoder framework. The susceptibility map pre-estimated from simple thresholding is the constant input of the module, thereby enhancing the network stability and accelerating training convergence. As another input, multi-level features generated by a cross-fusion block integrate the information of field and magnitude data effectively. Experimental results on in vivo human brain data, challenge data, clinical data and synthetic data demonstrate that the proposed method LCMnet can achieve outstanding performance on accurate susceptibility measurement and also excellent generalization.

A Causality-Driven Graph Convolutional Network for Postural Abnormality Diagnosis in Parkinsonians

Abnormal posture is a common movement disorder in the progress of Parkinson's disease (PD), and this abnormality can increase the risk of falls or even disabilities. The conventional assessment approach depends on the judgment of well-trained experts via canonical scales. However, this approach requires extensive clinical expertise and is highly subjective. Considering the potential of quantitative susceptibility mapping (QSM) in PD diagnosis, this study explored the QSM-based method for the automated classification between PD patients with and without postural abnormalities. Nevertheless, a major challenge is that unstable non-causal features typically lead to less reliable performance. Therefore, we propose a causality-driven graph-convolutional-network framework based on multi-instance learning, where performance stability is enhanced through the invariant prediction principle and causal interventions. Specifically, we adopt an intervention strategy that combines a non-causal intervenor with causal prediction. A stability constraint is proposed to ensure robust integrated prediction under different interventions. Moreover, an intra-class homogeneity constraint is enforced for each individually-learned causality scoring module to promote the extraction of group-level general features, and hence achieve a balance between subject-specific and group-level features. The proposed method demonstrated promising performance through extensive experiments on a real clinical dataset. Also, the features extracted by our method coincide with those reported in previous medical studies on PD posture abnormalities. In general, our work provides a clinically-valuable approach for automated, objective, and reliable diagnosis of postural abnormalities in Parkinsonians. Our source code is publicly available at https://github.com/SJTUBME-QianLab/CausalGCN-PDPA.

Renal MRI: From Nephron to NMR Signal

Renal diseases pose a significant socio-economic burden on healthcare systems. The development of better diagnostics and prognostics is well-recognized as a key strategy to resolve these challenges. Central to these developments are MRI biomarkers, due to their potential for monitoring of early pathophysiological changes, renal disease progression or treatment effects. The surge in renal MRI involves major cross-domain initiatives, large clinical studies, and educational programs. In parallel with these translational efforts, the need for greater (patho)physiological specificity remains, to enable engagement with clinical nephrologists and increase the associated health impact. The ISMRM 2022 Member Initiated Symposium (MIS) on renal MRI spotlighted this issue with the goal of inspiring more solutions from the ISMRM community. This work is a summary of the MIS presentations devoted to: 1) educating imaging scientists and clinicians on renal (patho)physiology and demands from clinical nephrologists, 2) elucidating the connection of MRI parameters with renal physiology, 3) presenting the current state of leading MR surrogates in assessing renal structure and functions as well as their next generation of innovation, and 4) describing the potential of these imaging markers for providing clinically meaningful renal characterization to guide or supplement clinical decision making. We hope to continue momentum of recent years and introduce new entrants to the development process, connecting (patho)physiology with (bio)physics, and conceiving new clinical applications. We envision this process to benefit from cross-disciplinary collaboration and analogous efforts in other body organs, but also to maximally leverage the unique opportunities of renal physiology. LEVEL OF EVIDENCE: 1 TECHNICAL EFFICACY STAGE: 2.

Histogram analysis of quantitative susceptibility mapping and apparent diffusion coefficient for identifying isocitrate dehydrogenase genotypes and tumor subtypes of adult-type diffuse gliomas

Background

Accurate preoperative identification of isocitrate dehydrogenase (IDH) genotypes and tumor subtypes is highly important for proper treatment planning and prognosis evaluation in patients with glioma. This study aimed to differentiate IDH genotypes and tumor subtypes of adult-type diffuse gliomas using histogram features of quantitative susceptibility mapping (QSM) and apparent diffusion coefficient (ADC).

Methods

This prospective study enrolled patients with suspected gliomas between March 2019 and January 2022 in a random series. Histogram features of QSM and ADC were extracted from the tumor parenchyma. The Mann-Whitney U test was used to compare the difference in histogram features between different IDH genotypes and among tumor subtypes. Receiver operating characteristic (ROC) curves were constructed to assess the corresponding diagnostic performance.

Results

This study included 47 patients with histopathologically confirmed adult-type diffuse gliomas. Totals of seven QSM features including 10th percentile (P10), 90th percentile (P90), interquartile range (IQR), maximum, mean absolute deviation (MAD), root mean squared (RMS), and variance, and five ADC features including P10, mean, median, RMS, and skewness exhibited significant differences between different IDH genotypes (P<0.05 for all), with the IQR of QSM demonstrating the highest area under curve (AUC) of 0.774 [95% confidence interval (CI): 0.635-0.913]. For separating tumor subtypes, the IQR of QSM also showed the highest AUC of 0.745 (95% CI: 0.566-0.924) for glioblastoma (GBM) versus astrocytoma and 0.848 (95% CI: 0.706-0.989) for GBM versus oligodendroglioma, but none of the features could discriminate astrocytoma from oligodendroglioma. The combination of the IQR of QSM, P10 of ADC, and age achieved the highest AUC of 0.910 (95% CI: 0.826-0.994) for IDH genotypes, and 0.939 (95% CI: 0.859-1.000) and 0.967 (95% CI: 0.904-1.000) for GBM versus astrocytoma and GBM versus oligodendroglioma, respectively.

Conclusions

QSM and ADC histogram features may serve as potential imaging markers for noninvasively assessing IDH genotypes and tumor subtypes of adult-type diffuse gliomas. Combining significant features may enhance the diagnostic performance substantially.

Quantitative susceptibility mapping shows lower brain iron content in children with childhood epilepsy with centrotemporal spikes

PURPOSE

The dysregulation of brain iron homeostasis is closely relevant to a multitude of chronic neurological disorders. This study employed quantitative susceptibility mapping (QSM) to detect and compare whole-brain iron content between childhood epilepsy with centrotemporal spikes (CECTS) children and typically developing children.

MATERIALS AND METHODS

32 children with CECTS and 25 age- and gender-matched healthy children were enrolled. All participants were imaged with 3.0-T MRI to acquire the structural and susceptibility-weighted data. The susceptibility-weighted data were processed using STISuite toolbox to obtain QSM. The magnetic susceptibility difference between the two groups was compared using voxel-wise and region of interest methods. Multivariable linear regression, controlling for age, were employed to investigate the associations between the brain magnetic susceptibility and age at onset.

RESULTS

Lower magnetic susceptibility was mainly observed in sensory- and motor-related brain regions in children with CECTS, including bilateral middle frontal gyrus, supplementary motor area, midcingulate cortex, paracentral lobule and precentral gyrus, the magnetic susceptibility of right paracentral lobule, right precuneus and left supplementary motor area were found to have positive correlation with the age at onset.

CONCLUSIONS

This study suggests that the potential iron deficiency in certain brain regions is associated with CECTS, which might be helpful for further illumination of potential pathogenesis mechanism of CECTS.

Magnetic Susceptibility in Progressive Supranuclear Palsy Variants, Parkinson's Disease, and Corticobasal Syndrome

BACKGROUND

Previous studies have shown that magnetic susceptibility is increased in several subcortical regions in progressive supranuclear palsy (PSP). However, it is still unclear how subcortical and cortical susceptibilities vary across different PSP variants, Parkinson's disease (PD), and corticobasal syndrome (CBS).

OBJECTIVE

This study aims to clarify the susceptibility profiles in the subcortical and cortical regions in different PSP variants, PD, and CBS.

METHODS

Sixty-four patients, 20 PSP-Richardson syndrome (PSP-RS), 9 PSP-parkinsonism (PSP-P), 7 PSP-progressive gait freezing, 4 PSP-postural instability, 11 PD, and 13 CBS, and 20 cognitively normal control subjects underwent a 3-Tesla magnetic resonance imaging scan to reconstruct quantitative susceptibility maps. Region-of-interest analysis was performed to obtain susceptibility in several subcortical and cortical regions. Bayesian linear mixed effect models were used to estimate susceptibility within group and differences between groups.

RESULTS

In the subcortical regions, patients with PSP-RS and PSP-P showed greater susceptibility than control subjects in the pallidum, substantia nigra, red nucleus, and cerebellar dentate (P < 0.05). Patients with PSP-RS also showed greater susceptibility than patients with PSP-progressive gait freezing, PD, and CBS in the red nucleus and cerebellar dentate, and patients with PSP-P showed greater susceptibility than PD in the red nucleus. Patients with PSP-postural instability and CBS showed greater susceptibility than control subjects in the pallidum and substantia nigra. No significant differences were observed in any cortical region.

CONCLUSIONS

The PSP variants and CBS had different patterns of magnetic susceptibility in the subcortical regions. The findings will contribute to our understanding about iron profiles and pathophysiology of PSP and may provide a potential biomarker to differentiate PSP variants, PD, and CBS. © 2023 International Parkinson and Movement Disorder Society.

Quantitative susceptibility mapping of brain iron in healthy aging and cognition

Quantitative susceptibility mapping (QSM) is a magnetic resonance imaging (MRI) technique that can assess the magnetic properties of cerebral iron in vivo. Although brain iron is necessary for basic neurobiological functions, excess iron content disrupts homeostasis, leads to oxidative stress, and ultimately contributes to neurodegenerative disease. However, some degree of elevated brain iron is present even among healthy older adults. To better understand the topographical pattern of iron accumulation and its relation to cognitive aging, we conducted an integrative review of 47 QSM studies of healthy aging, with a focus on five distinct themes. The first two themes focused on age-related increases in iron accumulation in deep gray matter nuclei versus the cortex. The overall level of iron is higher in deep gray matter nuclei than in cortical regions. Deep gray matter nuclei vary with regard to age-related effects, which are most prominent in the putamen, and age-related deposition of iron is also observed in frontal, temporal, and parietal cortical regions during healthy aging. The third theme focused on the behavioral relevance of iron content and indicated that higher iron in both deep gray matter and cortical regions was related to decline in fluid (speed-dependent) cognition. A handful of multimodal studies, reviewed in the fourth theme, suggest that iron interacts with imaging measures of brain function, white matter degradation, and the accumulation of neuropathologies. The final theme concerning modifiers of brain iron pointed to potential roles of cardiovascular, dietary, and genetic factors. Although QSM is a relatively recent tool for assessing cerebral iron accumulation, it has significant promise for contributing new insights into healthy neurocognitive aging.

A new 3D phase unwrapping method by region partitioning and local polynomial modeling in abdominal quantitative susceptibility mapping

Background

Accurate phase unwrapping is a critical prerequisite for successful applications in phase-related MRI, including quantitative susceptibility mapping (QSM) and susceptibility weighted imaging. However, many existing 3D phase unwrapping algorithms face challenges in the presence of severe noise, rapidly changing phase, and open-end cutline.

Methods

In this study, we introduce a novel 3D phase unwrapping approach utilizing region partitioning and a local polynomial model. Initially, the method leverages phase partitioning to create initial regions. Noisy voxels connecting areas within these regions are excluded and grouped into residual voxels. The connected regions within the region of interest are then reidentified and categorized into blocks and residual voxels based on voxel count thresholds. Subsequently, the method sequentially performs inter-block and residual voxel phase unwrapping using the local polynomial model. The proposed method was evaluated on simulation and in vivo abdominal QSM data, and was compared with the classical Region-growing, Laplacian_based, Graph-cut, and PRELUDE methods.

Results

Simulation experiments, conducted under different signal-to-noise ratios and phase change levels, consistently demonstrate that the proposed method achieves accurate unwrapping results, with mean error ratios not exceeding 0.01%. In contrast, the error ratios of Region-growing (N/A, 84.47%), Laplacian_based (20.65%, N/A), Graph-cut (2.26%, 20.71%), and PRELUDE (4.28%, 10.33%) methods are all substantially higher than those of the proposed method. In vivo abdominal QSM experiments further confirm the effectiveness of the proposed method in unwrapping phase data and successfully reconstructing susceptibility maps, even in scenarios with significant noise, rapidly changing phase, and open-end cutline in a large field of view.

Conclusion

The proposed method demonstrates robust and accurate phase unwrapping capabilities, positioning it as a promising option for abdominal QSM applications.

Straightforward Magnetic Resonance Temperature Measurements Combined with High Frame Rate and Magnetic Susceptibility Correction

Proton resonance frequency shift (PRFS) is an MRI-based simple temperature mapping method that exhibits higher spatial and temporal resolution than temperature mapping methods based on T1 relaxation time and diffusion. PRFS temperature measurements are validated against fiber-optic thermal sensors (FOSs). However, the use of FOSs may introduce temperature errors, leading to both underestimation and overestimation of PRFS measurements, primarily due to material susceptibility changes caused by the thermal sensors. In this study, we demonstrated susceptibility-corrected PRFS (scPRFS) with a high frame rate and accuracy for suitably distributed temperatures. A single-echo-based background removal technique was employed for phase variation correction, primarily owing to magnetic susceptibility, which enabled fast temperature mapping. The scPRFS was used to validate the temperature fidelity by comparing the temperatures of fiber-optic sensors and conventional PRFS through phantom-mimicked human and ex vivo experiments. This study demonstrates that scPRFS measurements in agar-gel are in good agreement with the thermal sensor readings, with a root mean square error (RMSE) of 0.33-0.36 °C in the phantom model and 0.12-0.16 °C in the ex vivo experiment. These results highlight the potential of scPRFS for precise thermal monitoring and ablation in both low- and high-temperature non-invasive therapies.

Free-breathing high isotropic resolution quantitative susceptibility mapping (QSM) of liver using 3D multi-echo UTE cones acquisition and respiratory motion-resolved image reconstruction

PURPOSE

To enable free-breathing and high isotropic resolution liver quantitative susceptibility mapping (QSM) using 3D multi-echo UTE cones acquisition and respiratory motion-resolved image reconstruction.

METHODS

Using 3D multi-echo UTE cones MRI, a respiratory motion was estimated from the k-space center of the imaging data. After sorting the k-space data with estimated motion, respiratory motion state-resolved reconstruction was performed for multi-echo data followed by nonlinear least-squares fitting for proton density fat fraction (PDFF),R 2 * $$ {\mathrm{R}}_2^{\ast } $$ , and fat-corrected B0 field maps. PDFF and B0 field maps were subsequently used for QSM reconstruction. The proposed method was compared with motion-averaged (gridding) reconstruction and conventional 3D multi-echo Cartesian MRI in moving gadolinium phantom and in vivo studies. Region of interest (ROI)-based linear regression analysis was performed on these methods to investigate correlations between gadolinium concentration and QSM in the phantom study and betweenR 2 * $$ {\mathrm{R}}_2^{\ast } $$ and QSM in in vivo study.

RESULTS

Cones with motion-resolved reconstruction showed sharper image quality compared to motion-averaged reconstruction with a substantial reduction of motion artifacts in both moving phantom and in vivo studies. For ROI-based linear regression analysis of the phantom study, susceptibility values from cones with motion-resolved reconstruction (QSM ppm $$ {\mathrm{QSM}}_{\mathrm{ppm}} $$  = 0.31 × gadolinium mM + $$ \times {\mathrm{gadolinium}}_{\mathrm{mM}}+ $$ 0.05,R 2 $$ {R}^2 $$  = 0.999) and Cartesian without motion (QSM ppm $$ {\mathrm{QSM}}_{\mathrm{ppm}} $$  = 0.32× gadolinium mM + $$ \times {\mathrm{gadolinium}}_{\mathrm{mM}}+ $$ 0.04,R 2 $$ {R}^2 $$  = 1.000) showed linear relationships with gadolinium concentrations and showed good agreement with each other. For in vivo, motion-resolved reconstruction showed higher goodness of fit (QSM ppm $$ {\mathrm{QSM}}_{\mathrm{ppm}} $$  = 0.00261 × R 2 s - 1 * - $$ \times {\mathrm{R}}_{2_{{\mathrm{s}}^{-1}}}^{\ast }- $$ 0.524,R 2 $$ {R}^2 $$  = 0.977) compared to motion-averaged reconstruction (QSM ppm $$ {\mathrm{QSM}}_{\mathrm{ppm}} $$ = 0.0021 × R 2 s - 1 * - $$ \times {\mathrm{R}}_{2_{{\mathrm{s}}^{-1}}}^{\ast }- $$ 0.572,R 2 $$ {R}^2 $$ = 0.723) in ROI-based linear regression analysis betweenR 2 * $$ {\mathrm{R}}_2^{\ast } $$ and QSM.

CONCLUSION

Feasibility of free-breathing liver QSM was demonstrated with motion-resolved 3D multi-echo UTE cones MRI, achieving high isotropic resolution currently unachievable in conventional Cartesian MRI.

Quantitative susceptibility mapping identifies hippocampal and other subcortical grey matter tissue composition changes in temporal lobe epilepsy

Temporal lobe epilepsy (TLE) is associated with widespread brain alterations. Using quantitative susceptibility mapping (QSM) alongside transverse relaxation rate (R 2 * ), we investigated regional brain susceptibility changes in 36 patients with left-sided (LTLE) or right-sided TLE (RTLE) secondary to hippocampal sclerosis, and 27 healthy controls (HC). We compared three susceptibility calculation methods to ensure image quality. Correlations of susceptibility andR 2 * with age of epilepsy onset, frequency of focal-to-bilateral tonic-clonic seizures (FBTCS), and neuropsychological test scores were examined. Weak-harmonic QSM (WH-QSM) successfully reduced noise and removed residual background field artefacts. Significant susceptibility increases were identified in the left putamen in the RTLE group compared to the LTLE group, the right putamen and right thalamus in the RTLE group compared to HC, and a significant susceptibility decrease in the left hippocampus in LTLE versus HC. LTLE patients who underwent epilepsy surgery showed significantly lower left-versus-right hippocampal susceptibility. SignificantR 2 * changes were found between TLE and HC groups in the amygdala, putamen, thalamus, and in the hippocampus. Specifically, decreased R2 * was found in the left and right hippocampus in LTLE and RTLE, respectively, compared to HC. Susceptibility andR 2 * were significantly correlated with cognitive test scores in the hippocampus, globus pallidus, and thalamus. FBTCS frequency correlated positively with ipsilateral thalamic and contralateral putamen susceptibility and withR 2 * in bilateral globi pallidi. Age of onset was correlated with susceptibility in the hippocampus and putamen, and withR 2 * in the caudate. Susceptibility andR 2 * changes observed in TLE groups suggest selective loss of low-myelinated neurons alongside iron redistribution in the hippocampi, predominantly ipsilaterally, indicating QSM's sensitivity to local pathology. Increased susceptibility andR 2 * in the thalamus and putamen suggest increased iron content and reflect disease severity.

Evaluating the role of Ginkgo biloba extract in the secondary prevention of acute ischemic stroke with cerebral microbleeds by quantitative susceptibility mapping (QSM)

BACKGROUND

EGb 761, a standardized dry extract of Ginkgo biloba leaves, has certain anti-inflammatory and thrombotic effects and can be used to treat cerebrovascular diseases.

METHODS

A total of 49 patients were randomly assigned to the Aspirin group (24 cases in Controlled group) and the Aspirin + Ginkgo biloba group (25 cases in Treatment group). The quantitative magnetic sensitivity and venous oxygen saturation of cerebral microbleeds were analyzed at admission, discharge, and after follow-up for 3 and 6 months.

RESULTS

The demographic details age, gender, and admission to NIHSS were not significantly different between the two groups (p < 0.05). Quantitative susceptibility mapping (QSM) showed that the magnetic sensitivity of patients in both groups remained stable after 3 and 6 months of follow-up, while the venous oxygen saturation of the Treatment group increased. The venous oxygen saturation at 3 and 6 months of follow-up was negatively correlated with the modified mRS grade score.

CONCLUSIONS

QSM can be used as a quantitative follow-up tool in monitoring both oxygen saturation and Magnetic susceptibility of microbleeds noninvasively in ischemic stroke patients with cerebral microbleeds. EGB combined with Aspirin can improve blood oxygen saturation in those patients and this effect is particularly significant in the long-term efficacy of secondary prevention.

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.

Robust quantitative susceptibility mapping via approximate message passing with parameter estimation

PURPOSE

For quantitative susceptibility mapping (QSM), the lack of ground-truth in clinical settings makes it challenging to determine suitable parameters for the dipole inversion. We propose a probabilistic Bayesian approach for QSM with built-in parameter estimation, and incorporate the nonlinear formulation of the dipole inversion to achieve a robust recovery of the susceptibility maps.

THEORY

From a Bayesian perspective, the image wavelet coefficients are approximately sparse and modeled by the Laplace distribution. The measurement noise is modeled by a Gaussian-mixture distribution with two components, where the second component is used to model the noise outliers. Through probabilistic inference, the susceptibility map and distribution parameters can be jointly recovered using approximate message passing (AMP).

METHODS

We compare our proposed AMP with built-in parameter estimation (AMP-PE) to the state-of-the-art L1-QSM, FANSI, and MEDI approaches on the simulated and in vivo datasets, and perform experiments to explore the optimal settings of AMP-PE. Reproducible code is available at: https://github.com/EmoryCN2L/QSM_AMP_PE.

RESULTS

On the simulated Sim2Snr1 dataset, AMP-PE achieved the lowest NRMSE, deviation from calcification moment and the highest SSIM, while MEDI achieved the lowest high-frequency error norm. On the in vivo datasets, AMP-PE is robust and successfully recovers the susceptibility maps using the estimated parameters, whereas L1-QSM, FANSI and MEDI typically require additional visual fine-tuning to select or double-check working parameters.

CONCLUSION

AMP-PE provides automatic and adaptive parameter estimation for QSM and avoids the subjectivity from the visual fine-tuning step, making it an excellent choice for the clinical setting.

Microstructure-Informed Myelin Mapping (MIMM) from Gradient Echo MRI using Stochastic Matching Pursuit

Quantification of the myelin content of the white matter is important for studying demyelination in neurodegenerative diseases such as Multiple Sclerosis (MS), particularly for longitudinal monitoring. A novel noninvasive MRI method, called Microstructure-Informed Myelin Mapping (MIMM), is developed to quantify the myelin volume fraction (MVF) by utilizing a multi gradient echo sequence (mGRE) and a detailed biophysical model of tissue microstructure. Myelin is modeled as anisotropic negative susceptibility source based on the Hollow Cylindrical Fiber Model (HCFM), and iron as isotropic positive susceptibility source in the extracellular region. Voxels with a range of biophysical parameters are simulated to create a dictionary of MR echo time magnitude signals and total susceptibility values. MRI signals measured using a mGRE sequence are then matched voxel-by-voxel to the created dictionary to obtain the spatial distributions of myelin and iron. Three different MIMM versions are presented to deal with the fiber orientation dependent susceptibility effects of the myelin sheaths: a basic variation, which assumes fiber orientation is an unknown to fit, two orientation informed variations, which assume the fiber orientation distribution is available either from a separate diffusion tensor imaging (DTI) acquisition or from a DTI atlas based fiber orientation map. While all showed a significant linear correlation with the reference method based on T2-relaxometry (p < 0.0001), DTI orientation informed and atlas orientation informed variations reduced overestimation at white matter tracts compared to the basic variation. Finally, the implications and usefulness of attaining an additional iron susceptibility distribution map are discussed.

Highlights

novel stochastic matching pursuit algorithm called microstructure-informed myelin mapping (MIMM) is developed to quantify Myelin Volume Fraction (MVF) using Magnetic Resonance Imaging (MRI) and microstructural modeling.utilizes a detailed biophysical model to capture the susceptibility effects on both magnitude and phase to quantify myelin and iron.matter fiber orientation effects are considered for the improved MVF quantification in the major fiber tracts.acquired myelin and iron maps may be utilized to monitor longitudinal disease progress.

Quantitative Susceptibility Mapping MRI in Deep-Brain Nuclei in First-Episode Psychosis

BACKGROUND

Psychosis is related to neurochemical changes in deep-brain nuclei, particularly suggesting dopamine dysfunctions. We used an magnetic resonance imaging-based technique called quantitative susceptibility mapping (QSM) to study these regions in psychosis. QSM quantifies magnetic susceptibility in the brain, which is associated with iron concentrations. Since iron is a cofactor in dopamine pathways and co-localizes with inhibitory neurons, differences in QSM could reflect changes in these processes.

METHODS

We scanned 83 patients with first-episode psychosis and 64 healthy subjects. We reassessed 22 patients and 21 control subjects after 3 months. Mean susceptibility was measured in 6 deep-brain nuclei. Using linear mixed models, we analyzed the effect of case-control differences, region, age, gender, volume, framewise displacement (FD), treatment duration, dose, laterality, session, and psychotic symptoms on QSM.

RESULTS

Patients showed a significant susceptibility reduction in the putamen and globus pallidus externa (GPe). Patients also showed a significant R2* reduction in GPe. Age, gender, FD, session, group, and region are significant predictor variables for QSM. Dose, treatment duration, and volume were not predictor variables of QSM.

CONCLUSIONS

Reduction in QSM and R2* suggests a decreased iron concentration in the GPe of patients. Susceptibility reduction in putamen cannot be associated with iron changes. Since changes observed in putamen and GPe were not associated with symptoms, dose, and treatment duration, we hypothesize that susceptibility may be a trait marker rather than a state marker, but this must be verified with long-term studies.