Nonlinear dipole inversion (NDI) enables robust quantitative susceptibility mapping (QSM)

Comparison study of quantitative susceptibility mapping with GRAPPA and wave-CAIPI: reproducibility, consistency, and microbleeds detection

PURPOSE

We compared quantitative susceptibility mapping (QSM) with wave-CAIPI 9 × (QSM_WC9 ×) with reference standard QSM with GRAPPA 2 × (QSM_G2 ×) in two MR scanners. We also compared detectability of microbleeds in both QSMs to demonstrate clinical feasibility of both QSMs.

MATERIALS AND METHODS

This prospective study was approved by the institutional review board and written informed consent was obtained from each subject. Healthy subjects were recruited to evaluate intra-scanner reproducibility, inter-scanner consistency, and inter-sequence consistency of QSM_G2 × and QSM_WC9 × at 2 MR scanners. Susceptibility values measured with volume of interests (VOIs) were evaluated. Patients who were requested for susceptibility weighted imaging were also recruited in this study to measure microbleeds on QSM_G2 × and QSM_WC9 × . The number of microbleeds was compared between two QSMs.

RESULTS

Total 55 healthy subjects (male 34, female 21, 38.3 years [23-79]) were included in this study. We investigated reproducibility and consistency of QSM_WC9 × by comparing reference standard QSM_G2 × in two MR scanners in this study, and high correlation (ρ, 0.93-0.97) and high intraclass correlation coefficient (ICC) (0.97-0.99) were obtained. Sixty patients (male 30, female 30; age, 55.4 years [21-85]) were finally enrolled in this prospective study. The ICC of the detected number of microbleeds between QSM_G2 × and QSM_WC9 × was 0.99 (0.98-0.99).

CONCLUSION

QSM_WC9 × and reference standard QSM_G2 × in two MR scanners showed good reproducibility and consistency in estimating magnetic susceptibilities. QSM_WC9 × and QSM_G2 × were also comparable in terms of microbleeds detection with good agreement of raters and high ICC.

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.

ComBating inter-site differences in field strength: harmonizing preclinical traumatic brain injury MRI data

Integrating datasets from multiple sites and scanners can increase statistical power for neuroimaging studies but can also introduce significant inter-site confounds. We evaluated the effectiveness of ComBat, an empirical Bayes approach, to combine longitudinal preclinical MRI data acquired at 4.7 or 9.4 T at two different sites in Australia. Male Sprague Dawley rats underwent MRI on Days 2, 9, 28, and 150 following moderate/severe traumatic brain injury (TBI) or sham injury as part of Project 1 of the NIH/NINDS-funded Centre Without Walls EpiBioS4Rx project. Diffusion-weighted and multiple-gradient-echo images were acquired, and outcomes included QSM, FA, and ADC. Acute injury measures including apnea and self-righting reflex were consistent between sites. Mixed-effect analysis of ipsilateral and contralateral corpus callosum (CC) summary values revealed a significant effect of site on FA and ADC values, which was removed following ComBat harmonization. Bland-Altman plots for each metric showed reduced variability across sites following ComBat harmonization, including for QSM, despite appearing to be largely unaffected by inter-site differences and no effect of site observed. Following harmonization, the combined inter-site data revealed significant differences in the imaging metrics consistent with previously reported outcomes. TBI resulted in significantly reduced FA and increased susceptibility in the ipsilateral CC, and significantly reduced FA in the contralateral CC compared with sham-injured rats. Additionally, TBI rats also exhibited a reversal in ipsilateral CC ADC values over time with significantly reduced ADC at Day 9, followed by increased ADC 150 days after injury. Our findings demonstrate the need for harmonizing multi-site preclinical MRI data and show that this can be successfully achieved using ComBat while preserving phenotypical changes due to TBI.

Spinet-QSM: model-based deep learning with schatten p-norm regularization for improved quantitative susceptibility mapping

OBJECTIVE

Quantitative susceptibility mapping (QSM) provides an estimate of the magnetic susceptibility of tissue using magnetic resonance (MR) phase measurements. The tissue magnetic susceptibility (source) from the measured magnetic field distribution/local tissue field (effect) inherent in the MR phase images is estimated by numerically solving the inverse source-effect problem. This study aims to develop an effective model-based deep-learning framework to solve the inverse problem of QSM.

MATERIALS AND METHODS

This work proposes a Schatten p -norm-driven model-based deep learning framework for QSM with a learnable norm parameter p to adapt to the data. In contrast to other model-based architectures that enforce the l2 -norm or l1 -norm for the denoiser, the proposed approach can enforce any p -norm ( 0 < p ≤ 2 ) on a trainable regulariser.

RESULTS

The proposed method was compared with deep learning-based approaches, such as QSMnet, and model-based deep learning approaches, such as learned proximal convolutional neural network (LPCNN). Reconstructions performed using 77 imaging volumes with different acquisition protocols and clinical conditions, such as hemorrhage and multiple sclerosis, showed that the proposed approach outperformed existing state-of-the-art methods by a significant margin in terms of quantitative merits.

CONCLUSION

The proposed SpiNet-QSM showed a consistent improvement of at least 5% in terms of the high-frequency error norm (HFEN) and normalized root mean squared error (NRMSE) over other QSM reconstruction methods with limited training data.

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.

Incorporating a-priori information in deep learning models for quantitative susceptibility mapping via adaptive convolution

Quantitative susceptibility mapping (QSM) has attracted considerable interest for tissue characterization (e.g., iron and calcium accumulation, myelination, venous vasculature) in the human brain and relies on extensive data processing of gradient-echo MRI phase images. While deep learning-based field-to-susceptibility inversion has shown great potential, the acquisition parameters applied in clinical settings such as image resolution or image orientation with respect to the magnetic field have not been fully accounted for. Furthermore, the lack of comprehensive training data covering a wide range of acquisition parameters further limits the current QSM deep learning approaches. Here, we propose the integration of a priori information of imaging parameters into convolutional neural networks with our approach, adaptive convolution, that learns the mapping between the additional presented information (acquisition parameters) and the changes in the phase images associated with these varying acquisition parameters. By associating a-priori information with the network parameters itself, the optimal set of convolution weights is selected based on data-specific attributes, leading to generalizability towards changes in acquisition parameters. Moreover, we demonstrate the feasibility of pre-training on synthetic data and transfer learning to clinical brain data to achieve substantial improvements in the computation of susceptibility maps. The adaptive convolution 3D U-Net demonstrated generalizability in acquisition parameters on synthetic and in-vivo data and outperformed models lacking adaptive convolution or transfer learning. Further experiments demonstrate the impact of the side information on the adaptive model and assessed susceptibility map computation on simulated pathologic data sets and measured phase data.

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.

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.

Stuttering associated with a pathogenic variant in the chaperone protein cyclophilin 40

Stuttering is a common speech disorder that interrupts speech fluency and tends to cluster in families. Typically, stuttering is characterized by speech sounds, words or syllables which may be repeated or prolonged and speech that may be further interrupted by hesitations or 'blocks'. Rare variants in a small number of genes encoding lysosomal pathway proteins have been linked to stuttering. We studied a large four-generation family in which persistent stuttering was inherited in an autosomal dominant manner with disruption of the cortico-basal-ganglia-thalamo-cortical network found on imaging. Exome sequencing of three affected family members revealed the PPID c.808C>T (p.Pro270Ser) variant that segregated with stuttering in the family. We generated a Ppid p.Pro270Ser knock-in mouse model and performed ex vivo imaging to assess for brain changes. Diffusion-weighted MRI in the mouse revealed significant microstructural changes in the left corticospinal tract, as previously implicated in stuttering. Quantitative susceptibility mapping also detected changes in cortico-striatal-thalamo-cortical loop tissue composition, consistent with findings in affected family members. This is the first report to implicate a chaperone protein in the pathogenesis of stuttering. The humanized Ppid murine model recapitulates network findings observed in affected family members.

BFRnet: A deep learning-based MR background field removal method for QSM of the brain containing significant pathological susceptibility sources

INTRODUCTION

Background field removal (BFR) is a critical step required for successful quantitative susceptibility mapping (QSM). However, eliminating the background field in brains containing significant susceptibility sources, such as intracranial hemorrhages, is challenging due to the relatively large scale of the field induced by these pathological susceptibility sources.

METHOD

This study proposes a new deep learning-based method, BFRnet, to remove the background field in healthy and hemorrhagic subjects. The network is built with the dual-frequency octave convolutions on the U-net architecture, trained with synthetic field maps containing significant susceptibility sources. The BFRnet method is compared with three conventional BFR methods and one previous deep learning method using simulated and in vivo brains from 4 healthy and 2 hemorrhagic subjects. Robustness against acquisition field-of-view (FOV) orientation and brain masking are also investigated.

RESULTS

For both simulation and in vivo experiments, BFRnet led to the best visually appealing results in the local field and QSM results with the minimum contrast loss and the most accurate hemorrhage susceptibility measurements among all five methods. In addition, BFRnet produced the most consistent local field and susceptibility maps between different sizes of brain masks, while conventional methods depend drastically on precise brain extraction and further brain edge erosions. It is also observed that BFRnet performed the best among all BFR methods for acquisition FOVs oblique to the main magnetic field.

CONCLUSION

The proposed BFRnet improved the accuracy of local field reconstruction in the hemorrhagic subjects compared with conventional BFR algorithms. The BFRnet method was effective for acquisitions of tilted orientations and retained whole brains without edge erosion as often required by traditional BFR methods.

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.

DeepSTI: Towards tensor reconstruction using fewer orientations in susceptibility tensor imaging

Susceptibility tensor imaging (STI) is an emerging magnetic resonance imaging technique that characterizes the anisotropic tissue magnetic susceptibility with a second-order tensor model. STI has the potential to provide information for both the reconstruction of white matter fiber pathways and detection of myelin changes in the brain at mm resolution or less, which would be of great value for understanding brain structure and function in healthy and diseased brain. However, the application of STI in vivo has been hindered by its cumbersome and time-consuming acquisition requirement of measuring susceptibility induced MR phase changes at multiple head orientations. Usually, sampling at more than six orientations is required to obtain sufficient information for the ill-posed STI dipole inversion. This complexity is enhanced by the limitation in head rotation angles due to physical constraints of the head coil. As a result, STI has not yet been widely applied in human studies in vivo. In this work, we tackle these issues by proposing an image reconstruction algorithm for STI that leverages data-driven priors. Our method, called DeepSTI, learns the data prior implicitly via a deep neural network that approximates the proximal operator of a regularizer function for STI. The dipole inversion problem is then solved iteratively using the learned proximal network. Experimental results using both simulation and in vivo human data demonstrate great improvement over state-of-the-art algorithms in terms of the reconstructed tensor image, principal eigenvector maps and tractography results, while allowing for tensor reconstruction with MR phase measured at much less than six different orientations. Notably, promising reconstruction results are achieved by our method from only one orientation in human in vivo, and we demonstrate a potential application of this technique for estimating lesion susceptibility anisotropy in patients with multiple sclerosis.

Quantitative susceptibility mapping using multi-channel convolutional neural networks with dipole-adaptive multi-frequency inputs

Quantitative susceptibility mapping (QSM) quantifies the distribution of magnetic susceptibility and shows great potential in assessing tissue contents such as iron, myelin, and calcium in numerous brain diseases. The accuracy of QSM reconstruction was challenged by an ill-posed field-to-susceptibility inversion problem, which is related to the impaired information near the zero-frequency response of the dipole kernel. Recently, deep learning methods demonstrated great capability in improving the accuracy and efficiency of QSM reconstruction. However, the construction of neural networks in most deep learning-based QSM methods did not take the intrinsic nature of the dipole kernel into account. In this study, we propose a dipole kernel-adaptive multi-channel convolutional neural network (DIAM-CNN) method for the dipole inversion problem in QSM. DIAM-CNN first divided the original tissue field into high-fidelity and low-fidelity components by thresholding the dipole kernel in the frequency domain, and it then inputs the two components as additional channels into a multichannel 3D Unet. QSM maps from the calculation of susceptibility through multiple orientation sampling (COSMOS) were used as training labels and evaluation reference. DIAM-CNN was compared with two conventional model-based methods [morphology enabled dipole inversion (MEDI) and improved sparse linear equation and least squares (iLSQR) and one deep learning method (QSMnet)]. High-frequency error norm (HFEN), peak signal-to-noise-ratio (PSNR), normalized root mean squared error (NRMSE), and the structural similarity index (SSIM) were reported for quantitative comparisons. Experiments on healthy volunteers demonstrated that the DIAM-CNN results had superior image quality to those of the MEDI, iLSQR, or QSMnet results. Experiments on data with simulated hemorrhagic lesions demonstrated that DIAM-CNN produced fewer shadow artifacts around the bleeding lesion than the compared methods. This study demonstrates that the incorporation of dipole-related knowledge into the network construction has a potential to improve deep learning-based QSM reconstruction.

Toward a realistic in silico abdominal phantom for QSM

PURPOSE

QSM outside the brain has recently gained interest, particularly in the abdominal region. However, the absence of reliable ground truths makes difficult to assess reconstruction algorithms, whose quality is already compromised by additional signal contributions from fat, gases, and different kinds of motion. This work presents a realistic in silico phantom for the development, evaluation and comparison of abdominal QSM reconstruction algorithms.

METHODS

Synthetic susceptibility andR 2 * $$ {R}_2^{\ast } $$ maps were generated by segmenting and postprocessing the abdominal 3T MRI data from a healthy volunteer. Susceptibility andR 2 * $$ {R}_2^{\ast } $$ values in different tissues/organs were assigned according to literature and experimental values and were also provided with realistic textures. The signal was simulated using as input the synthetic QSM andR 2 * $$ {R}_2^{\ast } $$ maps and fat contributions. Three susceptibility scenarios and two acquisition protocols were simulated to compare different reconstruction algorithms.

RESULTS

QSM reconstructions show that the phantom allows to identify the main strengths and limitations of the acquisition approaches and reconstruction algorithms, such as in-phase acquisitions, water-fat separation methods, and QSM dipole inversion algorithms.

CONCLUSION

The phantom showed its potential as a ground truth to evaluate and compare reconstruction pipelines and algorithms. The publicly available source code, designed in a modular framework, allows users to easily modify the susceptibility,R 2 * $$ {R}_2^{\ast } $$ and TEs, and thus creates different abdominal scenarios.

Affine transformation edited and refined deep neural network for quantitative susceptibility mapping

Deep neural networks have demonstrated great potential in solving dipole inversion for Quantitative Susceptibility Mapping (QSM). However, the performances of most existing deep learning methods drastically degrade with mismatched sequence parameters such as acquisition orientation and spatial resolution. We propose an end-to-end AFfine Transformation Edited and Refined (AFTER) deep neural network for QSM, which is robust against arbitrary acquisition orientation and spatial resolution up to 0.6 mm isotropic at the finest. The AFTER-QSM neural network starts with a forward affine transformation layer, followed by a Unet for dipole inversion, then an inverse affine transformation layer, followed by a Residual Dense Network (RDN) for QSM refinement. Simulation and in-vivo experiments demonstrated that the proposed AFTER-QSM network architecture had excellent generalizability. It can successfully reconstruct susceptibility maps from highly oblique and anisotropic scans, leading to the best image quality assessments in simulation tests and suppressed streaking artifacts and noise levels for in-vivo experiments compared with other methods. Furthermore, ablation studies showed that the RDN refinement network significantly reduced image blurring and susceptibility underestimation due to affine transformations. In addition, the AFTER-QSM network substantially shortened the reconstruction time from minutes using conventional methods to only a few seconds.

Multi-echo dipole inversion for magnetic susceptibility mapping

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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

Incomplete spectrum QSM using support information

Introduction

Reconstructing a bounded object from incomplete k-space data is a well posed problem, and it was recently shown that this incomplete spectrum approach can be used to reconstruct undersampled MRI images with similar quality to compressed sensing approaches. Here, we apply this incomplete spectrum approach to the field-to-source inverse problem encountered in quantitative magnetic susceptibility mapping (QSM). The field-to-source problem is an ill-posed problem because of conical regions in frequency space where the dipole kernel is zero or very small, which leads to the kernel's inverse being ill-defined. These "ill-posed" regions typically lead to streaking artifacts in QSM reconstructions. In contrast to compressed sensing, our approach relies on knowledge of the image-space support, more commonly referred to as the mask, of our object as well as the region in k-space with ill-defined values. In the QSM case, this mask is usually available, as it is required for most QSM background field removal and reconstruction methods.

Methods

We tuned the incomplete spectrum method (mask and band-limit) for QSM on a simulated dataset from the most recent QSM challenge and validated the QSM reconstruction results on brain images acquired in five healthy volunteers, comparing incomplete spectrum QSM to current state-of-the art-methods: FANSI, nonlinear dipole inversion, and conventional thresholded k-space division.

Results

Without additional regularization, incomplete spectrum QSM performs slightly better than direct QSM reconstruction methods such as thresholded k-space division (PSNR of 39.9 vs. 39.4 of TKD on a simulated dataset) and provides susceptibility values in key iron-rich regions similar or slightly lower than state-of-the-art algorithms, but did not improve the PSNR in comparison to FANSI or nonlinear dipole inversion. With added (ℓ1-wavelet based) regularization the new approach produces results similar to compressed sensing based reconstructions (at sufficiently high levels of regularization).

Discussion

Incomplete spectrum QSM provides a new approach to handle the "ill-posed" regions in the frequency-space data input to QSM.

[A multi-channel input convolutional neural network for artifact reduction in quantitative susceptibility mapping]

OBJECTIVE

To develop a deep learning-based QSM reconstruction method for reducing artifacts to improve the accuracy of magnetic susceptibility results.

METHODS

To eliminate artifacts caused by susceptibility interfaces with gigantic differences, we propose a multi-channel input convolutional neural network for artifact reduction (MAR-CNN) for solving the dipole inversion problem in QSM. In this neural network, the original tissue field was first separated into two components, which were subsequently imported as additional channels into a multi-channel 3D U-Net. MAR-CNN was compared with 3 conventional model-based methods, namely truncated k-space deconvolution (TKD), morphology enabled dipole inversion (MEDI), and improved sparse linear equation and least squares method (iLSQR), and with a deep learning method (QSMnet). High-frequency error norm, peak signal-to-noise ratio, normalized root mean squared error, and structure similarity index were reported for quantitative comparisons.

RESULTS

Experiments on healthy volunteers demonstrated that the results obtained using MAR-CNN had superior peak signal-to-noise ratio (43.12±1.19) and normalized root mean squared error (51.98± 3.65) to those of TKD, MEDI, iLSQR and QSMnet. MAR-CNN outperformed QSMnet reconstruction on all the 4 quantitative metrics with significant differences (P < 0.05). Experiment on data of simulated hemorrhagic lesion demonstrated that MAR-CNN produced less shadow artifacts around the bleeding lesion than the other 4 methods.

CONCLUSION

The proposed MAR-CNN for artifact reduction is capable of improving the accuracy of deep learning- based QSM reconstruction to effectively reduce artifacts.

Sub-acute Changes on MRI Measures of Cerebral Blood Flow and Venous Oxygen Saturation in Concussed Australian Rules Footballers

Background

Sports-related concussion (SRC) is common in collision sport athletes. There is growing evidence that repetitive SRC can have serious neurological consequences, particularly when the repetitive injuries occur when the brain has yet to fully recover from the initial injury. Hence, there is a need to identify biomarkers that are capable of determining SRC recovery so that they can guide clinical decisions pertaining to return-to-play. Cerebral venous oxygen saturation (SvO₂) and cerebral blood flow (CBF) can be measured using magnetic resonance imaging (MRI) and may provide insights into changing energy demands and recovery following SRC.

Results

In this study we therefore investigated SvO₂ and CBF in a cohort of concussed amateur Australian Football athletes (i.e., Australia’s most participated collision sport). Male and female Australian footballers (n = 13) underwent MRI after being cleared to return to play following a mandatory 13-day recovery period and were compared to a group of control Australian footballers (n = 16) with no recent history of SRC (i.e., > 3 months since last SRC). Despite the concussed Australian footballers being cleared to return to play at the time of MRI, we found evidence of significantly increased susceptibility in the global white matter (p = 0.020) and a trend (F_5,21 = 2.404, p = 0.071) for reduced relative CBF (relCBF) compared to the control group. Further, there was evidence of an interaction between sex and injury in straight sinus susceptibility values (F_1,25 = 3.858, p = 0.061) which were decreased in female SRC athletes (p = 0.053). Of note, there were significant negative correlations between straight sinus susceptibility and relCBF suggesting impaired metabolic function after SRC.

Conclusions

These findings support the use of quantitative susceptibility mapping (QSM) and relCBF as sensitive indicators of SRC, and raise further concerns related to SRC guidelines that allow for return-to-play in less than two weeks.

Quantitative susceptibility mapping as an imaging biomarker for Alzheimer’s disease: The expectations and limitations

Alzheimer’s disease (AD) is the most common type of dementia and a distressing diagnosis for individuals and caregivers. Researchers and clinical trials have mainly focused on β-amyloid plaques, which are hypothesized to be one of the most important factors for neurodegeneration in AD. Meanwhile, recent clinicopathological and radiological studies have shown closer associations of tau pathology rather than β-amyloid pathology with the onset and progression of Alzheimer’s symptoms. Toward a biological definition of biomarker-based research framework for AD, the 2018 National Institute on Aging–Alzheimer’s Association working group has updated the ATN classification system for stratifying disease status in accordance with relevant pathological biomarker profiles, such as cerebral β-amyloid deposition, hyperphosphorylated tau, and neurodegeneration. In addition, altered iron metabolism has been considered to interact with abnormal proteins related to AD pathology thorough generating oxidative stress, as some prior histochemical and histopathological studies supported this iron-mediated pathomechanism. Quantitative susceptibility mapping (QSM) has recently become more popular as a non-invasive magnetic resonance technique to quantify local tissue susceptibility with high spatial resolution, which is sensitive to the presence of iron. The association of cerebral susceptibility values with other pathological biomarkers for AD has been investigated using various QSM techniques; however, direct evidence of these associations remains elusive. In this review, we first briefly describe the principles of QSM. Second, we focus on a large variety of QSM applications, ranging from common applications, such as cerebral iron deposition, to more recent applications, such as the assessment of impaired myelination, quantification of venous oxygen saturation, and measurement of blood– brain barrier function in clinical settings for AD. Third, we mention the relationships among QSM, established biomarkers, and cognitive performance in AD. Finally, we discuss the role of QSM as an imaging biomarker as well as the expectations and limitations of clinically useful diagnostic and therapeutic implications for AD.

Towards in vivo ground truth susceptibility for single-orientation deep learning QSM: a multi-orientation gradient-echo MRI dataset

Recently, deep neural networks have shown great potential for solving dipole inversion of quantitative susceptibility mapping (QSM) with improved results. However, these studies utilized their limited dataset for network training and inference, which may lead to untrustworthy conclusions. Thus, a common dataset is needed for a fair comparison between different QSM reconstruction networks. Additionally, finding an in vivo reference susceptibility map that matches acquired single-orientation phase data remains an open problem. Susceptibility tensor imaging (STI) χ₃₃ and Calculation of Susceptibility through Multiple Orientation Sampling (COSMOS) are considered reference susceptibility candidates. However, a large number of multi-orientation GRE data for both STI and COSMOS reconstruction are now unavailable for training supervised neural networks for QSM. In this study, we reported the largest multi-orientation dataset, to the best of our knowledge in the QSM research field, with a total of 144 scans from 8 healthy subjects collected using a 3D GRE sequence from the same MR scanner. In addition, the parcellation of deep gray matter is also provided for automatically extracting susceptibility values. Five recently developed deep neural networks, i.e., xQSM, QSMnet, autoQSM, LPCNN, and MoDL-QSM were performed on this dataset. This potential data source could provide a common framework and labels to test the accuracy and robustness of deep neural networks for QSM reconstruction. This dataset has the potential to provide a benchmark of reference susceptibility for the deep learning-based QSM methods. Additionally, the trained COSMOS-labeled and χ₃₃-labeled networks were tested on the pathological data to explore their potential applications. The data together with deep gray matter parcellation maps are now publicly available via an open repository at https://osf.io/yfms7/, and the raw multi-orientation GRE data were also available at https://osf.io/y6rc3/.

Instant tissue field and magnetic susceptibility mapping from MRI raw phase using Laplacian enhanced deep neural networks

Quantitative susceptibility mapping (QSM) is an MRI post-processing technique that produces spatially resolved magnetic susceptibility maps from phase data. However, the traditional QSM reconstruction pipeline involves multiple non-trivial steps, including phase unwrapping, background field removal, and dipole inversion. These intermediate steps not only increase the reconstruction time but accumulates errors. This study aims to overcome existing limitations by developing a Laplacian-of-Trigonometric-functions (LoT) enhanced deep neural network for near-instant quantitative field and susceptibility mapping (i.e., iQFM and iQSM) from raw MRI phase data. The proposed iQFM and iQSM methods were compared with established reconstruction pipelines on simulated and in vivo datasets. In addition, experiments on patients with intracranial hemorrhage and multiple sclerosis were also performed to test the generalization of the proposed neural networks. The proposed iQFM and iQSM methods in healthy subjects yielded comparable results to those involving the intermediate steps while dramatically improving reconstruction accuracies on intracranial hemorrhages with large susceptibilities. High susceptibility contrast between multiple sclerosis lesions and healthy tissue was also achieved using the proposed methods. Comparative studies indicated that the most significant contributor to iQFM and iQSM over conventional multi-step methods was the elimination of traditional Laplacian unwrapping. The reconstruction time on the order of minutes for traditional approaches was shortened to around 0.1 seconds using the trained iQFM and iQSM neural networks.

Blip up-down acquisition for spin- and gradient-echo imaging (BUDA-SAGE) with self-supervised denoising enables efficient T2 , T2 *, para- and dia-magnetic susceptibility mapping

PURPOSE

To rapidly obtain high resolution T₂ , T₂ *, and quantitative susceptibility mapping (QSM) source separation maps with whole-brain coverage and high geometric fidelity.

METHODS

We propose Blip Up-Down Acquisition for Spin And Gradient Echo imaging (BUDA-SAGE), an efficient EPI sequence for quantitative mapping. The acquisition includes multiple T₂ *-, T₂ '-, and T₂ -weighted contrasts. We alternate the phase-encoding polarities across the interleaved shots in this multi-shot navigator-free acquisition. A field map estimated from interim reconstructions was incorporated into the joint multi-shot EPI reconstruction with a structured low rank constraint to eliminate distortion. A self-supervised neural network (NN), MR-Self2Self (MR-S2S), was used to perform denoising to boost SNR. Using Slider encoding allowed us to reach 1 mm isotropic resolution by performing super-resolution reconstruction on volumes acquired with 2 mm slice thickness. Quantitative T₂ (=1/R₂ ) and T₂ * (=1/R₂ *) maps were obtained using Bloch dictionary matching on the reconstructed echoes. QSM was estimated using nonlinear dipole inversion on the gradient echoes. Starting from the estimated R₂ /R₂ * maps, R₂ ' information was derived and used in source separation QSM reconstruction, which provided additional para- and dia-magnetic susceptibility maps.

RESULTS

In vivo results demonstrate the ability of BUDA-SAGE to provide whole-brain, distortion-free, high-resolution, multi-contrast images and quantitative T₂ /T₂ * maps, as well as yielding para- and dia-magnetic susceptibility maps. Estimated quantitative maps showed comparable values to conventional mapping methods in phantom and in vivo measurements.

CONCLUSION

BUDA-SAGE acquisition with self-supervised denoising and Slider encoding enables rapid, distortion-free, whole-brain T₂ /T₂ * mapping at 1 mm isotropic resolution under 90 s.

Streaking artifact suppression of quantitative susceptibility mapping reconstructions via L1-norm data fidelity optimization (L1-QSM)

PURPOSE

The presence of dipole-inconsistent data due to substantial noise or artifacts causes streaking artifacts in quantitative susceptibility mapping (QSM) reconstructions. Often used Bayesian approaches rely on regularizers, which in turn yield reduced sharpness. To overcome this problem, we present a novel L1-norm data fidelity approach that is robust with respect to outliers, and therefore prevents streaking artifacts.

METHODS

QSM functionals are solved with linear and nonlinear L1-norm data fidelity terms using functional augmentation, and are compared with equivalent L2-norm methods. Algorithms were tested on synthetic data, with phase inconsistencies added to mimic lesions, QSM Challenge 2.0 data, and in vivo brain images with hemorrhages.

RESULTS

The nonlinear L1-norm-based approach achieved the best overall error metric scores and better streaking artifact suppression. Notably, L1-norm methods could reconstruct QSM images without using a brain mask, with similar regularization weights for different data fidelity weighting or masking setups.

CONCLUSION

The proposed L1-approach provides a robust method to prevent streaking artifacts generated by dipole-inconsistent data, renders brain mask calculation unessential, and opens novel challenging clinical applications such asassessing brain hemorrhages and cortical layers.

Clinical Application of MPRAGE Wave Controlled Aliasing in Parallel Imaging (Wave-CAIPI): A Comparative Study with MPRAGE GRAPPA

Purpose: To compare reliability and elucidate clinical application of magnetization-prepared rapid gradient-echo (MPRAGE) with 9-fold acceleration by using wave-controlled aliasing in parallel imaging (Wave-CAIPI 3 × 3) in comparison to conventional MPRAGE accelerated by using generalized autocalibrating partially parallel acquisition (GRAPPA) 2 × 1.

Methods: A total of 26 healthy volunteers and 33 patients were included in this study. Subjects were scanned with two MPRAGEs, GRAPPA 2 × 1 and Wave-CAIPI 3 × 3 acquired in 5 min 21 s and 1 min 42 s, respectively, on a 3T MR scanner. Healthy volunteers underwent additional two MPRAGEs (CAIPI 3 × 3 and GRAPPA 3 × 3). The image quality of the four MPRAGEs was visually evaluated with a 5-point scale in healthy volunteers, and the SNR of four MPRAGEs was also calculated by measuring the phantom 10 times with each MPRAGE. Based on the results of the visual evaluation, voxel-based morphometry (VBM) analyses, including subfield analysis, were performed only for GRAPPA 2 × 1 and Wave-CAIPI 3 × 3. Correlation of segmentation results between GRAPPA 2 × 1 and Wave-CAIPI 3 × 3 was assessed.

Results: In visual evaluations, scores for MPRAGE GRAPPA 2 × 1 (mean rank: 4.00) were significantly better than those for Wave-CAIPI 3 × 3 (mean rank: 3.00), CAIPI 3 × 3 (mean rank: 1.83), and GRAPPA 3 × 3 (mean rank: 1.17), and scores for Wave-CAIPI 3×3 were significantly better than those for CAIPI 3 × 3 and GRAPPA 3 × 3. Image noise was evident at the center for additional MPRAGE CAIPI 3 × 3 and GRAPPA 3 × 3. The correlation of segmentation results between GRAPPA 2 × 1 and Wave-CAIPI 3 × 3 was higher than 0.85 in all VOIs except globus pallidus. Subfield analysis of hippocampus also showed a high correlation between GRAPPA 2 × 1 and Wave-CAIPI 3 × 3.

Conclusion: MPRAGE Wave-CAIPI 3 × 3 shows relatively better contrast, despite of its short scan time of 1 min 42 s. The volumes derived from automated segmentation of MPRAGE Wave-CAIPI are considered to be reliable measures.

On the separation of susceptibility sources in quantitative susceptibility mapping: Theory and phantom validation with an in vivo application to multiple sclerosis lesions of different age

PURPOSE

In biological tissue, phase contrast is determined by multiple substances such as iron, myelin or calcifications. Often, these substances occur co-located within the same measurement volume. However, quantitative susceptibility mapping can solely measure the average susceptibility per voxel. To provide new insight in disease progression and mechanisms in neurological diseases, where multiple processes such as demyelination and iron accumulation occur simultaneously in the same location, a separation of susceptibility sources is desirable to disentangle the underlying susceptibility proportions.

METHODS

The basic concept of separating the susceptibility effects from sources with different sign within one voxel is to include information on relaxation rate ΔR₂^∗ in the quantitative susceptibility mapping reconstruction pipeline. The presented reconstruction algorithm is implemented as a constrained minimization problem and solved using conjugate gradients. The algorithm is evaluated using a software phantom and validated in MRI measurements with a phantom containing mixtures of microscopic positive and negative susceptibility sources. Data from three multiple sclerosis patients are used to show in vivo feasibility.

RESULTS

In numerical simulations, the feasibility of disentangling susceptibility sources within the same voxel was confirmed provided the critera of the static dephasing regime were fulfilled. In phantom experiments, the magnitude decay kernel, which is an essential reconstruction parameter of the algorithm, was determined to be D_m=194.5T^-1s^-1ppm^-1, and susceptibility sources could be separated in MRI measurement data.

CONCLUSIONS

In conclusion, in this study a detailed description of the implementation of an algorithm for the separation of positive and negative susceptibility sources within the same volume element as well as its limitations is presented and validated quantitatively in both simulation and phantom experiments for the first time. An application to multiple sclerosis lesions shows promising results for in vivo usability.

Deep grey matter quantitative susceptibility mapping from small spatial coverages using deep learning

INTRODUCTION

Quantitative Susceptibility Mapping (QSM) is generally acquired with full brain coverage, even though many QSM brain-iron studies focus on the deep grey matter (DGM) region only. Reducing the spatial coverage to the DGM vicinity can substantially shorten the scan time or enhance the spatial resolution without increasing scan time; however, this may lead to significant DGM susceptibility underestimation.

METHOD

A recently proposed deep learning-based QSM method, namely xQSM, is investigated to assess the accuracy of dipole inversion on reduced brain coverages. The xQSM method is compared with two conventional dipole inversion methods using simulated and in vivo experiments from 4 healthy subjects at 3T. Pre-processed magnetic field maps are extended symmetrically from the centre of globus pallidus in the coronal plane to simulate QSM acquisitions of difference spatial coverages, ranging from 100% (∼32mm) to 400% (∼128mm) of the actual DGM physical size.

RESULTS

The proposed xQSM network led to the lowest DGM contrast loss in both simulated and in vivo subjects, with the smallest susceptibility variation range across all spatial coverages. For the digital brain phantom simulation, xQSM improved the DGM susceptibility underestimation more than 20% in small spatial coverages, as compared to conventional methods. For the in vivo acquisition, less than 5% DGM susceptibility error was achieved in 48mm axial slabs using the xQSM network, while a minimum of 112mm coverage was required for conventional methods. It is also shown that the background field removal process performed worse in reduced brain coverages, which further deteriorated the subsequent dipole inversion.

CONCLUSION

The recently proposed deep learning-based xQSM method significantly improves the accuracy of DGM QSM from small spatial coverages as compared with conventional QSM algorithms, which can shorten DGM QSM acquisition time substantially.

MoDL-QSM: Model-based deep learning for quantitative susceptibility mapping

Quantitative susceptibility mapping (QSM) has demonstrated great potential in quantifying tissue susceptibility in various brain diseases. However, the intrinsic ill-posed inverse problem relating the tissue phase to the underlying susceptibility distribution affects the accuracy for quantifying tissue susceptibility. Recently, deep learning has shown promising results to improve accuracy by reducing the streaking artifacts. However, there exists a mismatch between the observed phase and the theoretical forward phase estimated by the susceptibility label. In this study, we proposed a model-based deep learning architecture that followed the STI (susceptibility tensor imaging) physical model, referred to as MoDL-QSM. Specifically, MoDL-QSM accounts for the relationship between STI-derived phase contrast induced by the susceptibility tensor terms (χ₁₃, χ₂₃ and χ₃₃) and the acquired single-orientation phase. The convolutional neural networks are embedded into the physical model to learn a regularization term containing prior information. χ₃₃ and phase induced by χ₁₃ and χ₂₃ terms were used as the labels for network training. Quantitative evaluation metrics were compared with recently developed deep learning QSM methods. The results showed that MoDL-QSM achieved superior performance, demonstrating its potential for future applications.

Single-step calculation of susceptibility through multiple orientation sampling

Quantitative susceptibility mapping (QSM) was developed to estimate the spatial distribution of magnetic susceptibility from MR signal phase acquired using a gradient echo (GRE) sequence. The field-to-susceptibility inversion in QSM is known to be ill-posed and needs numerical stabilization through either regularization or data oversampling. The calculation of susceptibility through the multiple orientation sampling (COSMOS) method uses phase data acquired at three or more head orientations to achieve a well-conditioned field-to-susceptibility inversion and is often considered the gold standard for in vivo QSM. However, the conventional COSMOS approach, here named multistep COSMOS (MSCOSMOS), solves the dipole inversion from the local field derived from raw GRE phase through multiple steps of phase preprocessing. Error propagations between these consecutive phase processing steps can thus affect the final susceptibility quantification. On the other hand, recently proposed single-step QSM (SSQSM) methods aim to solve an integrated inversion from unprocessed or total phase to mitigate such error propagations but have been limited to single orientation QSM. This study therefore aimed to test the feasibility of using single-step COSMOS (SSCOSMOS) to jointly perform background field removal and dipole inversion with multiple orientation sampling, which could serve as a better standard for gauging SSQSM methods. We incorporated multiple spherical mean value (SMV) kernels of various radii with the dipole inversion in SSCOSMOS. QSM reconstructions with SSCOSMOS and MSCOSMOS were compared using both simulations with a numerical head phantom and in vivo human brain data. SSCOSMOS permitted integrated background removal and dipole inversion without the need to adjust any regularization parameters. In addition, with sufficiently large SMV kernels, SSCOSMOS performed consistently better than MSCOSMOS in all the tested error metrics in our simulations, giving better susceptibility quantification and smaller reconstruction error. Consistent tissue susceptibility values were obtained between SSCOSMOS and MSCOSMOS.

xQSM: quantitative susceptibility mapping with octave convolutional and noise-regularized neural networks

Quantitative susceptibility mapping (QSM) provides a valuable MRI contrast mechanism that has demonstrated broad clinical applications. However, the image reconstruction of QSM is challenging due to its ill-posed dipole inversion process. In this study, a new deep learning method for QSM reconstruction, namely xQSM, was designed by introducing noise regularization and modified octave convolutional layers into a U-net backbone and trained with synthetic and in vivo datasets, respectively. The xQSM method was compared with two recent deep learning (QSMnet⁺ and DeepQSM) and two conventional dipole inversion (MEDI and iLSQR) methods, using both digital simulations and in vivo experiments. Reconstruction error metrics, including peak signal-to-noise ratio, structural similarity, normalized root mean squared error and deep gray matter susceptibility measurements, were evaluated for comparison of the different methods. The results showed that the proposed xQSM network trained with in vivo datasets achieved the best reconstructions of all the deep learning methods. In particular, it led to, on average, 32.3%, 25.4% and 11.7% improvement in the accuracy of globus pallidus susceptibility estimation for digital simulations and 39.3%, 21.8% and 6.3% improvements for in vivo acquisitions compared with DeepQSM, QSMnet⁺ and iLSQR, respectively. It also exhibited the highest linearity against different susceptibility intensity scales and demonstrated the most robust generalization capability to various spatial resolutions of all the deep learning methods. In addition, the xQSM method also substantially shortened the reconstruction time from minutes using MEDI to only a few seconds.

Improved depiction of subthalamic nucleus and globus pallidus internus with optimized high-resolution quantitative susceptibility mapping at 7 T

The subthalamic nucleus (STN) and globus pallidus internus (GPi) are commonly used targets in deep-brain stimulation (DBS) surgery for the treatment of movement disorders. The success of DBS critically depends on the spatial precision of stimulation. By taking advantage of good contrast between iron-rich deep-brain nuclei and surrounding tissues, quantitative susceptibility mapping (QSM) has shown promise in differentiating the STN and GPi from the adjacent substantia nigra and globus pallidus externus, respectively. Nonlinear morphology-enabled dipole inversion (NMEDI) is a widely used QSM algorithm, but the image quality of reconstructed susceptibility maps relies on the regularization parameter selection. To date, few studies have systematically optimized the regularization parameter at the ultra-high field of 7 T. In this study, we optimized the regularization parameter in NMEDI to improve the depiction of STN and GPi at different spatial resolutions at both 3 T and 7 T. The optimized QSM images were further compared with other susceptibility-based images, including T2*-weighted (T2*w), R2*, susceptibility-weighted, and phase images. QSM showed better depiction of deep-brain nuclei with clearer boundaries compared with the other methods, and 7 T QSM at 0.35 × 0.35 × 1.0 mm³ demonstrated superior performance to the others. Our findings suggest that optimized high-resolution QSM at 7 T allows for improved delineation of deep-brain nuclei with clear and sharp borders between nuclei, which may become a promising tool for DBS nucleus preoperative localization.

Learned Proximal Networks for Quantitative Susceptibility Mapping

Quantitative Susceptibility Mapping (QSM) estimates tissue magnetic susceptibility distributions from Magnetic Resonance (MR) phase measurements by solving an ill-posed dipole inversion problem. Conventional single orientation QSM methods usually employ regularization strategies to stabilize such inversion, but may suffer from streaking artifacts or over-smoothing. Multiple orientation QSM such as calculation of susceptibility through multiple orientation sampling (COSMOS) can give well-conditioned inversion and an artifact free solution but has expensive acquisition costs. On the other hand, Convolutional Neural Networks (CNN) show great potential for medical image reconstruction, albeit often with limited interpretability. Here, we present a Learned Proximal Convolutional Neural Network (LP-CNN) for solving the ill-posed QSM dipole inversion problem in an iterative proximal gradient descent fashion. This approach combines the strengths of data-driven restoration priors and the clear interpretability of iterative solvers that can take into account the physical model of dipole convolution. During training, our LP-CNN learns an implicit regularizer via its proximal, enabling the decoupling between the forward operator and the data-driven parameters in the reconstruction algorithm. More importantly, this framework is believed to be the first deep learning QSM approach that can naturally handle an arbitrary number of phase input measurements without the need for any ad-hoc rotation or re-training. We demonstrate that the LP-CNN provides state-of-the-art reconstruction results compared to both traditional and deep learning methods while allowing for more flexibility in the reconstruction process.