Nonexponential T₂ decay in white matter

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.

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

PURPOSE

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

THEORY AND METHODS

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

RESULTS

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

CONCLUSION

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

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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

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

PURPOSE

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

THEORY AND METHODS

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

RESULTS

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

CONCLUSION

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

Rapid Whole-Knee Quantification of Cartilage Using T1, T2*, and TRAFF2 Mapping With Magnetic Resonance Fingerprinting

OBJECTIVE

Quantitative Magnetic Resonance Imaging (MRI) holds great promise for the early detection of cartilage deterioration. Here, a Magnetic Resonance Fingerprinting (MRF) framework is proposed for comprehensive and rapid quantification of T1, T2*, and TRAFF2 with whole-knee coverage.

METHODS

A MRF framework was developed to achieve quantification of Relaxation Along a Fictitious Field in the 2nd rotating frame of reference ( TRAFF2) along with T1 and T2*. The proposed sequence acquires 65 measurements of 25 high-resolution slices, interleaved with 7 inversion pulses and 40 RAFF2 trains, for whole-knee quantification in a total acquisition time of 3:25 min. Comparison with reference T1, T2*, and TRAFF2 methods was performed in phantom and in seven healthy subjects at 3 T. Repeatability (test-retest) with and without repositioning was also assessed.

RESULTS

Phantom measurements resulted in good agreement between MRF and the reference with mean biases of -54, 2, and 5 ms for T1, T2*, and TRAFF2, respectively. Complete characterization of the whole-knee cartilage was achieved for all subjects, and, for the femoral and tibial compartments, a good agreement between MRF and reference measurements was obtained. Across all subjects, the proposed MRF method yielded acceptable repeatability without repositioning ( R2 ≥ 0.94) and with repositioning ( R2 ≥ 0.57) for T1, T2*, and TRAFF2.

SIGNIFICANCE

The short scan time combined with the whole-knee coverage makes the proposed MRF framework a promising candidate for the early assessment of cartilage degeneration with quantitative MRI, but further research may be warranted to improve repeatability after repositioning and assess clinical value in patients.

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.

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.

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.

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.

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.

Adaptive slice-specific z-shimming for 2D spoiled gradient-echo sequences

Purpose

To reduce the misbalance between compensation gradients and macroscopic field gradients, we introduce an adaptive slice‐specific z‐shimming approach for 2D spoiled multi‐echo gradient‐echoe sequences in combination with modeling of the signal decay.

Methods

Macroscopic field gradients were estimated for each slice from a fast prescan (15 seconds) and then used to calculate slice‐specific compensation moments along the echo train. The coverage of the compensated field gradients was increased by applying three positive and three negative moments. With a forward model, which considered the effect of the slice profile, the z‐shim moment, and the field gradient, R2∗ maps were estimated. The method was evaluated in phantom and in vivo measurements at 3 T and compared with a spoiled multi‐echo gradient‐echo and a global z‐shimming approach without slice‐specific compensation.

Results

The proposed method yielded higher SNR in R2∗ maps due to a broader range of compensated macroscopic field gradients compared with global z‐shimming. In global white matter, the mean interquartile range, proxy for SNR, could be decreased to 3.06 s⁻¹ with the proposed approach, compared with 3.37 s⁻¹ for global z‐shimming and 3.52 s⁻¹ for uncompensated multi‐echo gradient‐echo.

Conclusion

Adaptive slice‐specific compensation gradients between echoes substantially improved the SNR of R2∗ maps, and the signal could also be rephased in anatomical areas, where it has already been completely dephased.

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.

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.

Longitudinal evaluation of demyelinated lesions in a multiple sclerosis model using ultrashort echo time magnetization transfer (UTE-MT) imaging

Alterations in myelin integrity are involved in many neurological disorders and demyelinating diseases, such as multiple sclerosis (MS). Although magnetic resonance imaging (MRI) is the gold standard method to diagnose and monitor MS patients, clinically available MRI protocols show limited specificity for myelin detection, notably in cerebral grey matter areas. Ultrashort echo time (UTE) MRI has shown great promise for direct imaging of lipids and myelin sheaths, and thus holds potential to improve lesion detection. In this study, we used a sequence combining magnetization transfer (MT) with UTE ("UTE-MT", TE ​= ​76 ​μs) and with short TE ("STE-MT", TE ​= ​3000 ​μs) to evaluate spatial and temporal changes in brain myelin content in the cuprizone mouse model for MS on a clinical 7 ​T scanner. During demyelination, UTE-MT ratio (UTE-MTR) and STE-MT ratio (STE-MTR) values were significantly decreased in most white matter and grey matter regions. However, only UTE-MTR detected cortical changes. After remyelination in subcortical and cortical areas, UTE-MTR values remained lower than baseline values, indicating that UTE-MT, but not STE-MT, imaging detected long-lasting changes following a demyelinating event. Next, we evaluated the potential correlations between imaging values and underlying histopathological markers. The strongest correlation was observed between UTE-MTR and percent coverage of myelin basic protein (MBP) immunostaining (r² ​= ​0.71). A significant, although lower, correlation was observed between STE-MTR and MBP (r² ​= ​0.48), and no correlation was found between UTE-MTR or STE-MTR and gliosis immunostaining. Interestingly, correlations varied across brain substructures. Altogether, our results demonstrate that UTE-MTR values significantly correlate with myelin content as measured by histopathology, not only in white matter, but also in subcortical and cortical grey matter regions in the cuprizone mouse model for MS. Readily implemented on a clinical 7 ​T system, this approach thus holds great potential for detecting demyelinating/remyelinating events in both white and grey matter areas in humans. When applied to patients with neurological disorders, including MS patient populations, UTE-MT methods may improve the non-invasive longitudinal monitoring of brain lesions, not only during disease progression but also in response to next generation remyelinating therapies.

Frequency difference mapping applied to the corpus callosum at 7T

Purpose

Frequency difference mapping (FDM) is a phase processing technique which characterizes the nonlinear temporal evolution of the phase of gradient echo (GE) signals. Here, a novel FDM‐processing algorithm is introduced, which is shown to reveal information about white matter microstructure. Unlike some other phase‐processing techniques, the FDM algorithm presented here does not require the use of phase unwrapping or sophisticated image processing. It uses a series of scaled complex divisions to unwrap phase and remove background fields.

Methods

Ten healthy subjects underwent a series of single‐slice, sagittal multi‐echo GE scans at 7T with the slice positioned at the midline. Phase data were processed with the novel FDM algorithm, and the temporal evolution of the magnitude signal and frequency difference was examined in 5 regions of the corpus callosum (CC; genu, anterior body, middle body, posterior body, and splenium).

Results

Consistent frequency difference contrast relative to surrounding tissue was observed in all subjects in the CC and in other white matter regions where the nerve fibers run perpendicular to B0, such as the superior cerebellar peduncle. Examination of the frequency difference curves shows distinct variations over the CC, with the genu and splenium displaying larger frequency differences than the other regions (in addition to a faster decay of signal magnitude).

Conclusion

The novel FDM algorithm presented here yields images sensitive to tissue microstructure and microstructural differences over the CC in a simple manner, without the requirement for phase unwrapping or sophisticated image processing.

Diffusion effects in myelin sheath free induction decay

Myelin sheath microstructure and composition produce MR signal decay characteristics that can be used to evaluate status and outcome of demyelinating disease. We extend a recently proposed model of neuronal magnetic susceptibility, that accounts for both the structural and inherent anisotropy of the myelin sheath, by including the whole dynamic range of diffusion effects. The respective Bloch-Torrey equation for local spin dephasing is solved with a uniformly convergent perturbation expansion method, and the resulting magnetization decay is validated with a numerical solution based on a finite difference method. We show that a variation of diffusion strengths can lead to substantially different MR signal decay curves. Our results may be used to adjust or control simulations for water diffusion in neuronal structures.

Magnetic resonance fingerprinting with quadratic RF phase for measurement of T2 * simultaneously with δf , T1 , and T2

PURPOSE

This study explores the possibility of using a gradient moment balanced sequence with a quadratically varied RF excitation phase in the magnetic resonance fingerprinting (MRF) framework to quantify T₂ ^* in addition to δ f , T₁ , and T₂ tissue properties.

METHODS

The proposed quadratic RF phase-based MRF method (qRF-MRF) combined a varied RF excitation phase with the existing balanced SSFP (bSSFP)-based MRF method to generate signals that were uniquely sensitive to δ f , T₁ , T₂ , as well as the distribution width of intravoxel frequency dispersion, Γ . A dictionary, generated through Bloch simulation, containing possible signal evolutions within the physiological range of δ f , T₁ , T₂ , and Γ , was used to perform parameter estimation. The estimated T₂ and Γ were subsequently used to estimate T₂ ^* . The proposed method was evaluated in phantom experiments and healthy volunteers (N = 5).

RESULTS

The T₁ and T₂ values from the phantom by qRF-MRF demonstrated good agreement with values obtained by traditional gold standard methods (r² = 0.995 and 0.997, respectively; concordance correlation coefficient = 0.978 and 0.995, respectively). The T₂ ^* values from the phantom demonstrated good agreement with values obtained through the multi-echo gradient-echo method (r² = 0.972, concordance correlation coefficient = 0.983). In vivo qRF-MRF-measured T₁ , T₂ , and T₂ ^* values were compared with measurements by existing methods and literature values.

CONCLUSION

The proposed qRF-MRF method demonstrated the potential for simultaneous quantification of δ f , T₁ , T₂ , and T₂ ^* tissue properties.

The effects of variations in tissue microstructure from postmortem rat brain on the asymmetry of the water proton resonance

PURPOSE

This work was performed to investigate the effects of tissue microstructure from postmortem rat brain on the shape of the water proton spectrum.

METHODS

Perfusion-fixed, resected rat brains (N = 4) were imaged at 9.4T. 3D DTI and 3D echo-planar spectroscopic imaging (EPSI) data were acquired with 150 μm isotropic resolution. DTI data were acquired over 60 directions with b = 3000 s/mm² . Water spectra were produced from EPSI data acquired over 128 echoes, with 2.9 Hz spectral resolution. A voxel-wise metric reflecting spectral asymmetry about the peak of the resonance was computed and compared with orientation estimates from DTI data by fitting data with the susceptibility anisotropy model.

RESULTS

Asymmetric broadening of the water resonance was computed for mixed populations of grey and/or white matter as determined by thresholding the fractional anisotropy. Asymmetry was shown to be differentially affected by tract orientation relative to B₀ in high FA voxels, whereas low FA voxels exhibited little sensitivity. Anatomic structures in the hippocampus were also found to produce distinct changes in the water resonance.

CONCLUSION

Present results demonstrate that structural variations in tissue architecture cause characteristic, reproducible changes in the water resonance shape. This suggests that water spectra are sensitive to cytoarchitectural variations in brain tissue.

Signal compartments in ultra-high field multi-echo gradient echo MRI reflect underlying tissue microstructure in the brain

Gradient recalled echo magnetic resonance imaging (GRE-MRI) at ultra-high field holds great promise for new contrast mechanisms and delineation of putative tissue compartments that contribute to the multi-echo GRE-MRI signal may aid structural characterization. Several studies have adopted the three water-pool compartment model to study white matter brain regions, associating individual compartments with myelin, axonal and extracellular water. However, the number and identifiability of GRE-MRI signal compartments has not been fully explored. We undertook this task for human brain imaging data. Multiple echo time GRE-MRI data were acquired in five healthy participants, specific anatomical structures were segmented in each dataset (substantia nigra, caudate, insula, putamen, thalamus, fornix, internal capsule, corpus callosum and cerebrospinal fluid), and the signal fitted with models comprising one to six signal compartments using a complex-valued plane wave formulation. Information criteria and cluster analysis methods were used to ascertain the number of distinct compartments within the signal from each structure and to determine their respective frequency shifts. We identified five principal signal compartments with different relative contributions to each structure's signal. Voxel-based maps of the volume fraction of each of these compartments were generated and demonstrated spatial correlation with brain anatomy.

MRI gradient-echo phase contrast of the brain at ultra-short TE with off-resonance saturation

Larmor-frequency shift or image phase measured by gradient-echo sequences has provided a new source of MRI contrast. This contrast is being used to study both the structure and function of the brain. So far, phase images of the brain have been largely obtained at long echo times as maximum phase signal-to-noise ratio (SNR) is achieved at TE = T2* (∼40 ms at 3T). The structures of the brain, however, are compartmentalized and complex with a wide range of signal relaxation times. At such long TE, the short-T2 components are largely attenuated and contribute minimally to phase contrast. The purpose of this study was to determine whether proton gradient-echo images of the brain exhibit phase contrast at ultra-short TE (UTE). Our data showed that UTE images acquired at 7 T without off-resonance saturation do not contain significant phase contrast between gray and white matter. However, UTE images of the brain can attain strong phase contrast even at a nominal TE of 106 μs by using off-resonance RF saturation pulses, which provide direct saturation of ultra-short-T2 components and indirect saturation of longer-T2 components via magnetization transfer. In addition, phase contrast between gray and white matter acquired at UTE with off-resonance saturation is reversed compared to that of the long-T2 signals acquired at long TEs. This finding opens up a potential new way to manipulate image phase contrast of the brain. By accessing short and ultra-short-T2 species, MRI phase images may further improve the characterization of tissue microstructure in the brain.

Studying brain microstructure with magnetic susceptibility contrast at high-field

A rapidly developing application of high field MRI is the study of brain anatomy and function with contrast based on the magnetic susceptibility of tissues. To study the subtle variations in susceptibility contrast between tissues and with changes in brain activity, dedicated scan techniques such as susceptibility-weighted MRI and blood-oxygen level dependent functional MRI have been developed. Particularly strong susceptibility contrast has been observed with systems that operate at 7T and above, and their recent widespread use has led to an improved understanding of contributing sources and mechanisms. To interpret magnetic susceptibility contrast, analysis approaches have been developed with the goal of extracting measures that report on local tissue magnetic susceptibility, a physical quantity that, under certain conditions, allows estimation of blood oxygenation, local tissue iron content, and quantification of its changes with disease. Interestingly, high field studies have also brought to light that not only the makeup of tissues affects MRI susceptibility contrast, but that also a tissue's sub-voxel structure at scales all the way down to the molecular level plays an important role as well. In this review, various ways will be discussed by which sub-voxel structure can affect the MRI signal in general, and magnetic susceptibility in particular, sometimes in a complex fashion. In the light of this complexity, it appears likely that accurate, brain-wide quantification of iron will require the combination of multiple contrasts that may include diffusion and magnetization transfer information with susceptibility-weighted contrast. On the other hand, this complexity also offers opportunities to use magnetic susceptibility contrast to inform about specific microstructural aspects of brain tissue. Details and several examples will be presented in this review.

Inferring brain tissue composition and microstructure via MR relaxometry

MRI relaxometry is sensitive to a variety of tissue characteristics in a complex manner, which makes it both attractive and challenging for characterizing tissue. This article reviews the most common water proton relaxometry measures, T₁, T₂, and T₂^*, and reports on their development and current potential to probe the composition and microstructure of brain tissue. The development of these relaxometry measures is challenged by the need for suitably accurate tissue models, as well as robust acquisition and analysis methodologies. MRI relaxometry has been established as a tool for characterizing neural tissue, particular with respect to myelination, and the potential for further development exists.

Investigating the spatiotemporal characteristics of the deoxyhemoglobin-related and deoxyhemoglobin-unrelated functional hemodynamic response across cortical layers in awake marmosets

Blood oxygenation level dependent (BOLD) functional magnetic resonance imaging (fMRI) has become a major tool to map neural activity. However, the spatiotemporal characteristics of the BOLD functional hemodynamic response across the cortical layers remain poorly understood. While human fMRI studies suffer from low spatiotemporal resolution, the use of anesthesia in animal models introduces confounding factors. Additionally, inflow contributions to the fMRI signal become non-negligible when short repetition times (TRs) are used. In the present work, we mapped the BOLD fMRI response to somatosensory stimulation in awake marmosets. To address the above technical concerns, we used a dual-echo gradient-recalled echo planar imaging (GR-EPI) sequence to separate the deoxyhemoglobin-related response (absolute T2* differences) from the deoxyhemoglobin-unrelated response (relative S0 changes). We employed a spatial saturation pulse to saturate incoming arterial spins and reduce inflow effects. Functional GR-EPI images were obtained from a single coronal slice with two different echo times (13.5 and 40.5ms) and TR=0.2s. BOLD, T2*, and S0 images were calculated and their functional responses were detected in both hemispheres of primary somatosensory cortex, from which five laminar regions (L1+2, L3, L4, L5, and L6) were derived. The spatiotemporal distribution of the BOLD response across the cortical layers was heterogeneous, with the middle layers having the highest BOLD amplitudes and shortest onset times. ΔT2* also showed a similar trend. However, functional S0 changes were detected only in L1+2, with a fast onset time. Because inflow effects were minimized, the source of S0 functional changes in L1+2 could be attributed to a reduction of cerebrospinal fluid volume fraction due to the functional increase in cerebral blood volume and to unmodeled T2* changes in the extra- and intra-venous compartments. Caution should be exercised when interpreting laminar BOLD fMRI changes in superficial layers as surrogates of underlying neural activity.

In vivo characterization of brain ultrashort-T2 components

PURPOSE

Recent nuclear magnetic resonance and MRI studies have measured a fast-relaxing signal component with T2∗<1 ms in white matter and myelin extracts. In ex vivo studies, evidence suggests that a large fraction of this component directly arises from bound protons in the myelin phospholipid membranes. Based on these results, this ultrashort-T₂ component in nervous tissue is a new potential imaging biomarker of myelination, which plays a critical role in neuronal signal conduction across the brain and loss or degradation of myelin is a key feature of many neurological disorders. The goal of this work was to characterize the relaxation times and frequency shifts of ultrashort-T₂ components in the human brain.

METHODS

This required development of an ultrashort echo time relaxometry acquisition strategy and fitting procedure for robust measurements in the presence of ultrashort T2∗ relaxation times and large frequency shifts.

RESULTS

We measured an ultrashort-T₂ component in healthy volunteers with a median T2∗ between 0.5-0.7 ms at 3T and 0.2-0.3 ms at 7T as well as an approximately -3 ppm frequency shift from water.

CONCLUSION

To our knowledge, this is the first time a chemical shift of the ultrashort-T₂ brain component has been measured in vivo. This chemical shift, at around 1.7 ppm, is similar to the primary resonance of most lipids, indicating that much of the ultrashort-T₂ component observed in vivo arises from bound protons in the myelin phospholipid membranes. Magn Reson Med 80:726-735, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Assessment of microstructural signal compartments across the corpus callosum using multi-echo gradient recalled echo at 7 T

Quantitative assessment of tissue microstructure is important in studying human brain diseases and disorders in which white matter is implicated, as it has been linked to demyelination, re-myelination, and axonal damage in clinical conditions. Ultra-high field magnetic resonance imaging data obtained using a multi-echo gradient echo sequence has been shown to contain information on myelin, axonal and extracellular compartments in white matter. In this study, we aimed to assess the sensitivity of a three-compartment model to estimate the variation of corresponding compartment parameters (water fraction, relaxation time and frequency shift) of the corpus callosum sub-regions, which are known to have different tissue structure. Additionally, we computed the g-ratio using myelin and axonal water fractions and performed a voxel-by-voxel analysis in the corpus callosum. Based on data acquired for ten participants, we show that the myelin compartment water fraction and T₂^∗ is consistent across the corpus callosum sub-regions, whilst myelin frequency shift varies. The results show that the variation in water fraction, T₂^∗ and frequency shift for the myelin signal compartment across the corpus callosum is smaller than for the axonal and extracellular signal compartments. The computed g-ratio was comparable to previously published studies in the corpus callosum. Our study suggests that a multi-echo GRE approach in vivo combined with a complex three-compartment model is sensitive to microstructural parameter variations across the human corpus callosum.

Magnetic resonance microdynamic imaging reveals distinct tissue microenvironments

Magnetic resonance imaging (MRI) provides a powerful set of tools with which to investigate biological tissues noninvasively and in vivo. Tissues are heterogeneous in nature; an imaging voxel contains an ensemble of different cells and extracellular matrix components. A long-standing challenge has been to infer the content of and interactions among these microscopic tissue components within a macroscopic imaging voxel. Spatially resolved multidimensional relaxation-diffusion correlation (REDCO) spectroscopy holds the potential to deliver such microdynamic information. However, to date, vast data requirements have mostly relegated these type of measurements to nuclear magnetic resonance applications and prevented them from being widely and successfully used in conjunction with imaging. By using a novel data acquisition and processing strategy in this study, spatially resolved REDCO could be performed in reasonable scanning times with excellent prospects for clinical applications. This new MR imaging framework-which we term "magnetic resonance microdynamic imaging (MRMI)"-permits the simultaneous noninvasive and model-free quantification of multiple subcellular, cellular, and interstitial tissue microenvironments within a voxel. MRMI is demonstrated with a fixed spinal cord specimen, enabling the quantification of microscopic tissue components with unprecedented specificity. Tissue components, such as axons, neuronal and glial soma, and myelin were identified on the basis of their multispectral signature within individual imaging voxels. These tissue elements could then be composed into images and be correlated with immunohistochemistry findings. MRMI provides novel image contrasts of tissue components and a new family of microdynamic biomarkers that could lead to new diagnostic imaging approaches to probe biological tissue alterations accompanied by pathological or developmental changes.

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.

PLANET: An ellipse fitting approach for simultaneous T1 and T2 mapping using phase-cycled balanced steady-state free precession

Purpose

To demonstrate the feasibility of a novel, ellipse fitting approach, named PLANET, for simultaneous estimation of relaxation times T₁ and T₂ from a single 3D phase‐cycled balanced steady‐state free precession (bSSFP) sequence.

Methods

A method is presented in which the elliptical signal model is used to describe the phase‐cycled bSSFP steady‐state signal. The fitting of the model to the acquired data is reformulated into a linear convex problem, which is solved directly by a linear least squares method, specific to ellipses. Subsequently, the relaxation times T₁ and T₂, the banding free magnitude, and the off‐resonance are calculated from the fitting results.

Results

Maps of T₁ and T₂, as well as an off‐resonance and a banding free magnitude can be simultaneously, quickly, and robustly estimated from a single 3D phase‐cycled bSSFP sequence. The feasibility of the method was demonstrated in a phantom and in the brain of healthy volunteers on a clinical MR scanner. The results were in good agreement for the phantom, but a systematic underestimation of T₁ was observed in the brain.

Conclusion

The presented method allows for accurate mapping of relaxation times and off‐resonance, and for the reconstruction of banding free magnitude images at realistic signal‐to‐noise ratios. Magn Reson Med 79:711–722, 2018. © 2017 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution‐NonCommercial‐NoDerivs License, which permits use and distribution in any medium, provided the original work is properly cited, the use is non‐commercial and no modifications or adaptations are made.

The larmor frequency shift in magnetically heterogeneous media depends on their mesoscopic structure

PURPOSE

Recent studies have addressed the determination of the NMR precession frequency in biological tissues containing magnetic susceptibility differences between cell types. The purpose of this study is to investigate the dependence of the precession frequency on medium microstructure using a simple physical model.

THEORY

This dependence is governed by diffusion of NMR-visible molecules in magnetically heterogeneous microenvironments. In the limit of fast diffusion, the precession frequency is determined by the average susceptibility-induced magnetic field, whereas in the limit of slow diffusion it is determined by the average local phase factor of precessing spins.

METHODS

The main method used is Monte Carlo simulation of isotropic suspensions of impermeable magnetized spheres. In addition, NMR spectroscopy was performed in aqueous suspensions of polystyrene microbeads.

RESULTS

The precession frequency depends on the structural organization of magnetized objects in the medium. Monte Carlo simulations demonstrated a nonmonotonic transition between the regimes of fast and slow diffusion. NMR experiments confirmed the transition, but were unable to confirm its precise form. Results for a given pattern of structural organization obey a scaling law.

CONCLUSION

The NMR precession frequency exhibits a complex dependence on medium structure. Our results suggest that the commonly assumed limit of fast water diffusion holds for biological tissues with small cells. Magn Reson Med 79:1101-1110, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Contributions to magnetic susceptibility of brain tissue

This review discusses the major contributors to the subtle magnetic properties of brain tissue and how they affect MRI contrast. With the increased availability of high-field scanners, the use of magnetic susceptibility contrast for the study of human brain anatomy and function has increased dramatically. This has not only led to novel applications, but has also improved our understanding of the complex relationship between MRI contrast and magnetic susceptibility. Chief contributors to the magnetic susceptibility of brain tissue have been found to include myelin as well as iron. In the brain, iron exists in various forms with diverse biological roles, many of which are now only starting to be uncovered. An interesting aspect of magnetic susceptibility contrast is its sensitivity to the microscopic distribution of iron and myelin, which provides opportunities to extract information at spatial scales well below MRI resolution. For example, in white matter, the myelin sheath that surrounds the axons can provide tissue contrast that is dependent on the axonal orientation and reflects the relative size of intra- and extra-axonal water compartments. The extraction of such ultrastructural information, together with quantitative information about iron and myelin concentrations, is an active area of research geared towards the characterization of brain structure and function, and their alteration in disease. Copyright © 2016 John Wiley & Sons, Ltd.

Contributions to magnetic susceptibility of brain tissue

This review discusses the major contributors to the subtle magnetic properties of brain tissue and how they affect MRI contrast. With the increased availability of high-field scanners, the use of magnetic susceptibility contrast for the study of human brain anatomy and function has increased dramatically. This has not only led to novel applications, but has also improved our understanding of the complex relationship between MRI contrast and magnetic susceptibility. Chief contributors to the magnetic susceptibility of brain tissue have been found to include myelin as well as iron. In the brain, iron exists in various forms with diverse biological roles, many of which are now only starting to be uncovered. An interesting aspect of magnetic susceptibility contrast is its sensitivity to the microscopic distribution of iron and myelin, which provides opportunities to extract information at spatial scales well below MRI resolution. For example, in white matter, the myelin sheath that surrounds the axons can provide tissue contrast that is dependent on the axonal orientation and reflects the relative size of intra- and extra-axonal water compartments. The extraction of such ultrastructural information, together with quantitative information about iron and myelin concentrations, is an active area of research geared towards the characterization of brain structure and function, and their alteration in disease. Copyright © 2016 John Wiley & Sons, Ltd.