Highlights • The functional characterization of callosal connections is informed by anatomical data.• Callosal connections play a conditional driving role depending on the brain state and behavioral demands.• Callosal connections play a modulatory function, in addition to a driving role.• The corpus callosum participates in learning and interhemispheric transfer of sensorimotor habits.• The corpus callosum contributes to language processing and cognitive functions.

The limited availability of large image datasets, mainly due to data privacy and differences in acquisition protocols or hardware, is a significant issue in the development of accurate and generalizable machine learning methods in medicine. This is especially the case for Magnetic Resonance (MR) images, where different MR scanners introduce a bias that limits the performance of a machine learning model. We present a novel method that learns to ignore the scanner-related features present in MR images, by introducing specific additional constraints on the latent space. We focus on a real-world classification scenario, where only a small dataset provides images of all classes. Our method \textit{Learn to Ignore (L2I)} outperforms state-of-the-art domain adaptation methods on a multi-site MR dataset for a classification task between multiple sclerosis patients and healthy controls.

Introduction: Graph theory has been applied to study the pathophysiology of multiple sclerosis (MS) since it provides global and focal measures of brain network properties that are affected by MS. Typically, the connection strength and, consequently, the network properties are computed by counting the number of streamlines (NOS) connecting couples of gray matter regions. However, recent studies have shown that this method is not quantitative. Methods: We evaluated diffusion-based microstructural measures extracted from three different models to assess the network properties in a group of 66 MS patients and 64 healthy subjects. Besides, we assessed their correlation with patients' disability and with a biological measure of neuroaxonal damage. Results: Graph metrics extracted from connectomes weighted by intra-axonal microstructural components were the most sensitive to MS pathology and the most related to clinical disability. In contrast, measures of network segregation extracted from the connectomes weighted by maps describing extracellular diffusivity were the most related to serum concentration of neurofilament light chain. Network properties assessed with NOS were neither sensitive to MS pathology nor correlated with clinical and pathological measures of disease impact in MS patients. Conclusion: Using tractometry-derived graph measures in MS patients, we identified a set of metrics based on microstructural components that are highly sensitive to the disease and that provide sensitive correlates of clinical and biological deterioration in MS patients. Impact statement Graph theory has been widely used to study the alterations in the structural connectivity of multiple sclerosis (MS) patients. Usually, brain graphs used for the extraction of network metrics are created by counting the number of streamlines connecting gray matter regions, however, this method is not quantitative. In this study, we used tractometry to average the values of diffusion-based microstructural maps along the reconstructed streamlines. Our results show that network metrics extracted from the connectomes weighted on microstructural maps provide sensitive information to MS pathology, which correlate with clinical and biological measures of disease impact.

In the central nervous system of primates, several pathways are characterized by different spectra of axon diameters. In vivo methods, based on diffusion-weighted magnetic resonance imaging, can provide axon diameter index estimates non-invasively. However, such methods report voxel-wise estimates, which vary from voxel-to-voxel for the same white matter bundle due to partial volume contributions from other pathways having different microstructure properties. Here, we propose a novel microstructure-informed tractography approach, COMMITAxSize, to resolve axon diameter index estimates at the streamline level, thus making the estimates invariant along trajectories. Compared to previously proposed voxel-wise methods, our formulation allows the estimation of a distinct axon diameter index value for each streamline, directly, furnishing a complementary measure to the existing calculation of the mean value along the bundle. We demonstrate the favourable performance of our approach comparing our estimates with existing histologically-derived measurements performed in the corpus callosum and the posterior limb of the internal capsule. Overall, our method provides a more robust estimation of the axon diameter index of pathways by jointly estimating the microstructure properties of the tissue and the macroscopic organisation of the white matter connectivity.

We aimed to determine in relapsing multiple sclerosis (MS) whether intrathecal synthesis of immunoglobulin (Ig) M and IgG is associated with outcomes reflecting inflammatory activity and chronic worsening.

The anisotropy of brain white matter microstructure manifests itself in orientational-dependence of various MRI contrasts, and can result in significant quantification biases if ignored. Understanding the origins of this orientation-dependence could enhance the interpretation of MRI signal changes in development, ageing and disease and ultimately improve clinical diagnosis. Using a novel experimental setup, this work studies the contributions of the intra- and extra-axonal water to the orientation-dependence of one of the most clinically-studied parameters, apparent transverse relaxation T2. Specifically, a tiltable receive coil is interfaced with an ultra-strong gradient MRI scanner to acquire multidimensional MRI data with an unprecedented range of acquisition parameters. Using this setup, compartmental T2 can be disentangled based on differences in diffusional-anisotropy, and its orientation-dependence further elucidated by re-orienting the head with respect to the main magnetic field B→0. A dependence of (compartmental) T2 on the fibre orientation w.r.t. B→0 was observed, and further quantified using characteristic representations for susceptibility- and magic angle effects. Across white matter, anisotropy effects were dominated by the extra-axonal water signal, while the intra-axonal water signal decay varied less with fibre-orientation. Moreover, the results suggest that the stronger extra-axonal T2 orientation-dependence is dominated by magnetic susceptibility effects (presumably from the myelin sheath) while the weaker intra-axonal T2 orientation-dependence may be driven by a combination of microstructural effects. Even though the current design of the tiltable coil only offers a modest range of angles, the results demonstrate an overall effect of tilt and serve as a proof-of-concept motivating further hardware development to facilitate experiments that explore orientational anisotropy. These observations have the potential to lead to white matter microstructural models with increased compartmental sensitivity to disease, and can have direct consequences for longitudinal and group-wise T2- and diffusion-MRI data analysis, where the effect of head-orientation in the scanner is commonly ignored.

Conventional magnetic resonance imaging (cMRI) in multiple sclerosis (MS) patients provides measures of focal brain damage and activity, which are fundamental for disease diagnosis, prognosis, and the evaluation of response to therapy. However, cMRI is insensitive to the damage to the microenvironment of the brain tissue and the heterogeneity of MS lesions. In contrast, the damaged tissue can be characterized by mathematical models on multishell diffusion imaging data, which measure different compartmental water diffusion. In this work, we obtained 12 diffusion measures from eight diffusion models, and we applied a deep-learning attention-based convolutional neural network (CNN) (GAMER-MRI) to select the most discriminating measures in the classification of MS lesions and the perilesional tissue by attention weights. Furthermore, we provided clinical and biological validation of the chosen metrics—and of their most discriminative combinations—by correlating their respective mean values in MS patients with the corresponding Expanded Disability Status Scale (EDSS) and the serum level of neurofilament light chain (sNfL), which are measures of disability and neuroaxonal damage. Our results show that the neurite density index from neurite orientation and dispersion density imaging (NODDI), the measures of the intra-axonal and isotropic compartments from microstructural Bayesian approach, and the measure of the intra-axonal compartment from the spherical mean technique NODDI were the most discriminating (respective attention weights were 0.12, 0.12, 0.15, and 0.13). In addition, the combination of the neurite density index from NODDI and the measures for the intra-axonal and isotropic compartments from the microstructural Bayesian approach exhibited a stronger correlation with EDSS and sNfL than the individual measures. This work demonstrates that the proposed method might be useful to select the microstructural measures that are most discriminative of focal tissue damage and that may also be combined to a unique contrast to achieve stronger correlations to clinical disability and neuroaxonal damage.

Background: MR imaging of the spinal cord (SC) gray matter (GM) at the cervical and lumbar enlargements' level may be particularly informative in lower motor neuron disorders, e. g., spinal muscular atrophy, but also in other neurodegenerative or autoimmune diseases affecting the SC. Radially sampled averaged magnetization inversion recovery acquisition (rAMIRA) is a novel approach to perform SC imaging in clinical settings with favorable contrast and is well-suited for SC GM quantitation. However, before applying rAMIRA in clinical studies, it is important to understand (i) the sources of inter-subject variability of total SC cross-sectional areas (TCA) and GM area (GMA) measurements in healthy subjects and (ii) their relation to age and sex to facilitate the detection of pathology-associated changes. In this study, we aimed to develop normalization strategies for rAMIRA-derived SC metrics using skull and spine-based metrics to reduce anatomical variability. Methods: Sixty-one healthy subjects (age range 11–93 years, 37.7% women) were investigated with axial two-dimensional rAMIRA imaging at 3T MRI. Cervical and thoracic levels including the level of the cervical (C4/C5) and lumbar enlargements (Tmax) were examined. SC T2-weighted sagittal images and high-resolution 3D whole-brain T1-weighted images were acquired. TCA and GMAs were quantified. Anatomical variables with associations of |r| > 0.30 in univariate association with SC areas, and age and sex were used to construct normalization models using backward selection with TCAC4/C5 as outcome. The effect of the normalization was assessed by % relative standard deviation (RSD) reductions. Results: Mean inter-individual variability and the SD of the SC area metrics were considerable: TCAC4/5: 8.1%/9.0; TCATmax: 8.9%/6.5; GMAC4/C5: 8.6%/2.2; GMATmax: 12.2%/3.8. Normalization based on sex, brain WM volume, and spinal canal area resulted in RSD reductions of 23.7% for TCAs and 12.0% for GM areas at C4/C5. Normalizations based on the area of spinal canal alone resulted in RSD reductions of 10.2% for TCAs and 9.6% for GM areas at C4/C5, respectively. Discussion: Anatomic inter-individual variability of SC areas is substantial. This study identified effective normalization models for inter-subject variability reduction in TCA and SC GMA in healthy subjects based on rAMIRA imaging.

Rahmanzadeh et al. quantify—for the first time in vivo—the relative damage to myelin and axons across the multiple sclerosis spectrum, in different types of lesions and in normal-appearing tissue. The data confirm neuropathological findings and extend them by revealing the complexity of the disease in living patients.

Diffusion MRI (dMRI) represents one of the few methods for mapping brain fiber orientations non-invasively. Unfortunately, dMRI fiber mapping is an indirect method that relies on inference from measured diffusion patterns. Comparing dMRI results with other modalities is a way to improve the interpretation of dMRI data and help advance dMRI technologies. Here, we present methods for comparing dMRI fiber orientation estimates with optical imaging of fluorescently labeled neurofilaments and vasculature in 3D human and primate brain tissue cuboids cleared using CLARITY. The recent advancements in tissue clearing provide a new opportunity to histologically map fibers projecting in 3D, which represents a captivating complement to dMRI measurements. In this work, we demonstrate the capability to directly compare dMRI and CLARITY in the same human brain tissue and assess multiple approaches for extracting fiber orientation estimates from CLARITY data. We estimate the three-dimensional neuronal fiber and vasculature orientations from neurofilament and vasculature stained CLARITY images by calculating the tertiary eigenvector of structure tensors. We then extend CLARITY orientation estimates to an orientation distribution function (ODF) formalism by summing multiple sub-voxel structure tensor orientation estimates. In a sample containing part of the human thalamus, there is a mean angular difference of 19° ± 15° between the primary eigenvectors of the dMRI tensors and the tertiary eigenvectors from the CLARITY neurofilament stain. We also demonstrate evidence that vascular compartments do not affect the dMRI orientation estimates by showing an apparent lack of correspondence (mean angular difference = 49° ± 23°) between the orientation of the dMRI tensors and the structure tensors in the vasculature stained CLARITY images. In a macaque brain dataset, we examine how the CLARITY feature extraction depends on the chosen feature extraction parameters. By varying the volume of tissue over which the structure tensor estimates are derived, we show that orientation estimates are noisier with more spurious ODF peaks for sub-voxels below 30 μm3 and that, for our data, the optimal gray matter sub-voxel size is between 62.5 μm 3 and 125 μm 3. The example experiments presented here represent an important advancement towards robust multi-modal MRI-CLARITY comparisons.

At the typical spatial resolution of MRI in the human brain, approximately 60–90% of voxels contain multiple fiber populations. Quantifying microstructural properties of distinct fiber populations within a voxel is therefore challenging but necessary. While progress has been made for diffusion and T1-relaxation properties, how to resolve intra-voxel T2 heterogeneity remains an open question. Here a novel framework, named COMMIT-T2, is proposed that uses tractography-based spatial regularization with diffusion-relaxometry data to estimate multiple intra-axonal T2 values within a voxel. Unlike previously-proposed voxel-based T2 estimation methods, which (when applied in white matter) implicitly assume just one fiber bundle in the voxel or the same T2 for all bundles in the voxel, COMMIT-T2 can recover specific T2 values for each unique fiber population passing through the voxel. In this approach, the number of recovered unique T2 values is not determined by a number of model parameters set a priori, but rather by the number of tractography-reconstructed streamlines passing through the voxel. Proof-of-concept is provided in silico and in vivo, including a demonstration that distinct tract-specific T2 profiles can be recovered even in the three-way crossing of the corpus callosum, arcuate fasciculus, and corticospinal tract. We demonstrate the favourable performance of COMMIT-T2 compared to that of voxelwise approaches for mapping intra-axonal T2 exploiting diffusion, including a direction-averaged method and AMICO-T2, a new extension to the previously-proposed Accelerated Microstructure Imaging via Convex Optimization (AMICO) framework.

White matter bundle segmentation using diffusion MRI fiber tractography has become the method of choice to identify white matter fiber pathways in vivo in human brains. However, like other analyses of complex data, there is considerable variability in segmentation protocols and techniques. This can result in different reconstructions of the same intended white matter pathways, which directly affects tractography results, quantification, and interpretation. In this study, we aim to evaluate and quantify the variability that arises from different protocols for bundle segmentation. Through an open call to users of fiber tractography, including anatomists, clinicians, and algorithm developers, 42 independent teams were given processed sets of human wholebrain streamlines and asked to segment 14 white matter fascicles on six subjects. In total, we received 57 different bundle segmentation protocols, which enabled detailed volume-based and streamline-based analyses of agreement and disagreement among protocols for each fiber pathway. Results show that even when given the exact same sets of underlying streamlines, the variability across protocols for bundle segmentation is greater than all other sources of variability in the virtual dissection process, including variability within protocols and variability across subjects. In order to foster the use of tractography bundle dissection in routine clinical settings, and as a fundamental analytical tool, future endeavors must aim to resolve and reduce this heterogeneity. Although external validation is needed to verify the anatomical accuracy of bundle dissections, reducing heterogeneity is a step towards reproducible research and may be achieved through the use of standard nomenclature and definitions of white matter bundles and well-chosen constraints and decisions in the dissection process.