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

Evaluation of tractography-based myelin-weighted connectivity across the lifespan

Introduction

Recent studies showed that the myelin of the brain changes in the life span, and demyelination contributes to the loss of brain plasticity during normal aging. Diffusion-weighted magnetic resonance imaging (dMRI) allows studying brain connectivity in vivo by mapping axons in white matter with tractography algorithms. However, dMRI does not provide insight into myelin; thus, combining tractography with myelin-sensitive maps is necessary to investigate myelin-weighted brain connectivity. Tractometry is designated for this purpose, but it suffers from some serious limitations. Our study assessed the effectiveness of the recently proposed Myelin Streamlines Decomposition (MySD) method in estimating myelin-weighted connectomes and its capacity to detect changes in myelin network architecture during the process of normal aging. This approach opens up new possibilities compared to traditional Tractometry.

Methods

In a group of 85 healthy controls aged between 18 and 68 years, we estimated myelin-weighted connectomes using Tractometry and MySD, and compared their modulation with age by means of three well-known global network metrics.

Results

Following the literature, our results show that myelin development continues until brain maturation (40 years old), after which degeneration begins. In particular, mean connectivity strength and efficiency show an increasing trend up to 40 years, after which the process reverses. Both Tractometry and MySD are sensitive to these changes, but MySD turned out to be more accurate.

Conclusion

After regressing the known predictors, MySD results in lower residual error, indicating that MySD provides more accurate estimates of myelin-weighted connectivity than Tractometry.

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.

Oligodendrocyte Progenitor Cell Transplantation Ameliorates Preterm Infant Cerebral White Matter Injury in Rats Model

Background

Cerebral white matter injury (WMI) is the most common brain injury in preterm infants, leading to motor and developmental deficits often accompanied by cognitive impairment. However, there is no effective treatment. One promising approach for treating preterm WMI is cell replacement therapy, in which lost cells can be replaced by exogenous oligodendrocyte progenitor cells (OPCs).

Methods

This study developed a method to differentiate human neural stem cells (hNSCs) into human OPCs (hOPCs). The preterm WMI animal model was established in rats on postnatal day 3, and OLIG2+/NG2+/PDGFRα+/O4+ hOPCs were enriched and transplanted into the corpus callosum on postnatal day 10. Then, histological analysis and electron microscopy were used to detect lesion structure; behavioral assays were performed to detect cognitive function.

Results

Transplanted hOPCs survived and migrated throughout the major white matter tracts. Morphological differentiation of transplanted hOPCs was observed. Histological analysis revealed structural repair of lesioned areas. Re-myelination of the axons in the corpus callosum was confirmed by electron microscopy. The Morris water maze test revealed cognitive function recovery.

Conclusion

Our study showed that exogenous hOPCs could differentiate into CC1+ OLS in the brain of WMI rats, improving their cognitive functions.

Fiber-orientation independent component of R2* obtained from single-orientation MRI measurements in simulations and a post-mortem human optic chiasm

The effective transverse relaxation rate (R2*) is sensitive to the microstructure of the human brain like the g-ratio which characterises the relative myelination of axons. However, the fibre-orientation dependence of R2* degrades its reproducibility and any microstructural derivative measure. To estimate its orientation-independent part (R2,iso*) from single multi-echo gradient-recalled-echo (meGRE) measurements at arbitrary orientations, a second-order polynomial in time model (hereafter M2) can be used. Its linear time-dependent parameter, β1, can be biophysically related to R2,iso* when neglecting the myelin water (MW) signal in the hollow cylinder fibre model (HCFM). Here, we examined the performance of M2 using experimental and simulated data with variable g-ratio and fibre dispersion. We found that the fitted β1 can estimate R2,iso* using meGRE with long maximum-echo time (TEmax ≈ 54 ms), but not accurately captures its microscopic dependence on the g-ratio (error 84%). We proposed a new heuristic expression for β1 that reduced the error to 12% for ex vivo compartmental R2 values. Using the new expression, we could estimate an MW fraction of 0.14 for fibres with negligible dispersion in a fixed human optic chiasm for the ex vivo compartmental R2 values but not for the in vivo values. M2 and the HCFM-based simulations failed to explain the measured R2*-orientation-dependence around the magic angle for a typical in vivo meGRE protocol (with TEmax ≈ 18 ms). In conclusion, further validation and the development of movement-robust in vivo meGRE protocols with TEmax ≈ 54 ms are required before M2 can be used to estimate R2,iso* in subjects.

Maternal antibiotic administration during gestation can affect the memory and brain structure in mouse offspring

Maternal antibiotics administration (MAA) is among the widely used therapeutic approaches in pregnancy. Although published evidence demonstrates that infants exposed to antibiotics immediately after birth have altered recognition memory responses at one month of age, very little is known about in utero effects of antibiotics on the neuronal function and behavior of children after birth. Therefore, this study aimed to evaluate the impact of MAA at different periods of pregnancy on memory decline and brain structural alterations in young mouse offspring after their first month of life. To study the effects of MAA on 4-week-old offspring, pregnant C57BL/6J mouse dams (2-3-month-old; n = 4/group) were exposed to a cocktail of amoxicillin (205 mg/kg/day) and azithromycin (51 mg/kg/day) in sterile drinking water (daily/1 week) during either the 2nd or 3rd week of pregnancy and stopped after delivery. A control group of pregnant dams was exposed to sterile drinking water alone during all three weeks of pregnancy. Then, the 4-week-old offspring mice were first evaluated for behavioral changes. Using the Morris water maze assay, we revealed that exposure of pregnant mice to antibiotics at the 2nd and 3rd weeks of pregnancy significantly altered spatial reference memory and learning skills in their offspring compared to those delivered from the control group of dams. In contrast, no significant difference in long-term associative memory was detected between offspring groups using the novel object recognition test. Then, we histologically evaluated brain samples from the same offspring individuals using conventional immunofluorescence and electron microscopy assays. To our knowledge, we observed a reduction in the density of the hippocampal CA1 pyramidal neurons and hypomyelination in the corpus callosum in groups of mice in utero exposed to antibiotics at the 2nd and 3rd weeks of gestation. In addition, offspring exposed to antibiotics at the 2nd or 3rd week of gestation demonstrated a decreased astrocyte cell surface area and astrocyte territories or depletion of neurogenesis in the dentate gyrus and hippocampal synaptic loss, respectively. Altogether, this study shows that MAA at different times of pregnancy can pathologically alter cognitive behavior and brain development in offspring at an early age after weaning.

The Potential of Myelin-Sensitive Imaging: Redefining Spatiotemporal Patterns of Myeloarchitecture

Recent advances in magnetic resonance imaging (MRI) have paved the way for approximation of myelin content in vivo. In this review, our main goal was to determine how to best capitalize on myelin-sensitive imaging. First, we briefly overview the theoretical and empirical basis for the myelin sensitivity of different MRI markers and, in doing so, highlight how multimodal imaging approaches are important for enhancing specificity to myelin. Then, we discuss recent studies that have probed the nonuniform distribution of myelin across cortical layers and along white matter tracts. These approaches, collectively known as myelin profiling, have provided detailed depictions of myeloarchitecture in both the postmortem and living human brain. Notably, MRI-based profiling studies have recently focused on investigating whether it can capture interindividual variability in myelin characteristics as well as trajectories across the lifespan. Finally, another line of recent evidence emphasizes the contribution of region-specific myelination to large-scale organization, demonstrating the impact of myelination on global brain networks. In conclusion, we suggest that combining well-validated MRI markers with profiling techniques holds strong potential to elucidate individual differences in myeloarchitecture, which has important implications for understanding brain function and disease.

Towards automated in vivo parcellation of the human cerebral cortex using supervised classification of magnetic resonance fingerprinting residuals

BACKGROUND

Accurate parcellation of the cerebral cortex in an individual is a guide to its underlying organization. The most promising in vivo quantitative magnetic resonance (MR)-based microstructural cortical mapping methods are yet to achieve a level of parcellation accuracy comparable to quantitative histology.

METHODS

We scanned 6 participants using a 3D echo-planar imaging MR fingerprinting (EPI-MRF) sequence on a 7T Siemens scanner. After projecting MRF signals to the individual-specific inflated model of the cortical surface, normalized autocorrelations of MRF residuals of vertices of 8 microstructurally distinct areas (BA1, BA2, BA4a, BA6, BA44, BA45, BA17, and BA18) from 3 cortical regions were used as feature vector inputs into linear support vector machine (SVM), radial basis function SVM (RBF-SVM), random forest, and k-nearest neighbors supervised classification algorithms. The algorithms' prediction performance was compared using: (i) features from each vertex or (ii) features from neighboring vertices.

RESULTS

The neighborhood-based RBF-SVM classifier achieved the highest prediction score of 0.85 for classification of MRF residuals in the central region from a held-out participant.

CONCLUSIONS

We developed an automated method of cortical parcellation using a combination of MR fingerprinting residual analysis and machine learning classification. Our findings provide the basis for employing unsupervised learning algorithms for whole-cortex structural parcellation in individuals.

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.

Longitudinal microstructural MRI markers of demyelination and neurodegeneration in early relapsing-remitting multiple sclerosis: Magnetisation transfer, water diffusion and g-ratio

INTRODUCTION

Quantitative microstructural MRI, such as myelin-sensitive magnetisation transfer ratio (MTR) or saturation (MTsat), axon-sensitive water diffusion Neurite Orientation Dispersion and Density Imaging (NODDI), and the aggregate g-ratio, may provide more specific markers of white matter integrity than conventional MRI for early patient stratification in relapsing-remitting multiple sclerosis (RRMS). The aim of this study was to determine the sensitivity of such markers to longitudinal pathological change within cerebral white matter lesions (WML) and normal-appearing white matter (NAWM) in recently diagnosed RRMS.

METHODS

Seventy-nine people with recently diagnosed RRMS, from the FutureMS longitudinal cohort, were recruited to an extended MRI protocol at baseline and one year later. Twelve healthy volunteers received the same MRI protocol, repeated within two weeks. Ethics approval and written informed consent were obtained. 3T MRI included magnetisation transfer, and multi-shell diffusion-weighted imaging. NAWM and whole brain were segmented from 3D T1-weighted MPRAGE, and WML from T2-weighted FLAIR. MTR, MTsat, NODDI isotropic (ISOVF) and intracellular (ICVF) volume fractions, and g-ratio (calculated from MTsat and NODDI data) were measured within WML and NAWM. Brain parenchymal fraction (BPF) was also calculated. Longitudinal change in BPF and microstructural metrics was assessed with paired t-tests (α = 0.05) and linear mixed models, adjusted for confounding factors with False Discovery Rate (FDR) correction for multiple comparisons. Longitudinal changes were compared with test-retest Bland-Altman limits of agreement from healthy control white matter. The influence of longitudinal change on g-ratio was explored through post-hoc analysis in silico by computing g-ratio with realistic simulated MTsat and NODDI values.

RESULTS

In NAWM, g-ratio and ICVF increased, and MTsat decreased over one year (adjusted mean difference = 0.007, 0.005, and -0.057 respectively, all FDR-corrected p < 0.05). There was no significant change in MTR, ISOVF, or BPF. In WML, MTsat, NODDI ICVF and ISOVF increased over time (adjusted mean difference = 0.083, 0.024 and 0.016, respectively, all FDR-corrected p < 0.05). Group-level longitudinal changes exceeded test-retest limits of agreement for NODDI ISOVF and ICVF in WML only. In silico analysis showed g-ratio may increase due to a decrease in MTsat or ISOVF, or an increase in ICVF.

DISCUSSION

G-ratio and MTsat changes in NAWM over one year may indicate subtle myelin loss in early RRMS, which were not apparent with BPF or NAWM MTR. Increases in NAWM and WML NODDI ICVF were not anticipated, and raise the possibility of axonal swelling or morphological change. Increases in WML MTsat may reflect myelin repair. Changes in NODDI ISOVF are more likely to reflect alterations in water content. Competing MTsat and ICVF changes may account for the absence of g-ratio change in WML. Longitudinal changes in microstructural measures are significant at a group level, however detection in individual patients in early RRMS is limited by technique reproducibility.

CONCLUSION

MTsat and g-ratio are more sensitive than MTR to early pathological changes in RRMS, but complex dependence of g-ratio on NODDI parameters limit the interpretation of aggregate measures in isolation. Improvements in technique reproducibility and validation of MRI biophysical models across a range of pathological tissue states are needed.

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.

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

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

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.

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.

Advanced Diffusion MR Imaging for Multiple Sclerosis in the Brain and Spinal Cord

Diffusion tensor imaging (DTI) has been established its usefulness in evaluating normal-appearing white matter (NAWM) and other lesions that are difficult to evaluate with routine clinical MRI in the evaluation of the brain and spinal cord lesions in multiple sclerosis (MS), a demyelinating disease. With the recent advances in the software and hardware of MRI systems, increasingly complex and sophisticated MRI and analysis methods, such as q-space imaging, diffusional kurtosis imaging, neurite orientation dispersion and density imaging, white matter tract integrity, and multiple diffusion encoding, referred to as advanced diffusion MRI, have been proposed. These are capable of capturing in vivo microstructural changes in the brain and spinal cord in normal and pathological states in greater detail than DTI.

This paper reviews the current status of recent advanced diffusion MRI for assessing MS in vivo as part of an issue celebrating two decades of magnetic resonance in medical sciences (MRMS), an official journal of the Japanese Society of Magnetic Resonance in Medicine.

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.

Rationale and design of the brain magnetic resonance imaging protocol for FutureMS: a longitudinal multi-centre study of newly diagnosed patients with relapsing-remitting multiple sclerosis in Scotland

Introduction: Multiple sclerosis (MS) is a chronic neuroinflammatory and neurodegenerative disease. MS prevalence varies geographically and is notably high in Scotland. Disease trajectory varies significantly between individuals and the causes for this are largely unclear. Biomarkers predictive of disease course are urgently needed to allow improved stratification for current disease modifying therapies and future targeted treatments aimed at neuroprotection and remyelination. Magnetic resonance imaging (MRI) can detect disease activity and underlying damage non-invasively in vivo at the micro and macrostructural level. FutureMS is a prospective Scottish longitudinal multi-centre cohort study, which focuses on deeply phenotyping patients with recently diagnosed relapsing-remitting MS (RRMS). Neuroimaging is a central component of the study and provides two main primary endpoints for disease activity and neurodegeneration. This paper provides an overview of MRI data acquisition, management and processing in FutureMS. FutureMS is registered with the Integrated Research Application System (IRAS, UK) under reference number 169955. Methods and analysis: MRI is performed at baseline (N=431) and 1-year follow-up, in Dundee, Glasgow and Edinburgh (3T Siemens) and in Aberdeen (3T Philips), and managed and processed in Edinburgh. The core structural MRI protocol comprises T1-weighted, T2-weighted, FLAIR and proton density images. Primary imaging outcome measures are new/enlarging white matter lesions (WML) and reduction in brain volume over one year. Secondary imaging outcome measures comprise WML volume as an additional quantitative structural MRI measure, rim lesions on susceptibility-weighted imaging, and microstructural MRI measures, including diffusion tensor imaging and neurite orientation dispersion and density imaging metrics, relaxometry, magnetisation transfer (MT) ratio, MT saturation and derived g-ratio measures. Conclusions: FutureMS aims to reduce uncertainty around disease course and allow for targeted treatment in RRMS by exploring the role of conventional and advanced MRI measures as biomarkers of disease severity and progression in a large population of RRMS patients in Scotland.

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

BACKGROUND

Wave-CAIPI Visualization of Short Transverse relaxation time component (ViSTa) is a recently developed, short-T₁-sensitized MRI method for fast quantification of myelin water fraction (MWF) in the human brain. It represents a promising technique for the evaluation of subtle, early signals of demyelination in the cerebral white matter of multiple sclerosis (MS) patients. Currently however, few studies exist that robustly assess the utility of ViSTa MWF measures of myelin compared to more conventional MRI measures of myelin in the brain of MS patients. Moreover, there are no previous studies evaluating the sensitivity of ViSTa MWF for the non-invasive detection of subtle tissue damage in both normal-appearing white matter (NAWM) and white matter lesions of MS patients. As a result, a central purpose of this study was to systematically evaluate the relationship between myelin sensitivity of T₁-based ViSTa MWF mapping and a more generally recognized metric, Magnetization Transfer Saturation (MT_sat), in healthy control and MS brain white matter.

METHODS

ViSTa MWF and MT_sat values were evaluated in automatically-classified normal appearing white matter (NAWM), white matter (WM) lesion tissue, cortical gray matter, and deep gray matter of 29 MS patients and 10 healthy controls using 3T MRI. MWF and MT _sat were also assessed in a tract-specific manner using the Johns Hopkins University WM atlas. MRI-derived measures of cerebral myelin content were uniquely compared by employing non-normal distribution-specific measures of median, interquartile range and skewness. Separate analyses of variance were applied to test tissue-specific differences in MT_sat and ViSTa MWF distribution metrics. Non-parametric tests were utilized when appropriate. All tests were corrected for multiple comparisons using the False Discovery Rate method at the level, α=0.05.

RESULTS

Differences in whole NAWM MS tissue damage were detected with a higher effect size when using ViSTa MWF (q = 0.0008; ƞ² = 0.34) compared to MT_sat (q = 0.02; ƞ²= 0.24). We also observed that, as a possible measure of WM pathology, ViSTa-derived NAWM MWF voxel distributions of MS subjects were consistently skewed towards lower MWF values, while MT_sat voxel distributions showed reduced skewness values. We further identified tract-specific reductions in mean ViSTa MWF of MS patients compared to controls that were not observed with MT_sat. However, MT_sat (q = 1.4 × 10^-21; ƞ² = 0.88) displayed higher effect sizes when differentiating NAWM and MS lesion tissue. Using regression analysis at the group level, we identified a linear relationship between MT_sat and ViSTa MWF in NAWM (R² = 0.46; p = 7.8 × 10^-4) lesions (R² = 0.30; p = 0.004), and with all tissue types combined (R² = 0.71; p = 8.4 × 10^-45). The linear relationship was also observed in most of the WM tracts we investigated. ViSTa MWF in NAWM of MS patients correlated with both disease duration (p = 0.02; R² = 0.27) and WM lesion volume (p = 0.002; R² = 0.34).

CONCLUSION

Because ViSTa MWF and MT_sat metrics exhibit differential sensitivities to tissue damage in MS white matter, they can be collected in combination to provide an efficient, comprehensive measure of myelin water and macromolecular pool proton signals. These complementary measures may offer a more sensitive, non-invasive biopsy of early precursor signals in NAWM that occur prior to lesion formation. They may also aid in monitoring the efficacy of remyelination therapies.

Reducing Susceptibility Distortion Related Image Blurring in Diffusion MRI EPI Data

Diffusion magnetic resonance imaging (MRI) is an increasingly popular technique in basic and clinical neuroscience. One promising application is to combine diffusion MRI with myelin maps from complementary MRI techniques such as multi-parameter mapping (MPM) to produce g-ratio maps that represent the relative myelination of axons and predict their conduction velocity. Statistical Parametric Mapping (SPM) can process both diffusion data and MPMs, making SPM the only widely accessible software that contains all the processing steps required to perform group analyses of g-ratio data in a common space. However, limitations have been identified in its method for reducing susceptibility-related distortion in diffusion data. More generally, susceptibility-related image distortion is often corrected by combining reverse phase-encoded images (blip-up and blip-down) using the arithmetic mean (AM), however, this can lead to blurred images. In this study we sought to (1) improve the susceptibility-related distortion correction for diffusion MRI data in SPM; (2) deploy an alternative approach to the AM to reduce image blurring in diffusion MRI data when combining blip-up and blip-down EPI data after susceptibility-related distortion correction; and (3) assess the benefits of these changes for g-ratio mapping. We found that the new processing pipeline, called consecutive Hyperelastic Susceptibility Artefact Correction (HySCO) improved distortion correction when compared to the standard approach in the ACID toolbox for SPM. Moreover, using a weighted average (WA) method to combine the distortion corrected data from each phase-encoding polarity achieved greater overlap of diffusion and more anatomically faithful structural white matter probability maps derived from minimally distorted multi-parameter maps as compared to the AM. Third, we showed that the consecutive HySCO WA performed better than the AM method when combined with multi-parameter maps to perform g-ratio mapping. These improvements mean that researchers can conveniently access a wide range of diffusion-related analysis methods within one framework because they are now available within the open-source ACID toolbox as part of SPM, which can be easily combined with other SPM toolboxes, such as the hMRI toolbox, to facilitate computation of myelin biomarkers that are necessary for g-ratio mapping.

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.

Decoding the microstructural properties of white matter using realistic models

Multi-echo gradient echo (ME-GRE) magnetic resonance signal evolution in white matter has a strong dependence on the orientation of myelinated axons with respect to the main static field. Although analytical solutions have been able to predict some of the white matter (WM) signal behaviour of the hollow cylinder model, it has been shown that realistic models of WM offer a better description of the signal behaviour observed. In this work, we present a pipeline to (i) generate realistic 2D WM models with their microstructure based on real axon morphology with adjustable fiber volume fraction (FVF) and g-ratio. We (ii) simulate their interaction with the static magnetic field to be able to simulate their MR signal. For the first time, we (iii) demonstrate that realistic 2D WM models can be used to simulate a MR signal that provides a good approximation of the signal obtained from a real 3D WM model derived from electron microscopy. We then (iv) demonstrate in silico that 2D WM models can be used to predict microstructural parameters in a robust way if ME-GRE multi-orientation data is available and the main fiber orientation in each pixel is known using DTI. A deep learning network was trained and characterized in its ability to recover the desired microstructural parameters such as FVF, g-ratio, free and bound water transverse relaxation and magnetic susceptibility. Finally, the network was trained to recover these micro-structural parameters from an ex vivo dataset acquired in 9 orientations with respect to the magnetic field and 12 echo times. We demonstrate that this is an overdetermined problem and that as few as 3 orientations can already provide comparable results for some of the decoded metrics.

A comprehensive library of human transcription factors for cell fate engineering

Human pluripotent stem cells (hPSCs) offer an unprecedented opportunity to model diverse cell types and tissues. To enable systematic exploration of the programming landscape mediated by transcription factors (TFs), we present the Human TFome, a comprehensive library containing 1,564 TF genes and 1,732 TF splice isoforms. By screening the library in three hPSC lines, we discovered 290 TFs, including 241 that were previously unreported, that induce differentiation in 4 days without alteration of external soluble or biomechanical cues. We used four of the hits to program hPSCs into neurons, fibroblasts, oligodendrocytes and vascular endothelial-like cells that have molecular and functional similarity to primary cells. Our cell-autonomous approach enabled parallel programming of hPSCs into multiple cell types simultaneously. We also demonstrated orthogonal programming by including oligodendrocyte-inducible hPSCs with unmodified hPSCs to generate cerebral organoids, which expedited in situ myelination. Large-scale combinatorial screening of the Human TFome will complement other strategies for cell engineering based on developmental biology and computational systems biology.

Image-based biophysical modeling predicts cortical potentials evoked with subthalamic deep brain stimulation

BACKGROUND

Subthalamic deep brain stimulation (DBS) is an effective surgical treatment for Parkinson's disease and continues to advance technologically with an enormous parameter space. As such, in-silico DBS modeling systems have become common tools for research and development, but their underlying methods have yet to be standardized and validated.

OBJECTIVE

Evaluate the accuracy of patient-specific estimates of neural pathway activations in the subthalamic region against intracranial, cortical evoked potential (EP) recordings.

METHODS

Pathway activations were modeled in eleven patients using the latest advances in connectomic modeling of subthalamic DBS, focusing on the hyperdirect pathway (HDP) and corticospinal/bulbar tract (CSBT) for their relevance in human research studies. Correlations between pathway activations and respective EP amplitudes were quantified.

RESULTS

Good model performance required accurate lead localization and image fusions, as well as appropriate selection of fiber diameter in the biophysical model. While optimal model parameters varied across patients, good performance could be achieved using a global set of parameters that explained 60% and 73% of electrophysiologic activations of CSBT and HDP, respectively. Moreover, restricted models fit to only EP amplitudes of eight standard (monopolar and bipolar) electrode configurations were able to extrapolate variation in EP amplitudes across other directional electrode configurations and stimulation parameters, with no significant reduction in model performance across the cohort.

CONCLUSIONS

Our findings demonstrate that connectomic models of DBS with sufficient anatomical and electrical details can predict recruitment dynamics of white matter. These results will help to define connectomic modeling standards for preoperative surgical targeting and postoperative patient programming applications.

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.

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, Trembler^J 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.

Parkinson's disease: deep learning with a parameter-weighted structural connectome matrix for diagnosis and neural circuit disorder investigation

Purpose

To investigate whether Parkinson’s disease (PD) can be differentiated from healthy controls and to identify neural circuit disorders in PD by applying a deep learning technique to parameter-weighted and number of streamlines (NOS)–based structural connectome matrices calculated from diffusion-weighted MRI.

Methods

In this prospective study, 115 PD patients and 115 healthy controls were enrolled. NOS-based and parameter-weighted connectome matrices were calculated from MRI images obtained with a 3-T MRI unit. With 5-fold cross-validation, diagnostic performance of convolutional neural network (CNN) models using those connectome matrices in differentiating patients with PD from healthy controls was evaluated. To identify the important brain connections for diagnosing PD, gradient-weighted class activation mapping (Grad-CAM) was applied to the trained CNN models.

Results

CNN models based on some parameter-weighted structural matrices (diffusion kurtosis imaging (DKI)–weighted, neurite orientation dispersion and density imaging (NODDI)–weighted, and g-ratio-weighted connectome matrices) showed moderate performance (areas under the receiver operating characteristic curve (AUCs) = 0.895, 0.801, and 0.836, respectively) in discriminating PD patients from healthy controls. The DKI-weighted connectome matrix performed significantly better than the conventional NOS-based matrix (AUC = 0.761) (DeLong’s test, p < 0.0001). Alterations of neural connections between the basal ganglia and cerebellum were indicated by applying Grad-CAM to the NODDI- and g-ratio-weighted matrices.

Conclusion

Patients with PD can be differentiated from healthy controls by applying the deep learning technique to the parameter-weighted connectome matrices, and neural circuit disorders including those between the basal ganglia on one side and the cerebellum on the contralateral side were visualized.

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.

The relationship between hippocampal-dependent task performance and hippocampal grey matter myelination and iron content

Individual differences in scene imagination, autobiographical memory recall, future thinking and spatial navigation have long been linked with hippocampal structure in healthy people, although evidence for such relationships is, in fact, mixed. Extant studies have predominantly concentrated on hippocampal volume. However, it is now possible to use quantitative neuroimaging techniques to model different properties of tissue microstructure in vivo such as myelination and iron. Previous work has linked such measures with cognitive task performance, particularly in older adults. Here we investigated whether performance on scene imagination, autobiographical memory, future thinking and spatial navigation tasks was associated with hippocampal grey matter myelination or iron content in young, healthy adult participants. Magnetic resonance imaging data were collected using a multi-parameter mapping protocol (0.8 mm isotropic voxels) from a large sample of 217 people with widely-varying cognitive task scores. We found little evidence that hippocampal grey matter myelination or iron content were related to task performance. This was the case using different analysis methods (voxel-based quantification, partial correlations), when whole brain, hippocampal regions of interest, and posterior:anterior hippocampal ratios were examined, and across different participant sub-groups (divided by gender and task performance). Variations in hippocampal grey matter myelin and iron levels may not, therefore, help to explain individual differences in performance on hippocampal-dependent tasks, at least in young, healthy individuals.

Myelin and Axonal Damage in Normal-Appearing White Matter in Patients with Moyamoya Disease

BACKGROUND AND PURPOSE

Although chronic ischemia is known to induce myelin and axonal damage in animal models, knowledge regarding patients with Moyamoya disease is limited. We aimed to investigate the presence of myelin and axonal damage in Moyamoya disease and their relationship with cognitive performance.

MATERIALS AND METHODS

Eighteen patients with Moyamoya disease (16-55 years of age) and 18 age- and sex-matched healthy controls were evaluated with myelin-sensitive MR imaging based on magnetization transfer saturation imaging and 2-shell diffusion MR imaging. The myelin volume fraction, which reflects the amount of myelin sheath; the g-ratio, which represents the ratio of the inner (axon) to the outer (axon plus myelin) diameter of the fiber; and the axon volume fraction, which reflects axonal components, were calculated and compared between the patients and controls. In the patients with Moyamoya disease, the relationship between these parameters and cognitive task-measuring performance speed was also evaluated.

RESULTS

Compared with the healthy controls, the patients with Moyamoya disease showed a significant decrease in the myelin and axon volume fractions (P < .05) in many WM regions, while the increases in the g-ratio values were not statistically significant. Correlations with cognitive performance were most frequently observed with the axon volume fraction (r = 0.52-0.54; P < .03 in the right middle and posterior cerebral artery areas) and were the strongest with the g-ratio values in the right posterior cerebral artery region (r = 0.64; P = .004).

CONCLUSIONS

Myelin-sensitive MR imaging and diffusion MR imaging revealed that myelin and axonal damage exist in patients with Moyamoya disease. The relationship with cognitive performance might be stronger with axonal damage than with myelin damage.

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

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

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.

Action potential propagation and synchronisation in myelinated axons

With the advent of advanced MRI techniques it has become possible to study axonal white matter non-invasively and in great detail. Measuring the various parameters of the long-range connections of the brain opens up the possibility to build and refine detailed models of large-scale neuronal activity. One particular challenge is to find a mathematical description of action potential propagation that is sufficiently simple, yet still biologically plausible to model signal transmission across entire axonal fibre bundles. We develop a mathematical framework in which we replace the Hodgkin-Huxley dynamics by a spike-diffuse-spike model with passive sub-threshold dynamics and explicit, threshold-activated ion channel currents. This allows us to study in detail the influence of the various model parameters on the action potential velocity and on the entrainment of action potentials between ephaptically coupled fibres without having to recur to numerical simulations. Specifically, we recover known results regarding the influence of axon diameter, node of Ranvier length and internode length on the velocity of action potentials. Additionally, we find that the velocity depends more strongly on the thickness of the myelin sheath than was suggested by previous theoretical studies. We further explain the slowing down and synchronisation of action potentials in ephaptically coupled fibres by their dynamic interaction. In summary, this study presents a solution to incorporate detailed axonal parameters into a whole-brain modelling framework.

With more and more data becoming available on white-matter tracts, the need arises to develop modelling frameworks that incorporate these data at the whole-brain level. This requires the development of efficient mathematical schemes to study parameter dependencies that can then be matched with data, in particular the speed of action potentials that cause delays between brain regions. Here, we develop a method that describes the formation of action potentials by threshold activated currents, often referred to as spike-diffuse-spike modelling. A particular focus of our study is the dependence of the speed of action potentials on structural parameters. We find that the diameter of axons and the thickness of the myelin sheath have a strong influence on the speed, whereas the length of myelinated segments and node of Ranvier length have a lesser effect. In addition to examining single axons, we demonstrate that action potentials between nearby axons can synchronise and slow down their propagation speed.

MR g-ratio-weighted connectome analysis in patients with multiple sclerosis

Multiple sclerosis (MS) is a brain network disconnection syndrome. Although the brain network topology in MS has been evaluated using diffusion MRI tractography, the mechanism underlying disconnection in the disorder remains unclear. In this study, we evaluated the brain network topology in MS using connectomes with connectivity strengths based on the ratio of the inner to outer myelinated axon diameter (i.e., g-ratio), thereby providing enhanced sensitivity to demyelination compared with the conventional measures of connectivity. We mapped g-ratio-based connectomes in 14 patients with MS and compared them with those of 14 age- and sex-matched healthy controls. For comparison, probabilistic tractography was also used to map connectomes based on the number of streamlines (NOS). We found that g-ratio- and NOS-based connectomes comprised significant connectivity reductions in patients with MS, predominantly in the motor, somatosensory, visual, and limbic regions. However, only the g-ratio-based connectome enabled detection of significant increases in nodal strength in patients with MS. Finally, we found that the g-ratio-weighted nodal strength in motor, visual, and limbic regions significantly correlated with inter-individual variation in measures of disease severity. The g-ratio-based connectome can serve as a sensitive biomarker for diagnosing and monitoring disease progression.

hMRI - A toolbox for quantitative MRI in neuroscience and clinical research

Neuroscience and clinical researchers are increasingly interested in quantitative magnetic resonance imaging (qMRI) due to its sensitivity to micro-structural properties of brain tissue such as axon, myelin, iron and water concentration. We introduce the hMRI-toolbox, an open-source, easy-to-use tool available on GitHub, for qMRI data handling and processing, presented together with a tutorial and example dataset. This toolbox allows the estimation of high-quality multi-parameter qMRI maps (longitudinal and effective transverse relaxation rates R1 and R2⋆, proton density PD and magnetisation transfer MT saturation) that can be used for quantitative parameter analysis and accurate delineation of subcortical brain structures. The qMRI maps generated by the toolbox are key input parameters for biophysical models designed to estimate tissue microstructure properties such as the MR g-ratio and to derive standard and novel MRI biomarkers. Thus, the current version of the toolbox is a first step towards in vivo histology using MRI (hMRI) and is being extended further in this direction. Embedded in the Statistical Parametric Mapping (SPM) framework, it benefits from the extensive range of established SPM tools for high-accuracy spatial registration and statistical inferences and can be readily combined with existing SPM toolboxes for estimating diffusion MRI parameter maps. From a user's perspective, the hMRI-toolbox is an efficient, robust and simple framework for investigating qMRI data in neuroscience and clinical research.

Evolution of white matter tract microstructure across the life span

The human brain undergoes dramatic structural change over the life span. In a large imaging cohort of 801 individuals aged 7-84 years, we applied quantitative relaxometry and diffusion microstructure imaging in combination with diffusion tractography to investigate tissue property dynamics across the human life span. Significant nonlinear aging effects were consistently observed across tracts and tissue measures. The age at which white matter (WM) fascicles attain peak maturation varies substantially across tissue measurements and tracts. These observations of heterochronicity and spatial heterogeneity of tract maturation highlight the importance of using multiple tissue measurements to investigate each region of the WM. Our data further provide additional quantitative evidence in support of the last-in-first-out retrogenesis hypothesis of aging, demonstrating a strong correlational relationship between peak maturational timing and the extent of quadratic measurement differences across the life span for the most myelin sensitive measures. These findings present an important baseline from which to assess divergence from normative aging trends in developmental and degenerative disorders, and to further investigate the mechanisms connecting WM microstructure to cognition.

Diffusivity and quantitative T1 profile of human visual white matter tracts after retinal ganglion cell damage

In patients with retinal ganglion cell diseases, recent diffusion tensor imaging (DTI) studies have revealed structural abnormalities in visual white matter tracts such as the optic tract, and optic radiation. However, the microstructural origin of these diffusivity changes is unknown as DTI metrics involve multiple biological factors and do not correlate directly with specific microstructural properties. In contrast, recent quantitative T1 (qT1) mapping methods provide tissue property measurements relatively specific to myelin volume fractions in white matter. This study aims to improve our understanding of microstructural changes in visual white matter tracts following retinal ganglion cell damage in Leber's hereditary optic neuropathy (LHON) patients by combining DTI and qT1 measurements. We collected these measurements from seven LHON patients and twenty age-matched control subjects. For all individuals, we identified the optic tract and the optic radiation using probabilistic tractography, and evaluated diffusivity and qT1 profiles along them. Both diffusivity and qT1 measurements in the optic tract differed significantly between LHON patients and controls. In the optic radiation, these changes were observed in diffusivity but were not evident in qT1 measurements. This suggests that myelin loss may not explain trans-synaptic diffusivity changes in the optic radiation as a consequence of retinal ganglion cell disease.

•

Retinal ganglion cell damage affects diffusivity and T1 along visual pathways.

•

DTI metric identified white matter change in both optic tract and optic radiation.

•

T1 measurement in optic radiation did not exhibit abnormality, unlike DTI metric.

•

Myelin loss may not be a major cause of diffusivity change along optic radiation.

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.

Age-related changes in myelin of axons of the corpus callosum and cognitive decline in common marmosets

Executive control is a higher-level cognitive function that involves a range of different processes that are involved in the planning, coordination, execution, and inhibition of responses. Many of the processes associated with executive control, such as response inhibition and mental flexibility, decline with age. Degeneration of white matter architecture is considered to be the one of the key factors underlying cognitive decline associated with aging. Here we investigated how white matter changes of the corpus callosum were related to cognitive aging in common marmosets (Callithrix jacchus). We hypothesized that reduction in myelin thickness, myelin density, and myelin fraction of axonal fibers in the corpus callosum would be associated with performance on a task of executive function in a small sample of geriatric marmosets (n = 4) and young adult marmosets (n = 2). Our results indicated declines in myelin thickness, density, and myelin fraction with age. Considerable variability was detected on these characteristics of myelin and cognitive performance assessed via the detoured reach task. Age-related changes in myelin in Region II of the corpus callosum were predictive of cognitive performance on the detoured reach task. Thus the detoured reach task appears to also measure aspects of corticostriatal function in addition to prefrontal cortical function.

The development of brain white matter microstructure

Throughout infancy, childhood, and adolescence, our brains undergo remarkable changes. Processes including myelination and synaptogenesis occur rapidly across the first 2-3 years of life, and ongoing brain remodeling continues into young adulthood. Studies have sought to characterize the patterns of structural brain development, and early studies predominately relied upon gross anatomical measures of brain structure, morphology, and organization. MRI offers the ability to characterize and quantify a range of microstructural aspects of brain tissue that may be more closely related to fundamental neurodevelopmental processes. Techniques such as diffusion, magnetization transfer, relaxometry, and myelin water imaging provide insight into changing cyto- and myeloarchitecture, neuronal density, and structural connectivity. In this review, we focus on the growing body of literature exploiting these MRI techniques to better understand the microstructural changes that occur in brain white matter during maturation. Our review focuses on studies of normative brain development from birth to early adulthood (∼25 years), and places particular emphasis on longitudinal studies and newer techniques that are being used to study microstructural white matter development. All imaging methods demonstrate consistent, rapid microstructural white matter development over the first 3 years of life, suggesting increased myelination and axonal packing. Diffusion studies clearly demonstrate continued white matter maturation during later childhood and adolescence, though the lack of consistent findings in other modalities suggests changes may be mainly due to axonal packing. An emerging literature details differential microstructural development in boys and girls, and connects developmental trajectories to cognitive abilities, behaviour, and/or environmental factors, though the nature of these relationships remains unclear. Future research will need to focus on newer imaging techniques and longitudinal studies to provide more detailed information about microstructural white matter development, particularly in the childhood years.

Developing 3D microscopy with CLARITY on human brain tissue: Towards a tool for informing and validating MRI-based histology

Recent breakthroughs in magnetic resonance imaging (MRI) enabled quantitative relaxometry and diffusion-weighted imaging with sub-millimeter resolution. Combined with biophysical models of MR contrast the emerging methods promise in vivo mapping of cyto- and myelo-architectonics, i.e., in vivo histology using MRI (hMRI) in humans. The hMRI methods require histological reference data for model building and validation. This is currently provided by MRI on post mortem human brain tissue in combination with classical histology on sections. However, this well established approach is limited to qualitative 2D information, while a systematic validation of hMRI requires quantitative 3D information on macroscopic voxels. We present a promising histological method based on optical 3D imaging combined with a tissue clearing method, Clear Lipid-exchanged Acrylamide-hybridized Rigid Imaging compatible Tissue hYdrogel (CLARITY), adapted for hMRI validation. Adapting CLARITY to the needs of hMRI is challenging due to poor antibody penetration into large sample volumes and high opacity of aged post mortem human brain tissue. In a pilot experiment we achieved transparency of up to 8 mm-thick and immunohistochemical staining of up to 5 mm-thick post mortem brain tissue by a combination of active and passive clearing, prolonged clearing and staining times. We combined 3D optical imaging of the cleared samples with tailored image processing methods. We demonstrated the feasibility for quantification of neuron density, fiber orientation distribution and cell type classification within a volume with size similar to a typical MRI voxel. The presented combination of MRI, 3D optical microscopy and image processing is a promising tool for validation of MRI-based microstructure estimates.

Linking structural and effective brain connectivity: structurally informed Parametric Empirical Bayes (si-PEB)

Despite the potential for better understanding functional neuroanatomy, the complex relationship between neuroimaging measures of brain structure and function has confounded integrative, multimodal analyses of brain connectivity. This is particularly true for task-related effective connectivity, which describes the causal influences between neuronal populations. Here, we assess whether measures of structural connectivity may usefully inform estimates of effective connectivity in larger scale brain networks. To this end, we introduce an integrative approach, capitalising on two recent statistical advances: Parametric Empirical Bayes, which provides group-level estimates of effective connectivity, and Bayesian model reduction, which enables rapid comparison of competing models. Crucially, we show that structural priors derived from high angular resolution diffusion imaging on a dynamic causal model of a 12-region network-based on functional MRI data from the same subjects-substantially improve model evidence (posterior probability 1.00). This provides definitive evidence that structural and effective connectivity depend upon each other in mediating distributed, large-scale interactions in the brain. Furthermore, this work offers novel perspectives for understanding normal brain architecture and its disintegration in clinical conditions.

White Matter Tissue Quantification at Low b-Values Within Constrained Spherical Deconvolution Framework

In the last decades, a number of Diffusion Weighted Imaging (DWI) based techniques have been developed to study non-invasively human brain tissues, especially white matter (WM). In this context, Constrained Spherical Deconvolution (CSD) is recognized as being able to accurately characterize water molecules displacement, as they emerge from the observation of MR diffusion weighted (MR-DW) images. CSD is suggested to be applied on MR-DW datasets consisting of b-values around 3,000 s/mm² and at least 45 unique diffusion weighting directions. Below such technical requirements, Diffusion Tensor Imaging (DT) remains the most widely accepted model. Unlike CSD, DTI is unable to resolve complex fiber geometries within the brain, thus affecting related tissues quantification. In addition, thanks to CSD, an index called Apparent Fiber Density (AFD) can be measured to estimate intra-axonal volume fraction within WM. In standard clinical settings, diffusion based acquisitions are well below such technical requirements. Therefore, in this study we wanted to extensively compare CSD and DTI model outcomes on really low demanding MR-DW datasets, i.e., consisting of a single shell (b-value = 1,000 s/mm²) and only 30 unique diffusion encoding directions. To this end, we performed deterministic and probabilistic tractographic reconstruction of two major WM pathways, namely the Corticospinal Tract and the Arcuate Fasciculus. We estimated and analyzed tensor based features as well as, for the first time, AFD interpretability in our data. By performing multivariate statistics and tract-based ROI analysis, we demonstrate that WM quantification is affected by both the diffusion model and threshold applied to noisy tractographic maps. Consistently with existing literature, we showed that CSD outperforms DTI even in our scenario. Most importantly, for the first time we address the problem of accuracy and interpretation of AFD in a low-demanding DW setup, and show that it is still a biological meaningful measure for the analysis of intra-axonal volume even in clinical settings.

Myelin Measurement: Comparison Between Simultaneous Tissue Relaxometry, Magnetization Transfer Saturation Index, and T1w/T2w Ratio Methods

Magnetization transfer (MT) imaging has been widely used for estimating myelin content in the brain. Recently, two other approaches, namely simultaneous tissue relaxometry of R1 and R2 relaxation rates and proton density (SyMRI) and the ratio of T1-weighted to T2-weighted images (T1w/T2w ratio), were also proposed as methods for measuring myelin. SyMRI and MT imaging have been reported to correlate well with actual myelin by histology. However, for T1w/T2w ratio, such evidence is limited. In 20 healthy adults, we examined the correlation between these three methods, using MT saturation index (MTsat) for MT imaging. After calibration, white matter (WM) to gray matter (GM) contrast was the highest for SyMRI among these three metrics. Even though SyMRI and MTsat showed strong correlation in the WM (r = 0.72), only weak correlation was found between T1w/T2w and SyMRI (r = 0.45) or MTsat (r = 0.38) (correlation coefficients significantly different from each other, with p values < 0.001). In subcortical and cortical GM, these measurements showed moderate to strong correlations to each other (r = 0.54 to 0.78). In conclusion, the high correlation between SyMRI and MTsat indicates that both methods are similarly suited to measure myelin in the WM, whereas T1w/T2w ratio may be less optimal.

Application of Quantitative Microstructural MR Imaging with Atlas-based Analysis for the Spinal Cord in Cervical Spondylotic Myelopathy

Mapping of MR fiber g-ratio, which is the ratio of the diameter of the axon to the diameter of the neuronal fiber, is introduced in this article. We investigated the MR fiber g-ratio, the axon volume fraction (AVF) and the myelin volume fraction (MVF) to evaluate microstructural changes in the spinal cord in patients with cervical spondylotic myelopathy (CSM) in vivo, using atlas-based analysis. We used diffusion MRI data acquired with a new simultaneous multi-slice accelerated readout-segmented echo planar imaging sequence for diffusion analysis for AVF calculation and magnetization transfer saturation imaging for MVF calculation. The AVFs of fasciculus gracilis in the affected side spinal cord, fasciculus cuneatus and lateral corticospinal tracts (LSCT) in the affected and unaffected side spinal cord were significantly lower (P = 0.019, 0.001, 0019, 0.000, and 0.002, respectively) than those of normal controls. No difference was found in the MVFs. The fiber g-ratio of LSCT was significantly lower (P = 0.040) in the affected side spinal cords than in the normal controls. The pathological microstructural changes in the spinal cord in patients with CSM, presumably partial axonal degenerations with preserved myelin. This technique has the potential to be a clinical biomarker in patients with CSM in vivo.

Computational neuroanatomy of human stratum proprium of interparietal sulcus

Recent advances in diffusion-weighted MRI (dMRI) and tractography have enabled identification of major long-range white matter tracts in the human brain. Yet, our understanding of shorter tracts, such as those within the parietal lobe, remains limited. Over a century ago, a tract connecting the superior and inferior parts of the parietal cortex was identified in a post-mortem study: stratum proprium of interparietal sulcus (SIPS; Sachs, Das hemisphärenmark des menschlichen grosshirns. Verlag von georg thieme, Leipzig, 1892). The tract has since been replicated in another fibre dissection study (Vergani et al., Cortex 56:145-156, 2014), however, it has not been fully investigated in the living human brain and its precise anatomical properties are yet to be described. We used dMRI and tractography to identify and characterise SIPS in vivo, and explored its spatial proximity to the cortical areas associated with optic-flow processing using fMRI. SIPS was identified bilaterally in all subjects, and its anatomical position and trajectory are consistent with previous post-mortem studies. Subsequent evaluation of the tractography results using the linear fascicle evaluation and virtual lesion analysis yielded strong statistical evidence for SIPS. We also found that the SIPS endpoints are adjacent to the optic-flow selective areas. In sum, we show that SIPS is a short-range tract connecting the superior and inferior parts of the parietal cortex, wrapping around the intraparietal sulcus, and that it may be a crucial anatomy underlying optic-flow processing. In vivo identification and characterisation of SIPS will facilitate further research on SIPS in relation to cortical functions, their development, and diseases that affect them.

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.

Quantitative MRI provides markers of intra-, inter-regional, and age-related differences in young adult cortical microstructure

Measuring the structural composition of the cortex is critical to understanding typical development, yet few investigations in humans have charted markers in vivo that are sensitive to tissue microstructural attributes. Here, we used a well-validated quantitative MR protocol to measure four parameters (R₁, MT, R₂*, PD*) that differ in their sensitivity to facets of the tissue microstructural environment (R₁, MT: myelin, macromolecular content; R₂*: myelin, paramagnetic ions, i.e., iron; PD*: free water content). Mapping these parameters across cortical regions in a young adult cohort (18–39 years, N = 93) revealed expected patterns of increased macromolecular content as well as reduced tissue water content in primary and primary adjacent cortical regions. Mapping across cortical depth within regions showed decreased expression of myelin and related processes – but increased tissue water content – when progressing from the grey/white to the grey/pial boundary, in all regions. Charting developmental change in cortical microstructure cross-sectionally, we found that parameters with sensitivity to tissue myelin (R₁ & MT) showed linear increases with age across frontal and parietal cortex (change 0.5–1.0% per year). Overlap of robust age effects for both parameters emerged in left inferior frontal, right parietal and bilateral pre-central regions. Our findings afford an improved understanding of ontogeny in early adulthood and offer normative quantitative MR data for inter- and intra-cortical composition, which may be used as benchmarks in further studies.

•

We mapped multi-parameter maps (MPMs) across and within cortical regions.

•

We charted age effects (ages 18–39) on myelin and related processes.

•

MPMs sensitive to myelin (R₁, MT) showed elevated values in primary areas over most cortical depths.

•

R₂* map foci tended to overlap MPMs sensitive to myelin (R₁, MT).

•

R₁ and MT increased with age (0.5–1.0% per year) at mid-depth in frontal and parietal cortex.

Modelling and interpretation of magnetization transfer imaging in the brain

Magnetization transfer contrast has yielded insight into brain tissue microstructure changes across the lifespan and in a range of disorders. This progress has been aided by the development of quantitative magnetization transfer imaging techniques able to extract intrinsic properties of the tissue that are independent of the specifics of the data acquisition. While the tissue properties extracted by these techniques do not map directly onto specific cellular structures or pathological processes, a growing body of work from animal models and histopathological correlations aids the in vivo interpretation of magnetization transfer properties of tissue. This review examines the biophysical models that have been developed to describe magnetization transfer contrast in tissue as well as the experimental evidence for the biological interpretation of magnetization transfer data in health and disease.

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

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