Myelin Water Fraction Mapping Research Articles

Quantifying Brain Myelin Water Fraction in a Guinea Pig Model of Spontaneous Intrauterine Growth Restriction.

Intrauterine growth restriction (IUGR) is an obstetrical condition where a fetus has not achieved its genetic potential. A consequence of IUGR is a decrease in brain myelin content. Myelin water imaging (MWI) has been used to assess fetal myelin water fraction (MWF) and might potentially assess myelination changes associated with IUGR. To quantify and compare the MWF of non-IUGR and IUGR fetal guinea pigs (GPs) in late gestation. Prospective animal model. Dunkin-Hartley GP model of spontaneous IUGR (mean ± SD: 60 ± 1.2 days gestation; 19 IUGR, 52 control). Eight spoiled gradient-recalled (SPGR) gradient echo volumes (flip angles [α]: 2°-16°), and two sets of eight balanced steady-state free precession (bSSFP) gradient echo volumes (α: 8° - 64°), at 0° and 180° phase increments, at 3.0 T. MWF maps were generated for each fetal GP brain using multicomponent driven equilibrium single pulse observation of T1 /T2 (mcDESPOT). MWF was quantified in the fetal corpus callosum (CC), fornix (FOR), parasagittal white matter (PSW), and whole fetal brain. Linear regression was performed between five fetal IUGR markers (body volume, body-to-pregnancy volume ratio, brain-to-liver volume ratio, brain-to-placenta volume ratio, and brain-to-body volume ratio) and MWF (coefficient of determination, R2 ). A t-test with a linear mixed model compared the MWF of non-IUGR and IUGR fetal GPs (significance was determined at α < 0.05). The MWF of the control fetuses are (mean ± SD): 0.23 ± 0.02 (CC), 0.31 ± 0.02 (FOR), 0.28 ± 0.02 (PSW), and 0.20 ± 0.01 (whole brain). The MWF of the IUGR fetuses are (mean ± SD): 0.19 ± 0.02 (CC), 0.27 ± 0.01 (FOR), 0.24 ± 0.03 (PSW), and 0.16 ± 0.01 (whole brain). Significant differences in MWF were found between the non-IUGR and IUGR fetuses in every comparison. The mean MWF of IUGR fetal GPs is significantly lower than non-IUGR fetal GPs. 2 TECHNICAL EFFICACY: Stage 1.

Magnetic resonance fingerprinting-based myelin water fraction mapping for the assessment of white matter maturation and integrity in typical development and leukodystrophies.

A quantitative biomarker for myelination, such as myelin water fraction (MWF), would boost the understanding of normative and pathological neurodevelopment, improving patients' diagnosis and follow-up. We quantified the fraction of a rapidly relaxing pool identified as MW using multicomponent three-dimensional (3D) magnetic resonance fingerprinting (MRF) to evaluate white matter (WM) maturation in typically developing (TD) children and alterations in leukodystrophies (LDs). We acquired DTI and 3D MRF-based R1, R2 and MWF data of 15 TD children and 17 LD patients (9 months-12.5 years old) at 1.5 T. We computed normative maturation curves in corpus callosum and corona radiata and performed WM tract profile analysis, comparing MWF with R1, R2 and fractional anisotropy (FA). Normative maturation curves demonstrated a steep increase for all tissue parameters in the first 3 years of age, followed by slower growth for MWF while R1, R2R2 and FA reached a plateau. Unlike FA, MWF values were similar for regions of interest (ROIs) with different degrees of axonal packing, suggesting independence from fiber bundle macro-organization and higher myelin specificity. Tract profile analysis indicated a specific spatial pattern of myelination in the major fiber bundles, consistent across subjects. LD were better distinguished from TD by MWF rather than FA, showing reduced MWF with respect to age-matched controls in both ROI-based and tract analysis. In conclusion, MRF-based MWF provides myelin-specific WM maturation curves and is sensitive to alteration due to LDs, suggesting its potential as a biomarker for WM disorders. As MRF allows fast simultaneous acquisition of relaxometry and MWF, it can represent a valuable diagnostic tool to study and follow up developmental WM disorders in children.

High-resolution myelin-water fraction and quantitative relaxation mapping using 3D ViSTa-MR fingerprinting.

This study aims to develop a high-resolution whole-brain multi-parametric quantitative MRI approach for simultaneous mapping of myelin-water fraction (MWF), T1, T2, and proton-density (PD), all within a clinically feasible scan time. We developed 3D visualization of short transverse relaxation time component (ViSTa)-MRF, which combined ViSTa technique with MR fingerprinting (MRF), to achieve high-fidelity whole-brain MWF and T1/T2/PD mapping on a clinical 3T scanner. To achieve fast acquisition and memory-efficient reconstruction, the ViSTa-MRF sequence leverages an optimized 3D tiny-golden-angle-shuffling spiral-projection acquisition and joint spatial-temporal subspace reconstruction with optimized preconditioning algorithm. With the proposed ViSTa-MRF approach, high-fidelity direct MWF mapping was achieved without a need for multicompartment fitting that could introduce bias and/or noise from additional assumptions or priors. The in vivo results demonstrate the effectiveness of the proposed acquisition and reconstruction framework to provide fast multi-parametric mapping with high SNR and good quality. The in vivo results of 1 mm- and 0.66 mm-isotropic resolution datasets indicate that the MWF values measured by the proposed method are consistent with standard ViSTa results that are 30× slower with lower SNR. Furthermore, we applied the proposed method to enable 5-min whole-brain 1 mm-iso assessment of MWF and T1/T2/PD mappings for infant brain development and for post-mortem brain samples. In this work, we have developed a 3D ViSTa-MRF technique that enables the acquisition of whole-brain MWF, quantitative T1, T2, and PD maps at 1 and 0.66 mm isotropic resolution in 5 and 15 min, respectively. This advancement allows for quantitative investigations of myelination changes in the brain.

Metabolic Syndrome Score is Associated with Lower Myelin in Cognitively Healthy Adults

AbstractBackgroundMetabolic syndrome (MetS) is a cluster of cardiovascular risk factors, including abdominal obesity, glucose intolerance, elevated blood pressure (BP), and dyslipidemia. MetS is a known risk factor for cognitive dysfunction and dementia (Yaffe et al., 2004. YAMA), and has been linked to brain atrophy and altered WM microstructure measured with diffusion tensor MR (Alfaro et al., 2018, Metabolism). The metabolic cost of myelin maintenance makes WM particularly vulnerable to metabolic dysfunction. Here, we aim to explore the relationships between MetS and myelin, as well as test whether WC‐myelin association (Bouhrara et al., 2021, Int J Obes) is independent from systolic BP.MethodSample included 90 cognitively and neurologically healthy adults (Bookheimer et al., 2019, Neuroimage), aged 20 to 79 (M = 53, SD = 17) with 76% women. Additional exclusion criteria included systolic blood pressure &gt;140, BMI&gt;33, and &gt;2 on modified Hachinski Ischemic Scale. Myelin water imaging was collected using a 3T Siemens Skyra with a 32‐channel coil and 3D‐gradient and spin‐echo (GRASE) sequence (1.9×1.9×4mm, TR = 1000ms, echoes = 32, first echo = 10ms, echo spacing = 10ms). Myelin Water Fraction maps were nonlinearly aligned to FMRIB58_FA using TBSS in FSL. MetS score was calculated as a sum of z‐scored waist circumference (WC), systolic BP, BMI, and self‐reports history of elevated cholesterol and glucose.ResultA linear regression analysis with MWF as the dependent variable, adjusting for sex, race, ethnicity, age, and age2, income and education, indicated that greater MetS score was associated with lower MWF in several regions including the genu corpus callosum, anterior corona radiata, and external capsule. WC was negatively associated with whole WM MWF (B = ‐.265, t(79,9) = ‐2.150, p = 0.035) after controlling for the above covariates and systolic BP.ConclusionAssociations between MetS with myelin content were detected especially in the late‐myelinating regions. This aligns with the development‐to‐degeneration hypothesis of WM aging, with lower oligodendrocyte‐to‐axon ratio as a possible factor contributing to metabolic vulnerability of late‐myelinating tracts. We also replicated the WC‐myelin association, although in fewer regions and with smaller effects sizes. Together, our results suggest that different metabolic syndrome components may exert both additive/synergistic, and independent effects on myelin integrity across adulthood.

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

To investigate the prospects of a multigradient-echo (mGRE) acquisition for in vivo myelin water imaging at 0.55 T. 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. 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. 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.

Improved myelin water fraction mapping with deep neural networks using synthetically generated 3D data

We introduce a generative model for synthesis of large scale 3D datasets for quantitative parameter mapping of myelin water fraction (MWF). Our model combines a MR physics signal decay model with an accurate probabilistic multi-component parametric T2 model. We synthetically generate a wide variety of high quality signals and corresponding parameters from a wide range of naturally occurring prior parameter values. To capture spatial variation, the generative signal decay model is combined with a generative spatial model conditioned on generic tissue segmentations. Synthesized 3D datasets can be used to train any convolutional neural network (CNN) based architecture for MWF estimation. Our source code is available at: https://github.com/quin-med-harvard-edu/synthmapReduction of acquisition time at the expense of lower SNR, as well as accuracy and repeatability of MWF estimation techniques, are key factors that affect the adoption of MWF mapping in clinical practice. We demonstrate that the synthetically trained CNN provides superior accuracy over the competing methods under the constraints of naturally occurring noise levels as well as on the synthetically generated images at low SNR levels. Normalized root mean squared error (nRMSE) is less than 7% on synthetic data, which is significantly lower than competing methods. Additionally, the proposed method yields a coefficient of variation (CoV) that is at least 4x better than the competing method on intra-session test-retest reference dataset.

REUSED: A deep neural network method for rapid whole-brain high-resolution myelin water fraction mapping from extremely under-sampled MRI.

Changes in myelination are a cardinal feature of brain development and the pathophysiology of several central nervous system diseases, including multiple sclerosis and dementias. Advanced magnetic resonance imaging (MRI) methods have been developed to probe myelin content through the measurement of myelin water fraction (MWF). However, the prolonged data acquisition and post-processing times of current MWF mapping methods pose substantial hurdles to their clinical implementation. Recently, fast steady-state MRI sequences have been implemented to produce high-spatial resolution whole-brain MWF mapping within ∼20min. Despite the subsequent significant advances in the inversion algorithm to derive MWF maps from steady-state MRI, the high-dimensional nature of such inversion does not permit further reduction of the acquisition time by data under-sampling. In this work, we present an unprecedented reduction in the computation (∼30s) and the acquisition time (∼7min) required for whole-brain high-resolution MWF mapping through a new Neural Network (NN)-based approach, named NN-Relaxometry of Extremely Under-SamplEd Data (NN-REUSED). Our analyses demonstrate virtually similar accuracy and precision in derived MWF values using NN-REUSED compared to results derived from the fully sampled reference method. The reduction in the acquisition and computation times represents a breakthrough toward clinically practical MWF mapping.

Myelin Water Fraction of Slowly Evolving Lesions is Associated With Disability in Multiple Sclerosis (P10-3.017)

Myelin Water Fraction of Slowly Evolving Lesions is Associated With Disability in Multiple Sclerosis (P10-3.017)

Mapping myelin in white matter with T1-weighted/T2-weighted maps: discrepancy with histology and other myelin MRI measures

The ratio of T1-weighted/T2-weighted magnetic resonance images (T1w/T2w MRI) has been successfully applied at the cortical level since 2011 and is now one of the most used myelin mapping methods. However, no reports have explored the histological validity of T1w/T2w myelin mapping in white matter. Here we compare T1w/T2w with ex vivo postmortem histology and in vivo MRI methods, namely quantitative susceptibility mapping (QSM) and multi-echo T2 myelin water fraction (MWF) mapping techniques. We report a discrepancy between T1w/T2w myelin maps of the human corpus callosum and the histology and analyse the putative causes behind such discrepancy. T1w/T2w does not positively correlate with Luxol Fast Blue (LFB)-Optical Density but shows a weak to moderate, yet significant, negative correlation. On the contrary, MWF is strongly and positively correlated with LFB, whereas T1w/T2w and MWF maps are weakly negatively correlated. The discrepancy between T1w/T2w MRI maps, MWF and histological myelin maps suggests caution in using T1w/T2w as a white matter mapping method at the callosal level. While T1w/T2w imaging may correlate with myelin content at the cortical level, it is not a specific method to map myelin density in white matter.

Feasibility of MRI Quantification of Myelin Water Fraction in the Fetal Guinea Pig Brain.

Fetal myelination assessment is important for understanding neurodevelopment and neurodegeneration. Myelin water imaging (MWI) quantifies myelin water fraction (MWF), a validated marker for myelin content, and has been used to assess brain myelin in children and neonates. To demonstrate that MWI can quantify MWF in fetal guinea pigs (GPs). Animal model. Nine pregnant, Dunkin-Hartley GPs with 31 fetuses (mean ± standard deviation=60 ± 1.5 days gestation). 3D spoiled gradient echo and balanced steady-state free precession sequences at 3.0 T. MWF maps were reconstructed for maternal and fetal GP brains using the multicomponent driven equilibrium single pulse observation of T1 and T2 (mcDESPOT) approach. Myelin basic protein (MBP) stain provided histological validation of the MWF. Regions of interest were placed in the maternal corpus callosum (CC), maternal fornix (FOR), fetal CC, and fetal FOR in MWF maps and MBP stains. Linear regression between MWF and MBP stain intensity (SI) of all four regions (coefficient of determination, R2 ). A paired t-test compared the MWF of maternal and mean fetal CC, MBP SI of maternal and mean fetal CC, MWF of maternal and mean fetal FOR, MBP SI of maternal and mean fetal FOR. A paired t-test with a linear mixed model compared the MWF of fetal CC and fetal FOR, and MBP SI of fetal CC and fetal FOR. A P value < 0.0083 was considered statistically significant. The mean MWF of the analyzed regions are as follows (mean ± standard deviation): 0.338 + 0.016 (maternal CC), 0.340 ± 0.017 (maternal FOR), 0.214 ± 0.016 (fetal CC), and 0.305 ± 0.025 (fetal FOR). MWF correlated with MBP SI in all regions (R2 =0.81). Significant differences were found between MWF and MBP SI of maternal and fetal CC, and MWF and MBP SI of fetal CC and fetal FOR. This study demonstrated the feasibility of MWI in assessing fetal brain myelin content. 2 Technical Efficacy: Stage 1.

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

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. 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. 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. 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 myelin water fraction mapping using complex-valued neural network analysis.

Previously, an artificial neural network method was introduced to estimate quantitative myelin water fraction (MWF) using multi-echo gradient-echo data. However, the fiber orientation of white matter with respect to B0 could bias the quantification of MWF. Here, we developed an advanced workflow for MWF estimation that could improve the quantification of MWF. To adopt fiber orientation effects, a complex-valued neural network with complex-valued operation was used. In addition, to compensate for the bias from different scan parameters, a signal model incorporating the T1 value was devised for training data generation. At the testing stage, a voxel-spread function approach was utilized for spatial B0 artifact correction. Finally, dropout-based variational inference was implemented for uncertainty estimates on the network model to provide a confidence interpretation of the output. According to simulation and in vivo analysis, the proposed method suggests improved quality of MWF estimation by correcting the bias and artifacts. The proposed complex-valued neural network approach can alleviate the dependency of fiber orientation effects compared to previous artificial neural network method. Uncertainty estimates provides information different from fitting error that can be used as a confidence level of the resulting MWF values. An improved MWF mapping using complex-valued neural network analysis has been proposed.

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

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 T2 -prep data and to evaluate its accuracy in comparison with multilayer perceptron (MLP) network. 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. 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. Subsecond whole-brain MWF extraction from fast acquisition with spiral trajectory and T2 -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.

Exploring in vivo lesion myelination dynamics: Longitudinal Myelin Water Imaging in early Multiple Sclerosis

BackgroundMultiple Sclerosis (MS) lesions are pathologically heterogeneous and the temporal behavior in terms of growth and myelination status of individual lesions is highly variable, especially in the early phase of the disease. Thus, monitoring the development of individual lesion myelination by using quantitative magnetic resonance myelin water imaging (MWI) could be valuable to capture the variability of disease pathology and get an individual insight into the subclinical disease activity. ObjectiveThe goal of this work was (1) to observe the variation and longitudinal change of in vivo lesion myelination by means of MWI and its parameter Myelin Water Fraction (MWF), and, (2) to identify individual lesion myelination patterns in early MS. MethodsIn this study n = 12 patients obtained conventional MRI and quantitative MWI derived from multi-component driven equilibrium single pulse observation of T1 and T2 (mcDESPOT) within four weeks after presenting a clinically isolated syndrome and remained within the study if clinically definitive MS was diagnosed within the 12 months study period. Four MRI sessions were acquired at baseline, 3, 6, and 12 months. The short-term and long-term variability of MWF maps was evaluated by scan-rescan measures and the coefficient of variation was determined in four healthy controls. Tracking of individual lesions was performed using the Automatic Follow-up of Individual Lesions (AFIL) algorithm. Lesion volume and MWF were evaluated for every individual lesion in all patients. Median lesion MWF change was used to define lesion categories as decreasing, varying, increasing and invariant for MWF variation. ResultsIn total n = 386 T2 lesions were detected with a subset of n = 225 permanent lesions present at all four time-points. Among those, a heterogeneous lesion MWF reduction was found, with the majority of lesions bearing only mild MWF reduction, approximately a third with an intermediate MWF decrease and highest MWF reduction in acute-inflammatory active lesions. A moderate negative correlation was determined between individual lesion volumes and median MWF consistent across all time-points. Permanent lesions featured variable temporal dynamics with the majority of varying MWF (58 %), however decreasing (16 %), increasing (15 %) and invariant (11 %) subgroups could be identified resembling demyelinating activity and post-demyelinating inactivity known from histopathology studies. Inflammatory-active enhancing lesions showed a distinct pattern of MWF reduction followed by partial recovery after 3 months. This was similar in new enhancing lesions and those with a non-enhancing precursor lesion. ConclusionThis work provides in vivo evidence for an individual evolution of early demyelinated MS lesions measured by means of MWF imaging. Our results support the hypothesis, that MS lesions undergo multiple demyelination and remyelination episodes in the early acute phase. The in vivo MRI surrogate of myelin turnover bears capacity as a novel biomarker to select and potentially monitor personalized MS treatment.

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. T1 and T2 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 T1 /T2 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.

Assessing the differential sensitivities of wave-CAIPI ViSTa myelin water fraction and magnetization transfer saturation for efficiently quantifying tissue damage in MS

BackgroundWave-CAIPI Visualization of Short Transverse relaxation time component (ViSTa) is a recently developed, short-T1-sensitized MRI method for fast quantification of myelin water fraction (MWF) in the human brain. It represents a promising technique for the evaluation of subtle, early signals of demyelination in the cerebral white matter of multiple sclerosis (MS) patients. Currently however, few studies exist that robustly assess the utility of ViSTa MWF measures of myelin compared to more conventional MRI measures of myelin in the brain of MS patients. Moreover, there are no previous studies evaluating the sensitivity of ViSTa MWF for the non-invasive detection of subtle tissue damage in both normal-appearing white matter (NAWM) and white matter lesions of MS patients. As a result, a central purpose of this study was to systematically evaluate the relationship between myelin sensitivity of T1-based ViSTa MWF mapping and a more generally recognized metric, Magnetization Transfer Saturation (MTsat), in healthy control and MS brain white matter. MethodsViSTa MWF and MTsat values were evaluated in automatically-classified normal appearing white matter (NAWM), white matter (WM) lesion tissue, cortical gray matter, and deep gray matter of 29 MS patients and 10 healthy controls using 3T MRI. MWF and MT sat were also assessed in a tract-specific manner using the Johns Hopkins University WM atlas. MRI-derived measures of cerebral myelin content were uniquely compared by employing non-normal distribution-specific measures of median, interquartile range and skewness. Separate analyses of variance were applied to test tissue-specific differences in MTsat and ViSTa MWF distribution metrics. Non-parametric tests were utilized when appropriate. All tests were corrected for multiple comparisons using the False Discovery Rate method at the level, α=0.05. ResultsDifferences in whole NAWM MS tissue damage were detected with a higher effect size when using ViSTa MWF (q = 0.0008; ƞ2 = 0.34) compared to MTsat (q = 0.02; ƞ2= 0.24). We also observed that, as a possible measure of WM pathology, ViSTa-derived NAWM MWF voxel distributions of MS subjects were consistently skewed towards lower MWF values, while MTsat voxel distributions showed reduced skewness values. We further identified tract-specific reductions in mean ViSTa MWF of MS patients compared to controls that were not observed with MTsat. However, MTsat (q = 1.4 × 10−21; ƞ2 = 0.88) displayed higher effect sizes when differentiating NAWM and MS lesion tissue. Using regression analysis at the group level, we identified a linear relationship between MTsat and ViSTa MWF in NAWM (R2 = 0.46; p = 7.8 × 10−4) lesions (R2 = 0.30; p = 0.004), and with all tissue types combined (R2 = 0.71; p = 8.4 × 10−45). The linear relationship was also observed in most of the WM tracts we investigated. ViSTa MWF in NAWM of MS patients correlated with both disease duration (p = 0.02; R2 = 0.27) and WM lesion volume (p = 0.002; R2 = 0.34). ConclusionBecause ViSTa MWF and MTsat metrics exhibit differential sensitivities to tissue damage in MS white matter, they can be collected in combination to provide an efficient, comprehensive measure of myelin water and macromolecular pool proton signals. These complementary measures may offer a more sensitive, non-invasive biopsy of early precursor signals in NAWM that occur prior to lesion formation. They may also aid in monitoring the efficacy of remyelination therapies.

Revisiting the T2 spectrum imaging inverse problem: Bayesian regularized non-negative least squares

Multi-echo T2 magnetic resonance images contain information about the distribution of T2 relaxation times of compartmentalized water, from which we can estimate relevant brain tissue properties such as the myelin water fraction (MWF). Regularized non-negative least squares (NNLS) is the tool of choice for estimating non-parametric T2 spectra. However, the estimation is ill-conditioned, sensitive to noise, and highly affected by the employed regularization weight. The purpose of this study is threefold: first, we want to underline that the apparently innocuous use of two alternative parameterizations for solving the inverse problem, which we called the standard and alternative regularization forms, leads to different solutions; second, to assess the performance of both parameterizations; and third, to propose a new Bayesian regularized NNLS method (BayesReg). The performance of BayesReg was compared with that of two conventional approaches (L-curve and Chi-square (X2) fitting) using both regularization forms. We generated a large dataset of synthetic data, acquired in vivo human brain data in healthy participants for conducting a scan-rescan analysis, and correlated the myelin content derived from histology with the MWF estimated from ex vivo data. Results from synthetic data indicate that BayesReg provides accurate MWF estimates, comparable to those from L-curve and X2, and with better overall stability across a wider signal-to-noise range. Notably, we obtained superior results by using the alternative regularization form. The correlations reported in this study are higher than those reported in previous studies employing the same ex vivo and histological data. In human brain data, the estimated maps from L-curve and BayesReg were more reproducible. However, the T2 spectra produced by BayesReg were less affected by over-smoothing than those from L-curve. These findings suggest that BayesReg is a good alternative for estimating T2 distributions and MWF maps.

Temperature dependence and histological correlation of inhomogeneous magnetization transfer and myelin water imaging in ex vivo brain

PurposeThe promise of inhomogeneous magnetization transfer (ihMT) as a new myelin imaging method was studied in ex vivo human brain tissue and in relation to myelin water fraction (MWF). The temperature dependence of both methods was characterized, as well as their correspondence with a histological measure of myelin content. Unfiltered and filtered ihMT protocols were studied by adjusting the saturation scheme to preserve or attenuate signal from tissue with short dipolar relaxation time T1D. MethodsihMT ratio (ihMTR) and MWF maps were acquired at 7 T from formalin-fixed human brain samples at 22.5 °C, 30 °C and 37 °C. The impact of temperature on unfiltered ihMTR, filtered ihMTR and MWF was investigated and compared to myelin basic protein staining. ResultsUnfiltered ihMTR exhibited no temperature dependence, whereas filtered ihMTR increased with increasing temperature. MWF decreased at higher temperature, with an increasing prevalence of areas where the myelin water signal was unreliably determined, likely related to a reduction in T2 peak separability at higher temperatures ex vivo. MWF and ihMTR showed similar per-sample correlation with myelin staining at room temperature. At 37 °C, filtered ihMTR was more strongly correlated with myelin staining and had increased dynamic range compared to unfiltered ihMTR. ConclusionsGiven the temperature dependence of filtered ihMT, increased dynamic range, and strong myelin specificity that persists at higher temperatures, we recommend carefully controlled temperatures close to 37 °C for filtered ihMT acquisitions. Unfiltered ihMT may also be useful, due to its independence from temperature, higher amplitude values, and sensitivity to short T1D components. Ex vivo myelin water imaging should be performed at room temperature, to avoid fitting issues found at higher temperatures.

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.