Detection of cortical transition regions utilizing statistical analyses of excess masses

Deep learning networks reflect cytoarchitectonic features used in brain mapping

The distribution of neurons in the cortex (cytoarchitecture) differs between cortical areas and constitutes the basis for structural maps of the human brain. Deep learning approaches provide a promising alternative to overcome throughput limitations of currently used cytoarchitectonic mapping methods, but typically lack insight as to what extent they follow cytoarchitectonic principles. We therefore investigated in how far the internal structure of deep convolutional neural networks trained for cytoarchitectonic brain mapping reflect traditional cytoarchitectonic features, and compared them to features of the current grey level index (GLI) profile approach. The networks consisted of a 10-block deep convolutional architecture trained to segment the primary and secondary visual cortex. Filter activations of the networks served to analyse resemblances to traditional cytoarchitectonic features and comparisons to the GLI profile approach. Our analysis revealed resemblances to cellular, laminar- as well as cortical area related cytoarchitectonic features. The networks learned filter activations that reflect the distinct cytoarchitecture of the segmented cortical areas with special regard to their laminar organization and compared well to statistical criteria of the GLI profile approach. These results confirm an incorporation of relevant cytoarchitectonic features in the deep convolutional neural networks and mark them as a valid support for high-throughput cytoarchitectonic mapping workflows.

Architectonic mapping of the medial region of the human orbitofrontal cortex by density profiles

The study of architectonic differentiation in the cortex is advanced by the articulation of objective definitions based on clear features of the cortical architecture. We adopted areal density profiles as a means of sampling the cortex. On the profiles, we isolated and quantified the density of individual cortical layers. These features may serve as criteria in objective definitions in a way that builds on qualitative observations found in the classical literature. A preprocessing procedure was introduced to overcome artefacts in the density profiles caused by the partial overlap of neighboring neuronal layers and cortical folding. We applied this method to the medial half of the orbital frontal cortex in specimens drawn from 10 human postmortem brain hemispheres. The measurements successfully confirmed the existence of several qualitatively observed areas (architectonic areas 14c, 14r, 11m, 11 and 13). The selection of specific sampling parameters was justified on the basis of simultaneous measurements of the cortical morphology which demonstrate its influence on the appearance of the cortical layers. We also examined the robustness of the measuring procedure by analyzing the outcome of varying systematically the sampling parameters. We describe here a novel method of sampling the cortex for architectonic analysis and demonstrate its application on histological sections obtained from the medial half of the human orbitofrontal cortex.

Quantitative architectural analysis: a new approach to cortical mapping

Recent progress in anatomical and functional MRI has revived the demand for a reliable, topographic map of the human cerebral cortex. Till date, interpretations of specific activations found in functional imaging studies and their topographical analysis in a spatial reference system are, often, still based on classical architectonic maps. The most commonly used reference atlas is that of Brodmann and his successors, despite its severe inherent drawbacks. One obvious weakness in traditional, architectural mapping is the subjective nature of localising borders between cortical areas, by means of a purely visual, microscopical examination of histological specimens. To overcome this limitation, more objective, quantitative mapping procedures have been established in the past years. The quantification of the neocortical, laminar pattern by defining intensity line profiles across the cortical layers, has a long tradition. During the last years, this method has been extended to enable a reliable, reproducible mapping of the cortex based on image analysis and multivariate statistics. Methodological approaches to such algorithm-based, cortical mapping were published for various architectural modalities. In our contribution, principles of algorithm-based mapping are described for cyto- and receptorarchitecture. In a cytoarchitectural parcellation of the human auditory cortex, using a sliding window procedure, the classical areal pattern of the human superior temporal gyrus was modified by a replacing of Brodmann's areas 41, 42, 22 and parts of area 21, with a novel, more detailed map. An extension and optimisation of the sliding window procedure to the specific requirements of receptorarchitectonic mapping, is also described using the macaque central sulcus and adjacent superior parietal lobule as a second, biologically independent example. Algorithm-based mapping procedures, however, are not limited to these two architectural modalities, but can be applied to all images in which a laminar cortical pattern can be detected and quantified, e.g. myeloarchitectonic and in vivo high resolution MR imaging. Defining cortical borders, based on changes in cortical lamination in high resolution, in vivo structural MR images will result in a rapid increase of our knowledge on the structural parcellation of the human cerebral cortex.

Three-dimensional cytoarchitectonic analysis of the posterior bank of the human precentral sulcus

Studies employing functional magnetic resonance imaging have identified the human frontal eye field as being in the anterior and partly in the posterior wall, as well as at the base of the precentral sulcus. Moreover, it is known that the frontal eye field extends rostrally to the superior frontal sulcus. According to Brodmann's cytoarchitectonic map, this region belongs to the dysgranular Brodmann area 6 of the premotor cortex. However, the frontal eye field in non-human primates has been located within the arcuate sulcus in Brodmann area 8, generating considerable debate about where to locate exactly the frontal eye field in humans. Functional studies of the primate frontal eye field have revealed a principal homology of voluntary saccadic control systems in human and old-world monkeys, especially the macaque. But these homologies seem to be contradicted by the reported topographic localization at the cytoarchitectonic level. Therefore, we studied the cytoarchitectonic structure of the posterior bank of the precentral sulcus of a human brain, employing newly developed spatial mapping techniques to provide data about whether or not this region should be considered cytoarchitecturally homogeneous or heterogeneous. We used functional magnetic resonance imaging results, as an initial guide in localizing a region which is activated by saccadic tasks. A maximum of activation was detected around the junction of the superior frontal sulcus and the precentral sulcus extending 1.5 cm along the precentral sulcus in direction of the lateral sulcus. Here, one human brain has been analyzed to obtain preliminary data about the cytoarchitectonical changes of a part of area 6. Statistical analysis of the three-dimensional architectonic data from this region allowed us to identify a zone at the posterior bank, which in other studies has been associated with a functional region that controls pursuit eye movements and performs sensory-to-motor transformations. We found two significant sectors along the ventral part of the posterior bank of the precentral sulcus. The caudal transition region coincides partly with a region that integrates retinal and eye position signals for target location, arm, and axial movements. The second more ventrally located region is attributed to process oral-facial movements. The caudal transition region coincides with our functional magnetic resonance imaging investigation. It was revealed that this region lies at the inferior frontal eye field, where a pronounced activation over a larger region can be stimulated. Currently, more studies are needed to combine functional magnetic resonance imaging data of maximal activation with data from whole histologic brain sections of more individuals and to quantify the variability of this region and its sub-regions by means of a standardized brain atlas.

Videomicroscopy, image processing, and analysis of whole histologic sections of the human brain

Serial histologic sections of a whole human brain may have extensions of up to 130 x 130 mm within the coronal plane around the temporal lobe. To date, however, technology has not provided a bright field microscope that is able to shift the object holder continuously in the x- and y-direction over such distances and still possess the same optical capabilities as comparable devices. We developed a new light microscope to continuously quantify such sections. We also developed the computing environment for controlling the device and for analyzing the data produced. In principle, we are now able to quantify each neuron of a human brain. The data ultimately will provide the most detailed structural information about the human brain ascertained thus far. Such detailed information of the spatial distribution of neurons is essential to develop realistic models for simulation of large-scale neuronal networks and to investigate the significance of neuronal arrangements with respect to neuronal signal processing in the CNS. After preprocessing of the data produced by the new microscope, we are able to detect lamination patterns in the spatial distribution of gravity centers of cells. Furthermore, morphological features like size of the projection area and mean staining intensity are visualized as a particle process. The particle process presents the sizes and staining intensity of perikaryons and allows a distinction of gray matter and white matter. These results provide evidence that the system works correctly and can be applied to a systematic analysis of a larger sequence of serial histologic sections. The objective of this study is to introduce the very large section analyzing microscope (VLSAM) and to present the initial data produced by the system. Moreover, we will discuss workload and future developments of the parallel image analysis system that are associated with the microscope.

Analysis of nerve fibers and their distribution in histologic sections of the human brain

The field of quantitative analysis and subsequent mapping of the cerebral cortex has developed rapidly. New powerful tools have been applied to investigate large regions of complex folded gyrencephalic cortices in order to detect structural transition regions that might partition different cortical fields of disjunct neuronal functions. We have developed a new mapping approach based on axoarchitectonics, a method of cortical visualization that previously has been used only indirectly with regard to myeloarchitectonics. Myeloarchitectonic visualization has the disadvantage of producing strong agglomerative effects of closely neighbored nerve fibers. Therefore, single and neurofunctional-relevant parameters such as axonal branchings, axon areas, and axon numbers have not been determinable with satisfying precision. As a result, different staining techniques had to be explored in order to achieve a suitable histologic staining for axon visualization. The best results were obtained after modifying the Naoumenko-Feigin staining for axons. From these contrast-rich stained histologic sections, videomicroscopic digital image tiles were generated and analyzed using a new fiber analysis framework. Finally, the analysis of histologic images provided topologic ordered parameters of axons that were transferred into parameter maps. The axon parameter maps were analyzed further via a recently developed traverse generating algorithm that calculated test lines oriented perpendicular to the cortical surface and white matter border. The gray value coded parameters of the parameter maps were then transferred into profile arrays. These profile arrays were statistically analyzed by a reliable excess mass approach we recently developed. We found that specific axonal parameters are preferentially distributed throughout granular and agranular types of cortex. Furthermore, our new procedure detected transition regions originally defined by changes of cytoarchitectonic layering. Statistically significant inhomogeneities of the distribution of certain axon quantities were shown to indicate a subparcellation of areas 4 and 6. The quantification techniques established here for the analysis of spatial axon distributions within larger regions of the cerebral cortex are suitable to detect inhomogeneities of laminar axon patterns. Hence, these techniques can be recommended for systematic and observer-supported cortical area mapping and parcellation studies.