Superficial white matter bundle atlas based on hierarchical fiber clustering over probabilistic tractography data.

The study of short association fibers is still an incomplete task due to their higher inter-subject variability and the smaller size of this kind of fibers in comparison to known long association bundles. However, their description is essential to understand human brain dysfunction and better characterize the human brain connectome. In this work, we present a multi-subject atlas of short association fibers, which was computed using a superficial white matter bundle identification method based on fiber clustering. To create the atlas, we used probabilistic tractography from one hundred subjects from the HCP database, aligned with non-linear registration. The method starts with an intra-subject clustering of short fibers (30-85 mm). Based on a cortical atlas, the intra-subject cluster centroids from all subjects are segmented to identify the centroids connecting each region of interest (ROI) of the atlas. To reduce computational load, the centroids from each ROI group are randomly separated into ten subgroups. Then, an inter-subject hierarchical clustering is applied to each centroid subgroup, followed by a second level of clustering to select the most-reproducible clusters across subjects for each ROI group. Finally, the clusters are labeled according to the regions that they connect, and clustered to create the final bundle atlas. The resulting atlas is composed of 525 bundles of superficial short association fibers along the whole brain, with 384 bundles connecting pairs of different ROIs and 141 bundles connecting portions of the same ROI. The reproducibility of the bundles was verified using automatic segmentation on three different tractogram databases. Results for deterministic and probabilistic tractography data show high reproducibility, especially for probabilistic tractography in HCP data. In comparison to previous work, our atlas features a higher number of bundles and greater cortical surface coverage.

Superficial white matter: A review on the dMRI analysis methods and applications.

The mapping of human brain connections is still an on going task. Unlike deep white matter (DWM), which has been extensively studied and well documented, superficial white matter (SWM) has been often left aside. Improving our understanding of the SWM is an important goal for a better understanding of the brain network and its relation to several pathologies. The shape and localization of these short bundles present a high variability across subjects. Furthermore, the small diameter of most superficial bundles and partial volume effects induced by their proximity to the cortex leads to complex tratography issues. Therefore, the mapping of SWM bundles and the use of the resulting atlases for clinical studies requiere dedicated methodologies that are reviewed in this paper.

Reproducibility of superficial white matter tracts using diffusion-weighted imaging tractography.

Human brain connection map is far from being complete. In particular the study of the superficial white matter (SWM) is an unachieved task. Its description is essential for the understanding of human brain function and the study of pathogenesis triggered by abnormal connectivity. In this work we automatically created a multi-subject atlas of SWM diffusion-based bundles of the whole brain. For each subject, the complete cortico-cortical tractogram is first split into sub-tractograms connecting pairs of gyri. Then intra-subject shape-based fiber clustering performs compression of each sub-tractogram into a set of bundles. Proceeding further with shape-based clustering provides a match of the bundles across subjects. Bundles found in most of the subjects are instantiated in the atlas. To increase robustness, this procedure was performed with two independent groups of subjects, in order to discard bundles without match across the two independent atlases. Finally, the resulting intersection atlas was projected on a third independent group of subjects in order to filter out bundles without reproducible and reliable projection. The final multi-subject diffusion-based U-fiber atlas is composed of 100 bundles in total, 50 per hemisphere, from which 35 are common to both hemispheres.

Short fiber bundle filtering and test-retest reproducibility of the Superficial White Matter.

In recent years, there has been a growing interest in studying the Superficial White Matter (SWM). The SWM consists of short association fibers connecting near giry of the cortex, with a complex organization due to their close relationship with the cortical folding patterns. Therefore, their segmentation from dMRI tractography datasets requires dedicated methodologies to identify the main fiber bundle shape and deal with spurious fibers. This paper presents an enhanced short fiber bundle segmentation based on a SWM bundle atlas and the filtering of noisy fibers. The method was tuned and evaluated over HCP test-retest probabilistic tractography datasets (44 subjects). We propose four fiber bundle filters to remove spurious fibers. Furthermore, we include the identification of the main fiber fascicle to obtain well-defined fiber bundles. First, we identified four main bundle shapes in the SWM atlas, and performed a filter tuning in a subset of 28 subjects. The filter based on the Convex Hull provided the highest similarity between corresponding test-retest fiber bundles. Subsequently, we applied the best filter in the 16 remaining subjects for all atlas bundles, showing that filtered fiber bundles significantly improve test-retest reproducibility indices when removing between ten and twenty percent of the fibers. Additionally, we applied the bundle segmentation with and without filtering to the ABIDE-II database. The fiber bundle filtering allowed us to obtain a higher number of bundles with significant differences in fractional anisotropy, mean diffusivity, and radial diffusivity of Autism Spectrum Disorder patients relative to controls.

Phybers: a package for brain tractography analysis.

We present a Python library (Phybers) for analyzing brain tractography data. Tractography datasets contain streamlines (also called fibers) composed of 3D points representing the main white matter pathways. Several algorithms have been proposed to analyze this data, including clustering, segmentation, and visualization methods. The manipulation of tractography data is not straightforward due to the geometrical complexity of the streamlines, the file format, and the size of the datasets, which may contain millions of fibers. Hence, we collected and structured state-of-the-art methods for the analysis of tractography and packed them into a Python library, to integrate and share tools for tractography analysis. Due to the high computational requirements, the most demanding modules were implemented in C/C++. Available functions include brain Bundle Segmentation (FiberSeg), Hierarchical Fiber Clustering (HClust), Fast Fiber Clustering (FFClust), normalization to a reference coordinate system, fiber sampling, calculation of intersection between sets of brain fibers, tools for cluster filtering, calculation of measures from clusters, and fiber visualization. The library tools were structured into four principal modules: Segmentation, Clustering, Utils, and Visualization (Fibervis). Phybers is freely available on a GitHub repository under the GNU public license for non-commercial use and open-source development, which provides sample data and extensive documentation. In addition, the library can be easily installed on both Windows and Ubuntu operating systems through the pip library.

Fiber Clustering Acceleration With a Modified Kmeans++ Algorithm Using Data Parallelism.

Fiber clustering methods are typically used in brain research to study the organization of white matter bundles from large diffusion MRI tractography datasets. These methods enable exploratory bundle inspection using visualization and other methods that require identifying brain white matter structures in individuals or a population. Some applications, such as real-time visualization and inter-subject clustering, need fast and high-quality intra-subject clustering algorithms. This work proposes a parallel algorithm using a General Purpose Graphics Processing Unit (GPGPU) for fiber clustering based on the FFClust algorithm. The proposed GPGPU implementation exploits data parallelism using both multicore and GPU fine-grained parallelism present in commodity architectures, including current laptops and desktop computers. Our approach implements all FFClust steps in parallel, improving execution times in all of them. In addition, our parallel approach includes a parallel Kmeans++ algorithm implementation and defines a new variant of Kmeans++ to reduce the impact of choosing outliers as initial centroids. The results show that our approach provides clustering quality results very similar to FFClust, and it requires an execution time of 3.5 s for processing about a million fibers, achieving a speedup of 11.5 times compared to FFClust.

A stringent fiber distance measure for dMRI tractography clustering and segmentation.

Most analysis and segmentation methods for diffusion MRI tractography datasets require a fiber distance measure able to determine the similarity between a pair of fibers. We present a stringent fiber distance measure able to perform a good discrimination between fiber shapes and lengths. It uses three terms: (i) a fiber maximum Euclidean distance, (ii) a fiber shape distance, and (iii) a fiber length distance. The distance was evaluated applying a hierarchical clustering of fibers connecting the pre-and post-central gyri of a subject. Results where compared with other known fiber distance measures. A better sensitivity to differences in fiber shape and length was found for the proposed distance. This will be very useful for the detailed study and description of white matter bundles. Known bundles will be better decomposed into sub-bundles, with more precision on the bundle shape and on the regions connected by the fibers. For short association bundles, this distance will be a real improvement, as even the most stringent distance used until now shows some limitations when evaluating the similarity of these fibers.

Clustering of Whole-Brain White Matter Short Association Bundles Using HARDI Data.

Human brain connectivity is extremely complex and variable across subjects. While long association and projection bundles are stable and have been deeply studied, short association bundles present higher intersubject variability, and few studies have been carried out to adequately describe the structure, shape, and reproducibility of these bundles. However, their analysis is crucial to understand brain function and better characterize the human connectome. In this study, we propose an automatic method to identify reproducible short association bundles of the superficial white matter, based on intersubject hierarchical clustering. The method is applied to the whole brain and finds representative clusters of similar fibers belonging to a group of subjects, according to a distance metric between fibers. We experimented with both affine and non-linear registrations and, due to better reproducibility, chose the results obtained from non-linear registration. Once the clusters are calculated, our method performs automatic labeling of the most stable connections based on individual cortical parcellations. We compare results between two independent groups of subjects from a HARDI database to generate reproducible connections for the creation of an atlas. To perform a better validation of the results, we used a bagging strategy that uses pairs of groups of 27 subjects from a database of 74 subjects. The result is an atlas with 44 bundles in the left hemisphere and 49 in the right hemisphere, of which 33 bundles are found in both hemispheres. Finally, we use the atlas to automatically segment 78 new subjects from a different HARDI database and to analyze stability and lateralization results.