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.
Traditional Diffusion Tensor Imaging (DTI) metrics are affected by crossing fibers and lesions. Most of the previous tractometry works use the single diffusion tensor, which leads to limited sensitivity and challenging interpretation of the results in crossing fiber regions. In this work, we propose a tractometry pipeline that combines white matter tractography with multi-tensor fixel-based metrics. These multi-tensors are estimated using the stable, accurate and robust to noise Multi-Resolution Discrete Search method (MRDS). The spatial coherence of the multi-tensor field estimated with MRDS, which includes up to three anisotropic and one isotropic tensors, is tractography-regularized using the Track Orientation Density Imaging method. Our end-to-end tractometry pipeline goes from raw data to track-specific multi-tensor-metrics tract profiles that are robust to noise and crossing fibers. A comprehensive evaluation conducted in a phantom simulating healthy and damaged tissue with the standard model, as well as in a healthy cohort of 20 individuals scanned along 5 time points, demonstrates the advantages of using multi-tensor metrics over traditional single-tensor metrics in tractometry. Qualitative assessment in a cohort of patients with relapsing-remitting multiple sclerosis reveals that the pipeline effectively detects white matter anomalies in the presence of crossing fibers and lesions.
Axon Diameter Distributions (ADDs) change during brain development and are altered in several brain pathologies. Mapping ADDs non-invasively using dMRI could provide a useful biomarker, but existing methods are either parametric, orientation dependent, surmmarize the whole ADD as a single measure or use non-standard protocols. We propose to estimate the ADD from an orientation-invariant PGSE protocol optimized for axon diameter sensitivity, using a discrete linear model with smoothness and sparsity regularization. To our knowledge, we are the first to report orientationally invarant ADD estimates from dMRI data.
Abstract Diffusion MRI fiber tractography is widely used to probe the structural connectivity of thebrain, with a range of applications in both clinical and basic neuroscience. Despite widespread use, tractography has well-known pitfalls that limits the anatomical accuracy of this technique. Numerous modern methods have been developed to address these shortcomings through advances in acquisition, modeling, and computation. To test whether these advances improve tractography accuracy, we organized the ISBI 2018 3D Validation of Tractography with Experimental MRI (3D-VoTEM) challenge. We made available three unique independent tractography validation datasets – a physical phantom and two ex vivo brain specimens - resulting in 176 distinct submissions from 9 research groups. By comparing results over a wide range of fiber complexities and algorithmic strategies, this challenge provides a more comprehensive assessment of tractography’s inherent limitations than has been reported previously. The central results were consistent across all sub-challenges in that, despite advances in tractography methods, the anatomical accuracy of tractography has not dramatically improved in recent years. Taken together, our results independently confirm findings from decades of tractography validation studies, demonstrate inherent limitations in reconstructing white matter pathways using diffusion MRI data alone, and highlight the need for alternative or combinatorial strategies to accurately map the fiber pathways of the brain.
Abstract Tractography is a family of algorithms that use diffusion-weighted magnetic resonance imaging data to reconstruct the white matter pathways of the brain. Although it has been proven to be particularly effective for studying non-invasively the neuronal architecture of the brain, recent studies have highlighted that the large incidence of false positive connections retrieved by these techniques can significantly bias any connectivity analysis. Some solutions have been proposed to overcome this issue and the ones relying on convex optimization framework showed a significant improvement. Here we propose an evolution of the Convex Optimization Modeling for Microstructure Informed Tractography (COMMIT) framework, that combines basic prior knowledge about brain anatomy with group-sparsity regularization into the optimization problem. We show that the new formulation dramatically reduces the incidence of false positives in synthetic DW-MRI data.
The presence of cortical lesions in multiple sclerosis patients has emerged as an important biomarker of the disease. They appear in the earliest stages of the illness and have been shown to correlate with the severity of clinical symptoms. However, cortical lesions are hardly visible in conventional magnetic resonance imaging (MRI) at 3T, and thus their automated detection has been so far little explored. In this study, we propose a fully-convolutional deep learning approach, based on the 3D U-Net, for the automated segmentation of cortical and white matter lesions at 3T. For this purpose, we consider a clinically plausible MRI setting consisting of two MRI contrasts only: one conventional T2-weighted sequence (FLAIR), and one specialized T1-weighted sequence (MP2RAGE). We include 90 patients from two different centers with a total of 728 and 3856 gray and white matter lesions, respectively. We show that two reference methods developed for white matter lesion segmentation are inadequate to detect small cortical lesions, whereas our proposed framework is able to achieve a detection rate of 76% for both cortical and white matter lesions with a false positive rate of 29% in comparison to manual segmentation. Further results suggest that our framework generalizes well for both types of lesion in subjects acquired in two hospitals with different scanners.
Purpose Non‐invasive axon diameter distribution (ADD) mapping using diffusion MRI is an ill‐posed problem. Current ADD mapping methods require knowledge of axon orientation before performing the acquisition. Instead, ActiveAx uses a 3D sampling scheme to estimate the orientation from the signal, providing orientationally invariant estimates. The mean diameter is estimated instead of the distribution for the solution to be tractable. Here, we propose an extension (ActiveAx ADD ) that provides non‐parametric and orientationally invariant estimates of the whole distribution. Theory The accelerated microstructure imaging with convex optimization (AMICO) framework accelerates mean diameter estimation using a linear formulation combined with Tikhonov regularization to stabilize the solution. Here, we implement a new formulation (ActiveAx ADD ) that uses Laplacian regularization to provide robust estimates of the whole ADD. Methods The performance of ActiveAx ADD was evaluated using Monte Carlo simulations on synthetic white matter samples mimicking axon distributions reported in histological studies. Results ActiveAx ADD provided robust ADD reconstructions when considering the isolated intra‐axonal signal. However, our formulation inherited some common microstructure imaging limitations. When accounting for the extra axonal compartment, estimated ADDs showed spurious peaks and increased variability because of the difficulty of disentangling intra and extra axonal contributions. Conclusion Laplacian regularization solves the ill‐posedness regarding the intra axonal compartment. ActiveAx ADD can potentially provide non‐parametric and orientationally invariant ADDs from isolated intra‐axonal signals. However, further work is required before ActiveAx ADD can be applied to real data containing extra‐axonal contributions, as disentangling the 2 compartment appears to be an overlooked challenge that affects microstructure imaging methods in general.
The choroid plexus has been shown to play a crucial role in CNS inflammation. Previous studies found larger choroid plexus in multiple sclerosis (MS) compared with healthy controls. However, it is not clear whether the choroid plexus is similarly involved in MS and in neuromyelitis optica spectrum disorder (NMOSD). Thus, the aim of this study was to compare the choroid plexus volume in MS and NMOSD.In this retrospective, cross-sectional study, patients were included by convenience sampling from 4 international MS centers. The choroid plexus of the lateral ventricles was segmented fully automatically on T1-weighted MRI sequences using a deep learning algorithm (Multi-Dimensional Gated Recurrent Units). Uni- and multivariable linear models were applied to investigate associations between the choroid plexus volume, clinically meaningful disease characteristics, and MRI parameters.We studied 180 patients with MS and 98 patients with NMOSD. In total, 94 healthy individuals and 47 patients with migraine served as controls. The choroid plexus volume was larger in MS (median 1,690 µL, interquartile range [IQR] 648 µL) than in NMOSD (median 1,403 µL, IQR 510 µL), healthy individuals (median 1,533 µL, IQR 570 µL), and patients with migraine (median 1,404 µL, IQR 524 µL; all p < 0.001), whereas there was no difference between NMOSD, migraine, and healthy controls. This was also true when adjusted for age, sex, and the intracranial volume. In contrast to NMOSD, the choroid plexus volume in MS was associated with the number of T2-weighted lesions in a linear model adjusted for age, sex, total intracranial volume, disease duration, relapses in the year before MRI, disease course, Expanded Disability Status Scale score, disease-modifying treatment, and treatment duration (beta 4.4; 95% CI 0.78-8.1; p = 0.018).This study supports an involvement of the choroid plexus in MS in contrast to NMOSD and provides clues to better understand the respective pathogenesis.
We performed an extensive assessment of the clinical relevance of a method that we had previously developed, which provides personalized quantitative MRI abnormality maps of individual multiple sclerosis (MS) patients. Specifically, we assessed the relationships between quantitative T1 (qT1), myelin water fraction (MWF), neurite density index (NDI), magnetization transfer saturation (MTsat) abnormality maps and clinical disability in a cohort of 102 MS patients and 98 healthy subjects. We found that qT1 and NDI alterations in white matter lesions were strongly related to patients' clinical disability, supporting the use of those personalized maps for patient stratification and follow-up in clinical practice.
In-vivo quantification of axon diameter is an attractive and debated topic in the MRI community. The possibility to resolve submicrometric axon diameters non-invasively yields the potential to push further the boundaries in research and clinics but yet, further work is needed to better explore and validate the existing approaches to estimate the inner axon diameter. Recently, the feasibility of estimating the axon diameter from the intra-axonal transverse relaxation time has been investigated combining a diffusion-relaxation protocol and histological data. In the present study, we apply this approach in a larger in vivo population to assess variability across participants.