ABLE: Automated Brain Lines Extraction Based on Laplacian Surface Collapse.

The archetypical folded shape of the human cortex has been a long-standing topic for neuroscientific research. Nevertheless, the accurate neuroanatomical segmentation of sulci remains a challenge. Part of the problem is the uncertainty of where a sulcus transitions into a gyrus and vice versa. This problem can be avoided by focusing on sulcal fundi and gyral crowns, which represent the topological opposites of cortical folding. We present Automated Brain Lines Extraction (ABLE), a method based on Laplacian surface collapse to reliably segment sulcal fundi and gyral crown lines. ABLE is built to work on standard FreeSurfer outputs and eludes the delineation of anastomotic sulci while maintaining sulcal fundi lines that traverse the regions with the highest depth and curvature. First, it segments the cortex into gyral and sulcal surfaces; then, each surface is spatially filtered. A Laplacian-collapse-based algorithm is applied to obtain a thinned representation of the surfaces. This surface is then used for careful detection of the endpoints of the lines. Finally, sulcal fundi and gyral crown lines are obtained by eroding the surfaces while preserving the connectivity between the endpoints. The method is validated by comparing ABLE with three other sulcal extraction methods using the Human Connectome Project (HCP) test-retest database to assess the reproducibility of the different tools. The results confirm ABLE as a reliable method for obtaining sulcal lines with an accurate representation of the sulcal topology while ignoring anastomotic branches and the overestimation of the sulcal fundi lines. ABLE is publicly available via https://github.com/HGGM-LIM/ABLE .This work was supported by the project exAScale ProgramIng models for extreme Data procEssing (ASPIDE), that has received funding from the European Union’s Horizon 2020 research and innovation program under grant agreement No 801091. This work has received funding from “la Caixa” Foundation under the project code LCF/PR/HR19/52160001. Susanna Carmona funded by Instituto de Salud Carlos III, co-funded by European Social Fund “Investing in your future” (Miguel Servet Type I research contract CP16/00096). The CNIC is supported by the Instituto de Salud Carlos III (ISCIII), the Ministerio de Ciencia e Innovación (MCIN) and the Pro CNIC Foundation, and is a Severo Ochoa Center of Excellence (SEV-2015-0505). Yasser Alemán-Gómez is supported by the Swiss National Science Foundation (185897) and the National Center of Competence in Research (NCCR) SYNAPSY - The Synaptic Bases of Mental Diseases, funded as well by the Swiss National Science Foundation (51AU40-1257).S

Surface-Based tools for Characterizing the Human Brain Cortical Morphology

Tesis por compendio de publicacionesThe cortex of the human brain is highly convoluted. These characteristic convolutions present advantages over lissencephalic brains. For instance, gyrification allows an expansion of cortical surface area without significantly increasing the cranial volume, thus facilitating the pass of the head through the birth channel. Studying the human brain’s cortical morphology and the processes leading to the cortical folds has been critical for an increased understanding of the pathological processes driving psychiatric disorders such as schizophrenia, bipolar disorders, autism, or major depression. Furthermore, charting the normal developmental changes in cortical morphology during adolescence or aging can be of great importance for detecting deviances that may be precursors for pathology. However, the exact mechanisms that push cortical folding remain largely unknown. The accurate characterization of the neurodevelopment processes is challenging. Multiple mechanisms co-occur at a molecular or cellular level and can only be studied through the analysis of ex-vivo samples, usually of animal models. Magnetic Resonance Imaging can partially fill the breach, allowing the portrayal of the macroscopic processes surfacing on in-vivo samples. Different metrics have been defined to measure cortical structure to describe the brain’s morphological changes and infer the associated microstructural events. Metrics such as cortical thickness, surface area, or cortical volume help establish a relation between the measured voxels on a magnetic resonance image and the underlying biological processes. However, the existing methods present limitations or room for improvement. Methods extracting the lines representing the gyral and sulcal morphology tend to over- or underestimate the total length. These lines can provide important information about how sulcal and gyral regions function differently due to their distinctive ontogenesis. Nevertheless, some methods label every small fold on the cortical surface as a sulcal fundus, thus losing the perspective of lines that travel through the deeper zones of a sulcal basin. On the other hand, some methods are too restrictive, labeling sulcal fundi only for a bunch of primary folds. To overcome this issue, we have proposed a Laplacian-collapse-based algorithm that can delineate the lines traversing the top regions of the gyri and the fundi of the sulci avoiding anastomotic sulci. For this, the cortex, represented as a 3D surface, is segmented into gyral and sulcal surfaces attending to the curvature and depth at every point of the mesh. Each resulting surface is spatially filtered, smoothing the boundaries. Then, a Laplacian-collapse-based algorithm is applied to obtain a thinned representation of the morphology of each structure. These thin curves are processed to detect where the extremities or endpoints lie. Finally, sulcal fundi and gyral crown lines are obtained by eroding the surfaces while preserving the structure topology and connectivity between the endpoints. The assessment of the presented algorithm showed that the labeled sulcal lines were close to the proposed ground truth length values while crossing through the deeper (and more curved) regions. The tool also obtained reproducibility scores better or similar to those of previous algorithms. A second limitation of the existing metrics concerns the measurement of sulcal width. This metric, understood as the physical distance between the points on opposite sulcal banks, can come in handy in detecting cortical flattening or complementing the information provided by cortical thickness, gyrification index, or such features. Nevertheless, existing methods only provided averaged measurements for different predefined sulcal regions, greatly restricting the possibilities of sulcal width and ignoring the intra-region variability. Regarding this, we developed a method that estimates the distance from each sulcal point in the cortex to its corresponding opposite, thus providing a per-vertex map of the physical sulcal distances. For this, the cortical surface is sampled at different depth levels, detecting the points where the sulcal banks change. The points corresponding to each sulcal wall are matched with the closest point on a different one. The distance between those points is the sulcal width. The algorithm was validated against a simulated sulcus that resembles a simple fold. Then the tool was used on a real dataset and compared against two widely-used sulcal width estimation methods, averaging the proposed algorithm’s values into the same region definition those reference tools use. The resulting values were similar for the proposed and the reference methods, thus demonstrating the algorithm’s accuracy. Finally, both algorithms were tested on a real aging population dataset to prove the methods’ potential in a use-case scenario. The main idea was to elucidate fine-grained morphological changes in the human cortex with aging by conducting three analyses: a comparison of the age-dependencies of cortical thickness in gyral and sulcal lines, an analysis of how the sulcal and gyral length changes with age, and a vertex-wise study of sulcal width and cortical thickness. These analyses showed a general flattening of the cortex with aging, with interesting findings such as a differential age-dependency of thickness thinning in the sulcal and gyral regions. By demonstrating that our method can detect this difference, our results can pave the way for future in vivo studies focusing on macro- and microscopic changes specific to gyri or sulci. Our method can generate new brain-based biomarkers specific to sulci and gyri, and these can be used on large samples to establish normative models to which patients can be compared. In parallel, the vertex-wise analyses show that sulcal width is very sensitive to changes during aging, independent of cortical thickness. This corroborates the concept of sulcal width as a metric that explains, in the least, the unique variance of morphology not fully captured by existing metrics. Our method allows for sulcal width vertex-wise analyses that were not possible previously, potentially changing our understanding of how changes in sulcal width shape cortical morphology. In conclusion, this thesis presents two new tools, open source and publicly available, for estimating cortical surface-based morphometrics. The methods have been validated and assessed against existing algorithms. They have also been tested on a real dataset, providing new, exciting insights into cortical morphology and showing their potential for defining innovative biomarkers.Programa de Doctorado en Ciencia y Tecnología Biomédica por la Universidad Carlos III de MadridPresidente: Juan Domingo Gispert López.- Secretario: Norberto Malpica González de Vega.- Vocal: Gemma Cristina Monté Rubi

A bipolar taxonomy of adult human brain sulcal morphology related to timing of fetal sulcation and trans-sulcal gene expression gradients

We developed a computational pipeline (now provided as a resource) for measuring morphological similarity between cortical surface sulci to construct a sulcal phenotype network (SPN) from each magnetic resonance imaging (MRI) scan in an adult cohort (N=34,725; 45-82 years). Networks estimated from pairwise similarities of 40 sulci on 5 morphological metrics comprised two clusters of sulci, represented also by the bipolar distribution of sulci on a linear-to-complex dimension. Linear sulci were more heritable and typically located in unimodal cortex; complex sulci were less heritable and typically located in heteromodal cortex. Aligning these results with an independent fetal brain MRI cohort (N=228; 21-36 gestational weeks), we found that linear sulci formed earlier, and the earliest and latest-forming sulci had the least between-adult variation. Using high-resolution maps of cortical gene expression, we found that linear sulcation is mechanistically underpinned by trans-sulcal gene expression gradients enriched for developmental processes.</p

Anatomo-functional correspondence in the superior temporal sulcus

The superior temporal sulcus (STS) is an intriguing region both for its complex anatomy and for the multiple functions that it hosts. Unfortunately, most studies explored either the functional organization or the anatomy of the STS only. Here, we link these two aspects by investigating anatomo-functional correspondences between the voice-sensitive cortex (Temporal Voice Areas) and the STS depth. To do so, anatomical and functional scans of 116 subjects were processed such as to generate individual surface maps on which both depth and functional voice activity can be analyzed. Individual depth profiles of manually drawn STS and functional profiles from a voice localizer (voice &gt; non-voice) maps were extracted and compared to assess anatomo-functional correspondences. Three major results were obtained: first, the STS exhibits a highly significant rightward depth asymmetry in its middle part. Second, there is an anatomo-functional correspondence between the location of the voice-sensitive peak and the deepest point inside this asymmetrical region bilaterally. Finally, we showed that this correspondence was independent of the gender and, using a machine learning approach, that it existed at the individual level. These findings offer new perspectives for the understanding of anatomo-functional correspondences in this complex cortical region

Cortical Surface Registration and Shape Analysis

A population analysis of human cortical morphometry is critical for insights into brain development or degeneration. Such an analysis allows for investigating sulcal and gyral folding patterns. In general, such a population analysis requires both a well-established cortical correspondence and a well-defined quantification of the cortical morphometry. The highly folded and convoluted structures render a reliable and consistent population analysis challenging. Three key challenges have been identified for such an analysis: 1) consistent sulcal landmark extraction from the cortical surface to guide better cortical correspondence, 2) a correspondence establishment for a reliable and stable population analysis, and 3) quantification of the cortical folding in a more reliable and biologically meaningful fashion. The main focus of this dissertation is to develop a fully automatic pipeline that supports a population analysis of local cortical folding changes. My proposed pipeline consists of three novel components I developed to overcome the challenges in the population analysis: 1) automatic sulcal curve extraction for stable/reliable anatomical landmark selection, 2) group-wise registration for establishing cortical shape correspondence across a population with no template selection bias, and 3) quantification of local cortical folding using a novel cortical-shape-adaptive kernel. To evaluate my methodological contributions, I applied all of them in an application to early postnatal brain development. I studied the human cortical morphological development using the proposed quantification of local cortical folding from neonate age to 1 year and 2 years of age, with quantitative developmental assessments. This study revealed a novel pattern of associations between the cortical gyrification and cognitive development.Doctor of Philosoph

Mindboggling morphometry of human brains

Mindboggle (http://mindboggle.info) is an open source brain morphometry platform that takes in preprocessed T1-weighted MRI data and outputs volume, surface, and tabular data containing label, feature, and shape information for further analysis. In this article, we document the software and demonstrate its use in studies of shape variation in healthy and diseased humans. The number of different shape measures and the size of the populations make this the largest and most detailed shape analysis of human brains ever conducted. Brain image morphometry shows great potential for providing much-needed biological markers for diagnosing, tracking, and predicting progression of mental health disorders. Very few software algorithms provide more than measures of volume and cortical thickness, while more subtle shape measures may provide more sensitive and specific biomarkers. Mindboggle computes a variety of (primarily surface-based) shapes: area, volume, thickness, curvature, depth, Laplace-Beltrami spectra, Zernike moments, etc. We evaluate Mindboggle’s algorithms using the largest set of manually labeled, publicly available brain images in the world and compare them against state-of-the-art algorithms where they exist. All data, code, and results of these evaluations are publicly available

Imaging local genetic influences on cortical folding

Recent progress in deciphering mechanisms of human brain cortical folding leave unexplained whether spatially patterned genetic influences contribute to this folding. High-resolution in vivo brain MRI can be used to estimate genetic correlations (covariability due to shared genetic factors) in interregional cortical thickness, and biomechanical studies predict an influence of cortical thickness on folding patterns. However, progress has been hampered because shared genetic influences related to folding patterns likely operate at a scale that is much more local (cm) than that addressed in prior imaging studies. Here, we develop methodological approaches to examine local genetic influences on cortical thickness and apply these methods to two large, independent samples. We find that such influences are markedly heterogeneous in strength, and in some cortical areas are notably stronger in specific orientations relative to gyri or sulci. The overall, phenotypic local correlation has a significant basis in shared genetic factors and is highly symmetric between left and right cortical hemispheres. Furthermore, the degree of local cortical folding relates systematically with the strength of local correlations, which tends to be higher in gyral crests and lower in sulcal fundi. The relationship between folding and local correlations is stronger in primary sensorimotor areas and weaker in association areas such as prefrontal cortex, consistent with reduced genetic constraints on the structural topology of association cortex. Collectively, our results suggest that patterned genetic influences on cortical thickness, measurable at the scale of in vivo MRI, may be a causal factor in the development of cortical folding

Characterization of U-shape streamline fibers: Methods and applications

Diffusion tensor imaging (DTI), high angular resolution diffusion imaging (HARDI), and Diffusion Spectrum Imaging (DSI) have been widely used in the neuroimaging field to examine the macro-scale fiber connection patterns in the cerebral cortex. However, the topographic and geometric relationships between diffusion imaging derived streamline fiber connection patterns and cortical folding patterns remain largely unknown. This paper specifically identifies and characterizes the U-shapes of diffusion imaging derived streamline fibers via a novel fiber clustering framework and examines their co-localization patterns with cortical sulci based on DTI, HARDI, and DSI datasets of human, chimpanzee and macaque brains. We verified the presence of these U-shaped streamline fibers that connect neighboring gyri by coursing around cortical sulci such as the central sulcus, pre-central sulcus, post-central sulcus, superior temporal sulcus, inferior frontal sulcus, and intra-parietal sulcus. This study also verified the existence of U-shape fibers across data modalities (DTI/HARDI/DSI) and primate species (macaque, chimpanzee and human), and suggests that the common pattern of U-shape fibers coursing around sulci is evolutionarily-preserved in cortical architectures