Plis de passage in the superior temporal sulcus: Morphology and local connectivity.

While there is a profusion of functional investigations involving the superior temporal sulcus (STS), our knowledge of the anatomy of this sulcus is still limited by a large individual variability. In particular, an accurate characterization of the "plis de passage" (PPs), annectant gyri inside the fold, is lacking to explain this variability. Performed on 90 subjects of the HCP database, our study revealed that PPs constitute landmarks that can be identified from the geometry of the STS walls. They were found associated with a specific U-shape white-matter connectivity between the two banks of the sulcus, the amount of connectivity being related to the depth of the PPs. These findings raise new hypotheses regarding the spatial organization of PPs, the relation between cortical anatomy and structural connectivity, as well as the possible role of PPs in the regional functional organization.

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 > 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.

Spatial normalization of brain images and beyond.

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.

Deep sulcal landmarks: algorithmic and conceptual improvements in the definition and extraction of sulcal pits.

Recent interest has been growing concerning points of maximum depth within folds, the sulcal pits, that can be used as reliable cortical landmarks. These remarkable points on the cortical surface are defined algorithmically as the outcome of an automatic extraction procedure. The influence of several crucial parameters of the reference technique (Im et al., 2010) has not been evaluated extensively, and no optimization procedure has been proposed so far. Designing an appropriate optimization framework for these parameters is mandatory to guarantee the reproducibility of results across studies and to ensure the feasibility of sulcal pit extraction and analysis on large cohorts. In this work, we propose a framework specifically dedicated to the optimization of the parameters of the method. This optimization framework relies on new measures for better quantifying the reproducibility of the number of sulcal pits per region across individuals, in line with the assumptions of one-to-one correspondence of sulcal roots across individuals which is an explicit aspect of the sulcal roots model (Régis et al., 2005). Our procedure benefits from a combination of improvements, including the use of a convenient sulcal depth estimation and is methodologically sound. Our experiments on two different groups of individuals, with a total of 137 subjects, show an increased reliability across subjects in deeper sulcal pits, as compared to the previous approach, and cover the entire cortical surface, including shallower and more variable folds that were not considered before. The effectiveness of our method ensures the feasibility of a systematic study of sulcal pits on large cohorts. On top of these methodological advances, we quantify the relationship between the reproducibility of the number of sulcal pits per region across individuals and their respective depth and demonstrate the relatively high reproducibility of several pits corresponding to shallower folds. Finally, we report new results regarding the local pit asymmetry, providing evidence that the algorithmic and conceptual approach defended here may contribute to better understanding of the key role of sulcal pits in neuroanatomy.

Benchmark data for sulcal pits extraction algorithms.

This article contains data related to the research article Auzias et al. (2015) [1]. This data can be used as a benchmark for quantitative evaluation of sulcal pits extraction algorithm. In particular, it allows a quantitative comparison with our method, and the assessment of the consistency of the sulcal pits extraction across two well-matched populations.

Model-driven harmonic parameterization of the cortical surface: HIP-HOP.

In the context of inter subject brain surface matching, we present a parameterization of the cortical surface constrained by a model of cortical organization. The parameterization is defined via an harmonic mapping of each hemisphere surface to a rectangular planar domain that integrates a representation of the model. As opposed to previous landmark-based registration methods we do not match folds between individuals but instead optimize the fit between cortical sulci and specific iso-coordinate axis in the model. This strategy overcomes some limitation to sulcus-based registration techniques such as topological variability in sulcal landmarks across subjects. Experiments on 62 subjects with manually traced sulci are presented and compared with the result of the Freesurfer software. The evaluation involves a measure of dispersion of sulci with both angular and area distortions. We show that the model-based strategy can lead to a natural, efficient and very fast (less than 5 min per hemisphere) method for defining inter subjects correspondences. We discuss how this approach also reduces the problems inherent to anatomically defined landmarks and open the way to the investigation of cortical organization through the notion of orientation and alignment of structures across the cortex.

Automatic sulcal line extraction on cortical surfaces using geodesic path density maps.

We present here a method that is designed to automatically extract sulcal lines on the mesh of any cortical surface. The method is based on the definition of a new function, the Geodesic Path Density Map (GPDM), within each sulcal basin (i.e. regions with a negative mean curvature). GPDM indicates at each vertex the likelihood that a shortest path between any two points of the basins boundary goes through that vertex. If the distance used to compute shortest path is anisotropic and constrained by a geometric information such as the depth, the GPDM indicates the likelihood that a vertex belongs to the sulcal line in the basin. An automatic GPDM adaptive thresholding procedure is proposed and sulcal lines are then defined. The process has been validated on a set of 25 subjects by comparing results to the manual segmentation from an expert and showed an average error below 2mm. It is also compared to our previous reference method in the context of inter-subject cortical surface registration and shows an significant improvement in performance.

Observer-independent characterization of sulcal landmarks and depth asymmetry in the central sulcus of the chimpanzee brain.

The central sulcus (CS) divides primary motor and sensory cortex in many mammalian brains. Recent studies have shown that experiential factors can influence the volume and lateralization of the CS in both human and nonhuman primates. In this study, we sought to define specific landmarks and the depth of the CS region corresponding to the motor-hand area of chimpanzees for comparison with humans using a novel, observer independent method applied to sample of 32 magnetic resonance images (MRI) scans. Our results showed that the dorsal-ventral location of the motor-hand region is comparable between humans and chimpanzees, though the depth of the CS was significantly greater in humans compared to chimpanzees. We further found that CS area corresponding to the motor-hand area was significantly larger in the hemisphere contralateral to the chimpanzees preferred hand. The methods employed here offer some potential advantages over traditional region-of-interest in the comparative study of cortical organization and gyrification in primates and are discussed.

Model-driven parameterization of the cortical surface for localization and inter-subject matching.

In this paper we present a generic and organized model of cortical folding, and a way to implement this model on any given cortical surface. This results in a model-driven parameterization, providing an anatomically meaningful coordinate system for cortical localization, and implicitly defining inter-subject surface matching without any deformation of surfaces. We present our cortical folding model and show how it naturally defines a parameterization of the cortex. The mapping of the model to any given cortical surface is detailed, leading to an anatomically invariant coordinate system. The process is evaluated on real data in terms of both anatomical and functional localization, and shows improved performance compared to a traditional volume-based normalization. It is fully automatic and available with the BrainVISA software platform.

Genetics of primary cerebral gyrification: Heritability of length, depth and area of primary sulci in an extended pedigree of Papio baboons.

Genetic control over morphological variability of primary sulci and gyri is of great interest in the evolutionary, developmental and clinical neurosciences. Primary structures emerge early in development and their morphology is thought to be related to neuronal differentiation, development of functional connections and cortical lateralization. We measured the proportional contributions of genetics and environment to regional variability, testing two theories regarding regional modulation of genetic influences by ontogenic and phenotypic factors. Our measures were surface area, and average length and depth of eleven primary cortical sulci from high-resolution MR images in 180 pedigreed baboons. Average heritability values for sulcal area, depth and length (h(2)(Area)=.38+/-.22; h(2)(Depth)=.42+/-.23; h(2)(Length)=.34+/-.22) indicated that regional cortical anatomy is under genetic control. The regional pattern of genetic contributions was complex and, contrary to previously proposed theories, did not depend upon sulcal depth, or upon the sequence in which structures appear during development. Our results imply that heritability of sulcal phenotypes may be regionally modulated by arcuate U-fiber systems. However, further research is necessary to unravel the complexity of genetic contributions to cortical morphology.

Projection of fMRI data onto the cortical surface using anatomically-informed convolution kernels.

As surface-based data analysis offer an attractive approach for intersubject matching and comparison, the projection of voxel-based 3D volumes onto the cortical surface is an essential problem. We present here a method that aims at producing representations of functional brain data on the cortical surface from functional MRI volumes. Such representations are for instance required for subsequent cortical-based functional analysis. We propose a projection technique based on the definition, around each node of the gray/white matter interface mesh, of convolution kernels whose shape and distribution rely on the geometry of the local anatomy. For one anatomy, a set of convolution kernels is computed that can be used to project any functional data registered with this anatomy. Therefore resulting in anatomically-informed projections of data onto the cortical surface, this kernel-based approach offers better sensitivity, specificity than other classical methods and robustness to misregistration errors. Influences of mesh and volumes spatial resolutions were also estimated for various projection techniques, using simulated functional maps.

Anatomically constrained surface parameterization for cortical localization.

We present here a method that aims at defining a surface-based coordinate system on the cortical surface. Such a system is needed for both cortical localization and intersubject matching in the framework of neuroimaging. We propose an automatic parameterization based on the spherical topology of the grey/white matter interface of each hemisphere and on the use of naturally organized and reproducible anatomical features. From those markers used as initial constraints, the coordinate system is propagated via a PDE solved on the cortical surface.

Cervical spinal cord MTR histogram analysis in multiple sclerosis using a 3D acquisition and a B-spline active surface segmentation technique.

The application of a three-dimensional magnetization transfer (MT) sequence and B-spline active surface segmentation method to produce MT histograms of the cervical spinal cord in a pilot study of controls and multiple sclerosis (MS) patients is presented. Subjects' cervical spinal cords were imaged with (a) a volume-acquired inversion-prepared fast spoiled gradient echo sequence and (b) a volume-acquired noninversion-prepared fast spoiled gradient echo MT sequence. The images were segmented using the B spline active surface technique and MT histograms were produced from the MT images. The method was sensitive enough to detect differences between seven MS patients and 10 controls in mean MT ratio (42.4 pu versus 44.0 pu, p = 0.03) and peak location (45.2 versus 46.8, p = 0.03). The spinal cord volumes obtained from the two sequences were associated with each other (parameter estimate 0.972, 95% confidence intervals 0.742, 1.202, p < 0.001).

Coordinate-based versus structural approaches to brain image analysis.

A basic issue in neurosciences is to look for possible relationships between brain architecture and cognitive models. The lack of architectural information in magnetic resonance images, however, has led the neuroimaging community to develop brain mapping strategies based on various coordinate systems without accurate architectural content. Therefore, the relationships between architectural and functional brain organizations are difficult to study when analyzing neuroimaging experiments. This paper advocates that the design of new brain image analysis methods inspired by the structural strategies often used in computer vision may provide better ways to address these relationships. The key point underlying this new framework is the conversion of the raw images into structural representations before analysis. These representations are made up of data-driven elementary features like activated clusters, cortical folds or fiber bundles. Two classes of methods are introduced. Inference of structural models via matching across a set of individuals is described first. This inference problem is illustrated by the group analysis of functional statistical parametric maps (SPMs). Then, the matching of new individual data with a priori known structural models is described, using the recognition of the cortical sulci as a prototypical example.

Diffusion tensor magnetic resonance image regularization.

As multi-dimensional complex data become more common, new regularization schemes tailored to those data are needed. In this paper we present a scheme for regularising diffusion tensor magnetic resonance (DT-MR) data, and more generally multi-dimensional data defined by a direction map and one or several magnitude maps. The scheme is divided in two steps. First, a variational method is proposed to restore direction fields with preservation of discontinuities. Its theoretical aspects are presented, as well as its application to the direction field that defines the main orientation of the diffusion tensors. The second step makes use of an anisotropic diffusion process to regularize the magnitude maps. The main idea is that for a range of data it is possible to use the restored direction as a prior to drive the regularization process in a way that preserves discontinuities and respects the local coherence of the magnitude map. We show that anisotropic diffusion is a convenient framework to implement that idea, and define a regularization process for the magnitude maps from our DT-MR data. Both steps are illustated on synthetic and real diffusion tensor magnetic resonance data.

Quantification of spinal cord atrophy from magnetic resonance images via a B-spline active surface model.

A method is presented that aims at segmenting and measuring the surface of the spinal cord from MR images in order to detect and quantify atrophy. A semiautomatic segmentation with very little intervention from an operator is proposed. It is based on the optimization of a B-spline active surface. The method allows for the computation of orthogonal cross-sections at any level along the cord, from which measurements are derived, such as cross-sectional area or curvature. An evaluation of the accuracy and reproducibility of the method is presented.

Structural group analysis of functional activation maps.

We present here a new method for cerebral activation detection over a group of subjects. This method is performed using individual activation maps of any sort. It aims at processing a group analysis while preserving individual information and at overcoming as far as possible limitations of the spatial normalization used to compare different subjects. We designed it such that it provides the individual occurrence of the activations detected at a group level. The localization can then be performed on the individual anatomy of each subject. The analysis starts with a hierarchical multiscale object-based description of each individual map. These descriptions are then compared, rather than comparing the images directly. The analysis is thus performed at an object level instead of voxel by voxel. It is made using a comparison graph, on which a labeling process is performed. The label field on the graph is modeled by a Markov random field, which allows us to introduce high-level rules of interrogation of the data. The process has been evaluated on simulated data and real data from a PET protocol.