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 whole-brain 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.

Diffusion-based biophysical models have been used in several recent works to study the microenvironment of brain tumours. While the pathophysiological interpretation of the parameters of these models remains unclear, their use as signal representations may yield useful biomarkers for monitoring the treatment and the progression of this complex and heterogeneous disease. Up to now, however, no study was devoted to assessing the mathematical stability of these approaches in cancerous brain regions. To this end, we analyzed in 11 brain tumour patients the fitting results of two microstructure models (Neurite Orientation Dispersion and Density Imaging and the Spherical Mean Technique) and of a signal representation (Diffusion Kurtosis Imaging) to compare the reliability of their parameter estimates in the healthy brain and in the tumoral lesion. The framework of our between-tissue analysis included the computation of 1) the residual sum of squares as a goodness-of-fit measure 2) the standard deviation of the models' derived metrics and 3) models' sensitivity functions to analyze the suitability of the employed protocol for parameter estimation in the different microenvironments. Our results revealed no issues concerning the fitting of the models in the tumoral lesion, with similar goodness of fit and parameter precisions occurring in normal appearing and pathological tissues. Lastly, with the aim of highlight possible biomarkers, in our analysis we briefly discuss the correlation between the metrics of the three techniques, identifying groups of indices which are significantly collinear in all tissues and thus provide no additional information when jointly used in data-driven analyses.

Diffusion magnetic resonance imaging (dMRI) tractography is an advanced imaging technique that enables in vivo reconstruction of the brain's white matter connections at macro scale. It provides an important tool for quantitative mapping of the brain's structural connectivity using measures of connectivity or tissue microstructure. Over the last two decades, the study of brain connectivity using dMRI tractography has played a prominent role in the neuroimaging research landscape. In this paper, we provide a high-level overview of how tractography is used to enable quantitative analysis of the brain's structural connectivity in health and disease. We focus on two types of quantitative analyses of tractography, including: 1) tract-specific analysis that refers to research that is typically hypothesis-driven and studies particular anatomical fiber tracts, and 2) connectome-based analysis that refers to research that is more data-driven and generally studies the structural connectivity of the entire brain. We first provide a review of methodology involved in three main processing steps that are common across most approaches for quantitative analysis of tractography, including methods for tractography correction, segmentation and quantification. For each step, we aim to describe methodological choices, their popularity, and potential pros and cons. We then review studies that have used quantitative tractography approaches to study the brain's white matter, focusing on applications in neurodevelopment, aging, neurological disorders, mental disorders, and neurosurgery. We conclude that, while there have been considerable advancements in methodological technologies and breadth of applications, there nevertheless remains no consensus about the "best" methodology in quantitative analysis of tractography, and researchers should remain cautious when interpreting results in research and clinical applications.

Neuropathological studies have shown that multiple sclerosis (MS) lesions are heterogeneous in terms of myelin/axon damage and repair as well as iron content. However, it remains a challenge to identify specific chronic lesion types, especially remyelinated lesions, in vivo in patients with MS.

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.

Although structural disconnection represents the hallmark of multiple sclerosis (MS) pathophysiology, classification attempts based on structural connectivity have achieved low accuracy levels. Here, we set out to fill this gap, exploring the performance of supervised classifiers on features derived from microstructure informed tractography and selected applying a novel robust approach.

Damage to the myelin sheath and the neuroaxonal unit is a cardinal feature of multiple sclerosis; however, a detailed characterization of the interaction between myelin and axon damage in vivo remains challenging. We applied myelin water and multi-shell diffusion imaging to quantify the relative damage to myelin and axons (i) among different lesion types; (ii) in normal-appearing tissue; and (iii) across multiple sclerosis clinical subtypes and healthy controls. We also assessed the relation of focal myelin/axon damage with disability and serum neurofilament light chain as a global biological measure of neuroaxonal damage. Ninety-one multiple sclerosis patients (62 relapsing-remitting, 29 progressive) and 72 healthy controls were enrolled in the study. Differences in myelin water fraction and neurite density index were substantial when lesions were compared to healthy control subjects and normal-appearing multiple sclerosis tissue: both white matter and cortical lesions exhibited a decreased myelin water fraction and neurite density index compared with healthy (P < 0.0001) and peri-plaque white matter (P < 0.0001). Periventricular lesions showed decreased myelin water fraction and neurite density index compared with lesions in the juxtacortical region (P < 0.0001 and P < 0.05). Similarly, lesions with paramagnetic rims showed decreased myelin water fraction and neurite density index relative to lesions without a rim (P < 0.0001). Also, in 75% of white matter lesions, the reduction in neurite density index was higher than the reduction in the myelin water fraction. Besides, normal-appearing white and grey matter revealed diffuse reduction of myelin water fraction and neurite density index in multiple sclerosis compared to healthy controls (P < 0.01). Further, a more extensive reduction in myelin water fraction and neurite density index in normal-appearing cortex was observed in progressive versus relapsing-remitting participants. Neurite density index in white matter lesions correlated with disability in patients with clinical deficits (P < 0.01, beta = -10.00); and neurite density index and myelin water fraction in white matter lesions were associated to serum neurofilament light chain in the entire patient cohort (P < 0.01, beta = -3.60 and P < 0.01, beta = 0.13, respectively). These findings suggest that (i) myelin and axon pathology in multiple sclerosis is extensive in both lesions and normal-appearing tissue; (ii) particular types of lesions exhibit more damage to myelin and axons than others; (iii) progressive patients differ from relapsing-remitting patients because of more extensive axon/myelin damage in the cortex; and (iv) myelin and axon pathology in lesions is related to disability in patients with clinical deficits and global measures of neuroaxonal damage.

Graph theory and network modelling have been previously applied to characterize motor network structural topology in multiple sclerosis (MS). However, between-group differences disclosed by graph analysis might be primarily driven by discrepancy in density, which is likely to be reduced in pathologic conditions as a consequence of macroscopic damage and fibre loss that may result in less streamlines properly traced. In this work, we employed the convex optimization modelling for microstructure informed tractography (COMMIT) framework, which, given a tractogram, estimates the actual contribution (or weight) of each streamline in order to optimally explain the diffusion magnetic resonance imaging signal, filtering out those that are implausible or not necessary. Then, we analysed the topology of this 'COMMIT-weighted sensory-motor network' in MS accounting for network density. By comparing with standard connectivity analysis, we also tested if abnormalities in network topology are still identifiable when focusing on more 'quantitative' network properties. We found that topology differences identified with standard tractography in MS seem to be mainly driven by density, which, in turn, is strongly influenced by the presence of lesions. We were able to identify a significant difference in density but also in network global and local properties when accounting for density discrepancy. Therefore, we believe that COMMIT may help characterize the structural organization in pathological conditions, allowing a fair comparison of connectomes which considers discrepancies in network density. Moreover, discrepancy-corrected network properties are clinically meaningful and may help guide prognosis assessment and treatment choice.

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.

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.

Limitations in the accuracy of brain pathways reconstructed by diffusion MRI (dMRI) tractography have received considerable attention. While the technical advances spearheaded by the Human Connectome Project (HCP) led to significant improvements in dMRI data quality, it remains unclear how these data should be analyzed to maximize tractography accuracy. Over a period of two years, we have engaged the dMRI community in the IronTract Challenge, which aims to answer this question by leveraging a unique dataset. Macaque brains that have received both tracer injections and ex vivo dMRI at high spatial and angular resolution allow a comprehensive, quantitative assessment of tractography accuracy on state-of-the-art dMRI acquisition schemes. We find that, when analysis methods are carefully optimized, the HCP scheme can achieve similar accuracy as a more time-consuming, Cartesian-grid scheme. Importantly, we show that simple pre- and post-processing strategies can improve the accuracy and robustness of many tractography methods. Finally, we find that fiber configurations that go beyond crossing (e.g., fanning, branching) are the most challenging for tractography. The IronTract Challenge remains open and we hope that it can serve as a valuable validation tool for both users and developers of dMRI analysis methods.

Diffusion magnetic resonance imaging is a noninvasive imaging modality that has been extensively used in the literature to study the neuronal architecture of the brain in a wide range of neurological conditions using tractography. However, recent studies highlighted that the anatomical accuracy of the reconstructions is inherently limited and challenged its appropriateness. Several solutions have been proposed to tackle this issue, but none of them proved effective to overcome this fundamental limitation. In this work, we present a novel processing framework to inject into the reconstruction problem basic prior knowledge about brain anatomy and its organization and evaluate its effectiveness using both simulated and real human brain data. Our results indicate that our proposed method dramatically increases the accuracy of the estimated brain networks and, thus, represents a major step forward for the study of connectivity.

Recent evidences have suggested the possible presence of an involvement of the extrapyramidal system in Fabry disease (FD), a rare X-linked lysosomal storage disorder. We aimed to investigate the microstructural integrity of the main tracts of the cortico-striatal-thalamo-cortical loop in FD patients.

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.

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.

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.

Soma and neurite density image (SANDI) is an advanced diffusion magnetic resonance imaging biophysical signal model devised to probe in vivo microstructural information in the gray matter (GM). This model requires acquisitions that include b values that are at least six times higher than those used in clinical practice. Such high b values are required to disentangle the signal contribution of water diffusing in soma from that diffusing in neurites and extracellular space, while keeping the diffusion time as short as possible to minimize potential bias due to water exchange. These requirements have limited the use of SANDI only to preclinical or cutting-edge human scanners. Here, we investigate the potential impact of neglecting water exchange in the SANDI model and present a 10-min acquisition protocol that enables to characterize both GM and white matter (WM) on 3 T scanners. We implemented analytical simulations to (i) evaluate the stability of the fitting of SANDI parameters when diminishing the number of shells; (ii) estimate the bias due to potential exchange between neurites and extracellular space in such reduced acquisition scheme, comparing it with the bias due to experimental noise. Then, we demonstrated the feasibility and assessed the repeatability and reproducibility of our approach by computing microstructural metrics of SANDI with AMICO toolbox and other state-of-the-art models on five healthy subjects. Finally, we applied our protocol to five multiple sclerosis patients. Results suggest that SANDI is a practical method to characterize WM and GM tissues in vivo on performant clinical scanners.