This paper presents an algorithm for fast segmentation of white matter bundles from massive dMRI tractography datasets using a multisubject atlas. We use a distance metric to compare streamlines in a subject dataset to labeled centroids in the atlas, and label them using a per-bundle configurable threshold. In order to reduce segmentation time, the algorithm first preprocesses the data using a simplified distance metric to rapidly discard candidate streamlines in multiple stages, while guaranteeing that no false negatives are produced. The smaller set of remaining streamlines is then segmented using the original metric, thus eliminating any false positives from the preprocessing stage. As a result, a single-thread implementation of the algorithm can segment a dataset of almost 9 million streamlines in less than 6 minutes. Moreover, parallel versions of our algorithm for multicore processors and graphics processing units further reduce the segmentation time to less than 22 seconds and to 5 seconds, respectively. This performance enables the use of the algorithm in truly interactive applications for visualization, analysis, and segmentation of large white matter tractography datasets.

Cortical folding is a hallmark of brain topography whose variability across individuals remains a puzzle. In this paper, we call for an effort to improve our understanding of the pli de passage phenomenon, namely annectant gyri buried in the depth of the main sulci. We suggest that plis de passage could become an interesting benchmark for models of the cortical folding process. As an illustration, we speculate on the link between modern biological models of cortical folding and the development of the Pli de Passage Frontal Moyen (PPFM) in the middle of the central sulcus. For this purpose, we have detected nine interrupted central sulci in the Human Connectome Project dataset, which are used to explore the organization of the hand sensorimotor areas in this rare configuration of the PPFM.

Schizophrenia (SZ) and bipolar disorder (BD) are often conceptualized as “disconnection syndromes,” with substantial evidence of abnormalities in deep white matter tracts, forming the substrates of long-range connectivity, seen in both disorders. However, the study of superficial white matter (SWM) U-shaped short-range tracts remained challenging until recently, although findings from postmortem studies suggest they are likely integral components of SZ and BD neuropathology. This diffusion weighted imaging (DWI) study aimed to investigate SWM microstructure in vivo in both SZ and BD for the first time. We performed whole brain tractography in 31 people with SZ, 32 people with BD and 54 controls using BrainVISA and Connectomist 2.0. Segmentation and labeling of SWM tracts were performed using a novel, comprehensive U-fiber atlas. Analysis of covariances yielded significant generalized fractional anisotropy (gFA) differences for 17 SWM bundles in frontal, parietal, and temporal cortices. Post hoc analyses showed gFA reductions in both patient groups as compared with controls in bundles connecting regions involved in language processing, mood regulation, working memory, and motor function (pars opercularis, insula, anterior cingulate, precentral gyrus). We also found increased gFA in SZ patients in areas overlapping the default mode network (inferior parietal, middle temporal, precuneus), supporting functional hyperconnectivity of this network evidenced in SZ. We thus illustrate that short U-fibers are vulnerable to the pathological processes in major psychiatric illnesses, encouraging improved understanding of their anatomy and function.

Robust spatial alignment of post mortem data and in vivo MRI acquisitions from different ages, especially from the early developmental stages, into standard spaces is still a bottleneck hampering easy comparison with the mainstream neuroimaging results. In this paper, we test a landmark-based spatial normalization strategy as a framework for the seamless integration of any macroscopic dataset in the context of the Human Brain Project (HBP). This strategy stems from an approach called DISCO embedding sulcal constraints in a registration framework used to initialize DARTEL, the widely used spatial normalization approach proposed in the SPM software. We show that this strategy is efficient with a heterogeneous dataset including challenging data as preterm newborns, infants, post mortem histological data and a synthetic atlas computed from averaging the ICBM database, as well as more commonly studied data acquired in vivo in adults. We then describe some perspectives for a research program aiming at improving folding pattern matching for atlas inference in the context of the future HBP's portal.

The deformable atlas paradigm has been at the core of computational anatomy during the last two decades. Spatial normalization is the variant endowing the atlas with a coordinate system used for voxel-based aggregation of images across subjects and studies. This framework has largely contributed to the success of brain mapping. Brain spatial normalization, however, is still ill-posed because of the complexity of the human brain architecture and the lack of architectural landmarks in standard morphological MRI. Multi-atlas strategies have been developed during the last decade to overcome some difficulties in the context of segmentation. A new generation of registration algorithms embedding architectural features inferred for instance from diffusion or functional MRI is on the verge to improve the architectural value of spatial normalization. A better understanding of the architectural meaning of the cortical folding pattern will lead to use some sulci as complementary constraints. Improving the architectural compliance of spatial normalization may impose to relax the diffeomorphic constraint usually underlying atlas warping. A two-level strategy could be designed: in each region, a dictionary of templates of incompatible folding patterns would be collected and matched in a way or another using rare architectural information, while individual subjects would be aligned using diffeomorphisms to the closest template. Manifold learning could help to aggregate subjects according to their morphology. Connectivity-based strategies could emerge as an alternative to deformation-based alignment leading to match the connectomes of the subjects rather than images.

We present an optimized algorithm that performs automatic classification of white matter fibers based on a multi-subject bundle atlas. We implemented a parallel algorithm that improves upon its previous version in both execution time and memory usage. Our new version uses the local memory of each processor, which leads to a reduction in execution time. Hence, it allows the analysis of bigger subject and/or atlas datasets. As a result, the segmentation of a subject of 4,145,000 fibers is reduced from about 14 minutes in the previous version to about 6 minutes, yielding an acceleration of 2.34. In addition, the new algorithm reduces the memory consumption of the previous version by a factor of 0.79.

Brain mapping studies often need to identify brain structures or functional circuits into a set of individual brains. To this end, multiple atlases have been published to represent such structures based on different modalities, subject sets, and techniques. The mainstream approach to exploit these atlases consists in spatially deforming each individual data onto a given atlas using dense deformation fields, which supposes the existence of a continuous mapping between atlases and individuals. However, this continuity is not always verified, and this “iconic” approach has limits. We present in this study an alternative, complementary, “structural” approach, which consists in extracting structures from the individual data, and comparing them without deformation. A “structural atlas” is thus a collection of annotated individual data with a common structure nomenclature. It may be used to characterize structure shape variability across individuals or species, or to train machine learning systems. This study exhibits Anatomist, a powerful structural 3D visualization software dedicated to building, exploring, and editing structural atlases involving a large number of subjects. It has been developed primarily to decipher the cortical folding variability; cortical sulci vary enormously in both size and shape, and some may be missing or have various topologies, which makes iconic approaches inefficient to study them. We, therefore, had to build structural atlases for cortical sulci, and use them to train sulci identification algorithms. Anatomist can display multiple subject data in multiple views, supports all kinds of neuroimaging data, including compound structural object graphs, handles arbitrary coordinate transformation chains between data, and has multiple display features. It is designed as a programming library in both C++ and Python languages, and may be extended or used to build dedicated custom applications. Its generic design makes all the display and structural aspects used to explore the variability of the cortical folding pattern work in other applications, for instance, to browse axonal fiber bundles, deep nuclei, functional activations, or other kinds of cortical parcellations. Multimodal, multi-individual, or inter-species display is supported, and adaptations to large scale screen walls have been developed. These very original features make it a unique viewer for structural atlas browsing.

This paper presents a parallel implementation of an algorithm for automatic segmentation of white matter fibers from tractography data. We execute the algorithm in parallel using a high-end video card with a Graphics Processing Unit (GPU) as a computation accelerator, using the CUDA language. By exploiting the parallelism and the properties of the memory hierarchy available on the GPU, we obtain a speedup in execution time of 33.6 with respect to an optimized sequential version of the algorithm written in C, and of 240 with respect to the original Python/C++ implementation. The execution time is reduced from more than two hours to only 35 seconds for a subject dataset of 800,000 fibers, thus enabling applications that use interactive segmentation and visualization of small to medium-sized tractography datasets.