Physiological noise compensation in gradient-echo myelin water imaging

Whole brain multiparametric mapping in two minutes using a dual-flip-angle stack-of-stars blipped multi-gradient-echo acquisition

A new MRI technique is presented for three-dimensional fast simultaneous whole brain mapping of myelin water fraction (MWF), T1, proton density (PD), R2*, magnetic susceptibility (QSM), and B1 transmit field (B1+). Phantom and human (N = 9) datasets were acquired using a dual-flip-angle blipped multi-gradient-echo (DFA-mGRE) sequence with a stack-of-stars (SOS) trajectory. Images were reconstructed using a subspace-based algorithm with a locally low-rank constraint. A novel joint-sparsity-constrained multicomponent T2*-T1 spectrum estimation (JMSE) algorithm is proposed to correct for the T1 saturation effect and B1+/B1- inhomogeneities in the quantification of MWF. A tissue-prior-based B1+ estimation algorithm was adapted for B1 correction in the mapping of T1 and PD. In the phantom study, measurements obtained at an acceleration factor (R) of 12 using prospectively under-sampled SOS showed good consistency (R2 > 0.997) with Cartesian reference for R2*/T1app/M0app. In the in vivo study, results of retrospectively under-sampled SOS with R = 6, 12, 18, showed good quality (structure similarity index measure > 0.95) compared with those of fully-sampled SOS. Besides, results of prospectively under-sampled SOS with R = 12 showed good consistency (intraclass correlation coefficient > 0.91) with Cartesian reference for T1/PD/B1+/MWF/QSM/R2*, and good reproducibility (coefficient of variation < 7.0 %) in the test-retest analysis for T1/PD/B1+/MWF/R2*. This study has demonstrated the feasibility of simultaneous whole brain multiparametric mapping with a two-minute scan using the DFA-mGRE SOS sequence, which may overcome a major obstacle for neurological applications of multiparametric MRI.

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.

JUST-Net: Jointly unrolled cross-domain optimization based spatio-temporal reconstruction network for accelerated 3D myelin water imaging

PURPOSE

We introduced a novel reconstruction network, jointly unrolled cross-domain optimization-based spatio-temporal reconstruction network (JUST-Net), aimed at accelerating 3D multi-echo gradient-echo (mGRE) data acquisition and improving the quality of resulting myelin water imaging (MWI) maps.

METHOD

An unrolled cross-domain spatio-temporal reconstruction network was designed. The main idea is to combine frequency and spatio-temporal image feature representations and to sequentially implement convolution layers in both domains. The k-space subnetwork utilizes shared information from adjacent frames, whereas the image subnetwork applies separate convolutions in both spatial and temporal dimensions. The proposed reconstruction network was evaluated for both retrospectively and prospectively accelerated acquisition. Furthermore, it was assessed in simulation studies and real-world cases with k-space corruptions to evaluate its potential for motion artifact reduction.

RESULTS

The proposed JUST-Net enabled highly reproducible and accelerated 3D mGRE acquisition for whole-brain MWI, reducing the acquisition time from fully sampled 15:23 to 2:22 min within a 3-min reconstruction time. The normalized root mean squared error of the reconstructed mGRE images increased by less than 4.0%, and the correlation coefficients for MWI showed a value of over 0.68 when compared to the fully sampled reference. Additionally, the proposed method demonstrated a mitigating effect on both simulated and clinical motion-corrupted cases.

CONCLUSION

The proposed JUST-Net has demonstrated the capability to achieve high acceleration factors for 3D mGRE-based MWI, which is expected to facilitate widespread clinical applications of MWI.

Quantitative magnetic resonance imaging reflects different levels of histologically determined myelin densities in multiple sclerosis, including remyelination in inactive multiple sclerosis lesions

Magnetic resonance imaging (MRI) of focal or diffuse myelin damage or remyelination may provide important insights into disease progression and potential treatment efficacy in multiple sclerosis (MS). We performed post-mortem MRI and histopathological myelin measurements in seven progressive MS cases to evaluate the ability of three myelin-sensitive MRI scans to distinguish different stages of MS pathology, particularly chronic demyelinated and remyelinated lesions. At 3 Tesla, we acquired two different myelin water imaging (MWI) scans and magnetisation transfer ratio (MTR) data. Histopathology included histochemical stainings for myelin phospholipids (LFB) and iron as well as immunohistochemistry for myelin proteolipid protein (PLP), CD68 (phagocytosing microglia/macrophages) and BCAS1 (remyelinating oligodendrocytes). Mixed-effects modelling determined which histopathological metric best predicted MWF and MTR in normal-appearing and diffusely abnormal white matter, active/inactive, inactive, remyelinated and ischemic lesions. Both MWI measures correlated well with each other and histology across regions, reflecting the different stages of MS pathology. MTR data showed a considerable influence of components other than myelin and a strong dependency on tissue storage duration. Both MRI and histology revealed increased myelin densities in inactive compared with active/inactive lesions. Chronic inactive lesions harboured single scattered myelin fibres indicative of low-level remyelination. Mixed-effects modelling showed that smaller differences between white matter areas were linked to PLP densities and only to a small extent confounded by iron. MWI reflects differences in myelin lipids and proteins across various levels of myelin densities encountered in MS, including low-level remyelination in chronic inactive lesions.

On the performance of multi-compartment relaxometry for myelin water imaging (MCR-MWI) - test-retest repeatability and inter-protocol reproducibility

In this study, we optimized the variable flip angle (VFA) acquisition scheme using numerical simulations to shorten the acquisition time of multicompartment relaxometry for myelin water imaging (MCR-MWI) to a clinically practical range in the absence of advanced image reconstruction methods. As the primary objective of this study, the test-retest repeatability of myelin water fraction (MWF) measurements of MCR-MWI is evaluated on three gradient echo (GRE) sequence settings using the optimized VFA schemes with different echo times and repetition times, emulating various scanner setups. The cross-protocol reproducibility of MCR-MWI and MCR with diffusion-informed myelin water imaging (MCR-DIMWI) is also examined. As a secondary objective, we explore the bundle-specific profiles of various microstructural parameters from MCR-(DI)MWI and their cross-correlations to determine if these parameters possess supplementary microstructure information beyond myelin concentration. Numerical simulations indicate that MCR-MWI can be performed with a minimum of three flip angles covering a wide range of T₁ weightings without adding significant bias. This is supported by the results of an in vivo experiment, allowing whole-brain 1.5 mm isotropic MWF maps to be acquired in 9 min, reducing the total scan time to 40% of the original implementation without significant quality degradation. Good test-retest repeatability is observed for MCR-MWI for all three GRE protocols. While good correlations can also be found in MWF across protocols, systematic differences are observed. Bundle-specific MWF analysis reveals that certain white matter bundles are similar in all participants. We also found that microstructure relaxation parameters have low linear correlations with MWF. MCR-MWI is a reproducible measure of myelin. However, attention should be paid to the protocol related MWF differences when comparing different studies, as the MWF bias up to 0.5% can be observed across the protocols examined in this work.

Mitigating the impact of flip angle and orientation dependence in single compartment R2* estimates via 2-pool modeling

PURPOSE

The effective transverse relaxation rate (R 2 * $$ {\mathrm{R}}_2^{\ast } $$ ) is influenced by biological features that make it a useful means of probing brain microstructure. However, confounding factors such as dependence on flip angle (α) and fiber orientation with respect to the main field (θ $$ \uptheta $$ ) complicate interpretation. The α- andθ $$ \uptheta $$ -dependence stem from the existence of multiple sub-voxel micro-environments (e.g., myelin and non-myelin water compartments). Ordinarily, it is challenging to quantify these sub-compartments; therefore, neuroscientific studies commonly make the simplifying assumption of a mono-exponential decay obtaining a singleR 2 * $$ {\mathrm{R}}_2^{\ast } $$ estimate per voxel. In this work, we investigated how the multi-compartment nature of tissue microstructure affects single compartmentR 2 * $$ {\mathrm{R}}_2^{\ast } $$ estimates.

METHODS

We used 2-pool (myelin and non-myelin water) simulations to characterize the bias in single compartmentR 2 * $$ {\mathrm{R}}_2^{\ast } $$ estimates. Based on our numeric observations, we introduced a linear model that partitionsR 2 * $$ {\mathrm{R}}_2^{\ast } $$ into α-dependent and α-independent components and validated this in vivo at 7T. We investigated the dependence of both components on the sub-compartment properties and assessed their robustness, orientation dependence, and reproducibility empirically.

RESULTS

R 2 * $$ {\mathrm{R}}_2^{\ast } $$ increased with myelin water fraction and residency time leading to a linear dependence on α. We observed excellent agreement between our numeric and empirical results. Furthermore, the α-independent component of the proposed linear model was robust to the choice of α and reduced dependence on fiber orientation, although it suffered from marginally higher noise sensitivity.

CONCLUSION

We have demonstrated and validated a simple approach that mitigates flip angle and orientation biases in single-compartmentR 2 * $$ {\mathrm{R}}_2^{\ast } $$ estimates.

Myelin Water Imaging of Nerve Recovery in Rehabilitating Stroke Patients

BACKGROUND

Myelin water imaging (MWI) using MRI has been introduced as a method to quantify the integrity of myelin in vivo. However, the investigation of its potential to probe myelin changes has been limited.

PURPOSE

To determine the myelin change using MWI in the corticospinal tract (CST) during the rehabilitation of stroke patients.

STUDY TYPE

Longitudinal.

POPULATION

A total of 24 stroke patients within 6 months from the onset (64.3 ± 16.1 years, 14 women, 10 men) and 10 healthy volunteers (27.0 ± 2.2 years, 2 women, 8 men).

FIELD STRENGTH/SEQUENCE

Three-dimensional multiecho gradient echo sequence and diffusion-weighted echoplanar imaging sequence at 3 T.

ASSESSMENT

The changes of myelin water fraction (MWF) and fractional anisotropy (FA) during rehabilitation were analyzed in the CST and other regions using tractography software and region of interest drawings by the radiologist.

STATISTICAL TESTS

A paired t-test was performed to investigate the change of MRI metrics during rehabilitation. In addition, an independent two-sample t-test was performed to investigate the effects of different rehabilitation protocols. A P-value <0.05 was considered significant.

RESULTS

In the CST, MWF significantly changed from 5.83 ± 0.91% to 6.23 ± 0.97% after rehabilitation while changes of FA (0.442 ± 0.038 to 0.443 ± 0.035) were not significant (P = 0.656). The rate of change in MWF and FA, which were 6.69% and 0.439% respectively, were significantly different. Other regions did not show significant changes (range of MWF change: -3.44% to -1.61%, range of FA change: -1.39% to 0.79%, and range of P-value: 0.144-0.761). Further analysis showed that those with additional robot-assisted rehabilitation had a significantly larger MWF change than those with conventional rehabilitation only (rate of change: 11.2% vs. 3.2%).

DATA CONCLUSION

The feasibility of using MWI to monitor myelin content was demonstrated by showing the MWF changes during rehabilitation.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY STAGE: 2.

Accelerated 3D myelin water imaging using joint spatio-temporal reconstruction

PURPOSE

To enable acceleration in 3D multi-echo gradient echo (mGRE) acquisition for myelin water imaging (MWI) by combining joint parallel imaging (JPI) and joint deep learning (JDL).

METHODS

We implemented a multistep reconstruction process using both advanced parallel imaging and deep learning network which can utilize joint spatiotemporal components between the multi-echo images to further accelerate 3D mGRE acquisition for MWI. In the first step, JPI was performed to estimate missing k-space lines. Next, JDL was implemented to reduce residual artifacts and produce high-fidelity reconstruction by using variable splitting optimization consisting of spatiotemporal denoiser block, data consistency block, and weighted average block. The proposed method was evaluated for MWI with 2D Cartesian uniform under-sampling for each echo, enabling scan times of up to approximately 2 min for 2 mm × 2 mm × 2 mm $2\ {\rm mm} \times 2\ {\rm mm} \times 2\ {\rm mm}$ 3D coverage.

RESULTS

The proposed method showed acceptable MWI quality with improved quantitative values compared to both JPI and JDL methods individually. The improved performance of the proposed method was demonstrated by the low normalized mean-square error and high-frequency error norm values of the reconstruction with high similarity to the fully sampled MWI.

CONCLUSION

Joint spatiotemporal reconstruction approach by combining JPI and JDL can achieve high acceleration factors for 3D mGRE-based MWI.

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.

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.

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.

Whole Brain Myelin Water Mapping in One Minute Using Tensor Dictionary Learning With Low-Rank Plus Sparse Regularization

The quantification of myelin water content in the brain can be obtained by the multi-echo [Formula: see text] weighted images ( [Formula: see text]WIs). To accelerate the long acquisition, a novel tensor dictionary learning algorithm with low-rank and sparse regularization (TDLLS) is proposed to reconstruct the [Formula: see text]WIs from the undersampled data. The proposed algorithm explores the local and nonlocal similarity and the global temporal redundancy in the real and imaginary parts of the complex relaxation signals. The joint application of the low-rank constraints on the dictionaries and the sparse constraints on the core coefficient tensors improves the performance of the tensor-based recovery. Parallel imaging is incorporated into the TDLLS algorithm (pTDLLS) for further acceleration. A pulse sequence is proposed to prospectively undersample the Ky-t space to obtain the whole brain high-quality myelin water fraction (MWF) maps within 1 minute at an undersampling rate (R) of 6.

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.

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.

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.

DeepResp: Deep learning solution for respiration-induced B0 fluctuation artifacts in multi-slice GRE

Respiration-induced B₀ fluctuation corrupts MRI images by inducing phase errors in k-space. A few approaches such as navigator have been proposed to correct for the artifacts at the expense of sequence modification. In this study, a new deep learning method, which is referred to as DeepResp, is proposed for reducing the respiration-artifacts in multi-slice gradient echo (GRE) images. DeepResp is designed to extract the respiration-induced phase errors from a complex image using deep neural networks. Then, the network-generated phase errors are applied to the k-space data, creating an artifact-corrected image. For network training, the computer-simulated images were generated using artifact-free images and respiration data. When evaluated, both simulated images and in-vivo images of two different breathing conditions (deep breathing and natural breathing) show improvements (simulation: normalized root-mean-square error (NRMSE) from 7.8 ± 5.2% to 1.3 ± 0.6%; structural similarity (SSIM) from 0.88 ± 0.08 to 0.99 ± 0.01; ghost-to-signal-ratio (GSR) from 7.9 ± 7.2% to 0.6 ± 0.6%; deep breathing: NRMSE from 13.9 ± 4.6% to 5.8 ± 1.4%; SSIM from 0.86 ± 0.03 to 0.95 ± 0.01; GSR 20.2 ± 10.2% to 5.7 ± 2.3%; natural breathing: NRMSE from 5.2 ± 3.3% to 4.0 ± 2.5%; SSIM from 0.94 ± 0.04 to 0.97 ± 0.02; GSR 5.7 ± 5.0% to 2.8 ± 1.1%). Our approach does not require any modification of the sequence or additional hardware, and may therefore find useful applications. Furthermore, the deep neural networks extract respiration-induced phase errors, which is more interpretable and reliable than results of end-to-end trained networks.

Probing myelin content of the human brain with MRI: A review

Rapid and efficient transmission of electric signals among neurons of vertebrates is ensured by myelin-insulating sheaths surrounding axons. Human cognition, sensation, and motor functions rely on the integrity of these layers, and demyelinating diseases often entail serious cognitive and physical impairments. Magnetic resonance imaging radically transformed the way these disorders are monitored, offering an irreplaceable tool to noninvasively examine the brain structure. Several advanced techniques based on MRI have been developed to provide myelin-specific contrasts and a quantitative estimation of myelin density in vivo. Here, the vast offer of acquisition strategies developed to date for this task is reviewed. Advantages and pitfalls of the different approaches are compared and discussed.

Artificial neural network for multi-echo gradient echo-based myelin water fraction estimation

PURPOSE

To demonstrate robust myelin water fraction (MWF) mapping using an artificial neural network (ANN) with multi-echo gradient-echo (GRE) signal.

METHODS

Multi-echo gradient-echo signals simulated with a three-pool exponential model were used to generate the training data set for the ANN, which was designed to yield the MWF. We investigated the performance of our proposed ANN for various conditions using both numerical simulations and in vivo data. Simulations were conducted with various SNRs to investigate the performance of the ANN. In vivo data with high spatial resolutions were applied in the analyses, and results were compared with MWFs derived by the nonlinear least-squares algorithm using a complex three-pool exponential model.

RESULTS

The network results for the simulations show high accuracies against noise compared with nonlinear least-squares MWFs: RMS-error value of 5.46 for the nonlinear least-squares MWF and 3.56 for the ANN MWF at an SNR of 150 (relative gain = 34.80%). These effects were also found in the in vivo data, with reduced SDs in the region-of-interest analyses. These effects of the ANN demonstrate the feasibility of acquiring high-resolution myelin water images.

CONCLUSION

The simulation results and in vivo data suggest that the ANN facilitates more robust MWF mapping in multi-echo gradient-echo sequences compared with the conventional nonlinear least-squares method.

Blind Source Separation for Myelin Water Fraction Mapping Using Multi-Echo Gradient Echo Imaging

In conventional gradient-echo myelin water imaging (GRE-MWI), myelin water fraction (MWF) is estimated by fitting the multi-echo gradient recalled echo (mGRE) signal to a pre-assumed numerical model (e.g., multi-component exponential curves or three component exponential curves). However, in mGRE, imaging artifacts (e.g., voxel spread function and physiological noise) and noise render the signal to deviate from the numerical model, leading to misfit of the model parameters. Here, as an alternative to the model-based GRE-MWI, a blind source separation (BSS) technique for the separation of multi-exponential mGRE signal is proposed. Among the various BSS techniques, a modified robust principal component analysis (rPCA) is presented to separate signal sources by enforcing the data-driven properties such as "low rankness" and "sparsity." Considering the signal evolution of T₂^∗ relaxation (i.e., non-negative exponential decay), low rankness of exponential decay was enforced by nonnegative matrix factorization (NMF) and hankelization. This method provides the separation of slow-decaying, fast-decaying exponential components and artifact components from mGRE images. After the separation, MWF map is reconstructed as the ratio of the fast-decaying component to the total decaying components. The proposed method was demonstrated in numerical simulations and in vivo scans. The method provided a robust estimation of MWF in the presence of statistical noise and imaging artifacts.

So You Want to Image Myelin Using MRI: An Overview and Practical Guide for Myelin Water Imaging

Myelin water imaging (MWI) is an MRI imaging biomarker for myelin. This method can generate an in vivo whole-brain myelin water fraction map in approximately 10 minutes. It has been applied in various applications including neurodegenerative disease, neurodevelopmental, and neuroplasticity studies. In this review we start with a brief introduction of myelin biology and discuss the contributions of myelin in conventional MRI contrasts. Then the MRI properties of myelin water and four different MWI methods, which are categorized as T₂ -, T₂ *-, T₁ -, and steady-state-based MWI, are summarized. After that, we cover more practical issues such as availability, interpretation, and validation of these methods. To illustrate the utility of MWI as a clinical research tool, MWI studies for two diseases, multiple sclerosis and neuromyelitis optica, are introduced. Additional topics about imaging myelin in gray matter and non-MWI methods for myelin imaging are also included. Although technical and physiological limitations exist, MWI is a potent surrogate biomarker of myelin that carries valuable and useful information of myelin. Evidence Level: 5 Technical Efficacy: 1 J. MAGN. RESON. IMAGING 2021;53:360-373.

Assessment and correction of macroscopic field variations in 2D spoiled gradient-echo sequences

Purpose

To model and correct the dephasing effects in the gradient‐echo signal for arbitrary RF excitation pulses with large flip angles in the presence of macroscopic field variations.

Methods

The dephasing of the spoiled 2D gradient‐echo signal was modeled using a numerical solution of the Bloch equations to calculate the magnitude and phase of the transverse magnetization across the slice profile. Additionally, regional variations of the transmit RF field and slice profile scaling due to macroscopic field gradients were included. Simulations, phantom, and in vivo measurements at 3 T were conducted for R2∗ and myelin water fraction (MWF) mapping.

Results

The influence of macroscopic field gradients on R2∗ and myelin water fraction estimation can be substantially reduced by applying the proposed model. Moreover, it was shown that the dephasing over time for flip angles of 60° or greater also depends on the polarity of the slice‐selection gradient because of phase variation along the slice profile.

Conclusion

Substantial improvements in R2∗ accuracy and myelin water fraction mapping coverage can be achieved using the proposed model if higher flip angles are required. In this context, we demonstrated that the phase along the slice profile and the polarity of the slice‐selection gradient are essential for proper modeling of the gradient‐echo signal in the presence of macroscopic field variations.

Improvement of reproducibility in quantitative susceptibility mapping (QSM) and transverse relaxation rates ( R 2 * ) after physiological noise correction

BACKGROUND

Numerous studies have suggested that quantitative susceptibility mapping (QSM) and transverse relaxation rates ( R 2 * ) are useful to monitor neurological diseases. For clinical use of QSM and R 2 * , reproducibility is an important feature. However, respiration-induced local magnetic field variation makes artifacts in gradient echo-based images and reduces the reproducibility of QSM and R 2 * .

PURPOSE

To investigate the improvement of reproducibility of QSM and R 2 * after the correction of respiration-induced field variation, and to assess the effect of varying types of the region of interest (ROI) analysis on reproducibility.

STUDY TYPE

Reproducibility study.

POPULATION

Ten controls.

FIELD STRENGTH/SEQUENCE

3T/multiecho gradient echo sequence.

ASSESSMENT

Intrascan reproducibility of QSM and R 2 * was investigated in ROIs before and after the respiration correction.

STATISTICAL TESTS

Reproducibility was obtained by the square of voxel-wise correlation coefficients between scans. A paired t-test was performed for comparison between before and after the respiration correction and between QSM and R 2 * .

RESULTS

Based on the ROI analysis, reproducibility increased after the respiration correction. Reproducibility in the white matter (11.89% increased in QSM and 23.38% in R 2 * , P = 0.009 and 0.024, respectively) and deep gray matter (5.50% increased in QSM and 13.96% in R 2 * , P = 0.024 and 0.019, respectively) increased significantly after the respiration correction. Reproducibility of R 2 * was higher than that of QSM in the whole brain and cortical gray matter, while QSM maps showed higher reproducibility than R 2 * in the white matter and deep gray matter.

DATA CONCLUSION

Respiration-induced error correction significantly improved reproducibility in QSM and R 2 * mapping. QSM and R 2 * mapping showed a different level of reproducibility depending on the types of ROI analysis.

LEVEL OF EVIDENCE

4 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2018.

Whole brain g-ratio mapping using myelin water imaging (MWI) and neurite orientation dispersion and density imaging (NODDI)

MR g-ratio, which measures the ratio of the aggregate volume of axons to that of fibers in a voxel, is a potential biomarker for white matter microstructures. In this study, a new approach for acquiring an in-vivo whole human brain g-ratio map is proposed. To estimate the g-ratio, myelin volume fraction and axonal volume fraction are acquired using multi-echo gradient echo myelin water imaging (GRE-MWI) and neurite orientation dispersion and density imaging (NODDI), respectively. In order to translate myelin water fraction measured in GRE-MWI into myelin volume fraction, a new scaling procedure is proposed and validated. This scaling approach utilizes geometric measures of myelin structure and, therefore, provides robustness over previous methods. The resulting g-ratio map reveals an expected range of g-ratios (0.71-0.85 in major fiber bundles) with a small inter-subject coefficient of variance (less than 2%). Additionally, a few fiber bundles (e.g. cortico-spinal tract and optic radiation) show different constituents of myelin volume fraction and axonal volume fraction, indicating potentials to utilize the measures for deciphering fiber tracking.

High resolution myelin water imaging incorporating local tissue susceptibility analysis

Quantitative myelin water imaging (MWI) from signal T₂* decay acquired with multiple Gradient-Recalled Echo (mGRE) sequence has been widely used since its first report. A recent study showed that with low resolution data (2mm isotropic voxels), direct application of complex fitting to a three-pool WM model with frequency shift terms could produce more stable parameter estimation for myelin water fraction mapping. MWI maps of higher spatial resolution resulting in more detailed tissue structures and reduced partial volume effects around white matter/gray matter (WM/GM) interface, however, is more desirable. Furthermore, as signal-to-noise ratio (SNR) of original images decreases due to reduced voxel size, the direct complex fitting procedure of myelin water imaging becomes more prone to systematic errors which severely compromised stability and reliability of the result. Instead of using the original part of T₂* decay, this work presents a new method based on the WM-induced phase from tissue susceptibility calculated with the same mGRE dataset, in a three-pool WM model (water of myelin, axonal and extracellular water), to improve high resolution MWI. Compared with direct complex fitting for the higher spatial resolution case, the proposed method is shown to provide a more stable and accurate estimation of MWI parameters, and finer details near WM/GM boundaries with greatly reduced partial volume effects.

Mechanisms of T2 * anisotropy and gradient echo myelin water imaging

In MRI, structurally aligned molecular or micro-organization (e.g. axonal fibers) can be a source of substantial signal variations that depend on the structural orientation and the applied magnetic field. This signal anisotropy gives us a unique opportunity to explore information that exists at a resolution several orders of magnitude smaller than that of typical MRI. In this review, one of the signal anisotropies, T₂ * anisotropy in white matter, and a related imaging method, gradient echo myelin water imaging (GRE-MWI), are explored. The T₂ * anisotropy has been attributed to isotropic and anisotropic magnetic susceptibility of myelin and compartmentalized microstructure of white matter fibers (i.e. axonal, myelin, and extracellular space). The susceptibility and microstructure create magnetic frequency shifts that change with the relative orientation of the fiber and the main magnetic field, generating the T₂ * anisotropy. The resulting multi-component magnitude decay and nonlinear phase evolution have been utilized for GRE-MWI, assisting in resolving the signal fraction of the multiple compartments in white matter. The GRE-MWI method has been further improved by signal compensation techniques including physiological noise compensation schemes. The T₂ * anisotropy and GRE-MWI provide microstructural information on a voxel (e.g. fiber orientation and tissue composition), and may serve as sensitive biomarkers for microstructural changes in the brain. Copyright © 2016 John Wiley & Sons, Ltd.

USPIO enhanced lymph node MRI using 3D multi-echo GRE in a rabbit model

Ultrasmallsuperparamagnetic iron oxide (USPIO) has been suggested to be a negative MR contrast agent to detect metastatic lymph nodes. Previously reported studies have evaluated the diagnostic performance of USPIO-enhanced MR lymph node imaging based on signal intensity. In this study, we investigate the specific performance of three different parametric approaches (normalized signal intensity, R₂ * and susceptibility) using 3D multi-echo gradient echo to quantify the USPIO particles in lymph nodes. Nine rabbits with VX2 tumor implants were scanned before and after USPIO injection. From 3D multi-echo GRE magnitude and phase data, we generated multi-echo combined T₂ *-weighted images, an R₂ * map, and a quantitative susceptibility map. Eighteen lymph nodes (nine reactive and nine metastatic) were evaluated and showed remarkable signal drops in the area of USPIO accumulation. On parametric analysis, the R₂ * difference before and after USPIO injection was significantly different (p < 0.05) between reactive and metastatic lymph nodes; in contrast, the normalized signal intensity and susceptibility were not significantly different between the nodes. Our study showed the potential utility of USPIO-enhanced MRI using R2* mapping from 3D multi-echo GRE for the detection of lymph node metastasis and parametric analysis of lymph node status in a rabbit model. Copyright © 2016 John Wiley & Sons, Ltd.

Imaging of nigrosome 1 in substantia nigra at 3T using multiecho susceptibility map-weighted imaging (SMWI)

PURPOSE

To enhance the visibility of nigrosome 1 in substantia nigra, which has recently been suggested as an imaging biomarker for Parkinson's disease (PD) at 3T magnetic resonance imaging (MRI).

MATERIALS AND METHODS

The substantia nigra structure was visualized at 3T MRI using multiecho susceptibility map-weighted imaging (SMWI) in 15 healthy volunteers and 6 patients with Parkinson's disease (PD). The visibility of nigrosome 1 was further enhanced by acquiring data in an oblique-coronal imaging plane at a high spatial resolution (0.5 × 0.5 × 1.0 mm³ ). To compare the visibility, the contrast-to-noise ratios (CNR) of the nigrosome 1 structure relative to the neighboring substantia nigra structure were evaluated in the SMWI and other conventional susceptibility contrast images (magnitude, frequency, quantitative susceptibility map [QSM] and susceptibility-weighted image).

RESULTS

In healthy volunteers, the CNRs of the nigrosome 1 structure were 1.04 ± 0.38, 0.84 ± 0.32, 1.04 ± 0.40, 0.86 ± 0.41, and 1.45 ± 0.48 for magnitude, frequency, quantitative susceptibility map, susceptibility-weighted image, and SMWI, respectively. Compared to conventional susceptibility contrast images, the SMWI method significantly improved the CNR of nigrosome 1 (P = 0.014 for magnitude, P = 0.030 for QSM, and P < 0.001 for frequency and SWI, respectively). The magnetic susceptibility difference between nigrosome 1 and neighboring substantia nigra structures was 0.037 ± 0.016 ppm (measured in QSM, P < 0.001) in healthy volunteers. In the PD patients, the visibility of the nigrosome 1 structures was reduced.

CONCLUSION

The SMWI method enhances the visibility of nigrosome 1 structures at 3T MRI when compared to conventional susceptibility contrast images.

LEVEL OF EVIDENCE

3 Technical Efficacy: Stage 2 J. MAGN. RESON. IMAGING 2017;46:528-536.

Rapid simultaneous high-resolution mapping of myelin water fraction and relaxation times in human brain using BMC-mcDESPOT

A number of central nervous system (CNS) diseases exhibit changes in myelin content and magnetic resonance longitudinal, T₁, and transverse, T₂, relaxation times, which therefore represent important biomarkers of CNS pathology. Among the methods applied for measurement of myelin water fraction (MWF) and relaxation times, the multicomponent driven equilibrium single pulse observation of T₁ and T₂ (mcDESPOT) approach is of particular interest. mcDESPOT permits whole brain mapping of multicomponent T₁ and T₂, with data acquisition accomplished within a clinically realistic acquisition time. Unfortunately, previous studies have indicated the limited performance of mcDESPOT in the setting of the modest signal-to-noise range of high-resolution mapping, required for the depiction of small structures and to reduce partial volume effects. Recently, we showed that a new Bayesian Monte Carlo (BMC) analysis substantially improved determination of MWF from mcDESPOT imaging data. However, our previous study was limited in that it did not discuss determination of relaxation times. Here, we extend the BMC analysis to the simultaneous determination of whole-brain MWF and relaxation times using the two-component mcDESPOT signal model. Simulation analyses and in-vivo human brain studies indicate the overall greater performance of this approach compared to the stochastic region contraction (SRC) algorithm, conventionally used to derive parameter estimates from mcDESPOT data. SRC estimates of the transverse relaxation time of the long T₂ fraction, T_2,l, and the longitudinal relaxation time of the short T₁ fraction, T_1,s, clustered towards the lower and upper parameter search space limits, respectively, indicating failure of the fitting procedure. We demonstrate that this effect is absent in the BMC analysis. Our results also showed improved parameter estimation for BMC as compared to SRC for high-resolution mapping. Overall we find that the combination of BMC analysis and mcDESPOT, BMC-mcDESPOT, shows excellent performance for accurate high-resolution whole-brain mapping of MWF and bi-component transverse and longitudinal relaxation times within a clinically realistic acquisition time.