Background field removal by solving the Laplacian boundary value problem

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.

The appearance of magnetic susceptibility objects in SWI phase depends on object size: Comparison with QSM and CT

PURPOSE

Tissue magnetic susceptibility sign can potentially be detected on susceptibility weighted imaging (SWI) phase (SW-P). This study aims to investigate its performance for depicting brain susceptibility structures.

METHODS

A simulation was performed to depict magnetic susceptibility structures of various geometries on SW-P and quantitative susceptibility mapping (QSM). Brain MRI was performed on 25 subjects using SWI on a 3 T MRI system. QSM was generated from the same data. SW-P and QSM were analyzed according to radiological assessment for depicting globus pallidus nuclei, optic radiation white matter tracts, and lateral ventricular choroid plexus calcifications. In 11 of these subjects, CT was available and correlated with SW-P and QSM to assess their performance in quantifying calcifications in the choroid plexus.

RESULTS

In simulation, the appearance of a sphere on SW-P ranged from centric nodule to mixed positive and negative values as the diameter increased. Large cylinders also appeared as mixed positive and negative values. In comparison, QSM correctly depicted the susceptibility distribution of all magnetic structures. On human brain images, SW-P depicted the globus pallidus and optic radiation with mixed positive and negative values, consistent with simulation, and small choroid plexus calcifications as either mixed positive and negative values or as centric nodules; QSM depicted all structures as solid structures with the expected signs. For measuring calcification in the choroid plexus, QSM vs CT linear regression had a higher coefficient of determination compared to SW-P vs CT and SW-P vs QSM.

CONCLUSION

Appearance of susceptibility sources on SW-P changes with object size. This problem can be overcome using QSM.

Delineation of prostatic calcification using quantitative susceptibility mapping: Spatial accuracy for magnetic resonance-only radiotherapy planning

To investigate the spatial accuracy of delineating prostatic calcifications by quantitative susceptibility mapping (QSM) in comparison with computed tomography (CT), we conducted phantom and human studies. Five differently‐sized spherical hydroxyapatites mimicking prostatic calcification (pseudo‐calcification) were arranged in the order of their sizes at the center of a plastic container filled with gelatin. This calcification phantom underwent magnetic resonance (MR) imaging, including the multiple spoiled gradient‐echo sequences (SPGR) for the QSM and CT as a reference. The volume of each pseudo‐calcification and center‐to‐center distance between the pseudo‐calcifications delineated by QSM and CT were measured. In the human study, eight patients with prostate cancer who underwent radiation therapy and had some prostatic calcifications were included. The patients underwent CT and SPGR and modified DIXON sequence for MR‐only simulation. The hybrid QSM processing combined with the complex signals in the SPGR and water and fat fraction maps estimated from the modified DIXON sequence were used to reconstruct the pelvic susceptibility map in humans. The threshold of CT numbers was set at 130 HU, while the QSM images were manually segmented in the calcification phantom and human studies. In the phantom study, there was an excellent agreement in the pseudo‐calcification volumes between QSM and CT (y = 1.02x – 7.38, R² = 0.99). The signal profiles had similar trends in CT and QSM. The center‐to‐center distances between the pseudo‐calcifications in the phantom were also identical in QSM and CT. The calcification volumes were almost identical between the QSM and CT in the human study (y = 0.95x – 9.32, R² = 1.00). QSM can offer geometric and volumetric accuracies to delineate prostatic calcifications, similar to CT. The prostatic calcification delineated by QSM may facilitate image‐guided radiotherapy in the MR‐only simulation workflow.

Investigating Optimal Echo Times for Quantitative Susceptibility Mapping of Basal Ganglia Nuclei in the Healthy Brain

BACKGROUND

Quantitative susceptibility mapping (QSM) technique had been used to measure the magnetic susceptibility of brain tissue in clinical practice. However, QSM presented echo-time (TE) dependence, and an appropriate number of echo-times (nTEs) for QSM became more important to obtain the reliable susceptibility value.

OBJECTIVE

The aim of the study was to explore the optimal nTEs for quantitative susceptibility mapping (QSM) measurements of basal ganglia nuclei in the healthy brain.

METHODS

3D multi-echo enhanced gradient recalled echo T2 star weighted angiography (ESWAN) sequence was acquired on a 3.0T MR scanner for QSM analysis. Regions of interests (ROIs) were drawn along the margin of the head of the caudate nucleus (HCN), putamen (Pu) and globus pallidus (GP). The mean susceptibility value and standard deviation of the ROIs were derived from the pixels within each region.

RESULTS

CV analysis demonstrated that TE6, TE8 and TE14 ESWAN sequences presented consistent lower CV value (< 1) for QSM measure of HCN, Pu and GP. ANOVA identified that susceptibility value showed no significant difference between TE6 and TE8 in HCN, Pu and GP (P > 0.05). ICC analysis demonstrated that the susceptibility value of TE6-TE8 had the highest ICC value as compared with TE6-TE14 and TE8-TE14 in HCN, Pu and GP. Combined with the timeefficiency of MRI scanning, TE6 sequence could not only provide the reliable QSM measurement but also short imaging time.

CONCLUSION

The current study identified that the optimal nTEs of ESWAN were 6 TEs (2.9ms ~ 80.9ms) for QSM measurement of basal ganglia nuclei in the healthy brain.

The central vein sign in multiple sclerosis lesions: Susceptibility relaxation optimization from a routine MRI multiecho gradient echo sequence

BACKGROUND AND PURPOSE

The white matter lesion central vein sign (CVS) is an emerging biomarker for multiple sclerosis (MS) differential diagnosis. Currently, CVS is detected on susceptibility weighted imaging (SWI) with suboptimal contrast. We developed an imaging method called susceptibility relaxation optimization (SRO) to improve CVS visualization.

METHODS

This was a retrospective study of MS patients who had MRI in June 2018 with routine 3D multiecho gradient echo (GRE) and T2-weighted fluid-attenuated inversion recovery (FLAIR) sequences. SRO and SWI images were reconstructed from GRE data. MS lesions were identified on FLAIR image. The CVS detection rate, the image quality score of CVS conspicuity (range 0-3), and central vein-to-lesion contrast were compared between SRO and SWI images.

RESULTS

In 20 MS patients (mean age 45 ± 9 years; 15 women), SRO significantly increased CVS detection rate compared to SWI (53.3%, 274/514 vs. 32.9%, 169/514; p<.001, McNemar's test). The median image quality score for SRO was 2 compared to 1 for SWI (p<.001, Wilcoxon signed-rank test). The median overall image quality score for SRO was 7 compared to 6 for SWI (p = .003; Wilcoxon signed-rank test). Central vein-to-lesion contrast was 0.12 ± 0.12 in SRO compared to 0.031 ± 0.075 in SWI (p<.001, t-test).

CONCLUSIONS

SRO yields better central vein contrast and increases CVS detection rate compared to SWI.

Variability in subthalamic nucleus targeting for deep brain stimulation with 3 and 7 Tesla magnetic resonance imaging

Deep brain stimulation (DBS) of the subthalamic nucleus (STN) is an effective surgical treatment for Parkinson's disease (PD). Side-effects may, however, be induced when the DBS lead is placed suboptimally. Currently, lower field magnetic resonance imaging (MRI) at 1.5 or 3 Tesla (T) is used for targeting. Ultra-high-field MRI (7 T and above) can obtain superior anatomical information and might therefore be better suited for targeting. This study aims to test whether optimized 7 T imaging protocols result in less variable targeting of the STN for DBS compared to clinically utilized 3 T images. Three DBS-experienced neurosurgeons determined the optimal STN DBS target site on three repetitions of 3 T-T2, 7 T-T2*, 7 T-R2* and 7 T-QSM images for five PD patients. The distance in millimetres between the three repetitive coordinates was used as an index of targeting variability and was compared between field strength, MRI contrast and repetition with a Bayesian ANOVA. Further, the target coordinates were registered to MNI space, and anatomical coordinates were compared between field strength, MRI contrast and repetition using a Bayesian ANOVA. The results indicate that the neurosurgeons are stable in selecting the DBS target site across MRI field strength, MRI contrast and repetitions. The analysis of the coordinates in MNI space however revealed that the actual selected location of the electrode is seemingly more ventral when using the 3 T scan compared to the 7 T scans.

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

OBJECTIVE

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

METHODS

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

RESULTS

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

CONCLUSION

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

Regional brain iron and gene expression provide insights into neurodegeneration in Parkinson's disease

The mechanisms responsible for the selective vulnerability of specific neuronal populations in Parkinson's disease are poorly understood. Oxidative stress secondary to brain iron accumulation is one postulated mechanism. We measured iron deposition in 180 cortical regions of 96 patients with Parkinson's disease and 35 control subjects using quantitative susceptibility mapping. We estimated the expression of 15 745 genes in the same regions using transcriptomic data from the Allen Human Brain Atlas. Using partial least squares regression, we then identified the profile of gene transcription in the healthy brain that underlies increased cortical iron in patients with Parkinson's disease relative to controls. Applying gene ontological tools, we investigated the biological processes and cell types associated with this transcriptomic profile and identified the sets of genes with spatial expression profiles in control brains that correlated significantly with the spatial pattern of cortical iron deposition in Parkinson's disease. Gene ontological analyses revealed that these genes were enriched for biological processes relating to heavy metal detoxification, synaptic function and nervous system development and were predominantly expressed in astrocytes and glutamatergic neurons. Furthermore, we demonstrated that the genes differentially expressed in Parkinson's disease are associated with the pattern of cortical expression identified in this study. Our findings provide mechanistic insights into regional selective vulnerabilities in Parkinson's disease, particularly the processes involving iron accumulation.

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.

Preconditioned water-fat total field inversion: Application to spine quantitative susceptibility mapping

PURPOSE

To (a) develop a preconditioned water-fat total field inversion (wfTFI) algorithm that directly estimates the susceptibility map from complex multi-echo gradient echo data for water-fat regions and to (b) evaluate the performance of the proposed wfTFI quantitative susceptibility mapping (QSM) method in comparison with a local field inversion (LFI) method and a linear total field inversion (TFI) method in the spine.

METHODS

Numerical simulations and in vivo spine multi-echo gradient echo measurements were performed to compare wfTFI to an algorithm based on disjoint background field removal (BFR) and LFI and to a formerly proposed TFI algorithm. The data from 1 healthy volunteer and 10 patients with metastatic bone disease were included in the analysis. Clinical routine computed tomography (CT) images were used as a reference standard to distinguish osteoblastic from osteolytic changes. The ability of the QSM methods to distinguish osteoblastic from osteolytic changes was evaluated.

RESULTS

The proposed wfTFI method was able to decrease the normalized root mean square error compared to the LFI and TFI methods in the simulation. The in vivo wfTFI susceptibility maps showed reduced BFR artifacts, noise amplification, and streaking artifacts compared to the LFI and TFI maps. wfTFI provided a significantly higher diagnostic confidence in differentiating osteolytic and osteoblastic lesions in the spine compared to the LFI method (p = .012).

CONCLUSION

The proposed wfTFI method can minimize BFR artifacts, noise amplification, and streaking artifacts in water-fat regions and can thus better differentiate between osteoblastic and osteolytic changes in patients with metastatic disease compared to LFI and the original TFI method.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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

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.

Brain iron assessment in patients with First-episode schizophrenia using quantitative susceptibility mapping

•

Patients with first-episode schizophrenia had significantly decreased QSM values in the bilateral substantia nigra, left red nucleus and left thalamus.

•

Patients with first-episode schizophrenia had significantly increased regional volumes in the bilateral putamen and bilateral substantia nigra.

•

QSM provides superior sensitivity over R₂* mapping in the evaluation of schizophrenia-related iron alterations.

•

QSM values in regions that showed intergroup differences did not exhibited significant correlations with PANSS scores.

Purpose

Decreased serum ferritin level was recently found in schizophrenia. Whether the brain iron concentration in schizophrenia exists abnormality is of research significance. Quantitative susceptibility mapping (QSM) was used in this study to assess brain iron changes in the grey matter nuclei of patients with first-episode schizophrenia.

Methods

The local ethics committee approved the study, and all subjects gave written informed consent. Thirty patients with first-episode schizophrenia and 30 age and gender-matched healthy controls were included in this study. QSM and effective transverse relaxation rate (R₂*) maps were reconstructed from a three-dimensional multi-echo gradient-echo sequence. The inter-group differences of regional QSM values, R₂* values and volumes were calculated in the grey matter nuclei, including bilateral caudate nucleus, putamen, globus pallidus, substantia nigra, red nucleus, and thalamus. The diagnostic performance of QSM and R₂* was evaluated using receiver operating characteristic curve. The correlations between regional iron variations and clinical PANSS (Positive and Negative Syndrome Scale) scores were assessed using partial correlation analysis.

Results

Compared to healthy controls, patients with first-episode schizophrenia had significantly decreased QSM values (less paramagnetic) in the bilateral substantia nigra, left red nucleus and left thalamus (p < 0.05, FDR correction). QSM proved more sensitive than R₂* regarding inter-group differences. The highest diagnostic performance for first-episode schizophrenia was observed in QSM value of the left substantia nigra (area under the curve, AUC = 0.718, p = 0.004). Regional volumes of bilateral putamen and bilateral substantia nigra were increased (p < 0.05, FDR correction) in first-episode schizophrenia. However, both QSM and R₂* values did not show significant correlations with PANSS scores (p > 0.05).

Conclusion

This study reveals decreased iron concentration in grey matter nuclei of patients with first-episode schizophrenia. QSM provides superior sensitivity over R₂* in the evaluation of schizophrenia-related brain iron changes. It demonstrated that QSM may be a potential biomarker for further understanding the pathophysiological mechanism of first-episode schizophrenia.

A role for magnetic susceptibility in synthetic computed tomography

PURPOSE

Radiotherapy treatment planning based on magnetic resonance imaging (MRI) benefits from increased soft-tissue contrast and functional imaging. MRI-only planning is attractive but limited by the lack of electron density information required for dose calculation, and the difficulty to differentiate air and bone. MRI can map magnetic susceptibility to separate bone from air. A method is introduced to produce synthetic CT (sCT) through automatic voxel-wise assignment of CT numbers from an MRI dataset processed that includes magnetic susceptibility mapping.

METHODS

Volumetric multi-echo gradient echo datasets were acquired in the heads of five healthy volunteers and fourteen patients with cancer using a 3 T MRI system. An algorithm for CT synthesis was designed using the volunteer data, based on fuzzy c-means clustering and adaptive thresholding of the MR data (magnitude, fat, water, and magnetic susceptibility). Susceptibility mapping was performed using a modified version of the iterative phase replacement algorithm. On patient data, the algorithm was assessed by direct comparison to X-ray computed tomography (CT) scans.

RESULTS

The skull, spine, teeth, and major sinuses were clearly distinguished in all sCT, from healthy volunteers and patients. The mean absolute CT number error between X-ray CT and sCT in patients ranged from 78 and 134 HU.

CONCLUSION

Susceptibility mapping using MRI can differentiate air and bone for CT synthesis. The proposed method is automated, fast, and based on a commercially available MRI pulse sequence. The method avoids registration errors and does not rely on a priori information, making it suitable for nonstandard anatomy.

Computed inverse MRI (CIMRI) for intrinsic brain magnetic susceptibility mapping

In magnetic resonance imaging (MRI), tissue magnetization in the main field B₀ is a necessary preparation for magnetic resonance signal formation that imposes an inherent dipole effect on MRI signals, which predisposes an artifact on tissue MRI. In the MRI principle, T2*-weighted MRI can be described by a cascade of data transformations: from the source of tissue magnetic susceptibility (denoted by χ) to the output of complex-valued T2* image (in a magnitude and phase pair). Under the linear approximation of the T2* phase MRI, we can computationally reconstruct the source χ by quantitative susceptibility mapping (QSM), which is an inverse solution that is modeled by computed inverse MRI (CIMRI). For a brain function study using MRI (fMRI), we can reconstruct a timeseries of brain χ images to represent the intrinsic brain function activity called functional QSM (fQSM). This intrinsic depiction is defined as the removal of the artifactual dipole effect and other MRI-introduced distortions from phase data through inverse mapping. With one high-resolution QSM experiment and one group (20 subjects) low-resolution fQSM experiment, we show that the dipole effect manifests as ripples around vessels and a spatial split at a local activation blob and that the dipole effect could be removed by CIMRI. In the context of inverse imaging or undoing MRI transformations (including dipole convolution), we computationally achieve brain intrinsic structural depiction by QSM and intrinsic functional depiction by fQSM.

Iterative Restoration of the Fringe Phase (REFRASE) for QSM

In quantitative susceptibility mapping (QSM), reconstructed results can be critically biased by misinterpreted or missing phase data near the edges of the brain support originating from the non-local relationship between field and susceptibility. These data either have to be excluded or corrected before further processing can take place. To address this, our iterative restoration of the fringe phase (REFRASE) approach simultaneously enhances the accuracy of multi-echo phase data QSM maps and the extent of the area available for evaluation. Data loss caused by strong local phase gradients near the surface of the brain support is recovered within the original phase data using harmonic and dipole-based fields extrapolated from a robust support region toward an extended brain mask. Over several iterations, phase data are rectified prior to the application of further QSM processing steps. The concept is successfully validated on numerical phantoms and brain scans from a cohort of volunteers. The increased extent of the mask and improved numerical stability within the segmented globus pallidus confirm the efficacy of the presented method in comparison to traditional evaluation.

The spatiotemporal changes in dopamine, neuromelanin and iron characterizing Parkinson's disease

In Parkinson's disease, there is a progressive reduction in striatal dopaminergic function, and loss of neuromelanin-containing dopaminergic neurons and increased iron deposition in the substantia nigra. We tested the hypothesis of a relationship between impairment of the dopaminergic system and changes in the iron metabolism. Based on imaging data of patients with prodromal and early clinical Parkinson's disease, we assessed the spatiotemporal ordering of such changes and relationships in the sensorimotor, associative and limbic territories of the nigrostriatal system. Patients with Parkinson's disease (disease duration < 4 years) or idiopathic REM sleep behaviour disorder (a prodromal form of Parkinson's disease) and healthy controls underwent longitudinal examination (baseline and 2-year follow-up). Neuromelanin and iron sensitive MRI and dopamine transporter single-photon emission tomography were performed to assess nigrostriatal levels of neuromelanin, iron, and dopamine. For all three functional territories of the nigrostriatal system, in the clinically most and least affected hemispheres separately, the following was performed: cross-sectional and longitudinal inter-group difference analysis of striatal dopamine and iron, and nigral neuromelanin and iron; in Parkinson's disease patients, exponential fitting analysis to assess the duration of the prodromal phase and the temporal ordering of changes in dopamine, neuromelanin or iron relative to controls; voxel-wise correlation analysis to investigate concomitant spatial changes in dopamine-iron, dopamine-neuromelanin and neuromelanin-iron in the substantia nigra pars compacta. The temporal ordering of dopaminergic changes followed the known spatial pattern of progression involving first the sensorimotor, then the associative and limbic striatal and nigral regions. Striatal dopaminergic denervation occurred first followed by abnormal iron metabolism and finally neuromelanin changes in the substantia nigra pars compacta, which followed the same spatial and temporal gradient observed in the striatum but shifted in time. In conclusion, dopaminergic striatal dysfunction and cell loss in the substantia nigra pars compacta are interrelated with increased nigral iron content.

Detection of Cerebral Microbleeds With Venous Connection at 7-Tesla MRI

OBJECTIVE

Cerebral microbleeds (MBs) are a common finding in patients with cerebral small vessel disease (CSVD) and Alzheimer disease as well as in healthy elderly people, but their pathophysiology remains unclear. To investigate a possible role of veins in the development of MBs, we performed an exploratory study, assessing in vivo presence of MBs with a direct connection to a vein.

METHODS

7-Tesla (7T) MRI was conducted and MBs were counted on quantitative susceptibility mapping (QSM). A submillimeter resolution QSM-based venogram allowed identification of MBs with a direct spatial connection to a vein.

RESULTS

A total of 51 people (mean age [SD] 70.5 [8.6] years, 37% female) participated in the study: 20 had CSVD (cerebral amyloid angiopathy [CAA] with strictly lobar MBs [n = 8], hypertensive arteriopathy [HA] with strictly deep MBs [n = 5], or mixed lobar and deep MBs [n = 7], 72.4 [6.1] years, 30% female) and 31 were healthy controls (69.4 [9.9] years, 42% female). In our cohort, we counted a total of 96 MBs with a venous connection, representing 14% of all detected MBs on 7T QSM. Most venous MBs (86%, n = 83) were observed in lobar locations and all of these were cortical. Patients with CAA showed the highest ratio of venous to total MBs (19%) (HA = 9%, mixed = 18%, controls = 5%).

CONCLUSION

Our findings establish a link between cerebral MBs and the venous vasculature, pointing towards a possible contribution of veins to CSVD in general and to CAA in particular. Pathologic studies are needed to confirm our observations.

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.

QSM in canine model of acute cerebral ischemia: A pilot study

PURPOSE

In the present study, we investigated the potential of QSM to assess the physiological state of cortical tissue in the middle cerebral artery occlusion canine model of a cerebral ischemia.

METHODS

Experiments were performed in 8 anesthetized canines. Gradient echo, perfusion, and DWI data of brains at normal and ischemic states were acquired. In the postprocessed susceptibility and quantitative cerebral blood flow maps, changes in values within the middle cerebral artery-fed cortical territories were quantified both on the ischemic and normal contralateral hemisphere side.

RESULTS

QSM values in critically ischemic tissue were significantly different from contralateral values-namely, susceptibility increase was observed in the cases in which cerebral perfusion was maintained above the threshold of neuronal death. Furthermore, the data indicates presence of a significant correlation between the changes in susceptibility values, cerebral perfusion, and the infarct volume and pial collateral scores. Additionally, our data suggests that difference in cortical susceptibility is prospectively indicative of the infarct growth rate.

CONCLUSION

In an experimental permanent middle cerebral artery occlusion model, QSM was shown to correlate with the functional parameters characterizing viability of ischemic tissue, thus warranting further research on its ability to provide complementary information during acute stroke MRI examinations in humans.

Susceptibility Weighted Imaging for evaluation of musculoskeletal lesions

The presence of blood or calcium in the musculoskeletal (MSK) system may be linked to specific pathological conditions. The ability of MRI for calcium detection is usually limited compared with other techniques such as CT. In a similar manner, the accuracy of MRI for detection and evaluation of hemorrhage in soft tissues is closely linked to the degree of degradation of blood products. Blood and calcium are substances that cause local inhomogeneity of the magnetic field resulting in susceptibility artifacts. To try to evaluate these substances, specific MRI sequences which are highly sensitive to these local magnetic field inhomogeneities such as Susceptibility Weighted Imaging (SWI) have been developed and successfully applied in the Central Nervous System, but scarcely used in MSK. SWI may increase the overall sensitivity of MRI to detect blood and calcium in several clinical scenarios such as degenerative joint disease or bone and soft tissue lesion assessment and discriminate between both compounds, something which is not always possible with conventional MRI approaches. In this paper, physical basis and technical adjustment for SWI acquisition at MSK are detailed reviewing the potential application of SWI in different MSK clinical scenarios.

The traveling heads 2.0: Multicenter reproducibility of quantitative imaging methods at 7 Tesla

OBJECT

This study evaluates inter-site and intra-site reproducibility at ten different 7 T sites for quantitative brain imaging.

MATERIAL AND METHODS

Two subjects - termed the "traveling heads" - were imaged at ten different 7 T sites with a harmonized quantitative brain MR imaging protocol. In conjunction with the system calibration, MP2RAGE, QSM, CEST and multi-parametric mapping/relaxometry were examined.

RESULTS

Quantitative measurements with MP2RAGE showed very high reproducibility across sites and subjects, and errors were in concordance with previous results and other field strengths. QSM had high inter-site reproducibility for relevant subcortical volumes. CEST imaging revealed systematic differences between the sites, but reproducibility was comparable to results in the literature. Relaxometry had also very high agreement between sites, but due to the high sensitivity, differences caused by different applications of the B1 calibration of the two RF coil types used were observed.

CONCLUSION

Our results show that quantitative brain imaging can be performed with high reproducibility at 7 T and with similar reliability as found at 3 T for multicenter studies of the supratentorial brain.

Iron accumulation in the oculomotor nerve of the progressive supranuclear palsy brain

Abnormal iron accumulation around the substantia nigra (SN) is a diagnostic indicator of Parkinsonism. This study aimed to identify iron-related microarchitectural changes around the SN of brains with progressive supranuclear palsy (PSP) via postmortem validations and in vivo magnetic resonance imaging (MRI). 7 T high-resolution MRI was applied to two postmortem brain tissues, from one normal brain and one PSP brain. Histopathological examinations were performed to demonstrate the molecular origin of the high-resolution postmortem MRI findings, by using ferric iron staining, myelin staining, and two-dimensional laser ablation-inductively coupled plasma-mass spectrometry (LA-ICP-MS) imaging. In vivo iron-related MRI was performed on five healthy controls, five patients with Parkinson’s disease (PD), and five patients with PSP. In the postmortem examination, excessive iron deposition along the myelinated fiber at the anterior SN and third cranial nerve (oculomotor nerve) fascicles of the PSP brain was verified by LA-ICP-MS. This region corresponded to those with high R₂^* values and positive susceptibility from quantitative susceptibility mapping (QSM), but was less sensitive in Perls’ Prussian blue staining. In in vivo susceptibility-weighted imaging, hypointense pixels were observed in the region between the SN and red nucleus (RN) in patients with PSP, but not in healthy controls and patients with PD. R₂^* and QSM values of such region were significantly higher in patients with PSP compared to those in healthy controls and patients with PD as well (vs. healthy control: p = 0.008; vs. PD: p = 0.008). Thus, excessive iron accumulation along the myelinated fibers at the anterior SN and oculomotor nerve fascicles may be a pathological characteristic and crucial MR biomarker in a brain with PSP.

SEPIA-Susceptibility mapping pipeline tool for phase images

Quantitative susceptibility mapping (QSM) is a physics-driven computational technique that has a high sensitivity in quantifying iron deposition based on MRI phase images. Furthermore, it has a unique ability to distinguish paramagnetic and diamagnetic contributions such as haemorrhage and calcification based on image contrast. These properties have contributed to a growing interest to use QSM not only in research but also in clinical applications. However, it is challenging to obtain high quality susceptibility map because of its ill-posed nature, especially for researchers who have less experience with QSM and the optimisation of its pipeline. In this paper, we present an open-source processing pipeline tool called SuscEptibility mapping PIpeline tool for phAse images (SEPIA) dedicated to the post-processing of MRI phase images and QSM. SEPIA connects various QSM toolboxes freely available in the field to offer greater flexibility in QSM processing. It also provides an interactive graphical user interface to construct and execute a QSM processing pipeline, simplifying the workflow in QSM research. The extendable design of SEPIA also allows developers to deploy their methods in the framework, providing a platform for developers and researchers to share and utilise the state-of-the-art methods in QSM.

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

High-quality Quantitative Susceptibility Mapping (QSM) with Nonlinear Dipole Inversion (NDI) is developed with pre-determined regularization while matching the image quality of state-of-the-art reconstruction techniques and avoiding over-smoothing that these techniques often suffer from. NDI is flexible enough to allow for reconstruction from an arbitrary number of head orientations and outperforms COSMOS even when using as few as 1-direction data. This is made possible by a nonlinear forward-model that uses the magnitude as an effective prior, for which we derived a simple gradient descent update rule. We synergistically combine this physics-model with a Variational Network (VN) to leverage the power of deep learning in the VaNDI algorithm. This technique adopts the simple gradient descent rule from NDI and learns the network parameters during training, hence requires no additional parameter tuning. Further, we evaluate NDI at 7 T using highly accelerated Wave-CAIPI acquisitions at 0.5 mm isotropic resolution and demonstrate high-quality QSM from as few as 2-direction data.

An optimized and highly repeatable MRI acquisition and processing pipeline for quantitative susceptibility mapping in the head-and-neck region

PURPOSE

Quantitative Susceptibility Mapping (QSM) is an emerging technique sensitive to disease-related changes including oxygenation. It is extensively used in brain studies and has increasing clinical applications outside the brain. Here we present the first MRI acquisition protocol and QSM pipeline optimized for the head-and-neck region together with a repeatability analysis performed in healthy volunteers.

METHODS

We investigated both the intrasession and the intersession repeatability of the optimized method in 10 subjects. We also implemented two, Tikhonov-regularisation-based susceptibility calculation techniques that were found to have higher contrast-to-noise than existing methods in the head-and-neck region. Repeatability was evaluated by calculating the distributions of susceptibility differences between repeated scans and the corresponding minimum detectable effect sizes (MDEs).

RESULTS

Deep brain regions had higher QSM repeatability than neck regions. As expected, intrasession repeatability was generally better than intersession repeatability. Susceptibility maps calculated using projection onto dipole fields for background field removal were more repeatable than using the Laplacian boundary value method in the head-and-neck region. Small (short-axis diameter <5 mm) lymph nodes had the lowest repeatability (MDE = 0.27 ppm) as imperfect segmentation included some of the surrounding paramagnetic fatty fascia, highlighting the importance of accurate region delineation. MDEs in the larger lymph nodes (0.16 ppm), submandibular glands (0.10 ppm), and especially the parotid glands (0.06 ppm) were much lower, comparable to those of the brain regions.

CONCLUSIONS

The high repeatability of the acquisition and pipeline optimized for QSM will facilitate clinical studies in the head-and-neck region.

Feasibility of quantitative susceptibility mapping (QSM) of the human kidney

Objective

To evaluate the feasibility of in-vivo quantitative susceptibility mapping (QSM) of the human kidney.

Methods

An axial single-breath-hold 3D multi-echo sequence (acquisition time 33 s) was completed on a 3 T-MRI-scanner (Magnetom Prisma, Siemens Healthineers, Erlangen, Germany) in 19 healthy volunteers. Graph-cut-based unwrapping combined with the T₂*-IDEAL approach was performed to remove the chemical shift of fat and to quantify QSM of the upper abdomen. Mean susceptibility values of the entire, renal cortex and medulla in both kidneys and the liver were determined and compared. Five subjects were measured twice to examine the reproducibility. One patient with severe renal fibrosis was included in the study to evaluate the potential clinical relevance of QSM.

Results

QSM was successful in 17 volunteers and the patient with renal fibrosis. Anatomical structures in the abdomen were clearly distinguishable by QSM and the susceptibility values obtained in the liver were comparable to those found in the literature. The results showed a good reproducibility. Besides, the mean renal QSM values obtained in healthy volunteers (0.04 ± 0.07 ppm for the right and − 0.06 ± 0.19 ppm for the left kidney) were substantially higher than that measured in the investigated fibrotic kidney (− 0.43 ± − 0.02 ppm).

Conclusion

QSM of the human kidney could be a promising approach for the assessment of information about microscopic renal tissue structure. Therefore, it might further improve functional renal MR imaging.

Electronic supplementary material

The online version of this article (10.1007/s10334-020-00895-9) contains supplementary material, which is available to authorized users.

Methods for quantitative susceptibility and R2* mapping in whole post-mortem brains at 7T applied to amyotrophic lateral sclerosis

Susceptibility weighted magnetic resonance imaging (MRI) is sensitive to the local concentration of iron and myelin. Here, we describe a robust image processing pipeline for quantitative susceptibility mapping (QSM) and R2* mapping of fixed post-mortem, whole-brain data. Using this pipeline, we compare the resulting quantitative maps in brains from patients with amyotrophic lateral sclerosis (ALS) and controls, with validation against iron and myelin histology. Twelve post-mortem brains were scanned with a multi-echo gradient echo sequence at 7T, from which susceptibility and R2* maps were generated. Semi-quantitative histological analysis for ferritin (the principal iron storage protein) and myelin proteolipid protein was performed in the primary motor, anterior cingulate and visual cortices. Magnetic susceptibility and R2* values in primary motor cortex were higher in ALS compared to control brains. Magnetic susceptibility and R2* showed positive correlations with both myelin and ferritin estimates from histology. Four out of nine ALS brains exhibited clearly visible hyperintense susceptibility and R2* values in the primary motor cortex. Our results demonstrate the potential for MRI-histology studies in whole, fixed post-mortem brains to investigate the biophysical source of susceptibility weighted MRI signals in neurodegenerative diseases like ALS.

3D magnetic resonance fingerprinting with quadratic RF phase

PURPOSE

To implement 3D magnetic resonance fingerprinting (MRF) with quadratic RF phase (qRF-MRF) for simultaneous quantification of T₁ , T₂ , ΔB₀ , and T 2 ∗ .

METHODS

3D MRF data with effective undersampling factor of 3 in the slice direction were acquired with quadratic RF phase patterns for T₁ , T₂ , and T 2 ∗ sensitivity. Quadratic RF phase encodes the off-resonance by modulating the on-resonance frequency linearly in time. Transition to 3D brings practical limitations for reconstruction and dictionary matching because of increased data and dictionary sizes. Randomized singular value decomposition (rSVD)-based compression in time and reduction in dictionary size with a quadratic interpolation method are combined to be able to process prohibitively large data sets in feasible reconstruction and matching times.

RESULTS

Accuracy of 3D qRF-MRF maps in various resolutions and orientations are compared to 3D fast imaging with steady-state precession (FISP) for T₁ and T₂ contrast and to 2D qRF-MRF for T 2 ∗ contrast and ΔB₀ . The precision of 3D qRF-MRF was 1.5-2 times higher than routine clinical scans. 3D qRF-MRF ΔB₀ maps were further processed to highlight the susceptibility contrast.

CONCLUSION

Natively co-registered 3D whole brain T₁ , T₂ , T 2 ∗ , ΔB₀ , and QSM maps can be acquired in as short as 5 min with 3D qRF-MRF.

Improved body quantitative susceptibility mapping by using a variable-layer single-min-cut graph-cut for field-mapping

PURPOSE

To develop a robust algorithm for field-mapping in the presence of water-fat components, large B 0 field inhomogeneities and MR signal voids and to apply the developed method in body applications of quantitative susceptibility mapping (QSM).

METHODS

A framework solving the cost-function of the water-fat separation problem in a single-min-cut graph-cut based on the variable-layer graph construction concept was developed. The developed framework was applied to a numerical phantom enclosing an MR signal void, an air bubble experimental phantom, 14 large field of view (FOV) head/neck region in vivo scans and to 6 lumbar spine in vivo scans. Field-mapping and subsequent QSM results using the proposed algorithm were compared to results using an iterative graph-cut algorithm and a formerly proposed single-min-cut graph-cut.

RESULTS

The proposed method was shown to yield accurate field-map and susceptibility values in all simulation and in vivo datasets when compared to reference values (simulation) or literature values (in vivo). The proposed method showed improved field-map and susceptibility results compared to iterative graph-cut field-mapping especially in regions with low SNR, strong field-map variations and high R 2 ∗ values.

CONCLUSIONS

A single-min-cut graph-cut field-mapping method with a variable-layer construction was developed for field-mapping in body water-fat regions, improving quantitative susceptibility mapping particularly in areas close to MR signal voids.

Different iron deposition patterns in Parkinson's disease subtypes: a quantitative susceptibility mapping study

Background

Parkinson's disease (PD) is a heterogeneous neurodegenerative disorder with various subtypes and uncertain pathogenesis. Iron deposition is considered to be involved in the pathological mechanisms of PD. The present study aimed to investigate the iron deposition in deep gray matter in patients with different PD subtypes using quantitative susceptibility mapping (QSM).

Methods

Forty-six PD patients and 22 healthy controls (HCs) were recruited for the study. PD patients were allocated to the tremor-dominant (TD) group (n=22), postural instability and gait disorder-dominant (PIGD) group (n=19), and intermediate group (n=5). Susceptibility values in deep gray matter nuclei measured by QSM among the PD-TD and PD-PIGD groups and the HCs, as well as the relationship between iron accumulation and clinical motor features, were investigated.

Results

Susceptibility values in the dentate nucleus (DN) were greater in the PD-TD (118.73±70.45) group than in the PD-PIGD (72.14±39.85, P=0.02) group and HCs (78.26±41.38, P=0.042). Further, a significant positive correlation was observed between the DN susceptibility values and tremor scores (r=0.324, P=0.028). Compared with the HCs (182.60±85.35), both the PD-TD (282.00±102.49, P=0.006) and PD-PIGD groups (284.91±118.54, P=0.007) exhibited greater susceptibility values in the substantia nigra (SN) pars reticulata. The susceptibility values in the SN pars compacta were also greater in the PD-PIGD group (164.51±89.44) than in the HCs (107.78±63.11, P=0.048).

Conclusions

The present study demonstrated various iron deposition patterns in different PD phenotypes. These findings give insight into the pathophysiology underlying different PD phenotypes, and potentially illustrate the involvement of iron deposition in the PD-TD and PD-PIGD subtypes.

Methodological Considerations for Neuroimaging in Deep Brain Stimulation of the Subthalamic Nucleus in Parkinson's Disease Patients

Deep brain stimulation (DBS) of the subthalamic nucleus is a neurosurgical intervention for Parkinson's disease patients who no longer appropriately respond to drug treatments. A small fraction of patients will fail to respond to DBS, develop psychiatric and cognitive side-effects, or incur surgery-related complications such as infections and hemorrhagic events. In these cases, DBS may require recalibration, reimplantation, or removal. These negative responses to treatment can partly be attributed to suboptimal pre-operative planning procedures via direct targeting through low-field and low-resolution magnetic resonance imaging (MRI). One solution for increasing the success and efficacy of DBS is to optimize preoperative planning procedures via sophisticated neuroimaging techniques such as high-resolution MRI and higher field strengths to improve visualization of DBS targets and vasculature. We discuss targeting approaches, MRI acquisition, parameters, and post-acquisition analyses. Additionally, we highlight a number of approaches including the use of ultra-high field (UHF) MRI to overcome limitations of standard settings. There is a trade-off between spatial resolution, motion artifacts, and acquisition time, which could potentially be dissolved through the use of UHF-MRI. Image registration, correction, and post-processing techniques may require combined expertise of traditional radiologists, clinicians, and fundamental researchers. The optimization of pre-operative planning with MRI can therefore be best achieved through direct collaboration between researchers and clinicians.

Multi-centre, multi-vendor reproducibility of 7T QSM and R2* in the human brain: Results from the UK7T study

Introduction

We present the reliability of ultra-high field T₂* MRI at 7T, as part of the UK7T Network's “Travelling Heads” study. T₂*-weighted MRI images can be processed to produce quantitative susceptibility maps (QSM) and R₂* maps. These reflect iron and myelin concentrations, which are altered in many pathophysiological processes. The relaxation parameters of human brain tissue are such that R₂* mapping and QSM show particularly strong gains in contrast-to-noise ratio at ultra-high field (7T) vs clinical field strengths (1.5–3T). We aimed to determine the inter-subject and inter-site reproducibility of QSM and R₂* mapping at 7T, in readiness for future multi-site clinical studies.

Methods

Ten healthy volunteers were scanned with harmonised single- and multi-echo T₂*-weighted gradient echo pulse sequences. Participants were scanned five times at each “home” site and once at each of four other sites. The five sites had 1× Philips, 2× Siemens Magnetom, and 2× Siemens Terra scanners. QSM and R₂* maps were computed with the Multi-Scale Dipole Inversion (MSDI) algorithm (https://github.com/fil-physics/Publication-Code). Results were assessed in relevant subcortical and cortical regions of interest (ROIs) defined manually or by the MNI152 standard space.

Results and Discussion

Mean susceptibility (χ) and R₂* values agreed broadly with literature values in all ROIs. The inter-site within-subject standard deviation was 0.001–0.005 ppm (χ) and 0.0005–0.001 ms⁻¹ (R₂*). For χ this is 2.1–4.8 fold better than 3T reports, and 1.1–3.4 fold better for R₂*. The median ICC from within- and cross-site R₂* data was 0.98 and 0.91, respectively. Multi-echo QSM had greater variability vs single-echo QSM especially in areas with large B₀ inhomogeneity such as the inferior frontal cortex. Across sites, R₂* values were more consistent than QSM in subcortical structures due to differences in B₀-shimming. On a between-subject level, our measured χ and R₂* cross-site variance is comparable to within-site variance in the literature, suggesting that it is reasonable to pool data across sites using our harmonised protocol.

Conclusion

The harmonized UK7T protocol and pipeline delivers on average a 3-fold improvement in the coefficient of reproducibility for QSM and R₂* at 7T compared to previous reports of multi-site reproducibility at 3T. These protocols are ready for use in multi-site clinical studies at 7T.

The 2016 QSM Challenge: Lessons learned and considerations for a future challenge design

PURPOSE

The 4th International Workshop on MRI Phase Contrast and QSM (2016, Graz, Austria) hosted the first QSM Challenge. A single-orientation gradient recalled echo acquisition was provided, along with COSMOS and the χ33 STI component as ground truths. The submitted solutions differed more than expected depending on the error metric used for optimization and were generally over-regularized. This raised (unanswered) questions about the ground truths and the metrics utilized.

METHODS

We investigated the influence of background field remnants by applying additional filters. We also estimated the anisotropic contributions from the STI tensor to the apparent susceptibility to amend the χ33 ground truth and to investigate the impact on the reconstructions. Lastly, we used forward simulations from the COSMOS reconstruction to investigate the impact noise had on the metric scores.

RESULTS

Reconstructions compared against the amended STI ground truth returned lower errors. We show that the background field remnants had a minor impact in the errors. In the absence of inconsistencies, all metrics converged to the same regularization weights, whereas structural similarity index metric was more insensitive to such inconsistencies.

CONCLUSION

There was a mismatch between the provided data and the ground truths due to the presence of unaccounted anisotropic susceptibility contributions and noise. Given the lack of reliable ground truths when using in vivo acquisitions, simulations are suggested for future QSM Challenges.

Decreasing iron susceptibility with temperature in quantitative susceptibility mapping: A phantom study

To clarify the temperature dependence of susceptibility estimated by quantitative susceptibility mapping (QSM) analysis, we investigated the relationship between temperature and susceptibility using a cylinder phantom with varying temperatures. Six solutions with various concentrations of superparamagnetic iron oxide (SPIO) nanoparticles were employed. These tubes were placed in a cylinder phantom and surrounded with water. The temperature of the circulated water was adjusted to change the temperature in the cylinder phantom from 25.8 °C to 42.5 °C. The cylinder phantom was scanned via a three-dimensional multiple spoiled gradient-echo sequence for R2* and QSM analyses with varying temperatures. The relationships between temperature, susceptibility, and R2* values were determined. Moreover, the temperature coefficients of susceptibility (χ-Tc) and (R2*-Tc) were calculated at each concentration and the linearities in these indices against each SPIO concentration were validated. Significant inverse correlations were found between temperature, susceptibility, and R2* values at each SPIO concentration due to the decrease in paramagnetic iron susceptibility that occurred with increasing temperature based on Curie's law. Moreover, although there were significant correlations between the susceptibility and R2* values at any temperature, the slopes of the regression lines grew in height with greater temperatures. The percentage of difference per Celsius degree in susceptibility in any SPIO concentration was lower than the corresponding finding among the R2* results. There were strong linearities between the SPIO concentration, χ-Tc (r = -0.994; p < 0.001), and R2*-Tc (r = -0.998; p < 0.001). The χ-Tc and R2*-Tc outcomes in a particular voxel varied considerably with the iron contents. Although there was an inverse correlation noted between temperature and susceptibility, the susceptibility analysis showed smaller temperature dependence relative to the R2* analysis. QSM analysis might be a more suitable option for magnetic resonance-based iron quantification in comparison with R2* relaxometry.

Simultaneous QSM and metabolic imaging of the brain using SPICE: Further improvements in data acquisition and processing

PURPOSE

To achieve high-resolution mapping of brain tissue susceptibility in simultaneous QSM and metabolic imaging.

METHODS

Simultaneous QSM and metabolic imaging was first achieved using SPICE (spectroscopic imaging by exploiting spatiospectral correlation), but the QSM maps thus obtained were at relatively low-resolution (2.0 × 3.0 × 3.0 mm³ ). We overcome this limitation using an improved SPICE data acquisition method with the following novel features: 1) sampling (k, t)-space in dual densities, 2) sampling central k-space fully to achieve nominal spatial resolution of 3.0 × 3.0 × 3.0 mm³ for metabolic imaging, and 3) sampling outer k-space sparsely to achieve spatial resolution of 1.0 × 1.0 × 1.9 mm³ for QSM. To keep the scan time short, we acquired spatiospectral encodings in echo-planar spectroscopic imaging trajectories in central k-space but in CAIPIRINHA (controlled aliasing in parallel imaging results in higher acceleration) trajectories in outer k-space using blipped phase encodings. For data processing and image reconstruction, a union-of-subspaces model was used, effectively incorporating sensitivity encoding, spatial priors, and spectral priors of individual molecules.

RESULTS

In vivo experiments were carried out to evaluate the feasibility and potential of the proposed method. In a 6-min scan, QSM maps at 1.0 × 1.0 × 1.9 mm³ resolution and metabolic maps at 3.0 × 3.0 × 3.0 mm³ nominal resolution were obtained simultaneously. Compared with the original method, the QSM maps obtained using the new method reveal fine-scale brain structures more clearly.

CONCLUSION

We demonstrated the feasibility of achieving high-resolution QSM simultaneously with metabolic imaging using a modified SPICE acquisition method. The improved capability of SPICE may further enhance its practical utility in brain mapping.

Cerebral oxygen extraction fraction: Comparison of dual-gas challenge calibrated BOLD with CBF and challenge-free gradient echo QSM+qBOLD

PURPOSE

To compare cortical gray matter oxygen extraction fraction (OEF) estimated from 2 MRI methods: (1) the quantitative susceptibility mapping (QSM) plus quantitative blood oxygen level dependent imaging (qBOLD) (QSM+qBOLD or QQ), and (2) the dual-gas calibrated-BOLD (DGCB) in healthy subjects; and to investigate the validity of iso-cerebral metabolic rate of oxygen consumption assumption during hypercapnia using QQ.

METHODS

In 10 healthy subjects, 3 tesla MRI including a multi-echo gradient echo sequence at baseline and hypercapnia for QQ, as well as an EPI dual-echo pseudo-continuous arterial spin labeling for DGCB, were performed under a hypercapnic and a hyperoxic condition. OEFs from QQ and DGCB were compared using region of interest analysis and paired t test. For QQ, cerebral metabolic rate of oxygen consumption = cerebral blood flow*OEF*arterial oxygen content was generated for both baseline and hypercapnia, which were compared.

RESULTS

Average OEF in cortical gray matter across 10 subjects from QQ versus DGCB was 35.5 ± 6.7% versus 38.0 ± 9.1% (P = .49) at baseline and 20.7 ± 4.4% versus 28.4 ± 7.6% (P = .02) in hypercapnia: OEF in cortical gray matter was significantly reduced as measured in QQ (P < .01) and in DGCB (P < .01). Cerebral metabolic rate of oxygen consumption (in μmol O₂ /min/100 g) was 168.2 ± 54.1 at baseline from DGCB and was 153.1 ± 33.8 at baseline and 126.4 ± 34.2 (P < .01) in hypercapnia from QQ.

CONCLUSION

The differences in OEF obtained from QQ and DGCB are small and nonsignificant at baseline but are statistically significant during hypercapnia. In addition, QQ shows a cerebral metabolic rate of oxygen consumption decrease (17.4%) during hypercapnia.

Quantitative Measurement of Metal Accumulation in Brain of Patients With Wilson's Disease

BACKGROUND

Currently, no study has evaluated metal accumulation in the brains of patients with Wilson's disease by using quantitative susceptibility mapping at 3T MRI. The objectives of this study were to qualitatively and quantitatively evaluate changes in magnetic susceptibility and R2* maps in deep gray matter nuclei to discriminate Wilson's disease patients from healthy controls and to evaluate their sensitivities in diagnosing Wilson's disease.

METHODS

Magnetic susceptibility and R2* maps and conventional T1-weighted, T2-weighted, and T2-weighted fluid-attenuated inversion recovery images were obtained from 17 Wilson's disease patients and 14 age-matched healthy controls on a 3T MRI scanner. Differences between Wilson's disease and healthy control groups in susceptibility and R2* values in multiple deep nuclei were evaluated using a Mann-Whitney U test and receiver operating characteristic curves. The correlations of susceptibility and R2* values with Unified Wilson's Disease Rating Scale score were also performed.

RESULTS

Magnetic susceptibility and R2* can effectively distinguish different types of signal abnormalities. Magnetic susceptibility and R2* values in multiple deep nuclei of Wilson's disease patients were significantly higher than those in healthy controls. Magnetic susceptibility value in the substantia nigra had the highest area under the curve (0.888). There were positive correlations of the Unified Wilson's Disease Rating Scale score with susceptibility values in the caudate nucleus (r = 0.757, P = 0.011), putamen (r = 0.679, P = 0.031), and red nucleus (r = 0.638, P = 0.047), as well as R2* values in the caudate nucleus (r = 0.754, P = 0.012).

CONCLUSIONS

Quantitative susceptibility mapping at 3T could be a useful tool to evaluate metal accumulation in deep gray matter nuclei of Wilson's disease patients. © 2020 International Parkinson and Movement Disorder Society.

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

PURPOSE

Multi-echo gradient-recalled echo acquisitions for QSM enable optimizing the SNR for several tissue types through multi-echo (TE) combination or investigating temporal variations in the susceptibility (potentially reflecting tissue microstructure) by calculating one QSM image at each TE (TE-dependent QSM). In contrast with multi-echo QSM, applying Laplacian-based methods (LBMs) for phase unwrapping and background field removal to single TEs could introduce nonlinear temporal variations (independent of tissue microstructure) into the measured susceptibility. Here, we aimed to compare the effect of LBMs on the QSM susceptibilities in TE-dependent versus multi-echo QSM.

METHODS

TE-dependent recalled echo data simulated in a numerical head phantom and gradient-recalled echo images acquired at 3 T in 10 healthy volunteers. Several QSM pipelines were tested, including four distinct LBMs: sophisticated harmonic artifact reduction for phase data (SHARP), variable-radius sophisticated harmonic artifact reduction for phase data (V-SHARP), Laplacian boundary value background field removal (LBV), and one-step total generalized variation (TGV). Results from distinct pipelines were compared using visual inspection, summary statistics of susceptibility in deep gray matter/white matter/venous regions of interest, and, in the healthy volunteers, regional susceptibility bias analysis and nonparametric tests.

RESULTS

Multi-echo versus TE-dependent QSM had higher regional accuracy, especially in high-susceptibility regions and at shorter TEs. Everywhere except in the veins, a processing pipeline incorporating TGV provided the most temporally stable TE-dependent QSM results with an accuracy similar to multi-echo QSM.

CONCLUSIONS

For TE-dependent QSM, carefully choosing LBMs can minimize the introduction of LBM-related nonlinear temporal susceptibility variations.

Correction of motion-induced susceptibility artifacts and B0 drift during proton resonance frequency shift-based MR thermometry in the pelvis with background field removal methods

Purpose

The linear change of the water proton resonance frequency shift (PRFS) with temperature is used to monitor temperature change based on the temporal difference of image phase. Here, the effect of motion‐induced susceptibility artifacts on the phase difference was studied in the context of mild radio frequency hyperthermia in the pelvis.

Methods

First, the respiratory‐induced field variations were disentangled from digestive gas motion in the pelvis. The projection onto dipole fields (PDF) as well as the Laplacian boundary value (LBV) algorithm were applied on the phase difference data to eliminate motion‐induced susceptibility artifacts. Both background field removal (BFR) algorithms were studied using simulations of susceptibility artifacts, a phantom heating experiment, and volunteer and patient heating data.

Results

Respiratory‐induced field variations were negligible in the presence of the filled water bolus. Even though LBV and PDF showed comparable results for most data, LBV seemed more robust in our data sets. Some data sets suggested that PDF tends to overestimate the background field, thus removing phase attributed to temperature. The BFR methods even corrected for susceptibility variations induced by a subvoxel displacement of the phantom. The method yielded successful artifact correction in 2 out of 4 patient treatment data sets during the entire treatment duration of mild RF heating of cervical cancer. The heating pattern corresponded well with temperature probe data.

Conclusion

The application of background field removal methods in PRFS‐based MR thermometry has great potential in various heating applications and body regions to reduce motion‐induced susceptibility artifacts that originate outside the region of interest, while conserving temperature‐induced PRFS. In addition, BFR automatically removes up to a first‐order spatial B₀ drift.

Simultaneous voxel-based magnetic susceptibility and morphometry analysis using magnetization-prepared spoiled turbo multiple gradient echo

This study aimed to develop and test a simultaneous acquisition and analysis pipeline for voxel-based magnetic susceptibility and morphometry (VBMSM) on a single dataset using young volunteers, elderly healthy volunteers, and an Alzheimer's disease (AD) group. 3D T₁ -weighted and multi-echo phase images for VBM and quantitative susceptibility mapping (QSM) were simultaneously acquired using a magnetization-prepared spoiled turbo multiple gradient echo sequence with inversion pulse for QSM (MP-QSM). The magnitude image was split into gray matter (GM) and white matter (WM) and was spatially normalized. The susceptibility map was reconstructed from the phase images. The segmented image and susceptibility map were compared with those obtained from conventional multiple spoiled gradient echo (mGRE) and MP-spoiled gradient echo (MP-GRE) in healthy volunteers to validate the availability of MP-QSM by numerical measurements. To assess the feasibility of the VBMSM analysis pipeline, voxel-based comparisons of susceptibility and morphometry in MP-QSM were conducted in volunteers with a bimodal age distribution, and in elderly volunteers and the AD group, using spatially normalized GM and WM volume images and a susceptibility map. GM/WM contrasts in MP-QSM, MP-GRE, and mGRE were 0.14 ± 0.011, 0.17 ± 0.015, and 0.045 ± 0.010, respectively. Segmented GM and WM volumes in the MP-QSM closely coincided with those in the MP-GRE. Region of interest analyses indicated that the mean susceptibility values in MP-QSM were completely in agreement with those in mGRE. In an evaluation of the aging effect, a significant increase and decrease in susceptibility and volume were found by VBMSM in deep GM and WM, respectively. Between the elderly volunteers and the AD group, the characteristic susceptibility and volume changes in GM and WM were observed. The proposed MP-QSM sequence makes it possible to acquire acceptable-quality images for simultaneous analysis and determine brain atrophy and susceptibility distribution without image registration by using voxel-based analyses.

Elevated brain oxygen extraction fraction measured by MRI susceptibility relates to perfusion status in acute ischemic stroke

Recent clinical trials of new revascularization therapies in acute ischemic stroke have highlighted the importance of physiological imaging to identify optimal treatments for patients. Oxygen extraction fraction (OEF) is a hallmark of at-risk tissue in stroke, and can be quantified from the susceptibility effect of deoxyhemoglobin molecules in venous blood on MRI phase scans. We measured OEF within cerebral veins using advanced quantitative susceptibility mapping (QSM) MRI reconstructions in 20 acute stroke patients. Absolute OEF was elevated in the affected (29.3 ± 3.4%) versus the contralateral hemisphere (25.5 ± 3.1%) of patients with large diffusion-perfusion lesion mismatch (P = 0.032). In these patients, OEF negatively correlated with relative CBF measured by dynamic susceptibility contrast MRI (P = 0.004), suggesting compensation for reduced flow. Patients with perfusion-diffusion match or no hypo-perfusion showed less OEF difference between hemispheres. Nine patients received longitudinal assessment and showed OEF ratio (affected to contralateral) of 1.2 ± 0.1 at baseline that normalized (decreased) to 1.0 ± 0.1 at follow-up three days later (P = 0.03). Our feasibility study demonstrates that QSM MRI can non-invasively quantify OEF in stroke patients, relates to perfusion status, and is sensitive to OEF changes over time. Clinical trial registration: Longitudinal MRI examinations of patients with brain ischemia and blood brain barrier permeability; clinicaltrials.org :NCT02077582.

Improved susceptibility-weighted imaging for high contrast and resolution thalamic nuclei mapping at 7T

PURPOSE

The thalamus is an important brain structure and neurosurgical target, but its constituting nuclei are challenging to image non-invasively. Recently, susceptibility-weighted imaging (SWI) at ultra-high field has shown promising capabilities for thalamic nuclei mapping. In this work, several methodological improvements were explored to enhance SWI quality and contrast, and specifically its ability for thalamic imaging.

METHODS

High-resolution SWI was performed at 7T in healthy participants, and the following techniques were applied: (a) monitoring and retrospective correction of head motion and B₀ perturbations using integrated MR navigators, (b) segmentation and removal of venous vessels on the SWI data using vessel enhancement filtering, and (c) contrast enhancement by tuning the parameters of the SWI phase-magnitude combination. The resulting improvements were evaluated with quantitative metrics of image quality, and by comparison to anatomo-histological thalamic atlases.

RESULTS

Even with sub-millimeter motion and natural breathing, motion and field correction produced clear improvements in both magnitude and phase data quality (76% and 41%, respectively). The improvements were stronger in cases of larger motion/field deviations, mitigating the dependence of image quality on subject performance. Optimizing the SWI phase-magnitude combination yielded substantial improvements in image contrast, particularly in the thalamus, well beyond previously reported SWI results. The atlas comparisons provided compelling evidence of anatomical correspondence between SWI features and several thalamic nuclei, for example, the ventral intermediate nucleus. Vein detection performed favorably inside the thalamus, and vein removal further improved visualization.

CONCLUSION

Altogether, the proposed developments substantially improve high-resolution SWI, particularly for thalamic nuclei imaging.

MRI T2 and T2* relaxometry to visualize neuromelanin in the dorsal substantia nigra pars compacta

Visualizing gradual changes in neuromelanin distribution within the substantia nigra is an important metric used to monitor the progression of Parkinsonism. This study aimed to identify the origin of the mismatch region between magnetic resonance transverse relaxation times (T₂ and T₂*) in the substantia nigra and investigate its feasibility and implications for in vivo detection of neuromelanin as a clinical biomarker. The relationships between neuromelanin distribution assessed by histological staining and the area of T₂ and T₂* mismatch determined by high- and low-resolution magnetic resonance relaxometry at 7T were directly compared in two normal and one depigmented substantia nigra collected at postmortem. In vivo feasibility of assessing T₂ and T₂* mismatch, clinically, was investigated using 3T magnetic resonance imaging. In the normal postmortem substantia nigra tissue, the T₂ and T₂* mismatch region exhibiting a linear pattern was strongly colocalized with neuromelanin distribution along the dorsal substantia nigra pars compacta, but a negligible amount of dorsal mismatch was observed in the depigmented brain. The regions of T₂ and T₂* mismatch from MRI, neuromelanin pigments from histology, and elevated iron signals from mass spectrometry were spatially overlapped for a normal postmortem brain. In preliminary in vivo studies, a similar, linear T₂ and T₂* mismatch region was observed in the dorsal area of the substantia nigra in eight normal subjects; this mismatch was significantly obscured in eight Parkinson's disease patients. The length of the dorsal linear mismatch line based on the T₂*-T₂ mask was significantly shorter in the Parkinson's disease patients compared to normal controls; this result was corroborated by reduced striatal uptake of [¹⁸F] FP-CIT dopamine transporters assessed by positron emission tomography scans. In conclusion, the measurement of T₂ and T₂* mismatch could serve as a complementary imaging biomarker to visualize the dorsal region of the substantia nigra pars compacta, which contains large amounts of neuromelanin.

Automated adaptive preconditioner for quantitative susceptibility mapping

PURPOSE

To develop an automated adaptive preconditioner for QSM reconstruction with improved susceptibility quantification accuracy and increased image quality.

THEORY AND METHODS

The total field was used to rapidly produce an approximate susceptibility map, which was then averaged and trended over R 2 ∗ binning to generate a spatially varying distribution of preconditioning values. This automated adaptive preconditioner was used to reconstruct QSM via total field inversion and was compared with its empirical counterparts in a numerical simulation, a brain experiment with 5 healthy subjects and 5 patients with intracerebral hemorrhage, and a cardiac experiment with 3 healthy subjects.

RESULTS

Among evaluated preconditioners, the automated adaptive preconditioner achieved the fastest convergence in reducing the RMSE of the QSM in the simulation, suppressed hemorrhage-associated artifacts while preserving surrounding brain tissue contrasts, and provided cardiac chamber oxygenation values consistent with those reported in the literature.

CONCLUSION

An automated adaptive preconditioner allows high-quality QSM from the total field in imaging various anatomies with dynamic susceptibility ranges.

Multi-site harmonization of 7 tesla MRI neuroimaging protocols

Increasing numbers of 7 T (7 T) magnetic resonance imaging (MRI) scanners are in research and clinical use. 7 T MRI can increase the scanning speed, spatial resolution and contrast-to-noise-ratio of many neuroimaging protocols, but technical challenges in implementation have been addressed in a variety of ways across sites. In order to facilitate multi-centre studies and ensure consistency of findings across sites, it is desirable that 7 T MRI sites implement common high-quality neuroimaging protocols that can accommodate different scanner models and software versions.

With the installation of several new 7 T MRI scanners in the United Kingdom, the UK7T Network was established with an aim to create a set of harmonized structural and functional neuroimaging sequences and protocols. The Network currently includes five sites, which use three different scanner platforms, provided by two different vendors.

Here we describe the harmonization of functional and anatomical imaging protocols across the three different scanner models, detailing the necessary changes to pulse sequences and reconstruction methods. The harmonized sequences are fully described, along with implementation details. Example datasets acquired from the same subject on all Network scanners are made available. Based on these data, an evaluation of the harmonization is provided. In addition, the implementation and validation of a common system calibration process is described.

Quantitative susceptibility mapping of the spine using in-phase echoes to initialize inhomogeneous field and R2* for the nonconvex optimization problem of fat-water separation

Quantitative susceptibility mapping (QSM) of human spinal vertebrae from a multi-echo gradient-echo (GRE) sequence is challenging, because comparable amounts of fat and water in the vertebrae make it difficult to solve the nonconvex optimization problem of fat-water separation (R2*-IDEAL) for estimating the magnetic field induced by tissue susceptibility. We present an in-phase (IP) echo initialization of R2*-IDEAL for QSM in the spinal vertebrae. Ten healthy human subjects were recruited for spine MRI. A 3D multi-echo GRE sequence was implemented to acquire out-phase and IP echoes. For the IP method, the R2* and field maps estimated by separately fitting the magnitude and phase of IP echoes were used to initialize gradient search R2*-IDEAL to obtain final R2*, field, water, and fat maps, and the final field map was used to generate QSM. The IP method was compared with the existing Zero method (initializing the field to zero), VARPRO-GC (variable projection using graphcuts but still initializing the field to zero), and SPURS (simultaneous phase unwrapping and removal of chemical shift using graphcuts for initialization) on both simulation and in vivo data. The single peak fat model was also compared with the multi-peak fat model. There was no substantial difference on QSM between the single peak and multi-peak fat models, but there were marked differences among different initialization methods. The simulations demonstrated that IP provided the lowest error in the field map. Compared to Zero, VARPRO-GC and SPURS, the proposed IP method provided substantially improved spine QSM in all 10 subjects.

MARBLES - Metal Artifact Based Landmark Enhanced Susceptibility Weighted Imaging For Interventional Device Localization In MRI

Susceptibility Weighted Imaging (SWI) is a method extensively studied for its application to improve contrast in MR imaging modality. The method enhances the visualization of magnetically susceptible content such as iron, calcium and zinc in the tissues by using the susceptibility differences in tissues to generate a unique image contrast. In this study, we propose an SWI based approach to improve the visualization of interventional devices in MRI data. Results obtained from two datasets (biopsy needle and brachytherapy seeds), indicate SWI to be suitable for visualization of the interventional devices, while also being computationally faster when compared with quantitative susceptibility mapping (QSM).