Inherent and unpredictable bias in multi-component DESPOT myelin water fraction estimation

Quantifying myelin in neonates using magnetic resonance imaging: a systematic literature review

This review aimed to assess the usefulness of various magnetic resonance imaging (MRI) techniques for the quantification of neonatal white matter myelination. The Scopus, PubMed, and Web of Science databases were searched to identify studies following the PRISMA (preferred reporting items for systematic reviews and meta-analyses) statement using quantitative MRI techniques to examine samples collected from neonates to quantify myelin. Twelve studies were ultimately included. The results demonstrated that in validation studies, relaxometry is the most frequently explored approach (83.33%), followed by magnetization transfer imaging (8.33%) and a new automatic segmentation technique (8.33%). Synthetic MRI is recommended for quantifying myelin in neonates because of several advantages that outweigh a few negligible limitations.

The myelin water imaging transcriptome: myelin water fraction regionally varies with oligodendrocyte-specific gene expression

Identifying sensitive and specific measures that can quantify myelin are instrumental in characterizing microstructural changes in neurological conditions. Neuroimaging transcriptomics is emerging as a valuable technique in this regard, offering insights into the molecular basis of promising candidates for myelin quantification, such as myelin water fraction (MWF). We aimed to demonstrate the utility of neuroimaging transcriptomics by validating MWF as a myelin measure. We utilized data from a normative MWF brain atlas, comprised of 50 healthy subjects (mean age = 25 years, range = 17-42 years) scanned at 3 Tesla. Magnetic resonance imaging data included myelin water imaging to extract MWF and T1 anatomical scans for image registration and segmentation. We investigated the inter-regional distributions of gene expression data from the Allen Human Brain Atlas in conjunction with inter-regional MWF distribution patterns. Pearson correlations were used to identify genes with expression profiles mirroring MWF. The Single Cell Type Atlas from the Human Protein Atlas was leveraged to classify genes into gene sets with high cell type specificity, and a control gene set with low cell type specificity. Then, we compared the Pearson correlation coefficients for each gene set to determine if cell type-specific gene expression signatures correlate with MWF. Pearson correlation coefficients between MWF and gene expression for oligodendrocytes and adipocytes were significantly higher than for the control gene set, whereas correlations between MWF and inhibitory/excitatory neurons were significantly lower. Our approach in integrating transcriptomics with neuroimaging measures supports an emerging technique for understanding and validating MRI-derived markers such as MWF.

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.

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

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

Inhomogeneous magnetization transfer imaging: Concepts and directions for further development

Off-resonance radio frequency irradiation can induce the ordering of proton spins in the dipolar fields of their neighbors, in molecules with restricted mobility. This dipolar order decays with a characteristic relaxation time, T_1D , that is very different from the T₁ and T₂ relaxation of the nuclear alignment with the main magnetic field. Inhomogeneous magnetization transfer (ihMT) imaging is a refinement of magnetization transfer (MT) imaging that isolates the MT signal dependence on dipolar order relaxation times within motion-constrained molecules. Because T_1D relaxation is a unique contrast mechanism, ihMT may enable improved characterization of tissue. Initial work has stressed the high correlation between ihMT signal and myelin density. Dipolar order relaxation appears to be much longer in membrane lipids than other molecules. Recent work has shown, however, that ihMT acquisitions may also be adjusted to emphasize different ranges of T_1D . These newer approaches may be sensitive to other microstructural components of tissue. Here, we review the concepts and history of ihMT and outline the requirements for further development to realize its full potential.

Quantitative MRI in leukodystrophies

Leukodystrophies constitute a large and heterogeneous group of genetic diseases primarily affecting the white matter of the central nervous system. Different disorders target different white matter structural components. Leukodystrophies are most often progressive and fatal. In recent years, novel therapies are emerging and for an increasing number of leukodystrophies trials are being developed. Objective and quantitative metrics are needed to serve as outcome measures in trials. Quantitative MRI yields information on microstructural properties, such as myelin or axonal content and condition, and on the chemical composition of white matter, in a noninvasive fashion. By providing information on white matter microstructural involvement, quantitative MRI may contribute to the evaluation and monitoring of leukodystrophies. Many distinct MR techniques are available at different stages of development. While some are already clinically applicable, others are less far developed and have only or mainly been applied in healthy subjects. In this review, we explore the background, current status, potential and challenges of available quantitative MR techniques in the context of leukodystrophies.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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

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

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

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

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

High-resolution quantitative MRI of multiple sclerosis spinal cord lesions

PURPOSE

Validation of quantitative MR measures for myelin imaging in the postmortem multiple sclerosis spinal cord.

METHODS

Four fixed spinal cord samples were imaged first with a 3T clinical MR scanner to identify areas of interest for scanning, and then with a 7T small bore scanner using a multicomponent-driven equilibrium single-pulse observation of T₁ and T₂ protocol to produce apparent proton density, T₁ , T₂ , myelin water, intracellular water, and free-water fraction maps. After imaging, the cords were sectioned and stained with histological markers (hematoxylin and eosin, myelin basic protein, and neurofilament protein), which were quantitatively compared with the MR maps.

RESULTS

Excellent correspondence was found between high-resolution MR parameter maps and histology, particularly for apparent proton density MRI and myelin basic protein staining.

CONCLUSION

High-resolution quantitative MRI of the spinal cord provides biologically meaningful measures, and could be beneficial to diagnose and track multiple sclerosis lesions in the spinal cord.

Data-driven algorithm for myelin water imaging: Probing subvoxel compartmentation based on identification of spatially global tissue features

PURPOSE

Multicomponent analysis of MRI T₂ relaxation time (mcT₂ ) is commonly used for estimating myelin content by separating the signal at each voxel into its underlying distribution of T₂ values. This voxel-based approach is challenging due to the large ambiguity in the multi-T₂ space and the low SNR of MRI signals. Herein, we present a data-driven mcT₂ analysis, which utilizes the statistical strength of identifying spatially global mcT₂ motifs in white matter segments before deconvolving the local signal at each voxel.

METHODS

Deconvolution is done using a tailored optimization scheme, which incorporates the global mcT₂ motifs without additional prior assumptions regarding the number of microscopic components. The end results of this process are voxel-wise myelin water fraction maps.

RESULTS

Validations are shown for computer-generated signals, uniquely designed subvoxel mcT₂ phantoms, and in vivo human brain. Results demonstrated excellent fitting accuracy, both for the numerical and the physical mcT₂ phantoms, exhibiting excellent agreement between calculated myelin water fraction and ground truth. Proof-of-concept in vivo validation is done by calculating myelin water fraction maps for white matter segments of the human brain. Interscan stability of myelin water fraction values was also estimated, showing good correlation between scans.

CONCLUSION

We conclude that studying global tissue motifs prior to performing voxel-wise mcT₂ analysis stabilizes the optimization scheme and efficiently overcomes the ambiguity in the T₂ space. This new approach can improve myelin water imaging and the investigation of microstructural compartmentation in general.

Investigation of the association between cerebral iron content and myelin content in normative aging using quantitative magnetic resonance neuroimaging

Myelin loss and iron accumulation are cardinal features of aging and various neurodegenerative diseases. Oligodendrocytes incorporate iron as a metabolic substrate for myelin synthesis and maintenance. An emerging hypothesis in Alzheimer’s disease research suggests that myelin breakdown releases substantial stores of iron that may accumulate, leading to further myelin breakdown and neurodegeneration. We assessed associations between iron content and myelin content in critical brain regions using quantitative magnetic resonance imaging (MRI) on a cohort of cognitively unimpaired adults ranging in age from 21 to 94 years. We measured whole-brain myelin water fraction (MWF), a surrogate of myelin content, using multicomponent relaxometry, and whole-brain iron content using susceptibility weighted imaging in all individuals. MWF was negatively associated with iron content in most brain regions evaluated indicating that lower myelin content corresponds to higher iron content. Moreover, iron content was significantly higher with advanced age in most structures, with men exhibiting a trend towards higher iron content as compared to women. Finally, relationship between MWF and age, in all brain regions investigated, suggests that brain myelination continues until middle age, followed by degeneration at older ages. This work establishes a foundation for further investigations of the etiology and sequelae of myelin breakdown and iron accumulation in neurodegeneration and may lead to new imaging markers for disease progression and treatment.

Myelination may be impaired in neonates following birth asphyxia

•

Myelination is a developmental process that intensifies after birth during the first years of life.

•

We used a T2* mapping sequence to assess myelination in healthy and critically ill neonates with neonatal encephalopathy.

•

Birth asphyxia, in addition to causing the previously well-described direct injury to the brain, may impair myelination.

Background

Myelination is a developmental process that begins during the end of gestation, intensifies after birth over the first years of life, and continues well into adolescence. Any event leading to brain injury around the time of birth and during the perinatal period, such as birth asphyxia, may impair this critical process. Currently, the impact of such brain injury related to birth asphyxia on the myelination process is unknown.

Objective

To assess the myelination pattern over the first month of life in neonates with neonatal encephalopathy (NE) developing brain injury, compared to neonates without injury (i.e., healthy neonates and neonates with NE who do not develop brain injury).

Methods

Brain magnetic resonance imaging (MRI) was performed around day of life 2, 10, and 30 in healthy neonates and near-term/term neonates with NE who were treated with hypothermia. We evaluated myelination in various regions of interest using a T2* mapping sequence. In each region of interest, we compared the T2* values of the neonates with NE with brain injury to the values of the neonates without injury, according to the MRI timing, by using a repeated measures generalized linear mixed model.

Results

We obtained 74 MRI scans over the first month of life for 6 healthy neonates, 17 neonates with NE who were treated with hypothermia and did not develop brain injury, and 16 neonates with NE who were treated with hypothermia and developed brain injury. The T2* values significantly increased in the neonates with NE who developed injury in the posterior limbs of the internal capsule (day 2: p < 0.001; day 10: p < 0.001; and day 30: p < 0.001), the thalami (day 2: p = 0.001; day 10: p = 0.006; and day 30: p = 0.016), the lentiform nuclei (day 2: p = 0.005), the anterior white matter (day 2: p = 0.002; day 10: p = 0.006; and day 30: p = 0.002), the posterior white matter (day 2: p = 0.001; day 10: p = 0.008; and day 30: p = 0.03), the genu of the corpus callosum (day 2: p = 0.01; and day 10: p = 0.006), and the optic radiations (day 30: p < 0.001).

Conclusion

In the neonates with NE who were treated with hypothermia and developed brain injury, birth asphyxia impaired myelination in the regions that are myelinated at birth or soon after birth (the posterior limbs of internal capsule, the thalami, and the lentiform nuclei), in the regions where the myelination process begins only after the perinatal period (optic radiations), and in the regions where this process does not occur until months after birth (anterior/posterior white matter), which suggests that birth asphyxia, in addition to causing the previously well-described direct injury to the brain, may impair myelination.

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

PURPOSE

The 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 T_1D.

METHODS

ihMT 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.

RESULTS

Unfiltered 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 T₂ 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.

CONCLUSIONS

Given 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 T_1D components. Ex vivo myelin water imaging should be performed at room temperature, to avoid fitting issues found at higher temperatures.

Decreased myelin content of the fornix predicts poorer memory performance beyond vascular risk, hippocampal volume, and fractional anisotropy in nondemented older adults

Alterations to cerebral white matter tracts have been associated with cognitive decline in aging and Alzheimer's disease (AD). In particular, the fornix has been implicated as especially vulnerable given that it represents the primary outflow tract of the hippocampus. Despite this, little work has focused on the fornix using a potential early marker of white matter degeneration-myelin water fraction (MWF; an in vivo marker of myelin content). Therefore, we sought to (1) clarify associations between MWF in the fornix and memory functioning, and (2) examine whether fornix MWF relates to memory performance above and beyond hippocampal volume and conventional imaging measures of white matter that may not be as specific to alterations in myelin content. Forty nondemented older adults (mean age = 72.9 years) underwent an MRI exam and neuropsychological assessment. Multicomponent driven equilibrium single pulse observation of T1 and T2 (mcDESPOT) was used to quantify fornix MWF and diffusion tensor imaging (DTI) was used to measure fornix fractional anisotropy (FA). Adjusting for age, sex, education, and vascular risk factors, linear regression models revealed that, lower fornix MWF was significantly associated with poorer memory functioning (β = 0.405, p = .007) across our sample of older adults. Notably, fornix MWF remained a significant predictor of memory functioning (β = 0.380, p = .015) even after adjusting for fornix DTI FA and hippocampal volume (in addition to the above covariates). Given the observed associations between myelin and memory in older adults without dementia, MWF may be a useful early marker of dementia risk.

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

BACKGROUND

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

NEW METHOD

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

RESULTS

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

COMPARISON WITH EXISTING METHODS

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

CONCLUSIONS

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

Comparison of multi echo T2 relaxation and steady state approaches for myelin imaging in the central nervous system

The traditional approach for measuring myelin-associated water with quantitative magnetic resonance imaging (MRI) uses multi-echo T₂ relaxation data to calculate the myelin water fraction (MWF). A fundamentally different approach, abbreviated “mcDESPOT”, uses a more efficient steady-state acquisition to generate an equivalent metric (f_M). Although previous studies have demonstrated inherent instability and bias in the complex mcDESPOT analysis procedure, f_M has often been used as a surrogate for MWF. We produced and compared multivariate atlases of MWF and f_M in healthy human brain and cervical spinal cord (available online) and compared their ability to detect multiple sclerosis pathology. A significant bias was found in all regions (p < 10^–5), albeit reversed for spinal cord (f_M-MWF =  − 3.4%) compared to brain (+ 6.2%). MWF and f_M followed an approximately linear relationship for regions with MWF <  ~ 10%. For MWF >  ~ 10%, the relationship broke down and f_M no longer increased in tandem with MWF. For multiple sclerosis patients, MWF and f_M Z score maps showed overlapping areas of low Z score and similar trends between patients and brain regions, although those of f_M generally had greater spatial extent and magnitude of severity. These results will guide future choice of myelin-sensitive quantitative MRI and improve interpretation of studies using either myelin imaging approach.

An atlas for human brain myelin content throughout the adult life span

Myelin water imaging is a quantitative neuroimaging technique that provides the myelin water fraction (MWF), a metric highly specific to myelin content, and the intra-/extra-cellular T₂ (IET2), which is related to water and iron content. We coupled high-resolution data from 100 adults with gold-standard methodology to create an optimized anatomical brain template and accompanying MWF and IET2 atlases. We then used the MWF atlas to characterize how myelin content relates to demographic factors. In most brain regions, myelin content followed a quadratic pattern of increase during the third decade of life, plateau at a maximum around the fifth decade, then decrease during later decades. The ranking of mean myelin content between brain regions remained consistent across age groups. These openly available normative atlases can facilitate evaluation of myelin imaging results on an individual basis and elucidate the distribution of myelin content between brain regions and in the context of aging.

Nonlinear associations of neurite density and myelin content with age revealed using multicomponent diffusion and relaxometry magnetic resonance imaging

Most magnetic resonance imaging (MRI) studies investigating the relationship between regional brain myelination or axonal density and aging have relied upon nonspecific methods to probe myelin and axonal content, including diffusion tensor imaging and relaxation time mapping. While these studies have provided pivotal insights into changes in cerebral architecture with aging and pathology, details of the underlying microstructural alterations have not been fully elucidated. In the current study, we used the BMC-mcDESPOT analysis, a direct and specific multicomponent relaxometry method for imaging of myelin water fraction (MWF), a marker of myelin content, and NODDI, an emerging multicomponent diffusion technique, for neurite density index (NDI) imaging, a proxy of axonal density. We investigated age-related differences in MWF and NDI in several white matter brain regions in a cohort of cognitively unimpaired participants over a wide age range. Our results indicate a quadratic, inverted U-shape, relationship between MWF and age in all brain regions investigated, suggesting that myelination continues until middle age followed by a decrease at older ages, in agreement with previous work. We found a similarly complex regional association between NDI and age, with several cerebral structures also exhibiting a quadratic, inverted U-shape, relationship. This novel observation suggests an increase in axonal density until the fourth decade of age followed by a rapid loss at older ages. We also observed that these age-related differences in MWF and NDI vary across different brain regions, as expected. Finally, our study indicates no significant association between MWF and NDI in most cerebral structures investigated, although this association approached significance in a limited number of brain regions, indicating the complementary nature of their information and encouraging further investigation. Overall, we find evidence of nonlinear associations between age and myelin or axonal density in a sample of well-characterized adults, using direct myelin and axonal content imaging methods.

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

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

Myelin quantification with MRI: A systematic review of accuracy and reproducibility

OBJECTIVES

Currently, multiple sclerosis is treated with anti-inflammatory therapies, but these treatments lack efficacy in progressive disease. New treatment strategies aim to repair myelin damage and efficacy evaluation of such new therapies would benefit from validated myelin imaging techniques. Several MRI methods for quantification of myelin density are available now. This systematic review aims to analyse the performance of these MRI methods.

METHODS

Studies comparing myelin quantification by MRI with histology, the current gold standard, or assessing reproducibility were retrieved from PubMed/MEDLINE and Embase (until December 2019). Included studies assessed both myelin histology and MRI quantitatively. Correlation or variance measurements were extracted from the studies. Non-parametric tests were used to analyse differences in study methodologies.

RESULTS

The search yielded 1348 unique articles. Twenty-two animal studies and 13 human studies correlated myelin MRI with histology. Eighteen clinical studies analysed the reproducibility. Overall bias risk was low or unclear. All MRI methods performed comparably, with a mean correlation between MRI and histology of R²=0.54 (SD=0.30) for animal studies, and R²=0.54 (SD=0.18) for human studies. Reproducibility for the MRI methods was good (ICC=0.75-0.93, R²=0.90-0.98, COV=1.3-27%), except for MTR (ICC=0.05-0.51).

CONCLUSIONS

Overall, MRI-based myelin imaging methods show a fairly good correlation with histology and a good reproducibility. However, the amount of validation data is too limited and the variability in performance between studies is too large to select the optimal MRI method for myelin quantification yet.

Myelin water fraction estimation using small-tip fast recovery MRI

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Microstructure of Human Corpus Callosum across the Lifespan: Regional Variations in Axon Caliber, Density, and Myelin Content

The myeloarchitecture of the corpus callosum (CC) is characterized as a mosaic of distinct differences in fiber density of small- and large-diameter axons along the anterior-posterior axis; however, regional and age differences across the lifespan are not fully understood. Using multiecho T2 magnetic resonance imaging combined with multi-T2 fitting, the myelin water fraction (MWF) and geometric-mean of the intra-/extracellular water T2 (geomT2IEW) in 395 individuals (7-85 years; 41% males) were examined. The approach was validated where regional patterns along the CC closely resembled the histology; MWF matched mean axon diameter and geomT2IEW mirrored the density of large-caliber axons. Across the lifespan, MWF exhibited a quadratic association with age in all 10 CC regions with evidence of a positive linear MWF-age relationship among younger participants and minimal age differences in the remainder of the lifespan. Regarding geomT2IEW, a significant linear age × region interaction reflected positive linear age dependence mostly prominent in the regions with the highest density of small-caliber fibers-genu and splenium. In all, these two indicators characterize distinct attributes that are consistent with histology, which is a first. In addition, these results conform to rapid developmental progression of CC myelination leveling in middle age as well as age-related degradation of axon sheaths in older adults.

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.

Longitudinal white matter and cognitive development in pediatric carriers of the apolipoprotein ε4 allele

We have previously demonstrated cross-sectional differences in magnetic resonance imaging (MRI) measurements of white matter myelin and gray matter in infants with or without the apolipoprotein ε4 allele, a major genetic risk factor for late-onset Alzheimer's disease (AD). In this study, we sought to compare longitudinal MRI white matter myelin and cognitive-behavioral changes in infants and young children with and without this allele. Serial MRI and cognitive tests were obtained on 223 infants and young children, including 74 ε4 carriers and 149 non-carriers, 2–68 months of age, matched for age, gestational duration, birth weight, sex ratio, maternal age, education, and socioeconomic status. Automated brain mapping algorithms and non-linear mixed models were used to characterize and compare trajectories of white matter myelin and cognitive-behavioral test scores. The APOE ε4 carriers had statistically significant differences in white matter myelin development, in the uncinate fasciculus, temporal lobe, internal capsule and occipital lobe. Additionally, ε4 carriers had a slightly greater rate of development in early learning composite a surrogate measure of IQ representative of expressive language, receptive language, fine motor, and visual skills, but displayed slightly lower non verbal development quotient scores a composite measure of fine motor and visual skills across the entire age range. This study supports the possibility that ε4 carriers have slightly altered rates of white matter and cognitive development in childhood. It continues to raise questions about the role of APOE in human brain development and the relevance of these developmental differences to the predisposition to AD.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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

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

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

PTSD, but not history of mTBI, is associated with altered myelin in combat-exposed Iraq and Afghanistan Veterans

OBJECTIVE

To investigate the biological, cognitive, and psychological presentations of combat-exposed Veterans with a history of mild traumatic brain injury (mTBI) and/or posttraumatic stress disorder (PTSD) using a novel white matter imaging technique and comprehensive neuropsychological assessment.

METHOD

74 Iraq and Afghanistan Veterans (mean age 33.89, 90.5% male) with history of mTBI (average 7.25 years since injury), PTSD, both, or neither underwent magnetic resonance imaging (MRI) exams including acquisition of a novel imaging technique, multicomponent-driven equilibrium single-pulse observation of T1/T2 (mcDESPOT) to quantify myelin water fraction (MWF), a surrogate measure of myelin content. Participants also underwent comprehensive neuropsychological assessment and three cognitive composite scores (memory, working memory/processing speed, and executive functioning) were created.

RESULTS

There were no significant group differences on the neuropsychological composite scores. ANCOVAs revealed a main effect of PTSD across all a priori regions of interest (ROI) in which PTSD was associated with higher MWF. There was no main effect of mTBI history or TBI by PTSD interaction on any ROI. Significant positive associations were observed between myelin and PTSD symptoms, but no significant associations were found between myelin and neurobehavioral symptoms. No significant associations were found between myelin in the a priori ROIs and the cognitive composite scores.

CONCLUSION

This study did not find neuropsychological or MWF differences in combat Veterans with a remote history of mTBI but did find myelin alterations related to PTSD. Psychological trauma should be a primary target for intervention in Veterans with comorbid PTSD and mTBI reporting subjective complaints, given its salience.

Brain Myelin Water Fraction and Diffusion Tensor Imaging Atlases for 9-10 Year-Old Children

BACKGROUND AND PURPOSE

Myelin water imaging (MWI) and diffusion tensor imaging (DTI) provide information about myelin and axon-related brain microstructure, which can be useful for investigating normal brain development and many childhood brain disorders. While pediatric DTI atlases exist, there are no pediatric MWI atlases available for the 9-10 years old age group. As myelination and structural development occurs throughout childhood and adolescence, studies of pediatric brain pathologies must use age-specific MWI and DTI healthy control data. We created atlases of myelin water fraction (MWF) and DTI metrics for healthy children aged 9-10 years for use as normative data in pediatric neuroimaging studies.

METHODS

3D-T₁ , DTI, and MWI scans were acquired from 20 healthy children (mean age: 9.6 years, range: 9.2-10.3 years, 4 females). ANTs and FSL registration were used to create quantitative MWF and DTI atlases. Region of interest (ROI) analysis in nine white matter regions was used to compare pediatric MWF with adult MWF values from a recent study and to investigate the correlation between pediatric MWF and DTI metrics.

RESULTS

Adults had significantly higher MWF than the pediatric cohort in seven of the nine white matter ROIs, but not in the genu of the corpus callosum or the cingulum. In the pediatric data, MWF correlated significantly with mean diffusivity, but not with axial diffusivity, radial diffusivity, or fractional anisotropy.

CONCLUSIONS

Normative MWF and DTI metrics from a group of 9-10 year old healthy children provide a resource for comparison to pathologies. The age-specific atlases are ready for use in pediatric neuroimaging research and can be accessed: https://sourceforge.net/projects/pediatric-mri-myelin-diffusion/.

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.

Estimating axon conduction velocity in vivo from microstructural MRI

The conduction velocity (CV) of action potentials along axons is a key neurophysiological property central to neural communication. The ability to estimate CV in humans in vivo from non-invasive MRI methods would therefore represent a significant advance in neuroscience. However, there are two major challenges that this paper aims to address: (1) Much of the complexity of the neurophysiology of action potentials cannot be captured with currently available MRI techniques. Therefore, we seek to establish the variability in CV that can be captured when predicting CV purely from parameters that have been reported to be estimatable from MRI: inner axon diameter (AD) and g-ratio. (2) errors inherent in existing MRI-based biophysical models of tissue will propagate through to estimates of CV, the extent to which is currently unknown. Issue (1) is investigated by performing a sensitivity analysis on a comprehensive model of axon electrophysiology and determining the relative sensitivity to various morphological and electrical parameters. The investigations suggest that 85% of the variance in CV is accounted for by variation in AD and g-ratio. The observed dependency of CV on AD and g-ratio is well characterised by the previously reported model by Rushton. Issue (2) is investigated through simulation of diffusion and relaxometry MRI data for a range of axon morphologies, applying models of restricted diffusion and relaxation processes to derive estimates of axon volume fraction (AVF), AD and g-ratio and estimating CV from the derived parameters. The results show that errors in the AVF have the biggest detrimental impact on estimates of CV, particularly for sparse fibre populations (AVF<0.3). For our equipment set-up and acquisition protocol, CV estimates are most accurate (below 5% error) where AVF is above 0.3, g-ratio is between 0.6 and 0.85 and AD is high (above 4μm). CV estimates are robust to errors in g-ratio estimation but are highly sensitive to errors in AD estimation, particularly where ADs are small. We additionally show CV estimates in human corpus callosum in a small number of subjects. In conclusion, we demonstrate accurate CV estimates are possible in regions of the brain where AD is sufficiently large. Problems with estimating ADs for smaller axons presents a problem for estimating CV across the whole CNS and should be the focus of further study.

Quantitative age-dependent differences in human brainstem myelination assessed using high-resolution magnetic resonance mapping

Previous in-vivo magnetic resonance imaging (MRI)-based studies of age-related differences in the human brainstem have focused on volumetric morphometry. These investigations have provided pivotal insights into regional brainstem atrophy but have not addressed microstructural age differences. However, growing evidence indicates the sensitivity of quantitative MRI to microstructural tissue changes in the brain. These studies have largely focused on the cerebrum, with very few MR investigations addressing age-dependent differences in the brainstem, in spite of its central role in the regulation of vital functions. Several studies indicate early brainstem alterations in a myriad of neurodegenerative diseases and dementias. The paucity of MR-focused investigations is likely due in part to the challenges imposed by the small structural scale of the brainstem itself as well as of substructures within, requiring accurate high spatial resolution imaging studies. In this work, we applied our recently developed approach to high-resolution myelin water fraction (MWF) mapping, a proxy for myelin content, to investigate myelin differences with normal aging within the brainstem. In this cross-sectional investigation, we studied a large cohort (n = 125) of cognitively unimpaired participants spanning a wide age range (21-94 years) and found a decrease in myelination with age in most brainstem regions studied, with several regions exhibiting a quadratic association between myelin and age. We believe that this study is the first investigation of MWF differences with normative aging in the adult brainstem. Further, our results provide reference MWF values.

Adult brain aging investigated using BMC-mcDESPOT-based myelin water fraction imaging

The relationship between regional brain myelination and aging has been the subject of intense study, with magnetic resonance imaging perhaps the most effective modality for elucidating this. However, most of these studies have used nonspecific methods to probe myelin content, including diffusion tensor imaging, magnetization transfer ratio, and relaxation times. In the present study, we used the BMC-mcDESPOT analysis, a direct and specific method for imaging of myelin water fraction (MWF), a surrogate of myelin content. We investigated age-related differences in MWF in several brain regions in a large cohort of cognitively unimpaired participants, spanning a wide age range. Our results indicate a quadratic, inverted U-shape, relationship between MWF and age in all brain regions investigated, suggesting that myelination continues until middle age followed by decreases at older ages. We also observed that these age-related differences vary across different brain regions, as expected. Our results provide evidence for nonlinear associations between age and myelin in a large sample of well-characterized adults, using a direct myelin content imaging method.

Myelin Water Fraction and Intra/Extracellular Water Geometric Mean T2 Normative Atlases for the Cervical Spinal Cord from 3T MRI

BACKGROUND AND PURPOSE

Acquiring and interpreting quantitative myelin-specific MRI data at an individual level is challenging because of technical difficulties and natural myelin variation in the population. To overcome these challenges, we used multiecho T₂ myelin water imaging (MWI) to create T₂ metric healthy population atlases that depict the mean and variation of myelin water fraction (MWF), and intra- and extracellular water mobility as described by geometric mean T₂ (IEGMT₂ ).

METHODS

Cervical cord MWI was performed at 3T on 20 healthy individuals (10M/10F, mean age: 36 years) and 3 relapsing remitting multiple sclerosis (RRMS) participants (1M/2F, age: 39/42/37 years). Anatomical data were collected for the purpose of image segmentation and registration. Atlases were created by coregistering and averaging T₂ metrics from all controls. Voxel-wise z-score maps from 3 RRMS participants were produced to demonstrate the preliminary utility of the MWF and IEGMT₂ atlases.

RESULTS

The average MWF atlas provides a representation of myelin in the spinal cord consistent with well-known spinal cord anatomical characteristics. The IEGMT₂ atlas also depicted structural variations in the spinal cord. Z-score analysis illustrated distinct abnormalities in MWF and IEGMT₂ in the 3 RRMS cases.

CONCLUSIONS

Our findings highlight the potential for using a quantitative T₂ relaxation metric atlas to visualize and detect pathology in spinal cord. Our MWF and IEGMT₂ atlases (URL: https://sourceforge.net/projects/mwi-spinal-cord-atlases/) can serve as normative references in the cervical spinal cord for other studies.