Beyond Crossing Fibers: Bootstrap Probabilistic Tractography Using Complex Subvoxel Fiber Geometries

Diffusion magnetic resonance imaging fiber tractography is a powerful tool for investigating human white matter connectivity in vivo. However, it is prone to false positive and false negative results, making interpretation of the tractography result difficult. Optimal tractography must begin with an accurate description of the subvoxel white matter fiber structure, includes quantification of the uncertainty in the fiber directions obtained, and quantifies the confidence in each reconstructed fiber tract. This paper presents a novel and comprehensive pipeline for fiber tractography that meets the above requirements. The subvoxel fiber geometry is described in detail using a technique that allows not only for straight crossing fibers but for fibers that curve and splay. This technique is repeatedly performed within a residual bootstrap statistical process in order to efficiently quantify the uncertainty in the subvoxel geometries obtained. A robust connectivity index is defined to quantify the confidence in the reconstructed connections. The tractography pipeline is demonstrated in the human brain

On Quantifying Local Geometric Structures of Fiber Tracts

International audienceIn diffusion MRI, fiber tracts, represented by densely distributed 3D curves, can be estimated from diffusion weighted images using tractography. The spatial geometric structure of white matter fiber tracts is known to be complex in human brain, but it carries intrinsic information of human brain. In this paper, inspired by studies of liquid crystals, we propose tract-based director field analysis (tDFA) with total six rotationally invariant scalar indices to quantify local geometric structures of fiber tracts. The contributions of tDFA include: 1) We propose orientational order (OO) and orientational dispersion (OD) indices to quantify the degree of alignment and dispersion of fiber tracts; 2) We define the local orthogonal frame for a set of unoriented curves, which is proved to be a generalization of the Frenet frame defined for a single oriented curve; 3) With the local orthogonal frame, we propose splay, bend, and twist indices to quantify three types of orientational distortion of local fiber tracts, and a total distortion index to describe distortions of all three types. The proposed tDFA for fiber tracts is a generalization of the voxel-based DFA (vDFA) which was recently proposed for a spherical function field (i.e., an ODF field). To our knowledge, this is the first work to quantify orientational distortion (splay, bend, twist, and total distortion) of fiber tracts. Experiments show that the proposed scalar indices are useful descriptors of local geometric structures to visualize and analyze fiber tracts

Ball and rackets: inferring fiber fanning from diffusion-weighted MRI

A number of methods have been proposed for resolving crossing fibers from diffusion-weighted (DW) MRI. However, other complex fiber geometries have drawn minimal attention. In this study, we focus on fiber orientation dispersion induced by within-voxel fanning. We use a multi-compartment, model-based approach to estimate fiber dispersion. Bingham distributions are employed to represent continuous distributions of fiber orientations, centered around a main orientation, and capturing anisotropic dispersion. We evaluate the accuracy of the model for different simulated fanning geometries, under different acquisition protocols and we illustrate the high SNR and angular resolution needs. We also perform a qualitative comparison between our parametric approach and five popular non-parametric techniques that are based on orientation distribution functions (ODFs). This comparison illustrates how the same underlying geometry can be depicted by different methods. We apply the proposed model on high-quality, post-mortem macaque data and present whole-brain maps of fiber dispersion, as well as exquisite details on the local anatomy of fiber distributions in various white matter regions

Local white matter geometry from diffusion tensor gradients

We introduce a mathematical framework for computing geometrical properties of white matter fibres directly from diffusion tensor fields. The key idea is to isolate the portion of the gradient of the tensor field corresponding to local variation in tensor orientation, and to project it onto a coordinate frame of tensor eigenvectors. The resulting eigenframe-centered representation then makes it possible to define scalar indices (or measures) that describe the local white matter geometry directly from the diffusion tensor field and its gradient, without requiring prior tractography. We derive new scalar indices of (1) fibre dispersion and (2) fibre curving, and we demonstrate them on synthetic and in vivo data. Finally, we illustrate their applicability to a group study on schizophrenia

Localized abnormalities in the cingulum bundle in patients with schizophrenia: A Diffusion Tensor tractography study

The cingulum bundle (CB) connects gray matter structures of the limbic system and as such has been implicated in the etiology of schizophrenia. There is growing evidence to suggest that the CB is actually comprised of a conglomeration of discrete sub-connections. The present study aimed to use Diffusion Tensor tractography to subdivide the CB into its constituent sub-connections, and to investigate the structural integrity of these sub-connections in patients with schizophrenia and matched healthy controls. Diffusion Tensor Imaging scans were acquired from 24 patients diagnosed with chronic schizophrenia and 26 matched healthy controls. Deterministic tractography was used in conjunction with FreeSurfer-based regions-of-interest to subdivide the CB into 5 sub-connections (I1 to I5). The patients with schizophrenia exhibited subnormal levels of FA in two cingulum sub-connections, specifically the fibers connecting the rostral and caudal anterior cingulate gyrus (I1) and the fibers connecting the isthmus of the cingulate with the parahippocampal cortex (I4). Furthermore, while FA in the I1 sub-connection was correlated with the severity of patients' positive symptoms (specifically hallucinations and delusions), FA in the I4 sub-connection was correlated with the severity of patients' negative symptoms (specifically affective flattening and anhedonia/asociality). These results support the notion that the CB is a conglomeration of structurally interconnected yet functionally distinct sub-connections, of which only a subset are abnormal in patients with schizophrenia. Furthermore, while acknowledging the fact that the present study only investigated the CB, these results suggest that the positive and negative symptoms of schizophrenia may have distinct neurobiological underpinnings

Brain connectivity using geodesics in HARDI

International audienceWe develop an algorithm for brain connectivity assessment using geodesics in HARDI (high angular resolution diffusion imaging). We propose to recast the problem of finding fibers bundles and connectivity maps to the calculation of shortest paths on a Riemannian manifold defined from fiber ODFs computed from HARDI measurements. Several experiments on real data show that out method is able to segment fibers bundles that are not easily recovered by other existing methods

Improved tractography using asymmetric fibre orientation distributions

Diffusion MRI allows us to make inferences on the structural organisation of the brain by mapping water diffusion to white matter microstructure. However, such a mapping is generally ill-defined; for instance, diffusion measurements are antipodally symmetric (diffusion along x and –x are equal), whereas the distribution of fibre orientations within a voxel is generally not symmetric. Therefore, different sub-voxel patterns such as crossing, fanning, or sharp bending, cannot be distinguished by fitting a voxel-wise model to the signal. However, asymmetric fibre patterns can potentially be distinguished once spatial information from neighbouring voxels is taken into account. We propose a neighbourhood-constrained spherical deconvolution approach that is capable of inferring asymmetric fibre orientation distributions (A-fods). Importantly, we further design and implement a tractography algorithm that utilises the estimated A-fods, since the commonly used streamline tractography paradigm cannot directly take advantage of the new information. We assess performance using ultra-high resolution histology data where we can compare true orientation distributions against sub-voxel fibre patterns estimated from down-sampled data. Finally, we explore the benefits of A-fods-based tractography using in vivo data by evaluating agreement of tractography predictions with connectivity estimates made using different in-vivo modalities. The proposed approach can reliably estimate complex fibre patterns such as sharp bending and fanning, which voxel-wise approaches cannot estimate. Moreover, histology-based and in-vivo results show that the new framework allows more accurate tractography and reconstruction of maps quantifying (symmetric and asymmetric) fibre complexity

Structural connectivity of cytoarchitectonically distinct human left temporal pole subregions: a diffusion MRI tractography study

The temporal pole (TP) is considered one of the major paralimbic cortical regions, and is involved in a variety of functions such as sensory perception, emotion, semantic processing, and social cognition. Based on differences in cytoarchitecture, the TP can be further subdivided into smaller regions (dorsal, ventrolateral and ventromedial), each forming key nodes of distinct functional networks. However, the brain structural connectivity profile of TP subregions is not fully clarified. Using diffusion MRI data in a set of 31 healthy subjects, we aimed to elucidate the comprehensive structural connectivity of three cytoarchitectonically distinct TP subregions. Diffusion tensor imaging (DTI) analysis suggested that major association fiber pathways such as the inferior longitudinal, middle longitudinal, arcuate, and uncinate fasciculi provide structural connectivity to the TP. Further analysis suggested partially overlapping yet still distinct structural connectivity patterns across the TP subregions. Specifically, the dorsal subregion is strongly connected with wide areas in the parietal lobe, the ventrolateral subregion with areas including constituents of the default-semantic network, and the ventromedial subregion with limbic and paralimbic areas. Our results suggest the involvement of the TP in a set of extensive but distinct networks of cortical regions, consistent with its functional roles

Unsupervised White Matter Fiber Clustering and Tract Probability Map Generation: Applications of a Gaussian Process framework for White Matter Fibers

With the increasing importance of fiber tracking in diffusion tensor images for clinical needs, there has been a growing demand for an objective mathematical framework to perform quantitative analysis of white matter fiber bundles incorporating their underlying physical significance. This paper presents such a novel mathematical framework that facilitates mathematical operations between tracts using an inner product based on Gaussian processes, between fibers which span a metric space. This metric facilitates combination of fiber tracts, rendering operations like tract membership to a bundle or bundle similarity simple. Based on this framework, we have designed an automated unsupervised atlas-based clustering method that does not require manual initialization nor an a priori knowledge of the number of clusters. Quantitative analysis can now be performed on the clustered tract volumes across subjects thereby avoiding the need for point parametrization of these fibers, or the use of medial or envelope representations as in previous work. Experiments on synthetic data demonstrate the mathematical operations. Subsequently, the applicability of the unsupervised clustering framework has been demonstrated on a 21 subject dataset

Gray matter alterations in early aging: A diffusion magnetic resonance imaging study

Many studies have observed altered neurofunctional and structural organization in the aging brain. These observations from functional neuroimaging studies show a shift in brain activity from the posterior to the anterior regions with aging (PASA model), as well as a decrease in cortical thickness, which is more pronounced in the frontal lobe followed by the parietal, occipital, and temporal lobes (retrogenesis model). However, very little work has been done using diffusion MRI (dMRI) with respect to examining the structural tissue alterations underlying these neurofunctional changes in the gray matter. Thus, for the first time, we propose to examine gray matter changes using diffusion MRI in the context of aging. In this work, we propose a novel dMRI based measure of gray matter “heterogeneity” that elucidates these functional and structural models (PASA and retrogenesis) of aging from the viewpoint of diffusion MRI. In a cohort of 85 subjects (all males, ages 15–55 years), we show very high correlation between age and “heterogeneity” (a measure of structural layout of tissue in a region-of-interest) in specific brain regions. We examine gray matter alterations by grouping brain regions into anatomical lobes as well as functional zones. Our findings from dMRI data connects the functional and structural domains and confirms the “retrogenesis” hypothesis of gray matter alterations while lending support to the neurofunctional PASA model of aging in addition to showing the preservation of paralimbic areas during healthy aging