Improved estimation of myelin water fraction using complex model fitting

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.

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.

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.

Conversion map from quantitative parameter mapping to myelin water fraction: comparison with R1·R2* and myelin water fraction in white matter

OBJECTIVE

To clarify the relationship between myelin water fraction (MWF) and R1⋅R2* and to develop a method to calculate MWF directly from parameters derived from QPM, i.e., MWF converted from QPM (MWFQPM).

MATERIALS AND METHODS

Subjects were 12 healthy volunteers. On a 3 T MR scanner, dataset was acquired using spoiled gradient-echo sequence for QPM. MWF and R1⋅R2* maps were derived from the multi-gradient-echo (mGRE) dataset. Volume-of-interest (VOI) analysis using the JHU-white matter (WM) atlas was performed. All the data in the 48 WM regions measured VOI were plotted, and quadratic polynomial approximations of each region were derived from the relationship between R1·R2* and the two-pool model-MWF. The R1·R2* map was converted to MWFQPM map. MWF atlas template was generated using converted to MWF from R1·R2* per WM region.

RESULTS

The mean MWF and R1·R2* values for the 48 WM regions were 11.96 ± 6.63%, and 19.94 ± 4.59 s-2, respectively. A non-linear relationship in 48 regions of the WM between MWF and R1·R2* values was observed by quadratic polynomial approximation (R2 ≥ 0.963, P < 0.0001).

DISCUSSION

MWFQPM map improved image quality compared to the mGRE-MWF map. Myelin water atlas template derived from MWFQPM may be generated with combined multiple WM regions.

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.

Myelin water imaging at 0.55 T using a multigradient-echo sequence

PURPOSE

To investigate the prospects of a multigradient-echo (mGRE) acquisition for in vivo myelin water imaging at 0.55 T.

METHODS

Scans were performed on the brain of four healthy volunteers at 0.55 and 3 T, using a 3D mGRE sequence. The myelin water fraction (MWF) was calculated for both field strengths using a nonnegative least squares (NNLS) algorithm, implemented in the qMRLab suite. The quality of these maps as well as single-voxel fits were compared visually for 0.55 and 3 T.

RESULTS

The obtained MWF values at 0.55 T are consistent with previously reported ones at higher field strengths. The MWF maps are a considerable improvement over the ones at 3 T. Example fits show that 0.55 T data is better described by an exponential model than 3 T data, making the assumed multi-exponential model of the NNLS algorithm more accurate.

CONCLUSION

This first assessment shows that mGRE myelin water imaging at 0.55 T is feasible and has the potential to yield better results than at higher fields.

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.

Quantitative myelin water imaging using short TR adiabatic inversion recovery prepared echo-planar imaging (STAIR-EPI) sequence

Introduction

Numerous techniques for myelin water imaging (MWI) have been devised to specifically assess alterations in myelin. The biomarker employed to measure changes in myelin content is known as the myelin water fraction (MWF). The short TR adiabatic inversion recovery (STAIR) sequence has recently been identified as a highly effective method for calculating MWF. The purpose of this study is to develop a new clinical transitional myelin water imaging (MWI) technique that combines STAIR preparation and echo-planar imaging (EPI) (STAIR-EPI) sequence for data acquisition.

Methods

Myelin water (MW) in the brain has shorter T1 and T2 relaxation times than intracellular and extracellular water. In the proposed STAIR-EPI sequence, a short TR (e.g., ≤300 ms) together with an optimized inversion time enable robust long T1 water suppression with a wide range of T1 values [i.e., (600, 2,000) ms]. The EPI allows fast data acquisition of the remaining MW signals. Seven healthy volunteers and seven patients with multiple sclerosis (MS) were recruited and scanned in this study. The apparent myelin water fraction (aMWF), defined as the signal ratio of MW to total water, was measured in the lesions and normal-appearing white matter (NAWM) in MS patients and compared with those measured in the normal white matter (NWM) in healthy volunteers.

Results

As seen in the STAIR-EPI images acquired from MS patients, the MS lesions show lower signal intensities than NAWM do. The aMWF measurements for both MS lesions (3.6 ± 1.3%) and NAWM (8.6 ± 1.2%) in MS patients are significantly lower than NWM (10 ± 1.3%) in healthy volunteers (P < 0.001).

Discussion

The proposed STAIR-EPI technique, which can be implemented in MRI scanners from all vendors, is able to detect myelin loss in both MS lesions and NAWM in MS patients.

Microstructure-Informed Myelin Mapping (MIMM) from Gradient Echo MRI using Stochastic Matching Pursuit

Quantification of the myelin content of the white matter is important for studying demyelination in neurodegenerative diseases such as Multiple Sclerosis (MS), particularly for longitudinal monitoring. A novel noninvasive MRI method, called Microstructure-Informed Myelin Mapping (MIMM), is developed to quantify the myelin volume fraction (MVF) by utilizing a multi gradient echo sequence (mGRE) and a detailed biophysical model of tissue microstructure. Myelin is modeled as anisotropic negative susceptibility source based on the Hollow Cylindrical Fiber Model (HCFM), and iron as isotropic positive susceptibility source in the extracellular region. Voxels with a range of biophysical parameters are simulated to create a dictionary of MR echo time magnitude signals and total susceptibility values. MRI signals measured using a mGRE sequence are then matched voxel-by-voxel to the created dictionary to obtain the spatial distributions of myelin and iron. Three different MIMM versions are presented to deal with the fiber orientation dependent susceptibility effects of the myelin sheaths: a basic variation, which assumes fiber orientation is an unknown to fit, two orientation informed variations, which assume the fiber orientation distribution is available either from a separate diffusion tensor imaging (DTI) acquisition or from a DTI atlas based fiber orientation map. While all showed a significant linear correlation with the reference method based on T2-relaxometry (p < 0.0001), DTI orientation informed and atlas orientation informed variations reduced overestimation at white matter tracts compared to the basic variation. Finally, the implications and usefulness of attaining an additional iron susceptibility distribution map are discussed.

Highlights

novel stochastic matching pursuit algorithm called microstructure-informed myelin mapping (MIMM) is developed to quantify Myelin Volume Fraction (MVF) using Magnetic Resonance Imaging (MRI) and microstructural modeling.utilizes a detailed biophysical model to capture the susceptibility effects on both magnitude and phase to quantify myelin and iron.matter fiber orientation effects are considered for the improved MVF quantification in the major fiber tracts.acquired myelin and iron maps may be utilized to monitor longitudinal disease progress.

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.

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.

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.

Self-labelled encoder-decoder (SLED) for multi-echo gradient echo-based myelin water imaging

PURPOSE

Reconstruction of high quality myelin water imaging (MWI) maps is challenging, particularly for data acquired using multi-echo gradient echo (mGRE) sequences. A non-linear least squares fitting (NLLS) approach has often been applied for MWI. However, this approach may produce maps with limited detail and, in some cases, sub-optimal signal to noise ratio (SNR), due to the nature of the voxel-wise fitting. In this study, we developed a novel, unsupervised learning method called self-labelled encoder-decoder (SLED) to improve gradient echo-based MWI data fitting.

METHODS

Ultra-high resolution, MWI data was collected from five mouse brains with variable levels of myelination, using a mGRE sequence. Imaging data was acquired using a 7T preclinical MRI system. A self-labelled, encoder-decoder network was implemented in TensorFlow for calculation of myelin water fraction (MWF) based on the mGRE signal decay. A simulated MWI phantom was also created to evaluate the performance of MWF estimation.

RESULTS

Compared to NLLS, SLED demonstrated improved MWF estimation, in terms of both stability and accuracy in phantom tests. In addition, SLED produced less noisy MWF maps from high resolution MR microscopy images of mouse brain tissue. It specifically resulted in lower noise amplification for all mouse genotypes that were imaged and yielded mean MWF values in white matter ROIs that were highly correlated with those derived from standard NLLS fitting. Lastly, SLED also exhibited higher tolerance to low SNR data.

CONCLUSION

Due to its unsupervised and self-labeling nature, SLED offers a unique alternative to analyze gradient echo-based MWI data, providing accurate and stable MWF estimations.

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.

Myelin water imaging using a short-TR adiabatic inversion-recovery (STAIR) sequence

PURPOSE

To develop a new myelin water imaging (MWI) technique using a short-TR adiabatic inversion-recovery (STAIR) sequence on a clinical 3T MR scanner.

METHODS

Myelin water (MW) in the brain has both a much shorter T₁ and a much shorter T₂ * than intracellular/extracellular water. A STAIR sequence with a short TR was designed to efficiently suppress long T₁ signals from intracellular/extracellular water, and therefore allow selective imaging of MW, which has a much shorter T₁ . Numerical simulation and phantom studies were performed to investigate the effectiveness of long T₁ signal suppression. TheT₂ * in white matter (WM) was measured with STAIR and compared with T₂ * measured with a conventional gradient recall echo in in vivo study. Four healthy volunteers and 4 patients with multiple sclerosis were recruited for qualitative and quantitative MWI. Apparent MW fraction was generated to compare MW in normal WM in volunteers to MW in lesions in patients with multiple sclerosis.

RESULTS

Both simulation and phantom studies showed that when TR was sufficiently short (eg, 250 ms), the STAIR sequence effectively suppressed long T₁ signals from tissues with a broad range of T₁ s using a single TR/TI combination. The volunteer study showed a short T₂ * of 9.5 ± 1.7 ms in WM, which is similar to reported values for MW. Lesions in patients with multiple sclerosis showed a significantly lower apparent MW fraction (4.5% ± 1.0%) compared with that of normal WM (9.2% ± 1.5%) in healthy volunteers (p < 0.05).

CONCLUSIONS

The STAIR sequence provides selective MWI in brain and can quantify reductions in MW content in patients with multiple sclerosis.

Mapping the human connectome using diffusion MRI at 300 mT/m gradient strength: Methodological advances and scientific impact

Tremendous efforts have been made in the last decade to advance cutting-edge MRI technology in pursuit of mapping structural connectivity in the living human brain with unprecedented sensitivity and speed. The first Connectom 3T MRI scanner equipped with a 300 mT/m whole-body gradient system was installed at the Massachusetts General Hospital in 2011 and was specifically constructed as part of the Human Connectome Project. Since that time, numerous technological advances have been made to enable the broader use of the Connectom high gradient system for diffusion tractography and tissue microstructure studies and leverage its unique advantages and sensitivity to resolving macroscopic and microscopic structural information in neural tissue for clinical and neuroscientific studies. The goal of this review article is to summarize the technical developments that have emerged in the last decade to support and promote large-scale and scientific studies of the human brain using the Connectom scanner. We provide a brief historical perspective on the development of Connectom gradient technology and the efforts that led to the installation of three other Connectom 3T MRI scanners worldwide - one in the United Kingdom in Cardiff, Wales, another in continental Europe in Leipzig, Germany, and the latest in Asia in Shanghai, China. We summarize the key developments in gradient hardware and image acquisition technology that have formed the backbone of Connectom-related research efforts, including the rich array of high-sensitivity receiver coils, pulse sequences, image artifact correction strategies and data preprocessing methods needed to optimize the quality of high-gradient strength diffusion MRI data for subsequent analyses. Finally, we review the scientific impact of the Connectom MRI scanner, including advances in diffusion tractography, tissue microstructural imaging, ex vivo validation, and clinical investigations that have been enabled by Connectom technology. We conclude with brief insights into the unique value of strong gradients for diffusion MRI and where the field is headed in the coming years.

Data-driven myelin water imaging based on T1 and T2 relaxometry

Long acquisition times preclude the application of multiecho spin echo (MESE) sequences for myelin water fraction (MWF) mapping in daily clinical practice. In search of alternative methods, previous studies of interest explored the biophysical modeling of MWF from measurements of different tissue properties that can be obtained in scan times shorter than those required for the MESE. In this work, a novel data-driven estimation of MWF maps from fast relaxometry measurements is proposed and investigated. T₁ and T₂ relaxometry maps were acquired in a cohort of 20 healthy subjects along with a conventional MESE sequence. Whole-brain quantitative mapping was achieved with a fast protocol in 6 min 24 s. Reference MWF maps were derived from the MESE sequence (TA = 11 min 17 s) and their data-driven estimation from relaxometry measurements was investigated using three different modeling strategies: two general linear models (GLMs) with linear and quadratic regressors, respectively; a random forest regression model; and two deep neural network architectures, a U-Net and a conditional generative adversarial network (cGAN). Models were validated using a 10-fold crossvalidation. The resulting maps were visually and quantitatively compared by computing the root mean squared error (RMSE) between the estimated and reference MWF maps, the intraclass correlation coefficients (ICCs) between corresponding MWF values in different brain regions, and by performing Bland-Altman analysis. Qualitatively, the estimated maps appear to generally provide a similar, yet more blurred MWF contrast in comparison with the reference, with the cGAN model best capturing MWF variabilities in small structures. By estimating the average adjusted coefficient of determination of the GLM with quadratic regressors, we showed that 87% of the variability in the MWF values can be explained by relaxation times alone. Further quantitative analysis showed an average RMSE smaller than 0.1% for all methods. The ICC was greater than 0.81 for all methods, and the bias smaller than 2.19%. It was concluded that this work confirms the notion that relaxometry parameters contain a large part of the information on myelin water and that MWF maps can be generated from T₁ /T₂ data with minimal error. Among the investigated modeling approaches, the cGAN provided maps with the best trade-off between accuracy and blurriness. Fast relaxometry, like the 6 min 24 s whole-brain protocol used in this work in conjunction with machine learning, may thus have the potential to replace time-consuming MESE acquisitions.

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.

Susceptibility source separation from gradient echo data using magnitude decay modeling

BACKGROUND AND PURPOSE

The objective is to demonstrate feasibility of separating magnetic sources in quantitative susceptibility mapping (QSM) by incorporating magnitude decay rates R 2 ∗ $R_2^{\rm{*}}$ in gradient echo (GRE) MRI.

METHODS

Magnetic susceptibility source separation was developed using R 2 ∗ $R_2^{\rm{*}}$ and compared with a prior method using R 2 ' = R 2 ∗ - R 2 ${R^{\prime}_2} = R_2^* - {R_2}$ that required an additional sequence to measure the transverse relaxation rate R₂ . Both susceptibility separation methods were compared in multiple sclerosis (MS) patients (n = 17). Susceptibility values of negative sources estimated with R 2 ∗ $R_2^{\rm{*}}$ -based source separation in a set of enhancing MS lesions (n = 44) were correlated against longitudinal myelin water fraction (MWF) changes.

RESULTS

In in vivo data, linear regression of the estimated χ + ${\chi}^{+}$ and χ - ${\chi}^{-}$ susceptibility values between the R 2 ∗ $R_2^*$ - and the R 2 ' ${R^{\prime}_2}$ -based separation methods performed across 182 segmented lesions revealed correlation coefficient r = .96 and slope close .99. Correlation analysis in enhancing lesions revealed a significant positive association between the χ - ${\chi}^{-}$ increase at 1-year post-onset relative to 0 year and the MWF increase at 1 year relative to 0 year (β = -0.144, 95% confidence interval: [-0.199, -0.1], p = .0008) and good agreement between R 2 ' ${R^{\prime}_2}$ and R 2 ∗ $R_2^*$ methods (r = .79, slope = .95).

CONCLUSIONS

Separation of magnetic sources based solely on GRE complex data is feasible by combining magnitude decay rate modeling and phase-based QSM and χ - ${\chi}^{-}$ change may serve as a biomarker for myelin recovery or damage in acute MS lesions.

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.

Insights into human cerebral white matter maturation and degeneration across the adult lifespan

White matter (WM) microstructural properties change across the adult lifespan and with neuronal diseases. Understanding microstructural changes due to aging is paramount to distinguish them from neuropathological changes. Conducted on a large cohort of 147 cognitively unimpaired subjects, spanning a wide age range of 21 to 94 years, our study evaluated sex- and age-related differences in WM microstructure. Specifically, we used diffusion tensor imaging (DTI) magnetic resonance imaging (MRI) indices, sensitive measures of myelin and axonal density in WM, and myelin water fraction (MWF), a measure of the fraction of the signal of water trapped within the myelin sheets, to probe these differences. Furthermore, we examined regional correlations between MWF and DTI indices to evaluate whether the DTI metrics provide information complementary to MWF. While sexual dimorphism was, overall, nonsignificant, we observed region-dependent differences in MWF, that is, myelin content, and axonal density with age and found that both exhibit nonlinear, but distinct, associations with age. Furthermore, DTI indices were moderately correlated with MWF, indicating their good sensitivity to myelin content as well as to other constituents of WM tissue such as axonal density. The microstructural differences captured by our MRI metrics, along with their weak to moderate associations with MWF, strongly indicate the potential value of combining these outcome measures in a multiparametric approach. Furthermore, our results support the last-in-first-out and the gain-predicts-loss hypotheses of WM maturation and degeneration. Indeed, our results indicate that the posterior WM regions are spared from neurodegeneration as compared to anterior regions, while WM myelination follows a temporally symmetric time course across the adult life span.

Subsecond accurate myelin water fraction reconstruction from FAST-T2 data with 3D UNET

PURPOSE

To develop a 3D UNET convolutional neural network for rapid extraction of myelin water fraction (MWF) maps from six-echo fast acquisition with spiral trajectory and T₂ -prep data and to evaluate its accuracy in comparison with multilayer perceptron (MLP) network.

METHODS

The MWF maps were extracted from 138 patients with multiple sclerosis using an iterative three-pool nonlinear least-squares algorithm (NLLS) without and with spatial regularization (srNLLS), which were used as ground-truth labels to train, validate, and test UNET and MLP networks as a means to accelerate data fitting. Network testing was performed in 63 patients with multiple sclerosis and a numerically simulated brain phantom at SNR of 200, 100 and 50.

RESULTS

Simulations showed that UNET reduced the MWF mean absolute error by 30.1% to 56.4% and 16.8% to 53.6% over the whole brain and by 41.2% to 54.4% and 21.4% to 49.4% over the lesions for predicting srNLLS and NLLS MWF, respectively, compared to MLP, with better performance at lower SNRs. UNET also outperformed MLP for predicting srNLLS MWF in the in vivo multiple-sclerosis brain data, reducing mean absolute error over the whole brain by 61.9% and over the lesions by 67.5%. However, MLP yielded 41.1% and 51.7% lower mean absolute error for predicting in vivo NLLS MWF over the whole brain and the lesions, respectively, compared with UNET. The whole-brain MWF processing time using a GPU was 0.64 seconds for UNET and 0.74 seconds for MLP.

CONCLUSION

Subsecond whole-brain MWF extraction from fast acquisition with spiral trajectory and T₂ -prep data using UNET is feasible and provides better accuracy than MLP for predicting MWF output of srNLLS algorithm.

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.

Brain myelin water fraction is associated with APOE4 allele status in patients with cognitive impairment

BACKGROUND AND PURPOSE

Apolipoprotein E4 (APOE4) is a major genetic risk factor for Alzheimer's disease. However, the effect of APOE4 status on myelin remains unclear. This study investigated the effect of APOE4 on myelin content in cognitively impaired individuals using T2* gradient echo (GRE)-based myelin water fraction (MWF) imaging.

METHODS

Between August 2017 and January 2019, we evaluated 39 cognitively impaired patients (median age, 75 years; male:female = 8:31; Alzheimer's disease: mild cognitive impairment = 11:28). We obtained brain MWF values from white matter hyperintensities (WMHs) and normal-appearing white matter (NAWM). Linear regression analysis was performed to investigate the relationship between the APOE4 status and MWF and cognitive function and MWF.

RESULTS

Among the 39 cognitively impaired patients, nine (23.1%) were APOE4 carriers and 30 (76.9%) were noncarriers. APOE4 carriers had a lower hippocampal volume than noncarriers (p = .045), but other brain volume parameters were not differed. After age adjustment, the APOE4 status was significantly associated with reduced MWF in NAWM (β = -0.310 per allele; p = .049) but not in WMH (β = -0.258 per allele; p = .113). After age adjustment, MWF in NAWM was significantly associated with Mini-Mental State Examination score (β = 0.313, p = .031).

CONCLUSIONS

T2* GRE-based MWF imaging can reveal myelin loss, particularly in NAWM, in cognitively impaired patients among APOE4 carriers. In vivo MWF in NAWM might be a novel imaging marker of Alzheimer's disease, for clarifying the interactions between the white matter and cognitive dysfunction with respect to the APOE4 status.

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.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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

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.

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.

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.

Comparison of non-parametric T2 relaxometry methods for myelin water quantification

Multi-component T₂ relaxometry allows probing tissue microstructure by assessing compartment-specific T₂ relaxation times and water fractions, including the myelin water fraction. Non-negative least squares (NNLS) with zero-order Tikhonov regularization is the conventional method for estimating smooth T₂ distributions. Despite the improved estimation provided by this method compared to non-regularized NNLS, the solution is still sensitive to the underlying noise and the regularization weight. This is especially relevant for clinically achievable signal-to-noise ratios. In the literature of inverse problems, various well-established approaches to promote smooth solutions, including first-order and second-order Tikhonov regularization, and different criteria for estimating the regularization weight have been proposed, such as L-curve, Generalized Cross-Validation, and Chi-square residual fitting. However, quantitative comparisons between the available reconstruction methods for computing the T₂ distribution, and between different approaches for selecting the optimal regularization weight, are lacking. In this study, we implemented and evaluated ten reconstruction algorithms, resulting from the individual combinations of three penalty terms with three criteria to estimate the regularization weight, plus non-regularized NNLS. Their performance was evaluated both in simulated data and real brain MRI data acquired from healthy volunteers through a scan-rescan repeatability analysis. Our findings demonstrate the need for regularization. As a result of this work, we provide a list of recommendations for selecting the optimal reconstruction algorithms based on the acquired data. Moreover, the implemented methods were packaged in a freely distributed toolbox to promote reproducible research, and to facilitate further research and the use of this promising quantitative technique in clinical practice.

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.

Myelin water fraction estimation using small-tip fast recovery MRI

PURPOSE

To demonstrate the feasibility of an optimized set of small-tip fast recovery (STFR) MRI scans for rapidly estimating myelin water fraction (MWF) in the brain.

METHODS

We optimized a set of STFR scans to minimize the Cramér-Rao Lower Bound of MWF estimates. We evaluated the RMSE of MWF estimates from the optimized scans in simulation. We compared STFR-based MWF estimates (both modeling exchange and not modeling exchange) to multi-echo spin echo (MESE)-based estimates. We used the optimized scans to acquire in vivo data from which a MWF map was estimated. We computed the STFR-based MWF estimates using PERK, a recently developed kernel regression technique, and the MESE-based MWF estimates using both regularized non-negative least squares (NNLS) and PERK.

RESULTS

In simulation, the optimized STFR scans led to estimates of MWF with low RMSE across a range of tissue parameters and across white matter and gray matter. The STFR-based MWF estimates that modeled exchange compared well to MESE-based MWF estimates in simulation. When the optimized scans were tested in vivo, the MWF map that was estimated using a 3-compartment model with exchange was closer to the MESE-based MWF map.

CONCLUSIONS

The optimized STFR scans appear to be well suited for estimating MWF in simulation and in vivo when we model exchange in training. In this case, the STFR-based MWF estimates are close to the MESE-based estimates.

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.

Fast and high-resolution myelin water imaging: Accelerating multi-echo GRASE with CAIPIRINHA

PURPOSE

Although several MRI methods have been explored to achieve in vivo myelin quantification, imaging the whole brain in clinically acceptable times and sufficiently high resolution remains challenging. To address this problem, this work investigates the acceleration of multi-echo T₂ acquisitions based on the multi-echo gradient and spin echo (GRASE) sequence using CAIPIRINHA undersampling and adapted k-space reordering patterns.

METHODS

A prototype multi-echo GRASE sequence supporting CAIPIRINHA parallel imaging was implemented. Multi-echo T₂ data were acquired from 12 volunteers using the implemented sequence (1.6 × 1.6 × 1.6 mm³ , 84 slices, acquisition time [TA] = 10:30 min) and a multi-echo spin echo (MESE) sequence as reference (1.6 × 1.6 × 3.2 mm³ , single-slice, TA = 5:41 min). Myelin water fraction (MWF) maps derived from both acquisitions were compared via correlation and Bland-Altman analyses. In addition, scan-rescan datasets were acquired to evaluate the repeatability of the derived maps.

RESULTS

Resulting maps from the MESE and multi-echo GRASE sequences were found to be correlated (r = 0.83). The Bland-Altman analysis revealed a mean bias of -0.2% (P = .24) with the limits of agreement ranging from -3.7% to 3.3%. The Pearson's correlation coefficient among MWF values obtained from the scan-rescan datasets was found to be 0.95 and the mean bias equal to 0.11% (P = .32), indicating good repeatability of the retrieved maps.

CONCLUSION

By combining a 3D multi-echo GRASE sequence with CAIPIRINHA sampling, whole-brain MWF maps were obtained in 10:30 min with 1.6 mm isotropic resolution. The good correlation with conventional MESE-based maps demonstrates that the implemented sequence may be a promising alternative to time-consuming MESE acquisitions.

Sensitivity to myelin using model-free analysis of the water resonance line-shape in postmortem mouse brain

PURPOSE

Dysmyelinating diseases are characterized by abnormal myelin formation and function. Such microstructural abnormalities in myelin have been demonstrated to produce measurable effects on the MR signal. This work examines these effects on measurements of voxel-wise, high-resolution water spectra acquired using a 3D echo-planar spectroscopic imaging (EPSI) pulse sequence from both postmortem fixed control mouse brains and a dysmyelination mouse brain model.

METHODS

Perfusion fixed, resected control (n = 5) and shiverer (n = 4) mouse brains were imaged using 3D-EPSI with 100 µm isotropic resolution. The free induction decay (FID) was sampled every 2.74 ms over 192 echoes, for a total sampling duration of 526.08 ms. Voxel-wise FIDs were Fourier transformed to produce water spectra with 1.9 Hz resolution. Spectral asymmetry was computed and compared between the two tissue types.

RESULTS

The water resonance is more asymmetrically broadened in the white matter of control mouse brain compared with dysmyelinated white matter. In control brain, this is modulated by and consistent with previously reported orientationally dependent effects of white matter relative to B₀ . Similar sensitivity to orientation is observed in dysmyelinated white matter as well; however, the magnitude of the resonance asymmetry is much lower across all directions.

CONCLUSION

Results demonstrate that components of the spectra are specifically differentially affected by myelin concentration. This suggests that water proton spectra may be sensitive to the presence of myelin, and as such, could serve as a MRI-based biomarker of dysmyelinating disease, free of mathematical models.

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.

Field strength influences on gradient recalled echo MRI signal compartment frequency shifts

Different trends of echo time dependent gradient recalled echo MRI signals in different brain regions have been attributed to signal compartments in image voxels. It remains unclear how variations in gradient recalled echo MRI signals change as a function of MRI field strength, and how data processing may impact signal compartment parameters. We used two popular quantitative susceptibility mapping methods of processing raw phase images (Laplacian and path-based unwrapping with V-SHARP) and expressed values in the form of induced frequency shifts (in Hz) in six specific brain regions at 3T and 7T. We found the frequency shift curves to vary with echo time, and a good overlap between 3T and 7T mean frequency shift curves was present. However, the amount of variation across participants was greater at 3T, and we were able to obtain better compartment model fits of the signal at 7T. We also found the temporal trends in the signal and compartment frequency shifts to change with the method used to process images. The inter-participant averaged trends were consistent between 3T and 7T for each quantitative susceptibility pipeline. However, signal compartment frequency shifts generated using different pipelines may not be comparable.

Advances in MRI of the myelin bilayer

Myelin plays a key role in the function of the central nervous system and is involved in many neurodegenerative diseases. Hence, depiction of myelin is desired for both research and diagnosis. However, MRI of the lipid bilayer constituting the myelin membrane is hampered by extremely rapid signal decay and cannot be accomplished with conventional sequences. Dedicated short-T₂ techniques have therefore been employed, yet with extended sequence timings not well matched to the rapid transverse relaxation in the bilayer, which leads to signal loss and blurring. In the present work, capture and encoding of the ultra-short-T₂ signals in the myelin bilayer is considerably improved by employing advanced short-T₂ methodology and hardware, in particular a high-performance human-sized gradient insert. The approach is applied to tissue samples excised from porcine brain and in vivo in a human volunteer. It is found that the rapidly decaying non-aqueous components in the brain can indeed be depicted with MRI at useful resolution. As a considerable fraction of these signals is related to the myelin bilayer, the presented approach has strong potential to contribute to myelin research and diagnosis.

Influence of 7T GRE-MRI Signal Compartment Model Choice on Tissue Parameters

Quantitative assessment of tissue microstructure is important in studying human brain diseases and disorders. Ultra-high field magnetic resonance imaging (MRI) data obtained using a multi-echo gradient echo sequence have been shown to contain information on myelin, axonal, and extracellular compartments in tissue. Quantitative assessment of water fraction, relaxation time (T₂*), and frequency shift using multi-compartment models has been shown to be useful in studying white matter properties via specific tissue parameters. It remains unclear how tissue parameters vary with model selection based on 7T multiple echo time gradient-recalled echo (GRE) MRI data. We applied existing signal compartment models to the corpus callosum and investigated whether a three-compartment model can be reduced to two compartments and still resolve white matter parameters [i.e., myelin water fraction (MWF) and g-ratio]. We show that MWF should be computed using a three-compartment model in the corpus callosum, and the g-ratios obtained using three compartment models are consistent with previous reports. We provide results for other parameters, such as signal compartment frequency shifts.

Strong diffusion gradients allow the separation of intra- and extra-axonal gradient-echo signals in the human brain

The quantification of brain white matter properties is a key area of application of Magnetic Resonance Imaging (MRI), with much effort focused on using MR techniques to quantify tissue microstructure. While diffusion MRI probes white matter (WM) microstructure by characterising the sensitivity of Brownian motion of water molecules to anisotropic structures, susceptibility-based techniques probe the tissue microstructure by observing the effect of interaction between the tissue and the magnetic field. Here, we unify these two complementary approaches by combining ultra-strong (300 mT/m) gradients with a novel Diffusion-Filtered Asymmetric Spin Echo (D-FASE) technique. Using D-FASE we can separately assess the evolution of the intra- and extra-axonal signals under the action of susceptibility effects, revealing differences in the behaviour in different fibre tracts. We observed that the effective relaxation rate of the ASE signal in the corpus callosum decreases with increasing b-value in all subjects (from 17.1 ± 0.7 s⁻¹ at b = 0 s/mm² to 14.6 ± 0.7 s⁻¹ at b = 4800 s/mm²), while this dependence on b in the corticospinal tract is less pronounced (from 12.0± 1.1 s⁻¹ at b = 0s/mm² to 10.7 ± 0.5 s⁻¹ at b = 4800 s/mm²). Voxelwise analysis of the signal evolution with respect to b-factor and acquisition delay using a microscopic model demonstrated differences in gradient echo signal evolution between the intra- and extra-axonal pools.

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.

Multi-Exponential Relaxometry Using l1 -Regularized Iterative NNLS (MERLIN) With Application to Myelin Water Fraction Imaging

A new parameter estimation algorithm, MERLIN, is presented for accurate and robust multi-exponential relaxometry using magnetic resonance imaging, a tool that can provide valuable insight into the tissue microstructure of the brain. Multi-exponential relaxometry is used to analyze the myelin water fraction and can help to detect related diseases. However, the underlying problem is ill-conditioned, and as such, is extremely sensitive to noise and measurement imperfections, which can lead to less precise and more biased parameter estimates. MERLIN is a fully automated, multi-voxel approach that incorporates state-of-the-art l₁ -regularization to enforce sparsity and spatial consistency of the estimated distributions. The proposed method is validated in simulations and in vivo experiments, using a multi-echo gradient-echo (MEGE) sequence at 3 T. MERLIN is compared to the conventional single-voxel l₂ -regularized NNLS (rNNLS) and a multi-voxel extension with spatial priors (rNNLS + SP), where it consistently showed lower root mean squared errors of up to 70 percent for all parameters of interest in these simulations.