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.

Quantitative magnetic resonance imaging biomarkers for cortical pathology in multiple sclerosis at 7 T.

Substantial cortical gray matter tissue damage, which correlates with clinical disease severity, has been revealed in multiple sclerosis (MS) using advanced magnetic resonance imaging (MRI) methods at 3 T and the use of ultra-high field, as well as in histopathology studies. While clinical assessment mainly focuses on lesions using T 1 - and T 2 -weighted MRI, quantitative MRI (qMRI) methods are capable of uncovering subtle microstructural changes. The aim of this ultra-high field study is to extract possible future MR biomarkers for the quantitative evaluation of regional cortical pathology. Because of their sensitivity to iron, myelin, and in part specifically to cortical demyelination, T 1 , T 2 , R 2 * , and susceptibility mapping were performed including two novel susceptibility markers; in addition, cortical thickness as well as the volumes of 34 cortical regions were computed. Data were acquired in 20 patients and 16 age- and sex-matched healthy controls. In 18 cortical regions, large to very large effect sizes (Cohen's d ≥ 1) and statistically significant differences in qMRI values between patients and controls were revealed compared with only four regions when using more standard MR measures, namely, volume and cortical thickness. Moreover, a decrease in all susceptibility contrasts ( χ , χ + , χ - ) and R 2 * values indicates that the role of cortical demyelination might outweigh inflammatory processes in the form of iron accumulation in cortical MS pathology, and might also indicate iron loss. A significant association between susceptibility contrasts as well as R 2 * of the caudal middle frontal gyrus and disease duration was found (adjusted R2 : 0.602, p = 0.0011). Quantitative MRI parameters might be more sensitive towards regional cortical pathology compared with the use of conventional markers only and therefore may play a role in early detection of tissue damage in MS in the future.

QSM Throughout the Body.

Magnetic materials in tissue, such as iron, calcium, or collagen, can be studied using quantitative susceptibility mapping (QSM). To date, QSM has been overwhelmingly applied in the brain, but is increasingly utilized outside the brain. QSM relies on the effect of tissue magnetic susceptibility sources on the MR signal phase obtained with gradient echo sequence. However, in the body, the chemical shift of fat present within the region of interest contributes to the MR signal phase as well. Therefore, correcting for the chemical shift effect by means of water-fat separation is essential for body QSM. By employing techniques to compensate for cardiac and respiratory motion artifacts, body QSM has been applied to study liver iron and fibrosis, heart chamber blood and placenta oxygenation, myocardial hemorrhage, atherosclerotic plaque, cartilage, bone, prostate, breast calcification, and kidney stone.

A novel gradient echo data based vein segmentation algorithm and its application for the detection of regional cerebral differences in venous susceptibility.

Accurate segmentation of cerebral venous vasculature from gradient echo data is of central importance in several areas of neuroimaging such as for the susceptibility-based assessment of brain oxygenation or planning of electrode placement in deep brain stimulation. In this study, a vein segmentation algorithm for single- and multi-echo gradient echo data is proposed. First, susceptibility maps, true susceptibility-weighted images, and, in the multi-echo case, R2* maps were generated from the gradient echo data. These maps were filtered with an inverted Hamming filter to suppress background contrast as well as artifacts from field inhomogeneities at the brain boundaries. A shearlet-based scale-wise representation was generated to calculate a vesselness function and to generate segmentations based on local thresholding. The accuracy of the proposed algorithm was evaluated for different echo times and image resolutions using a manually generated reference segmentation and two vein segmentation algorithms (Frangi vesselness-based, recursive vesselness filter) as a reference with the Dice and Cohen's coefficients as well as the modified Hausdorff distance. The Frangi-based and recursive vesselness filter methods were significantly outperformed with regard to all error metrics. Applying the algorithm, susceptibility differences likely related to differences in blood oxygenation between superficial and deep venous territories could be demonstrated.

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.

Toward quantitative neuroimaging biomarkers for Friedreich's ataxia at 7 Tesla: Susceptibility mapping, diffusion imaging, R2 and R1 relaxometry.

Friedreich's ataxia (FRDA) is a rare genetic disorder leading to degenerative processes. So far, no effective treatment has been found. Therefore, it is important to assist the development of medication with imaging biomarkers reflecting disease status and progress. Ten FRDA patients (mean age 37 ± 14 years; four female) and 10 age- and sex-matched controls were included. Acquisition of magnetic resonance imaging (MRI) data for quantitative susceptibility mapping, R1 , R2 relaxometry and diffusion imaging was performed at 7 Tesla. Results of volume of interest (VOI)-based analyses of the quantitative data were compared with a voxel-based morphometry (VBM) evaluation. Differences between patients and controls were assessed using the analysis of covariance (ANCOVA; p < 0.01) with age and sex as covariates, effect size of group differences, and correlations with disease characteristics with Spearman correlation coefficient. For the VBM analysis, a statistical threshold of 0.001 for uncorrected and 0.05 for corrected p-values was used. Statistically significant differences between FRDA patients and controls were found in five out of twelve investigated structures, and statistically significant correlations with disease characteristics were revealed. Moreover, VBM revealed significant white matter atrophy within regions of the brainstem, and the cerebellum. These regions overlapped partially with brain regions for which significant differences between healthy controls and patients were found in the VOI-based quantitative MRI evaluation. It was shown that two independent analyses provided overlapping results. Moreover, positive results on correlations with disease characteristics were found, indicating that these quantitative MRI parameters could provide more detailed information and assist the search for effective treatments.

Mapping the human brainstem: Brain nuclei and fiber tracts at 3 T and 7 T.

Structural high-resolution imaging of the brainstem can be of high importance in clinical practice. However, ultra-high field magnetic resonance imaging (MRI) is still restricted in use due to limited availability. Therefore, quantitative MRI techniques (quantitative susceptibility mapping [QSM], relaxation measurements [ R2* , R1 ], diffusion tensor imaging [DTI]) and T2 - and proton density (PD)-weighted imaging in the human brainstem at 3 T and 7 T are compared. Five healthy volunteers (mean age: 21.5 ± 1.9 years) were measured at 3 T and 7 T using multi-echo gradient echo sequences for susceptibility mapping and R2* relaxometry, magnetization-prepared 2 rapid acquisition gradient echo sequences for R1 relaxometry, turbo-spin echo sequences for PD- and T2 -weighted imaging and readout-segmented echo planar sequences for DTI. Susceptibility maps were computed using Laplacian-based phase unwrapping, V-SHARP for background field removal and the streaking artifact reduction for QSM algorithm for dipole inversion. Contrast-to-noise ratios (CNRs) were determined at 3 T and 7 T in ten volumes of interest (VOIs). Data acquired at 7 T showed higher CNR. However, in four VOIs, lower CNR was observed for R2* at 7 T. QSM was shown to be the contrast with which the highest number of structures could be identified. The depiction of very fine tracts such as the medial longitudinal fasciculus throughout the brainstem was only possible in susceptibility maps acquired at 7 T. DTI effectively showed the main tracts (crus cerebri, transverse pontine fibers, corticospinal tract, middle and superior cerebellar peduncle, pontocerebellar tract, and pyramid) at both field strengths. Assessing the brainstem with quantitative MRI methods such as QSM, R2* , as well as PD- and T2 -weighted imaging with great detail, is also possible at 3 T, especially when using susceptibility mapping calculated from a gradient echo sequence with a wide range of echo times from 10.5 to 52.5 ms. However, tracing smallest structures strongly benefits from imaging at ultra-high field.

Rapid and accurate dictionary-based T2 mapping from multi-echo turbo spin echo data at 7 Tesla.

BACKGROUND: Using lower refocusing flip angles in multi-echo turbo spin echo (ME-TSE) sequences at ultra-high magnetic field leads to non-monoexponential signal decay and overestimation of T2 values due to stimulated and secondary echoes. PURPOSE: To investigate the feasibility of a fast and accurate reconstruction of quantitative T2 values using an ME-TSE sequence with reduced refocusing flip angles at 7 Tesla, a dictionary-based reconstruction method was developed and is presented in this work. STUDY TYPE: Prospective. SUBJECTS: Phantom measurements with relaxation phantom, four healthy volunteers. FIELD STRENGTH/SEQUENCE: 7 Tesla MRI, multi-echo turbo spin echo (ME-TSE), spin echo (SE), and B1 mapping. ASSESSMENT: Based on Bloch simulations and the extended phase graph model, signal decay curves were calculated to account for nonrectangular slice profile, B1 inhomogeneity, and reduced refocusing flip angles and stored in a dictionary. Data obtained with an ME-TSE sequence at 7 Tesla were matched to this dictionary to obtain T2 values. To compare the proposed method to reference T2 values, a spin echo sequence with different echo times was used. STATISTICAL TESTS: Welch's t-test was used to compare T2 values in phantom measurements. RESULTS: T2 values obtained with the proposed ME-TSE method coincided with the T2 values from the spin echo experiment in phantom measurements (P = 0.89 for 120° flip angle, P = 0.75 for 180° flip angle). Only for very low B1 transmit fields, a slight overestimation of T2 values was observed. In vivo measurements showed lower T2 values in gray matter (55 ± 2 millisecond) and white matter (39 ± 5 millisecond) compared with literature values of 3 Tesla data. DATA CONCLUSIONS: The proposed dictionary-based ME-TSE approach provided accurate T2 values in short measurement time at 7 Tesla with low specific absorption rate burden due to the reduction of refocusing flip angles. Therefore, it can provide new opportunities in clinical high-field MRI to further improve radiographic diagnosis by using quantitative imaging. LEVEL OF EVIDENCE: 1 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2019;49:1253-1262.

Quantitative susceptibility mapping and 23 Na imaging-based in vitro characterization of blood clotting kinetics.

Blood clotting is a fundamental biochemical process in post-hemorrhagic hemostasis. Although the varying appearance of coagulating blood in T1 - and T2 -weighted images is widely used to qualitatively determine bleeding age, the technique permits only a rough discrimination of coagulation stages, and it remains difficult to distinguish acute and chronic hemorrhagic stages because of low T1 - and T2 -weighted signal intensities in both instances. To investigate new biomedical parameters for magnetic resonance imaging-based characterization of blood clotting kinetics, sodium imaging and quantitative susceptibility mapping (QSM) were compared with conventional T1 - and T2 -weighted imaging, as well as with biochemical hemolysis parameters. For this purpose, a blood-filled spherical agar phantom was investigated daily for 14 days, as well as after 24 days at 7 T after initial preparation with fresh blood. T1 - and T2 -weighted sequences, a three-dimensional (3D) gradient echo sequence and a density-adapted 3D radial projection reconstruction pulse sequence for 23 Na imaging were applied. For hemolysis estimations, free hemoglobin and free potassium concentrations were measured photometrically and with the direct ion-selective electrode method, respectively, in separate heparinized whole-blood samples along the same timeline. Initial mean susceptibility was low (0.154 ± 0.020 ppm) and increased steadily during the course of coagulation to reach up to 0.570 ± 0.165 ppm. The highest total sodium (NaT) values (1.02 ± 0.06 arbitrary units) in the clot were observed initially, dropped to 0.69 ± 0.13 arbitrary units after one day and increased again to initial values. Compartmentalized sodium (NaS) showed a similar signal evolution, and the NaS/NaT ratio steadily increased over clot evolution. QSM depicts clot evolution in vitro as a process associated with hemoglobin accumulation and transformation, and enables the differentiation of the acute and chronic coagulation stages. Sodium imaging visualizes clotting independent of susceptibility and seems to correspond to clot integrity. A combination of QSM and sodium imaging may enhance the characterization of hemorrhage.

Suitable reference tissues for quantitative susceptibility mapping of the brain.

PURPOSE: Since quantitative susceptibility mapping (QSM) quantifies magnetic susceptibility relative to a reference value, a suitable reference tissue has to be available to compare different subjects and stages of disease. METHODS: To find such a suitable reference tissue for QSM of the brain, melanoma patients with and without brain metastases were measured. Twelve reference regions were chosen and assessed for stability of susceptibility values with respect to multiple intra-individual and inter-individual measurements, age, and stage of disease. RESULTS: Cerebrospinal fluid (CSF), the internal capsule and one region in the splenium of the corpus callosum are the regions with the smallest standard deviations of the mean susceptibility value. The mean susceptibility is 0.010 ± 0.014 ppm for CSF in the atrium of the lateral ventricles (csfpost ), -0.060 ± 0.019 ppm for the posterior limb of the internal capsule (ci2), and -0.008 ± 0.019 ppm for the splenium of the corpus callosum. csfpost and ci2 show nearly no dependence on age or stage of disease, whereas some other regions, e.g., the red nucleus, show moderate dependence on age or disease. CONCLUSION: The internal capsule and CSF appear to be the most suitable reference regions for QSM of the brain in the melanoma patients studied. Both showed virtually no dependence on age or disease and small variations among patients. Magn Reson Med 78:204-214, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Potential of quantitative susceptibility mapping for detection of prostatic calcifications.

PURPOSE: To evaluate whether quantitative susceptibility (QSM) may be used as an alternative to computed tomography (CT) to detect calcification in prostate cancer patients. MATERIALS AND METHODS: Susceptibility map calculation was performed using 3D gradient echo magnetic resonance imaging (MRI) data from 26 patients measured at 3T who previously received a planning CT of the prostate. Phase images were unwrapped using Laplacian-based phase unwrapping, the background field was removed with the V-SHARP method, and susceptibility maps were calculated with the iLSQR method. Two blinded readers were asked to identify peri- and intraprostatic calcifications. RESULTS: Average mean and minimum susceptibility values (referenced to iliopsoas muscle) of calcifications were -0.249 ± 0.179 ppm and -0.551 ± 0.323 ppm, and average mean and maximum intensities in CT images were 319 ± 164 HU and 679 ± 392 HU. Twenty-one and 17 out of 22 prostatic calcifications were identified using susceptibility maps and magnitude images, respectively, as well as more than half of periprostatic phleboliths depicted by CT. Calcifications in the prostate and its periphery were quantitatively differentiable from noncalcified prostate tissue in CT (mean values for calcifications / for noncalcified tissue: 71 to 649 / -1 to 83 HU) and in QSM (mean values for calcifications / for noncalcified tissue: -0.641 to 0.063 / -0.046 to 0.181 ppm). Moreover, there was a significant correlation between susceptibility values and CT image intensities for calcifications (P < 0.004). CONCLUSION: Prostatic calcifications could be well identified with QSM. Susceptibility maps can be easily obtained from clinical prostate MR protocols that include a 3D gradient echo sequence, rendering it a promising technique for detection and quantification of intraprostatic calcifications. LEVEL OF EVIDENCE: 1 J. Magn. Reson. Imaging 2017;45:889-898.

On contrast mechanisms in p-space imaging.

PURPOSE: The effects of microstructure on susceptibility mapping have recently received increasing attention. In this work, the capability of p-space imaging to resolve subvoxel structure and susceptibility is assessed. METHODS: In a simulation study, the p-space contrast of axon bundles comprised of hollow cylinders is investigated. Various axon and susceptibility distributions are considered and compared with the corresponding case of a voxel with homogeneous substructure of the order of the voxel size. RESULTS: MR signal behavior for p-space imaging of a voxel containing axon bundles and a voxel containing homogeneous substructure are nearly identical. CONCLUSION: p-Space imaging resolves subvoxel structure of the order of the voxel size. Due to dephasing effects, strong susceptibility variations alter p-space contrast. However, p-space contrast is not directly linked to the susceptibility of the axon compartment. Magn Reson Med 75:2526-2533, 2016. © 2015 Wiley Periodicals, Inc.

Diagnostic utility of 7T neuromelanin imaging of the substantia nigra in Parkinson's disease.

Parkinson's disease (PD) is a prevalent neurodegenerative disorder that presents a diagnostic challenge due to symptom overlap with other disorders. Neuromelanin (NM) imaging is a promising biomarker for PD, but adoption has been limited, in part due to subpar performance at standard MRI field strengths. We aimed to evaluate the diagnostic utility of ultra-high field 7T NM-sensitive imaging in the diagnosis of PD versus controls and essential tremor (ET), as well as NM differences among PD subtypes. A retrospective case-control study was conducted including PD patients, ET patients, and controls. 7T NM-sensitive 3D-GRE was acquired, and substantia nigra pars compacta (SNpc) volumes, contrast ratios, and asymmetry indices were calculated. Statistical analyses, including general linear models and ROC curves, were employed. Twenty-one PD patients, 13 ET patients, and 18 controls were assessed. PD patients exhibited significantly lower SNpc volumes compared to non-PD subjects. SNpc total volume showed 100% sensitivity and 96.8% specificity (AUC = 0.998) for differentiating PD from non-PD and 100% sensitivity and 95.2% specificity (AUC = 0.996) in differentiating PD from ET. Contrast ratio was not significantly different between PD and non-PD groups (p = 0.07). There was also significantly higher asymmetry index in SNpc volume in PD compared to non-PD cohorts (p < 0.001). NM signal loss in PD predominantly involved the inferior, posterior, and lateral aspects of SNpc. Akinetic-rigid subtype showed more significant NM signal loss compared to tremor dominant subtype (p < 0.001). 7T NM imaging demonstrates potential as a diagnostic tool for PD, including potential distinction between subtypes, allowing improved understanding of disease progression and subtype-related characteristics.

7 T Lesion-Attenuated Magnetization-Prepared Gradient Echo Acquisition for Detection of Posterior Fossa Demyelinating Lesions in Multiple Sclerosis.

OBJECTIVES: Detection of infratentorial demyelinating lesions in multiple sclerosis (MS) presents a challenge in magnetic resonance imaging (MRI), a difficulty that is further heightened in 7 T MRI. This study aimed to assess the efficacy of a novel MRI approach, lesion-attenuated magnetization-prepared gradient echo acquisition (LAMA), for detecting demyelinating lesions within the posterior fossa and upper cervical spine on 7 T MRI and contrast its performance with conventional double-inversion recovery (DIR) and T2-weighted turbo spin echo sequences. MATERIALS AND METHODS: We conducted a retrospective cross-sectional study in 42 patients with a confirmed diagnosis of MS. All patients had 7 T MRI that incorporated LAMA, 3D DIR, and 2D T2-weighted turbo spin echo sequences. Three readers assessed lesion count in the brainstem, cerebellum, and upper cervical spinal cord using both DIR and T2-weighted images in one session. In a separate session, LAMA was analyzed alone. Contrast-to-noise ratio was also compared between LAMA and the conventional sequences. Lesion counts between methods were assessed using nonparametric Wilcoxon signed rank test. Interrater agreement in lesion detection was estimated by intraclass correlation coefficients. RESULTS: LAMA identified a significantly greater number of lesions than DIR + T2 (mean 6.4 vs 3.0; P < 0.001). LAMA also exhibited better interrater agreement (intraclass correlation coefficient [95% confidence interval], 0.75 [0.41-0.88] vs 0.61 [0.35-0.78]). The contrast-to-noise ratio for LAMA (3.7 ± 0.9) significantly exceeded that of DIR (1.94 ± 0.7) and T2 (1.2 ± 0.7) (all P 's < 0.001). In cases with no lesions detected using DIR + T2, at least 1 lesion was identified in 83.3% with LAMA. Across all analyzed brain regions, LAMA consistently detected more lesions than DIR + T2. CONCLUSIONS: LAMA significantly improves the detection of infratentorial demyelinating lesions in MS patients compared with traditional methods. Integrating LAMA with standard magnetization-prepared 2 rapid acquisition gradient echo acquisition provides a valuable tool for accurately characterizing the extent of MS disease.

Enhancing outcomes in deep brain stimulation: a comparative study of direct targeting using 7T versus 3T MRI.

OBJECTIVE: The aim of this study was to compare outcomes of direct targeting in deep brain stimulation (DBS) for essential tremor using 7T MRI versus 3T MRI. The authors hypothesized that 7T MRI direct targeting would be noninferior to 3T MRI in early tremor outcomes. METHODS: A retrospective study was conducted on patients undergoing unilateral thalamic DBS for essential tremor between 2021 and 2023. Two matched cohorts were assessed, one using 7T MRI and the other using 3T MRI for surgical planning. The primary endpoint was the percentage improvement in the Fahn-Tolosa-Marin Tremor Rating Scale (TRS) scores. Additionally, the authors assessed optimized programming settings and variance in electrode position on postoperative imaging. Demographic and clinical data were compared using the nonparametric Mann-Whitney U-test. The squared Euclidean distance of each contact from the group mean centroid was calculated and averaged across the entire cohort to provide the variance (i.e., the mean squared distance) of electrode contact position. RESULTS: A total of 34 patients were analyzed, with 17 in each cohort. There were no significant differences in demographic information or mean surgical dates between the groups. There were no differences in intraoperative target repositioning or adverse events. The 7T group had a significantly greater TRS improvement than the 3T group (64.9% ± 11.4% vs 50.9% ± 16.4%, p = 0.004). Patients in the 7T cohort also had a lower mean stimulation current compared with those in the 3T cohort (2.0 ± 0.8 mA vs 2.7 ± 0.9 mA, p = 0.01). Image evaluation revealed that although the mean electrode position was comparable between 7T and 3T, the 7T electrode positioning was more clustered, indicating a lower variance in the final electrode location. The mean Euclidean distance between the individual electrode tips and the group centroid was significantly less at 7T than at 3T (1.82 ± 0.68 mm vs 2.75 ± 0.81 mm, p = 0.001). CONCLUSIONS: Despite concerns for increased artifacts and distortions at 7T, the authors show that these effects can be mitigated with an appropriate workflow, leading to improved surgical outcomes with direct targeting using 7T MRI. Their results suggest similar accuracy but greater precision in targeting with 7T MRI compared with 3T MRI, resulting in lower stimulation currents and improved tremor reduction. Future studies are needed to assess outcomes related to 7T MRI in targeting other subcortical structures.

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.

Central Vein Sign in Multiple Sclerosis: A Comparison Study of the Diagnostic Performance of 3T versus 7T MRI.

BACKGROUND AND PURPOSE: An early and accurate diagnosis of multiple sclerosis remains challenging in clinical neurology. Established diagnostic methods have less than desirable sensitivity and specificity. An accurate, noninvasive diagnostic test for MS could have a major impact on diagnostic criteria. We compared the frequency of detection of the central vein sign (CVS) in white matter lesions of MS and controls on 7T T2*-weighted and SWI to 3T SWI. Additionally, we assessed the diagnostic performance of 7T T2*, 7T SWI, and 3T SWI for MS. MATERIALS AND METHODS: A retrospective case-control study was performed of patients with MS having both 7T MRI and 3T MRI. A control group of patients without MS was selected. Diagnosis of MS was established by board-certified neurologists with fellowship training in autoimmune neurology in line with the 2017 McDonald criteria. Percentage of lesions with a CVS was blindly measured for each technique. Diagnostic performance was computed by sensitivity, specificity, and positive and negative likelihood ratios (LRs). RESULTS: Sixty-one patients with MS (903 lesions) and 39 controls (1088 lesions) were included. 7T T2* showed significantly more CVS (87%) than both 7T SWI (73%) and 3T SWI (31%) (all P < .001). CVS was identified in the control group in ≤6% of lesions on all sequences. Using a threshold of >40% of lesions with CVS on 7T T2* and >15% on 7T SWI, both sequences had an accuracy = 100%, sensitivity = 100%, specificity = 100%, infinite positive LR, and zero negative LR. Using an optimal threshold of >12%, 3T SWI had an accuracy = 96.0%, sensitivity = 93.4%, specificity = 100%, infinite positive LR, and negative LR = 0.066. CONCLUSIONS: 7T MRI had 100% sensitivity and specificity for the diagnosis of MS and is superior to 3T. Future revisions to MS diagnostic criteria may consider recommendations for 7T MRI and inclusion of CVS as a biomarker.

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 ΔR2∗ 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 Dm=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.

A novel phantom with dia- and paramagnetic substructure for quantitative susceptibility mapping and relaxometry.

PURPOSE: A phantom is presented in this study that allows for an experimental evaluation of QSM reconstruction algorithms. The phantom contains susceptibility producing particles with dia- and paramagnetic properties embedded in an MRI visible medium and is suitable to assess the performance of algorithms that attempt to separate isotropic dia- and paramagnetic susceptibility at the sub-voxel level. METHODS: The phantom was built from calcium carbonate (diamagnetic) and tungsten carbide particles (paramagnetic) embedded in gelatin and surrounded by agarose gel. Different mass fractions and mixing ratios of both susceptibility sources were used. Gradient echo data were acquired at 1.5 T, 3 T and 7 T. Susceptibility maps were calculated using the MEDI toolbox and relaxation rates ΔR2∗ were determined using exponential fitting. RESULTS: Relaxation rates as well as susceptibility values generally coincide with the theoretical values for particles fulfilling the assumptions of the the static dephasing regime with stronger deviations for relaxation rates at higher field strength and for high susceptibility values. MRI raw data are available for free academic use as supplementary material. CONCLUSIONS: In this study, a susceptibility phantom is presented that might find its application in the development and quantitative validation of current and future QSM reconstruction algorithms which aim to separate the influence of isotropic dia- and paramagnetic substructure in quantitative susceptibility mapping.

Multiparametric MRI for Characterization of the Basal Ganglia and the Midbrain.

Objectives To characterize subcortical nuclei by multi-parametric quantitative magnetic resonance imaging. Materials and Methods: The following quantitative multiparametric MR data of five healthy volunteers were acquired on a 7T MRI system: 3D gradient echo (GRE) data for the calculation of quantitative susceptibility maps (QSM), GRE sequences with and without off-resonant magnetic transfer pulse for magnetization transfer ratio (MTR) calculation, a magnetization-prepared 2 rapid acquisition gradient echo sequence for T1 mapping, and (after a coil change) a density-adapted 3D radial pulse sequence for 23Na imaging. First, all data were co-registered to the GRE data, volumes of interest (VOIs) for 21 subcortical structures were drawn manually for each volunteer, and a combined voxel-wise analysis of the four MR contrasts (QSM, MTR, T1, 23Na) in each structure was conducted to assess the quantitative, MR value-based differentiability of structures. Second, a machine learning algorithm based on random forests was trained to automatically classify the groups of multi-parametric voxel values from each VOI according to their association to one of the 21 subcortical structures. Results The analysis of the integrated multimodal visualization of quantitative MR values in each structure yielded a successful classification among nuclei of the ascending reticular activation system (ARAS), the limbic system and the extrapyramidal system, while classification among (epi-)thalamic nuclei was less successful. The machine learning-based approach facilitated quantitative MR value-based structure classification especially in the group of extrapyramidal nuclei and reached an overall accuracy of 85% regarding all selected nuclei. Conclusion Multimodal quantitative MR enabled excellent differentiation of a wide spectrum of subcortical nuclei with reasonable accuracy and may thus enable sensitive detection of disease and nucleus-specific MR-based contrast alterations in the future.