Micro-compartment specific T2* relaxation in the brain

Quantitative BOLD (qBOLD) imaging of oxygen metabolism and blood oxygenation in the human body: A scoping review

PURPOSE

There are many approaches to the quantitative BOLD (qBOLD) technique described in the literature, differing in pulse sequences, MRI parameters and data processing. Thus, in this review, we summarized the acquisition methods, approaches used for oxygenation quantification and clinical populations investigated.

METHODS

Three databases were systematically searched (Medline, Embase, and Web of Science) for published research that used qBOLD methods for quantification of oxygen metabolism. Data extraction and synthesis were performed by one author and reviewed by a second author.

RESULTS

A total of 93 relevant papers were identified. Acquisition strategies were summarized, and oxygenation parameters were found to have been investigated in many pathologies such as steno-occlusive diseases, stroke, glioma, and multiple sclerosis disease.

CONCLUSION

A summary of qBOLD approaches for oxygenation measurements and applications could help researchers to identify good practice and provide objective information to inform the development of future consensus recommendations.

The Larmor frequency shift of a white matter magnetic microstructure model with multiple sources

Magnetic susceptibility imaging may provide valuable information about chemical composition and microstructural organization of tissue. However, its estimation from the MRI signal phase is particularly difficult as it is sensitive to magnetic tissue properties ranging from the molecular to the macroscopic scale. The MRI Larmor frequency shift measured in white matter (WM) tissue depends on the myelinated axons and other magnetizable sources such as iron-filled ferritin. We have previously derived the Larmor frequency shift arising from a dense medium of cylinders with scalar susceptibility and arbitrary orientation dispersion. Here, we extend our model to include microscopic WM susceptibility anisotropy as well as spherical inclusions with scalar susceptibility to represent subcellular structures, biologically stored iron, and so forth. We validate our analytical results with computer simulations and investigate the feasibility of estimating susceptibility using simple iterative linear least squares without regularization or preconditioning. This is done in a digital brain phantom synthesized from diffusion MRI measurements of an ex vivo mouse brain at ultra-high field.

Harnessing myelin water fraction as an imaging biomarker of human cerebral aging, neurodegenerative diseases, and risk factors influencing myelination: A review

Myelin water fraction (MWF) imaging has emerged as a promising magnetic resonance imaging (MRI) biomarker for investigating brain function and composition. This comprehensive review synthesizes the current state of knowledge on MWF as a biomarker of human cerebral aging, neurodegenerative diseases, and risk factors influencing myelination. The databases used include Web of Science, Scopus, Science Direct, and PubMed. We begin with a brief discussion of the theoretical foundations of MWF imaging, including its basis in MR physics and the mathematical modeling underlying its calculation, with an overview of the most adopted MRI methods of MWF imaging. Next, we delve into the clinical and research applications that have been explored to date, highlighting its advantages and limitations. Finally, we explore the potential of MWF to serve as a predictive biomarker for neurological disorders and identify future research directions for optimizing MWF imaging protocols and interpreting MWF in various contexts. By harnessing the power of MWF imaging, we may gain new insights into brain health and disease across the human lifespan, ultimately informing novel diagnostic and therapeutic strategies.

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

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

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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

Quantum dipole interactions and transient hydrogen bond orientation order in cells, cellular membranes and myelin sheath: Implications for MRI signal relaxation, anisotropy, and T1 magnetic field dependence

PURPOSE

Despite significant impact on the study of human brain, MRI lacks a theory of signal formation that integrates quantum interactions involving proton dipoles (a primary MRI signal source) with brain intricate cellular environment. The purpose of the present study is developing such a theory.

METHODS

We introduce the Transient Hydrogen Bond (THB) model, where THB-mediated quantum dipole interactions between water and protons of hydrophilic heads of amphipathic biomolecules forming cells, cellular membranes and myelin sheath serve as a major source of MR signal relaxation.

RESULTS

The THB theory predicts the existence of a hydrogen-bond-driven structural order of dipole-dipole connections within THBs as a primary factor for the anisotropy observed in MRI signal relaxation. We have also demonstrated that the conventional Lorentzian spectral density function decreases too fast at high frequencies to adequately capture the field dependence of brain MRI signal relaxation. To bridge this gap, we introduced a stretched spectral density function that surpasses the limitations of Lorentzian dispersion. In human brain, our findings reveal that at any time point only about 4% to 7% of water protons are engaged in quantum encounters within THBs. These ultra-short (2 to 3 ns), but frequent quantum spin exchanges lead to gradual recovery of magnetization toward thermodynamic equilibrium, that is, relaxation of MRI signal.

CONCLUSION

By incorporating quantum proton interactions involved in brain imaging, the THB approach introduces new insights on the complex relationship between brain tissue cellular structure and MRI measurements, thus offering a promising new tool for better understanding of brain microstructure in health and disease.

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.

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.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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

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.

Empirical relationship between TEM-derived myelin volume fraction and MRI-R2 values in aging ex vivo rat corpus callosum

Ex vivo ratiometric measurements of short- and long-T2 components using the multiple spin echo sequence of MRI are often employed to evaluate alterations in myelin content in the white matter (WM) of the brain. However, the relationship between absolute MRI-T2 values (long-T2 component) and myelin volumetric information in aged ex vivo rodent WM appears to be influenced by factors such as animal species, field strength, and fixation durations/washing. Here, multiple spin echo sequence-based MRI-R2 (the reciprocal of T2) values were measured in the corpus callosum (CC) region in the post-mortem rat brains (n = 9) of different age groups with common fixation techniques without washing at 7 T. Transmission electron microscopy (TEM)-based quantification of myelin volume fraction (MVF) and corresponding Monte-Carlo simulation to estimate relaxation rates (R2,IE) due to diffusion in the presence of inhomogeneous magnetic field perturbation in intra- and extra-cellular (IE) spaces were respectively performed. To determine whether the short-T2 components originating from myelin water were mixed with long-T2 components from IE water or were undetectable, the MVF values obtained from TEM results were respectively compared with MRI-R2 and R2,IE values. A significant correlation (Pearson's correlation coefficient r = 0.8763; p < 0.01) of average MRI-R2 and MVF values was observed. Estimated R2,IE values from Monte-Carlo simulations in IE water signals were also positively correlated (r = 0.8281; p < 0.01) with MVF values. However, the magnitudes of R2,IE values were much smaller than those observed for MRI-R2 values, indicating that changes in R2 related MVF are likely dominated by myelin water components. Such comparisons between independent parameters from MRI, TEM, and simulations support the suggestion that myelin water signals were indistinguishably mixed to exhibit mono-exponential T2 relaxation, and multiple spin echo sequence-based MRI-R2 values in aging ex vivo rat CC without prolonged washing still reflect the volumetric information of myelin, likely due to enhanced water exchange across the myelin.

The impact of head orientation with respect to B0 on diffusion tensor MRI measures

Diffusion tensor MRI (DT-MRI) remains the most commonly used approach to characterise white matter (WM) anisotropy. However, DT estimates may be affected by tissue orientation w.r.t. B → 0 due to local gradients and intrinsic T 2 orientation dependence induced by the microstructure. This work aimed to investigate whether and how diffusion tensor MRI-derived measures depend on the orientation of the head with respect to the static magnetic field, B → 0 . By simulating WM as two compartments, we demonstrated that compartmental T 2 anisotropy can induce the dependence of diffusion tensor measures on the angle between WM fibres and the magnetic field. In in vivo experiments, reduced radial diffusivity and increased axial diffusivity were observed in white matter fibres perpendicular to B → 0 compared to those parallel to B → 0 . Fractional anisotropy varied by up to 20 % as a function of the angle between WM fibres and the orientation of the main magnetic field. To conclude, fibre orientation w.r.t. B → 0 is responsible for up to 7 % variance in diffusion tensor measures across the whole brain white matter from all subjects and head orientations. Fibre orientation w.r.t. B → 0 may introduce additional variance in clinical research studies using diffusion tensor imaging, particularly when it is difficult to control for (e.g., fetal or neonatal imaging, or when the trajectories of fibres change due to, e.g., space occupying lesions).

Improved test-retest reliability of R2* and susceptibility quantification using multi-shot multi-echo 3D EPI

This study aimed to evaluate the potential of 3D echo-planar imaging (EPI) for improving the reliability of T2*-weighted (T2*w) data and quantification of R2* decay rate and susceptibility (χ) compared to conventional gradient echo (GRE)-based acquisition. Eight healthy subjects in a wide age range were recruited. Each subject received repeated scans for both GRE and EPI acquisitions with an isotropic 1 mm resolution at 3 T. Maps of R2* and χ were quantified and compared using their inter-scan difference to evaluate the test-retest reliability. Inter-protocol differences of R2* and χ between GRE and EPI were also measured voxel by voxel and in selected ROIs to test the consistency between the two acquisition methods. The quantifications of R2* and χ using EPI protocols showed increased test-retest reliability with higher EPI factors up to 5 as performed in the experiment and were consistent with those based on GRE. This result suggested multi-shot multi-echo 3D EPI can be a useful alternative acquisition method for T2*w MRI and quantification of R2* and χ with reduced scan time, improved test-retest reliability and similar accuracy compared to commonly used 3D GRE.

Recommended Implementation of Quantitative Susceptibility Mapping for Clinical Research in The Brain: A Consensus of the ISMRM Electro-Magnetic Tissue Properties Study Group

This article provides recommendations for implementing quantitative susceptibility mapping (QSM) for clinical brain research. It is a consensus of the ISMRM Electro-Magnetic Tissue Properties Study Group. While QSM technical development continues to advance rapidly, the current QSM methods have been demonstrated to be repeatable and reproducible for generating quantitative tissue magnetic susceptibility maps in the brain. However, the many QSM approaches available give rise to the need in the neuroimaging community for guidelines on implementation. This article describes relevant considerations and provides specific implementation recommendations for all steps in QSM data acquisition, processing, analysis, and presentation in scientific publications. We recommend that data be acquired using a monopolar 3D multi-echo GRE sequence, that phase images be saved and exported in DICOM format and unwrapped using an exact unwrapping approach. Multi-echo images should be combined before background removal, and a brain mask created using a brain extraction tool with the incorporation of phase-quality-based masking. Background fields should be removed within the brain mask using a technique based on SHARP or PDF, and the optimization approach to dipole inversion should be employed with a sparsity-based regularization. Susceptibility values should be measured relative to a specified reference, including the common reference region of whole brain as a region of interest in the analysis, and QSM results should be reported with - as a minimum - the acquisition and processing specifications listed in the last section of the article. These recommendations should facilitate clinical QSM research and lead to increased harmonization in data acquisition, analysis, and reporting.

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

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

Effect of motion, cortical orientation and spatial resolution on quantitative imaging of cortical R2* and magnetic susceptibility at 0.3 mm in-plane resolution at 7 T

MR images of the effective relaxation rate R2* and magnetic susceptibility χ derived from multi-echo T2*-weighted (T2*w) MRI can provide insight into iron and myelin distributions in the brain, with the potential of providing biomarkers for neurological disorders. Quantification of R2* and χ at submillimeter resolution in the cortex in vivo has been difficult because of challenges such as head motion, limited signal to noise ratio, long scan time, and motion related magnetic field fluctuations. This work aimed to improve the robustness for quantifying intracortical R2* and χ and analyze the effects from motion, spatial resolution, and cortical orientation. T2*w data was acquired with a spatial resolution of 0.3 × 0.3 × 0.4 mm3 at 7 T and downsampled to various lower resolutions. A combined correction for motion and B0 changes was deployed using volumetric navigators. Such correction improved the T2*w image quality rated by experienced image readers and test-retest reliability of R2* and χ quantification with reduced median inter-scan differences up to 10 s-1 and 5 ppb, respectively. R2* and χ near the line of Gennari, a cortical layer high in iron and myelin, were as much as 10 s-1 and 10 ppb higher than the region at adjacent cortical depth. In addition, a significant effect due to the cortical orientation relative to the static field (B0) was observed in χ with a peak-to-peak amplitude of about 17 ppb. In retrospectively downsampled data, the capability to distinguish different cortical depth regions based on R2* or χ contrast remained up to isotropic 0.5 mm resolution. This study highlights the unique characteristics of R2* and χ along the cortical depth at submillimeter resolution and the need for motion and B0 corrections for their robust quantification in vivo.

Larmor frequency shift from magnetized cylinders with arbitrary orientation distribution

The magnetic susceptibility of tissue can provide valuable information about its chemical composition and microstructural organization. However, the relation between the magnetic microstructure and the measurable Larmor frequency shift is understood only for a few idealized cases. Here we analyze the microstructure formed by magnetized, NMR-invisible infinite cylinders suspended in an NMR-reporting fluid. Through simulations, we scrutinize various geometries of mesoscopic Lorentz cavities and inclusions, and show that the cavity size should be approximately one order of magnitude larger than the width of the inclusions. We also analytically derive the Larmor frequency shift for a population of cylinders with arbitrary orientation dispersion and show that it is determined by the l = 2 Laplace expansion coefficients p 2 m of the cylinders' orientation distribution function. Our work underscores the need to account for microstructural organization when estimating magnetic tissue properties.

Ultra-short T2 components imaging of the whole brain using 3D dual-echo UTE MRI with rosette k-space pattern

PURPOSE

This study aimed to develop a new 3D dual-echo rosette k-space trajectory, specifically designed for UTE MRI applications. The imaging of the ultra-short transverse relaxation time (uT₂ ) of brain was acquired to test the performance of the proposed UTE sequence.

THEORY AND METHODS

The rosette trajectory was developed based on rotations of a "petal-like" pattern in the k_x -k_y plane, with oscillated extensions in the k_z -direction for 3D coverage. Five healthy volunteers underwent 10 dual-echo 3D rosette UTE scans with various TEs. Dual-exponential complex model fitting was performed on the magnitude data to separate uT₂ signals, with the output of uT₂ fraction, uT₂ value, and long-T₂ value.

RESULTS

The 3D rosette dual-echo UTE sequence showed better performance than a 3D radial UTE acquisition. More significant signal intensity decay in white matter than gray matter was observed along with the TEs. The white matter regions had higher uT₂ fraction values than gray matter (10.9% ± 1.9% vs. 5.7% ± 2.4%). The uT₂ value was approximately 0.10 ms in white matter .

CONCLUSION

The higher uT₂ fraction value in white matter compared to gray matter demonstrated the ability of the proposed sequence to capture rapidly decaying signals.

Diagnostic accuracy of quantitative susceptibility mapping in multiple system atrophy: The impact of echo time and the potential of histogram analysis

•

We performed histogram analysis on χ maps at different TEs on MSA patients and HC.

•

We found altered χ distribution in Pu, SN, GP, CN for MSAp and in SN, DN for MSAc.

•

QSM diagnostic accuracy is TE-dependent and is enhanced at short TEs.

•

Short TEs capture rapidly-decaying contributions of high χ sources.

•

Histogram features detect χ spatial heterogeneities improving diagnostic accuracy.

The non-invasive quantification of iron stores via Quantitative Susceptibility Mapping (QSM) could play an important role in the diagnosis and the differential diagnosis of atypical Parkinsonisms. However, the susceptibility (χ) values measured via QSM depend on echo time (TE). This effect relates to the microstructural organization within the voxel, whose composition can be altered by the disease. Moreover, pathological iron deposition in a brain area may not be spatially uniform, and conventional Region of Interest (ROI)-based analysis may fail in detecting alterations. Therefore, in this work we evaluated the impact of echo time on the diagnostic accuracy of QSM on a population of patients with Multiple System Atrophy (MSA) of either Parkinsonian (MSAp) or cerebellar (MSAc) phenotypes. In addition, we tested the potential of histogram analysis to improve QSM classification accuracy.

We enrolled 32 patients (19 MSAp and 13 MSAc) and 16 healthy controls, who underwent a 7T MRI session including a gradient-recalled multi-echo sequence for χ mapping. Nine histogram features were extracted from the χ maps computed for each TE in atlas-based ROIs covering deep brain nuclei, and compared among groups.

Alterations of susceptibility distribution were found in the Putamen, Substantia Nigra, Globus Pallidus and Caudate Nucleus for MSAp and in the Substantia Nigra and Dentate Nucleus for MSAc. Increased iron deposition was observed in a larger number of ROIs for the two shortest TEs and the standard deviation, the 75^th and the 90^th percentile were the most informative features yielding excellent diagnostic accuracy with area under the ROC curve > 0.9.

In conclusion, short TEs may enhance QSM diagnostic performances, as they can capture variations in rapidly-decaying contributions of high χ sources. The analysis of histogram features allowed to reveal fine heterogeneities in the spatial distribution of susceptibility alteration, otherwise undetected by a simple evaluation of ROI χ mean values.

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.

Learning-based motion artifact removal networks for quantitative R 2 ∗ mapping

PURPOSE

To introduce two novel learning-based motion artifact removal networks (LEARN) for the estimation of quantitative motion- and B 0 -inhomogeneity-corrected R 2 ∗ maps from motion-corrupted multi-Gradient-Recalled Echo (mGRE) MRI data.

METHODS

We train two convolutional neural networks (CNNs) to correct motion artifacts for high-quality estimation of quantitative B 0 -inhomogeneity-corrected R 2 ∗ maps from mGRE sequences. The first CNN, LEARN-IMG, performs motion correction on complex mGRE images, to enable the subsequent computation of high-quality motion-free quantitative R 2 ∗ (and any other mGRE-enabled) maps using the standard voxel-wise analysis or machine learning-based analysis. The second CNN, LEARN-BIO, is trained to directly generate motion- and B 0 -inhomogeneity-corrected quantitative R 2 ∗ maps from motion-corrupted magnitude-only mGRE images by taking advantage of the biophysical model describing the mGRE signal decay.

RESULTS

We show that both CNNs trained on synthetic MR images are capable of suppressing motion artifacts while preserving details in the predicted quantitative R 2 ∗ maps. Significant reduction of motion artifacts on experimental in vivo motion-corrupted data has also been achieved by using our trained models.

CONCLUSION

Both LEARN-IMG and LEARN-BIO can enable the computation of high-quality motion- and B 0 -inhomogeneity-corrected R 2 ∗ maps. LEARN-IMG performs motion correction on mGRE images and relies on the subsequent analysis for the estimation of R 2 ∗ maps, while LEARN-BIO directly performs motion- and B 0 -inhomogeneity-corrected R 2 ∗ estimation. Both LEARN-IMG and LEARN-BIO jointly process all the available gradient echoes, which enables them to exploit spatial patterns available in the data. The high computational speed of LEARN-BIO is an advantage that can lead to a broader clinical application.

Improved magnetic resonance myelin water imaging using multi-channel denoising convolutional neural networks (MCDnCNN)

BACKGROUND

Myelin water imaging (MWI) is powerful and important for studying and diagnosing neurological and psychiatric diseases. In particular, myelin water fraction (MWF) is derived from MWI data for quantifying myelination. However, MWF estimation is typically sensitive to noise. Improving the accuracy of MWF estimation based on WMI data acquired using a magnetic resonance (MR) multiple gradient recalled echo (mGRE) imaging sequence is desired.

METHODS

The proposed method employs a recently introduced the multi-channel denoising convolutional neural networks (MCDnCNN). Five different MCDnCNN models, denoted as Delevel1, Delevel2, Delevel3, Delevel4 and DelevelMix corresponding to five noise levels (Level1, Level2, Level3, Level4 and LevelMix), were trained using the data of the first echo of the mGRE brain images acquired from 15 healthy human subjects. Using simulated noisy data that employed a hollow cylinder model, we first evaluated the improvement in estimating MWF based on data denoised by the five different MCDnCNNs, by comparing the MWF maps calculated from the denoised data with ground truth. Next, we again evaluated the improvement using real-world in vivo datasets of 11 human participants acquired using the mGRE sequence. The datasets were first denoised by five different MCDnCNNs (Delevel1, 2, 3, 4 and DelevelMix), and subsequently their MWF maps were calculated and compared with the MWF maps directly calculated from the raw mGRE images without being denoised.

RESULTS

Experiments using the simulation data denoised by the appropriate MCDnCNN models showed that the standard deviation (SD) of the absolute error (AE) of the derived MWF results was significantly reduced (maximal reduction =15.5%, Level3 simulated noisy data, orientation angle =0, all the five MCDnCNN models). In the test using in vivo data, estimating MWF based on data particularly denoised by the appropriate MCDnCNN models was found to be the best, compared to otherwise not using the appropriate models. The results demonstrated that the appropriate MCDnCNN models may permit high-quality MWF mapping, i.e., substantial reduction of random variation in estimating MWF-maps while preserving accuracy and structural details.

CONCLUSIONS

Appropriate MCDnCNN models as proposed may improve both the accuracy and precision in estimating MWF maps, thereby making it a more clinically feasible alternative.

Quantitative susceptibility mapping and R1 measurement: Determination of the myelin volume fraction in the aging ex vivo rat corpus callosum

In studies of the white matter (WM) in aging brains, both quantitative susceptibility mapping (QSM) and direct R₁ measurement offer potentially useful ex vivo MRI tools that allow volumetric characterization of myelin content changes. Despite the technical importance of such MRI methods in numerous age-related diseases, the supposed linear relationship between the estimates of either the QSM or R₁ method and age-affected myelin contents has not been validated. In this study, the absolute myelin volume fraction (MVF) was determined by transmission electron microscopy (TEM) as a gold standard measure for comparison with the values obtained by the aforementioned MR methods. To theoretically evaluate and understand the MR signal characteristics, QSM simulations were performed using the finite perturber method (FPM). Specifically, the simulation geometry modeling was based on TEM-derived structures aligned orthogonally to the main magnetic field, the construct of which was used to estimate the magnetic field shift (ΔB) changes arising from the conjectured myelin structures. Experimentally, ex vivo corpus callosum (CC) samples from rat brains obtained at 6 weeks (n = 3), 4 months (n = 3), and 20 months (n = 3) after birth were used to establish the relationship between changes quantified by either QSM or R₁ with the absolute MVF by TEM. From the ex vivo brain samples, the scatterplot of mean MVF versus R₁ was fitted to a linear equation, where R_1mean  = 0.7948 × MVF_mean  + 0.8118 (Pearson's correlation coefficient r = 0.9138; p < 0.01), while the scatterplot of mean MVF versus MRI-derived magnetic susceptibility (χ) was also fitted to a line where χ_{measured,mean}  = -0.1218 × MVF_mean  - 0.006345 (r = -0.8435; p < 0.01). As a result of the FPM-based QSM simulations, a linearly proportional relationship between the simulated magnetic susceptibility, χ_{simulated,mean} , and MVF (r = -0.9648; p < 0.01) was established. Such a statistically significant linear correlation between MRI-derived values by the QSM (or R₁ ) method and MVF demonstrated that variable myelin contents in the WM (i.e., CC) can be quantified across multiple stages of aging. These findings further support that both techniques based on QSM and R₁ provide an efficient means of studying the brain-aging process with accurate volumetric quantification of the myelin content in WM.

Improved Multi-Echo Gradient-Echo-Based Myelin Water Fraction Mapping Using Dimensionality Reduction

Multi-echo gradient-echo (mGRE)-based myelin water fraction (MWF) mapping is a promising myelin water imaging (MWI) modality but is vulnerable to noise and artifact corruption. The linear dimensionality reduction (LDR) method has recently shown improvements with regard to these challenges. However, the magnitude value based low rank operators have been shown to misestimate the MWF for regions with [Formula: see text] anisotropy. This paper presents a nonlinear dimensionality reduction (NLDR) method to estimate the MWF map better by encouraging nonlinear low dimensionality of mGRE signal sources. Specifically, we implemented a fully connected deep autoencoder to extract the low-dimensional features of complex-valued signals and incorporated a sparse regularization to separate the anomaly sources that do not reside in the low-dimensional manifold. Simulations and in vivo experiments were performed to evaluate the accuracy of the MWF map under various situations. The proposed NLDR-based MWF improves the accuracy of the MWF map over the conventional nonlinear least-squares method and the LDR-based MWF and maintains robustness against noise and artifact corruption.

Nonparametric D-R1-R2 distribution MRI of the living human brain

Diffusion-relaxation correlation NMR can simultaneously characterize both the microstructure and the local chemical composition of complex samples that contain multiple populations of water. Recent developments on tensor-valued diffusion encoding and Monte Carlo inversion algorithms have made it possible to transfer diffusion-relaxation correlation NMR from small-bore scanners to clinical MRI systems. Initial studies on clinical MRI systems employed 5D D-R₁ and D-R₂ correlation to characterize healthy brain in vivo. However, these methods are subject to an inherent bias that originates from not including R₂ or R₁ in the analysis, respectively. This drawback can be remedied by extending the concept to 6D D-R₁-R₂ correlation. In this work, we present a sparse acquisition protocol that records all data necessary for in vivo 6D D-R₁-R₂ correlation MRI across 633 individual measurements within 25 min-a time frame comparable to previous lower-dimensional acquisition protocols. The data were processed with a Monte Carlo inversion algorithm to obtain nonparametric 6D D-R₁-R₂ distributions. We validated the reproducibility of the method in repeated measurements of healthy volunteers. For a post-therapy glioblastoma case featuring cysts, edema, and partially necrotic remains of tumor, we present representative single-voxel 6D distributions, parameter maps, and artificial contrasts over a wide range of diffusion-, R₁-, and R₂-weightings based on the rich information contained in the D-R₁-R₂ distributions.

Combined diffusion-relaxometry microstructure imaging: Current status and future prospects

Microstructure imaging seeks to noninvasively measure and map microscopic tissue features by pairing mathematical modeling with tailored MRI protocols. This article reviews an emerging paradigm that has the potential to provide a more detailed assessment of tissue microstructure-combined diffusion-relaxometry imaging. Combined diffusion-relaxometry acquisitions vary multiple MR contrast encodings-such as b-value, gradient direction, inversion time, and echo time-in a multidimensional acquisition space. When paired with suitable analysis techniques, this enables quantification of correlations and coupling between multiple MR parameters-such as diffusivity, T 1 , T 2 , and T 2 ∗ . This opens the possibility of disentangling multiple tissue compartments (within voxels) that are indistinguishable with single-contrast scans, enabling a new generation of microstructural maps with improved biological sensitivity and specificity.

Pure balanced steady-state free precession imaging (pure bSSFP)

Purpose

To show that for tissues the conspicuous asymmetries in the frequency response function of bSSFP can be mitigated by using a short enough TR.

Theory and Methods

Configuration theory indicates that bSSFP becomes apparently “pure” (i.e., exhibiting a symmetric profile) in the limit of TR →0. To this end, the frequency profile of bSSFP was measured as a function of the TR using a manganese‐doped aqueous probe, as well as brain tissue that was shown to exhibit a pronounced asymmetry due to its microstructure. The frequency response function was sampled using N=72 (phantom) and N=36 (in vivo) equally distributed linear RF phase increments in the interval [0,2π). Imaging was performed with 2.0 mm isotropic resolution over a TR range of 1.5–8 ms at 3 and 1.5 T.

Results

As expected, pure substances showed a symmetric TR‐independent frequency profile, whereas brain tissue revealed a pronounced asymmetry. The observed asymmetry for the tissue, however, decreases with decreasing TR and gives strong evidence that the frequency response function of bSSFP becomes symmetric in the limit of TR →0, in agreement with theory. The limit of apparently pure bSSFP imaging can thus be achieved for a TR ∼ 1.5 ms at 1.5 T, whereas at 3 T, tissues still show some residual asymmetry.

Conclusion

In the limit of short enough TR, tissues become apparently pure for bSSFP. This limit can be reached for brain tissue at 1.5 T with TR ∼ 1–2 ms at clinically relevant resolutions.

Quantifying myelin in crossing fibers using diffusion-prepared phase imaging: Theory and simulations

PURPOSE

Myelin has long been the target of neuroimaging research. However, most available techniques can only provide a voxel-averaged estimate of myelin content. In the human brain, white matter fiber pathways connecting different brain areas and carrying different functions often cross each other in the same voxel. A measure that can differentiate the degree of myelination of crossing fibers would provide a more specific marker of myelination.

THEORY AND METHODS

One MRI signal property that is sensitive to myelin is the phase accumulation. This sensitivity is used by measuring the phase accumulation of the signal remaining after diffusion-weighting, which is called diffusion-prepared phase imaging (DIPPI). Including diffusion-weighting before estimating the phase accumulation has two distinct advantages for estimating the degree of myelination: (1) It increases the relative contribution of intra-axonal water, whose phase is related linearly to the thickness of the surrounding myelin (in particular the log g-ratio); and (2) it gives directional information, which can be used to distinguish between crossing fibers. Here the DIPPI sequence is described, an approach is proposed to estimate the log g-ratio, and simulations are used and DIPPI data acquired in an isotropic phantom to quantify other sources of phase accumulation.

RESULTS

The expected bias is estimated in the log g-ratio for reasonable in vivo acquisition parameters caused by eddy currents (~4%-10%), remaining extra-axonal signal (~15%), and gradients in the bulk off-resonance field (<10% for most of the brain).

CONCLUSION

This new sequence may provide a g-ratio estimate per fiber population crossing within a voxel.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Measuring compartmental T2-orientational dependence in human brain white matter using a tiltable RF coil and diffusion-T2 correlation MRI

The anisotropy of brain white matter microstructure manifests itself in orientational-dependence of various MRI contrasts, and can result in significant quantification biases if ignored. Understanding the origins of this orientation-dependence could enhance the interpretation of MRI signal changes in development, ageing and disease and ultimately improve clinical diagnosis. Using a novel experimental setup, this work studies the contributions of the intra- and extra-axonal water to the orientation-dependence of one of the most clinically-studied parameters, apparent transverse relaxation T2. Specifically, a tiltable receive coil is interfaced with an ultra-strong gradient MRI scanner to acquire multidimensional MRI data with an unprecedented range of acquisition parameters. Using this setup, compartmental T2 can be disentangled based on differences in diffusional-anisotropy, and its orientation-dependence further elucidated by re-orienting the head with respect to the main magnetic field B→0. A dependence of (compartmental) T2 on the fibre orientation w.r.t. B→0 was observed, and further quantified using characteristic representations for susceptibility- and magic angle effects. Across white matter, anisotropy effects were dominated by the extra-axonal water signal, while the intra-axonal water signal decay varied less with fibre-orientation. Moreover, the results suggest that the stronger extra-axonal T2 orientation-dependence is dominated by magnetic susceptibility effects (presumably from the myelin sheath) while the weaker intra-axonal T2 orientation-dependence may be driven by a combination of microstructural effects. Even though the current design of the tiltable coil only offers a modest range of angles, the results demonstrate an overall effect of tilt and serve as a proof-of-concept motivating further hardware development to facilitate experiments that explore orientational anisotropy. These observations have the potential to lead to white matter microstructural models with increased compartmental sensitivity to disease, and can have direct consequences for longitudinal and group-wise T2- and diffusion-MRI data analysis, where the effect of head-orientation in the scanner is commonly ignored.

In vivo evaluation of heme and non-heme iron content and neuronal density in human basal ganglia

Non-heme iron is an important element supporting the structure and functioning of biological tissues. Imbalance in non-heme iron can lead to different neurological disorders. Several MRI approaches have been developed for iron quantification relying either on the relaxation properties of MRI signal or measuring tissue magnetic susceptibility. Specific quantification of the non-heme iron can, however, be constrained by the presence of the heme iron in the deoxygenated blood and contribution of cellular composition. The goal of this paper is to introduce theoretical background and experimental MRI method allowing disentangling contributions of heme and non-heme irons simultaneously with evaluation of tissue neuronal density in the iron-rich basal ganglia. Our approach is based on the quantitative Gradient Recalled Echo (qGRE) MRI technique that allows separation of the total R2* metric characterizing decay of GRE signal into tissue-specific (R2t*) and the baseline blood oxygen level-dependent (BOLD) contributions. A combination with the QSM data (also available from the qGRE signal phase) allowed further separation of the tissue-specific R2t* metric in a cell-specific and non-heme-iron-specific contributions. It is shown that the non-heme iron contribution to R2t* relaxation can be described with the previously developed Gaussian Phase Approximation (GPA) approach. qGRE data were obtained from 22 healthy control participants (ages 26-63 years). Results suggest that the ferritin complexes are aggregated in clusters with an average radius about 100nm comprising approximately 2600 individual ferritin units. It is also demonstrated that the concentrations of heme and non-heme iron tend to increase with age. The strongest age effect was seen in the pallidum region, where the highest age-related non-heme iron accumulation was observed.

χ-separation: Magnetic susceptibility source separation toward iron and myelin mapping in the brain

Obtaining a histological fingerprint from the in-vivo brain has been a long-standing target of magnetic resonance imaging (MRI). In particular, non-invasive imaging of iron and myelin, which are involved in normal brain functions and are histopathological hallmarks in neurodegenerative diseases, has practical utilities in neuroscience and medicine. Here, we propose a biophysical model that describes the individual contribution of paramagnetic (e.g., iron) and diamagnetic (e.g., myelin) susceptibility sources to the frequency shift and transverse relaxation of MRI signals. Using this model, we develop a method, χ-separation, that generates the voxel-wise distributions of the two sources. The method is validated using computer simulation and phantom experiments, and applied to ex-vivo and in-vivo brains. The results delineate the well-known histological features of iron and myelin in the specimen, healthy volunteers, and multiple sclerosis patients. This new technology may serve as a practical tool for exploring the microstructural information of the brain.

Improved estimation of myelin water fractions with learned parameter distributions

PURPOSE

To improve estimation of myelin water fraction (MWF) in the brain from multi-echo gradient-echo imaging data.

METHODS

A systematic sensitivity analysis was first conducted to characterize the conventional exponential models used for MWF estimation. A new estimation method was then proposed for improved estimation of MWF from practical gradient-echo imaging data. The proposed method uses an extended signal model that includes a finite impulse response filter to compensate for practical signal variations. This new model also enables the use of prelearned parameter distributions as well as low-rank signal structures to improve parameter estimation. The resulting parameter estimation problem was solved optimally in the Bayesian sense.

RESULTS

Our sensitivity analysis results showed that the conventional exponential models were very sensitive to measurement noise and modeling errors. Our simulation and experimental results showed that our proposed method provided a substantial improvement in reliability, reproducibility, and robustness of MWF estimates over the conventional methods. Clinical results obtained from stroke patients indicated that the proposed method, with its improved capability, could reveal the loss of myelin in lesions, demonstrating its translational potentials.

CONCLUSION

This paper addressed the problem of robust MWF estimation from gradient-echo imaging data. A new method was proposed to provide improved MWF estimation in the presence of significant noise and modeling errors. The performance of the proposed method has been evaluated using both simulated and experimental data, showing significantly improved robustness over the existing methods. The proposed method may prove useful for quantitative myelin imaging in clinical applications.

Single-step calculation of susceptibility through multiple orientation sampling

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

Sensitivity of fiber orientation dependent R 2 ∗ to temperature and post mortem interval

PURPOSE

R 2 ∗ imaging of brain white matter is well known for being sensitive to the orientation of nerve fibers with respect to the B₀ field of the MRI scanner. The goal of this study was to evaluate whether and to which extent fiber orientation dependent R 2 ∗ differs between in vivo and post mortem in situ examinations, and to investigate the influence of varying temperatures and post mortem intervals (PMI).

METHODS

Post mortem in situ and in vivo MRI scans were conducted at 3T. R 2 ∗ was acquired with a multi-echo gradient-echo sequence, and the orientation of white matter fibers was computed using diffusion tensor imaging (DTI). Fitting of the measured fiber orientation dependent R 2 ∗ was performed using three different formulations of a previously proposed model.

RESULTS

R 2 ∗ increased with increasing fiber angle for in vivo and post mortem in situ examinations, whereby the orientation dependency was lower post mortem. The different formulations of the fiber orientation model resulted in an identical fit, but showed large variations of the estimated parameters. The higher order orientation dependent R 2 ∗ components significantly decreased with decreasing temperature, while the orientation independent R 2 ∗ components showed no significant correlation with either temperature or PMI.

CONCLUSION

Although the mean diffusivity is strongly reduced post mortem, we could successfully estimate the fiber angle using DTI. Due to the strong correlation of the higher order orientation dependent R 2 ∗ components with temperature, the decreased R 2 ∗ fiber orientation dependency post mortem in situ might primarily be attributed to the lower brain temperature.

Decoding the microstructural properties of white matter using realistic models

Multi-echo gradient echo (ME-GRE) magnetic resonance signal evolution in white matter has a strong dependence on the orientation of myelinated axons with respect to the main static field. Although analytical solutions have been able to predict some of the white matter (WM) signal behaviour of the hollow cylinder model, it has been shown that realistic models of WM offer a better description of the signal behaviour observed. In this work, we present a pipeline to (i) generate realistic 2D WM models with their microstructure based on real axon morphology with adjustable fiber volume fraction (FVF) and g-ratio. We (ii) simulate their interaction with the static magnetic field to be able to simulate their MR signal. For the first time, we (iii) demonstrate that realistic 2D WM models can be used to simulate a MR signal that provides a good approximation of the signal obtained from a real 3D WM model derived from electron microscopy. We then (iv) demonstrate in silico that 2D WM models can be used to predict microstructural parameters in a robust way if ME-GRE multi-orientation data is available and the main fiber orientation in each pixel is known using DTI. A deep learning network was trained and characterized in its ability to recover the desired microstructural parameters such as FVF, g-ratio, free and bound water transverse relaxation and magnetic susceptibility. Finally, the network was trained to recover these micro-structural parameters from an ex vivo dataset acquired in 9 orientations with respect to the magnetic field and 12 echo times. We demonstrate that this is an overdetermined problem and that as few as 3 orientations can already provide comparable results for some of the decoded metrics.

Multi-compartment analysis of the complex gradient-echo signal quantifies myelin breakdown in premanifest Huntington's disease

White matter (WM) alterations have been identified as a relevant pathological feature of Huntington's disease (HD). Increasing evidence suggests that WM changes in this disorder are due to alterations in myelin-associated biological processes. Multi-compartmental analysis of the complex gradient-echo MRI signal evolution in WM has been shown to quantify myelin in vivo, therefore pointing to the potential of this technique for the study of WM myelin changes in health and disease. This study first characterized the reproducibility of metrics derived from the complex multi-echo gradient-recalled echo (mGRE) signal across the corpus callosum in healthy participants, finding highest reproducibility in the posterior callosal segment. Subsequently, the same analysis pipeline was applied in this callosal region in a sample of premanifest HD patients (n = 19) and age, sex and education matched healthy controls (n = 21). In particular, we focused on two myelin-associated derivatives: i. the myelin water signal fraction (fm), a parameter dependent on myelin content; and ii. The difference in frequency between myelin and intra-axonal water pools (Δω), a parameter dependent on the ratio between the inner and the outer axonal radii. fm was found to be lower in HD patients (β = -0.13, p = 0.03), while Δω did not show a group effect. Performance in tests of working memory, executive function, social cognition and movement was also assessed, and a greater age-related decline in executive function was detected in HD patients (β = -0.06, p = 0.006), replicating previous evidence of executive dysfunction in HD. Finally, the correlation between fm, executive function, and proximity to disease onset was explored in patients, and a positive correlation between executive function and fm was detected (r = 0.542; p = 0.02). This study emphasises the potential of complex mGRE signal analysis for aiding understanding of HD pathogenesis and progression. Moreover, expanding on evidence from pathology and animal studies, it provides novel in vivo evidence supporting myelin breakdown as an early feature of HD.

Variable flip angle echo planar time-resolved imaging (vFA-EPTI) for fast high-resolution gradient echo myelin water imaging

Myelin water imaging techniques based on multi-compartment relaxometry have been developed as an important tool to measure myelin concentration in vivo, but are limited by the long scan time of multi-contrast multi-echo acquisition. In this work, a fast imaging technique, termed variable flip angle Echo Planar Time-Resolved Imaging (vFA-EPTI), is developed to acquire multi-echo and multi-flip-angle gradient-echo data with significantly reduced acquisition time, providing rich information for multi-compartment analysis of gradient-echo myelin water imaging (GRE-MWI). The proposed vFA-EPTI method achieved 26 folds acceleration with good accuracy by utilizing an efficient continuous readout, optimized spatiotemporal encoding across echoes and flip angles, as well as a joint subspace reconstruction. An approach to estimate off-resonance field changes between different flip-angle acquisitions was also developed to ensure high-quality joint reconstruction across flip angles. The accuracy of myelin water fraction (MWF) estimate under high acceleration was first validated by a retrospective undersampling experiment using a lengthy fully-sampled data as reference. Prospective experiments were then performed where whole-brain MWF and multi-compartment quantitative maps were obtained in 5 min at 1.5 mm isotropic resolution and 24 min at 1 mm isotropic resolution at 3T. Additionally, ultra-high resolution data at 600 μm isotropic resolution were acquired at 7T, which show detailed structures within the cortex such as the line of Gennari, demonstrating the ability of the proposed method for submillimeter GRE-MWI that can be used to study cortical myeloarchitecture in vivo.

Probing the myelin water compartment with a saturation-recovery, multi-echo gradient-recalled echo sequence

PURPOSE

To investigate the effect of varying levels of T 1 -weighting on the evolution of the complex signal from white matter in a multi-echo gradient-recalled echo (mGRE) saturation-recovery sequence.

THEORY AND METHODS

Analysis of the complex signal evolution in an mGRE sequence allows the contributions from short- and long- T 2 ∗ components to be separated, thus providing a measure of the relative strength of signals from the myelin water, and the external and intra-axonal compartments. Here we evaluated the effect of different levels of T 1 -weighting on these signals, expecting that the previously reported, short T 1 of the myelin water would lead to a relative enhancement of the myelin water signal in the presence of signal saturation. Complex, saturation-recovery mGRE data from the splenium of the corpus callosum from 5 healthy volunteers were preprocessed using a frequency difference mapping (FDM) approach and analyzed using the 3-pool model of complex signal evolution in white matter.

RESULTS

An increase in the apparent T 1 as a function of echo time was demonstrated, but this increase was an order of magnitude smaller than that expected from previously reported myelin water T 1 -values. This suggests the presence of magnetization transfer and exchange effects which counteract the T 1 -weighting.

CONCLUSION

Variation of the B 1 + amplitude in a saturation-recovery mGRE sequence can be used to modulate the relative strength of signals from the different compartments in white matter, but the modulation is less than predicted from previously reported T 1 -values.

Understanding aqueous and non-aqueous proton T1 relaxation in brain

A full picture of longitudinal relaxation in complex heterogeneous environments like white matter brain tissue remains elusive. In tissue, successive approximations, from the solvation layer model to the two pool model, have highlighted how longitudinal magnetization evolution depends on both inter-compartmental exchange and spin-lattice relaxation. In white matter, however, these models fail to capture the behaviour of the two distinct aqueous pools, myelin water and intra/extra-cellular water. A challenge with testing more comprehensive multi-pool models lies in directly observing all pools, both aqueous and non-aqueous. In this work, we advance these efforts by integrating three main experimental and analytical elements: direct observation of the longitudinal relaxation of both the aqueous and the non-aqueous protons in white matter, a wide range of different initial conditions, and application of an analysis pipeline which includes lineshape, CPMG, and fitting of a four pool model. An eigenvector interpretation of the four pool model highlights how longitudinal relaxation in white matter depends on initial conditions. We find that a single set of model parameters is able to describe the entire range of relaxation behaviour observed in all the separable aqueous and non-aqueous pools in experiments involving six different initial conditions. Understanding of the nature and connectedness of the tissue components is crucial in the design and interpretation of many MRI measurements, especially those based on magnetization transfer and longitudinal relaxation. In particular, the dependency of relaxation behaviour on initial conditions is likely the basis for understanding method-dependent discrepancies in in vivo T₁.

A Physical Perspective to the Inductive Function of Myelin-A Missing Piece of Neuroscience

Starting from the inductance in neurons, two physical origins are discussed, which are the coil inductance of myelin and the piezoelectric effect of the cell membrane. The direct evidence of the coil inductance of myelin is the opposite spiraling phenomenon between adjacent myelin sheaths confirmed by previous studies. As for the piezoelectric effect of the cell membrane, which has been well-known in physics, the direct evidence is the mechanical wave accompany with action potential. Therefore, a more complete physical nature of neural signals is provided. In conventional neuroscience, the neural signal is a pure electrical signal. In our new theory, the neural signal is an energy pulse containing electrical, magnetic, and mechanical components. Such a physical understanding of the neural signal and neural systems significantly improve the knowledge of the neurons. On the one hand, we achieve a corrected neural circuit of an inductor-capacitor-capacitor (LCC) form, whose frequency response and electrical characteristics have been validated by previous studies and the modeling fitting of artifacts in our experiments. On the other hand, a number of phenomena observed in neural experiments are explained. In particular, they are the mechanism of magnetic nerve stimulations and ultrasound nerve stimulations, the MRI image contrast issue and Anode Break Excitation. At last, the biological function of myelin is summarized. It is to provide inductance in the process of neural signal, which can enhance the signal speed in peripheral nervous systems and provide frequency modulation function in central nervous systems.

Towards in vivo g-ratio mapping using MRI: Unifying myelin and diffusion imaging

BACKGROUND

The g-ratio, quantifying the comparative thickness of the myelin sheath encasing an axon, is a geometrical invariant that has high functional relevance because of its importance in determining neuronal conduction velocity. Advances in MRI data acquisition and signal modelling have put in vivo mapping of the g-ratio, across the entire white matter, within our reach. This capacity would greatly increase our knowledge of the nervous system: how it functions, and how it is impacted by disease.

NEW METHOD

This is the second review on the topic of g-ratio mapping using MRI.

RESULTS

This review summarizes the most recent developments in the field, while also providing methodological background pertinent to aggregate g-ratio weighted mapping, and discussing pitfalls associated with these approaches.

COMPARISON WITH EXISTING METHODS

Using simulations based on recently published data, this review reveals caveats to the state-of-the-art calibration methods that have been used for in vivo g-ratio mapping. It highlights the need to estimate both the slope and offset of the relationship between these MRI-based markers and the true myelin volume fraction if we are really to achieve the goal of precise, high sensitivity g-ratio mapping in vivo. Other challenges discussed in this review further evidence the need for gold standard measurements of human brain tissue from ex vivo histology.

CONCLUSIONS

We conclude that the quest to find the most appropriate MRI biomarkers to enable in vivo g-ratio mapping is ongoing, with the full potential of many novel techniques yet to be investigated.

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.

Parsimonious discretization for characterizing multi-exponential decay in magnetic resonance

We address the problem of analyzing noise-corrupted magnetic resonance transverse decay signals as a superposition of underlying independently decaying monoexponentials of positive amplitude. First, we indicate the manner in which this is an ill-conditioned inverse problem, rendering the analysis unstable with respect to noise. Second, we define an approach to this analysis, stabilized solely by the nonnegativity constraint without regularization. This is made possible by appropriate discretization, which is coarser than that often used in practice. Thirdly, we indicate further stabilization by inspecting the plateaus of cumulative distributions. We demonstrate our approach through analysis of simulated myelin water fraction measurements, and compare the accuracy with more conventional approaches. Finally, we apply our method to brain imaging data obtained from a human subject, showing that our approach leads to maps of the myelin water fraction which are much more stable with respect to increasing noise than those obtained with conventional approaches.

Non-negative least squares computation for in vivo myelin mapping using simulated multi-echo spin-echo T2 decay data

Multi-compartment T₂ mapping has gained particular relevance for the study of myelin water in the brain. As a facilitator of rapid saltatory axonal signal transmission, myelin is a cornerstone indicator of white matter development and function. Regularized non-negative least squares fitting of multi-echo T₂ data has been widely employed for the computation of the myelin water fraction (MWF), and the obtained MWF maps have been histopathologically validated. MWF measurements depend upon the quality of the data acquisition, B₁⁺ homogeneity and a range of fitting parameters. In this special issue article, we discuss the relevance of these factors for the accurate computation of multi-compartment T₂ and MWF maps. We generated multi-echo spin-echo T₂ decay curves following the Carr-Purcell-Meiboom-Gill approach for various myelin concentrations and myelin T₂ scenarios by simulating the evolution of the magnetization vector between echoes based on the Bloch equations. We demonstrated that noise and imperfect refocusing flip angles yield systematic underestimations in MWF and intra-/extracellular water geometric mean T₂ (gmT₂ ). MWF estimates were more stable than myelin water gmT₂ time across different settings of the T₂ analysis. We observed that the lower limit of the T₂ distribution grid should be slightly shorter than TE₁ . Both TE₁ and the acquisition echo spacing also have to be sufficiently short to capture the rapidly decaying myelin water T₂ signal. Among all parameters of interest, the estimated MWF and intra-/extracellular water gmT₂ differed by approximately 0.13-4 percentage points and 3-4 ms, respectively, from the true values, with larger deviations observed in the presence of greater B₁⁺ inhomogeneities and at lower signal-to-noise ratio. Tailoring acquisition strategies may allow us to better characterize the T₂ distribution, including the myelin water, in vivo.

Microscopic susceptibility anisotropy imaging

PURPOSE

The gradient-echo MR signal in brain white matter depends on the orientation of the fibers with respect to the external magnetic field. To map microstructure-specific magnetic susceptibility in orientationally heterogeneous material, it is thus imperative to regress out unwanted orientation effects.

METHODS

This work introduces a novel framework, referred to as microscopic susceptibility anisotropy imaging, that disentangles the 2 principal effects conflated in gradient-echo measurements, (a) the susceptibility properties of tissue microenvironments, especially the myelin microstructure, and (b) the axon orientation distribution relative to the magnetic field. Specifically, we utilize information about the orientational tissue structure inferred from diffusion MRI data to factor out the B 0 -direction dependence of the frequency difference signal.

RESULTS

A human pilot study at 3 T demonstrates proxy maps of microscopic susceptibility anisotropy unconfounded by fiber crossings and orientation dispersion as well as magnetic field direction. The developed technique requires only a dual-echo gradient-echo scan acquired at 1 or 2 head orientations with respect to the magnetic field and a 2-shell diffusion protocol achievable on standard scanners within practical scan times.

CONCLUSIONS

The quantitative recovery of microscopic susceptibility features in the presence of orientational heterogeneity potentially improves the assessment of microstructural tissue integrity.

Unraveling the contributions to the neuromelanin-MRI contrast

The Locus Coeruleus (LC) and the Substantia Nigra (SN) are small brainstem nuclei that change with aging and may be involved in the development of various neurodegenerative and psychiatric diseases. Magnetization Transfer (MT) MRI has been shown to facilitate LC and the SN visualization, and the observed contrast is assumed to be related to neuromelanin accumulation. Imaging these nuclei may have predictive value for the progression of various diseases, but interpretation of previous studies is hindered by the fact that the precise biological source of the contrast remains unclear, though several hypotheses have been put forward. To inform clinical studies on the possible biological interpretation of the LC- and SN contrast, we examined an agar-based phantom containing samples of natural Sepia melanin and synthetic Cys-Dopa-Melanin and compared this to the in vivo human LC and SN. T₁ and T₂* maps, MT spectra and relaxation times of the phantom, the LC and the SN were measured, and a two-pool MT model was fitted. Additionally, Bloch simulations and a transient MT experiment were conducted to confirm the findings. Overall, our results indicate that Neuromelanin-MRI contrast in the LC likely results from a lower macromolecular fraction, thus facilitating interpretation of results in clinical populations. We further demonstrate that in older individuals T₁ lengthening occurs in the LC.