Myelin volume fraction imaging with MRI

Permittivity tensor imaging: modular label-free imaging of 3D dry mass and 3D orientation at high resolution

The dry mass and the orientation of biomolecules can be imaged without a label by measuring their permittivity tensor (PT), which describes how biomolecules affect the phase and polarization of light. Three-dimensional (3D) imaging of PT has been challenging. We present a label-free computational microscopy technique, PT imaging (PTI), for the 3D measurement of PT. PTI encodes the invisible PT into images using oblique illumination, polarization-sensitive detection and volumetric sampling. PT is decoded from the data with a vectorial imaging model and a multi-channel inverse algorithm, assuming uniaxial symmetry in each voxel. We demonstrate high-resolution imaging of PT of isotropic beads, anisotropic glass targets, mouse brain tissue, infected cells and histology slides. PTI outperforms previous label-free imaging techniques such as vector tomography, ptychography and light-field imaging in resolving the 3D orientation and symmetry of organelles, cells and tissue. We provide open-source software and modular hardware to enable the adoption of the method.

Insights and improvements in correspondence between axonal volume fraction measured with diffusion-weighted MRI and electron microscopy

Biophysical diffusion-weighted imaging (DWI) models are increasingly used in neuroscience to estimate the axonal water fraction (f AW ), which in turn is key for noninvasive estimation of the axonal volume fraction (f A ). These models require thorough validation by comparison with a reference method, for example, electron microscopy (EM). While EM studies often neglect the unmyelinated axons and solely report the fraction of myelinated axons, in DWI both myelinated and unmyelinated axons contribute to the DWI signal. However, DWI models often include simplifications, for example, the neglect of differences in the compartmental relaxation times or fixed diffusivities, which in turn might affect the estimation off AW . We investigate whether linear calibration parameters (scaling and offset) can improve the comparability between EM- and DWI-based metrics off A . To this end, we (a) used six DWI models based on the so-called standard model of white matter (WM), including two models with fixed compartmental diffusivities (e.g., neurite orientation dispersion and density imaging, NODDI) and four models that fitted the compartmental diffusivities (e.g., white matter tract integrity, WMTI), and (b) used a multimodal data set including ex vivo diffusion DWI and EM data in mice with a broad dynamic range of fibre volume metrics. We demonstrated that the offset is associated with the volume fraction of unmyelinated axons and the scaling factor is associated with different compartmentalT 2 and can substantially enhance the comparability between EM- and DWI-based metrics off A . We found that DWI models that fitted compartmental diffusivities provided the most accurate estimates of the EM-basedf A . Finally, we introduced a more efficient hybrid calibration approach, where only the offset is estimated but the scaling is fixed to a theoretically predicted value. Using this approach, a similar one-to-one correspondence to EM was achieved for WMTI. The method presented can pave the way for use of validated DWI-based models in clinical research and neuroscience.

Characterization of hyperpolarization-activated cyclic nucleotide-gated channels in oligodendrocytes

Mature oligodendrocytes (OLG) are the myelin-forming cells of the central nervous system. Recent work has shown a dynamic role for these cells in the plasticity of neural circuits, leading to a renewed interest in voltage-sensitive currents in OLG. Hyperpolarization-activated cyclic nucleotide-gated (HCN) channels and their respective current (Ih) were recently identified in mature OLG and shown to play a role in regulating myelin length. Here we provide a biochemical and electrophysiological characterization of HCN channels in cells of the oligodendrocyte lineage. We observed that mice with a nonsense mutation in the Hcn2 gene (Hcn2ap/ap) have less white matter than their wild type counterparts with fewer OLG and fewer oligodendrocyte progenitor cells (OPCs). Hcn2ap/ap mice have severe motor impairments, although these deficits were not observed in mice with HCN2 conditionally eliminated only in oligodendrocytes (Cnpcre/+; Hcn2F/F). However, Cnpcre/+; Hcn2F/F mice develop motor impairments more rapidly in response to experimental autoimmune encephalomyelitis (EAE). We conclude that HCN2 channels in OLG may play a role in regulating metabolism.

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/synthmap Reduction 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.

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.

Empirical relationship between TEM-derived myelin volume fraction and MRI-R2 values in aging ex vivo rat corpus callosum

Ex vivo ratiometric measurements of short- and long-T2 components using the multiple spin echo sequence of MRI are often employed to evaluate alterations in myelin content in the white matter (WM) of the brain. However, the relationship between absolute MRI-T2 values (long-T2 component) and myelin volumetric information in aged ex vivo rodent WM appears to be influenced by factors such as animal species, field strength, and fixation durations/washing. Here, multiple spin echo sequence-based MRI-R2 (the reciprocal of T2) values were measured in the corpus callosum (CC) region in the post-mortem rat brains (n = 9) of different age groups with common fixation techniques without washing at 7 T. Transmission electron microscopy (TEM)-based quantification of myelin volume fraction (MVF) and corresponding Monte-Carlo simulation to estimate relaxation rates (R2,IE) due to diffusion in the presence of inhomogeneous magnetic field perturbation in intra- and extra-cellular (IE) spaces were respectively performed. To determine whether the short-T2 components originating from myelin water were mixed with long-T2 components from IE water or were undetectable, the MVF values obtained from TEM results were respectively compared with MRI-R2 and R2,IE values. A significant correlation (Pearson's correlation coefficient r = 0.8763; p < 0.01) of average MRI-R2 and MVF values was observed. Estimated R2,IE values from Monte-Carlo simulations in IE water signals were also positively correlated (r = 0.8281; p < 0.01) with MVF values. However, the magnitudes of R2,IE values were much smaller than those observed for MRI-R2 values, indicating that changes in R2 related MVF are likely dominated by myelin water components. Such comparisons between independent parameters from MRI, TEM, and simulations support the suggestion that myelin water signals were indistinguishably mixed to exhibit mono-exponential T2 relaxation, and multiple spin echo sequence-based MRI-R2 values in aging ex vivo rat CC without prolonged washing still reflect the volumetric information of myelin, likely due to enhanced water exchange across the myelin.

Quantitative MRI reveals widespread, network-specific myelination change during generalized epilepsy progression

Activity-dependent myelination is a fundamental mode of brain plasticity which significantly influences network function. We recently discovered that absence seizures, which occur in multiple forms of generalized epilepsy, can induce activity-dependent myelination, which in turn promotes further progression of epilepsy. Structural alterations of myelin are likely to be widespread, given that absence seizures arise from an extensive thalamocortical network involving frontoparietal regions of the bilateral hemispheres. However, the temporal course and spatial extent of myelin plasticity is unknown, due to limitations of gold-standard histological methods such as electron microscopy (EM). In this study, we leveraged magnetization transfer and diffusion MRI for estimation of g-ratios across major white matter tracts in a mouse model of generalized epilepsy with progressive absence seizures. EM was performed on the same brains after MRI. After seizure progression, we found increased myelination (decreased g-ratios) throughout the anterior portion (genu-to-body) of the corpus callosum but not in the posterior portion (body-splenium) nor in the fornix or the internal capsule. Curves obtained from averaging g-ratio values at every longitudinal point of the corpus callosum were statistically different with p<0.001. Seizure-associated myelin differences found in the corpus callosum body with MRI were statistically significant (p = 0.0027) and were concordant with EM in the same region (p = 0.01). Notably, these differences were not detected by diffusion tensor imaging. This study reveals widespread myelin structural change that is specific to the absence seizure network. Furthermore, our findings demonstrate the potential utility and importance of MRI-based g-ratio estimation to non-invasively detect myelin plasticity.

Associations between myelin water imaging and measures of fall risk and functional mobility in multiple sclerosis

BACKGROUND AND PURPOSE

Myelin water fraction (MWF) deficits as measured by myelin water imaging (MWI) have been related to worse motor function in persons with multiple sclerosis (PwMS). However, it is unknown if measures from MWI metrics in motor areas relate to fall risk measures in PwMS. The objective of this study was to examine the relationship between MWI measures in motor areas to performance on clinical measures of fall risk and disability in PwMS.

METHODS

Sixteen individuals with relapsing-remitting MS participated (1 male, 15 female; age 47.1 years [12.3]; Expanded Disability Status Scale 4.0 [range 0-6.5]) and completed measures of walking and fall risk (Timed 25 Foot Walk [T25FW] and Timed Up and Go). MWF and the geometric mean of the intra-/extracellular water T₂ (geomT_2IEW ) values reflecting myelin content and contribution of large-diameter axons/density, respectively, were assessed in three motor-related regions.

RESULTS

The geomT_2IEW of the corticospinal tract (r = -.599; p = .018) and superior cerebellar peduncles (r = -.613; p = .015) demonstrated significant inverse relationships with T25FW, suggesting that decreased geomT_2IEW was related to slower walking. Though not significant, MWF in the corticospinal tract and superior cerebellar peduncles also demonstrated fair relationships with the T25FW, suggesting that worse performance on the T25FW was associated with lower MWF values.

CONCLUSIONS

MWI of key motor regions was associated with walking performance in PwMS. Further MWI studies are needed to identify relationships between pathology and clinical function in PwMS to guide targeted rehabilitation therapies aimed at preventing falls.

Comparing myelin-sensitive magnetic resonance imaging measures and resulting g-ratios in healthy and multiple sclerosis brains

The myelin concentration and the degree of myelination of nerve fibers can provide valuable information on the integrity of human brain tissue. Magnetic resonance imaging (MRI) of myelin-sensitive parameters can help to non-invasively evaluate demyelinating diseases such as multiple sclerosis (MS). Several different myelin-sensitive MRI methods have been proposed to determine measures of the degree of myelination, in particular the g-ratio. However, variability in underlying physical principles and different biological models influence measured myelin concentrations, and consequently g-ratio values. We therefore investigated similarities and differences between five different myelin-sensitive MRI measures and their effects on g-ratio mapping in the brains of both MS patients and healthy volunteers. We compared two different estimates of the myelin water fraction (MWF) as well as the inhomogeneous magnetization transfer ratio (ihMTR), magnetization transfer saturation (MTsat), and macromolecular tissue volume (MTV) in 13 patients with MS and 14 healthy controls. In combination with diffusion-weighted imaging, we derived g-ratio parameter maps for each of the five different myelin measures. The g-ratio values calculated from different myelin measures varied strongly, especially in MS lesions. While, compared to normal-appearing white matter, MTsat and one estimate of the MWF resulted in higher g-ratio values within lesions, ihMTR, MTV, and the second MWF estimate resulted in lower lesion g-ratio values. As myelin-sensitive measures provide rough estimates of myelin content rather than absolute myelin concentrations, resulting g-ratio values strongly depend on the utilized myelin measure and model used for g-ratio mapping. When comparing g-ratio values, it is, thus, important to utilize the same MRI methods and models or to consider methodological differences. Particular caution is necessary in pathological tissue such as MS lesions.

Brain myelination at 7 months of age predicts later language development

Between 6 and 12 months of age there are dramatic changes in infants' processing of language. The neurostructural underpinnings of these changes are virtually unknown. The objectives of this study were to (1) examine changes in brain myelination during this developmental period and (2) examine the relationship between myelination during this period and later language development. Macromolecular proton fraction (MPF) was used as a marker of myelination. Whole-brain MPF maps were obtained with 1.25 mm3 isotropic spatial resolution from typically developing children at 7 and 11 months of age. Effective myelin density was calculated from MPF based on a linear relationship known from the literature. Voxel-based analyses were used to identify longitudinal changes in myelin density and to calculate correlations between myelin density at these ages and later language development. Increases in myelin density were more predominant in white matter than in gray matter. A strong predictive relationship was found between myelin density at 7 months of age, language production at 24 and 30 months of age, and rate of language growth. No relationships were found between myelin density at 11 months, or change in myelin density between 7 and 11 months of age, and later language measures. Our findings suggest that critical changes in brain structure may precede periods of pronounced change in early language skills.

Evaluating brain damage in multiple sclerosis with simultaneous multi-angular-relaxometry of tissue

OBJECTIVE

Multiple sclerosis (MS) is a common demyelinating central nervous system disease. MRI methods that can quantify myelin loss are needed for trials of putative remyelinating agents. Quantitative magnetization transfer MRI introduced the macromolecule proton fraction (MPF), which correlates with myelin concentration. We developed an alternative approach, Simultaneous-Multi-Angular-Relaxometry-of-Tissue (SMART) MRI, to generate MPF. Our objective was to test SMART-derived MPF metric as a potential imaging biomarker of demyelination.

METHODS

Twenty healthy control (HC), 11 relapsing-remitting MS (RRMS), 22 progressive MS (PMS), and one subject with a biopsied tumefactive demyelinating lesion were scanned at 3T using SMART MRI. SMART-derived MPF metric was determined in normal-appearing cortical gray matter (NAGM), normal-appearing subcortical white matter (NAWM), and demyelinating lesions. MPF metric was evaluated for correlations with physical and cognitive test scores. Comparisons were made between HC and MS and between MS subtypes. Furthermore, correlations were determined between MPF and neuropathology in the biopsied person.

RESULTS

SMART-derived MPF in NAGM and NAWM were lower in MS than HC (p < 0.001). MPF in NAGM, NAWM and lesions differentiated RRMS from PMS (p < 0.01, p < 0.001, p < 0.001, respectively), whereas lesion volumes did not. MPF in NAGM, NAWM and lesions correlated with the Expanded Disability Status Scale (p < 0.01, p < 0.001, p < 0.001, respectively) and nine-hole peg test (p < 0.001, p < 0.001, p < 0.01, respectively). MPF was lower in the histopathologically confirmed inflammatory demyelinating lesion than the contralateral NAWM and increased in the biopsied lesion over time, mirroring improved clinical performance.

INTERPRETATION

SMART-derived MPF metric holds potential as a quantitative imaging biomarker of demyelination and remyelination.

Inconsistencies in atlas-based volumetric measures of the human nucleus basalis of Meynert: A need for high-resolution alternatives

The nucleus basalis of Meynert (nbM) is the major source of cortical acetylcholine (ACh) and has been related to cognitive processes and to neurological disorders. However, spatially delineating the human nbM in MRI studies remains challenging. Due to the absence of a functional localiser for the human nbM, studies to date have localised it using nearby neuroanatomical landmarks or using probabilistic atlases. To understand the feasibility of MRI of the nbM we set our four goals; our first goal was to review current human nbM region-of-interest (ROI) selection protocols used in MRI studies, which we found have reported highly variable nbM volume estimates. Our next goal was to quantify and discuss the limitations of existing atlas-based volumetry of nbM. We found that the identified ROI volume depends heavily on the atlas used and on the probabilistic threshold set. In addition, we found large disparities even for data/studies using the same atlas and threshold. To test whether spatial resolution contributes to volume variability, as our third goal, we developed a novel nbM mask based on the normalized BigBrain dataset. We found that as long as the spatial resolution of the target data was 1.3 mm isotropic or above, our novel nbM mask offered realistic and stable volume estimates. Finally, as our last goal we tried to discern nbM using publicly available and novel high resolution structural MRI ex vivo MRI datasets. We find that, using an optimised 9.4T quantitative T₂^⁎ ex vivo dataset, the nbM can be visualised using MRI. We conclude caution is needed when applying the current methods of mapping nbM, especially for high resolution MRI data. Direct imaging of the nbM appears feasible and would eliminate the problems we identify, although further development is required to allow such imaging using standard (f)MRI scanning.

Error quantification in multi-parameter mapping facilitates robust estimation and enhanced group level sensitivity

Multi-Parameter Mapping (MPM) is a comprehensive quantitative neuroimaging protocol that enables estimation of four physical parameters (longitudinal and effective transverse relaxation rates R₁ and R₂^*, proton density PD, and magnetization transfer saturation MT_sat) that are sensitive to microstructural tissue properties such as iron and myelin content. Their capability to reveal microstructural brain differences, however, is tightly bound to controlling random noise and artefacts (e.g. caused by head motion) in the signal. Here, we introduced a method to estimate the local error of PD, R₁, and MT_sat maps that captures both noise and artefacts on a routine basis without requiring additional data. To investigate the method's sensitivity to random noise, we calculated the model-based signal-to-noise ratio (mSNR) and showed in measurements and simulations that it correlated linearly with an experimental raw-image-based SNR map. We found that the mSNR varied with MPM protocols, magnetic field strength (3T vs. 7T) and MPM parameters: it halved from PD to R₁ and decreased from PD to MT_sat by a factor of 3-4. Exploring the artefact-sensitivity of the error maps, we generated robust MPM parameters using two successive acquisitions of each contrast and the acquisition-specific errors to down-weight erroneous regions. The resulting robust MPM parameters showed reduced variability at the group level as compared to their single-repeat or averaged counterparts. The error and mSNR maps may better inform power-calculations by accounting for local data quality variations across measurements. Code to compute the mSNR maps and robustly combined MPM maps is available in the open-source hMRI toolbox.

Bundle myelin fraction (BMF) mapping of different white matter connections using microstructure informed tractography

To date, we have scarce information about the relative myelination level of different fiber bundles in the human brain. Indirect evidence comes from postmortem histology data but histological stainings are unable to follow a specific bundle and determine its intrinsic myelination. In this context, quantitative MRI, and diffusion MRI tractography may offer a viable solution by providing, respectively, voxel-wise myelin sensitive maps and the pathways of the major tracts of the brain. Then, "tractometry" can be used to combine these two pieces of information by averaging tissue features (obtained from any voxel-wise map) along the streamlines recovered with diffusion tractography. Although this method has been widely used in the literature, in cases of voxels containing multiple fiber populations (each with different levels of myelination), tractometry provides biased results because the same value will be attributed to any bundle passing through the voxel. To overcome this bias, we propose a new method - named "myelin streamline decomposition" (MySD) - which extends convex optimization modeling for microstructure informed tractography (COMMIT) allowing the actual value measured by a microstructural map to be deconvolved on each individual streamline, thereby recovering unique bundle-specific myelin fractions (BMFs). We demonstrate the advantage of our method with respect to tractometry in well-studied bundles and compare the cortical projection of the obtained bundle-wise myelin values of both methods. We also prove the stability of our approach across different subjects and different MRI sensitive myelin mapping approaches. This work provides a proof-of-concept of in vivo investigations of entire neuronal pathways that, to date, are not possible.

Signatures of life course socioeconomic conditions in brain anatomy

Socioeconomic status (SES) plays a significant role in health and disease. At the same time, early-life conditions affect neural function and structure, suggesting the brain may be a conduit for the biological embedding of SES. Here, we investigate the brain anatomy signatures of SES in a large-scale population cohort aged 45-85 years. We assess both gray matter morphometry and tissue properties indicative of myelin content. Higher life course SES is associated with increased volume in several brain regions, including postcentral and temporal gyri, cuneus, and cerebellum. We observe more widespread volume differences and higher myelin content in the sensorimotor network but lower myelin content in the temporal lobe associated with childhood SES. Crucially, childhood SES differences persisted in adult brains even after controlling for adult SES, highlighting the unique contribution of early-life conditions to brain anatomy, independent of later changes in SES. These findings inform on the biological underpinnings of social inequality, particularly as they pertain to early-life conditions.

Macromolecular Proton Fraction as a Myelin Biomarker: Principles, Validation, and Applications

Macromolecular proton fraction (MPF) is a quantitative MRI parameter describing the magnetization transfer (MT) effect and defined as a relative amount of protons bound to biological macromolecules with restricted molecular motion, which participate in magnetic cross-relaxation with water protons. MPF attracted significant interest during past decade as a biomarker of myelin. The purpose of this mini review is to provide a brief but comprehensive summary of MPF mapping methods, histological validation studies, and MPF applications in neuroscience. Technically, MPF maps can be obtained using a variety of quantitative MT methods. Some of them enable clinically reasonable scan time and resolution. Recent studies demonstrated the feasibility of MPF mapping using standard clinical MRI pulse sequences, thus substantially enhancing the method availability. A number of studies in animal models demonstrated strong correlations between MPF and histological markers of myelin with a minor influence of potential confounders. Histological studies validated the capability of MPF to monitor both demyelination and re-myelination. Clinical applications of MPF have been mainly focused on multiple sclerosis where this method provided new insights into both white and gray matter pathology. Besides, several studies used MPF to investigate myelin role in other neurological and psychiatric conditions. Another promising area of MPF applications is the brain development studies. MPF demonstrated the capabilities to quantitatively characterize the earliest stage of myelination during prenatal brain maturation and protracted myelin development in adolescence. In summary, MPF mapping provides a technically mature and comprehensively validated myelin imaging technology for various preclinical and clinical neuroscience applications.

Rapid simultaneous acquisition of macromolecular tissue volume, susceptibility, and relaxometry maps

PURPOSE

A major obstacle to the clinical implementation of quantitative MR is the lengthy acquisition time required to derive multi-contrast parametric maps. We sought to reduce the acquisition time for QSM and macromolecular tissue volume by acquiring both contrasts simultaneously by leveraging their redundancies. The joint virtual coil concept with GRAPPA (JVC-GRAPPA) was applied to reduce acquisition time further.

METHODS

Three adult volunteers were imaged on a 3 Tesla scanner using a multi-echo 3D GRE sequence acquired at 3 head orientations. Macromolecular tissue volume, QSM, R2∗ , T₁ , and proton density maps were reconstructed. The same sequence (GRAPPA R = 4) was performed in subject 1 with a single head orientation for comparison. Fully sampled data was acquired in subject 2, from which retrospective undersampling was performed (R = 6 GRAPPA and R = 9 JVC-GRAPPA). Prospective undersampling was performed in subject 3 (R = 6 GRAPPA and R = 9 JVC-GRAPPA) using gradient blips to shift k-space sampling in later echoes.

RESULTS

Subject 1's multi-orientation and single-orientation macromolecular tissue volume maps were not significantly different based on RMSE. For subject 2, the retrospectively undersampled JVC-GRAPPA and GRAPPA generated similar results as fully sampled data. This approach was validated with the prospectively undersampled images in subject 3. Using QSM, R2∗ , and macromolecular tissue volume, the contributions of myelin and iron content to susceptibility were estimated.

CONCLUSION

We have developed a novel strategy to simultaneously acquire data for the reconstruction of 5 intrinsically coregistered 1-mm isotropic resolution multi-parametric maps, with a scan time of 6 min using JVC-GRAPPA.

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.

Direct imaging of white matter ultrashort T2∗ components at 7 Tesla

PURPOSE

To demonstrate direct imaging of the white matter ultrashort T₂^∗ components at 7 Tesla using inversion recovery (IR)-enhanced ultrashort echo time (UTE) MRI. To investigate its characteristics, potentials and limitations, and to establish a clinical protocol.

MATERIAL AND METHODS

The IR UTE technique suppresses long T₂^∗ signals within white matter by using adiabatic inversion in combination with dual-echo difference imaging. Artifacts arising at 7 T from long T₂^∗ scalp fat components were reduced by frequency shifting the IR pulse such that those frequencies were inverted likewise. For 8 healthy volunteers, the T₂^∗ relaxation times of white matter were then quantified. In 20 healthy volunteers, the UTE difference and fraction contrast were evaluated. Finally, in 6 patients with multiple sclerosis (MS), the performance of the technique was assessed.

RESULTS

A frequency shift of -1.2 ppm of the IR pulse (i.e. towards the fat frequency) provided a good suppression of artifacts. With this, an ultrashort compartment of (68 ± 6) % with a T₂^∗ time of (147 ± 58) μs was quantified with a chemical shift of (-3.6 ± 0.5) ppm from water. Within healthy volunteers' white matter, a stable ultrashort T₂^∗ fraction contrast was calculated. For the MS patients, a significant fraction reduction in the identified lesions as well as in the normal-appearing white matter was observed.

CONCLUSIONS

The quantification results indicate that the observed ultrashort components arise primarily from myelin tissue. Direct IR UTE imaging of the white matter ultrashort T₂^∗ components is thus feasible at 7 T with high quantitative inter-subject repeatability and good detection of signal loss in MS.

Quantitative Macromolecular Proton Fraction Mapping Reveals Altered Cortical Myelin Profile in Schizophrenia Spectrum Disorders

Myelin abnormalities have been reported in schizophrenia spectrum disorders (SSD) in white matter. However, in vivo examinations of cortical myeloarchitecture in SSD, especially those using quantitative measures, are limited. Here, we employed macromolecular proton fraction (MPF) obtained from quantitative magnetization transfer imaging to characterize intracortical myelin organization in 30 SSD patients versus 34 healthy control (HC) participants. We constructed cortical myelin profiles by extracting MPF values at various cortical depths and quantified their shape using a nonlinearity index (NLI). To delineate the association of illness duration with myelin changes, SSD patients were further divided into 3 duration groups. Between-group comparisons revealed reduced NLI in the SSD group with the longest illness duration (>5.5 years) compared with HC predominantly in bilateral prefrontal areas. Within the SSD group, cortical NLI decreased with disease duration and was positively associated with a measure of spatial working memory capacity as well as with cortical thickness (CT). Layer-specific analyses suggested that NLI decreases in the long-duration SSD group may arise in part from significantly increased MPF values in the midcortical layers. The current study reveals cortical myelin profile changes in SSD with illness progression, which may reflect an abnormal compensatory mechanism of the disorder.

Can MRI measure myelin? Systematic review, qualitative assessment, and meta-analysis of studies validating microstructural imaging with myelin histology

Recent years have seen an increased understanding of the importance of myelination in healthy brain function and neuropsychiatric diseases. Non-invasive microstructural magnetic resonance imaging (MRI) holds the potential to expand and translate these insights to basic and clinical human research, but the sensitivity and specificity of different MR markers to myelination is a subject of debate. To consolidate current knowledge on the topic, we perform a systematic review and meta-analysis of studies that validate microstructural imaging by combining it with myelin histology. We find meta-analytic evidence for correlations between various myelin histology metrics and markers from different MRI modalities, including fractional anisotropy, radial diffusivity, macromolecular pool, magnetization transfer ratio, susceptibility and longitudinal relaxation rate, but not mean diffusivity. Meta-analytic correlation effect sizes range widely, between R2 = 0.26 and R2 = 0.82. However, formal comparisons between MRI-based myelin markers are limited by methodological variability, inconsistent reporting and potential for publication bias, thus preventing the establishment of a single most sensitive strategy to measure myelin with MRI. To facilitate further progress, we provide a detailed characterisation of the evaluated studies as an online resource. We also share a set of 12 recommendations for future studies validating putative MR-based myelin markers and deploying them in vivo in humans.

In vivo tensor-valued diffusion MRI of focal demyelination in white and deep grey matter of rodents

•

We performed in-vivo tensor-valued diffusion MRI in demyelinating rodents.

•

Lysolecithin was injected in white and deep grey matter to cause focal demyelination.

•

Focal demyelination reduced microscopic fractional anisotropy (µFA).

•

Isotropic kurtosis may be particularly sensitive to deep grey matter lesions.

Background

Multiple sclerosis (MS) is a chronic inflammatory demyelinating disease leading to damage of white matter (WM) and grey matter (GM). Magnetic resonance imaging (MRI) is the modality of choice to assess brain damage in MS, but there is an unmet need in MRI for achieving higher sensitivity and specificity to MS-related microstructural alterations in WM and GM.

Objective

To explore whether tensor-valued diffusion MRI (dMRI) can yield sensitive microstructural read-outs for focal demyelination in cerebral WM and deep GM (DGM).

Methods

Eight rats underwent L-α-Lysophosphatidylcholine (LPC) injections in the WM and striatum to introduce focal demyelination. Multimodal MRI was performed at 7 Tesla after 7 days. Tensor-valued dMRI was complemented by diffusion tensor imaging, quantitative MRI and proton magnetic resonance spectroscopy (MRS).

Results

Quantitative MRI and MRS confirmed that LPC injections caused inflammatory demyelinating lesions in WM and DGM. Tensor-valued dMRI illustrated a significant decline of microscopic fractional anisotropy (µFA) in both LPC-treated WM and DGM (P < 0.005) along with a marked increase of isotropic kurtosis (MK_I) in DGM (P < 0.0001).

Conclusion

Tensor-valued dMRI bears considerable potential for microstructural imaging in MS, suggesting a regional µFA decrease may be a sensitive indicator of MS lesions, while a regional MK_I increase may be particularly sensitive in detecting DGM lesions of MS.

A GRM7 mutation associated with developmental delay reduces mGlu7 expression and produces neurological phenotypes

The metabotropic glutamate receptor 7 (mGlu₇) is a G protein–coupled receptor that has been recently linked to neurodevelopmental disorders. This association is supported by the identification of GRM7 variants in patients with autism spectrum disorder, attention deficit hyperactivity disorder, and severe developmental delay. One GRM7 mutation previously reported in 2 patients results in a single amino acid change, I154T, within the mGlu₇ ligand-binding domain. Here, we report 2 new patients with this mutation who present with severe developmental delay and epilepsy. Functional studies of the mGlu₇-I154T mutant reveal that this substitution resulted in significant loss of mGlu₇ protein expression in HEK293A cells and in mice. We show that this occurred posttranscriptionally at the level of protein expression and trafficking. Similar to mGlu₇–global KO mice, mGlu₇-I154T animals exhibited reduced motor coordination, deficits in contextual fear learning, and seizures. This provides functional evidence that a disease-associated mutation affecting the mGlu₇ receptor was sufficient to cause neurological dysfunction in mice and further validates GRM7 as a disease-causing gene in the human population.

Understanding aqueous and non-aqueous proton T1 relaxation in brain

A full picture of longitudinal relaxation in complex heterogeneous environments like white matter brain tissue remains elusive. In tissue, successive approximations, from the solvation layer model to the two pool model, have highlighted how longitudinal magnetization evolution depends on both inter-compartmental exchange and spin-lattice relaxation. In white matter, however, these models fail to capture the behaviour of the two distinct aqueous pools, myelin water and intra/extra-cellular water. A challenge with testing more comprehensive multi-pool models lies in directly observing all pools, both aqueous and non-aqueous. In this work, we advance these efforts by integrating three main experimental and analytical elements: direct observation of the longitudinal relaxation of both the aqueous and the non-aqueous protons in white matter, a wide range of different initial conditions, and application of an analysis pipeline which includes lineshape, CPMG, and fitting of a four pool model. An eigenvector interpretation of the four pool model highlights how longitudinal relaxation in white matter depends on initial conditions. We find that a single set of model parameters is able to describe the entire range of relaxation behaviour observed in all the separable aqueous and non-aqueous pools in experiments involving six different initial conditions. Understanding of the nature and connectedness of the tissue components is crucial in the design and interpretation of many MRI measurements, especially those based on magnetization transfer and longitudinal relaxation. In particular, the dependency of relaxation behaviour on initial conditions is likely the basis for understanding method-dependent discrepancies in in vivo T₁.

Sensitivity analysis of effective transverse shear viscoelastic and diffusional properties of myelinated white matter

Motivated by the need to interpret the results from a combined use of in vivo brain Magnetic Resonance Elastography (MRE) and Diffusion Tensor Imaging (DTI), we developed a computational framework to study the sensitivity of single-frequency MRE and DTI metrics to white matter microstructure and cell-level mechanical and diffusional properties. White matter was modeled as a triphasic unidirectional composite, consisting of parallel cylindrical inclusions (axons) surrounded by sheaths (myelin), and embedded in a matrix (glial cells plus extracellular matrix). Only 2D mechanics and diffusion in the transverse plane (perpendicular to the axon direction) was considered, and homogenized (effective) properties were derived for a periodic domain containing a single axon. The numerical solutions of the MRE problem were performed with ABAQUS and by employing a sophisticated boundary-conforming grid generation scheme. Based on the linear viscoelastic response to harmonic shear excitation and steady-state diffusion in the transverse plane, a systematic sensitivity analysis of MRE metrics (effective transverse shear storage and loss moduli) and DTI metric (effective radial diffusivity) was performed for a wide range of microstructural and intrinsic (phase-based) physical properties. The microstructural properties considered were fiber volume fraction, and the myelin sheath/axon diameter ratio. The MRE and DTI metrics are very sensitive to the fiber volume fraction, and the intrinsic viscoelastic moduli of the glial phase. The MRE metrics are nonlinear functions of the fiber volume fraction, but the effective diffusion coefficient varies linearly with it. Finally, the transverse metrics of both MRE and DTI are insensitive to the axon diameter in steady state. Our results are consistent with the limited anisotropic MRE and co-registered DTI measurements, mainly in the corpus callosum, available in the literature. We conclude that isotropic MRE and DTI constitutive models are good approximations for myelinated white matter in the transverse plane. The unidirectional composite model presented here is used for the first time to model harmonic shear stress under MRE-relevant frequency on the cell level. This model can be extended to 3D in order to inform the solution of the inverse problem in MRE, establish the biological basis of MRE metrics, and integrate MRE/DTI with other modalities towards increasing the specificity of neuroimaging.

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.

Perilesional neurodegenerative injury in multiple sclerosis: Relation to focal lesions and impact on disability

BACKGROUND

Axonal injury is the primary source of irreversible neurological decline in persons with multiple sclerosis (pwMS). Identifying and quantifying myelin and axonal loss in lesional and perilesional tissue in vivo is fundamental for a better understanding of multiple sclerosis (MS) outcomes and patient impairment. Using advanced magnetic resonance imaging (MRI) methods, consisting of selective inversion recovery quantitative magnetization transfer imaging (SIR-qMT) and multi-compartment diffusion MRI with the spherical mean technique (SMT), we conducted a cross-sectional pilot study to assess myelin and axonal damage in the normal appearing white matter (NAWM) surrounding chronic black holes (cBHs) and how this pathology correlates with disability in vivo. We hypothesized that lesional axonal transection propagates tissue injury in the surrounding NAWM and that the degree of this injury is related to patient disability.

METHODS

Eighteen pwMS underwent a 3.0 Tesla conventional clinical MRI, inclusive of T1 and T2 weighted protocols, as well as SIR-qMT and SMT. Regions of interests (ROIs) were manually delineated in cBHs, NAWM neighboring cBHs (perilesional NAWM), distant ipsilateral NAWM and contra-lateral distant NAWM. SIR-qMT-derived macromolecular-to-free pool size ratio (PSR) and SMT-derived apparent axonal volume fraction (V_ax) were extracted to infer on myelin and axonal content, respectively. Group differences were assessed using mixed-effects regression models and correlation analyses were obtained by bootstrapping 95% confidence interval.

RESULTS

In comparison to perilesional NAWM, both PSR and V_ax values were reduced in cBHs (p < 0.0001) and increased in distant contra-lateral NAWM ROIs (p < 0.001 for PSR and p < 0.0001 for V_ax) but not ipsilateral NAWM (p = 0.176 for PSR and p = 0.549 for V_ax). V_ax values measured in cBHs correlated with those in perilesional NAWM (Pearson rho = 0.63, p < 0.001). No statistically relevant associations were seen between PSR/V_ax values and clinical and/or MRI metrics of the disease with the exception of cBH PSR values, which correlated with the Expanded Disability Status Scale (Pearson rho = -0.63, p = 0.03).

CONCLUSIONS

Our results show that myelin and axonal content, detected by PSR and V_ax, are reduced in perilesional NAWM, as a function of the degree of focal cBH axonal injury. This finding is indicative of an ongoing anterograde/retrograde degeneration and suggests that treatment prevention of cBH development is a key factor for preserving NAWM integrity in surrounding tissue. It also suggests that measuring changes in perilesional areas over time may be a useful measure of outcome for proof-of-concept clinical trials on neuroprotection and repair. PSR and V_ax largely failed to capture associations with clinical and MRI characteristics, likely as a result of the small sample size and cross-sectional design, however, longitudinal assessment of a larger cohort may unravel the impact of this pathology on disease progression.

An interactive meta-analysis of MRI biomarkers of myelin

Several MRI measures have been proposed as in vivo biomarkers of myelin, each with applications ranging from plasticity to pathology. Despite the availability of these myelin-sensitive modalities, specificity and sensitivity have been a matter of discussion. Debate about which MRI measure is the most suitable for quantifying myelin is still ongoing. In this study, we performed a systematic review of published quantitative validation studies to clarify how different these measures are when compared to the underlying histology. We analyzed the results from 43 studies applying meta-analysis tools, controlling for study sample size and using interactive visualization (https://neurolibre.github.io/myelin-meta-analysis). We report the overall estimates and the prediction intervals for the coefficient of determination and find that MT and relaxometry-based measures exhibit the highest correlations with myelin content. We also show which measures are, and which measures are not statistically different regarding their relationship with histology.

Myelin water fraction estimation using small-tip fast recovery MRI

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

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

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

Traumatic and nontraumatic spinal cord injury: pathological insights from neuroimaging

Pathophysiological changes in the spinal cord white and grey matter resulting from injury can be observed with MRI techniques. These techniques provide sensitive markers of macrostructural and microstructural tissue integrity, which correlate with histological findings. Spinal cord MRI findings in traumatic spinal cord injury (tSCI) and nontraumatic spinal cord injury - the most common form of which is degenerative cervical myelopathy (DCM) - have provided important insights into the pathophysiological processes taking place not just at the focal injury site but also rostral and caudal to the spinal injury. Although tSCI and DCM have different aetiologies, they show similar degrees of spinal cord pathology remote from the injury site, suggesting the involvement of similar secondary degenerative mechanisms. Advanced quantitative MRI protocols that are sensitive to spinal cord pathology have the potential to improve diagnosis and, more importantly, predict outcomes in patients with tSCI or nontraumatic spinal cord injury. This Review describes the insights into tSCI and DCM that have been revealed by neuroimaging and outlines current activities and future directions for the field.

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.

White matter mean diffusivity correlates with myelination in tuberous sclerosis complex

Objective

Diffusion tensor imaging (DTI) of the white matter is a biomarker for neurological disease burden in tuberous sclerosis complex (TSC). To clarify the basis of abnormal diffusion in TSC, we correlated ex vivo high‐resolution diffusion imaging with histopathology in four tissue types: cortex, tuber, perituber, and white matter.

Methods

Surgical specimens of three children with TSC were scanned in a 3T or 7T MRI with a structural image isotropic resolution of 137–300 micron, and diffusion image isotropic resolution of 270‐1,000 micron. We stained for myelin (luxol fast blue, LFB), gliosis (glial fibrillary acidic protein, GFAP), and neurons (NeuN) and registered the digitized histopathology slides (0.686 micron resolution) to MRI for visual comparison. We then performed colocalization analysis in four tissue types in each specimen. Finally, we applied a linear mixed model (LMM) for pooled analysis across the three specimens.

Results

In white matter and perituber regions, LFB optical density measures correlated with fractional anisotropy (FA) and inversely with mean diffusivity (MD). In white matter only, GFAP correlated with MD, and inversely with FA. In tubers and in the cortex, there was little variation in mean LFB and GFAP signal intensity, and no correlation with MRI metrics. Neuronal density correlated with MD. In the analysis of the combined specimens, the most robust correlation was between white matter MD and LFB metrics.

Interpretation

In TSC, diffusion imaging abnormalities in microscopic tissue types correspond to specific histopathological markers. Across all specimens, white matter diffusivity correlates with myelination.

Evaluation of principal component analysis image denoising on multi-exponential MRI relaxometry

PURPOSE

Multi-exponential relaxometry is a powerful tool for characterizing tissue, but generally requires high image signal-to-noise ratio (SNR). This work evaluates the use of principal-component-analysis (PCA) denoising to mitigate these SNR demands and improve the precision of relaxometry measures.

METHODS

PCA denoising was evaluated using both simulated and experimental MRI data. Bi-exponential transverse relaxation signals were simulated for a wide range of acquisition and sample parameters, and experimental data were acquired from three excised and fixed mouse brains. In both cases, standard relaxometry analysis was performed on both original and denoised image data, and resulting estimated signal parameters were compared.

RESULTS

Denoising reduced the root-mean-square-error of parameters estimated from multi-exponential relaxometry by factors of ≈3×, for typical acquisition and sample parameters. Denoised images and subsequent parameter maps showed little or no signs of spatial artifact or loss of resolution.

CONCLUSION

Experimental studies and simulations demonstrate that PCA denoising of MRI relaxometry data is an effective method of improving parameter precision without sacrificing image resolution. This simple yet important processing step thus paves the way for broader applicability of multi-exponential MRI relaxometry.

Formalin tissue fixation biases myelin-sensitive MRI

PURPOSE

Chemical fixatives such as formalin form cross-links between proteins and affect the relaxation times and diffusion properties of tissue. These fixation-induced changes likely also affect myelin density measurements produced by quantitative magnetization transfer and myelin water imaging. In this work, we evaluate these myelin-sensitive MRI methods for fixation-induced biases.

METHODS

We perform quantitative magnetization transfer, myelin water imaging, and deuterium oxide-exchanged zero TE imaging on unfixed human spinal cord tissue at 9.4 Tesla and repeat these measurements after 1 day and 31 days of formalin fixation.

RESULTS

The quantitative magnetization-transfer bound pool fraction increased by 30.7% ± 21.1% after 1 day of fixation and by 42.6% ± 33.9% after 31 days of fixation. Myelin water fraction increased by 39.7% ± 15.5% and 37.0% ± 15.9% at these same time points, and mean T₂ of the myelin water pool nearly doubled. Reference-normalized deuterium oxide-exchanged zero TE signal intensity increased by 8.17% ± 6.03% after 31 days of fixation but did not change significantly after 1 day of fixation. After fixation, specimen cross-sectional area decreased by approximately 5%; after correction for shrinkage, changes in deuterium oxide-exchanged zero TE intensity were nearly eliminated.

CONCLUSION

Bound pool fraction and myelin water fraction are significantly increased by formalin fixation, whereas deuterium oxide-exchanged zero TE intensity is minimally affected. Changes in quantitative magnetization transfer and myelin water imaging may be due in part to delamination and formation of vacuoles in the myelin sheath. Deuterium oxide-exchanged signal intensity may be altered by fixation-induced changes in myelin lipid solid-state ¹ H T₁ . We urge caution in the comparison of these measurements across subjects or specimens in different states, especially unfixed versus fixed tissue.

Iron, Myelin, and the Brain: Neuroimaging Meets Neurobiology

Although iron is crucial for neuronal functioning, many aspects of cerebral iron biology await clarification. The ability to quantify specific iron forms in the living brain would open new avenues for diagnosis, therapeutic monitoring, and understanding pathogenesis of diseases. A modality that allows assessment of brain tissue composition in vivo, in particular of iron deposits or myelin content on a submillimeter spatial scale, is magnetic resonance imaging (MRI). Multimodal strategies combining MRI with complementary analytical techniques ex vivo have emerged, which may lead to improved specificity. Interdisciplinary collaborations will be key to advance beyond simple correlative analyses in the biological interpretation of MRI data and to gain deeper insights into key factors leading to iron accumulation and/or redistribution associated with neurodegeneration.

Modeling conduction delays in the corpus callosum using MRI-measured g-ratio

Conduction of action potentials along myelinated axons is affected by their structural features, such as the axonal g-ratio, the ratio between the inner and outer diameters of the myelin sheath surrounding the axon. The effect of g-ratio variance on conduction properties has been quantitatively evaluated using single-axon models. It has recently become possible to estimate a g-ratio weighted measurement in vivo using quantitative MRI. Nevertheless, it is still unclear whether the variance in the g-ratio in the healthy human brain leads to significant differences in conduction velocity. In this work we tested whether the g-ratio MRI measurement can be used to predict conduction delays in the corpus callosum. We present a comprehensive framework in which the structural properties of fibers (i.e. length and g-ratio, measured using MRI), are incorporated in a biophysical model of axon conduction, to model conduction delays of long-range white matter fibers. We applied this framework to the corpus callosum, and found conduction delay estimates that are compatible with previously estimated values of conduction delays. We account for the variance in the velocity given the axon diameter distribution in the splenium, mid-body and genu, to further compare the fibers within the corpus callosum. Conduction delays have been suggested to increase with age. Therefore, we investigated whether there are differences in the g-ratio and the fiber length between young and old adults, and whether this leads to a difference in conduction speed and delays. We found very small differences between the predicted delays of the two groups in the motor fibers of the corpus callosum. We also found that the motor fibers of the corpus callosum have the fastest conduction estimates. Using the axon diameter distributions, we found that the occipital fibers have the slowest estimations, while the frontal and motor fiber tracts have similar estimates. Our study provides a framework for predicting conduction latencies in vivo. The framework could have major implications for future studies of white matter diseases and large range network computations. Our results highlight the need for improving additional in vivo measurements of white matter microstructure.

Arterial stiffness and white matter integrity in the elderly: A diffusion tensor and magnetization transfer imaging study

BACKGROUND AND PURPOSE

The stiffness of large arteries and increased pulsatility can have an impact on the brain white matter (WM) microstructure, however those mechanisms are still poorly understood. The aim of this study was to investigate the association between central artery stiffness, axonal and myelin integrity in 54 cognitively unimpaired elderly subjects (65-75 years old).

METHODS

The neuronal fiber integrity of brain WM was assessed using diffusion tensor metrics and magnetization transfer imaging as measures of axonal organization (Fractional anisotropy, Radial diffusivity) and state of myelination (Myelin volume fraction). Central artery stiffness was measured by carotid-femoral pulse wave velocity (cfPWV). Statistical analyses included 4 regions (the corpus callosum, the internal capsule, the corona radiata and the superior longitudinal fasciculus) which have been previously denoted as vulnerable to increased central artery stiffness.

RESULTS

cfPWV was significantly associated with fractional anisotropy and radial diffusivity (p < 0.05, corrected for multiple comparisons) but not with myelin volume fraction. Findings from this study also show that improved executive function performance correlates with Fractional anisotropy positively (p < 0.05 corrected) as well as with myelin volume fraction and radial diffusivity negatively (p < 0.05 corrected).

CONCLUSIONS

These findings suggest that arterial stiffness is associated with axon degeneration rather than demyelination. Controlling arterial stiffness may play a role in maintaining the health of WM axons in the aging brain.

On modeling

Mapping tissue microstructure with MRI holds great promise as a noninvasive window into tissue organization at the cellular level. Having originated within the realm of diffusion NMR in the late 1970s, this field is experiencing an exponential growth in the number of publications. At the same time, model-based approaches are also increasingly incorporated into advanced MRI acquisition and reconstruction techniques. However, after about two decades of intellectual and financial investment, microstructural mapping has yet to find a single commonly accepted clinical application. Here, we suggest that slow progress in clinical translation may signify unresolved fundamental problems. We outline such problems and related practical pitfalls, as well as review strategies for developing and validating tissue microstructure models, to provoke a discussion on how to bridge the gap between our scientific aspirations and the clinical reality. We argue for recalibrating the efforts of our community toward a more systematic focus on fundamental research aimed at identifying relevant degrees of freedom affecting the measured MR signal. Such a focus is essential for realizing the truly revolutionary potential of noninvasive three-dimensional in vivo microstructural mapping.

SyMRI of the Brain: Rapid Quantification of Relaxation Rates and Proton Density, With Synthetic MRI, Automatic Brain Segmentation, and Myelin Measurement

Conventional magnetic resonance images are usually evaluated using the image signal contrast between tissues and not based on their absolute signal intensities. Quantification of tissue parameters, such as relaxation rates and proton density, would provide an absolute scale; however, these methods have mainly been performed in a research setting. The development of rapid quantification, with scan times in the order of 6 minutes for full head coverage, has provided the prerequisites for clinical use. The aim of this review article was to introduce a specific quantification method and synthesis of contrast-weighted images based on the acquired absolute values, and to present automatic segmentation of brain tissues and measurement of myelin based on the quantitative values, along with application of these techniques to various brain diseases. The entire technique is referred to as “SyMRI” in this review. SyMRI has shown promising results in previous studies when used for multiple sclerosis, brain metastases, Sturge-Weber syndrome, idiopathic normal pressure hydrocephalus, meningitis, and postmortem imaging.

Microstructural imaging in the spinal cord and validation strategies

In vivo histology using magnetic resonance imaging (MRI) is a newly emerging research field that aims to non-invasively characterize tissue microstructure. The implications of in vivo histology are many, from discovering novel biomarkers to studying human development, to providing tools for disease diagnosis and monitoring the effects of novel treatments on tissue. This review focuses on quantitative MRI (qMRI) techniques that are used to map spinal cord microstructure. Opening with a rationale for non-invasive imaging of the spinal cord, this article continues with a brief overview of the existing MRI techniques for axon and myelin imaging, followed by the specific challenges and potential solutions for acquiring and processing such data. The final part of this review focuses on histological validation, with suggested tissue preparation, acquisition and processing protocols for large-scale microscopy.

Inferring brain tissue composition and microstructure via MR relaxometry

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

Experimental studies of g-ratio MRI in ex vivo mouse brain

This study aimed to experimentally evaluate a previously proposed MRI method for mapping axonal g-ratio (ratio of axon diameters, measured to the inner and outer boundary of myelin). MRI and electron microscopy were used to study excised and fixed brains of control mice and three mouse models of abnormal white matter. The results showed that g-ratio measured with MRI correlated with histological measures of myelinated axon g-ratio, but with a bias that is likely due to the presence of non-myelinated axons. The results also pointed to cases where the MRI g-ratio model simplifies to be primarily a function of total myelin content.

Abnormal social behavior in mice with tyrosinemia type I is associated with an increase of myelin in the cerebral cortex

Hereditary tyrosinemia type I (HT1) is caused by mutations in the fumarylacetoacetate hydrolase (FAH) gene, the template for the final enzyme in the tyrosine catabolism pathway. If left untreated this deficiency of functional FAH leads to a buildup of toxic metabolites that can cause liver disease, kidney dysfunction and high mortality. The current treatment with the drug NTBC prevents the production of these metabolites and has consequently increased the survival rate in HT1 children. As a result of this increased survival, long term complications of HT1 are now being observed, including slower learning, impaired cognition and altered social behavior. We studied a mouse model of HT1 to gain insight into the effects of HT1 and treatment with NTBC on social behavior in mice. We showed that mice with HT1 display abnormal social behavior in that they spend more time in the absence of another mouse and do not discriminate between a novel mouse and an already familiar mouse. This altered behavior was due to HT1 and not treatment with NTBC. Quantification of cerebral cortex myelin in mice with HT1 showed a two to threefold increase in myelin expression. Our findings suggest that absence of FAH expression in the brain produces an altered brain biochemistry resulting in increased expression of myelin. This increase in myelination could lead to abnormal action potential velocity and altered neuronal connections that provide a mechanism for the altered learning, social behavior and cognitive issues recently seen in HT1.

Loss of mTORC2 signaling in oligodendrocyte precursor cells delays myelination

Myelin abnormalities are increasingly being recognized as an important component of a number of neurologic developmental disorders. The integration of many signaling pathways and cell types are critical for correct myelinogenesis. The PI3-K and mechanistic target of rapamycin (mTOR) pathways have been found to play key roles. mTOR is found within two distinct complexes, mTORC1 and mTORC2. mTORC1 activity has been shown to play a major role during myelination, while the role of mTORC2 is not yet well understood. To determine the role of mTORC2 signaling in myelinogenesis, we generated a mouse lacking the critical mTORC2 component Rictor in oligodendrocyte precursors (OPCs). Targeted deletion of Rictor in these cells decreases and delays the expression of myelin related proteins and reduces the size of cerebral white matter tracts. This is developmentally manifest as a transient reduction in myelinated axon density and g-ratio. OPC cell number is reduced at birth without detectable change in proliferation with proportional reductions in mature oligodendrocyte number at P15. The total number of oligodendrocytes as well as extent of myelination, does improve over time. Adult conditional knock-out (CKO) animals do not demonstrate a behavioral phenotype likely due in part to preserved axonal conduction velocities. These data support and extend prior studies demonstrating an important but transient contribution of mTORC2 signaling to myelin development.

Analysis of White Matter Damage in Patients with Multiple Sclerosis via a Novel In Vivo MR Method for Measuring Myelin, Axons, and G-Ratio

BACKGROUND AND PURPOSE

Myelin and axon volume fractions can now be estimated via MR imaging in vivo, as can the g-ratio, which equals the ratio of the inner to the outer diameter of a nerve fiber. The purpose of this study was to evaluate WM damage in patients with MS via this novel MR imaging technique.

MATERIALS AND METHODS

Twenty patients with relapsing-remitting MS with a combined total of 149 chronic plaques were analyzed. Myelin volume fraction was calculated based on simultaneous tissue relaxometry. Intracellular and CSF compartment volume fractions were quantified via neurite orientation dispersion and density imaging. Axon volume fraction and g-ratio were calculated by combining these measurements. Myelin and axon volume fractions and g-ratio were measured in plaques, periplaque WM, and normal-appearing WM.

RESULTS

All metrics differed significantly across the 3 groups (P < .001, except P = .027 for g-ratio between periplaque WM and normal-appearing WM). Those in plaques differed most from those in normal-appearing WM. The percentage changes in plaque and periplaque WM metrics relative to normal-appearing WM were significantly larger in absolute value for myelin volume fraction than for axon volume fraction and g-ratio (P < .001, except P = .033 in periplaque WM relative to normal-appearing WM for comparison between myelin and axon volume fraction).

CONCLUSIONS

In this in vivo MR imaging study, the myelin of WM was more damaged than axons in plaques and periplaque WM of patients with MS. Myelin and axon volume fractions and g-ratio may potentially be useful for evaluating WM damage in patients with MS.

Promise and pitfalls of g-ratio estimation with MRI

The fiber g-ratio is the ratio of the inner to the outer diameter of the myelin sheath of a myelinated axon. It has a limited dynamic range in healthy white matter, as it is optimized for speed of signal conduction, cellular energetics, and spatial constraints. In vivo imaging of the g-ratio in health and disease would greatly increase our knowledge of the nervous system and our ability to diagnose, monitor, and treat disease. MRI based g-ratio imaging was first conceived in 2011, and expanded to be feasible in full brain white matter with preliminary results in 2013. This manuscript reviews the growing g-ratio imaging literature and speculates on future applications. It details the methodology for imaging the g-ratio with MRI, and describes the known pitfalls and challenges in doing so.

Evaluating g-ratio weighted changes in the corpus callosum as a function of age and sex

Recent years have seen a growing interest in relating MRI measurements to the structural-biophysical properties of white matter fibers. The fiber g-ratio, defined as the ratio between the inner and outer radii of the axon myelin sheath, is an important structural property of white matter, affecting signal conduction. Recently proposed modeling methods that use a combination of quantitative-MRI signals, enable a measurement of the fiber g-ratio in vivo. Here we use an MRI-based g-ratio estimation to observe the variance of the g-ratio within the corpus callosum, and evaluate sex and age related differences. To estimate the g-ratio we used a model (Stikov et al., 2011; Duval et al., 2017) based on two different WM microstructure parameters: the relative amounts of myelin (myelin volume fraction, MVF) and fibers (fiber volume fraction, FVF) in a voxel. We derived the FVF from the fractional anisotropy (FA), and estimated the MVF by using the lipid and macromolecular tissue volume (MTV), calculated from the proton density (Mezer et al., 2013). In comparison to other methods of estimating the MVF, MTV represents a stable parameter with a straightforward route of acquisition. To establish our model, we first compared histological MVF measurements (West et al., 2016) with the MRI derived MTV. We then implemented our model on a large database of 92 subjects (44 males), aged 7 to 81, in order to evaluate age and sex related changes within the corpus callosum. Our results show that the MTV provides a good estimation of MVF for calculating g-ratio, and produced values from the corpus callosum that correspond to those found in animals ex vivo and are close to the theoretical optimum, as well as to published in vivo data. Our results demonstrate that the MTV derived g-ratio provides a simple and reliable in vivo g-ratio-weighted (GR*) measurement in humans. In agreement with theoretical predictions, and unlike other tissue parameters measured with MRI, the g-ratio estimations were found to be relatively stable with age, and we found no support for a significant sexual dimorphism with age.