Quantification of voxel-wise total fibre density: Investigating the problems associated with track-count mapping

MVComp toolbox: MultiVariate Comparisons of brain MRI features accounting for common information across metrics

Multivariate approaches have recently gained in popularity to address the physiological unspecificity of neuroimaging metrics and to better characterize the complexity of biological processes underlying behavior. However, commonly used approaches are biased by the intrinsic associations between variables, or they are computationally expensive and may be more complicated to implement than standard univariate approaches. Here, we propose using the Mahalanobis distance (D2), an individual-level measure of deviation relative to a reference distribution that accounts for covariance between metrics. To facilitate its use, we introduce an open-source python-based tool for computing D2 relative to a reference group or within a single individual: the MultiVariate Comparison (MVComp) toolbox. The toolbox allows different levels of analysis (i.e., group- or subject-level), resolutions (e.g., voxel-wise, ROI-wise) and dimensions considered (e.g., combining MRI metrics or WM tracts). Several example cases are presented to showcase the wide range of possible applications of MVComp and to demonstrate the functionality of the toolbox. The D2 framework was applied to the assessment of white matter (WM) microstructure at 1) the group-level, where D2 can be computed between a subject and a reference group to yield an individualized measure of deviation. We observed that clustering applied to D2 in the corpus callosum yields parcellations that highly resemble known topography based on neuroanatomy, suggesting that D2 provides an integrative index that meaningfully reflects the underlying microstructure. 2) At the subject level, D2 was computed between voxels to obtain a measure of (dis)similarity. The loadings of each MRI metric (i.e., its relative contribution to D2) were then extracted in voxels of interest to showcase a useful option of the MVComp toolbox. These relative contributions can provide important insights into the physiological underpinnings of differences observed. Integrative multivariate models are crucial to expand our understanding of the complex brain-behavior relationships and the multiple factors underlying disease development and progression. Our toolbox facilitates the implementation of a useful multivariate method, making it more widely accessible.

Tremor suppression following treatment with MRgFUS: skull density ratio consistency and degree of posterior dentatorubrothalamic tract lesioning predicts long-term clinical outcomes in essential tremor

Objectives

Magnetic resonance-guided focussed ultrasound (MRgFUS) is an incisionless ablative procedure, widely used for treatment of Parkinsonian and Essential Tremor (ET). Enhanced understanding of the patient- and treatment-specific factors that influence sustained long-term tremor suppression could help clinicians achieve superior outcomes via improved patient screening and treatment strategy.

Methods

We retrospectively analysed data from 31 subjects with ET, treated with MRgFUS at a single centre. Tremor severity was assessed with parts A, B and C of the Clinical Rating Scale for Tremor (CRST) as well as the combined CRST. Tremor in the dominant and non-dominant hand was assessed with Hand Tremor Scores (HTS), derived from the CRST. Pre- and post-treatment imaging data were analysed to determine ablation volume overlap with automated thalamic segmentations, and the dentatorubrothalamic tract (DRTT) and compared with percentage change in CRST and HTS following treatment.

Results

Tremor symptoms were significantly reduced following treatment. Combined pre-treatment CRST (mean: 60.7 ± 17.3) and HTS (mean: 19.2 ± 5.7) improved by an average of 45.5 and 62.6%, respectively. Percentage change in CRST was found to be significantly negatively associated with age (β = -0.375, p = 0.015), and SDR standard deviation (SDRSD; β = -0.324, p = 0.006), and positively associated with ablation overlap with the posterior DRTT (β = 0.535, p < 0.001). Percentage HTS improvement in the dominant hand decreased significantly with older age (β = -0.576, p < 0.01).

Conclusion

Our results suggest that increased lesioning of the posterior region of the DRTT could result in greater improvements in combined CRST and non-dominant hand HTS, and that subjects with lower SDR standard deviation tended to experience greater improvement in combined CRST.

Investigating Tissue-Specific Abnormalities in Alzheimer's Disease with Multi-Shell Diffusion MRI

BACKGROUND

Most studies using diffusion-weighted MRI (DW-MRI) in Alzheimer's disease (AD) have focused their analyses on white matter (WM) microstructural changes using the diffusion (kurtosis) tensor model. Although recent works have addressed some limitations of the tensor model, such as the representation of crossing fibers and partial volume effects with cerebrospinal fluid (CSF), the focus remains in modeling and analyzing the WM.

OBJECTIVE

In this work, we present a brain analysis approach for DW-MRI that disentangles multiple tissue compartments as well as micro- and macroscopic effects to investigate differences between groups of subjects in the AD continuum and controls.

METHODS

By means of the multi-tissue constrained spherical deconvolution of multi-shell DW-MRI, underlying brain tissue is modeled with a WM fiber orientation distribution function along with the contributions of gray matter (GM) and CSF to the diffusion signal. From this multi-tissue model, a set of measures capturing tissue diffusivity properties and morphology are extracted. Group differences were interrogated following fixel-, voxel-, and tensor-based morphometry approaches while including strong FWE control across multiple comparisons.

RESULTS

Abnormalities related to AD stages were detected in WM tracts including the splenium, cingulum, longitudinal fasciculi, and corticospinal tract. Changes in tissue composition were identified, particularly in the medial temporal lobe and superior longitudinal fasciculus.

CONCLUSION

This analysis framework constitutes a comprehensive approach allowing simultaneous macro and microscopic assessment of WM, GM, and CSF, from a single DW-MRI dataset.

Gradual changes in microarchitectural properties of cortex and juxtacortical white matter: Observed by anatomical and diffusion MRI

PURPOSE

Characterization of cerebral cortex is challenged by the complexity and heterogeneity of its cyto- and myeloarchitecture. This study evaluates quantitative MRI metrics, measured across two cortical depths and in subcortical white matter (WM) adjacent to cortex (juxtacortical WM), indicative of myelin content, neurite density, and diffusion microenvironment, for a comprehensive characterization of cortical microarchitecture.

METHODS

High-quality structural and diffusion MRI data (N = 30) from the Human Connectome Project were processed to compute myelin index, neurite density index, fractional anisotropy, mean diffusivity, axial diffusivity, and radial diffusivity from superficial cortex, deep cortex, and juxtacortical WM. The distributional patterns of these metrics were analyzed individually, correlated to one another, and were compared to established parcellations.

RESULTS

Our results supported that myeloarchitectonic and the coexisting cytoarchitectonic structures influence the diffusion properties of water molecules residing in cortex. Full cortical thickness showed myelination patterns similar to those previously observed in humans. Higher myelin indices with similar distributional patterns were observed in deep cortex whereas lower myelin indices were observed in superficial cortex. Neurite density index and other diffusion MRI derived parameters provided complementary information to myelination. Reliable and reproducible correlations were identified among the cortical microarchitectural properties and fiber distributional patterns in proximal WM structures.

CONCLUSION

We demonstrated gradual changes across the cortical sheath by assessing depth-specific cortical micro-architecture using anatomical and diffusion MRI. Mutually independent but coexisting features of cortical layers and juxtacortical WM provided new insights towards structural organizational units and variabilities across cortical regions and through depth.

Leveraging manifold learning techniques to explore white matter anomalies: An application of the TractLearn pipeline in epilepsy

An accurate description of brain white matter anatomy in vivo remains a challenge. However, technical progress allows us to analyze structural variations in an increasingly sophisticated way. Current methods of processing diffusion MRI data now make it possible to correct some limiting biases. In addition, the development of statistical learning algorithms offers the opportunity to analyze the data from a new perspective. We applied newly developed tractography models to extract quantitative white matter parameters in a group of patients with chronic temporal lobe epilepsy. Furthermore, we implemented a statistical learning workflow optimized for the MRI diffusion data - the TractLearn pipeline - to model inter-individual variability and predict structural changes in patients. Finally, we interpreted white matter abnormalities in the context of several other parameters reflecting clinical status, as well as neuronal and cognitive functioning for these patients. Overall, we show the relevance of such a diffusion data processing pipeline for the evaluation of clinical populations. The "global to fine scale" funnel statistical approach proposed in this study also contributes to the understanding of neuroplasticity mechanisms involved in refractory epilepsy, thus enriching previous findings.

MRI-Visible Anatomy of the Brainstem

Human brainstem internal anatomy is intricate, complex, and essential to normal brain function. The brainstem is affected by stroke, multiple sclerosis, and most neurodegenerative diseases-a 1-mm focus of pathologic condition can have profound clinical consequences. Unfortunately, detailed internal brainstem anatomy is difficult to see with conventional MRI sequences. We review normal brainstem anatomy visualized on widely available clinical 3-T MRI scanners using fast gray matter acquisition T1 inversion recovery, probabilistic diffusion tractography, neuromelanin, and susceptibility-weighted imaging. Better anatomic localization using these recent innovations improves our ability to diagnose, localize, and treat brainstem diseases. We aim to provide an accessible review of the most clinically relevant brainstem neuroanatomy.

MRI-Visible Anatomy of the Basal Ganglia and Thalamus

Conventional MR imaging does not discriminate basal ganglia and thalamic internal anatomy well. Radiology reports describe anatomic locations but not specific functional structures. Functional neurosurgery uses indirect targeting based on commissural coordinates or atlases that do not fully account for individual variability. We describe innovative MR imaging sequences that improve the visualization of normal anatomy in this complex brain region and may increase our understanding of basal ganglia and thalamic function. Better visualization also may improve treatments for movement disorders and other emerging functional neurosurgery targets. We aim to provide an accessible review of the most clinically-relevant neuroanatomy within the thalamus and basal ganglia.

The impact of levamisole and alcohol on white matter microstructure in adult chronic cocaine users

Previous brain imaging studies with chronic cocaine users (CU) using diffusion tensor imaging (DTI) mostly focused on fractional anisotropy to investigate white matter (WM) integrity. However, a quantitative interpretation of fractional anisotropy (FA) alterations is often impeded by the inherent limitations of the underlying tensor model. A more fine-grained measure of WM alterations could be achieved by measuring fibre density (FD). This study investigates this novel DTI metric comparing 23 chronic CU and 32 healthy subjects. Quantitative hair analysis was used to determine intensity of cocaine and levamisole exposure-a cocaine adulterant with putative WM neurotoxicity. We first assessed the impact of cocaine use, levamisole exposure and alcohol use on group differences in WM integrity. Compared with healthy controls, all models revealed cortical reductions of FA and FD in CU. At the within-patient group level, we found that alcohol use and levamisole exposure exhibited regionally different FA and FD alterations than cocaine use. We found mostly negative correlations of tract-based WM associated with levamisole and weekly alcohol use. Specifically, levamisole exposure was linked with stronger WM reductions in the corpus callosum than alcohol use. Cocaine use duration correlated negatively with FA and FD in some regions. Yet, most of these correlations did not survive a correction for multiple testing. Our results suggest that chronic cocaine use, levamisole exposure and alcohol use were all linked to significant WM impairments in CU. We conclude that FD could be a sensitive marker to detect the impact of the use of multiple substances on WM integrity in cocaine but also other substance use disorders.

Structural brain connectivity correlates with outcome in mild traumatic brain injury

We investigated the topology of structural brain connectivity networks and its association to outcome following mild traumatic brain injury, a major cause of permanent disability. Eighty-five patients with mild traumatic brain injury underwent MRI twice, about three weeks and eight months after injury, and 30 age-matched orthopedic trauma control subjects were scanned. Outcome was assessed with Extended Glasgow Outcome Scale on average eight months after injury. We performed constrained spherical deconvolution based probabilistic streamlines tractography on diffusion MRI data and parcellated cortical and subcortical gray matter into 84 regions based on T1-weighted data to reconstruct structural brain connectivity networks weighted by the number of streamlines. Graph theoretical methods were employed to measure network properties in both patients and controls, and correlations between these properties and outcome were calculated. We found no global differences in the network properties between patients with mild traumatic brain injury and orthopedic control subjects at either stage. However, we found significantly increased betweenness centrality of the right pars opercularis in the chronic stage compared to control subjects. Furthermore, both global and local network properties correlated significantly with outcome. Higher normalized global efficiency, degree, and strength as well as lower small-worldness were associated with better outcome. Correlations between the outcome and the local network properties were the most prominent in the left putamen and the left postcentral gyrus. Our results indicate that both global and local network properties provide valuable information about the outcome already in the acute/subacute stage, and therefore, are promising biomarkers for prognostic purposes in mild traumatic brain injury.

TractLearn: A geodesic learning framework for quantitative analysis of brain bundles

Deep learning-based convolutional neural networks have recently proved their efficiency in providing fast segmentation of major brain fascicles structures, based on diffusion-weighted imaging. The quantitative analysis of brain fascicles then relies on metrics either coming from the tractography process itself or from each voxel along the bundle. Statistical detection of abnormal voxels in the context of disease usually relies on univariate and multivariate statistics models, such as the General Linear Model (GLM). Yet in the case of high-dimensional low sample size data, the GLM often implies high standard deviation range in controls due to anatomical variability, despite the commonly used smoothing process. This can lead to difficulties to detect subtle quantitative alterations from a brain bundle at the voxel scale. Here we introduce TractLearn, a unified framework for brain fascicles quantitative analyses by using geodesic learning as a data-driven learning task. TractLearn allows a mapping between the image high-dimensional domain and the reduced latent space of brain fascicles using a Riemannian approach. We illustrate the robustness of this method on a healthy population with test-retest acquisition of multi-shell diffusion MRI data, demonstrating that it is possible to separately study the global effect due to different MRI sessions from the effect of local bundle alterations. We have then tested the efficiency of our algorithm on a sample of 5 age-matched subjects referred with mild traumatic brain injury. Our contributions are to propose: 1/ A manifold approach to capture controls variability as standard reference instead of an atlas approach based on a Euclidean mean. 2/ A tool to detect global variation of voxels' quantitative values, which accounts for voxels' interactions in a structure rather than analyzing voxels independently. 3/ A ready-to-plug algorithm to highlight nonlinear variation of diffusion MRI metrics. With this regard, TractLearn is a ready-to-use algorithm for precision medicine.

Investigation of Cerebral White Matter Changes After Spinal Cord Injury With a Measure of Fiber Density

Remote neurodegenerative changes in supraspinal white matter (WM) can manifest after central lesions such as spinal cord injury (SCI). The majority of diffusion tensor imaging (DTI) studies use traditional metrics such as fractional anisotropy (FA) and mean diffusivity (MD) to investigate microstructural changes in cerebral WM after SCI. However, interpretation of FA readouts is often challenged by inherent limitations of the tensor model. Recent developments in novel diffusion markers, such as fiber density (FD), allows more accurate depictions of WM pathways and has shown more reliable quantification of WM alterations compared to FA in recent studies of neurological diseases. This study investigated if FD provides useful characterization of supraspinal WM integrity after SCI in addition to the traditional DTI readouts. FA, MD, and FD maps were derived from diffusion datasets of 20 patients with chronic SCI and compared with 19 healthy controls (HC). Group differences were investigated across whole brain WM using tract-based spatial statistics and averaged diffusion values of the corticospinal tract (CST) and thalamic radiation (TR) were extracted for comparisons between HC and SCI subgroups. We also related diffusion readouts of the CST and TR with clinical scores of sensorimotor function. To investigate which diffusion markers of the CST and TR delineate HC and patients with SCI a receiver operating characteristic (ROC) analysis was performed. Overall, patients with an SCI showed decreased FA of the TR and CST. ROC analysis differentiated HC and SCI based on diffusion markers of large WM tracts including FD of the TR. Furthermore, patients' motor function was positively correlated with greater microstructural integrity of the CST. While FD showed the strongest correlation, motor function was also associated with FA and MD of the CST. In summary, microstructural changes of supraspinal WM in patients with SCI can be detected using FD as a complementary marker to traditional DTI readouts and correlates with their clinical characteristics. Future DTI studies may benefit from utilizing this novel marker to investigate complex large WM tracts in patient cohorts with varying presentations of SCI or neurodegenerative diseases.

Reduced apparent fiber density in the white matter of premature-born adults

Premature-born adults exhibit lasting white matter alterations as demonstrated by widespread reduction in fractional anisotropy (FA) based on diffusion-weighted imaging (DWI). FA reduction, however, is non-specific for microscopic underpinnings such as aberrant myelination or fiber density (FD). Using recent advances in DWI, we tested the hypothesis of reduced FD in premature-born adults and investigated its link with the degree of prematurity and cognition. 73 premature- and 89 mature-born adults aged 25-27 years underwent single-shell DWI, from which a FD measure was derived using convex optimization modeling for microstructure informed tractography (COMMIT). Premature-born adults exhibited lower FD in numerous tracts including the corpus callosum and corona radiata compared to mature-born adults. These FD alterations were associated with both the degree of prematurity, as assessed via gestational age and birth weight, as well as with reduced cognition as measured by full-scale IQ. Finally, lower FD overlapped with lower FA, suggesting lower FD underlie unspecific FA reductions. Results provide evidence that premature birth leads to lower FD in adulthood which links with lower full-scale IQ. Data suggest that lower FD partly underpins FA reductions of premature birth but that other processes such as hypomyelination might also take place.

Technical and radiographic considerations for magnetic resonance imaging-guided focused ultrasound capsulotomy

OBJECTIVE

Magnetic resonance imaging-guided focused ultrasound (MRgFUS) is an emerging treatment modality that enables incisionless ablative neurosurgical procedures. Bilateral MRgFUS capsulotomy has recently been demonstrated to be safe and effective in treating obsessive-compulsive disorder (OCD) and major depressive disorder (MDD). Preliminary evidence has suggested that bilateral MRgFUS capsulotomy can present increased difficulties in reaching lesional temperatures as compared to unilateral thalamotomy. The authors of this article aimed to study the parameters associated with successful MRgFUS capsulotomy lesioning and to present longitudinal radiographic findings following MRgFUS capsulotomy.

METHODS

Using data from 22 attempted MRgFUS capsulotomy treatments, the authors investigated the relationship between various sonication parameters and the maximal temperature achieved at the intracranial target. Lesion volume and morphology were analyzed longitudinally using structural and diffusion tensor imaging. A retreatment procedure was attempted in one patient, and their postoperative imaging is presented.

RESULTS

Skull density ratio (SDR), skull thickness, and angle of incidence were significantly correlated with the maximal temperature achieved. MRgFUS capsulotomy lesions appeared similar to those following MRgFUS thalamotomy, with three concentric zones observed on MRI. Lesion volumes regressed substantially over time following MRgFUS. Fractional anisotropy analysis revealed a disruption in white matter integrity, followed by a gradual return to near-baseline levels concurrent with lesion regression. In the patient who underwent retreatment, successful bilateral lesioning was achieved, and there were no adverse clinical or radiographic events.

CONCLUSIONS

With the current iteration of MRgFUS technology, skull-related parameters such as SDR, skull thickness, and angle of incidence should be considered when selecting patients suitable for MRgFUS capsulotomy. Lesions appear to follow morphological patterns similar to what is seen following MRgFUS thalamotomy. Retreatment appears to be safe, although additional cases will be necessary to further evaluate the associated safety profile.

The Seven Deadly Sins of Measuring Brain Structural Connectivity Using Diffusion MRI Streamlines Fibre-Tracking

There is great interest in the study of brain structural connectivity, as white matter abnormalities have been implicated in many disease states. Diffusion magnetic resonance imaging (MRI) provides a powerful means to characterise structural connectivity non-invasively, by using a fibre-tracking algorithm. The most widely used fibre-tracking strategy is based on the step-wise generation of streamlines. Despite their popularity and widespread use, there are a number of practical considerations that must be taken into account in order to increase the robustness of streamlines tracking results, particularly when these methods are used to study brain structural connectivity, and the connectome. This review article describes what we consider the ‘seven deadly sins’ of mapping structural connections using diffusion MRI streamlines fibre-tracking, with particular emphasis on ‘sins’ that can be practically avoided and they can have an important impact in the results. It is shown that there are important ‘deadly sins’ to be avoided at every step of the pipeline, such as during data acquisition, during data modelling to estimate local fibre architecture, during the fibre-tracking process itself, and during quantification of the tracking results. The recommendations here are intended to inform users on potential important shortcomings of their current tracking protocols, as well as to guide future users on some of the key issues and decisions that must be faced when designing their processing pipelines.

Diffusion MRI fiber tractography of the brain

The ability of fiber tractography to delineate non-invasively the white matter fiber pathways of the brain raises possibilities for clinical applications and offers enormous potential for neuroscience. In the last decade, fiber tracking has become the method of choice to investigate quantitative MRI parameters in specific bundles of white matter. For neurosurgeons, it is quickly becoming an invaluable tool for the planning of surgery, allowing for visualization and localization of important white matter pathways before and even during surgery. Fiber tracking has also claimed a central role in the field of "connectomics," a technique that builds and studies comprehensive maps of the complex network of connections within the brain, and to which significant resources have been allocated worldwide. Despite its unique abilities and exciting applications, fiber tracking is not without controversy, in particular when it comes to its interpretation. As neuroscientists are eager to study the brain's connectivity, the quantification of tractography-derived "connection strengths" between distant brain regions is becoming increasingly popular. However, this practice is often frowned upon by fiber-tracking experts. In light of this controversy, this paper provides an overview of the key concepts of tractography, the technical considerations at play, and the different types of tractography algorithm, as well as the common misconceptions and mistakes that surround them. We also highlight the ongoing challenges related to fiber tracking. While recent methodological developments have vastly increased the biological accuracy of fiber tractograms, one should be aware that, even with state-of-the-art techniques, many issues that severely bias the resulting structural "connectomes" remain unresolved.

Modelling white matter with spherical deconvolution: How and why?

Since the realization that diffusion MRI can probe the microstructural organization and orientation of biological tissue in vivo and non-invasively, a multitude of diffusion imaging methods have been developed and applied to study the living human brain. Diffusion tensor imaging was the first model to be widely adopted in clinical and neuroscience research, but it was also clear from the beginning that it suffered from limitations when mapping complex configurations, such as crossing fibres. In this review, we highlight the main steps that have led the field of diffusion imaging to move from the tensor model to the adoption of diffusion and fibre orientation density functions as a more effective way to describe the complexity of white matter organization within each brain voxel. Among several techniques, spherical deconvolution has emerged today as one of the main approaches to model multiple fibre orientations and for tractography applications. Here we illustrate the main concepts and the reasoning behind this technique, as well as the latest developments in the field. The final part of this review provides practical guidelines and recommendations on how to set up processing and acquisition protocols suitable for spherical deconvolution.

Subtle white matter alterations in schizophrenia identified with a new measure of fiber density

Altered cerebral connectivity is one of the core pathophysiological mechanism underlying the development and progression of information-processing deficits in schizophrenia. To date, most diffusion tensor imaging (DTI) studies used fractional anisotropy (FA) to investigate disrupted white matter connections. However, a quantitative interpretation of FA changes is often impeded by the inherent limitations of the underlying tensor model. A more fine-grained measure of white matter alterations could be achieved by measuring fiber density (FD) - a novel non-tensor-derived diffusion marker. This study investigates, for the first time, FD alterations in schizophrenia patients. FD and FA maps were derived from diffusion data of 25 healthy controls (HC) and 21 patients with schizophrenia (SZ). Using tract-based spatial statistics (TBSS), group differences in FD and FA were investigated across the entire white matter. Furthermore, we performed a region of interest (ROI) analysis of frontal fasciculi to detect potential correlations between FD and positive symptoms. As a result, whole brain TBSS analysis revealed reduced FD in SZ patients compared to HC in several white matter tracts including the left and right thalamic radiation (TR), superior longitudinal fasciculus (SLF), corpus callosum (CC), and corticospinal tract (CST). In contrast, there were no significant FA differences between groups. Further, FD values in the TR were negatively correlated with the severity of positive symptoms and medication dose in SZ patients. In summary, a novel diffusion-weighted data analysis approach enabled us to identify widespread FD changes in SZ patients with most prominent white matter alterations in the frontal and subcortical regions. Our findings suggest that the new FD measure may be more sensitive to subtle changes in the white matter microstructure compared to FA, particularly in the given population. Therefore, investigating FD may be a promising approach to detect subtle changes in the white matter microstructure of altered connectivity in schizophrenia.

Reproducibility of multi-shell diffusion tractography on traveling subjects: A multicenter study prospective

Reproducibility of multicenter diffusion magnetic resonance imaging has drawn more attention recently due to rapidly increasing need for large-size brain imaging studies. Advanced multi-shell diffusion models are recommended for their potentials to provide variety of physio-pathological information. While previous studies have investigated the consistency of single-shell diffusion acquisition from various hardware and protocols, a well-controlled study with multi-shell acquisition would be necessary to understand the inherent factors of reproducibility from new complexity of such acquisition protocol. In this study, three traveling subjects were scanned at eight imaging centers equipped with the same type of scanners using the same multi-shell diffusion imaging protocol. Track density imaging and structure connectomes were investigated in local-scale distribution and in distal-scale connectivity, respectively. With evaluations of the coefficient of variation and the intra-class correlation coefficient, our results indicated: 1) similar to single-shell schemes, the intra-center reproducibility of multi-shell is higher than inter-center; 2) multi-shell schemes produce higher reproducibility and precision among centers compared to the single-shell schemes; and 3) in addition to the diffusion schemes, image quality and the presence of complex fiber structure could also associated with multicenter reproducibility.

Multivariate pattern classification of brain white matter connectivity predicts classic trigeminal neuralgia

Trigeminal neuralgia (TN) is a severe form of chronic facial neuropathic pain. Increasing interest in the neuroimaging of pain has highlighted changes in the root entry zone in TN, but also group-level central nervous system gray and white matter (WM) abnormalities. Group differences in neuroimaging data are frequently evaluated with univariate statistics; however, this approach is limited because it is based on single, or clusters of, voxels. By contrast, multivariate pattern analyses consider all the model's neuroanatomical features to capture a specific distributed spatial pattern. This approach has potential use as a prediction tool at the individual level. We hypothesized that a multivariate pattern classification method can distinguish specific patterns of abnormal WM connectivity of classic TN from healthy controls (HCs). Diffusion-weighted scans in 23 right-sided TN and matched controls were processed to extract whole-brain interregional streamlines. We used a linear support vector machine algorithm to differentiate interregional normalized streamline count between TN and HC. This algorithm successfully differentiated between TN and HC with an accuracy of 88%. The structural pattern emphasized WM connectivity of regions that subserve sensory, affective, and cognitive dimensions of pain, including the insula, precuneus, inferior and superior parietal lobules, and inferior and medial orbital frontal gyri. Normalized streamline counts were associated with longer pain duration and WM metric abnormality between the connections. This study demonstrates that machine-learning algorithms can detect characteristic patterns of structural alterations in TN and highlights the role of structural brain imaging for identification of neuroanatomical features associated with neuropathic pain disorders.

Bundle-specific tractography with incorporated anatomical and orientational priors

Anatomical white matter bundles vary in shape, size, length, and complexity, making diffusion MRI tractography reconstruction of some bundles more difficult than others. As a result, bundles reconstruction often suffers from a poor spatial extent recovery. To fill-up the white matter volume as much and as best as possible, millions of streamlines can be generated and filtering techniques applied to address this issue. However, well-known problems and biases are introduced such as the creation of a large number of false positives and over-representation of easy-to-track parts of bundles and under-representation of hard-to-track. To address these challenges, we developed a Bundle-Specific Tractography (BST) algorithm. It incorporates anatomical and orientational prior knowledge during the process of streamline tracing to increase reproducibility, sensitivity, specificity and efficiency when reconstructing certain bundles of interest. BST outperforms classical deterministic, probabilistic, and global tractography methods. The increase in anatomically plausible streamlines, with larger spatial coverage, helps to accurately represent the full shape of bundles, which could greatly enhance and robustify tract-based and connectivity-based neuroimaging studies.

Linking spherical mean diffusion weighted signal with intra-axonal volume fraction

Diffusion MRI has been widely used to assess brain tissue microstructure. However, the conventional diffusion tensor imaging (DTI) is inadequate for characterizing fiber direction or fiber density in voxels with crossing fibers in brain white matter. The constrained spherical deconvolution (CSD) technique has been proposed to measure the complex fiber orientation distribution (FOD) using a single high b-value (b ≥ 3000 s/mm²) to derive the intra-axonal volume fraction (V_in) from the calculated FOD. Recently, the spherical mean technique (SMT) was developed to fit V_in directly from a multi-compartment model with multi-shell b-values. Although different numbers of b-values are needed in the two techniques, both methods have been suggested to be related to the spherical mean diffusion weighted signal (S¯). The current study compared the two techniques on the same high-quality Human Connectome Project diffusion data and investigated the relation between S¯ and V_in systematically. At high b-values (b ≥ 3000 s/mm²), S¯ is linearly related to V_in, and S¯ provides similar contrast with V_in in white matter. At low b-values (b ~ 1000 s/mm²), the linear relation between S¯ and V_in is sensitive to the variations of intrinsic diffusivity. These results demonstrate that S¯ measured with the typical b-value of 1000 s/mm² is not an indicator of V_in, and previous DTI studies acquired with b = 1000 s/mm² cannot be re-analyzed to provide V_in-weighted contrast.

Minimal number of gradient directions for robust measurement of spherical mean diffusion weighted signal

PURPOSE

Determination of the minimum number of gradient directions (N_min) for robust measurement of spherical mean diffusion weighted signal (S¯).

METHODS

Computer simulations were employed to characterize the relative standard deviation (RSD) of the measured spherical mean signal as a function of the number of gradient directions (N). The effects of diffusion weighting b-value and signal-to-noise ratio (SNR) were investigated. Multi-shell high angular resolution Human Connectome Project diffusion data were analyzed to support the simulation results.

RESULTS

RSD decreases with increasing N, and the minimum number of N needed for RSD ≤ 5% is referred to as N_min. At high SNRs, N_min increases with increasing b-value to achieve sufficient sampling. Simulations showed that N_min is linearly dependent on the b-value. At low SNRs, N_min increases with increasing b-value to reduce the noise. RSD can be estimated as σS¯N, where σ = 1/SNR is the noise level. The experimental results were in good agreement with the simulation results. The spherical mean signal can be measured accurately with a subset of gradient directions.

CONCLUSION

As N_min is affected by b-value and SNR, we recommend using 10 × b / b₁ (b₁ = 1 ms/μm²) uniformly distributed gradient directions for typical human diffusion studies with SNR ~ 20 for robust spherical mean signal measurement.

White matter degeneration profile in the cognitive cortico-subcortical tracts in Parkinson's disease

BACKGROUND

In Parkinson's disease cognitive impairment is an early nonmotor feature, but it is still unclear why some patients are able to maintain their cognitive performance at normal levels, as quantified by neuropsychological tests, whereas others cannot. The objectives of this study were to perform a cross-sectional study and analyze the white matter changes in the cognitive and motor bundles in patients with Parkinson's disease.

METHODS

Sixteen Parkinson's disease patients with normal cognitive performance, 19 with mild cognitive impairment (based on their performance of 1.5 standard deviations below the healthy population mean), and 16 healthy controls were compared with respect to their tractography patterns between the cortical cognitive / motor regions and subcortical structures, using high angular resolution diffusion imaging and constrained spherical deconvolution computation.

RESULTS

Motor bundles showed decreased apparent fiber density in both PD groups, associated with a significant increase in diffusivity metrics, number of reconstructed streamlines, and track volumes, compared with healthy controls. By contrast, in the cognitive bundles, decreased fiber density in both Parkinson's groups was compounded by the absence of changes in diffusivity in patients with normal cognition, whereas patients with cognitive impairment had increased diffusivity metrics, lower numbers of reconstructed streamlines, and lower track volumes.

CONCLUSIONS

Both PD groups showed similar patterns of white matter neurodegeneration in the motor bundles, whereas cognitive bundles showed a distinct pattern: Parkinson's patients with normal cognition had white matter diffusivity metrics similar to healthy controls, whereas in patients with cognitive impairment white matter showed a neurodegeneration pattern. © 2018 International Parkinson and Movement Disorder Society.

Uniformity and Deviation of Intra-axonal Cross-sectional Area Coverage of the Gray-to-White Matter Interface

Diffusion magnetic resonance imaging (dMRI) is a compelling tool for investigating the structure and geometry of brain tissue based on indirect measurement of the diffusion anisotropy of water. Recent developments in global top-down tractogram optimizations enable the estimation of streamline weights, which characterize the connection between gray matter areas. In this work, the intra-axonal cross-sectional area coverage of the gray-to-white matter interface was examined by intersecting tractography streamlines with cortical regions of interest. The area coverage is the ratio of streamline weights divided by the surface area at the gray-to-white matter interface and assesses the estimated percentage which is covered by intra-axonal space. A high correlation (r = 0.935) between streamline weights and the cortical surface area was found across all regions of interest in all subjects. The variance across different cortical regions exhibits similarities to myelin maps. Additionally, we examined the effect of different diffusion gradient subsets at a lower, clinically feasible spatial resolution. Subsampling of the initial high-resolution diffusion dataset did not alter the tendency of the area coverage at the gray-to-white matter interface across cortical areas and subjects. However, single-shell acquisition schemes with lower b-values lead to a steady increase in area coverage in comparison to the full acquisition scheme at high resolution.

Super-Resolution Track-Density Imaging Reveals Fine Anatomical Features in Tree Shrew Primary Visual Cortex and Hippocampus

Diffusion-weighted magnetic resonance imaging (dMRI) is widely used to study white and gray matter (GM) micro-organization and structural connectivity in the brain. Super-resolution track-density imaging (TDI) is an image reconstruction method for dMRI data, which is capable of providing spatial resolution beyond the acquired data, as well as novel and meaningful anatomical contrast that cannot be obtained with conventional reconstruction methods. TDI has been used to reveal anatomical features in human and animal brains. In this study, we used short track TDI (stTDI), a variation of TDI with enhanced contrast for GM structures, to reconstruct direction-encoded color maps of fixed tree shrew brain. The results were compared with those obtained with the traditional diffusion tensor imaging (DTI) method. We demonstrated that fine microstructures in the tree shrew brain, such as Baillarger bands in the primary visual cortex and the longitudinal component of the mossy fibers within the hippocampal CA3 subfield, were observable with stTDI, but not with DTI reconstructions from the same dMRI data. The possible mechanisms underlying the enhanced GM contrast are discussed.

The role of whole-brain diffusion MRI as a tool for studying human in vivo cortical segregation based on a measure of neurite density

PURPOSE

To investigate whether diffusion MRI can be used to study cortical segregation based on a contrast related to neurite density, thus providing a complementary tool to myelin-based MRI techniques used for myeloarchitecture.

METHODS

Several myelin-sensitive MRI methods (e.g., based on T₁ , T₂ , and T2*) have been proposed to parcellate cortical areas based on their myeloarchitecture. Recent improvements in hardware, acquisition, and analysis methods have opened the possibility of achieving a more robust characterization of cortical microstructure using diffusion MRI. High-quality diffusion MRI data from the Human Connectome Project was combined with recent advances in fiber orientation modeling. The orientational average of the fiber orientation distribution was used as a summary parameter, which was displayed as inflated brain surface views.

RESULTS

Diffusion MRI identifies cortical patterns consistent with those previously seen by MRI methods used for studying myeloarchitecture, which have shown patterns of high myelination in the sensorimotor strip, visual cortex, and auditory areas and low myelination in frontal and anterior temporal areas.

CONCLUSION

In vivo human diffusion MRI provides a useful complementary noninvasive approach to myelin-based methods used to study whole-brain cortical parcellation, by exploiting a contrast based on tissue microstructure related to neurite density, rather than myelin itself. Magn Reson Med 79:2738-2744, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

A test-retest study on Parkinson's PPMI dataset yields statistically significant white matter fascicles

In this work, we propose a diffusion MRI protocol for mining Parkinson's disease diffusion MRI datasets and recover robust disease-specific biomarkers. Using advanced high angular resolution diffusion imaging (HARDI) crossing fiber modeling and tractography robust to partial volume effects, we automatically dissected 50 white matter (WM) fascicles. These fascicles connect deep nuclei (thalamus, putamen, pallidum) to different cortical functional areas (associative, motor, sensorimotor, limbic), basal forebrain and substantia nigra. Then, among these 50 candidate WM fascicles, only the ones that passed a test-retest reproducibility procedure qualified for further tractometry analysis. Leveraging the unique 2-timepoints test-retest Parkinson's Progression Markers Initiative (PPMI) dataset of over 600 subjects, we found statistically significant differences in tract profiles along the subcortico-cortical pathways between Parkinson's disease patients and healthy controls. In particular, significant increases in FA, apparent fiber density, tract-density and generalized FA were detected in some locations of the nigro-subthalamo-putaminal-thalamo-cortical pathway. This connection is one of the major motor circuits balancing the coordination of motor output. Detailed and quantifiable knowledge on WM fascicles in these areas is thus essential to improve the quality and outcome of Deep Brain Stimulation, and to target new WM locations for investigation.

Visualization, Interaction and Tractometry: Dealing with Millions of Streamlines from Diffusion MRI Tractography

Recently proposed tractography and connectomics approaches often require a very large number of streamlines, in the order of millions. Generating, storing and interacting with these datasets is currently quite difficult, since they require a lot of space in memory and processing time. Compression is a common approach to reduce data size. Recently such an approach has been proposed consisting in removing collinear points in the streamlines. Removing points from streamlines results in files that cannot be robustly post-processed and interacted with existing tools, which are for the most part point-based. The aim of this work is to improve visualization, interaction and tractometry algorithms to robustly handle compressed tractography datasets. Our proposed improvements are threefold: (i) An efficient loading procedure to improve visualization (reduce memory usage up to 95% for a 0.2 mm step size); (ii) interaction techniques robust to compressed tractograms; (iii) tractometry techniques robust to compressed tractograms to eliminate biased in tract-based statistics. The present work demonstrates the need of correctly handling compressed streamlines to avoid biases in future tractometry and connectomics studies.

Improvement in White Matter Tract Reconstruction with Constrained Spherical Deconvolution and Track Density Mapping in Low Angular Resolution Data: A Pediatric Study and Literature Review

INTRODUCTION

Diffusion-weighted magnetic resonance imaging (DW-MRI) allows noninvasive investigation of brain structure in vivo. Diffusion tensor imaging (DTI) is a frequently used application of DW-MRI that assumes a single main diffusion direction per voxel, and is therefore not well suited for reconstructing crossing fiber tracts. Among the solutions developed to overcome this problem, constrained spherical deconvolution with probabilistic tractography (CSD-PT) has provided superior quality results in clinical settings on adult subjects; however, it requires particular acquisition parameters and long sequences, which may limit clinical usage in the pediatric age group. The aim of this work was to compare the results of DTI with those of track density imaging (TDI) maps and CSD-PT on data from neonates and children, acquired with low angular resolution and low b-value diffusion sequences commonly used in pediatric clinical MRI examinations.

MATERIALS AND METHODS

We analyzed DW-MRI studies of 50 children (eight neonates aged 3-28 days, 20 infants aged 1-8 months, and 22 children aged 2-17 years) acquired on a 1.5 T Philips scanner using 34 gradient directions and a b-value of 1,000 s/mm². Other sequence parameters included 60 axial slices; acquisition matrix, 128 × 128; average scan time, 5:34 min; voxel size, 1.75 mm × 1.75 mm × 2 mm; one b = 0 image. For each subject, we computed principal eigenvector (EV) maps and directionally encoded color TDI maps (DEC-TDI maps) from whole-brain tractograms obtained with CSD-PT; the cerebellar-thalamic, corticopontocerebellar, and corticospinal tracts were reconstructed using both CSD-PT and DTI. Results were compared by two neuroradiologists using a 5-point qualitative score.

RESULTS

The DEC-TDI maps obtained presented higher anatomical detail than EV maps, as assessed by visual inspection. In all subjects, white matter (WM) tracts were successfully reconstructed using both tractography methodologies. The mean qualitative scores of all tracts obtained with CSD-PT were significantly higher than those obtained with DTI (p-value < 0.05 for all comparisons).

CONCLUSION

CSD-PT can be successfully applied to DW-MRI studies acquired at 1.5 T with acquisition parameters adapted for pediatric subjects, thus providing TDI maps with greater anatomical detail. This methodology yields satisfactory results for clinical purposes in the pediatric age group.

Semi-Automatic Segmentation of Optic Radiations and LGN, and Their Relationship to EEG Alpha Waves

At rest, healthy human brain activity is characterized by large electroencephalography (EEG) fluctuations in the 8-13 Hz range, commonly referred to as the alpha band. Although it is well known that EEG alpha activity varies across individuals, few studies have investigated how this may be related to underlying morphological variations in brain structure. Specifically, it is generally believed that the lateral geniculate nucleus (LGN) and its efferent fibres (optic radiation, OR) play a key role in alpha activity, yet it is unclear whether their shape or size variations contribute to its inter-subject variability. Given the widespread use of EEG alpha in basic and clinical research, addressing this is important, though difficult given the problems associated with reliably segmenting the LGN and OR. For this, we employed a multi-modal approach and combined diffusion magnetic resonance imaging (dMRI), functional magnetic resonance imaging (fMRI) and EEG in 20 healthy subjects to measure structure and function, respectively. For the former, we developed a new, semi-automated approach for segmenting the OR and LGN, from which we extracted several structural metrics such as volume, position and diffusivity. Although these measures corresponded well with known morphology based on previous post-mortem studies, we nonetheless found that their inter-subject variability was not significantly correlated to alpha power or peak frequency (p >0.05). Our results therefore suggest that alpha variability may be mediated by an alternative structural source and our proposed methodology may in general help in better understanding the influence of anatomy on function such as measured by EEG or fMRI.

New Clinically Feasible 3T MRI Protocol to Discriminate Internal Brain Stem Anatomy

Two new 3T MR imaging contrast methods, track density imaging and echo modulation curve T2 mapping, were combined with simultaneous multisection acquisition to reveal exquisite anatomic detail at 7 canonical levels of the brain stem. Compared with conventional MR imaging contrasts, many individual brain stem tracts and nuclear groups were directly visualized for the first time at 3T. This new approach is clinically practical and feasible (total scan time = 20 minutes), allowing better brain stem anatomic localization and characterization.

Whole-brain analytic measures of network communication reveal increased structure-function correlation in right temporal lobe epilepsy

The in vivo structure-function relationship is key to understanding brain network reorganization due to pathologies. This relationship is likely to be particularly complex in brain network diseases such as temporal lobe epilepsy, in which disturbed large-scale systems are involved in both transient electrical events and long-lasting functional and structural impairments. Herein, we estimated this relationship by analyzing the correlation between structural connectivity and functional connectivity in terms of analytical network communication parameters. As such, we targeted the gradual topological structure-function reorganization caused by the pathology not only at the whole brain scale but also both in core and peripheral regions of the brain. We acquired diffusion (dMRI) and resting-state fMRI (rsfMRI) data in seven right-lateralized TLE (rTLE) patients and fourteen healthy controls and analyzed the structure-function relationship by using analytical network communication metrics derived from the structural connectome. In rTLE patients, we found a widespread hypercorrelated functional network. Network communication analysis revealed greater unspecific branching of the shortest path (search information) in the structural connectome and a higher global correlation between the structural and functional connectivity for the patient group. We also found evidence for a preserved structural rich-club in the patient group. In sum, global augmentation of structure-function correlation might be linked to a smaller functional repertoire in rTLE patients, while sparing the central core of the brain which may represent a pathway that facilitates the spread of seizures.

Structural white-matter connections mediating distinct behavioral components of spatial neglect in right brain-damaged patients

Spatial neglect is a neuropsychological syndrome in which patients fail to perceive and orient to stimuli located in the space contralateral to the lesioned hemisphere. It is characterized by a wide heterogeneity in clinical symptoms which can be grouped into distinct behavioral components correlating with different lesion sites. Moreover, damage to white-matter (WM) fiber tracts has been suggested to disconnect brain networks that mediate different functions associated with spatial cognition and attention. However, it remains unclear what WM pathways are associated with functionally dissociable neglect components. In this study we examined nine patients with a focal right hemisphere stroke using a series of neuropsychological tests and diffusion tensor imaging (DTI) in order to disentangle the role of specific WM pathways in neglect symptoms. First, following previous work, the behavioral test scores of patients were factorized into three independent components reflecting perceptual, exploratory, and object-centered deficits in spatial awareness. We then examined the structural neural substrates of these components by correlating indices of WM integrity (fractional anisotropy) with the severity of deficits along each profile. Several locations in the right parietal and frontal WM correlated with neuropsychological scores. Fiber tracts projecting from these locations indicated that posterior parts of the superior longitudinal fasciculus (SLF), as well as nearby callosal fibers connecting ipsilateral and contralateral parietal areas, were associated with perceptual spatial deficits, whereas more anterior parts of SLF and inferior fronto-occipital fasciculus (IFOF) were predominantly associated with object-centered deficits. In addition, connections between frontal areas and superior colliculus were found to be associated with the exploratory deficits. Our results provide novel support to the view that neglect may result from disconnection lesions in distributed brain networks, but also extend these notions by highlighting the role of dissociable circuits in different functional components of the neglect syndrome. However these preliminary findings require replication with larger samples of patients.