What do brain endocasts tell us? A comparative analysis of the accuracy of sulcal identification by experts and perspectives in palaeoanthropology.

Palaeoneurology is a complex field as the object of study, the brain, does not fossilize. Studies rely therefore on the (brain) endocranial cast (often named endocast), the only available and reliable proxy for brain shape, size and details of surface. However, researchers debate whether or not specific marks found on endocasts correspond reliably to particular sulci and/or gyri of the brain that were imprinted in the braincase. The aim of this study is to measure the accuracy of sulcal identification through an experiment that reproduces the conditions that palaeoneurologists face when working with hominin endocasts. We asked 14 experts to manually identify well-known foldings in a proxy endocast that was obtained from an MRI of an actual in vivo Homo sapiens head. We observe clear differences in the results when comparing the non-corrected labels (the original labels proposed by each expert) with the corrected labels. This result illustrates that trying to reconstruct a sulcus following the very general known shape/position in the literature or from a mean specimen may induce a bias when looking at an endocast and trying to follow the marks observed there. We also observe that the identification of sulci appears to be better in the lower part of the endocast compared to the upper part. The results concerning specific anatomical traits have implications for highly debated topics in palaeoanthropology. Endocranial description of fossil specimens should in the future consider the variation in position and shape of sulci in addition to using models of mean brain shape. Moreover, it is clear from this study that researchers can perceive sulcal imprints with reasonably high accuracy, but their correct identification and labelling remains a challenge, particularly when dealing with extinct species for which we lack direct knowledge of the brain.

"Plis de passage" Deserve a Role in Models of the Cortical Folding Process.

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.

Browsing Multiple Subjects When the Atlas Adaptation Cannot Be Achieved via a Warping Strategy.

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.

Increased and Decreased Superficial White Matter Structural Connectivity in Schizophrenia and Bipolar Disorder.

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.

Fast Automatic Segmentation of White Matter Streamlines Based on a Multi-Subject Bundle Atlas.

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.

Interactive segmentation of white-matter fibers using a multi-subject atlas.

We present a fast algorithm for automatic segmentation of white matter fibers from tractography datasets based on a multi-subject bundle atlas. We describe a sequential version of the algorithm that runs on a desktop computer CPU, as well as a highly parallel version that uses a Graphics Processing Unit (GPU) as an accelerator. Our sequential implementation runs 270 times faster than a C++/Python implementation of a previous algorithm based on the same segmentation method, and 21 times faster than a highly optimized C version of the same previous algorithm. Our parallelized implementation exploits the multiple computation units and memory hierarchy of the GPU to further speed up the algorithm by a factor of 30 with respect to our sequential code. As a result, the time to segment a subject dataset of 800,000 fibers is reduced from more than 2.5 hours in the Python/C++ code, to less than one second in the GPU version.

GPU-based acceleration of an automatic white matter segmentation algorithm using CUDA.

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.