Diffusion MRI fiber tractography of the brain

A tractometry principal component analysis of white matter tract network structure and relationships with cognitive function in relapsing-remitting multiple sclerosis

Understanding the brain changes underlying cognitive dysfunction is a key priority in multiple sclerosis (MS) to improve monitoring and treatment of this debilitating symptom. Functional connectivity network changes are associated with cognitive dysfunction, but it is less well understood how changes in normal appearing white matter relate to cognitive symptoms. If white matter tracts have network structure it would be expected that tracts within a network share susceptibility to MS pathology. In the present study, we used a tractometry approach to explore patterns of variance in white matter metrics across white matter (WM) tracts, and assessed how such patterns relate to neuropsychological test performance across cognitive domains. A sample of 102 relapsing-remitting MS patients and 27 healthy controls underwent MRI and neuropsychological testing. Tractography was performed on diffusion MRI data to extract 40 WM tracts and microstructural measures were extracted from each tract. Principal component analysis (PCA) was used to decompose metrics from all tracts to assess the presence of any co-variance structure among the tracts. Similarly, PCA was applied to cognitive test scores to identify the main cognitive domains. Finally, we assessed the ability of tract co-variance patterns to predict test performance across cognitive domains. We found that a single co-variance pattern which captured microstructure across all tracts explained the most variance (65% variance explained) and that there was little evidence for separate, smaller network patterns of pathology. Variance in this pattern was explained by effects related to lesions, but one main co-variance pattern persisted after this effect was regressed out. This main WM tract co-variance pattern contributed to explaining a modest degree of variance in one of our four cognitive domains in MS. These findings highlight the need to investigate the relationship between the normal appearing white matter and cognitive impairment further and on a more granular level, to improve the understanding of the network structure of the brain in MS.

The Case for Optimized Edge-Centric Tractography at Scale

The anatomic validity of structural connectomes remains a significant uncertainty in neuroimaging. Edge-centric tractography reconstructs streamlines in bundles between each pair of cortical or subcortical regions. Although edge bundles provides a stronger anatomic embedding than traditional connectomes, calculating them for each region-pair requires exponentially greater computation. We observe that major speedup can be achieved by reducing the number of streamlines used by probabilistic tractography algorithms. To ensure this does not degrade connectome quality, we calculate the identifiability of edge-centric connectomes between test and re-test sessions as a proxy for information content. We find that running PROBTRACKX2 with as few as 1 streamline per voxel per region-pair has no significant impact on identifiability. Variation in identifiability caused by streamline count is overshadowed by variation due to subject demographics. This finding even holds true in an entirely different tractography algorithm using MRTrix. Incidentally, we observe that Jaccard similarity is more effective than Pearson correlation in calculating identifiability for our subject population.

Diffusion tractography of superior cerebellar peduncle and dentatorubrothalamic tracts in two autopsy confirmed progressive supranuclear palsy variants: Richardson syndrome and the speech-language variant

•

Different changes in DTI metrics in SCP and DRTT can be seen across PSP subtypes.

•

DRTT tractography reconstructions demonstrated specific changes in PSP-RS.

•

DTI and clinical PSP scores are specifically linked across each PSP variant.

Background

Progressive supranuclear palsy (PSP) is a 4-repeat tauopathy with neurodegeneration typically observed in the superior cerebellar peduncle (SCP) and dentatorubrothalamic tracts (DRTT). However, it is unclear how these tracts are differentially affected in different clinical variants of PSP.

Objectives

To determine whether diffusion tractography of the SCP and DRTT can differentiate autopsy-confirmed PSP with Richardson’s syndrome (PSP-RS) and PSP with predominant speech/language disorder (PSP-SL).

Methods

We studied 22 autopsy-confirmed PSP patients that included 12 with PSP-RS and 10 with PSP-SL. We compared these two groups to 11 patients with autopsy-confirmed Alzheimer’s disease with SL problems, i.e., logopenic progressive aphasia (AD-LPA) (disease controls) and 10 healthy controls. Whole brain tractography was performed to identify the SCP and DRTT, as well as the frontal aslant tract and superior longitudinal fasciculus. We assessed fractional anisotropy and mean diffusivity for each tract. Hierarchical linear modeling was used for statistical comparisons, and correlations were assessed with clinical disease severity, ocular motor impairment, and parkinsonism. DRTT connectomics matrix analysis was also performed across groups.

Results

The SCP showed decreased fractional anisotropy for PSP-RS and PSP-SL and increased mean diffusivity in PSP-RS, compared to controls and AD-LPA. Right DRTT fibers showed lower fractional anisotropy in PSP-RS and PSP-SL compared to controls and AD-LPA, with PSP-RS also showing lower values compared to PSP-SL. Reductions in connectivity were observed in infratentorial DRTT regions in PSP-RS vs cortical regions in PSP-SL. PSP-SL showed greater abnormalities in the frontal aslant tract and superior longitudinal fasciculus compared to controls, PSP-RS, and AD-LPA. Significant correlations were observed between ocular motor impairment and SCP in PSP-RS (p = 0.042), and DRTT in PSP-SL (p = 0.022). In PSP-SL, the PSP Rating Scale correlated with the SCP (p = 0.045) and DRTT (p = 0.008), and the Unified Parkinson’s Disease Rating Scale correlated with the DRTT (p = 0.014).

Conclusions

Degeneration of the SCP and DRTT are diagnostic features of both PSP-RS and PSP-SL and associations with clinical metrics validate the role of these tracts in PSP-related clinical features, particularly in PSP-SL.

Assessment of structural disconnections in gliomas: comparison of indirect and direct approaches

Gliomas are amongst the most common primary brain tumours in adults and are often associated with poor prognosis. Understanding the extent of white matter (WM) which is affected outside the tumoral lesion may be of paramount importance to explain cognitive deficits and the clinical progression of the disease. To this end, we explored both direct (i.e., tractography based) and indirect (i.e., atlas-based) approaches to quantifying WM structural disconnections in a cohort of 44 high- and low-grade glioma patients. While these methodologies have recently gained popularity in the context of stroke and other pathologies, to our knowledge, this is the first time they are applied in patients with brain tumours. More specifically, in this work, we present a quantitative comparison of the disconnection maps provided by the two methodologies by applying well-known metrics of spatial similarity, extension, and correlation. Given the important role the oedematous tissue plays in the physiopathology of tumours, we performed these analyses both by including and excluding it in the definition of the tumoral lesion. This was done to investigate possible differences determined by this choice. We found that direct and indirect approaches offer two distinct pictures of structural disconnections in patients affected by brain gliomas, presenting key differences in several regions of the brain. Following the outcomes of our analysis, we eventually discuss the strengths and pitfalls of these two approaches when applied in this critical field.

Imaging versus electrographic connectivity in human mood-related fronto-temporal networks

BACKGROUND

The efficacy of psychiatric DBS is thought to be driven by the connectivity of stimulation targets with mood-relevant fronto-temporal networks, which is typically evaluated using diffusion-weighted tractography.

OBJECTIVE

Leverage intracranial electrophysiology recordings to better predict the circuit-wide effects of neuromodulation to white matter targets. We hypothesize strong convergence between tractography-predicted structural connectivity and stimulation-induced electrophysiological responses.

METHODS

Evoked potentials were elicited by single-pulse stimulation to two common DBS targets for treatment-resistant depression - the subcallosal cingulate (SCC) and ventral capsule/ventral striatum (VCVS) - in two patients undergoing DBS with stereo-electroencephalographic (sEEG) monitoring. Evoked potentials were compared with predicted structural connectivity between DBS leads and sEEG contacts using probabilistic, patient-specific diffusion-weighted tractography.

RESULTS

Evoked potentials and tractography showed strong convergence in both patients in orbitofrontal, ventromedial prefrontal, and lateral prefrontal cortices for both SCC and VCVS stimulation targets. Low convergence was found in anterior cingulate (ACC), where tractography predicted structural connectivity from SCC targets but produced no evoked potentials during SCC stimulation. Further, tractography predicted no connectivity to ACC from VCVS targets, but VCVS stimulation produced robust evoked potentials.

CONCLUSION

The two connectivity methods showed significant convergence, but important differences emerged with respect to the ability of tractography to predict electrophysiological connectivity between SCC and VCVS to regions of the mood-related network. This multimodal approach raises intriguing implications for the use of tractography in surgical targeting and provides new data to enhance our understanding of the network-wide effects of neuromodulation.

GPU-accelerated connectome discovery at scale

Diffusion magnetic resonance imaging and tractography enable the estimation of anatomical connectivity in the human brain, in vivo. Yet, without ground-truth validation, different tractography algorithms can yield widely varying connectivity estimates. Although streamline pruning techniques mitigate this challenge, slow compute times preclude their use in big-data applications. We present 'Regularized, Accelerated, Linear Fascicle Evaluation' (ReAl-LiFE), a GPU-based implementation of a state-of-the-art streamline pruning algorithm (LiFE), which achieves >100× speedups over previous CPU-based implementations. Leveraging these speedups, we overcome key limitations with LiFE's algorithm to generate sparser and more accurate connectomes. We showcase ReAl-LiFE's ability to estimate connections with superlative test-retest reliability, while outperforming competing approaches. Moreover, we predicted inter-individual variations in multiple cognitive scores with ReAl-LiFE connectome features. We propose ReAl-LiFE as a timely tool, surpassing the state of the art, for accurate discovery of individualized brain connectomes at scale. Finally, our GPU-accelerated implementation of a popular non-negative least-squares optimization algorithm is widely applicable to many real-world problems.

Establishing the functional relevancy of white matter connections in the visual system and beyond

For over a century, researchers have examined the functional relevancy of white matter bundles. Consequently, many large-scale bundles spanning several centimeters have been associated in their entirety with specific brain functions, such as language or attention. However, these coarse structural-functional relationships are at odds with modern understanding of the fine-grained functional organization of human cortex, such as the mosaic of category-selective regions in ventral temporal cortex. Here, we review a multimodal approach that combines fMRI to define functional regions of interest within individual's brains with dMRI tractography to identify the white matter bundles of the same individual. Combining these data allows to determine which subsets of streamlines within a white matter bundle connect to specific functional regions in each individual. That is, this approach identifies the functionally defined white matter sub-bundles of the brain. We argue that this approach not only enhances the accuracy of interpreting the functional relevancy of white matter bundles, but also enables segmentation of these large-scale bundles into meaningful functional units, which can then be linked to behavior with enhanced precision. Importantly, this approach has the potential for making new discoveries of the fine-grained functional relevancy of white matter connections in the visual system and the brain more broadly, akin to the flurry of research that has identified functional regions in cortex.

Methodological evaluation of individual cognitive prediction based on the brain white matter structural connectome

An emerging trend is to use regression‐based machine learning approaches to predict cognitive functions at the individual level from neuroimaging data. However, individual prediction models are inherently influenced by the vast options for network construction and model selection in machine learning pipelines. In particular, the brain white matter (WM) structural connectome lacks a systematic evaluation of the effects of different options in the pipeline on predictive performance. Here, we focused on the methodological evaluation of brain structural connectome‐based predictions. For network construction, we considered two parcellation schemes for defining nodes and seven strategies for defining edges. For the regression algorithms, we used eight regression models. Four cognitive domains and brain age were targeted as predictive tasks based on two independent datasets (Beijing Aging Brain Rejuvenation Initiative [BABRI]: 633 healthy older adults; Human Connectome Projects in Aging [HCP‐A]: 560 healthy older adults). Based on the results, the WM structural connectome provided a satisfying predictive ability for individual age and cognitive functions, especially for executive function and attention. Second, different parcellation schemes induce a significant difference in predictive performance. Third, prediction results from different data sets showed that dMRI with distinct acquisition parameters may plausibly result in a preference for proper fiber reconstruction algorithms and different weighting options. Finally, deep learning and Elastic‐Net models are more accurate and robust in connectome‐based predictions. Together, significant effects of different options in WM network construction and regression algorithms on the predictive performances are identified in this study, which may provide important references and guidelines to select suitable options for future studies in this field.

Fibre tract segmentation for intraoperative diffusion MRI in neurosurgical patients using tract-specific orientation atlas and tumour deformation modelling

Purpose:

Intraoperative diffusion MRI could provide a means of visualising brain fibre tracts near a neurosurgical target after preoperative images have been invalidated by brain shift. We propose an atlas-based intraoperative tract segmentation method, as the standard preoperative method, streamline tractography, is unsuitable for intraoperative implementation.

Methods:

A tract-specific voxel-wise fibre orientation atlas is constructed from healthy training data. After registration with a target image, a radial tumour deformation model is applied to the orientation atlas to account for displacement caused by lesions. The final tract map is obtained from the inner product of the atlas and target image fibre orientation data derived from intraoperative diffusion MRI.

Results:

The simple tumour model takes only seconds to effectively deform the atlas into alignment with the target image. With minimal processing time and operator effort, maps of surgically relevant tracts can be achieved that are visually and qualitatively comparable with results obtained from streamline tractography.

Conclusion:

Preliminary results demonstrate feasibility of intraoperative streamline-free tract segmentation in challenging neurosurgical cases. Demonstrated results in a small number of representative sample subjects are realistic despite the simplicity of the tumour deformation model employed. Following this proof of concept, future studies will focus on achieving robustness in a wide range of tumour types and clinical scenarios, as well as quantitative validation of segmentations.

Quantitative mapping of the brain's structural connectivity using diffusion MRI tractography: A review

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

Learning white matter subject-specific segmentation from structural MRI

PURPOSE

Mapping brain white matter (WM) is essential for building an understanding of brain anatomy and function. Tractography-based methods derived from diffusion-weighted MRI (dMRI) are the principal tools for investigating WM. These procedures rely on time-consuming dMRI acquisitions that may not always be available, especially for legacy or time-constrained studies. To address this problem, we aim to generate WM tracts from structural magnetic resonance imaging (MRI) image by deep learning.

METHODS

Following recently proposed innovations in structural anatomical segmentation, we evaluate the feasibility of training multiply spatial localized convolution neural networks to learn context from fixed spatial patches from structural MRI on standard template. We focus on six widely used dMRI tractography algorithms (TractSeg, RecoBundles, XTRACT, Tracula, automated fiber quantification (AFQ), and AFQclipped) and train 125 U-Net models to learn these techniques from 3870 T1-weighted images from the Baltimore Longitudinal Study of Aging, the Human Connectome Project S1200 release, and scans acquired at Vanderbilt University.

RESULTS

The proposed framework identifies fiber bundles with high agreement against tractography-based pathways with a median Dice coefficient from 0.62 to 0.87 on a test cohort, achieving improved subject-specific accuracy when compared to population atlas-based methods. We demonstrate the generalizability of the proposed framework on three externally available datasets.

CONCLUSIONS

We show that patch-wise convolutional neural network can achieve robust bundle segmentation from T1w. We envision the use of this framework for visualizing the expected course of WM pathways when dMRI is not available.

Automation of Cranial Nerve Tractography by Filtering Tractograms for Skull Base Surgery

Fiber tractography enables the in vivo reconstruction of white matter fibers in 3 dimensions using data collected by diffusion tensor imaging, thereby helping to understand functional neuroanatomy. In a pre-operative context, it provides essential information on the trajectory of fiber bundles of medical interest, such as cranial nerves. However, the optimization of tractography parameters is a time-consuming process and requires expert neuroanatomical knowledge, making the use of tractography difficult in clinical routine. Tractogram filtering is a method used to isolate the most relevant fibers. In this work, we propose to use filtering as a post-processing of tractography to avoid the manual optimization of tracking parameters and therefore making a step forward automation of tractography. To question the feasibility of automated tractography of cranial nerves, we perform an analysis of main cranial nerves on a series of patients with skull base tumors. A quantitative evaluation of the filtering performance of two state-of-the-art and a new entropy-based methods is carried out on the basis of reference tractograms produced by experts. Our approach proves to be more stable in the selection of the optimal filtering threshold and turns out to be interesting in terms of computational time complexity.

The orbitofrontal cortex of the sheep. Topography, organization, neurochemistry, digital tensor imaging and comparison with the chimpanzee and human

Areas dedicated to higher brain functions such as the orbitofrontal cortex (OFC) are thought to be unique to hominidae. The OFC is involved in social behavior, reward and punishment encoding and emotional control. Here, we focused on the putative corresponding area in the sheep to assess its homology to the OFC in humans. We used classical histology in five sheep (Ovis aries) and four chimpanzees (Pan troglodytes) as a six-layered-cortex primate, and Diffusion Tensor Imaging (DTI) in three sheep and five human brains. Nissl’s staining exhibited a certain alteration in cortical lamination since no layer IV was found in the sheep. A reduction of the total cortical thickness was also evident together with a reduction of the prevalence of layer one and an increased layer two on the total thickness. Tractography of the sheep OFC, on the other hand, revealed similarities both with human tracts and those described in the literature, as well as a higher number of cortico-cortical fibers connecting the OFC with the visual areas in the right hemisphere. Our results evidenced the presence of the basic components necessary for complex abstract thought in the sheep and a pronounced laterality, often associated with greater efficiency of a certain function, suggested an evolutionary adaptation of this prey species.

Structural Connectivity of Human Inferior Colliculus Subdivisions Using in vivo and post mortem Diffusion MRI Tractography

Inferior colliculus (IC) is an obligatory station along the ascending auditory pathway that also has a high degree of top-down convergence via efferent pathways, making it a major computational hub. Animal models have attributed critical roles for the IC in in mediating auditory plasticity, egocentric selection, and noise exclusion. IC contains multiple functionally distinct subdivisions. These include a central nucleus that predominantly receives ascending inputs and external and dorsal nuclei that receive more heterogeneous inputs, including descending and multisensory connections. Subdivisions of human IC have been challenging to identify and quantify using standard brain imaging techniques such as MRI, and the connectivity of each of these subnuclei has not been identified in the human brain. In this study, we estimated the connectivity of human IC subdivisions with diffusion MRI (dMRI) tractography, using both anatomical-based seed analysis as well as unsupervised k-means clustering. We demonstrate sensitivity of tractography to overall IC connections in both high resolution post mortem and in vivo datasets. k-Means clustering of the IC streamlines in both the post mortem and in vivo datasets generally segregated streamlines based on their terminus beyond IC, such as brainstem, thalamus, or contralateral IC. Using fine-grained anatomical segmentations of the major IC subdivisions, the post mortem dataset exhibited unique connectivity patterns from each subdivision, including commissural connections through dorsal IC and lateral lemniscal connections to central and external IC. The subdivisions were less distinct in the context of in vivo connectivity, although lateral lemniscal connections were again highest to central and external IC. Overall, the unsupervised and anatomically driven methods provide converging evidence for distinct connectivity profiles for each of the IC subdivisions in both post mortem and in vivo datasets, suggesting that dMRI tractography with high quality data is sensitive to neural pathways involved in auditory processing as well as top-down control of incoming auditory information.

Normal development of the brain: a survey of joint structural-functional brain studies

Joint structural-functional (S-F) developmental studies present a novel approach to address the complex neuroscience questions on how the human brain works and how it matures. Joint S-F biomarkers have the inherent potential to model effectively the brain's maturation, fill the information gap in temporal brain atlases, and demonstrate how the brain's performance matures during the lifespan. This review presents the current state of knowledge on heterochronous and heterogeneous development of S-F links during the maturation period. The S-F relationship has been investigated in early-matured unimodal and prolonged-matured transmodal regions of the brain using a variety of structural and functional biomarkers and data acquisition modalities. Joint S-F unimodal studies have employed auditory and visual stimuli, while the main focus of joint S-F transmodal studies has been resting-state and cognitive experiments. However, nonsignificant associations between some structural and functional biomarkers and their maturation show that designing and developing effective S-F biomarkers is still a challenge in the field. Maturational characteristics of brain asymmetries have been poorly investigated by the joint S-F studies, and the results were partially inconsistent with previous nonjoint ones. The inherent complexity of the brain performance can be modeled using multifactorial and nonlinear techniques as promising methods to simulate the impact of age on S-F relations considering their analysis challenges.

Diffusion MRI Connections in the Octopus Brain

Using high angle resolution diffusion magnetic resonance imaging (HARDI) with fiber tractography analysis we map out a meso-scale connectome of the Octopus bimaculoides brain. The brain of this cephalopod has a qualitatively different organization than that of vertebrates, yet it exhibits complex behavior, an elaborate sensory system and high cognitive abilities. Over the last 60 years wide ranging and detailed studies of octopus brain anatomy have been undertaken, including classical histological sectioning/staining, electron microscopy and neuronal tract tracing with injected dyes. These studies have elucidated many neuronal connections within and among anatomical structures. Diffusion MRI based tractography utilizes a qualitatively different method of tracing connections within the brain and offers facile three-dimensional images of anatomy and connections that can be quantitatively analyzed. Twenty-five separate lobes of the brain were segmented in the 3D MR images of each of three samples, including all five sub-structures in the vertical lobe. These parcellations were used to assay fiber tracings between lobes. The connectivity matrix constructed from diffusion MRI data was largely in agreement with that assembled from earlier studies. The one major difference was that connections between the vertical lobe and more basal supra-esophageal structures present in the literature were not found by MRI. In all, 92 connections between the 25 different lobes were noted by diffusion MRI: 53 between supra-esophageal lobes and 26 between the optic lobes and other structures. These represent the beginnings of a mesoscale connectome of the octopus brain.

Dynamic primitives of brain network interaction Special Issue "Advances in Mapping the Connectome"

What dynamic processes underly functional brain networks? Functional connectivity (FC) and functional connectivity dynamics (FCD) are used to represent the patterns and dynamics of functional brain networks. FC(D) is related to the synchrony of brain activity: when brain areas oscillate in a coordinated manner this yields a high correlation between their signal time series. To explain the processes underlying FC(D) we review how synchronized oscillations emerge from coupled neural populations in brain network models (BNMs). From detailed spiking networks to more abstract population models, there is strong support for the idea that the brain operates near critical instabilities that give rise to multistable or metastable dynamics that in turn lead to the intermittently synchronized slow oscillations underlying FC(D). We explore further consequences from these fundamental mechanisms and how they fit with reality. We conclude by highlighting the need for integrative brain models that connect separate mechanisms across levels of description and spatiotemporal scales and link them with cognitive function.

Volumetric Segmentation of White Matter Tracts with Label Embedding

Convolutional neural networks have achieved state-of-the-art performance for white matter (WM) tract segmentation based on diffusion magnetic resonance imaging (dMRI). However, the segmentation can still be difficult for challenging WM tracts with thin bodies or complicated shapes; the segmentation is even more problematic in challenging scenarios with reduced data quality or domain shift between training and test data, which can be easily encountered in clinical settings. In this work, we seek to improve the segmentation of WM tracts, especially for challenging WM tracts in challenging scenarios. In particular, our method is based on volumetric WM tract segmentation, where voxels are directly labeled without performing tractography. To improve the segmentation, we exploit the characteristics of WM tracts that different tracts can cross or overlap and revise the network design accordingly. Specifically, because multiple tracts can co-exist in a voxel, we hypothesize that the different tract labels can be correlated. The tract labels at a single voxel are concatenated as a label vector, the length of which is the number of tract labels. Due to the tract correlation, this label vector can be projected into a lower-dimensional space-referred to as the embedded space-for each voxel, which allows the segmentation network to solve a simpler problem. By predicting the coordinate in the embedded space for the tracts at each voxel and subsequently mapping the coordinate to the label vector with a reconstruction module, the segmentation result can be achieved. To facilitate the learning of the embedded space, an auxiliary label reconstruction loss is integrated with the segmentation accuracy loss during network training, and network training and inference are end-to-end. Our method was validated on two dMRI datasets under various settings. The results show that the proposed method improves the accuracy of WM tract segmentation, and the improvement is more prominent for challenging tracts in challenging scenarios.

Potential Pitfalls of Using Fractional Anisotropy, Axial Diffusivity, and Radial Diffusivity as Biomarkers of Cerebral White Matter Microstructure

Fractional anisotropy (FA), axial diffusivity (AD), and radial diffusivity (RD) are commonly used as MRI biomarkers of white matter microstructure in diffusion MRI studies of neurodevelopment, brain aging, and neurologic injury/disease. Some of the more frequent practices include performing voxel-wise or region-based analyses of these measures to cross-sectionally compare individuals or groups, longitudinally assess individuals or groups, and/or correlate with demographic, behavioral or clinical variables. However, it is now widely recognized that the majority of cerebral white matter voxels contain multiple fiber populations with different trajectories, which renders these metrics highly sensitive to the relative volume fractions of the various fiber populations, the microstructural integrity of each constituent fiber population, and the interaction between these factors. Many diffusion imaging experts are aware of these limitations and now generally avoid using FA, AD or RD (at least in isolation) to draw strong reverse inferences about white matter microstructure, but based on the continued application and interpretation of these metrics in the broader biomedical/neuroscience literature, it appears that this has perhaps not yet become common knowledge among diffusion imaging end-users. Therefore, this paper will briefly discuss the complex biophysical underpinnings of these measures in the context of crossing fibers, provide some intuitive “thought experiments” to highlight how conventional interpretations can lead to incorrect conclusions, and suggest that future studies refrain from using (over-interpreting) FA, AD, and RD values as standalone biomarkers of cerebral white matter microstructure.

High fidelity fiber orientation density functions from fiber ball imaging

The fiber orientation density function (fODF) in white matter is a primary physical quantity that can be estimated with diffusion MRI. It has often been employed for fiber tracking and microstructural modeling. Requirements for the construction of high fidelity fODFs, in the sense of having good angular resolution, adequate data to avoid sampling errors, and minimal noise artifacts, are described for fODFs calculated with fiber ball imaging. A criterion is formulated for the number of diffusion encoding directions needed to achieve a given angular resolution. The advantages of using large b-values (≥6000 s/mm² ) are also discussed. For the direct comparison of different fODFs, a method is developed for defining a local frame of reference tied to each voxel's individual axonal structure. The Matusita anisotropy axonal is proposed as a scalar fODF measure for quantifying angular variability. Experimental results, obtained at 3 T from human volunteers, are used as illustrations.

Fluid mechanics approach to analyzing collagen fiber organization

SIGNIFICANCE

The spatial organization of collagen fibers has been used as a biomarker for assessing injury and disease progression. However, quantifying this organization for complex structures is challenging.

AIM

To quantify and classify complex collagen fiber organizations.

APPROACH

Using quantitative second-harmonic generation (SHG) microscopy, we show that collagen-fiber orientation can be viewed as pseudovector fields. Subsequently, we analyze them using fluid mechanic metrics, such as energy U, enstrophy E, and tortuosity τ.

RESULTS

We show that metrics used in fluid mechanics for analyzing fluid flow can be adapted to analyze complex collagen fiber organization. As examples, we consider SHG images of collagenous tissue for straight, wavy, and circular fiber structures.

CONCLUSIONS

The results of this study show the utility of the chosen metrics to distinguish diverse and complex collagen organizations. We find that the distribution of values for E and U increases with collagen fiber disorganization, where they divide between low and high values corresponding to uniformly aligned fibers and disorganized collagen fibers, respectively. We also confirm that the values of τ cluster around 1 when the fibers are straight, and the range increases up to 1.5 when wavier fibers are present.

Utilizing the TractSeg Tool for Automatic Corticospinal Tract Segmentation in Patients With Brain Pathology

Purpose: White-matter tract segmentation in patients with brain pathology can guide surgical planning and can be used for tissue integrity assessment. Recently, TractSeg was proposed for automatic tract segmentation in healthy subjects. The aim of this study was to assess the use of TractSeg for corticospinal-tract (CST) segmentation in a large cohort of patients with brain pathology and to evaluate its consistency in repeated measurements. Methods: A total of 649 diffusion-tensor-imaging scans were included, of them: 625 patients and 24 scans from 12 healthy controls (scanned twice for consistency assessment). Manual CST labeling was performed in all cases, and by 2 raters for the healthy subjects. Segmentation results were evaluated based on the Dice score. In order to evaluate consistency in repeated measurements, volume, Fractional Anisotropy (FA), and Mean Diffusivity (MD) values were extracted and correlated for the manual versus automatic methods. Results: For the automatic CST segmentation Dice scores of 0.63 and 0.64 for the training and testing datasets were obtained. Higher consistency between measurements was detected for the automatic segmentation, with between measurements correlations of volume = 0.92/0.65, MD = 0.94/0.75 for the automatic versus manual segmentation. Conclusions: The TractSeg method enables automatic CST segmentation in patients with brain pathology. Superior measurements consistency was detected for the automatic in comparison to manual fiber segmentation, which indicates an advantage when using this method for clinical and longitudinal studies.

Investigating the Anatomy and Microstructure of the Dentato-rubro-thalamic and Subthalamo-ponto-cerebellar Tracts in Parkinson's Disease

Cerebellar-thalamic connections play a central role in deep brain stimulation-based treatment of tremor syndromes. Here, we used diffusion Magnetic Resonance Imaging (MRI) tractography to delineate the main cerebellar peduncles as well as two main white matter tracts that connect the cerebellum with the thalamus, the dentato-rubro-thalamic tract (DRTT) and the subthalamo-ponto-cerebellar tract (SPCT). We first developed a reconstruction protocol in young healthy adults with high-resolution diffusion imaging data and then demonstrate feasibility of transferring this protocol to clinical studies using standard diffusion MRI data from a cohort of patients with Parkinson's disease (PD) and their matched healthy controls. The tracts obtained closely corresponded to the previously described anatomical pathways and features of the DRTT and the SPCT. Second, we investigated the microstructure of these tracts with fractional anisotropy (FA), radial diffusivity (RD), and hindrance modulated orientational anisotropy (HMOA) in patients with PD and healthy controls. By reducing dimensionality of both the microstructural metrics and the investigated cerebellar and cerebellar-thalamic tracts using principal component analyses, we found global differences between patients with PD and controls, suggestive of higher fractional anisotropy, lower radial diffusivity, and higher hindrance modulated orientational anisotropy in patients. However, separate analyses for each of the tracts did not yield any significant differences. Our findings contribute to the characterization of the distinct anatomical connections between the cerebellum and the diencephalon. Microstructural differences between patients and controls in the cerebellar pathways suggest involvement of these structures in PD, complementing previous functional and diffusion imaging studies.

Toward personalized medicine in connectomic deep brain stimulation

At the group-level, deep brain stimulation leads to significant therapeutic benefit in a multitude of neurological and neuropsychiatric disorders. At the single-patient level, however, symptoms may sometimes persist despite "optimal" electrode placement at established treatment coordinates. This may be partly explained by limitations of disease-centric strategies that are unable to account for heterogeneous phenotypes and comorbidities observed in clinical practice. Instead, tailoring electrode placement and programming to individual patients' symptom profiles may increase the fraction of top-responding patients. Here, we propose a three-step, circuit-based framework with the aim of developing patient-specific treatment targets that address the unique symptom constellation prevalent in each patient. First, we describe how a symptom network target library could be established by mapping beneficial or undesirable DBS effects to distinct circuits based on (retrospective) group-level data. Second, we suggest ways of matching the resulting symptom networks to circuits defined in the individual patient (template matching). Third, we introduce network blending as a strategy to calculate optimal stimulation targets and parameters by selecting and weighting a set of symptom-specific networks based on the symptom profile and subjective priorities of the individual patient. We integrate the approach with published literature and conclude by discussing limitations and future challenges.

Deterministic Tractography Analysis of Rat Brain Using SIGMA Atlas in 9.4T MRI

Preclinical studies using rodents have been the choice for many neuroscience researchers due totheir close reflection of human biology. In particular, research involving rodents has utilized MRI to accurately identify brain regions and characteristics by acquiring high resolution cavity images with different contrasts non-invasively, and this has resulted in high reproducibility and throughput. In addition, tractographic analysis using diffusion tensor imaging to obtain information on the neural structure of white matter has emerged as a major methodology in the field of neuroscience due to its contribution in discovering significant correlations between altered neural connections and various neurological and psychiatric diseases. However, unlike image analysis studies with human subjects where a myriad of human image analysis programs and procedures have been thoroughly developed and validated, methods for analyzing rat image data using MRI in preclinical research settings have seen significantly less developed. Therefore, in this study, we present a deterministic tractographic analysis pipeline using the SIGMA atlas for a detailed structural segmentation and structural connectivity analysis of the rat brain’s structural connectivity. In addition, the structural connectivity analysis pipeline presented in this study was preliminarily tested on normal and stroke rat models for initial observation.

Prevalence of white matter pathways coming into a single white matter voxel orientation: The bottleneck issue in tractography

Characterizing and understanding the limitations of diffusion MRI fiber tractography is a prerequisite for methodological advances and innovations which will allow these techniques to accurately map the connections of the human brain. The so-called "crossing fiber problem" has received tremendous attention and has continuously triggered the community to develop novel approaches for disentangling distinctly oriented fiber populations. Perhaps an even greater challenge occurs when multiple white matter bundles converge within a single voxel, or throughout a single brain region, and share the same parallel orientation, before diverging and continuing towards their final cortical or sub-cortical terminations. These so-called "bottleneck" regions contribute to the ill-posed nature of the tractography process, and lead to both false positive and false negative estimated connections. Yet, as opposed to the extent of crossing fibers, a thorough characterization of bottleneck regions has not been performed. The aim of this study is to quantify the prevalence of bottleneck regions. To do this, we use diffusion tractography to segment known white matter bundles of the brain, and assign each bundle to voxels they pass through and to specific orientations within those voxels (i.e. fixels). We demonstrate that bottlenecks occur in greater than 50-70% of fixels in the white matter of the human brain. We find that all projection, association, and commissural fibers contribute to, and are affected by, this phenomenon, and show that even regions traditionally considered "single fiber voxels" often contain multiple fiber populations. Together, this study shows that a majority of white matter presents bottlenecks for tractography which may lead to incorrect or erroneous estimates of brain connectivity or quantitative tractography (i.e., tractometry), and underscores the need for a paradigm shift in the process of tractography and bundle segmentation for studying the fiber pathways of the human brain.

Nonparametric D-R1-R2 distribution MRI of the living human brain

Diffusion-relaxation correlation NMR can simultaneously characterize both the microstructure and the local chemical composition of complex samples that contain multiple populations of water. Recent developments on tensor-valued diffusion encoding and Monte Carlo inversion algorithms have made it possible to transfer diffusion-relaxation correlation NMR from small-bore scanners to clinical MRI systems. Initial studies on clinical MRI systems employed 5D D-R₁ and D-R₂ correlation to characterize healthy brain in vivo. However, these methods are subject to an inherent bias that originates from not including R₂ or R₁ in the analysis, respectively. This drawback can be remedied by extending the concept to 6D D-R₁-R₂ correlation. In this work, we present a sparse acquisition protocol that records all data necessary for in vivo 6D D-R₁-R₂ correlation MRI across 633 individual measurements within 25 min-a time frame comparable to previous lower-dimensional acquisition protocols. The data were processed with a Monte Carlo inversion algorithm to obtain nonparametric 6D D-R₁-R₂ distributions. We validated the reproducibility of the method in repeated measurements of healthy volunteers. For a post-therapy glioblastoma case featuring cysts, edema, and partially necrotic remains of tumor, we present representative single-voxel 6D distributions, parameter maps, and artificial contrasts over a wide range of diffusion-, R₁-, and R₂-weightings based on the rich information contained in the D-R₁-R₂ distributions.

Multi-tissue spherical deconvolution of tensor-valued diffusion MRI

Multi-tissue constrained spherical deconvolution (MT-CSD) leverages the characteristic b-value dependency of each tissue type to estimate both the apparent tissue densities and the white matter fiber orientation distribution function from diffusion MRI data. In this work, we generalize MT-CSD to tensor-valued diffusion encoding with arbitrary b-tensor shapes. This enables the use of data encoded with mixed b-tensors, rather than being limited to the subset of linear (conventional) b-tensors. Using the complete set of data, including all b-tensor shapes, provides a categorical improvement in the estimation of apparent tissue densities, fiber ODF, and resulting tractography. Furthermore, we demonstrate that including multiple b-tensor shapes in the analysis provides improved contrast between tissue types, in particular between gray matter and white matter. We also show that our approach provides high-quality apparent tissue density maps and high-quality fiber tracking from data, even with sparse sampling across b-tensors that yield whole-brain coverage at 2 mm isotropic resolution in approximately 5:15 min.

Structural Brain Network Reproducibility: Influence of Different Diffusion Acquisition and Tractography Reconstruction Schemes on Graph Metrics

Background: Graph metrics of structural brain networks demonstrate to be a powerful tool for investigating brain topology at a large scale. However, the variability of the results related to applying different magnetic resonance acquisition schemes and tractography reconstruction techniques is not fully characterized. Materials and Methods: The present work aims to evaluate the influence of different combinations of diffusion acquisition schemes (single and multishell), diffusion models (tensor and spherical deconvolution), and tractography reconstruction approaches (deterministic and probabilistic) on the reproducibility of graph metrics derived from structural connectome on test/retest (TRT) data released by the Human Connectome Project. From each implemented experimental setup, both global and local graph metrics were evaluated and their reproducibility was estimated by the intraclass correlation coefficient (ICC). Moreover, the percentage relative standard deviation (pRSD) from the ICC values of local graph metrics was calculated to quantify how much the reproducibility varied across nodes within each experimental setup. Results: The presented results show that different combinations of diffusion acquisition schemes, diffusion models, and tractography algorithms can strongly affect the reproducibility of global and local graph metrics. The combination of constrained spherical deconvolution (CSD) and deterministic tractography gave generally high reproducibility (ICCs >0.75) and lowest pRSD for the considered graph metrics, meanwhile probabilistic CSD with a high b-value returned the highest reproducibility. Notably, the introduction of streamline selection filters on CSD can substantially affect the reproducibility. Discussion: This work demonstrates that the TRT reproducibility of graph metrics is generally high but can vary substantially with different combinations of acquisition and reconstruction schemes. Impact statement This work demonstrates the influence of different diffusion acquisition schemes, diffusion models, and tractography reconstruction approaches on the reproducibility of graph metrics derived from structural connectome. The presented findings impact on the choice of both acquisition protocol and processing pipeline for topological analyses to produce reproducible measurements for brain network studies.

A critical review of connectome validation studies

Diffusion MRI tractography is the most widely used macroscale method for mapping connectomes in vivo. However, tractography is prone to various errors and biases, and thus tractography-derived connectomes require careful validation. Here, we critically review studies that have developed or utilized phantoms and tracer maps to validate tractography-derived connectomes, either quantitatively or qualitatively. We identify key factors impacting connectome reconstruction accuracy, including streamline seeding, propagation and filtering methods, and consider the strengths and limitations of state-of-the-art connectome phantoms and associated validation studies. These studies demonstrate the inherent limitations of current fiber orientation models and tractography algorithms and their impact on connectome reconstruction accuracy. Reconstructing connectomes with both high sensitivity and high specificity is challenging, given that some tractography methods can generate an abundance of spurious connections, while others can overlook genuine fiber bundles. We argue that streamline filtering can minimize spurious connections and potentially improve the biological plausibility of connectomes derived from tractography. We find that algorithmic choices such as the tractography seeding methodology, angular threshold, and streamline propagation method can substantially impact connectome reconstruction accuracy. Hence, careful application of tractography is necessary to reconstruct accurate connectomes. Improvements in diffusion MRI acquisition techniques will not necessarily overcome current tractography limitations without accompanying modeling and algorithmic advances.

Development of the visual white matter pathways mediates development of electrophysiological responses in visual cortex

The latency of neural responses in the visual cortex changes systematically across the lifespan. Here, we test the hypothesis that development of visual white matter pathways mediates maturational changes in the latency of visual signals. Thirty-eight children participated in a cross-sectional study including diffusion magnetic resonance imaging (MRI) and magnetoencephalography (MEG) sessions. During the MEG acquisition, participants performed a lexical decision and a fixation task on words presented at varying levels of contrast and noise. For all stimuli and tasks, early evoked fields were observed around 100 ms after stimulus onset (M100), with slower and lower amplitude responses for low as compared to high contrast stimuli. The optic radiations and optic tracts were identified in each individual's brain based on diffusion MRI tractography. The diffusion properties of the optic radiations predicted M100 responses, especially for high contrast stimuli. Higher optic radiation fractional anisotropy (FA) values were associated with faster and larger M100 responses. Over this developmental window, the M100 responses to high contrast stimuli became faster with age and the optic radiation FA mediated this effect. These findings suggest that the maturation of the optic radiations over childhood accounts for individual variations observed in the developmental trajectory of visual cortex responses.

The role of spatial embedding in mouse brain networks constructed from diffusion tractography and tracer injections

Diffusion MRI tractography is the only noninvasive method to measure the structural connectome in humans. However, recent validation studies have revealed limitations of modern tractography approaches, which lead to significant mistracking caused in part by local uncertainties in fiber orientations that accumulate to produce larger errors for longer streamlines. Characterizing the role of this length bias in tractography is complicated by the true underlying contribution of spatial embedding to brain topology. In this work, we compare graphs constructed with ex vivo tractography data in mice and neural tracer data from the Allen Mouse Brain Connectivity Atlas to random geometric surrogate graphs which preserve the low-order distance effects from each modality in order to quantify the role of geometry in various network properties. We find that geometry plays a substantially larger role in determining the topology of graphs produced by tractography than graphs produced by tracers. Tractography underestimates weights at long distances compared to neural tracers, which leads tractography to place network hubs close to the geometric center of the brain, as do corresponding tractography-derived random geometric surrogates, while tracer graphs place hubs further into peripheral areas of the cortex. We also explore the role of spatial embedding in modular structure, network efficiency and other topological measures in both modalities. Throughout, we compare the use of two different tractography streamline node assignment strategies and find that the overall differences between tractography approaches are small relative to the differences between tractography- and tracer-derived graphs. These analyses help quantify geometric biases inherent to tractography and promote the use of geometric benchmarking in future tractography validation efforts.

The structural connectome in ADHD

Attention-deficit/hyperactivity disorder (ADHD) has been conceptualized as a brain dysconnectivity disorder. In the past decade, noninvasive diffusion magnetic resonance imaging (dMRI) studies have demonstrated that individuals with ADHD have alterations in the white matter structural connectome, and that these alterations are associated with core symptoms and cognitive deficits in patients. This review aims to summarize recent dMRI-based structural connectome studies in ADHD from voxel-, tractography-, and network-based perspectives. Voxel- and tractography-based studies have demonstrated disrupted microstructural properties predominantly located in the frontostriatal tracts, the corpus callosum, the corticospinal tracts, and the cingulum bundle in patients with ADHD. Network-based studies have suggested abnormal global and local efficiency as well as nodal properties in the prefrontal and parietal regions in the ADHD structural connectomes. The altered structural connectomes in those with ADHD provide significant signatures for prediction of symptoms and diagnostic classification. These studies suggest that abnormalities in the structural connectome may be one of the neural underpinnings of ADHD psychopathology and show potential for establishing imaging biomarkers in clinical evaluation. However, given that there are inconsistent findings across studies due to sample heterogeneity and analysis method variations, these ADHD-related white matter alterations are still far from informing clinical practice. Future studies with larger and more homogeneous samples are needed to validate the consistency of current results; advanced dMRI techniques can help to generate much more precise estimation of white matter pathways and assure specific fiber configurations; and finally, dimensional analysis frameworks can deepen our understanding of the neurobiology underlying ADHD.

Anatomy and physiology of word-selective visual cortex: from visual features to lexical processing

Over the past 2 decades, researchers have tried to uncover how the human brain can extract linguistic information from a sequence of visual symbols. The description of how the brain's visual system processes words and enables reading has improved with the progressive refinement of experimental methodologies and neuroimaging techniques. This review provides a brief overview of this research journey. We start by describing classical models of object recognition in non-human primates, which represent the foundation for many of the early models of visual word recognition in humans. We then review functional neuroimaging studies investigating the word-selective regions in visual cortex. This research led to the differentiation of highly specialized areas, which are involved in the analysis of different aspects of written language. We then consider the corresponding anatomical measurements and provide a description of the main white matter pathways carrying neural signals crucial to word recognition. Finally, in an attempt to integrate structural, functional, and electrophysiological findings, we propose a view of visual word recognition, accounting for spatial and temporal facets of word-selective neural processes. This multi-modal perspective on the neural circuitry of literacy highlights the relevance of a posterior-anterior differentiation in ventral occipitotemporal cortex for visual processing of written language and lexical features. It also highlights unanswered questions that can guide us towards future research directions. Bridging measures of brain structure and function will help us reach a more precise understanding of the transformation from vision to language.

Diffusion kurtosis imaging of white matter in bipolar disorder

White matter pathology likely contributes to the pathogenesis of bipolar disorder (BD). Most studies of white matter in BD have used diffusion tensor imaging (DTI), but the advent of more advanced multi-shell diffusion MRI imaging offers the possibility to investigate other aspects of white matter microstructure. Diffusion kurtosis imaging (DKI) extends the DTI model and provides additional measures related to diffusion restriction. Here, we investigated white matter in BD by applying whole-brain voxel-based analysis (VBA) and a network-based connectivity approach using constrained spherical deconvolution tractography to assess differences in DKI and DTI metrics between BD (n = 25) and controls (n = 24). The VBA showed lower mean kurtosis in the corona radiata and posterior association fibers in BD. Regional differences in connectivity were indicated by lower mean kurtosis and kurtosis anisotropy in streamlines traversing the temporal and occipital lobes, and lower mean axial kurtosis in the right cerebellar, thalamo-subcortical pathways in BD. Significant differences were not seen in DTI metrics following FDR-correction. The DKI findings indicate altered connectivity across cortical, subcortical and cerebellar areas in BD. DKI is sensitive to different microstructural properties and is a useful complementary technique to DTI to more fully investigate white matter in BD.

Parcellation-Free prediction of task fMRI activations from dMRI tractography

The relationship between brain structure and function plays a crucial role in cognitive and clinical neuroscience. We present a supervised machine learning based approach that captures this relationship by predicting the spatial extent of activations that are observed with task based functional Magnetic Resonance Imaging (fMRI) from the local white matter connectivity, as reflected in diffusion MRI (dMRI) tractography. In particular, we explore three different feature representations of local connectivity patterns that do not require a pre-defined parcellation of cortical and subcortical structures. Instead, they employ cluster-based Bag of Features, Gaussian Mixture Models, and Fisher vectors. We demonstrate that our framework can be used to test the statistical significance of structure-function relationships, compare it to parcellation-based and group-average benchmarks, and propose an algorithm for visualizing our chosen feature representations that permits a neuroanatomical interpretation of our results.

Fiber tractography bundle segmentation depends on scanner effects, vendor effects, acquisition resolution, diffusion sampling scheme, diffusion sensitization, and bundle segmentation workflow

When investigating connectivity and microstructure of white matter pathways of the brain using diffusion tractography bundle segmentation, it is important to understand potential confounds and sources of variation in the process. While cross-scanner and cross-protocol effects on diffusion microstructure measures are well described (in particular fractional anisotropy and mean diffusivity), it is unknown how potential sources of variation effect bundle segmentation results, which features of the bundle are most affected, where variability occurs, nor how these sources of variation depend upon the method used to reconstruct and segment bundles. In this study, we investigate six potential sources of variation, or confounds, for bundle segmentation: variation (1) across scan repeats, (2) across scanners, (3) across vendors (4) across acquisition resolution, (5) across diffusion schemes, and (6) across diffusion sensitization. We employ four different bundle segmentation workflows on two benchmark multi-subject cross-scanner and cross-protocol databases, and investigate reproducibility and biases in volume overlap, shape geometry features of fiber pathways, and microstructure features within the pathways. We find that the effects of acquisition protocol, in particular acquisition resolution, result in the lowest reproducibility of tractography and largest variation of features, followed by vendor-effects, scanner-effects, and finally diffusion scheme and b-value effects which had similar reproducibility as scan-rescan variation. However, confounds varied both across pathways and across segmentation workflows, with some bundle segmentation workflows more (or less) robust to sources of variation. Despite variability, bundle dissection is consistently able to recover the same location of pathways in the deep white matter, with variation at the gray matter/ white matter interface. Next, we show that differences due to the choice of bundle segmentation workflows are larger than any other studied confound, with low-to-moderate overlap of the same intended pathway when segmented using different methods. Finally, quantifying microstructure features within a pathway, we show that tractography adds variability over-and-above that which exists due to noise, scanner effects, and acquisition effects. Overall, these confounds need to be considered when harmonizing diffusion datasets, interpreting or combining data across sites, and when attempting to understand the successes and limitations of different methodologies in the design and development of new tractography or bundle segmentation methods.

Striatal topographical organization: Bridging the gap between molecules, connectivity and behavior

The striatum represents the major hub of the basal ganglia, receiving projections from the entire cerebral cortex and it is assumed to play a key role in a wide array of complex behavioral tasks. Despite being extensively investigated during the last decades, the topographical organization of the striatum is not well understood yet. Ongoing efforts in neuroscience are focused on analyzing striatal anatomy at different spatial scales, to understand how structure relates to function and how derangements of this organization are involved in various neuropsychiatric diseases. While being subdivided at the macroscale level into dorsal and ventral divisions, at a mesoscale level the striatum represents an anatomical continuum sharing the same cellular makeup. At the same time, it is now increasingly ascertained that different striatal compartments show subtle histochemical differences, and their neurons exhibit peculiar patterns of gene expression, supporting functional diversity across the whole basal ganglia circuitry. Such diversity is further supported by afferent connections which are heterogenous both anatomically, as they originate from distributed cortical areas and subcortical structures, and biochemically, as they involve a variety of neurotransmitters. Specifically, the cortico-striatal projection system is topographically organized delineating a functional organization which is maintained throughout the basal ganglia, subserving motor, cognitive and affective behavioral functions. While such functional heterogeneity has been firstly conceptualized as a tripartite organization, with sharply defined limbic, associative and sensorimotor territories within the striatum, it has been proposed that such territories are more likely to fade into one another, delineating a gradient-like organization along medio-lateral and ventro-dorsal axes. However, the molecular and cellular underpinnings of such organization are less understood, and their relations to behavior remains an open question, especially in humans. In this review we aimed at summarizing the available knowledge on striatal organization, especially focusing on how it links structure to function and its alterations in neuropsychiatric diseases. We examined studies conducted on different species, covering a wide array of different methodologies: from tract-tracing and immunohistochemistry to neuroimaging and transcriptomic experiments, aimed at bridging the gap between macroscopic and molecular levels.

Development of white matter microstructure and executive functions during childhood and adolescence: a review of diffusion MRI studies

Diffusion magnetic resonance imaging (dMRI) provides indirect measures of white matter microstructure that can be used to make inferences about structural connectivity within the brain. Over the last decade, a growing literature of cross-sectional and longitudinal studies have documented relationships between dMRI indices and cognitive development. In this review, we provide a brief overview of dMRI methods and how they can be used to study white matter and connectivity and review the extant literature examining the links between dMRI indices and executive functions during development. We explore the links between white matter microstructure and specific executive functions: inhibition, working memory and cognitive shifting, as well as performance on complex executive function tasks. Concordance in findings across studies are highlighted, and potential explanations for discrepancies between results, together with challenges with using dMRI in child and adolescent populations, are discussed. Finally, we explore future directions that are necessary to better understand the links between child and adolescent development of structural connectivity of the brain and executive functions.

Diffusion-Weighted Imaging: Recent Advances and Applications

Quantitative diffusion imaging techniques enable the characterization of tissue microstructural properties of the human brain "in vivo", and are widely used in neuroscientific and clinical contexts. In this review, we present the basic physical principles behind diffusion imaging and provide an overview of the current diffusion techniques, including standard and advanced techniques as well as their main clinical applications. Standard diffusion tensor imaging (DTI) offers sensitivity to changes in microstructure due to diseases and enables the characterization of single fiber distributions within a voxel as well as diffusion anisotropy. Nonetheless, its inability to represent complex intravoxel fiber topologies and the limited biological specificity of its metrics motivated the development of several advanced diffusion MRI techniques. For example, high-angular resolution diffusion imaging (HARDI) techniques enabled the characterization of fiber crossing areas and other complex fiber topologies in a single voxel and supported the development of higher-order signal representations aiming to decompose the diffusion MRI signal into distinct microstructure compartments. Biophysical models, often known by their acronym (e.g., CHARMED, WMTI, NODDI, DBSI, DIAMOND) contributed to capture the diffusion properties from each of such tissue compartments, enabling the computation of voxel-wise maps of axonal density and/or morphology that hold promise as clinically viable biomarkers in several neurological and neuroscientific applications; for example, to quantify tissue alterations due to disease or healthy processes. Current challenges and limitations of state-of-the-art models are discussed, including validation efforts. Finally, novel diffusion encoding approaches (e.g., b-tensor or double diffusion encoding) may increase the biological specificity of diffusion metrics towards intra-voxel diffusion heterogeneity in clinical settings, holding promise in neurological applications.

Functionnectome as a framework to analyse the contribution of brain circuits to fMRI

In recent years, the field of functional neuroimaging has moved away from a pure localisationist approach of isolated functional brain regions to a more integrated view of these regions within functional networks. However, the methods used to investigate functional networks rely on local signals in grey matter and are limited in identifying anatomical circuitries supporting the interaction between brain regions. Mapping the brain circuits mediating the functional signal between brain regions would propel our understanding of the brain's functional signatures and dysfunctions. We developed a method to unravel the relationship between brain circuits and functions: The Functionnectome. The Functionnectome combines the functional signal from fMRI with white matter circuits' anatomy to unlock and chart the first maps of functional white matter. To showcase this method's versatility, we provide the first functional white matter maps revealing the joint contribution of connected areas to motor, working memory, and language functions. The Functionnectome comes with an open-source companion software and opens new avenues into studying functional networks by applying the method to already existing datasets and beyond task fMRI.

Active Cortex Tractography

Most existing diffusion tractography algorithms are affected by gyral bias, causing the termination of streamlines at gyral crowns instead of sulcal banks. In this paper, we propose a tractography technique, called active cortex tractography (ACT), to overcome gyral bias by enabling fiber streamlines to curve naturally into the cortex. We show that the cortex can play an active role in cortical tractography by providing anatomical information to overcome orientation ambiguities as the streamlines enter the superficial white matter in gyral blades and approach the cortex. This is achieved by devising a direction scouting mechanism that takes into account the white matter surface normal vectors. The scouting mechanism allows probing of directions further in space to prepare the streamlines to turn at appropriate angles. The surface normal vectors guide the streamlines to turn into the cortex, perpendicular to the white-gray matter interface. Evaluation using synthetic, macaque and human data with different streamline seeding schemes demonstrates that ACT improves cortical tractography.

Ventral intermediate nucleus structural connectivity-derived segmentation: anatomical reliability and variability

The Ventral intermediate nucleus (Vim) of thalamus is the most targeted structure for the treatment of drug-refractory tremors. Since methodological differences across existing studies are remarkable and no gold-standard pipeline is available, in this study, we tested different parcellation pipelines for tractography-derived putative Vim identification. Thalamic parcellation was performed on a high quality, multi-shell dataset and a downsampled, clinical-like dataset using two different diffusion signal modeling techniques and two different voxel classification criteria, thus implementing a total of four parcellation pipelines. The most reliable pipeline in terms of inter-subject variability has been picked and parcels putatively corresponding to motor thalamic nuclei have been selected by calculating similarity with a histology-based mask of Vim. Then, spatial relations with optimal stimulation points for the treatment of essential tremor have been quantified. Finally, effect of data quality and parcellation pipelines on a volumetric index of connectivity clusters has been assessed. We found that the pipeline characterized by higher-order signal modeling and threshold-based voxel classification criteria was the most reliable in terms of inter-subject variability regardless data quality. The maps putatively corresponding to Vim were those derived by precentral and dentate nucleus-thalamic connectivity. However, tractography-derived functional targets showed remarkable differences in shape and sizes when compared to a ground truth model based on histochemical staining on seriate sections of human brain. Thalamic voxels connected to contralateral dentate nucleus resulted to be the closest to literature-derived stimulation points for essential tremor but at the same time showing the most remarkable inter-subject variability. Finally, the volume of connectivity parcels resulted to be significantly influenced by data quality and parcellation pipelines. Hence, caution is warranted when performing thalamic connectivity-based segmentation for stereotactic targeting.

Comparison of multiple tractography methods for reconstruction of the retinogeniculate visual pathway using diffusion MRI

The retinogeniculate visual pathway (RGVP) conveys visual information from the retina to the lateral geniculate nucleus. The RGVP has four subdivisions, including two decussating and two nondecussating pathways that cannot be identified on conventional structural magnetic resonance imaging (MRI). Diffusion MRI tractography has the potential to trace these subdivisions and is increasingly used to study the RGVP. However, it is not yet known which fiber tracking strategy is most suitable for RGVP reconstruction. In this study, four tractography methods are compared, including constrained spherical deconvolution (CSD) based probabilistic (iFOD1) and deterministic (SD-Stream) methods, and multi-fiber (UKF-2T) and single-fiber (UKF-1T) unscented Kalman filter (UKF) methods. Experiments use diffusion MRI data from 57 subjects in the Human Connectome Project. The RGVP is identified using regions of interest created by two clinical experts. Quantitative anatomical measurements and expert anatomical judgment are used to assess the advantages and limitations of the four tractography methods. Overall, we conclude that UKF-2T and iFOD1 produce the best RGVP reconstruction results. The iFOD1 method can better quantitatively estimate the percentage of decussating fibers, while the UKF-2T method produces reconstructed RGVPs that are judged to better correspond to the known anatomy and have the highest spatial overlap across subjects. Overall, we find that it is challenging for current tractography methods to both accurately track RGVP fibers that correspond to known anatomy and produce an approximately correct percentage of decussating fibers. We suggest that future algorithm development for RGVP tractography should take consideration of both of these two points.

Synchrotron X-ray micro-CT as a validation dataset for diffusion MRI in whole mouse brain

PURPOSE

To introduce synchrotron X-ray microcomputed tomography (microCT) and demonstrate its use as a natively isotropic, nondestructive, 3D validation modality for diffusion MRI in whole, fixed mouse brain.

METHODS

Postmortem diffusion MRI and microCT data were acquired of the same whole mouse brain. Diffusion data were processed using constrained spherical deconvolution. Synchrotron data were acquired at an isotropic voxel size of 1.17 μm. Structure tensor analysis was used to calculate fiber orientation distribution functions from the microCT data. A pipeline was developed to spatially register the 2 datasets in order to perform qualitative comparisons of fiber geometries represented by fiber orientation distribution functions. Fiber orientations from both modalities were used to perform whole-brain deterministic tractography to demonstrate validation of long-range white matter pathways.

RESULTS

Fiber orientation distribution functions were able to be extracted throughout the entire microCT dataset, with spatial registration to diffusion MRI simplified due to the whole brain extent of the microCT data. Fiber orientations and tract pathways showed good agreement between modalities.

CONCLUSION

Synchrotron microCT is a potentially valuable new tool for future multi-scale diffusion MRI validation studies, providing comparable value to optical histology validation methods while addressing some key limitations in data acquisition and ease of processing.

Diffusion MRI tractography for neurosurgery: the basics, current state, technical reliability and challenges

Diffusion magnetic resonance imaging (dMRI) tractography is currently the only imaging technique that allows for non-invasive delineation and visualisation of white matter (WM) tractsin vivo,prompting rapid advances in related fields of brain MRI research in recent years. One of its major clinical applications is for pre-surgical planning and intraoperative image guidance in neurosurgery, where knowledge about the location of WM tracts nearby the surgical target can be helpful to guide surgical resection and optimise post-surgical outcomes. Surgical injuries to these WM tracts can lead to permanent neurological and functional deficits, making the accuracy of tractography reconstructions paramount. The quality of dMRI tractography is influenced by many modifiable factors, ranging from MRI data acquisition through to the post-processing of tractography output, with the potential of error propagation based on decisions made at each and subsequent processing steps. Research over the last 25 years has significantly improved the anatomical accuracy of tractography. An updated review about tractography methodology in the context of neurosurgery is now timely given the thriving research activities in dMRI, to ensure more appropriate applications in the clinical neurosurgical realm. This article aims to review the dMRI physics, and tractography methodologies, highlighting recent advances to provide the key concepts of tractography-informed neurosurgery, with a focus on the general considerations, the current state of practice, technical challenges, potential advances, and future demands to this field.

Structure, function and connectivity fingerprints of the frontal eye field versus the inferior frontal junction: A comprehensive comparison

The human prefrontal cortex contains two prominent areas, the frontal eye field and the inferior frontal junction, that are crucially involved in the orchestrating functions of attention, working memory and cognitive control. Motivated by comparative evidence in non-human primates, we review the human neuroimaging literature, suggesting that the functions of these regions can be clearly dissociated. We found remarkable differences in how these regions relate to sensory domains and visual topography, top-down and bottom-up spatial attention, spatial versus non-spatial (i.e., feature- and object-based) attention and working memory and, finally, the multiple-demand system. Functional magnetic resonance imaging (fMRI) studies using multivariate pattern analysis reveal the selectivity of the frontal eye field and inferior frontal junction to spatial and non-spatial information, respectively. The analysis of functional and effective connectivity provides evidence of the modulation of the activity in downstream visual areas from the frontal eye field and inferior frontal junction and sheds light on their reciprocal influences. We therefore suggest that future studies should aim at disentangling more explicitly the role of these regions in the control of spatial and non-spatial selection. We propose that the analysis of the structural and functional connectivity (i.e., the connectivity fingerprints) of the frontal eye field and inferior frontal junction may be used to further characterize their involvement in a spatial ('where') and a non-spatial ('what') network, respectively, highlighting segregated brain networks that allow biasing visual selection and working memory performance to support goal-driven behaviour.

Enhancing the estimation of fiber orientation distributions using convolutional neural networks

Local fiber orientation distributions (FODs) can be computed from diffusion magnetic resonance imaging (dMRI). The accuracy and ability of FODs to resolve complex fiber configurations benefits from acquisition protocols that sample a high number of gradient directions, a high maximum b-value, and multiple b-values. However, acquisition time and scanners that follow these standards are limited in clinical settings, often resulting in dMRI acquired at a single shell (single b-value). In this work, we learn improved FODs from clinically acquired dMRI. We evaluate patch-based 3D convolutional neural networks (CNNs) on their ability to regress multi-shell FODs from single-shell FODs, using constrained spherical deconvolution (CSD). We evaluate U-Net and High-Resolution Network (HighResNet) 3D CNN architectures on data from the Human Connectome Project and an in-house dataset. We evaluate how well each CNN can resolve FODs 1) when training and testing on datasets with the same dMRI acquisition protocol; 2) when testing on a dataset with a different dMRI acquisition protocol than used to train the CNN; and 3) when testing on a dataset with a fewer number of gradient directions than used to train the CNN. This work is a step towards more accurate FOD estimation in time- and resource-limited clinical environments.