A revised model for estimating g-ratio from MRI

Modeling and simulation for prediction of multiple sclerosis progression

In light of extensive work that has created a wide range of techniques for predicting the course of multiple sclerosis (MS) disease, this paper attempts to provide an overview of these approaches and put forth an alternative way to predict the disease progression. For this purpose, the existing methods for estimating and predicting the course of the disease have been categorized into clinical, radiological, biological, and computational or artificial intelligence-based markers. Weighing the weaknesses and strengths of these prognostic groups is a profound method that is yet in need and works directly at the level of diseased connectivity. Therefore, we propose using the computational models in combination with established connectomes as a predictive tool for MS disease trajectories. The fundamental conduction-based Hodgkin-Huxley model emerged as promising from examining these studies. The advantage of the Hodgkin-Huxley model is that certain properties of connectomes, such as neuronal connection weights, spatial distances, and adjustments of signal transmission rates, can be taken into account. It is precisely these properties that are particularly altered in MS and that have strong implications for processing, transmission, and interactions of neuronal signaling patterns. The Hodgkin-Huxley (HH) equations as a point-neuron model are used for signal propagation inside a small network. The objective is to change the conduction parameter of the neuron model, replicate the changes in myelin properties in MS and observe the dynamics of the signal propagation across the network. The model is initially validated for different lengths, conduction values, and connection weights through three nodal connections. Later, these individual factors are incorporated into a small network and simulated to mimic the condition of MS. The signal propagation pattern is observed after inducing changes in conduction parameters at certain nodes in the network and compared against a control model pattern obtained before the changes are applied to the network. The signal propagation pattern varies as expected by adapting to the input conditions. Similarly, when the model is applied to a connectome, the pattern changes could give an insight into disease progression. This approach has opened up a new path to explore the progression of the disease in MS. The work is in its preliminary state, but with a future vision to apply this method in a connectome, providing a better clinical tool.

The effects of axonal beading and undulation on axonal diameter estimation from diffusion MRI: Insights from simulations in human axons segmented from three-dimensional electron microscopy

The increasing availability of high-performance gradient systems in human MRI scanners has generated great interest in diffusion microstructural imaging applications such as axonal diameter mapping. Practically, sensitivity to axon diameter in diffusion MRI is attained at strong diffusion weightings b , where the deviation from the expected 1 / b scaling in white matter yields a finite transverse diffusivity, which is then translated into an axon diameter estimate. While axons are usually modeled as perfectly straight, impermeable cylinders, local variations in diameter (caliber variation or beading) and direction (undulation) are known to influence axonal diameter estimates and have been observed in microscopy data of human axons. In this study, we performed Monte Carlo simulations of diffusion in axons reconstructed from three-dimensional electron microscopy of a human temporal lobe specimen using simulated sequence parameters matched to the maximal gradient strength of the next-generation Connectome 2.0 human MRI scanner ( ≲ 500 mT/m). We show that axon diameter estimation is accurate for nonbeaded, nonundulating fibers; however, in fibers with caliber variations and undulations, the axon diameter is heavily underestimated due to caliber variations, and this effect overshadows the known overestimation of the axon diameter due to undulations. This unexpected underestimation may originate from variations in the coarse-grained axial diffusivity due to caliber variations. Given that increased axonal beading and undulations have been observed in pathological tissues, such as traumatic brain injury and ischemia, the interpretation of axon diameter alterations in pathology may be significantly confounded.

Pathological ultrastructural alterations of myelinated axons in normal appearing white matter in progressive multiple sclerosis

Multiple sclerosis (MS) pathophysiology includes inflammation, demyelination and neurodegeneration, but the exact mechanisms of disease initiation and progression are unknown. A major feature of lesions is lack of myelin, which increases axonal energy demand and requires adaptation in number and size of mitochondria. Outside lesions, subtle and diffuse alterations are observed in normal appearing white matter (NAWM) and normal appearing grey matter (NAGM), including increased oxidative stress, reduced axon density and changes in myelin composition and morphology. On an ultrastructural level, only limited data is available on alterations in myelinated axons. We generated large scale 2D scanning transmission electron microscopy images ('nanotomy') of non-demyelinated brain tissue of control and progressive MS donors, accessible via an open-access online repository. We observed a reduced density of myelinated axons in NAWM, without a decrease in cross-sectional axon area. Small myelinated axons were less frequently and large myelinated axons were more frequently present in NAWM, while the g-ratio was similar. The correlation between axonal mitochondrial radius and g-ratio was lost in NAWM, but not in NAGM. Myelinated axons in control GM and NAGM had a similar g-ratio and radius distribution. We hypothesize that axonal loss in NAWM is likely compensated by swelling of the remaining myelinated axons and subsequent adjustment of myelin thickness to maintain their g-ratio. Failure of axonal mitochondria to adjust their size and fine-tuning of myelin thickness may render NAWM axons and their myelin more susceptible to injury.

Investigation of the association between central arterial stiffness and aggregate g-ratio in cognitively unimpaired adults

Stiffness of the large arteries has been shown to impact cerebral white matter (WM) microstructure in both younger and older adults. However, no study has yet demonstrated an association between arterial stiffness and aggregate g-ratio, a specific magnetic resonance imaging (MRI) measure of axonal myelination that is highly correlated with neuronal signal conduction speed. In a cohort of 38 well-documented cognitively unimpaired adults spanning a wide age range, we investigated the association between central arterial stiffness, measured using pulse wave velocity (PWV), and aggregate g-ratio, measured using our recent advanced quantitative MRI methodology, in several cerebral WM structures. After adjusting for age, sex, smoking status, and systolic blood pressure, our results indicate that higher PWV values, that is, elevated arterial stiffness, were associated with lower aggregate g-ratio values, that is, lower microstructural integrity of WM. Compared to other brain regions, these associations were stronger and highly significant in the splenium of the corpus callosum and the internal capsules, which have been consistently documented as very sensitive to elevated arterial stiffness. Moreover, our detailed analysis indicates that these associations were mainly driven by differences in myelination, measured using myelin volume fraction, rather than axonal density, measured using axonal volume fraction. Our findings suggest that arterial stiffness is associated with myelin degeneration, and encourages further longitudinal studies in larger study cohorts. Controlling arterial stiffness may represent a therapeutic target in maintaining the health of WM tissue in cerebral normative aging.

The influence of axonal beading and undulation on axonal diameter mapping

We consider the effect of non-cylindrical axonal shape on axonal diameter mapping with diffusion MRI. Practical sensitivity to axon diameter is attained at strong diffusion weightings b , where the deviation from the 1 / b scaling yields the finite transverse diffusivity, which is then translated into axon diameter. While axons are usually modeled as perfectly straight, impermeable cylinders, the local variations in diameter (caliber variation or beading) and direction (undulation) have been observed in microscopy data of human axons. Here we quantify the influence of cellular-level features such as caliber variation and undulation on axon diameter estimation. For that, we simulate the diffusion MRI signal in realistic axons segmented from 3-dimensional electron microscopy of a human brain sample. We then create artificial fibers with the same features and tune the amplitude of their caliber variations and undulations. Numerical simulations of diffusion in fibers with such tunable features show that caliber variations and undulations result in under- and over-estimation of axon diameters, correspondingly; this bias can be as large as 100%. Given that increased axonal beading and undulations have been observed in pathological tissues, such as traumatic brain injury and ischemia, the interpretation of axon diameter alterations in pathology may be significantly confounded.

Temporal trajectories of normal myelination and axonal development assessed by quantitative macromolecular and diffusion MRI: Ultrastructural and immunochemical validation in a rabbit model

INTRODUCTION

Quantitative and non-invasive measures of brain myelination and maturation during development are of great importance to both clinical and translational research communities. While the metrics derived from diffusion tensor imaging, are sensitive to developmental changes and some pathologies, they remain difficult to relate to the actual microstructure of the brain tissue. The advent of advanced model-based microstructural metrics requires histological validation. The purpose of the study was to validate novel, model-based MRI techniques, such as macromolecular proton fraction mapping (MPF) and neurite orientation and dispersion indexing (NODDI), against histologically derived indexes of myelination and microstructural maturation at various stages of development.

METHODS

New Zealand White rabbit kits underwent serial in-vivo MRI examination at postnatal days 1, 5, 11, 18, and 25, and as adults. Multi-shell, diffusion-weighted experiments were processed to fit NODDI model to obtain estimates, intracellular volume fraction (ICVF) and orientation dispersion index (ODI). Macromolecular proton fraction (MPF) maps were obtained from three source (MT-, PD-, and T1-weighted) images. After MRI sessions, a subset of animals was euthanized and regional samples of gray and white matter were taken for western blot analysis, to determine myelin basic protein (MBP), and electron microscopy, to estimate axonal, myelin fractions and g-ratio.

RESULTS

MPF of white matter regions showed a period of fast growth between P5 and P11 in the internal capsule, with a later onset in the corpus callosum. This MPF trajectory was in agreement with levels of myelination in the corresponding brain region, as assessed by western blot and electron microscopy. In the cortex, the greatest increase of MPF occurred between P18 and P26. In contrast, myelin, according to MBP western blot, saw the largest hike between P5 and P11 in the sensorimotor cortex and between P11 and P18 in the frontal cortex, which then seemingly plateaued after P11 and P18 respectively. G-ratio by MRI markers decreased with age in the white matter. However, electron microscopy suggest a relatively stable g-ratio throughout development.

CONCLUSION

Developmental trajectories of MPF accurately reflected regional differences of myelination rate in different cortical regions and white matter tracts. MRI-derived estimation of g-ratio was inaccurate during early development, likely due to the overestimation of axonal volume fraction by NODDI due to the presence of a large proportion of unmyelinated axons.

Conduction velocity along a key white matter tract is associated with autobiographical memory recall ability

Conduction velocity is the speed at which electrical signals travel along axons and is a crucial determinant of neural communication. Inferences about conduction velocity can now be made in vivo in humans using a measure called the magnetic resonance (MR) g-ratio. This is the ratio of the inner axon diameter relative to that of the axon plus the myelin sheath that encases it. Here, in the first application to cognition, we found that variations in MR g-ratio, and by inference conduction velocity, of the parahippocampal cingulum bundle were associated with autobiographical memory recall ability in 217 healthy adults. This tract connects the hippocampus with a range of other brain areas. We further observed that the association seemed to be with inner axon diameter rather than myelin content. The extent to which neurites were coherently organised within the parahippocampal cingulum bundle was also linked with autobiographical memory recall ability. Moreover, these findings were specific to autobiographical memory recall and were not apparent for laboratory-based memory tests. Our results offer a new perspective on individual differences in autobiographical memory recall ability, highlighting the possible influence of specific white matter microstructure features on conduction velocity when recalling detailed memories of real-life past experiences.

Myelin characteristics of the corpus callosum in capuchin monkeys (Sapajus [Cebus] apella) across the lifespan

The midsagittal area of the corpus callosum (CC) is frequently studied in relation to brain development, connectivity, and function. Here we quantify myelin characteristics from electron microscopy to understand more fully differential patterns of white matter development occurring within the CC. We subdivided midsagittal regions of the CC into: I-rostrum and genu, II-rostral body, III-anterior midbody, IV-posterior midbody, and V-isthmus and splenium. The sample represented capuchin monkeys ranging in age from 2 weeks to 35 years (Sapajus [Cebus] apella, n = 8). Measurements of myelin thickness, myelin fraction, and g-ratio were obtained in a systematic random fashion. We hypothesized there would be a period of rapid myelin growth within the CC in early development. Using a locally weighted regression analysis (LOESS), we found regional differences in myelin characteristics, with posterior regions showing more rapid increases in myelin thickness and sharper decreases in g-ratio in early development. The most anterior region showed the most sustained growth in myelin thickness. For all regions over the lifespan, myelin fraction increased, plateaued, and decreased. These results suggest differential patterns of nonlinear myelin growth occur early in development and well into adulthood in the CC of capuchin monkeys.

In vivo Estimation of Axonal Morphology From Magnetic Resonance Imaging and Electroencephalography Data

Purpose

We present a novel approach that allows the estimation of morphological features of axonal fibers from data acquired in vivo in humans. This approach allows the assessment of white matter microscopic properties non-invasively with improved specificity.

Theory

The proposed approach is based on a biophysical model of Magnetic Resonance Imaging (MRI) data and of axonal conduction velocity estimates obtained with Electroencephalography (EEG). In a white matter tract of interest, these data depend on (1) the distribution of axonal radius [P(r)] and (2) the g-ratio of the individual axons that compose this tract [g(r)]. P(r) is assumed to follow a Gamma distribution with mode and scale parameters, M and θ, and g(r) is described by a power law with parameters α and β.

Methods

MRI and EEG data were recorded from 14 healthy volunteers. MRI data were collected with a 3T scanner. MRI-measured g-ratio maps were computed and sampled along the visual transcallosal tract. EEG data were recorded using a 128-lead system with a visual Poffenberg paradigm. The interhemispheric transfer time and axonal conduction velocity were computed from the EEG current density at the group level. Using the MRI and EEG measures and the proposed model, we estimated morphological properties of axons in the visual transcallosal tract.

Results

The estimated interhemispheric transfer time was 11.72 ± 2.87 ms, leading to an average conduction velocity across subjects of 13.22 ± 1.18 m/s. Out of the 4 free parameters of the proposed model, we estimated θ – the width of the right tail of the axonal radius distribution – and β – the scaling factor of the axonal g-ratio, a measure of fiber myelination. Across subjects, the parameter θ was 0.40 ± 0.07 μm and the parameter β was 0.67 ± 0.02 μm^−α.

Conclusion

The estimates of axonal radius and myelination are consistent with histological findings, illustrating the feasibility of this approach. The proposed method allows the measurement of the distribution of axonal radius and myelination within a white matter tract, opening new avenues for the combined study of brain structure and function, and for in vivo histological studies of the human brain.

Resolution and b value dependent Structural Connectome in ex vivo Mouse Brain

Diffusion magnetic resonance imaging has been widely used in both clinical and preclinical studies to characterize tissue microstructure and structural connectivity. The diffusion MRI protocol for the Human Connectome Project (HCP) has been developed and optimized to obtain high-quality, high-resolution diffusion MRI (dMRI) datasets. However, such efforts have not been fully explored in preclinical studies, especially for rodents. In this study, high quality dMRI datasets of mouse brains were acquired at 9.4T system from two vendors. In particular, we acquired a high-spatial resolution dMRI dataset (25 μm isotropic with 126 diffusion encoding directions), which we believe to be the highest spatial resolution yet obtained; and a high-angular resolution dMRI dataset (50 μm isotropic with 384 diffusion encoding directions), which we believe to be the highest angular resolution compared to the dMRI datasets at the microscopic resolution. We systematically investigated the effects of three important parameters that affect the final outcome of the connectome: b value (1000s/mm² to 8000 s/mm²), angular resolution (10 to 126), and spatial resolution (25 μm to 200 μm). The stability of tractography and connectome increase with the angular resolution, where more than 50 angles is necessary to achieve consistent results. The connectome and quantitative parameters derived from graph theory exhibit a linear relationship to the b value (R² > 0.99); a single-shell acquisition with b value of 3000 s/mm² shows comparable results to the multi-shell high angular resolution dataset. The dice coefficient decreases and both false positive rate and false negative rate gradually increase with coarser spatial resolution. Our study provides guidelines and foundations for exploration of tradeoffs among acquisition parameters for the structural connectome in ex vivo mouse brain.

AAV1.NT-3 gene therapy for X-linked Charcot-Marie-Tooth neuropathy type 1

X-linked Charcot-Marie-Tooth neuropathy (CMTX) is caused by mutations in the gene encoding Gap Junction Protein Beta-1 (GJB1)/Connexin32 (Cx32) in Schwann cells. Neurotrophin-3 (NT-3) is an important autocrine factor supporting Schwann cell survival and differentiation and stimulating axon regeneration and myelination. Improvements in these parameters have been shown previously in a CMT1 model, TremblerJ mouse, with NT-3 gene transfer therapy. For this study, scAAV1.tMCK.NT-3 was delivered to the gastrocnemius muscle of 3-month-old Cx32 knockout (KO) mice. Measurable levels of NT-3 were found in the serum at 6-month post gene delivery. The outcome measures included functional, electrophysiological and histological assessments. At 9-months of age, NT-3 treated mice showed no functional decline with normalized compound muscle action potential amplitudes. Myelin thickness and nerve conduction velocity significantly improved compared with untreated cohort. A normalization toward age-matched wildtype histopathological parameters included increased number of Schmidt-Lanterman incisures, and muscle fiber diameter. Collectively, these findings suggest a translational application to CMTX1.

Towards a representative reference for MRI-based human axon radius assessment using light microscopy

Non-invasive assessment of axon radii via MRI bears great potential for clinical and neuroscience research as it is a main determinant of the neuronal conduction velocity. However, there is a lack of representative histological reference data at the scale of the cross-section of MRI voxels for validating the MRI-visible, effective radius (r_eff). Because the current gold standard stems from neuroanatomical studies designed to estimate the bulk-determined arithmetic mean radius (r_arith) on small ensembles of axons, it is unsuited to estimate the tail-weighted r_eff. We propose CNN-based segmentation on high-resolution, large-scale light microscopy (lsLM) data to generate a representative reference for r_eff. In a human corpus callosum, we assessed estimation accuracy and bias of r_arith and r_eff. Furthermore, we investigated whether mapping anatomy-related variation of r_arith and r_eff is confounded by low-frequency variation of the image intensity, e.g., due to staining heterogeneity. Finally, we analyzed the error due to outstandingly large axons in r_eff. Compared to r_arith, r_eff was estimated with higher accuracy (maximum normalized-root-mean-square-error of r_eff: 8.5 %; r_arith: 19.5 %) and lower bias (maximum absolute normalized-mean-bias-error of r_eff: 4.8 %; r_arith: 13.4 %). While r_arith was confounded by variation of the image intensity, variation of r_eff seemed anatomy-related. The largest axons contributed between 0.8 % and 2.9 % to r_eff. In conclusion, the proposed method is a step towards representatively estimating r_eff at MRI voxel resolution. Further investigations are required to assess generalization to other brains and brain areas with different axon radii distributions.

Cerebral aggregate g-ratio mapping using magnetic resonance relaxometry and diffusion tensor imaging to investigate sex and age-related differences in white matter microstructure

Axonal demyelination is a cardinal feature of aging and age-related diseases. The g-ratio, mathematically defined as the inner-to-outer diameter of a myelinated axon, is used as a structural index of optimal axonal myelination and has been shown to represent a sensitive imaging biomarker of microstructural integrity. Several magnetic resonance imaging (MRI) methods for whole-brain mapping of aggregate g-ratio have been introduced. Computation of the aggerate g-ratio requires estimates of the myelin volume fraction (MVF) and the axonal volume fraction (AVF). While accurate determinations of MVF and AVF can be obtained through multicomponent relaxometry or diffusion analyses, respectively, these methods require lengthy acquisition times making their implementation challenging in a clinical context. Therefore, any attempt to overcome this drawback is needed. Expanding on our previous work, we introduced a new MRI method for whole-brain mapping of aggregate g-ratio. This new approach is based on the use of a single-shell diffusion for AVF determination, reducing the acquisition time by approximately ~10 min from our recently introduced approach, while offering the possibility to investigate g-ratio differences in previous studies with existing data for MVF mapping and single-shell diffusion data for AVF mapping. Our comparison analysis indicates that our newly derived aggregate g-ratio values were similar to those derived from our previous method, which requires a longer acquisition time. Further, in agreement with our previous observations, we found quadratic U-shaped relationships between aggregate g-ratio and age in this much larger study cohort. However, our results show that sexual dimorphism in g-ratio was not significant in any brain region investigated.

Overview of myelin, major myelin lipids, and myelin-associated proteins

Myelin is a modified cell membrane that forms a multilayer sheath around the axon. It retains the main characteristics of biological membranes, such as lipid bilayer, but differs from them in several important respects. In this review, we focus on aspects of myelin composition that are peculiar to this structure and differentiate it from the more conventional cell membranes, with special attention to its constituent lipid components and several of the most common and important myelin proteins: myelin basic protein, proteolipid protein, and myelin protein zero. We also discuss the many-fold functions of myelin, which include reliable electrical insulation of axons to ensure rapid propagation of nerve impulses, provision of trophic support along the axon and organization of the unmyelinated nodes of Ranvier, as well as the relationship between myelin biology and neurologic disease such as multiple sclerosis. We conclude with a brief history of discovery in the field and outline questions for future research.

Motor Learning and Physical Exercise in Adaptive Myelination and Remyelination

The idea that myelination is driven by both intrinsic and extrinsic cues has gained much traction in recent years. Studies have demonstrated that myelination occurs in an intrinsic manner during early development and continues through adulthood in an activity-dependent manner called adaptive myelination. Motor learning, the gradual acquisition of a specific novel motor skill, promotes adaptive myelination in both the healthy and demyelinated central nervous system (CNS). On the other hand, exercise, a physical activity that involves planned, structured and repetitive bodily movements that expend energy and benefits one's fitness, promotes remyelination in pathology, but it is less clear whether it promotes adaptive myelination in healthy subjects. Studies on these topics have also investigated whether the timing of motor learning or physical exercise is important for successful addition of myelin. Here we review our current understanding of the relationship of motor skill learning and physical exercise on myelination.

Quantifying myelin in crossing fibers using diffusion-prepared phase imaging: Theory and simulations

PURPOSE

Myelin has long been the target of neuroimaging research. However, most available techniques can only provide a voxel-averaged estimate of myelin content. In the human brain, white matter fiber pathways connecting different brain areas and carrying different functions often cross each other in the same voxel. A measure that can differentiate the degree of myelination of crossing fibers would provide a more specific marker of myelination.

THEORY AND METHODS

One MRI signal property that is sensitive to myelin is the phase accumulation. This sensitivity is used by measuring the phase accumulation of the signal remaining after diffusion-weighting, which is called diffusion-prepared phase imaging (DIPPI). Including diffusion-weighting before estimating the phase accumulation has two distinct advantages for estimating the degree of myelination: (1) It increases the relative contribution of intra-axonal water, whose phase is related linearly to the thickness of the surrounding myelin (in particular the log g-ratio); and (2) it gives directional information, which can be used to distinguish between crossing fibers. Here the DIPPI sequence is described, an approach is proposed to estimate the log g-ratio, and simulations are used and DIPPI data acquired in an isotropic phantom to quantify other sources of phase accumulation.

RESULTS

The expected bias is estimated in the log g-ratio for reasonable in vivo acquisition parameters caused by eddy currents (~4%-10%), remaining extra-axonal signal (~15%), and gradients in the bulk off-resonance field (<10% for most of the brain).

CONCLUSION

This new sequence may provide a g-ratio estimate per fiber population crossing within a voxel.

Dissecting whole-brain conduction delays through MRI microstructural measures

Network models based on structural connectivity have been increasingly used as the blueprint for large-scale simulations of the human brain. As the nodes of this network are distributed through the cortex and interconnected by white matter pathways with different characteristics, modeling the associated conduction delays becomes important. The goal of this study is to estimate and characterize these delays directly from the brain structure. To achieve this, we leveraged microstructural measures from a combination of advanced magnetic resonance imaging acquisitions and computed the main determinants of conduction velocity, namely axonal diameter and myelin content. Using the model proposed by Rushton, we used these measures to calculate the conduction velocity and estimated the associated delays using tractography. We observed that both the axonal diameter and conduction velocity distributions presented a rather constant trend across different connection lengths, with resulting delays that scale linearly with the connection length. Relying on insights from graph theory and Kuramoto simulations, our results support the approximation of constant conduction velocity but also show path- and region-specific differences.

The Influence of Radio-Frequency Transmit Field Inhomogeneities on the Accuracy of G-ratio Weighted Imaging

G-ratio weighted imaging is a non-invasive, in-vivo MRI-based technique that aims at estimating an aggregated measure of relative myelination of axons across the entire brain white matter. The MR g-ratio and its constituents (axonal and myelin volume fraction) are more specific to the tissue microstructure than conventional MRI metrics targeting either the myelin or axonal compartment. To calculate the MR g-ratio, an MRI-based myelin-mapping technique is combined with an axon-sensitive MR technique (such as diffusion MRI). Correction for radio-frequency transmit (B1+) field inhomogeneities is crucial for myelin mapping techniques such as magnetization transfer saturation. Here we assessed the effect of B1+ correction on g-ratio weighted imaging. To this end, the B1+ field was measured and the B1+ corrected MR g-ratio was used as the reference in a Bland-Altman analysis. We found a substantial bias (≈-89%) and error (≈37%) relative to the dynamic range of g-ratio values in the white matter if the B1+ correction was not applied. Moreover, we tested the efficiency of a data-driven B1+ correction approach that was applied retrospectively without additional reference measurements. We found that it reduced the bias and error in the MR g-ratio by a factor of three. The data-driven correction is readily available in the open-source hMRI toolbox (www.hmri.info) which is embedded in the statistical parameter mapping (SPM) framework.

[Protective effect of transplantation of human oligodendrocyte precursor cells in a rat model of white matter injury]

OBJECTIVE

To study the effect of human oligodendrocyte precursor cell (hOPC) transplantation in the treatment of white matter injury (WMI).

METHODS

Neonatal rats were randomly divided into a sham-operation group, a model group, and a transplantation group (n=10 each). At the age of 3 days, the rats in the model group and the transplantation group were treated with right common carotid artery ligation, followed by hypoxia for 2 hours, to prepare a rat model of WMI. hOPCs were isolated from a spontaneously aborted human fetal brain at week 11 of gestation, and then hOPCs were cultured and transplanted into the rats with WMI. At 3 months after transplantation, the water maze test was performed to evaluate neurological function, and an electron microscope was used to observe myelin sheath thickness and proliferation.

RESULTS

The place navigation test using the Morris water maze showed that the model group had a significantly longer escape latency than the sham-operation group, and compared with the model group, the transplantation group had a significant reduction in escape latency (P < 0.05). To a certain degree, hOPC transplantation alleviated cognitive impairment in rats with WMI at the age of 90 days. The electron microscope images showed that hOPC transplantation promoted remyelination in the brain of WMI rats. Compared with the sham-operation group, the model group had a significant increase in the g-ratio (total axon diameter/total fiber diameter). Compared with the model group, the transplantation group had a significant reduction in the g-ratio (P < 0.05).

CONCLUSIONS

Intrathecal hOPC transplantation may alleviate neurological injury and promote remyelination in a rat model of WMI.

[Protective effect of transplantation of human oligodendrocyte precursor cells in a rat model of white matter injury]

OBJECTIVE

To study the effect of human oligodendrocyte precursor cell (hOPC) transplantation in the treatment of white matter injury (WMI).

METHODS

Neonatal rats were randomly divided into a sham-operation group, a model group, and a transplantation group (n=10 each). At the age of 3 days, the rats in the model group and the transplantation group were treated with right common carotid artery ligation, followed by hypoxia for 2 hours, to prepare a rat model of WMI. hOPCs were isolated from a spontaneously aborted human fetal brain at week 11 of gestation, and then hOPCs were cultured and transplanted into the rats with WMI. At 3 months after transplantation, the water maze test was performed to evaluate neurological function, and an electron microscope was used to observe myelin sheath thickness and proliferation.

RESULTS

The place navigation test using the Morris water maze showed that the model group had a significantly longer escape latency than the sham-operation group, and compared with the model group, the transplantation group had a significant reduction in escape latency (P < 0.05). To a certain degree, hOPC transplantation alleviated cognitive impairment in rats with WMI at the age of 90 days. The electron microscope images showed that hOPC transplantation promoted remyelination in the brain of WMI rats. Compared with the sham-operation group, the model group had a significant increase in the g-ratio (total axon diameter/total fiber diameter). Compared with the model group, the transplantation group had a significant reduction in the g-ratio (P < 0.05).

CONCLUSIONS

Intrathecal hOPC transplantation may alleviate neurological injury and promote remyelination in a rat model of WMI.

Age-related estimates of aggregate g-ratio of white matter structures assessed using quantitative magnetic resonance neuroimaging

The g‐ratio, defined as the inner‐to‐outer diameter of a myelinated axon, is associated with the speed of nerve impulse conduction, and represents an index of axonal myelination and integrity. It has been shown to be a sensitive and specific biomarker of neurodevelopment and neurodegeneration. However, there have been very few magnetic resonance imaging studies of the g‐ratio in the context of normative aging; characterizing regional and time‐dependent cerebral changes in g‐ratio in cognitively normal subjects will be a crucial step in differentiating normal from abnormal microstructural alterations. In the current study, we investigated age‐related differences in aggregate g‐ratio, that is, g‐ratio averaged over all fibers within regions of interest, in several white matter regions in a cohort of 52 cognitively unimpaired participants ranging in age from 21 to 84 years. We found a quadratic, U‐shaped, relationship between aggregate g‐ratio and age in most cerebral regions investigated, suggesting myelin maturation until middle age followed by a decrease at older ages. As expected, we observed that these age‐related differences vary across different brain regions, with the frontal lobes and parietal lobes exhibiting slightly earlier ages of minimum aggregate g‐ratio as compared to more posterior structures such as the occipital lobes and temporal lobes; this agrees with the retrogenesis paradigm. Our results provide evidence for a nonlinear association between age and aggregate g‐ratio in a sample of adults from a highly controlled population. Finally, sex differences in aggregate g‐ratio were observed in several cerebral regions, with women exhibiting overall lower values as compared to men; this likely reflects the greater myelin content in women's brain, in agreement with recent investigations.

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.

Realistic Microstructure Simulator (RMS): Monte Carlo simulations of diffusion in three-dimensional cell segmentations of microscopy images

BACKGROUND

Monte Carlo simulations of diffusion are commonly used as a model validation tool as they are especially suitable for generating the diffusion MRI signal in complicated tissue microgeometries.

NEW METHOD

Here we describe the details of implementing Monte Carlo simulations in three-dimensional (3d) voxelized segmentations of cells in microscopy images. Using the concept of the corner reflector, we largely reduce the computational load of simulating diffusion within and exchange between multiple cells. Precision is further achieved by GPU-based parallel computations.

RESULTS

Our simulation of diffusion in white matter axons segmented from a mouse brain demonstrates its value in validating biophysical models. Furthermore, we provide the theoretical background for implementing a discretized diffusion process, and consider the finite-step effects of the particle-membrane reflection and permeation events, needed for efficient simulation of interactions with irregular boundaries, spatially variable diffusion coefficient, and exchange.

COMPARISON WITH EXISTING METHODS

To our knowledge, this is the first Monte Carlo pipeline for MR signal simulations in a substrate composed of numerous realistic cells, accounting for their permeable and irregularly-shaped membranes.

CONCLUSIONS

The proposed RMS pipeline makes it possible to achieve fast and accurate simulations of diffusion in realistic tissue microgeometry, as well as the interplay with other MR contrasts. Presently, RMS focuses on simulations of diffusion, exchange, and T₁ and T₂ NMR relaxation in static tissues, with a possibility to straightforwardly account for susceptibility-induced T₂^* effects and flow.

The impact of realistic axonal shape on axon diameter estimation using diffusion MRI

To study axonal microstructure with diffusion MRI, axons are typically modeled as straight impermeable cylinders, whereby the transverse diffusion MRI signal can be made sensitive to the cylinder’s inner diameter. However, the shape of a real axon varies along the axon direction, which couples the longitudinal and transverse diffusion of the overall axon direction. Here we develop a theory of the intra-axonal diffusion MRI signal based on coarse-graining of the axonal shape by 3-dimensional diffusion. We demonstrate how the estimate of the inner diameter is confounded by the diameter variations (beading), and by the local variations in direction (undulations) along the axon. We analytically relate diffusion MRI metrics, such as time-dependent radial diffusivity D_⊥(t) and kurtosis K_⊥(t), to the axonal shape, and validate our theory using Monte Carlo simulations in synthetic undulating axons with randomly positioned beads, and in realistic axons reconstructed from electron microscopy images of mouse brain white matter. We show that (i) In the narrow pulse limit, the inner diameter from D_⊥(t) is overestimated by about twofold due to a combination of axon caliber variations and undulations (each contributing a comparable effect size); (ii) The narrow-pulse kurtosis K_⊥∣_t→∞ deviates from that in an ideal cylinder due to caliber variations; we also numerically calculate the fourth-order cumulant for an ideal cylinder in the wide pulse limit, which is relevant for inner diameter overestimation; (iii) In the wide pulse limit, the axon diameter overestimation is mainly due to undulations at low diffusion weightings b; and (iv) The effect of undulations can be considerably reduced by directional averaging of high-b signals, with the apparent inner diameter given by a combination of the axon caliber (dominated by the thickest axons), caliber variations, and the residual contribution of undulations.

Cognitive impairments induced by necrotizing enterocolitis can be prevented by inhibiting microglial activation in mouse brain

Necrotizing enterocolitis (NEC) is a severe gastrointestinal disease of the premature infant. One of the most important long-term complications observed in children who survive NEC early in life is the development of profound neurological impairments. However, the pathways leading to NEC-associated neurological impairments remain unknown, thus limiting the development of prevention strategies. We have recently shown that NEC development is dependent on the expression of the lipopolysaccharide receptor Toll-like receptor 4 (TLR4) on the intestinal epithelium, whose activation by bacteria in the newborn gut leads to mucosal inflammation. Here, we hypothesized that damage-induced production of TLR4 endogenous ligands in the intestine might lead to activation of microglial cells in the brain and promote cognitive impairments. We identified a gut-brain signaling axis in an NEC mouse model in which activation of intestinal TLR4 signaling led to release of high-mobility group box 1 in the intestine that, in turn, promoted microglial activation in the brain and neurological dysfunction. We further demonstrated that an orally administered dendrimer-based nanotherapeutic approach to targeting activated microglia could prevent NEC-associated neurological dysfunction in neonatal mice. These findings shed light on the molecular pathways leading to the development of NEC-associated brain injury, provide a rationale for early removal of diseased intestine in NEC, and indicate the potential of targeted therapies that protect the developing brain in the treatment of NEC in early childhood.

Along-axon diameter variation and axonal orientation dispersion revealed with 3D electron microscopy: implications for quantifying brain white matter microstructure with histology and diffusion MRI

Tissue microstructure modeling of diffusion MRI signal is an active research area striving to bridge the gap between macroscopic MRI resolution and cellular-level tissue architecture. Such modeling in neuronal tissue relies on a number of assumptions about the microstructural features of axonal fiber bundles, such as the axonal shape (e.g., perfect cylinders) and the fiber orientation dispersion. However, these assumptions have not yet been validated by sufficiently high-resolution 3-dimensional histology. Here, we reconstructed sequential scanning electron microscopy images in mouse brain corpus callosum, and introduced a random-walker (RaW)-based algorithm to rapidly segment individual intra-axonal spaces and myelin sheaths of myelinated axons. Confirmed by a segmentation based on human annotations initiated with conventional machine-learning-based carving, our semi-automatic algorithm is reliable and less time-consuming. Based on the segmentation, we calculated MRI-relevant estimates of size-related parameters (inner axonal diameter, its distribution, along-axon variation, and myelin g-ratio), and orientation-related parameters (fiber orientation distribution and its rotational invariants; dispersion angle). The reported dispersion angle is consistent with previous 2-dimensional histology studies and diffusion MRI measurements, while the reported diameter exceeds those in other mouse brain studies. Furthermore, we calculated how these quantities would evolve in actual diffusion MRI experiments as a function of diffusion time, thereby providing a coarse-graining window on the microstructure, and showed that the orientation-related metrics have negligible diffusion time-dependence over clinical and pre-clinical diffusion time ranges. However, the MRI-measured inner axonal diameters, dominated by the widest cross sections, effectively decrease with diffusion time by ~ 17% due to the coarse-graining over axonal caliber variations. Furthermore, our 3d measurement showed that there is significant variation of the diameter along the axon. Hence, fiber orientation dispersion estimated from MRI should be relatively stable, while the "apparent" inner axonal diameters are sensitive to experimental settings, and cannot be modeled by perfectly cylindrical axons.

Automated 3D Axonal Morphometry of White Matter

Axonal structure underlies white matter functionality and plays a major role in brain connectivity. The current literature on the axonal structure is based on the analysis of two-dimensional (2D) cross-sections, which, as we demonstrate, is precarious. To be able to quantify three-dimensional (3D) axonal morphology, we developed a novel pipeline, called ACSON (AutomatiC 3D Segmentation and morphometry Of axoNs), for automated 3D segmentation and morphometric analysis of the white matter ultrastructure. The automated pipeline eliminates the need for time-consuming manual segmentation of 3D datasets. ACSON segments myelin, myelinated and unmyelinated axons, mitochondria, cells and vacuoles, and analyzes the morphology of myelinated axons. We applied the pipeline to serial block-face scanning electron microscopy images of the corpus callosum of sham-operated (n = 2) and brain injured (n = 3) rats 5 months after the injury. The 3D morphometry showed that cross-sections of myelinated axons were elliptic rather than circular, and their diameter varied substantially along their longitudinal axis. It also showed a significant reduction in the myelinated axon diameter of the ipsilateral corpus callosum of rats 5 months after brain injury, indicating ongoing axonal alterations even at this chronic time-point.

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.

Myelin volume fraction imaging with MRI

MRI is a valuable tool to assess myelin during development and demyelinating disease processes. While multiexponential T2 and quantitative magnetization transfer measures correlate with myelin content, neither provides the total myelin volume fraction. In many cases correlative measures are adequate; but to assess microstructure of myelin, (e.g. calculate the g-ratio using MRI), an accurate measure of myelin volume fraction is imperative. Using a volumetric model of white matter, we relate MRI measures of myelin to absolute measures of myelin volume fraction and compare them to quantitative histology. We assess our approach in control mice along with two models of hypomyelination and one model of hypermyelination and find strong agreement between MRI and histology amongst models. This work investigates the sensitivities of MRI myelin measures to changes in axon geometry and displays promise for estimating g-ratio from MRI.

Four in vivo g-ratio-weighted imaging methods: Comparability and repeatability at the group level

A recent method, denoted in vivo g-ratio-weighted imaging, has related the microscopic g-ratio, only accessible by ex vivo histology, to noninvasive MRI markers for the fiber volume fraction (FVF) and myelin volume fraction (MVF). Different MRI markers have been proposed for g-ratio weighted imaging, leaving open the question which combination of imaging markers is optimal. To address this question, the repeatability and comparability of four g-ratio methods based on different combinations of, respectively, two imaging markers for FVF (tract-fiber density, TFD, and neurite orientation dispersion and density imaging, NODDI) and two imaging markers for MVF (magnetization transfer saturation rate, MT, and, from proton density maps, macromolecular tissue volume, MTV) were tested in a scan-rescan experiment in two groups. Moreover, it was tested how the repeatability and comparability were affected by two key processing steps, namely the masking of unreliable voxels (e.g., due to partial volume effects) at the group level and the calibration value used to link MRI markers to MVF (and FVF). Our data showed that repeatability and comparability depend largely on the marker for the FVF (NODDI outperformed TFD), and that they were improved by masking. Overall, the g-ratio method based on NODDI and MT showed the highest repeatability (90%) and lowest variability between groups (3.5%). Finally, our results indicate that the calibration procedure is crucial, for example, calibration to a lower g-ratio value (g = 0.6) than the commonly used one (g = 0.7) can change not only repeatability and comparability but also the reported dependency on the FVF imaging marker. Hum Brain Mapp 39:24-41, 2018. © 2017 Wiley Periodicals, Inc.

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.

Whole brain g-ratio mapping using myelin water imaging (MWI) and neurite orientation dispersion and density imaging (NODDI)

MR g-ratio, which measures the ratio of the aggregate volume of axons to that of fibers in a voxel, is a potential biomarker for white matter microstructures. In this study, a new approach for acquiring an in-vivo whole human brain g-ratio map is proposed. To estimate the g-ratio, myelin volume fraction and axonal volume fraction are acquired using multi-echo gradient echo myelin water imaging (GRE-MWI) and neurite orientation dispersion and density imaging (NODDI), respectively. In order to translate myelin water fraction measured in GRE-MWI into myelin volume fraction, a new scaling procedure is proposed and validated. This scaling approach utilizes geometric measures of myelin structure and, therefore, provides robustness over previous methods. The resulting g-ratio map reveals an expected range of g-ratios (0.71-0.85 in major fiber bundles) with a small inter-subject coefficient of variance (less than 2%). Additionally, a few fiber bundles (e.g. cortico-spinal tract and optic radiation) show different constituents of myelin volume fraction and axonal volume fraction, indicating potentials to utilize the measures for deciphering fiber tracking.

Introducing axonal myelination in connectomics: A preliminary analysis of g-ratio distribution in healthy subjects

Microstructural imaging and connectomics are two research areas that hold great potential for investigating brain structure and function. Combining these two approaches can lead to a better and more complete characterization of the brain as a network. The aim of this work is characterizing the connectome from a novel perspective using the myelination measure given by the g-ratio. The g-ratio is the ratio of the inner to the outer diameters of a myelinated axon, whose aggregated value can now be estimated in vivo using MRI. In two different datasets of healthy subjects, we reconstructed the structural connectome and then used the g-ratio estimated from diffusion and magnetization transfer data to characterize the network structure. Significant characteristics of g-ratio weighted graphs emerged. First, the g-ratio distribution across the edges of the graph did not show the power-law distribution observed using the number of streamlines as a weight. Second, connections involving regions related to motor and sensory functions were the highest in myelin content. We also observed significant differences in terms of the hub structure and the rich-club organization suggesting that connections involving hub regions present higher myelination than peripheral connections. Taken together, these findings offer a characterization of g-ratio distribution across the connectome in healthy subjects and lay the foundations for further investigating plasticity and pathology using a similar approach.

Mapping White Matter Microstructure in the One Month Human Brain

White matter microstructure, essential for efficient and coordinated transmission of neural communications, undergoes pronounced development during the first years of life, while deviations to this neurodevelopmental trajectory likely result in alterations of brain connectivity relevant to behavior. Hence, systematic evaluation of white matter microstructure in the normative brain is critical for a neuroscientific approach to both typical and atypical early behavioral development. However, few studies have examined the infant brain in detail, particularly in infants under 3 months of age. Here, we utilize quantitative techniques of diffusion tensor imaging and neurite orientation dispersion and density imaging to investigate neonatal white matter microstructure in 104 infants. An optimized multiple b-value diffusion protocol was developed to allow for successful acquisition during non-sedated sleep. Associations between white matter microstructure measures and gestation corrected age, regional asymmetries, infant sex, as well as newborn growth measures were assessed. Results highlight changes of white matter microstructure during the earliest periods of development and demonstrate differential timing of developing regions and regional asymmetries. Our results contribute to a growing body of research investigating the neurobiological changes associated with neurodevelopment and suggest that characteristics of white matter microstructure are already underway in the weeks immediately following birth.

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.

Modeling white matter microstructure

Quantitative magnetic resonance imaging can be combined with advanced biophysical models to measure microstructural features of white matter. Non-invasive microstructural imaging has the potential to revolutionize neuroscience, and acquiring these measures in clinically feasible times would greatly improve patient monitoring and clinical studies of drug efficacy. However, a good understanding of microstructural imaging techniques is essential to set realistic expectations and to prevent over-interpretation of results. This review explains the methodology behind microstructural modeling and imaging, and gives an overview of the breakthroughs and challenges associated with it.

g-Ratio weighted imaging of the human spinal cord in vivo

The fiber g-ratio is defined as the ratio of the inner to the outer diameter of the myelin sheath. This ratio provides a measure of the myelin thickness that complements axon morphology (diameter and density) for assessment of demyelination in diseases such as multiple sclerosis. Previous work has shown that an aggregate g-ratio map can be computed using a formula that combines axon and myelin density measured with quantitative MRI. In this work, we computed g-ratio weighted maps in the cervical spinal cord of nine healthy subjects. We utilized the 300mT/m gradients from the CONNECTOM scanner to estimate the fraction of restricted water (fr) with high accuracy, using the CHARMED model. Myelin density was estimated using the lipid and macromolecular tissue volume (MTV) method, derived from normalized proton density (PD) mapping. The variability across spinal level, laterality and subject were assessed using a three-way ANOVA. The average g-ratio value obtained in the white matter was 0.76+/-0.03, consistent with previous histology work. Coefficients of variation of fr and MTV were respectively 4.3% and 13.7%. fr and myelin density were significantly different across spinal tracts (p=3×10-7 and 0.004 respectively) and were positively correlated in the white matter (r=0.42), suggesting shared microstructural information. The aggregate g-ratio did not show significant differences across tracts (p=0.6). This study suggests that fr and myelin density can be measured in vivo with high precision and that they can be combined to produce a g-ratio-weighted map robust to free water pool contamination from cerebrospinal fluid or veins. Potential applications include the study of early demyelination in multiple sclerosis, and the quantitative assessment of remyelination drugs.

The role of myelination in measures of white matter integrity: Combination of diffusion tensor imaging and two-photon microscopy of CLARITY intact brains

Diffusion tensor imaging (DTI) is used extensively in neuroscience to noninvasively estimate white matter (WM) microarchitecture. However, the diffusion signal is inherently ambiguous because it infers WM structure from the orientation of water diffusion and cannot identify the biological sources of diffusion changes. To compare inferred WM estimates to directly labeled axonal elements, we performed a novel within-subjects combination of high-resolution ex vivo DTI with two-photon laser microscopy of intact mouse brains rendered optically transparent by Clear Lipid-exchanged, Anatomically Rigid, Imaging/immunostaining compatible, Tissue hYdrogel (CLARITY). We found that myelin basic protein (MBP) immunofluorescence significantly correlated with fractional anisotropy (FA), especially in WM regions with coherent fiber orientations and low fiber dispersion. Our results provide evidence that FA is particularly sensitive to myelination in WM regions with these characteristics. Furthermore, we found that radial diffusivity (RD) was only sensitive to myelination in a subset of WM tracts, suggesting that the association of RD with myelin should be used cautiously. This combined DTI-CLARITY approach illustrates, for the first time, a framework for using brain-wide immunolabeling of WM targets to elucidate the relationship between the diffusion signal and its biological underpinnings. This study also demonstrates the feasibility of a within-subject combination of noninvasive neuroimaging and tissue clearing techniques that has broader implications for neuroscience research.

Characterizing axonal myelination within the healthy population: a tract-by-tract mapping of effects of age and gender on the fiber g-ratio

The g-ratio, equal to the ratio of the inner-to-outer diameter of a myelinated axon, is associated with the speed of conduction, and thus reflects axonal function and integrity. It is now possible to estimate an “aggregate” g-ratio in vivo using MRI. The aim of this study was to assess the variation of the MRI-derived fiber g-ratio in the brain of healthy individuals, and to characterize its variation across the lifespan. Thirty-eight healthy participants, aged between 20 and 76, were recruited. Whole-brain g-ratio maps were computed and analyzed voxel-wise. Median tract g-ratio values were also extracted. No significant effect of gender was found, whereas age was found to be significantly associated with the g-ratio within the white matter. The tract-specific analysis showed this relationship to follow a nearly-linear increase, although the slope appears to slow down slightly after the 6th decade of life. The most likely interpretation is a subtle but consistent reduction in myelin throughout adulthood, with the density of axons beginning to decrease between the 4th and 5th decade.

Whole-Brain In-vivo Measurements of the Axonal G-Ratio in a Group of 37 Healthy Volunteers

The g-ratio, quantifying the ratio between the inner and outer diameters of a fiber, is an important microstructural characteristic of fiber pathways and is functionally related to conduction velocity. We introduce a novel method for estimating the MR g-ratio non-invasively across the whole brain using high-fidelity magnetization transfer (MT) imaging and single-shell diffusion MRI. These methods enabled us to map the MR g-ratio in vivo across the brain's prominent fiber pathways in a group of 37 healthy volunteers and to estimate the inter-subject variability. Effective correction of susceptibility-related distortion artifacts was essential before combining the MT and diffusion data, in order to reduce partial volume and edge artifacts. The MR g-ratio is in good qualitative agreement with histological findings despite the different resolution and spatial coverage of MRI and histology. The MR g-ratio holds promise as an important non-invasive biomarker due to its microstructural and functional relevance in neurodegeneration.