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

Advanced Quantitative MRI Unveils Microstructural Thalamic Changes Reflecting Disease Progression in Multiple Sclerosis

BACKGROUND AND OBJECTIVES

In patients with multiple sclerosis (PwMS), thalamic atrophy occurs during the disease course. However, there is little understanding of the mechanisms leading to volume loss and of the relationship between microstructural thalamic pathology and disease progression. This cross-sectional and longitudinal study aimed to comprehensively characterize in vivo pathologic changes within thalamic microstructure in PwMS using advanced multiparametric quantitative MRI (qMRI).

METHODS

Thalamic microstructural integrity was evaluated using quantitative T1, magnetization transfer saturation, multishell diffusion, and quantitative susceptibility mapping (QSM) in 183 PwMS and 105 healthy controls (HCs). The same qMRI protocol was available for 127 PwMS and 73 HCs after a 2-year follow-up period. Inclusion criteria for PwMS encompassed either an active relapsing-remitting MS (RRMS) or inactive progressive MS (PMS) disease course. Thalamic alterations were compared between PwMS and HCs and among disease phenotypes. In addition, the study investigated the relationship between thalamic damage and clinical and conventional MRI measures of disease severity.

RESULTS

Compared with HCs, PwMS exhibited substantial thalamic alterations, indicative of microstructural and macrostructural damage, demyelination, and disruption in iron homeostasis. These alterations extended beyond focal thalamic lesions, affecting normal-appearing thalamic tissue diffusely. Over the follow-up period, PwMS displayed an accelerated decrease in myelin volume fraction [mean difference in annualized percentage change (MD-ApC) = -1.50; p = 0.041] and increase in quantitative T1 (MD-ApC = 0.92; p < 0.0001) values, indicating heightened demyelinating and neurodegenerative processes. The observed differences between PwMS and HCs were substantially driven by the subgroup with PMS, wherein thalamic degeneration was significantly accelerated, even in comparison with patients with RRMS. Thalamic qMRI alterations showed extensive correlations with conventional MRI, clinical, and cognitive disease burden measures. Disability progression over follow-up was associated with accelerated thalamic degeneration, as reflected by enhanced diffusion (β = -0.067; p = 0.039) and QSM (β = -0.077; p = 0.027) changes. Thalamic qMRI metrics emerged as significant predictors of neurologic and cognitive disability even when accounting for other established markers including white matter lesion load and brain and thalamic atrophy.

DISCUSSION

These findings offer deeper insights into thalamic pathology in PwMS, emphasizing the clinical relevance of thalamic damage and its link to disease progression. Advanced qMRI biomarkers show promising potential in guiding interventions aimed at mitigating thalamic neurodegenerative processes.

Advanced MRI Measures of Myelin and Axon Volume Identify Repair in Multiple Sclerosis

OBJECTIVE

Pathological studies suggest that multiple sclerosis (MS) lesions endure multiple waves of damage and repair; however, the dynamics and characteristics of these processes are poorly understood in patients living with MS.

METHODS

We studied 128 MS patients (75 relapsing-remitting, 53 progressive) and 72 healthy controls who underwent advanced magnetic resonance imaging and clinical examination at baseline and 2 years later. Magnetization transfer saturation and multi-shell diffusion imaging were used to quantify longitudinal changes in myelin and axon volumes within MS lesions. Lesions were grouped into 4 classes (repair, damage, mixed repair damage, and stable). The frequency of each class was correlated to clinical measures, demographic characteristics, and levels of serum neurofilament light chain (sNfL).

RESULTS

Stable lesions were the most frequent (n = 2,276; 44%), followed by lesions with patterns of "repair" (n = 1,352; 26.2%) and damage (n = 1,214; 23.5%). The frequency of "repair" lesion was negatively associated with disability (β = -0.04; p < 0.001) and sNfL (β = -0.16; p < 0.001) at follow-up. The frequency of the "damage" class was higher in progressive than relapsing-remitting patients (p < 0.05) and was related to disability (baseline: β = -0.078; follow-up: β = -0.076; p < 0.001) and age (baseline: β = -0.078; p < 0.001). Stable lesions were more frequent in relapsing-remitting than in progressive patients (p < 0.05), and in younger patients versus older (β = -0.07; p < 0.001) at baseline. Further, "mixed" lesions were most frequent in older patients (β = 0.004; p < 0.001) at baseline.

INTERPRETATION

These findings show that repair and damage processes within MS lesions occur across the entire disease spectrum and that their frequency correlates with patients disability, age, disease duration, and extent of neuroaxonal damage. ANN NEUROL 2024.

A multi-modal, asymmetric, weighted, and signed description of anatomical connectivity

The macroscale connectome is the network of physical, white-matter tracts between brain areas. The connections are generally weighted and their values interpreted as measures of communication efficacy. In most applications, weights are either assigned based on imaging features-e.g. diffusion parameters-or inferred using statistical models. In reality, the ground-truth weights are unknown, motivating the exploration of alternative edge weighting schemes. Here, we explore a multi-modal, regression-based model that endows reconstructed fiber tracts with directed and signed weights. We find that the model fits observed data well, outperforming a suite of null models. The estimated weights are subject-specific and highly reliable, even when fit using relatively few training samples, and the networks maintain a number of desirable features. In summary, we offer a simple framework for weighting connectome data, demonstrating both its ease of implementation while benchmarking its utility for typical connectome analyses, including graph theoretic modeling and brain-behavior associations.

Pathologic classification of a late-onset peripheral neuropathy in a spontaneous Labrador retriever dog model

Late-onset peripheral neuropathy (LPN) is a heritable canine neuropathy commonly found in Labrador retrievers and is characterized by laryngeal paralysis and pelvic limb paresis. Our objective was to establish canine LPN as a model for human hereditary peripheral neuropathy by classifying it as either an axonopathy or myelinopathy and evaluating length-dependent degeneration. We conducted a motor nerve conduction study of the sciatic and ulnar nerves, electromyography (EMG) of appendicular and epaxial musculature, and histologic analysis of sciatic and recurrent laryngeal nerves in LPN-affected and control dogs. LPN-affected dogs exhibited significant decreases in compound muscle action potential (CMAP) amplitude, CMAP area, and pelvic limb latencies. However, no differences were found in motor nerve conduction velocity, residual latencies, or CMAP duration. Distal limb musculature showed greater EMG changes in LPN-affected dogs. Histologically, LPN-affected dogs exhibited a reduction in the number of large-diameter axons, especially in distal nerve regions. In conclusion, LPN in Labrador retrievers is a common, spontaneous, length-dependent peripheral axonopathy that is a novel animal model of age-related peripheral neuropathy that could be used for fundamental research and clinical trials.

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.

The human brain connectome weighted by the myelin content and total intra-axonal cross-sectional area of white matter tracts

A central goal in neuroscience is the development of a comprehensive mapping between structural and functional brain features, which facilitates mechanistic interpretation of brain function. However, the interpretability of structure-function brain models remains limited by a lack of biological detail. Here, we characterize human structural brain networks weighted by multiple white matter microstructural features including total intra-axonal cross-sectional area and myelin content. We report edge-weight-dependent spatial distributions, variance, small-worldness, rich club, hubs, as well as relationships with function, edge length, and myelin. Contrasting networks weighted by the total intra-axonal cross-sectional area and myelin content of white matter tracts, we find opposite relationships with functional connectivity, an edge-length-independent inverse relationship with each other, and the lack of a canonical rich club in myelin-weighted networks. When controlling for edge length, networks weighted by either fractional anisotropy, radial diffusivity, or neurite density show no relationship with whole-brain functional connectivity. We conclude that the co-utilization of structural networks weighted by total intra-axonal cross-sectional area and myelin content could improve our understanding of the mechanisms mediating the structure-function brain relationship.

Towards a biologically annotated brain connectome

The brain is a network of interleaved neural circuits. In modern connectomics, brain connectivity is typically encoded as a network of nodes and edges, abstracting away the rich biological detail of local neuronal populations. Yet biological annotations for network nodes - such as gene expression, cytoarchitecture, neurotransmitter receptors or intrinsic dynamics - can be readily measured and overlaid on network models. Here we review how connectomes can be represented and analysed as annotated networks. Annotated connectomes allow us to reconceptualize architectural features of networks and to relate the connection patterns of brain regions to their underlying biology. Emerging work demonstrates that annotated connectomes help to make more veridical models of brain network formation, neural dynamics and disease propagation. Finally, annotations can be used to infer entirely new inter-regional relationships and to construct new types of network that complement existing connectome representations. In summary, biologically annotated connectomes offer a compelling way to study neural wiring in concert with local biological features.

Multiparametric MRI: From Simultaneous Rapid Acquisition Methods and Analysis Techniques Using Scoring, Machine Learning, Radiomics, and Deep Learning to the Generation of Novel Metrics

ABSTRACT

With the recent advancements in rapid imaging methods, higher numbers of contrasts and quantitative parameters can be acquired in less and less time. Some acquisition models simultaneously obtain multiparametric images and quantitative maps to reduce scan times and avoid potential issues associated with the registration of different images. Multiparametric magnetic resonance imaging (MRI) has the potential to provide complementary information on a target lesion and thus overcome the limitations of individual techniques. In this review, we introduce methods to acquire multiparametric MRI data in a clinically feasible scan time with a particular focus on simultaneous acquisition techniques, and we discuss how multiparametric MRI data can be analyzed as a whole rather than each parameter separately. Such data analysis approaches include clinical scoring systems, machine learning, radiomics, and deep learning. Other techniques combine multiple images to create new quantitative maps associated with meaningful aspects of human biology. They include the magnetic resonance g-ratio, the inner to the outer diameter of a nerve fiber, and the aerobic glycolytic index, which captures the metabolic status of tumor tissues.

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.

Anti-Correlated Myelin-Sensitive MRI Levels in Humans Consistent with a Subcortical to Sensorimotor Regulatory Process-Multi-Cohort Multi-Modal Evidence

Differential axonal myelination synchronises signalling over different axon lengths. The consequences of myelination processes described at the cellular level for the regulation of myelination at the macroscopic level are unknown. We analysed multiple cohorts of myelin-sensitive brain MRI. Our aim was to (i) confirm a previous report of anti-correlation between myelination in subcortical and sensorimotor areas in healthy subjects, (ii) and thereby test our hypothesis for a regulatory interaction between them. We analysed nine image-sets across three different human cohorts using six MRI modalities. Each image-set contained healthy controls (HC) and ME/CFS subjects. Subcortical and Sensorimotor regions of interest (ROI) were optimised for the detection of anti-correlations and the same ROIs were used to test the HC in all image-sets. For each cohort, median MRI values were computed in both regions for each subject and their correlation across the cohort was computed. We confirmed negative correlations in healthy controls between subcortical and sensorimotor regions in six image-sets: three T1wSE (p = 5 × 10^-8, 5 × 10^-7, 0.002), T2wSE (p =2 × 10^-6), MTC (p = 0.01), and WM volume (p = 0.02). T1/T2 was the exception with a positive correlation (p = 0.01). This myelin regulation study is novel in several aspects: human subjects, cross-sectional design, ROI optimization, spin-echo MRI and reproducible across multiple independent image-sets. In multiple independent image-sets we confirmed an anti-correlation between subcortical and sensorimotor myelination which supports a previously unreported regulatory interaction. The subcortical region contained the brain's primary regulatory nuclei. We suggest a mechanism has evolved whereby relatively low subcortical myelination in an individual is compensated by upregulated sensorimotor myelination to maintain adequate sensorimotor performance.

Brain white matter morphological structure correlation with its optical properties estimated from optical coherence tomography (OCT) data

A pilot post-mortem study identifies a strong correlation between the attenuation coefficient estimated from the OCT data and some morphological features of the sample, namely the number of nuclei in the field of view of the histological image and the fiber structural parameter introduced in the study to quantify the difference in the myelinated fibers arrangements. The morphological features were identified from the histopathological images of the sample taken from the same locations as the OCT images and stained with the immunohistochemical (IHC) staining specific to the myelin. It was shown that the linear regression of the IHC quantitative characteristics allows adequate prediction of the attenuation coefficient of the sample. This discovery opens the opportunity for the usage of the OCT as a neuronavigation tool.

The Myelin-Weighted Connectome in Parkinson's Disease

BACKGROUND

Even though Parkinson's disease (PD) is typically viewed as largely affecting gray matter, there is growing evidence that there are also structural changes in the white matter. Traditional connectomics methods that study PD may not be specific to underlying microstructural changes, such as myelin loss.

OBJECTIVE

The primary objective of this study is to investigate the PD-induced changes in myelin content in the connections emerging from the basal ganglia and the brainstem. For the weighting of the connectome, we used the longitudinal relaxation rate as a biologically grounded myelin-sensitive metric.

METHODS

We computed the myelin-weighted connectome in 35 healthy control subjects and 81 patients with PD. We used partial least squares to highlight the differences between patients with PD and healthy control subjects. Then, a ring analysis was performed on selected brainstem and subcortical regions to evaluate each node's potential role as an epicenter for disease propagation. Then, we used behavioral partial least squares to relate the myelin alterations with clinical scores.

RESULTS

Most connections (~80%) emerging from the basal ganglia showed a reduced myelin content. The connections emerging from potential epicentral nodes (substantia nigra, nucleus basalis of Meynert, amygdala, hippocampus, and midbrain) showed significant decrease in the longitudinal relaxation rate (P < 0.05). This effect was not seen for the medulla and the pons.

CONCLUSIONS

The myelin-weighted connectome was able to identify alteration of the myelin content in PD in basal ganglia connections. This could provide a different view on the importance of myelination in neurodegeneration and disease progression. © 2021 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

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.

Repair of peripheral nerve injuries using a prevascularized cell-based tissue-engineered nerve conduit

One of the major challenges in the development of a larger and longer nerve conduit for peripheral nerve repair is the limitation in oxygen and nutrient diffusion within the tissue after transplantation preventing Schwann cell and axonal migration. This restriction is due to the slow neovascularization process of the graft starting from both nerve endings. To overcome this limitation, we propose the design of a living tissue-engineered nerve conduit made of an internal tube with a three-dimensional structure supporting axonal migration, which is inserted inside a hollow external tube that plays the role of an epineurium and is strong enough to be stitched to the severed nerve stumps. The internal tube is made of a rolled living fibroblast sheet and can be seeded with endothelial cells to promote the formation of a network containing capillary-like structures which allow rapid inosculation with the host nerve microvasculature after grafting. Human nerve conduits were grafted in immunodeficient rats to bridge a 15 mm sciatic nerve gap. Human capillaries within the pre-vascularized nerve conduit successfully connected to the host circulation 2 weeks after grafting. Twenty-two weeks after surgery, rats transplanted with the nerve conduits had a similar motor function recovery compared to the autograft group. By promoting rapid vascularization of the internal nerve tube from both ends of the nerve stumps, this endothelialized nerve conduit model displays a favorable environment to enhance axonal migration in both larger caliber and longer nerve grafts.

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.

Predicting MEG resting-state functional connectivity from microstructural information

Understanding how human brain microstructure influences functional connectivity is an important endeavor. In this work, magnetic resonance imaging data from 90 healthy participants were used to calculate structural connectivity matrices using the streamline count, fractional anisotropy, radial diffusivity, and a myelin measure (derived from multicomponent relaxometry) to assign connection strength. Unweighted binarized structural connectivity matrices were also constructed. Magnetoencephalography resting-state data from those participants were used to calculate functional connectivity matrices, via correlations of the Hilbert envelopes of beamformer time series in the delta, theta, alpha, and beta frequency bands. Nonnegative matrix factorization was performed to identify the components of the functional connectivity. Shortest path length and search-information analyses of the structural connectomes were used to predict functional connectivity patterns for each participant. The microstructure-informed algorithms predicted the components of the functional connectivity more accurately than they predicted the total functional connectivity. This provides a methodology to understand functional mechanisms better. The shortest path length algorithm exhibited the highest prediction accuracy. Of the weights of the structural connectivity matrices, the streamline count and the myelin measure gave the most accurate predictions, while the fractional anisotropy performed poorly. Overall, different structural metrics paint very different pictures of the structural connectome and its relationship to functional connectivity.

Toward nonparametric diffusion- T1 characterization of crossing fibers in the human brain

PURPOSE

To estimate T 1 for each distinct fiber population within voxels containing multiple brain tissue types.

METHODS

A diffusion- T 1 correlation experiment was carried out in an in vivo human brain using tensor-valued diffusion encoding and multiple repetition times. The acquired data were inverted using a Monte Carlo algorithm that retrieves nonparametric distributions P ( D , R 1 ) of diffusion tensors and longitudinal relaxation rates R 1 = 1 / T 1 . Orientation distribution functions (ODFs) of the highly anisotropic components of P ( D , R 1 ) were defined to visualize orientation-specific diffusion-relaxation properties. Finally, Monte Carlo density-peak clustering (MC-DPC) was performed to quantify fiber-specific features and investigate microstructural differences between white matter fiber bundles.

RESULTS

Parameter maps corresponding to P ( D , R 1 ) 's statistical descriptors were obtained, exhibiting the expected R 1 contrast between brain tissue types. Our ODFs recovered local orientations consistent with the known anatomy and indicated differences in R 1 between major crossing fiber bundles. These differences, confirmed by MC-DPC, were in qualitative agreement with previous model-based works but seem biased by the limitations of our current experimental setup.

CONCLUSIONS

Our Monte Carlo framework enables the nonparametric estimation of fiber-specific diffusion- T 1 features, thereby showing potential for characterizing developmental or pathological changes in T 1 within a given fiber bundle, and for investigating interbundle T 1 differences.

The R1-weighted connectome: complementing brain networks with a myelin-sensitive measure

Myelin plays a crucial role in how well information travels between brain regions. Complementing the structural connectome, obtained with diffusion MRI tractography, with a myelin-sensitive measure could result in a more complete model of structural brain connectivity and give better insight into white-matter myeloarchitecture. In this work we weight the connectome by the longitudinal relaxation rate (R1), a measure sensitive to myelin, and then we assess its added value by comparing it with connectomes weighted by the number of streamlines (NOS). Our analysis reveals differences between the two connectomes both in the distribution of their weights and the modular organization. Additionally, the rank-based analysis shows that R1 can be used to separate transmodal regions (responsible for higher-order functions) from unimodal regions (responsible for low-order functions). Overall, the R1-weighted connectome provides a different perspective on structural connectivity taking into account white matter myeloarchitecture.

In the present work, we show that by using a myelin-sensitive measure we can complement the diffusion MRI-based connectivity and provide a different picture of the brain organization. We show that the R1-weighted average distribution does not follow the same trend as the number of streamlines strength distribution, and the two connectomes exhibit different modular organization. We also show that unimodal cortical regions tend to be connected by more streamlines, but the connections exhibit a lower R1-weighted average, while the transmodal regions have higher R1-weighted average but fewer streamlines. This could imply that the unimodal regions require more connections with lower myelination, whereas the transmodal regions rely on connections with higher myelination.

Efficient whole-brain tract-specific T1 mapping at 3T with slice-shuffled inversion-recovery diffusion-weighted imaging

PURPOSE

Most voxels in white matter contain multiple fiber populations with different orientations and levels of myelination. Conventional T₁ mapping measures 1 T₁ value per voxel, representing a weighted average of the multiple tract T₁ times. Inversion-recovery diffusion-weighted imaging (IR-DWI) allows the T₁ times of multiple tracts in a voxel to be disentangled, but the scan time is prohibitively long. Recently, slice-shuffled IR-DWI implementations have been proposed to significantly reduce scan time. In this work, we demonstrate that we can measure tract-specific T₁ values in the whole brain using simultaneous multi-slice slice-shuffled IR-DWI at 3T.

METHODS

We perform simulations to evaluate the accuracy and precision of our crossing fiber IR-DWI signal model for various fiber parameters. The proposed sequence and signal model are tested in a phantom consisting of crossing asparagus pieces doped with gadolinium to vary T₁ , and in 2 human subjects.

RESULTS

Our simulations show that tract-specific T₁ times can be estimated within 5% of the nominal fiber T₁ values. Tract-specific T₁ values were resolved in subvoxel 2 fiber crossings in the asparagus phantom. Tract-specific T₁ times were resolved in 2 different tract crossings in the human brain where myelination differences have previously been reported; the crossing of the cingulum and genu of the corpus callosum and the crossing of the corticospinal tract and pontine fibers.

CONCLUSION

Whole-brain tract-specific T₁ mapping is feasible using slice-shuffled IR-DWI at 3T. This technique has the potential to improve the microstructural characterization of specific tracts implicated in neurodevelopment, aging, and demyelinating disorders.

Navigating the link between processing speed and network communication in the human brain

Processing speed on cognitive tasks relies upon efficient communication between widespread regions of the brain. Recently, novel methods of quantifying network communication like 'navigation efficiency' have emerged, which aim to be more biologically plausible compared to traditional shortest path length-based measures. However, it is still unclear whether there is a direct link between these communication measures and processing speed. We tested this relationship in forty-five healthy adults (27 females), where processing speed was defined as decision-making time and measured using drift rate from the hierarchical drift diffusion model. Communication measures were calculated from a graph theoretical analysis of the whole-brain structural connectome and of a task-relevant fronto-parietal structural subnetwork, using the large-scale Desikan-Killiany atlas. We found that faster processing speed on trials that require greater cognitive control are correlated with higher navigation efficiency (of both the whole-brain and the task-relevant subnetwork). In contrast, faster processing speed on trials that require more automatic processing are correlated with shorter path length within the task-relevant subnetwork. Our findings reveal that differences in the way communication is modelled between shortest path length and navigation may be sensitive to processing of automatic and controlled responses, respectively. Further, our findings suggest that there is a relationship between the speed of cognitive processing and the structural constraints of the human brain network.

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.

Evaluation of the diffusion MRI white matter tract integrity model using myelin histology and Monte-Carlo simulations

Quantitative evaluation of brain myelination has drawn considerable attention. Conventional diffusion-based magnetic resonance imaging models, including diffusion tensor imaging and diffusion kurtosis imaging (DKI),¹ have been used to infer the microstructure and its changes in neurological diseases. White matter tract integrity (WMTI) was proposed as a biophysical model to relate the DKI-derived metrics to the underlying microstructure. Although the model has been validated on ex vivo animal brains, it was not well evaluated with ex vivo human brains. In this study, histological samples (namely corpus callosum) from postmortem human brains have been investigated based on WMTI analyses on a clinical 3T scanner and comparisons with gold standard myelin staining in proteolipid protein and Luxol fast blue. In addition, Monte Carlo simulations were conducted to link changes from ex vivo to in vivo conditions based on the microscale parameters of water diffusivity and permeability. The results show that WMTI metrics, including axonal water fraction AWF, radial extra-axonal diffusivity D_e^⊥, and intra-axonal diffusivity D_awere needed to characterize myelin content alterations. Thus, WMTI model metrics are shown to be promising candidates as sensitive biomarkers of demyelination.

NODDI in clinical research

Diffusion MRI (dMRI) has proven to be a useful imaging approach for both clinical diagnosis and research investigating the microstructures of nervous tissues, and it has helped us to better understand the neurophysiological mechanisms of many diseases. Though diffusion tensor imaging (DTI) has long been the default tool to analyze dMRI data in clinical research, acquisition with stronger diffusion weightings beyond the DTI regimen is now possible with modern clinical scanners, potentially enabling even more detailed characterization of tissue microstructures. To take advantage of such data, neurite orientation dispersion and density imaging (NODDI) has been proposed as a way to relate the dMRI signal to tissue features via biophysically inspired modeling. The number of reports demonstrating the potential clinical utility of NODDI is rapidly increasing. At the same time, the pitfalls and limitations of NODDI, and general challenges in microstructure modeling, are becoming increasingly recognized by clinicians. dMRI microstructure modeling is a rapidly evolving field with great promise, where people from different scientific backgrounds, such as physics, medicine, biology, neuroscience, and statistics, are collaborating to build novel tools that contribute to improving human healthcare. Here, we review the applications of NODDI in clinical research and discuss future perspectives for investigations toward the implementation of dMRI microstructure imaging in clinical practice.

Constructing an Axonal-Specific Myelin Developmental Graph and its Application to Childhood Absence Epilepsy

BACKGROUND AND PURPOSE

The process of myelination starts in utero around 20 weeks of gestation and continues through adulthood. We first set out to characterize the maturation of the tract-specific myelin content in healthy subjects from childhood (7-12 years) into adulthood (18-32 years). Second, we apply the resulting development graph to children with childhood absence epilepsy (CAE), a pediatric epilepsy that was previously characterized by changes in myelin content.

METHODS

In a prospective cross-sectional study, 15 healthy children (7-12 years), 14 healthy adult participants (18-32 years) and 17 children with a clinical diagnosis of CAE (6-12 years) were included. For each participant, diffusion weighted images were acquired to reconstruct bundles of white matter tracts and multi-echo multi-slice GRASE images were acquired for myelin-water estimation. Subsequently, a tract-specific myelin development graph was constructed using the percentual difference in myelin-water content from childhood (12 year) to adulthood (25 year).

RESULTS

The graph revealed myelination patterns, where tracts in the central regions myelinate prior to peripheral tracts and intra-hemispheric tracts as well as tracts in the left hemisphere myelinate prior to inter-hemispheric tracts and tracts in the right hemisphere, respectively. No significant differences were found in myelin-water content between children with CAE and healthy children for neither the early developing tracts, nor the tracts that develop in a later stage. However, the difference between the myelin-water of late and early developing tracts is significantly smaller in the children with CAE.

CONCLUSION

These results indicate that CAE is associated with widespread neurodevelopmental myelin differences.

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.

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.

Microstructure-Informed Connectomics: Enriching Large-Scale Descriptions of Healthy and Diseased Brains

Rapid advances in neuroimaging and network science have produced powerful tools and measures to appreciate human brain organization at multiple spatial and temporal scales. It is now possible to obtain increasingly meaningful representations of whole-brain structural and functional brain networks and to formally assess macroscale principles of network topology. In addition to its utility in characterizing healthy brain organization, individual variability, and life span-related changes, there is high promise of network neuroscience for the conceptualization and, ultimately, management of brain disorders. In the current review, we argue for a science of the human brain that, while strongly embracing macroscale connectomics, also recommends awareness of brain properties derived from meso- and microscale resolutions. Such features include MRI markers of tissue microstructure, local functional properties, as well as information from nonimaging domains, including cellular, genetic, or chemical data. Integrating these measures with connectome models promises to better define the individual elements that constitute large-scale networks, and clarify the notion of connection strength among them. By enriching the description of large-scale networks, this approach may improve our understanding of fundamental principles of healthy brain organization. Notably, it may also better define the substrate of prevalent brain disorders, including stroke, autism, as well as drug-resistant epilepsies that are each characterized by intriguing interactions between local anomalies and network-level perturbations.

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.