Mapping Topographic Structure in White Matter Pathways with Level Set Trees

Fiber tractography on diffusion imaging data offers rich potential for describing white matter pathways in the human brain, but characterizing the spatial organization in these large and complex data sets remains a challenge. We show that level set trees---which provide a concise representation of the hierarchical mode structure of probability density functions---offer a statistically-principled framework for visualizing and analyzing topography in fiber streamlines. Using diffusion spectrum imaging data collected on neurologically healthy controls (N=30), we mapped white matter pathways from the cortex into the striatum using a deterministic tractography algorithm that estimates fiber bundles as dimensionless streamlines. Level set trees were used for interactive exploration of patterns in the endpoint distributions of the mapped fiber tracks and an efficient segmentation of the tracks that has empirical accuracy comparable to standard nonparametric clustering methods. We show that level set trees can also be generalized to model pseudo-density functions in order to analyze a broader array of data types, including entire fiber streamlines. Finally, resampling methods show the reliability of the level set tree as a descriptive measure of topographic structure, illustrating its potential as a statistical descriptor in brain imaging analysis. These results highlight the broad applicability of level set trees for visualizing and analyzing high-dimensional data like fiber tractography output

Mapping the Connectome: Multi-Level Analysis of Brain Connectivity

Background and scope The brain contains vast numbers of interconnected neurons that constitute anatomical and functional networks. Structural descriptions of neuronal network elements and connections make up the “connectome ” of the brain (Hagmann, 2005; Sporns et al., 2005; Sporns, 2011), and are important for understanding normal brain function and disease-related dysfunction. A long-standing ambition of the neuroscience community has been to achieve complete connectome maps for the human brain as well as the brains of non-human primates, rodents, and other species (Bohland et al., 2009; Hagmann et al., 2010; Van Essen and Ugurbil, 2012). A wide repertoire of experimental tools is currently available to map neural connectivity at multiple levels, from the tracing of mesoscopic axonal connections and the delineation of white matter tracts (Saleem et al., 2002; Van der Linden et al., 2002; Sporns et al., 2005; Schmahmann et al., 2007; Hagmann et al., 2010), the mappin

Automated multi-subject fiber clustering of mouse brain using dominant sets

Mapping of structural and functional connectivity may provide deeper understanding of brain function and disfunction. Diffusion Magnetic Resonance Imaging (DMRI) is a powerful technique to non-invasively delineate white matter (WM) tracts and to obtain a three-dimensional description of the structural architecture of the brain. However, DMRI tractography methods produce highly multi-dimensional datasets whose interpretation requires advanced analytical tools. Indeed, manual identification of specific neuroanatomical tracts based on prior anatomical knowledge is time-consuming and prone to operator-induced bias. Here we propose an automatic multi-subject fiber clustering method that enables retrieval of group-wise WM fiber bundles. In order to account for variance across subjects, we developed a multi-subject approach based on a method known as Dominant Sets algorithm, via an intra- and cross-subject clustering. The intra-subject step allows us to reduce the complexity of the raw tractography data, thus obtaining homogeneous neuroanatomically-plausible bundles in each diffusion space. The cross-subject step, characterized by a proper space-invariant metric in the original diffusion space, enables the identification of the same WM bundles across multiple subjects without any prior neuroanatomical knowledge. Quantitative analysis was conducted comparing our algorithm with spectral clustering and affinity propagation methods on synthetic dataset. We also performed qualitative analysis on mouse brain tractography retrieving significant WM structures. The approach serves the final goal of detecting WM bundles at a population level, thus paving the way to the study of the WM organization across groups.Mapping of structural and functional connectivity may provide deeper understanding of brain function and disfunction. Diffusion Magnetic Resonance Imaging (DMRI) is a powerful technique to non-invasively delineate white matter (WM) tracts and to obtain a three-dimensional description of the structural architecture of the brain. However, DMRI tractography methods produce highly multi-dimensional datasets whose interpretation requires advanced analytical tools. Indeed, manual identification of specific neuroanatomical tracts based on prior anatomical knowledge is time-consuming and prone to operator-induced bias. Here we propose an automatic multi-subject fiber clustering method that enables retrieval of group-wise WM fiber bundles. In order to account for variance across subjects, we developed a multi-subject approach based on a method known as Dominant Sets algorithm, via an intra-and cross-subject clustering. The intra-subject step allows us to reduce the complexity of the raw tractography data, thus obtaining homogeneous neuroanatomically-plausible bundles in each diffusion space. The cross-subject step, characterized by a proper space-invariant metric in the original diffusion space, enables the identification of the same WM bundles across multiple subjects without any prior neuroanatomical knowledge. Quantitative analysis was conducted comparing our algorithm with spectral clustering and affinity propagation methods on synthetic dataset. We also performed qualitative analysis on mouse brain tractography retrieving significant WM structures. The approach serves the final goal of detecting WM bundles at a population level, thus paving the way to the study of the WM organization across groups

Mapping changes of in vivo connectivity patterns in the human mediodorsal thalamus: correlations with higher cognitive and executive functions

The mediodorsal thalamic nucleus is recognized as an association hub mediating interconnections with mainly the prefrontal cortex. Tracer studies in primates and in vivo diffusion tensor tractography findings in both humans and monkeys confirm its role in relaying networks that connect to the dorsolateral prefrontal, orbitofrontal, frontal medial and cingulate cortex. Our study was designed to use in vivo probabilistic tractography to describe the pathways emerging from or projecting to the mediodorsal nucleus; moreover, to use such information to automatically define subdivisions based on the divergence of remote structural connections. Diffusion tensor MR imaging data of 156 subjects were utilized to perform connectivity-based segmentation of the mediodorsal nucleus by employing a k-means clustering algorithm. Two domains were revealed (medial and lateral) that are separated from each other by a sagittally oriented plane. For each subject, general assessment of cognitive performance by means of the Wechsler Abbreviated Scale of Intelligence and measures of Delis-Kaplan Executive Function System (D-KEFS) test was utilized. Inter-subject variability in terms of connectivity-based cluster sizes was discovered and the relative sizes of the lateral mediodorsal domain correlated with the individuals' performance in the D-KEFS Sorting test (r = 0.232, p = 0.004). Our results show that the connectivity-based parcellation technique applied to the mediodorsal thalamic nucleus delivers a single subject level descriptor of connectional topography; furthermore, we revealed a possible weak interaction between executive performance and the size of the thalamic area from which pathways converge to the lateral prefrontal corte

A Tutorial in Connectome Analysis: Topological and Spatial Features of Brain Networks

High-throughput methods for yielding the set of connections in a neural system, the connectome, are now being developed. This tutorial describes ways to analyze the topological and spatial organization of the connectome at the macroscopic level of connectivity between brain regions as well as the microscopic level of connectivity between neurons. We will describe topological features at three different levels: the local scale of individual nodes, the regional scale of sets of nodes, and the global scale of the complete set of nodes in a network. Such features can be used to characterize components of a network and to compare different networks, e.g. the connectome of patients and control subjects for clinical studies. At the global scale, different types of networks can be distinguished and we will describe Erd\"os-R\'enyi random, scale-free, small-world, modular, and hierarchical archetypes of networks. Finally, the connectome also has a spatial organization and we describe methods for analyzing wiring lengths of neural systems. As an introduction for new researchers in the field of connectome analysis, we discuss the benefits and limitations of each analysis approach.Comment: Neuroimage, in pres

Application of neuroanatomical features to tractography clustering: Neural Features to Fiber Clustering

Diffusion tensor imaging allows unprecedented insight into brain neural connectivity in vivo by allowing reconstruction of neuronal tracts via captured patterns of water diffusion in white matter microstructures. However, tractography algorithms often output hundreds of thousands of fibers, rendering subsequent data analysis intractable. As a remedy, fiber clustering techniques are able to group fibers into dozens of bundles and thus facilitate analyses. Most existing fiber clustering methods rely on geometrical information of fibers, by viewing them as curves in 3D Euclidean space. The important neuroanatomical aspect of fibers, however, is ignored. In this article, the neuroanatomical information of each fiber is encapsulated in the associativity vector, which functions as the unique “fingerprint” of the fiber. Specifically, each entry in the associativity vector describes the relationship between the fiber and a certain anatomical ROI in a fuzzy manner. The value of the entry approaches 1 if the fiber is spatially related to the ROI at high confidence; on the contrary, the value drops closer to 0. The confidence of the ROI is calculated by diffusing the ROI according to the underlying fibers from tractography. In particular, we have adopted the fast marching method for simulation of ROI diffusion. Using the associativity vectors of fibers, we further model fibers as observations sampled from multivariate Gaussian mixtures in the feature space. To group all fibers into relevant major bundles, an expectation-maximization clustering approach is employed. Experimental results indicate that our method results in anatomically meaningful bundles that are highly consistent across subjects