SEGUE: A Speedy rEgion-Growing Algorithm for Unwrapping Estimated Phase

SIPAS: A comprehensive susceptibility imaging process and analysis studio

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

How to accurately quantify brain magnetic susceptibility in the context of Parkinson's disease: Validation on phantoms and healthy volunteers at 1.5 and 3 T

Currently, brain iron content represents a new neuromarker for understanding the physiopathological mechanisms leading to Parkinson's disease (PD). In vivo quantification of biological iron is possible by reconstructing magnetic susceptibility maps obtained using quantitative susceptibility mapping (QSM). Applying QSM is challenging, as up to now, no standardization of acquisition protocols and phase image processing has emerged from referenced studies. Our objectives were to compare the accuracy and the sensitivity of 10 QSM pipelines built from algorithms from the literature, applied on phantoms data and on brain data. Two phantoms, with known magnetic susceptibility ranges, were created from several solutions of gadolinium chelate. Twenty healthy volunteers from two age groups were included. Phantoms and brain data were acquired at 1.5 and 3 T, respectively. Susceptibility-weighted images were obtained using a 3D multigradient-recalled-echo sequence. For brain data, 3D anatomical T1- and T2-weighted images were also acquired to segment the deep gray nuclei of interest. Concerning in vitro data, the linear dependence of magnetic susceptibility versus gadolinium concentration and deviations from the theoretically expected values were calculated. For brain data, the accuracy and sensitivity of the QSM pipelines were evaluated in comparison with results from the literature and regarding the expected magnetic susceptibility increase with age, respectively. A nonparametric Mann-Whitney U-test was used to compare the magnetic susceptibility quantification in deep gray nuclei between the two age groups. Our methodology enabled quantifying magnetic susceptibility in human brain and the results were consistent with those from the literature. Statistically significant differences were obtained between the two age groups in all cerebral regions of interest. Our results show the importance of optimizing QSM pipelines according to the application and the targeted magnetic susceptibility range, to achieve accurate quantification. We were able to define the optimal QSM pipeline for future applications on patients with PD.

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

INTRODUCTION

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

METHODS

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

RESULTS

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

CONCLUSION

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

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.

Single-shot frequency offset measurement with HASTE using the selective parity approach

Measurements of frequency offset are commonly required in MRI. The standard method measures the signal phase as a function of evolution time. Here we use a single shot turbo-spin-echo acquisition method to measure frequency offset at a single evolution time. After excitation the transverse magnetisation evolves during the evolution time, and is then repeatedly refocused. The phase is conjugated between alternate echoes. Using partial parallel acquisition techniques we obtain separate odd- and even- echo images. An iterative procedure ensures self-consistency between them. The difference in phase between the two images yields frequency offset maps. The technique was implemented at 3 Tesla and tested on a healthy human volunteer for a range of evolution times between 6 and 42 ms. A standard method using a similar readout train and multiple evolution times was used as a gold-standard measure. In a statistical comparison with the gold standard no evidence for bias or offset was found. There was no systematic variation in precision or accuracy as a function of evolution time. We conclude that the presented approach represents a viable method for the rapid generation of frequency offset maps with a high image quality and minimal distortion.

A predictor-corrector phase unwrapping algorithm for temporally undersampled gradient-echo MRI

PURPOSE

To develop a method for unwrapping temporally undersampled and nonlinear gradient recalled echo (GRE) phase.

THEORY AND METHODS

Temporal unwrapping is performed as a sequential one step prediction of the echo phase, followed by a correction to the nearest integer wrap-count. A spatio-temporal extension of the 1D predictor corrector unwrapping (PCU) algorithm improves the prediction accuracy, and thereby maintains spatial continuity. The proposed method is evaluated using numerical phantom, physical phantom, and in vivo brain data at both 3 T and 9.4 T. The unwrapping performance is compared with the state-of-the-art temporal and spatial unwrapping algorithms, and the spatio-temporal iterative virtual-echo based Nyquist sampled (iVENyS) algorithm.

RESULTS

Simulation results showed significant reduction in unwrapping errors at higher echoes compared with the state-of-the-art algorithms. Similar to the iVENyS algorithm, the PCU algorithm was able to generate spatially smooth phase images for in vivo data acquired at 3 T and 9.4 T, bypassing the use of additional spatial unwrapping step. A key advantage over iVENyS algorithm is the superior performance of PCU algorithm at higher echoes.

CONCLUSION

PCU algorithm serves as a robust phase unwrapping method for temporally undersampled and nonlinear GRE phase, particularly in the presence of high field gradients.

Neuroimaging of Parkinson's disease by quantitative susceptibility mapping

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

Incorporating the effect of white matter microstructure in the estimation of magnetic susceptibility in ex vivo mouse brain

PURPOSE

To extend quantitative susceptibility mapping to account for microstructure of white matter (WM) and demonstrate its effect on ex vivo mouse brain at 16.4T.

THEORY AND METHODS

Previous studies have shown that the MRI measured Larmor frequency also depends on local magnetic microstructure at the mesoscopic scale. Here, we include effects from WM microstructure using our previous results for the mesoscopic Larmor frequencyΩ ‾ Meso $$ {\overline{\Omega}}^{\mathrm{Meso}} $$ of cylinders with arbitrary orientations. We scrutinize the validity of our model and QSM in a digital brain phantom includingΩ ‾ Meso $$ {\overline{\Omega}}^{\mathrm{Meso}} $$ from a WM susceptibility tensor and biologically stored iron with scalar susceptibility. We also apply susceptibility tensor imaging to the phantom and investigate how the fitted tensors are biased fromΩ ‾ Meso $$ {\overline{\Omega}}^{\mathrm{Meso}} $$ . Last, we demonstrate how to combine multi-gradient echo and diffusion MRI images of ex vivo mouse brains acquired at 16.4T to estimate an apparent scalar susceptibility without sample rotations.

RESULTS

Our new model improves susceptibility estimation compared to QSM for the brain phantom. Applying susceptibility tensor imaging to the phantom withΩ ‾ Meso $$ {\overline{\Omega}}^{\mathrm{Meso}} $$ from WM axons with scalar susceptibility produces a highly anisotropic susceptibility tensor that mimics results from previous susceptibility tensor imaging studies. For the ex vivo mouse brain we find theΩ ‾ Meso $$ {\overline{\Omega}}^{\mathrm{Meso}} $$ due to WM microstructure to be substantial, changing susceptibility in WM up to 25% root-mean-squared-difference.

CONCLUSION

Ω ‾ Meso $$ {\overline{\Omega}}^{\mathrm{Meso}} $$ impacts susceptibility estimates and biases susceptibility tensor imaging fitting substantially. Hence, it should not be neglected when imaging structurally anisotropic tissue such as brain WM.

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

Background

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

Methods

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

Results

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

Conclusion

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

Regional cortical thinning, demyelination and iron loss in cerebral small vessel disease

The link between white matter hyperintensities (WMH) and cortical thinning is thought to be an important pathway by which WMH contributes to cognitive deficits in cerebral small vessel disease (SVD). However, the mechanism behind this association and the underlying tissue composition abnormalities are unclear. The objective of this study is to determine the association between WMH and cortical thickness, and the in vivo tissue composition abnormalities in the WMH-connected cortical regions. In this cross-sectional study, we included 213 participants with SVD who underwent standardized protocol including multimodal neuroimaging scans and cognitive assessment (i.e. processing speed, executive function and memory). We identified the cortex connected to WMH using probabilistic tractography starting from the WMH and defined the WMH-connected regions at three connectivity levels (low, medium and high connectivity level). We calculated the cortical thickness, myelin and iron of the cortex based on T1-weighted, quantitative R1, R2* and susceptibility maps. We used diffusion-weighted imaging to estimate the mean diffusivity of the connecting white matter tracts. We found that cortical thickness, R1, R2* and susceptibility values in the WMH-connected regions were significantly lower than in the WMH-unconnected regions (all Pcorrected < 0.001). Linear regression analyses showed that higher mean diffusivity of the connecting white matter tracts were related to lower thickness (β = -0.30, Pcorrected < 0.001), lower R1 (β = -0.26, Pcorrected = 0.001), lower R2* (β = -0.32, Pcorrected < 0.001) and lower susceptibility values (β = -0.39, Pcorrected < 0.001) of WMH-connected cortical regions at high connectivity level. In addition, lower scores on processing speed were significantly related to lower cortical thickness (β = 0.20, Pcorrected = 0.030), lower R1 values (β = 0.20, Pcorrected = 0.006), lower R2* values (β = 0.29, Pcorrected = 0.006) and lower susceptibility values (β = 0.19, Pcorrected = 0.024) of the WMH-connected regions at high connectivity level, independent of WMH volumes and the cortical measures of WMH-unconnected regions. Together, our study demonstrated that the microstructural integrity of white matter tracts passing through WMH is related to the regional cortical abnormalities as measured by thickness, R1, R2* and susceptibility values in the connected cortical regions. These findings are indicative of cortical thinning, demyelination and iron loss in the cortex, which is most likely through the disruption of the connecting white matter tracts and may contribute to processing speed impairment in SVD, a key clinical feature of SVD. These findings may have implications for finding intervention targets for the treatment of cognitive impairment in SVD by preventing secondary degeneration.

Reducing the ambiguity of field inhomogeneity and chemical shift effect for fat-water separation by field factor

PURPOSE

To reduce the ambiguity between chemical shift and field inhomogeneity with flexible TE combinations by introducing a variable (field factor).

THEORY AND METHODS

The ambiguity between chemical shift and field inhomogeneity can be eliminated directly from the multiple in-phase images acquired at different TEs; however, it is only applicable to few echo combinations. In this study, we accommodated such an implementation in flexible TE combinations by introducing a new variable (field factor). The effects of the chemical shift were removed from the field inhomogeneity in the candidate solutions, thus reducing the ambiguity problem. To validate this concept, multi-echo MRI data acquired from various anatomies with different imaging parameters were tested. The derived fat and water images were compared with those of the state-of-the-art fat-water separation algorithms.

RESULTS

Robust fat-water separation was achieved with the accurate solution of field inhomogeneity, and no apparent fat-water swap was observed. In addition to the good performance, the proposed method is applicable to various fat-water separation applications, including different sequence types and flexible TE choices.

CONCLUSION

We propose an algorithm to reduce the ambiguity of chemical shift and field inhomogeneity and achieved robust fat-water separation in various applications.

Investigating the effect of oblique image acquisition on the accuracy of QSM and a robust tilt correction method

PURPOSE

Quantitative susceptibility mapping (QSM) is used increasingly for clinical research where oblique image acquisition is commonplace, but its effects on QSM accuracy are not well understood.

THEORY AND METHODS

The QSM processing pipeline involves defining the unit magnetic dipole kernel, which requires knowledge of the direction of the main magnetic fieldB ^ 0 $$ {\hat{\boldsymbol{B}}}_{\mathbf{0}} $$ with respect to the acquired image volume axes. The direction ofB ^ 0 $$ {\hat{\boldsymbol{B}}}_{\mathbf{0}} $$ is dependent on the axis and angle of rotation in oblique acquisition. Using both a numerical brain phantom and in vivo acquisitions in 5 healthy volunteers, we analyzed the effects of oblique acquisition on magnetic susceptibility maps. We compared three tilt-correction schemes at each step in the QSM pipeline: phase unwrapping, background field removal and susceptibility calculation, using the RMS error and QSM-tuned structural similarity index.

RESULTS

Rotation of wrapped phase images gave severe artifacts. Background field removal with projection onto dipole fields gave the most accurate susceptibilities when the field map was first rotated into alignment withB ^ 0 $$ {\hat{\boldsymbol{B}}}_{\mathbf{0}} $$ . Laplacian boundary value and variable-kernel sophisticated harmonic artifact reduction for phase data background field removal methods gave accurate results without tilt correction. For susceptibility calculation, thresholded k-space division, iterative Tikhonov regularization, and weighted linear total variation regularization, all performed most accurately when local field maps were rotated into alignment withB ^ 0 $$ {\hat{\boldsymbol{B}}}_{\mathbf{0}} $$ before susceptibility calculation.

CONCLUSION

For accurate QSM, oblique acquisition must be taken into account. Rotation of images into alignment withB ^ 0 $$ {\hat{\boldsymbol{B}}}_{\mathbf{0}} $$ should be carried out after phase unwrapping and before background-field removal. We provide open-source tilt-correction code to incorporate easily into existing pipelines: https://github.com/o-snow/QSM_TiltCorrection.git.

Simulated and experimental phantom data for multi-center quality assurance of quantitative susceptibility maps at 3 T, 7 T and 9.4 T

PURPOSE

To develop methods for quality assurance of quantitative susceptibility mapping (QSM) using MRI at different magnetic field strengths, and scanners, using different MR-sequence protocols, and post-processing pipelines.

METHODS

We built a custom phantom based on iron in two forms: homogeneous susceptibility ('free iron') and with fine-scaled variations in susceptibility ('clustered iron') at different iron concentrations. The phantom was measured at 3.0 T (two scanners), 7.0 T and 9.4 T using multi-echo, gradient echo acquisition sequences. A digital phantom analogue to the iron-phantom, tailored to obtain similar results as in experimentation was developed, with similar geometry and susceptibility values. Morphology enabled dipole inversion was applied to the phase images to obtain QSM for experimental and simulated data using the MEDI + 0 approach for background regularization.

RESULTS

Across all scanners, QSM-values showed a linear increase with iron concentrations. The QSM-relaxivity was 0.231 ± 0.047 ppm/mM for free and 0.054 ± 0.013 ppm/mM for clustered iron, with adjusted determination coefficients (DoC) ≥ 0.87. Similarly, the simulations yielded linear increases (DoC ≥ 0.99). In both the experimental and digital phantoms, the estimated molar susceptibility was lower with clustered iron, because clustering led to highly localized field effects.

CONCLUSION

Our iron phantom can be used to evaluate the capability of QSM to detect local variations in susceptibility across different field strengths, when using different MR-sequence protocols. The devised simulation method captures the effect of iron clustering in QSM as seen experimentally and could be used in the future to optimize QSM processing pipelines and achieve higher accuracy for local field effects, as also seen in Alzheimer's beta-amyloid plaques.

Increased iron deposition in nucleus accumbens associated with disease progression and chronicity in migraine

BACKGROUND

Migraine is one of the world's most prevalent and disabling diseases. Despite huge advances in neuroimaging research, more valuable neuroimaging markers are still urgently needed to provide important insights into the brain mechanisms that underlie migraine symptoms. We therefore aim to investigate the regional iron deposition in subcortical nuclei of migraineurs as compared to controls and its association with migraine-related pathophysiological assessments.

METHODS

A total of 200 migraineurs (56 chronic migraine [CM], 144 episodic migraine [EM]) and 41 matched controls were recruited. All subjects underwent MRI and clinical variables including frequency/duration of migraine, intensity of migraine, 6-item Headache Impact Test (HIT-6), Migraine Disability Assessment (MIDAS), and Pittsburgh Sleep Quality Index (PSQI) were recorded. Quantitative susceptibility mapping was employed to quantify the regional iron content in subcortical regions. Associations between clinical variables and regional iron deposition were studied as well.

RESULTS

Increased iron deposition in the putamen, caudate, and nucleus accumbens (NAC) was observed in migraineurs more than controls. Meanwhile, patients with CM had a significantly higher volume of iron deposits compared to EM in multiple subcortical nuclei, especially in NAC. Volume of iron in NAC can be used to distinguish patients with CM from EM with a sensitivity of 85.45% and specificity of 71.53%. As the most valuable neuroimaging markers in all of the subcortical nuclei, higher iron deposition in NAC was significantly associated with disease progression, and higher HIT-6, MIDAS, and PSQI.

CONCLUSIONS

These findings provide evidence that iron deposition in NAC may be a biomarker for migraine chronicity and migraine-related dysfunctions, thus may help to understand the underlying vascular and neural mechanisms of migraine.

TRIAL REGISTRATION

ClinicalTrials.gov, number NCT04939922.

Repeatability of brain phase-based magnetic resonance electric properties tomography methods and effect of compressed SENSE and RF shimming

Magnetic resonance electrical properties tomography (MREPT) is an emerging imaging modality to noninvasively measure tissue conductivity and permittivity. Implementation of MREPT in the clinic requires repeatable measurements at a short scan time and an appropriate protocol. The aim of this study was to investigate the repeatability of conductivity measurements using phase-based MREPT and the effects of compressed SENSE (CS), and RF shimming on the precision of conductivity measurements. Conductivity measurements using turbo spin echo (TSE) and three-dimensional balanced fast field echo (bFFE) with CS factors were repeatable. Conductivity measurement using bFFE phase showed smaller mean and variance that those measured by TSE. The conductivity measurements using bFFE showed minimal deviation with CS factors up to 8, with deviation increasing at CS factors > 8. Subcortical structures produced less consistent measurements than cortical parcellations at higher CS factors. RF shimming using full slice coverage 2D dual refocusing echo acquisition mode (DREAM) and full coverage 3D dual TR approaches further improved measurement precision. BFFE is a more optimal sequence than TSE for phase-based MREPT in brain. Depending on the area of the brain being measured, the scan can be safely accelerated with compressed SENSE without sacrifice of precision, offering the potential to employ MREPT in clinical research and applications. RF shimming with better field mapping further improves precision of the conductivity measures.

Depth-wise profiles of iron and myelin in the cortex and white matter using χ-separation: A preliminary study

The in-vivo profiling of iron and myelin across cortical depths and underlying white matter has important implications for advancing knowledge about their roles in brain development and degeneration. Here, we utilize χ-separation, a recently-proposed advanced susceptibility mapping that creates positive (χ_pos) and negative (χ_neg) susceptibility maps, to generate the depth-wise profiles of χ_pos and χ_neg as surrogate biomarkers for iron and myelin, respectively. Two regional sulcal fundi of precentral and middle frontal areas are profiled and compared with findings from previous studies. The results show that the χ_pos profiles peak at superificial white matter (SWM), which is an area beneath cortical gray matter known to have the highest accumulation of iron within the cortex and white matter. On the other hand, the χ_neg profiles increase in SWM toward deeper white matter. These characteristics in the two profiles are in agreement with histological findings of iron and myelin. Furthermore, the χ_neg profiles report regional differences that agree with well-known distributions of myelin concentration. When the two profiles are compared with those of QSM and R₂*, different shapes and peak locations are observed. This preliminary study offers an insight into one of the possible applications of χ-separation for exploring microstructural information of the human brain, as well as clinical applications in monitoring changes of iron and myelin in related diseases.

A method for measuring B0 field inhomogeneity using quantitative double-echo in steady-state

PURPOSE

To develop and validate a method forB 0 $$ {B}_0 $$ mapping for knee imaging using the quantitative Double-Echo in Steady-State (qDESS) exploiting the phase difference (Δ θ $$ \Delta \theta $$ ) between the two echoes acquired. Contrary to a two-gradient-echo (2-GRE) method,Δ θ $$ \Delta \theta $$ depends only on the first echo time.

METHODS

Bloch simulations were applied to investigate robustness to noise of the proposed methodology and all imaging studies were validated with phantoms and in vivo simultaneous bilateral knee acquisitions. Two phantoms and five healthy subjects were scanned using qDESS, water saturation shift referencing (WASSR), and multi-GRE sequences.Δ B 0 $$ \Delta {B}_0 $$ maps were calculated with the qDESS and the 2-GRE methods and compared against those obtained with WASSR. The comparison was quantitatively assessed exploiting pixel-wise difference maps, Bland-Altman (BA) analysis, and Lin's concordance coefficient (ρ c $$ {\rho}_c $$ ). For in vivo subjects, the comparison was assessed in cartilage using average values in six subregions.

RESULTS

The proposed method for measuringΔ B 0 $$ \Delta {B}_0 $$ inhomogeneities from a qDESS acquisition providedΔ B 0 $$ \Delta {B}_0 $$ maps that were in good agreement with those obtained using WASSR.Δ B 0 $$ \Delta {B}_0 $$ ρ c $$ {\rho}_c $$ values were≥ $$ \ge $$ 0.98 and 0.90 in phantoms and in vivo, respectively. The agreement between qDESS and WASSR was comparable to that of a 2-GRE method.

CONCLUSION

The proposed method may allow B0 correction for qDESST 2 $$ {T}_2 $$ mapping using an inherently co-registeredΔ B 0 $$ \Delta {B}_0 $$ map without requiring an additional B0 measurement sequence. More generally, the method may help shorten knee imaging protocols that require an auxiliaryΔ B 0 $$ \Delta {B}_0 $$ map by exploiting a qDESS acquisition that also providesT 2 $$ {T}_2 $$ measurements and high-quality morphological imaging.

Unsupervised physiological noise correction of functional magnetic resonance imaging data using phase and magnitude information (PREPAIR)

Of the sources of noise affecting blood oxygen level-dependent functional magnetic resonance imaging (fMRI), respiration and cardiac fluctuations are responsible for the largest part of the variance, particularly at high and ultrahigh field. Existing approaches to removing physiological noise either use external recordings, which can be unwieldy and unreliable, or attempt to identify physiological noise from the magnitude fMRI data. Data-driven approaches are limited by sensitivity, temporal aliasing, and the need for user interaction. In the light of the sensitivity of the phase of the MR signal to local changes in the field stemming from physiological processes, we have developed an unsupervised physiological noise correction method using the information carried in the phase and the magnitude of echo-planar imaging data. Our technique, Physiological Regressor Estimation from Phase and mAgnItude, sub-tR (PREPAIR) derives time series signals sampled at the slice TR from both phase and magnitude images. It allows physiological noise to be captured without aliasing, and efficiently removes other sources of signal fluctuations not related to physiology, prior to regressor estimation. We demonstrate that the physiological signal time courses identified with PREPAIR agree well with those from external devices and retrieve challenging cardiac dynamics. The removal of physiological noise was as effective as that achieved with the most used approach based on external recordings, RETROICOR. In comparison with widely used recording-free physiological noise correction tools-PESTICA and FIX, both performed in unsupervised mode-PREPAIR removed significantly more respiratory and cardiac noise than PESTICA, and achieved a larger increase in temporal signal-to-noise-ratio at both 3 and 7 T.

Phase-based masking for quantitative susceptibility mapping of the human brain at 9.4T

PURPOSE

To develop improved tissue masks for QSM.

METHODS

Masks including voxels at the brain surface were automatically generated from the magnitude alone (MM) or combined with test functions from the first (PG) or second (PB) derivative of the sign of the wrapped phase. Phase images at 3T and 9.4T were simulated at different TEs and used to generate a mask, PItoh , with between-voxel phase differences less than π. MM, PG, and PB were compared with PItoh . QSM were generated from 3D multi-echo gradient-echo data acquired at 9.4T (21 subjects aged: 20-56y), and from the QSM2016 challenge 3T data using different masks, unwrapping, background removal, and dipole inversion algorithms. QSM contrast was quantified using age-based iron concentrations.

RESULTS

Close to air cavities, phase wraps became denser with increasing field and echo time, yielding increased values of the test functions. Compared with PItoh , PB had the highest Dice coefficient, while PG had the lowest and MM the highest percentage of voxels outside PItoh. Artifacts observed in QSM at 9.4T with MM were mitigated by stronger background filters but yielded a reduced QSM contrast. With PB, QSM contrast was greater and artifacts diminished. Similar results were obtained with challenge data, evidencing larger effects of mask close to air cavities.

CONCLUSION

Automatic, phase-based masking founded on the second derivative of the sign of the wrapped phase, including cortical voxels at the brain surface, was able to mitigate artifacts and restore QSM contrast across cortical and subcortical brain regions.

Phase unwrapping algorithm for images with local high-density noise

Due to undersampling and the local phase with local high-density noise, it is still difficult to develop a robust phase unwrapping algorithm. In order to resolve this issue, here, we propose what we believe to be a novel multiple path-following phase unwrapping (MPIPU) algorithm based on the shearing interference principle to recover the undersampling phase (non-noise). By calculating the unwrapping coefficient k, the phase iteration filling algorithm based on least-squares is developed for the high-density noise region in order to reconstruct the three-dimensional surface topography of interferometric synthetic aperture radar (InSAR) data. The proposed algorithm takes advantage of the MPIPU's ability to fill in the missing phase with fitting data and can successfully suppress the error transfer caused by the blocky noise phase iteration process. Several experiments are conducted using both simulated and actual InSAR image data. The experimental findings show that the proposed method can achieve robust phase unwrapping performance on a phase of local high-density noise.

Quantitative susceptibility mapping (QSM) and R2* of silent cerebral infarcts in sickle cell anemia

Silent cerebral infarction (SCI) is the most commonly reported radiological abnormality in patients with sickle cell anemia (SCA) and is associated with future clinical stroke risk. To date, there have been few histological and quantitative MRI studies of SCI and multiple radiological definitions exist. As a result, the tissue characteristics and composition of SCI remain elusive. The objective of this work was therefore to investigate the composition of segmented SCI lesions using quantitative MRI for R2 * and quantitative magnetic susceptibility mapping (QSM). 211 SCI lesions were segmented from 32 participants with SCA and 6 controls. SCI were segmented according to two definitions (FLAIR+/-T1w-based threshold) using a semi-automated pipeline. Magnetic susceptibility (χ) and R2 * maps were calculated from a multi-echo gradient echo sequence and mean SCI values were compared to an equivalent region of interest in normal appearing white matter (NAWM). SCI χ and R2 * were investigated as a function of SCI definition, patient demographics, anatomical location, and cognition. Compared to NAWM, SCI were significantly less diamagnetic (χ = -0.0067 ppm vs. -0.0153 ppm, p < 0.001) and had significantly lower R2 * (16.7 s-1 vs. 19.2 s-1, p < 0.001). SCI definition had a significant effect on the mean SCI χ and R2 * , with lesions becoming significantly less diamagnetic and having significantly lower R2 * after the application of a more stringent T1w-based threshold. SCI-NAWM R2 * decrease was significantly greater in patients with SCA compared with controls (-2.84 s-1 vs. -0.64 s-1, p < 0.0001). No significant association was observed between mean SCI-NAWM χ or R2* differences and subject age, lesion anatomical location, or cognition. The increased χ and decreased R2 * in SCI relative to NAWM observed in both patients and controls is indicative of lower myelin or increased water content within the segmented lesions. The significant SCI-NAWM R2 * differences observed between SCI in patients with SCA and controls suggests there may be differences in tissue composition relative to NAWM in SCI in the two populations. Quantitative MRI techniques such as QSM and R2 * mapping can be used to enhance our understanding of the pathophysiology and composition of SCI in patients with SCA as well as controls.

A Comparison of MRI Quantitative Susceptibility Mapping and TRUST-Based Measures of Brain Venous Oxygen Saturation in Sickle Cell Anaemia

In recent years, interest has grown in the potential for magnetic resonance imaging (MRI) measures of venous oxygen saturation (Yv) to improve neurological risk prediction. T2-relaxation-under-spin-tagging (TRUST) is an MRI technique which has revealed changes in Yv in patients with sickle cell anemia (SCA). However, prior studies comparing Yv in patients with SCA relative to healthy controls have reported opposing results depending on whether the calibration model, developed to convert blood T2 to Yv, is based on healthy human hemoglobin (HbA), bovine hemoglobin (HbBV) or sickle hemoglobin (HbS). MRI Quantitative Susceptibility Mapping (QSM) is an alternative technique that may hold promise for estimating Yv in SCA as blood magnetic susceptibility is linearly dependent upon Yv, and no significant difference has been found between the magnetic susceptibility of HbA and HbS. Therefore, the aim of this study was to compare estimates of Yv using QSM and TRUST with five published calibration models in healthy controls and patients with SCA. 17 patients with SCA and 13 healthy controls underwent MRI. Susceptibility maps were calculated from a multi-parametric mapping acquisition and Yv was calculated from the mean susceptibility in a region of interest in the superior sagittal sinus. TRUST estimates of T2, within a similar but much smaller region, were converted to Yv using five different calibration models. Correlation and Bland-Altman analyses were performed to compare estimates of Yv between TRUST and QSM methods. For each method, t-tests were also used to explore group-wise differences between patients with SCA and healthy controls. In healthy controls, significant correlations were observed between QSM and TRUST measures of Yv, while in SCA, there were no such correlations. The magnitude and direction of group-wise differences in Yv varied with method. The TRUST-HbBV and QSM methods suggested decreased Yv in SCA relative to healthy controls, while the TRUST-HbS (p &lt; 0.01) and TRUST-HbA models suggested increased Yv in SCA as in previous studies. Further validation of all MRI measures of Yv, relative to ground truth measures such as O15 PET and jugular vein catheterization, is required in SCA before QSM or TRUST methods can be considered for neurological risk prediction.

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

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

Shim Coils Tailored for Correcting B0 Inhomogeneity in the Human Brain (SCOTCH): Design Methodology and 48-Channel Prototype Assessment in 7-Tesla MRI

Increased static field inhomogeneities are a burden for human brain MRI at Ultra-High-Field. In particular they cause enhanced Echo-Planar image distortions and signal losses due to magnetic susceptibility gradients at air-tissue interfaces in the subject's head. In the past decade, Multi-Coil Arrays (MCA) have been proposed to shim the field in the brain better than the 2nd or 3rd order Spherical Harmonic (SH) coils usually offered by MRI manufacturers. Here we present a novel MCA, named SCOTCH, optimized for whole brain shimming. Based on a cylindrical structure, it features several layers of small coils whose shape, size and location are found from a principal component analysis of ideal stream functions computed from an internal 100-brain fieldmap database. From an Open-Access external database of 126 brains, our SCOTCH implementation is shown to be equivalent to a partial 7th-order SH system with unlimited power, outperforming all known existing MCA prototypes. This result is further confirmed by a low-cost  30-cm diameter SCOTCH prototype built with 48 coils on 3 layers, and tested on 7 volunteers at 7T with a parallel-transmit RF coil made to be inserted in SCOTCH. Echo-Planar images of the subject brains before and after SCOTCH shimming show large signal recoveries, especially in the prefrontal cortex.

Hybrid data fidelity term approach for quantitative susceptibility mapping

Purpose

Methods

Results

Conclusions

Imaging white matter microstructure with gradient-echo phase imaging: Is ex vivo imaging with formalin-fixed tissue a good approximation of the in vivo brain?

Purpose

Ex vivo imaging is a commonly used approach to investigate the biophysical mechanism of orientation‐dependent signal phase evolution in white matter. Yet, how phase measurements are influenced by the structural alteration in the tissue after formalin fixation is not fully understood. Here, we study the effects on magnetic susceptibility, microstructural compartmentalization, and chemical exchange measurement with a postmortem formalin‐fixed whole‐brain human tissue.

Methods

A formalin‐fixed, postmortem human brain specimen was scanned with multiple orientations to the main magnetic field direction for robust bulk magnetic susceptibility measurement with conventional quantitative susceptibility imaging models. White matter samples were subsequently excised from the whole‐brain specimen and scanned in multiple rotations on an MRI scanner to measure the anisotropic magnetic susceptibility and microstructure‐related contributions in the signal phase and to validate the findings of the whole‐brain data.

Results

The bulk isotropic magnetic susceptibility of ex vivo whole‐brain imaging is comparable to in vivo imaging, with noticeable enhanced nonsusceptibility contributions. The excised specimen experiment reveals that anisotropic magnetic susceptibility and compartmentalization phase effect were considerably reduced in the formalin‐fixed white matter specimens.

Conclusions

Formalin‐fixed postmortem white matter exhibits comparable isotropic magnetic susceptibility to previous in vivo imaging findings. However, the measured phase and magnitude data of the fixed white matter tissue shows a significantly weaker orientation dependency and compartmentalization effect. Alternatives to formalin fixation are needed to better reproduce the in vivo microstructural effects in postmortem samples.

Magnetic Resonance Imaging Versus Computed Tomography for Three-Dimensional Bone Imaging of Musculoskeletal Pathologies: A Review

Magnetic resonance imaging (MRI) is increasingly utilized as a radiation‐free alternative to computed tomography (CT) for the diagnosis and treatment planning of musculoskeletal pathologies. MR imaging of hard tissues such as cortical bone remains challenging due to their low proton density and short transverse relaxation times, rendering bone tissues as nonspecific low signal structures on MR images obtained from most sequences. Developments in MR image acquisition and post‐processing have opened the path for enhanced MR‐based bone visualization aiming to provide a CT‐like contrast and, as such, ease clinical interpretation. The purpose of this review is to provide an overview of studies comparing MR and CT imaging for diagnostic and treatment planning purposes in orthopedic care, with a special focus on selective bone visualization, bone segmentation, and three‐dimensional (3D) modeling. This review discusses conventional gradient‐echo derived techniques as well as dedicated short echo time acquisition techniques and post‐processing techniques, including the generation of synthetic CT, in the context of 3D and specific bone visualization. Based on the reviewed literature, it may be concluded that the recent developments in MRI‐based bone visualization are promising. MRI alone provides valuable information on both bone and soft tissues for a broad range of applications including diagnostics, 3D modeling, and treatment planning in multiple anatomical regions, including the skull, spine, shoulder, pelvis, and long bones.

Improved susceptibility weighted imaging at ultra-high field using bipolar multi-echo acquisition and optimized image processing: CLEAR-SWI

Purpose

Susceptibility Weighted Imaging (SWI) has become established in the clinical investigation of stroke, microbleeds, tumor vascularization, calcification and iron deposition, but suffers from a number of shortcomings and artefacts. The goal of this study was to reduce the sensitivity of SWI to strong B₁ and B₀ inhomogeneities at ultra-high field to generate homogeneous images with increased contrast and free of common artefacts. All steps in SWI processing have been addressed −coil combination, phase unwrapping, image combination over echoes, phase filtering and homogeneity correction −and applied to an efficient bipolar multi-echo acquisition to substantially improve the quality of SWI.

Principal results

Our findings regarding the optimal individual processing steps lead us to propose a Contrast-weighted, Laplace-unwrapped, bipolar multi-Echo, ASPIRE-combined, homogeneous, improved Resolution SWI, or CLEAR-SWI. CLEAR-SWI was compared to two other multi-echo SWI methods and standard, single-echo SWI with the same acquisition time at 7 T in 10 healthy volunteers and with single-echo SWI in 13 patients with brain tumors. CLEAR-SWI had improved contrast-to-noise and homogeneity, reduced signal dropout and was not compromised by the artefacts which affected standard SWI in 10 out of 13 cases close to tumors (as assessed by expert raters), as well as generating T2* maps and phase images which can be used for Quantitative Susceptibility Mapping. In a comparison with other multi-echo SWI methods, CLEAR-SWI had the fewest artefacts, highest SNR and generally higher contrast-to-noise.

Major conclusions

CLEAR-SWI eliminates the artefacts common in standard, single-echo SWI, reduces signal dropouts and improves image homogeneity and contrast-to-noise. Applied clinically, in a study of brain tumor patients, CLEAR-SWI was free of the artefacts which affected standard, single-echo SWI.

The PHU-NET: A robust phase unwrapping method for MRI based on deep learning

PURPOSE

This work was aimed at designing a deep-learning-based approach for MR image phase unwrapping to improve the robustness and efficiency of traditional methods.

METHODS

A deep learning network called PHU-NET was designed for MR phase unwrapping. In this network, a novel training data generation method was proposed to simulate the wrapping patterns in MR phase images. The wrapping boundary and wrapping counts were explicitly estimated and used for network training. The proposed method was quantitatively evaluated and compared to other methods using a number of simulated datasets with varying signal-to-noise ratio (SNR) and MR phase images from various parts of the human body.

RESULTS

The results showed that our method performed better in the simulated data even under an extremely low SNR. The proposed method had less residual wrapping in the images from various parts of human body and worked well in the presence of severe anatomical discontinuity. Our method was also advantageous in terms of computational efficiency compared to the traditional methods.

CONCLUSION

This work proposed a robust and computationally efficient MR phase unwrapping method based on a deep learning network, which has promising performance in applications using MR phase information.

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.

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.

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.

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.

Phase unwrapping with a rapid opensource minimum spanning tree algorithm (ROMEO)

PURPOSE

To develop a rapid and accurate MRI phase-unwrapping technique for challenging phase topographies encountered at high magnetic fields, around metal implants, or postoperative cavities, which is sufficiently fast to be applied to large-group studies including Quantitative Susceptibility Mapping and functional MRI (with phase-based distortion correction).

METHODS

The proposed path-following phase-unwrapping algorithm, ROMEO, estimates the coherence of the signal both in space-using MRI magnitude and phase information-and over time, assuming approximately linear temporal phase evolution. This information is combined to form a quality map that guides the unwrapping along a 3D path through the object using a computationally efficient minimum spanning tree algorithm. ROMEO was tested against the two most commonly used exact phase-unwrapping methods, PRELUDE and BEST PATH, in simulated topographies and at several field strengths: in 3T and 7T in vivo human head images and 9.4T ex vivo rat head images.

RESULTS

ROMEO was more reliable than PRELUDE and BEST PATH, yielding unwrapping results with excellent temporal stability for multi-echo or multi-time-point data. It does not require image masking and delivers results within seconds, even in large, highly wrapped multi-echo data sets (eg, 9 seconds for a 7T head data set with 31 echoes and a 208 × 208 × 96 matrix size).

CONCLUSION

Overall, ROMEO was both faster and more accurate than PRELUDE and BEST PATH, delivering exact results within seconds, which is well below typical image acquisition times, enabling potential on-console application.

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.

Multi-parametric artificial neural network fitting of phase-cycled balanced steady-state free precession data

PURPOSE

Standard relaxation time quantification using phase-cycled balanced steady-state free precession (bSSFP), eg, motion-insensitive rapid configuration relaxometry (MIRACLE), is subject to a considerable underestimation of tissue T₁ and T₂ due to asymmetric intra-voxel frequency distributions. In this work, an artificial neural network (ANN) fitting approach is proposed to simultaneously extract accurate reference relaxation times (T₁ , T₂ ) and robust field map estimates ( B 1 + , ΔB₀ ) from the bSSFP profile.

METHODS

Whole-brain bSSFP data acquired at 3T were used for the training of a feedforward ANN with N = 12, 6, and 4 phase-cycles. The magnitude and phase of the Fourier transformed complex bSSFP frequency response served as input and the multi-parametric reference set [T₁ , T₂ , B 1 + , ∆B₀ ] as target. The ANN predicted relaxation times were validated against the target and MIRACLE.

RESULTS

The ANN prediction of T₁ and T₂ for trained and untrained data agreed well with the reference, even for only four acquired phase-cycles. In contrast, relaxometry based on 4-point MIRACLE was prone to severe off-resonance-related artifacts. ANN predicted B 1 + and ∆B₀ maps showed the expected spatial inhomogeneity patterns in high agreement with the reference measurements for 12-point, 6-point, and 4-point bSSFP phase-cycling schemes.

CONCLUSION

ANNs show promise to provide accurate brain tissue T₁ and T₂ values as well as reliable field map estimates. Moreover, the bSSFP acquisition can be accelerated by reducing the number of phase-cycles while still delivering robust T₁ , T₂ , B 1 + , and ∆B₀ estimates.

Experimental setup for bulk susceptibility effect-minimized, multi-orientation MRI of ex vivo tissue samples

PURPOSE

Many conventional ex vivo magnetic resonance imaging (MRI) setups utilize cylindrical or other nonspherical tissue containers which can cause static field (B₀ ) inhomogeneity affecting the accuracy of the measurements in an orientation-dependent manner. In this work we demonstrate an experimental method to obtain MRI of ex vivo tissue samples held in a spherical container in order to minimize bulk susceptibility-induced B₀ inhomogeneity in arbitrary orientations.

METHODS

B₀ inhomogeneity caused by tissue-air susceptibility mismatch can be theoretically eliminated if the surface of susceptibility discontinuity is spherical. This situation can be approximated by putting a tissue sample in a spherical shell filled with materials with tissue-like magnetic susceptibility. We achieved this on an intact monkey brain by (a) holding the brain with a three-dimensional (3D)-printed holder with tissue-like (within 0.5 ppm) susceptibility, and (b) enclosing the brain and the holder in an acrylic spherical shell filled with diamagnetic liquid. Furthermore, the sphere and the radio-frequency coil for MRI were mounted on a 3D-printed frame designed to reduce B₀ inhomogeneity contributions. The sphere could be rotated freely without disturbing the RF coil to facilitate multi-orientation imaging. We verified our setup by mapping B₀ in the monkey brain at 13 different orientations in a human 7T scanner, and measuring orientation-dependent R 2 ∗ relaxation rates in the white and gray matters of the brain. The results were then compared with a setup where the brain was held inside a cylindrical container.

RESULTS

In all orientations, the B₀ standard deviation in the brain in the spherical setup (converted to Larmor frequency offset) was less than about 10 Hz. This corresponds to two-sigma deviation of B₀ of <0.07 ppm. The B₀ gradient was <9 Hz/mm in 95 % of the brain voxels in all orientations. In high-resolution imaging with e.g. voxel size <0.4 mm, this corresponds to voxel line broadening of <4 Hz (0.013 ppm). R 2 ∗ in the corpus callosum showed distinctly different orientation dependence compared to the gray matter. The B₀ uniformity and R 2 ∗ reliability were much reduced in the cylindrical container setup.

CONCLUSIONS

We have demonstrated an experimental method to effectively minimize bulk susceptibility-induced B₀ perturbation in multi-orientation ex vivo MRI. The method promises to benefit a range of tissue orientation-dependent MR property studies.