Search

The white matter structures of the human brain can be represented using diffusion-weighted MRI tractography. Unfortunately, tractography is prone to find false-positive streamlines causing a severe decline in its specificity and limiting its feasibility in accurate structural brain connectivity analyses. Filtering algorithms have been proposed to reduce the number of invalid streamlines but the currently available filtering algorithms are not suitable to process data that contains motion artefacts which are typical in clinical research. We augmented the Convex Optimization Modelling for Microstructure Informed Tractography (COMMIT) algorithm to adjust for these signals drop-out motion artefacts. We demonstrate with comprehensive Monte-Carlo whole brain simulations and in vivo infant data that our robust algorithm is capable of properly filtering tractography reconstructions despite these artefacts. We evaluated the results using parametric and non-parametric statistics and our results demonstrate that if not accounted for, motion artefacts can have severe adverse effects in human brain structural connectivity analyses as well as in microstructural property mappings. In conclusion, the usage of robust filtering methods to mitigate motion related errors in tractogram filtering is highly beneficial, especially in clinical studies with uncooperative patient groups such as infants. With our presented robust augmentation and open-source implementation, robust tractogram filtering is readily available., (Copyright © 2021. Published by Elsevier Inc.)

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., (© The Author(s) (2021). Published by Oxford University Press on behalf of the Guarantors of Brain.)

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 T 1 -relaxation properties, how to resolve intra-voxel T 2 heterogeneity remains an open question. Here a novel framework, named COMMIT-T 2 , is proposed that uses tractography-based spatial regularization with diffusion-relaxometry data to estimate multiple intra-axonal T 2 values within a voxel. Unlike previously-proposed voxel-based T 2 estimation methods, which (when applied in white matter) implicitly assume just one fiber bundle in the voxel or the same T 2 for all bundles in the voxel, COMMIT-T 2 can recover specific T 2 values for each unique fiber population passing through the voxel. In this approach, the number of recovered unique T 2 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 T 2 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-T 2 compared to that of voxelwise approaches for mapping intra-axonal T 2 exploiting diffusion, including a direction-averaged method and AMICO-T 2 , a new extension to the previously-proposed Accelerated Microstructure Imaging via Convex Optimization (AMICO) framework., (Copyright © 2020 The Authors. Published by Elsevier Inc. All rights reserved.)

Purpose: New analytical reconstruction techniques of diffusion weighted signal have been proposed. A previous work evidenced the exploitability of some indices derived from the simple harmonic oscillator-based reconstruction and estimation (3D-SHORE) model as numerical biomarkers of neural plasticity after stroke. Here, the analysis is extended to two additional indices: return to the plane/origin (RTPP/RTOP) probabilities. Moreover, several motor networks were introduced and the results were analyzed at different time scales., Methods: Ten patients underwent three diffusion spectrum imaging (DSI) scans [1 week (tp1), 1 month (tp2) and 6 months (tp3) after stroke]. Ten matched controls underwent two DSI scans 1 month apart. 3D-SHORE was used for reconstructing the signal and the microstructural indices were derived. Tract-based analysis was performed along motor cortical, subcortical and transcallosal networks in the contralesional area., Results: The optimal intra-class correlation coefficient (ICC) was obtained in the subcortical loop for propagator anisotropy (ICC [Formula: see text] 0.96), followed by generalized fractional anisotropy (ICC [Formula: see text] 0.94). The new indices reached the highest stability in the transcallosal network and performed well in the cortical and subcortical networks with the exception of RTOP in the cortical loop (ICC [Formula: see text] 0.59). They allowed discriminating patients from controls at the majority of the timescales. Finally, the regression model using indices calculated along the subcortical loop at tp1 resulted in the best prediction of clinical outcome., Conclusions: The whole set of microstructural indices provide measurements featuring high precision. The new indices allow discriminating patients from controls in all networks, except for RTPP in the cortical loop. Moreover, the 3D-SHORE indices in subcortical connections constitute a good regression model for predicting the clinical outcome at 6 months, supporting their suitability as numerical biomarkers for neuronal plasticity after stroke.

Diffusion weighted magnetic resonance signals convey information about tissue microstructure and cytoarchitecture. In the last years, many models have been proposed for recovering the diffusion signal and extracting information to constitute new families of numerical indices. Two main categories of reconstruction models can be identified in diffusion magnetic resonance imaging (DMRI): ensemble average propagator (EAP) models and compartmental models. From both, descriptors can be derived for elucidating the underlying microstructural architecture. While compartmental models indices directly quantify the fraction of different cell compartments in each voxel, EAP-derived indices are only a derivative measure and the effect of the different microstructural configurations on the indices is still unclear. In this paper, we analyze three EAP indices calculated using the 3D Simple Harmonic Oscillator based Reconstruction and Estimation (3D-SHORE) model and estimate their changes with respect to the principal microstructural configurations. We take advantage of the state of the art simulations to quantify the variations of the indices with the simulation parameters. Analysis of in-vivo data correlates the EAP indices with the microstructural parameters obtained from the Neurite Orientation Dispersion and Density Imaging (NODDI) model as a pseudo ground truth for brain data. Results show that the EAP derived indices convey information on the tissue microstructure and that their combined values directly reflect the configuration of the different compartments in each voxel., (Copyright © 2016 Elsevier B.V. All rights reserved.)

Objectives: In this study, we investigated the structural plasticity of the contralesional motor network in ischemic stroke patients using diffusion magnetic resonance imaging (MRI) and explored a model that combines a MRI-based metric of contralesional network integrity and clinical data to predict functional outcome at 6 months after stroke., Methods: MRI and clinical examinations were performed in 12 patients in the acute phase, at 1 and 6 months after stroke. Twelve age- and gender-matched controls underwent 2 MRIs 1 month apart. Structural remodeling after stroke was assessed using diffusion MRI with an automated measurement of generalized fractional anisotropy (GFA), which was calculated along connections between contralesional cortical motor areas. The predictive model of poststroke functional outcome was computed using a linear regression of acute GFA measures and the clinical assessment., Results: GFA changes in the contralesional motor tracts were found in all patients and differed significantly from controls (0.001 ≤ p < 0.05). GFA changes in intrahemispheric and interhemispheric motor tracts correlated with age (p ≤ 0.01); those in intrahemispheric motor tracts correlated strongly with clinical scores and stroke sizes (p ≤ 0.001). GFA measured in the acute phase together with a routine motor score and age were a strong predictor of motor outcome at 6 months (r(2) = 0.96, p = 0.0002)., Conclusion: These findings represent a proof of principle that contralesional diffusion MRI measures may provide reliable information for personalized rehabilitation planning after ischemic motor stroke.

Background: The role of the non-injured hemisphere in stroke recovery is poorly understood. In this pilot study, we sought to explore the presence of structural changes detectable by diffusion tensor imaging (DTI) in the contralesional hemispheres of patients who recovered well from ischemic stroke., Methods: We analyzed serial DTI data from 16 stroke patients who had moderate initial neurological deficits (NIHSS scores 3-12) and good functional outcome at 3-6 months (NIHSS score 0 or modified Rankin Score ≤1). We segmented the brain tissue in gray and white matter (GM and WM) and measured the apparent diffusion coefficient (ADC) and fractional anisotropy in the infarct, in the contralesional infarct mirror region as well as in concentrically expanding regions around them., Results: We found that GM and WM ADC significantly increased in the infarct region (p < 0.01) from acute to chronic time points, whereas in the infarct mirror region, GM and WM ADC increased (p < 0.01) and WM fractional anisotropy decreased (p < 0.05). No significant changes were detected in other regions., Conclusion: DTI-based metrics are sensitive to regional structural changes in the contralesional hemisphere during stroke recovery. Prospective studies in larger cohorts with varying levels of recovery are needed to confirm our findings., (Copyright © 2012 S. Karger AG, Basel.)

q-Space-based techniques such as diffusion spectrum imaging, q-ball imaging, and their variations have been used extensively in research for their desired capability to delineate complex neuronal architectures such as multiple fiber crossings in each of the image voxels. The purpose of this article was to provide an introduction to the q-space formalism and the principles of basic q-space techniques together with the discussion on the advantages as well as challenges in translating these techniques into the clinical environment. A review of the currently used q-space-based protocols in clinical research is also provided.

Background: Patients suffering from migraine with aura (MWA) and migraine without aura (MWoA) show abnormalities in visual motion perception during and between attacks. Whether this represents the consequences of structural changes in motion-processing networks in migraineurs is unknown. Moreover, the diagnosis of migraine relies on patient's history, and finding differences in the brain of migraineurs might help to contribute to basic research aimed at better understanding the pathophysiology of migraine., Methods and Findings: To investigate a common potential anatomical basis for these disturbances, we used high-resolution cortical thickness measurement and diffusion tensor imaging (DTI) to examine the motion-processing network in 24 migraine patients (12 with MWA and 12 MWoA) and 15 age-matched healthy controls (HCs). We found increased cortical thickness of motion-processing visual areas MT+ and V3A in migraineurs compared to HCs. Cortical thickness increases were accompanied by abnormalities of the subjacent white matter. In addition, DTI revealed that migraineurs have alterations in superior colliculus and the lateral geniculate nucleus, which are also involved in visual processing., Conclusions: A structural abnormality in the network of motion-processing areas could account for, or be the result of, the cortical hyperexcitability observed in migraineurs. The finding in patients with both MWA and MWoA of thickness abnormalities in area V3A, previously described as a source in spreading changes involved in visual aura, raises the question as to whether a "silent" cortical spreading depression develops as well in MWoA. In addition, these experimental data may provide clinicians and researchers with a noninvasively acquirable migraine biomarker.

IntroductionDuring the last decade, a multitude of novel quantitative and semiquantitative MRI techniques have provided new information about the pathophysiology of neurological diseases. Yet, selection of the most relevant contrasts for a given pathology remains challenging. In this work, we developed and validated a method, Gated-Attention MEchanism Ranking of multi-contrast MRI in brain pathology (GAMER MRI), to rank the relative importance of MR measures in the classification of well understood ischemic stroke lesions. Subsequently, we applied this method to the classification of multiple sclerosis (MS) lesions, where the relative importance of MR measures is less understood.MethodsGAMER MRI was developed based on the gated attention mechanism, which computes attention weights (AWs) as proxies of importance of hidden features in the classification. In the first two experiments, we used Trace-weighted (Trace), apparent diffusion coefficient (ADC), Fluid-Attenuated Inversion Recovery (FLAIR), and T1-weighted (T1w) images acquired in 904 acute/subacute ischemic stroke patients and in 6,230 healthy controls and patients with other brain pathologies to assess if GAMER MRI could produce clinically meaningful importance orders in two different classification scenarios. In the first experiment, GAMER MRI with a pretrained convolutional neural network (CNN) was used in conjunction with Trace, ADC, and FLAIR to distinguish patients with ischemic stroke from those with other pathologies and healthy controls. In the second experiment, GAMER MRI with a patch-based CNN used Trace, ADC and T1w to differentiate acute ischemic stroke lesions from healthy tissue. The last experiment explored the performance of patch-based CNN with GAMER MRI in ranking the importance of quantitative MRI measures to distinguish two groups of lesions with different pathological characteristics and unknown quantitative MR features. Specifically, GAMER MRI was applied to assess the relative importance of the myelin water fraction (MWF), quantitative susceptibility mapping (QSM), T1 relaxometry map (qT1), and neurite density index (NDI) in distinguishing 750 juxtacortical lesions from 242 periventricular lesions in 47 MS patients. Pair-wise permutation t-tests were used to evaluate the differences between the AWs obtained for each quantitative measure.ResultsIn the first experiment, we achieved a mean test AUC of 0.881 and the obtained AWs of FLAIR and the sum of AWs of Trace and ADC were 0.11 and 0.89, respectively, as expected based on previous knowledge. In the second experiment, we achieved a mean test F1 score of 0.895 and a mean AW of Trace = 0.49, of ADC = 0.28, and of T1w = 0.23, thereby confirming the findings of the first experiment. In the third experiment, MS lesion classification achieved test balanced accuracy = 0.777, sensitivity = 0.739, and specificity = 0.814. The mean AWs of T1map, MWF, NDI, and QSM were 0.29, 0.26, 0.24, and 0.22 (p

Pathological data showed focal inflammation and regions of diffuse neuronal loss in the cortex of people with multiple sclerosis (MS). In this work, we applied a novel model ("soma and neurite density imaging (SANDI)") to multishell diffusion-weighted MRI data acquired in healthy subjects and people with multiple sclerosis (pwMS), in order to investigate inflammation and degeneration-related changes in the cortical tissue of pwMS. We aimed to (i) establish whether SANDI is applicable in vivo clinical data; (ii) investigate inflammatory and degenerative changes using SANDI soma fraction (fsoma)—a marker of cellularity—in both cortical lesions and in the normal-appearing-cortex and (iii) correlate SANDI fsoma with clinical and biological measures in pwMS. We applied a simplified version of SANDI to a clinical scanners. We then provided evidence that pwMS exhibited an overall decrease in cortical SANDI fsoma compared to healthy subjects, suggesting global degenerative processes compatible with neuronal loss. On the other hand, we have found that progressive pwMS showed a higher SANDI fsoma in the outer part of the cortex compared to relapsing–remitting pwMS, possibly supporting current pathological knowledge of increased innate inflammatory cells in these regions. A similar finding was obtained in subpial lesions in relapsing–remitting patients, reflecting existing pathological data in these lesion types. A significant correlation was found between SANDI fsoma and serum neurofilament light chain—a biomarker of inflammatory axonal damage—suggesting a relationship between SANDI soma fraction and inflammatory processes in pwMS again. Overall, our data show that SANDI fsoma is a promising biomarker to monitor changes in cellularity compatible with neurodegeneration and neuroinflammation in the cortex of MS patients. [ABSTRACT FROM AUTHOR]

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 myelinweighted 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. [ABSTRACT FROM AUTHOR]

Axon radius is a potential biomarker for brain diseases and a crucial tissue microstructure parameter that determines the speed of action potentials. Diffusion MRI (dMRI) allows non-invasive estimation of axon radius, but accurately estimating the radius of axons in the human brain is challenging. Most axons in the brain have a radius below one micrometer, which falls below the sensitivity limit of dMRI signals even when using the most advanced human MRI scanners. Therefore, new MRI methods that are sensitive to small axon radii are needed. In this proof-of-concept investigation, we examine whether a surface-based axonal relaxation process could mediate a relationship between intra-axonal T2 and T1 times and inner axon radius, as measured using postmortem histology. A unique in vivo human diffusion-T1-T2 relaxation dataset was acquired on a 3T MRI scanner with ultra-strong diffusion gradients, using a strong diffusion-weighting (i.e., b = 6,000 s/mm²) and multiple inversion and echo times. A second reduced diffusion-T2 dataset was collected at various echo times to evaluate the model further. The intra-axonal relaxation times were estimated by fitting a diffusionrelaxation model to the orientation-averaged spherical mean signals. Our analysis revealed that the proposed surface-based relaxation model effectively explains the relationship between the estimated relaxation times and the histological axon radius measured in various corpus callosum regions. Using these histological values, we developed a novel calibration approach to predict axon radius in other areas of the corpus callosum. Notably, the predicted radii and those determined from histological measurements were in close agreement. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Quantitative MRI provides biophysical measures of the microstructural integrity of the CNS, which can be compared across CNS regions, patients, and centres. In patients with multiple sclerosis, quantitative MRI techniques such as relaxometry, myelin imaging, magnetization transfer, diffusion MRI, quantitative susceptibility mapping, and perfusion MRI, complement conventional MRI techniques by providing insight into disease mechanisms. These include: (i) presence and extent of diffuse damage in CNS tissue outside lesions (normal-appearing tissue); (ii) heterogeneity of damage and repair in focal lesions; and (iii) specific damage to CNS tissue components. This review summarizes recent technical advances in quantitative MRI, existing pathological validation of quantitative MRI techniques, and emerging applications of quantitative MRI to patients with multiple sclerosis in both research and clinical settings. The current level of clinical maturity of each quantitative MRI technique, especially regarding its integration into clinical routine, is discussed. We aim to provide a better understanding of how quantitative MRI may help clinical practice by improving stratification of patients with multiple sclerosis, and assessment of disease progression, and evaluation of treatment response. [ABSTRACT FROM AUTHOR]

• Projecting the voxel myelin index on bundles using tractometry provides biased results when more than one bundle crosses a voxel. • Microstructure informed tractography can deconvolve myelin markers derived from voxel-wise maps on bundles to provide Bundle-specific myelin fractions (BMFs). • Myelin streamline decomposition (MySD) derived values projected on cortex show a pattern in line with literature. • BMFs estimated with MySD are consistent across two different myelin sensitive maps. 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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]