BigBrain 3D atlas of cortical layers: Cortical and laminar thickness gradients diverge in sensory and motor cortices

Interpretable machine learning approach for neuron-centric analysis of human cortical cytoarchitecture

The complexity of the cerebral cortex underlies its function and distinguishes us as humans. Here, we present a principled veridical data science methodology for quantitative histology that shifts focus from image-level investigations towards neuron-level representations of cortical regions, with the neurons in the image as a subject of study, rather than pixel-wise image content. Our methodology relies on the automatic segmentation of neurons across whole histological sections and an extensive set of engineered features, which reflect the neuronal phenotype of individual neurons and the properties of neurons' neighborhoods. The neuron-level representations are used in an interpretable machine learning pipeline for mapping the phenotype to cortical layers. To validate our approach, we created a unique dataset of cortical layers manually annotated by three experts in neuroanatomy and histology. The presented methodology offers high interpretability of the results, providing a deeper understanding of human cortex organization, which may help formulate new scientific hypotheses, as well as to cope with systematic uncertainty in data and model predictions.

Layer-dependent stability of intracortical recordings and neuronal cell loss

Intracortical recordings can be used to voluntarily control external devices via brain-machine interfaces (BMI). Multiple factors, including the foreign body response (FBR), limit the stability of these neural signals over time. Current clinically approved devices consist of multi-electrode arrays with a single electrode site at the tip of each shank, confining the recording interface to a single layer of the cortex. Advancements in manufacturing technology have led to the development of high-density electrodes that can record from multiple layers. However, the long-term stability of neural recordings and the extent of neuronal cell loss around the electrode across different cortical depths have yet to be explored. To answer these questions, we recorded neural signals from rats chronically implanted with a silicon-substrate microelectrode array spanning the layers of the cortex. Our results show the long-term stability of intracortical recordings varies across cortical depth, with electrode sites around L4-L5 having the highest stability. Using machine learning guided segmentation, our novel histological technique, DeepHisto, revealed that the extent of neuronal cell loss varies across cortical layers, with L2/3 and L4 electrodes having the largest area of neuronal cell loss. These findings suggest that interfacing depth plays a major role in the FBR and long-term performance of intracortical neuroprostheses.

Three-Dimensional Digital Reconstruction of the Cerebellar Cortex: Lobule Thickness, Surface Area Measurements, and Layer Architecture

The cerebellum is ontogenetically one of the first structures to develop in the central nervous system; nevertheless, it has been only recently reconsidered for its significant neurobiological, functional, and clinical relevance in humans. Thus, it has been a relatively under-studied compared to the cerebrum. Currently, non-invasive imaging modalities can barely reach the necessary resolution to unfold its entire, convoluted surface, while only histological analyses can reveal local information at the micrometer scale. Herein, we used the BigBrain dataset to generate area and point-wise thickness measurements for all layers of the cerebellar cortex and for each lobule in particular. We found that the overall surface area of the cerebellar granular layer (including Purkinje cells) was 1,732 cm2 and the molecular layer was 1,945 cm2. The average thickness of the granular layer is 0.88 mm (± 0.83) and that of the molecular layer is 0.32 mm (± 0.08). The cerebellum (both granular and molecular layers) is thicker at the depth of the sulci and thinner at the crowns of the gyri. Globally, the granular layer is thicker in the lateral-posterior-inferior region than the medial-superior regions. The characterization of individual layers in the cerebellum achieved herein represents a stepping-stone for investigations interrelating structural and functional connectivity with cerebellar architectonics using neuroimaging, which is a matter of considerable relevance in basic and clinical neuroscience. Furthermore, these data provide templates for the construction of cerebellar topographic maps and the precise localization of structural and functional alterations in diseases affecting the cerebellum.

Characterization of regional differences in resting-state fMRI with a data-driven network model of brain dynamics

Model-based data analysis of whole-brain dynamics links the observed data to model parameters in a network of neural masses. Recently, studies focused on the role of regional variance of model parameters. Such analyses however necessarily depend on the properties of preselected neural mass model. We introduce a method to infer from the functional data both the neural mass model representing the regional dynamics and the region- and subject-specific parameters while respecting the known network structure. We apply the method to human resting-state fMRI. We find that the underlying dynamics can be described as noisy fluctuations around a single fixed point. The method reliably discovers three regional parameters with clear and distinct role in the dynamics, one of which is strongly correlated with the first principal component of the gene expression spatial map. The present approach opens a novel way to the analysis of resting-state fMRI with possible applications for understanding the brain dynamics during aging or neurodegeneration.

Mapping Cytoarchitectonics and Receptor Architectonics to Understand Brain Function and Connectivity

This review focuses on cytoarchitectonics and receptor architectonics as biological correlates of function and connectivity. It introduces the 3-dimensional cytoarchitectonic probabilistic maps of cortical areas and nuclei of the Julich-Brain Atlas, available at EBRAINS, to study structure-function relationships. The maps are linked to the BigBrain as microanatomical reference model and template space. The siibra software tool suite enables programmatic access to the maps and to receptor architectonic data that are anchored to brain areas. Such cellular and molecular data are tools for studying magnetic resonance connectivity including modeling and simulation. At the end, we highlight perspectives of the Julich-Brain as well as methodological considerations. Thus, microstructural maps as part of a multimodal atlas help elucidate the biological correlates of large-scale networks and brain function with a high level of anatomical detail, which provides a basis to study brains of patients with psychiatric disorders.

BrainStat: A toolbox for brain-wide statistics and multimodal feature associations

Analysis and interpretation of neuroimaging datasets has become a multidisciplinary endeavor, relying not only on statistical methods, but increasingly on associations with respect to other brain-derived features such as gene expression, histological data, and functional as well as cognitive architectures. Here, we introduce BrainStat - a toolbox for (i) univariate and multivariate linear models in volumetric and surface-based brain imaging datasets, and (ii) multidomain feature association of results with respect to spatial maps of post-mortem gene expression and histology, task-based fMRI meta-analysis, as well as resting-state fMRI motifs across several common surface templates. The combination of statistics and feature associations into a turnkey toolbox streamlines analytical processes and accelerates cross-modal research. The toolbox is implemented in both Python and MATLAB, two widely used programming languages in the neuroimaging and neuroinformatics communities. BrainStat is openly available and complemented by an expandable documentation.

Neurophysiological signatures of cortical micro-architecture

Systematic spatial variation in micro-architecture is observed across the cortex. These micro-architectural gradients are reflected in neural activity, which can be captured by neurophysiological time-series. How spontaneous neurophysiological dynamics are organized across the cortex and how they arise from heterogeneous cortical micro-architecture remains unknown. Here we extensively profile regional neurophysiological dynamics across the human brain by estimating over 6 800 timeseries features from the resting state magnetoencephalography (MEG) signal. We then map regional time-series profiles to a comprehensive multi-modal, multi-scale atlas of cortical micro-architecture, including microstructure, metabolism, neurotransmitter receptors, cell types and laminar differentiation. We find that the dominant axis of neurophysiological dynamics reflects characteristics of power spectrum density and linear correlation structure of the signal, emphasizing the importance of conventional features of electromagnetic dynamics while identifying additional informative features that have traditionally received less attention. Moreover, spatial variation in neurophysiological dynamics is colocalized with multiple micro-architectural features, including genomic gradients, intracortical myelin, neurotransmitter receptors and transporters, and oxygen and glucose metabolism. Collectively, this work opens new avenues for studying the anatomical basis of neural activity.

Gradual changes in microarchitectural properties of cortex and juxtacortical white matter: Observed by anatomical and diffusion MRI

PURPOSE

Characterization of cerebral cortex is challenged by the complexity and heterogeneity of its cyto- and myeloarchitecture. This study evaluates quantitative MRI metrics, measured across two cortical depths and in subcortical white matter (WM) adjacent to cortex (juxtacortical WM), indicative of myelin content, neurite density, and diffusion microenvironment, for a comprehensive characterization of cortical microarchitecture.

METHODS

High-quality structural and diffusion MRI data (N = 30) from the Human Connectome Project were processed to compute myelin index, neurite density index, fractional anisotropy, mean diffusivity, axial diffusivity, and radial diffusivity from superficial cortex, deep cortex, and juxtacortical WM. The distributional patterns of these metrics were analyzed individually, correlated to one another, and were compared to established parcellations.

RESULTS

Our results supported that myeloarchitectonic and the coexisting cytoarchitectonic structures influence the diffusion properties of water molecules residing in cortex. Full cortical thickness showed myelination patterns similar to those previously observed in humans. Higher myelin indices with similar distributional patterns were observed in deep cortex whereas lower myelin indices were observed in superficial cortex. Neurite density index and other diffusion MRI derived parameters provided complementary information to myelination. Reliable and reproducible correlations were identified among the cortical microarchitectural properties and fiber distributional patterns in proximal WM structures.

CONCLUSION

We demonstrated gradual changes across the cortical sheath by assessing depth-specific cortical micro-architecture using anatomical and diffusion MRI. Mutually independent but coexisting features of cortical layers and juxtacortical WM provided new insights towards structural organizational units and variabilities across cortical regions and through depth.

High-resolution cortical MAP-MRI reveals areal borders and laminar substructures observed with histological staining

The variations in cellular composition and tissue architecture measured with histology provide the biological basis for partitioning the brain into distinct cytoarchitectonic areas and for characterizing neuropathological tissue alterations. Clearly, there is an urgent need to develop whole-brain neuroradiological methods that can assess cortical cyto- and myeloarchitectonic features non-invasively. Mean apparent propagator (MAP) MRI is a clinically feasible diffusion MRI method that quantifies efficiently and comprehensively the net microscopic displacements of water molecules diffusing in tissues. We investigate the sensitivity of high-resolution MAP-MRI to detecting areal and laminar variations in cortical cytoarchitecture and compare our results with observations from corresponding histological sections in the entire brain of a rhesus macaque monkey. High-resolution images of MAP-derived parameters, in particular the propagator anisotropy (PA), non-gaussianity (NG), and the return-to-axis probability (RTAP) reveal cortical area-specific lamination patterns in good agreement with the corresponding histological stained sections. In a few regions, the MAP parameters provide superior contrast to the five histological stains used in this study, delineating more clearly boundaries and transition regions between cortical areas and laminar substructures. Throughout the cortex, various MAP parameters can be used to delineate transition regions between specific cortical areas observed with histology and to refine areal boundaries estimated using atlas registration-based cortical parcellation. Using surface-based analysis of MAP parameters we quantify the cortical depth dependence of diffusion propagators in multiple regions-of-interest in a consistent and rigorous manner that is largely independent of the cortical folding geometry. The ability to assess cortical cytoarchitectonic features efficiently and non-invasively, its clinical feasibility, and translatability make high-resolution MAP-MRI a promising 3D imaging tool for studying whole-brain cortical organization, characterizing abnormal cortical development, improving early diagnosis of neurodegenerative diseases, identifying targets for biopsies, and complementing neuropathological investigations.

Long-term learning transforms prefrontal cortex representations during working memory

The role of the lateral prefrontal cortex (lPFC) in working memory (WM) is debated. Non-human primate (NHP) electrophysiology shows that the lPFC stores WM representations, but human neuroimaging suggests that the lPFC controls WM content in sensory cortices. These accounts are confounded by differences in task training and stimulus exposure. We tested whether long-term training alters lPFC function by densely sampling WM activity using functional MRI. Over 3 months, participants trained on both a WM and serial reaction time (SRT) task, wherein fractal stimuli were embedded within sequences. WM performance improved for trained (but not novel) fractals and, neurally, delay activity increased in distributed lPFC voxels across learning. Item-level WM representations became detectable within lPFC patterns, and lPFC activity reflected sequence relationships from the SRT task. These findings demonstrate that human lPFC develops stimulus-selective responses with learning, and WM representations are shaped by long-term experience, which could reconcile competing accounts of WM functioning.

Evolution of cortical neurons supporting human cognition

Human cognitive abilities are generally thought to arise from cortical expansion over the course of human brain evolution. In addition to increased neuron numbers, this cortical expansion might be driven by adaptations in the properties of single neurons and their local circuits. We review recent findings on the distinct structural, functional, and transcriptomic features of human cortical neurons and their organization in cortical microstructure. We focus on the supragranular cortical layers, which showed the most prominent expansion during human brain evolution, and the properties of their principal cells: pyramidal neurons. We argue that the evolutionary adaptations in neuronal features that accompany the expansion of the human cortex partially underlie interindividual variability in human cognitive abilities.

Broad transcriptomic dysregulation occurs across the cerebral cortex in ASD

Neuropsychiatric disorders classically lack defining brain pathologies, but recent work has demonstrated dysregulation at the molecular level, characterized by transcriptomic and epigenetic alterations1-3. In autism spectrum disorder (ASD), this molecular pathology involves the upregulation of microglial, astrocyte and neural-immune genes, the downregulation of synaptic genes, and attenuation of gene-expression gradients in cortex1,2,4-6. However, whether these changes are limited to cortical association regions or are more widespread remains unknown. To address this issue, we performed RNA-sequencing analysis of 725 brain samples spanning 11 cortical areas from 112 post-mortem samples from individuals with ASD and neurotypical controls. We find widespread transcriptomic changes across the cortex in ASD, exhibiting an anterior-to-posterior gradient, with the greatest differences in primary visual cortex, coincident with an attenuation of the typical transcriptomic differences between cortical regions. Single-nucleus RNA-sequencing and methylation profiling demonstrate that this robust molecular signature reflects changes in cell-type-specific gene expression, particularly affecting excitatory neurons and glia. Both rare and common ASD-associated genetic variation converge within a downregulated co-expression module involving synaptic signalling, and common variation alone is enriched within a module of upregulated protein chaperone genes. These results highlight widespread molecular changes across the cerebral cortex in ASD, extending beyond association cortex to broadly involve primary sensory regions.

High-resolution atlasing and segmentation of the subcortex: Review and perspective on challenges and opportunities created by machine learning

This paper reviews almost three decades of work on atlasing and segmentation methods for subcortical structures in human brain MRI. In writing this survey, we have three distinct aims. First, to document the evolution of digital subcortical atlases of the human brain, from the early MRI templates published in the nineties, to the complex multi-modal atlases at the subregion level that are available today. Second, to provide a detailed record of related efforts in the automated segmentation front, from earlier atlas-based methods to modern machine learning approaches. And third, to present a perspective on the future of high-resolution atlasing and segmentation of subcortical structures in in vivo human brain MRI, including open challenges and opportunities created by recent developments in machine learning.

Mapping neurotransmitter systems to the structural and functional organization of the human neocortex

Neurotransmitter receptors support the propagation of signals in the human brain. How receptor systems are situated within macro-scale neuroanatomy and how they shape emergent function remain poorly understood, and there exists no comprehensive atlas of receptors. Here we collate positron emission tomography data from more than 1,200 healthy individuals to construct a whole-brain three-dimensional normative atlas of 19 receptors and transporters across nine different neurotransmitter systems. We found that receptor profiles align with structural connectivity and mediate function, including neurophysiological oscillatory dynamics and resting-state hemodynamic functional connectivity. Using the Neurosynth cognitive atlas, we uncovered a topographic gradient of overlapping receptor distributions that separates extrinsic and intrinsic psychological processes. Finally, we found both expected and novel associations between receptor distributions and cortical abnormality patterns across 13 disorders. We replicated all findings in an independently collected autoradiography dataset. This work demonstrates how chemoarchitecture shapes brain structure and function, providing a new direction for studying multi-scale brain organization.

neuromaps: structural and functional interpretation of brain maps

Imaging technologies are increasingly used to generate high-resolution reference maps of brain structure and function. Comparing experimentally generated maps to these reference maps facilitates cross-disciplinary scientific discovery. Although recent data sharing initiatives increase the accessibility of brain maps, data are often shared in disparate coordinate systems, precluding systematic and accurate comparisons. Here we introduce neuromaps, a toolbox for accessing, transforming and analyzing structural and functional brain annotations. We implement functionalities for generating high-quality transformations between four standard coordinate systems. The toolbox includes curated reference maps and biological ontologies of the human brain, such as molecular, microstructural, electrophysiological, developmental and functional ontologies. Robust quantitative assessment of map-to-map similarity is enabled via a suite of spatial autocorrelation-preserving null models. neuromaps combines open-access data with transparent functionality for standardizing and comparing brain maps, providing a systematic workflow for comprehensive structural and functional annotation enrichment analysis of the human brain.

Inconsistencies in atlas-based volumetric measures of the human nucleus basalis of Meynert: A need for high-resolution alternatives

The nucleus basalis of Meynert (nbM) is the major source of cortical acetylcholine (ACh) and has been related to cognitive processes and to neurological disorders. However, spatially delineating the human nbM in MRI studies remains challenging. Due to the absence of a functional localiser for the human nbM, studies to date have localised it using nearby neuroanatomical landmarks or using probabilistic atlases. To understand the feasibility of MRI of the nbM we set our four goals; our first goal was to review current human nbM region-of-interest (ROI) selection protocols used in MRI studies, which we found have reported highly variable nbM volume estimates. Our next goal was to quantify and discuss the limitations of existing atlas-based volumetry of nbM. We found that the identified ROI volume depends heavily on the atlas used and on the probabilistic threshold set. In addition, we found large disparities even for data/studies using the same atlas and threshold. To test whether spatial resolution contributes to volume variability, as our third goal, we developed a novel nbM mask based on the normalized BigBrain dataset. We found that as long as the spatial resolution of the target data was 1.3 mm isotropic or above, our novel nbM mask offered realistic and stable volume estimates. Finally, as our last goal we tried to discern nbM using publicly available and novel high resolution structural MRI ex vivo MRI datasets. We find that, using an optimised 9.4T quantitative T₂^⁎ ex vivo dataset, the nbM can be visualised using MRI. We conclude caution is needed when applying the current methods of mapping nbM, especially for high resolution MRI data. Direct imaging of the nbM appears feasible and would eliminate the problems we identify, although further development is required to allow such imaging using standard (f)MRI scanning.

The cortical spectrum: A robust structural continuum in primate cerebral cortex revealed by histological staining and magnetic resonance imaging

High-level characterizations of the primate cerebral cortex sit between two extremes: on one end the cortical mantle is seen as a mosaic of structurally and functionally unique areas, and on the other it is seen as a uniform six-layered structure in which functional differences are defined solely by extrinsic connections. Neither of these extremes captures the crucial neuroanatomical finding: that the cortex exhibits systematic gradations in architectonic structure. These gradations have been shown to predict cortico-cortical connectivity, which in turn suggests powerful ways to ground connectomics in anatomical structure, and by extension cortical function. A challenge to widespread use of this concept is the labor-intensive and invasive nature of histological staining, which is the primary means of recognizing anatomical gradations. Here we show that a novel computational analysis technique can provide a coarse-grained picture of cortical variation. For each of 78 cortical areas spanning the entire cortical mantle of the rhesus macaque, we created a high dimensional set of anatomical features derived from captured images of cortical tissue stained for myelin and SMI-32. The method involved semi-automated de-noising of images, and enabled comparison of brain areas without hand-labeling of features such as layer boundaries. We applied multidimensional scaling (MDS) to the dataset to visualize similarity among cortical areas. This analysis shows a systematic variation between weakly laminated (limbic) cortices and sharply laminated (eulaminate) cortices. We call this smooth continuum the “cortical spectrum”. We also show that this spectrum is visible within subsystems of the cortex: the occipital, parietal, temporal, motor, prefrontal, and insular cortices. We compared the MDS-derived spectrum with a spectrum produced using T1- and T2-weighted magnetic resonance imaging (MRI) data derived from macaque, and found close agreement of the two coarse-graining methods. This suggests that T1w/T2w data, routinely obtained in human MRI studies, can serve as an effective proxy for data derived from high-resolution histological methods. More generally, this approach shows that the cortical spectrum is robust to the specific method used to compare cortical areas, and is therefore a powerful tool to understand the principles of organization of the primate cortex.

A whole-brain 3D myeloarchitectonic atlas: Mapping the Vogt-Vogt legacy to the cortical surface

Building precise and detailed parcellations of anatomically and functionally distinct brain areas has been a major focus in Neuroscience. Pioneer anatomists parcellated the cortical manifold based on extensive histological studies of post-mortem brain, harnessing local variations in cortical cyto- and myeloarchitecture to define areal boundaries. Compared to the cytoarchitectonic field, where multiple neuroimaging studies have recently translated this old legacy data into useful analytical resources, myeloarchitectonics, which parcellate the cortex based on the organization of myelinated fibers, has received less attention. Here, we present the neocortical surface-based myeloarchitectonic atlas based on the histology-derived maps of the Vogt-Vogt school and its 2D translation by Nieuwenhuys. In addition to a myeloarchitectonic parcellation, our package includes intracortical laminar profiles of myelin content based on Vogt-Vogt-Hopf original publications. Histology-derived myelin density mapped on our atlas demonstrated a close overlap with in vivo quantitative MRI markers for myelin and relates to cytoarchitectural features. Complementing the existing battery of approaches for digital cartography, the whole-brain myeloarchitectonic atlas offers an opportunity to validate imaging surrogate markers of myelin in both health and disease.

Human electromagnetic and haemodynamic networks systematically converge in unimodal cortex and diverge in transmodal cortex

Whole-brain neural communication is typically estimated from statistical associations among electromagnetic or haemodynamic time-series. The relationship between functional network architectures recovered from these 2 types of neural activity remains unknown. Here, we map electromagnetic networks (measured using magnetoencephalography (MEG)) to haemodynamic networks (measured using functional magnetic resonance imaging (fMRI)). We find that the relationship between the 2 modalities is regionally heterogeneous and systematically follows the cortical hierarchy, with close correspondence in unimodal cortex and poor correspondence in transmodal cortex. Comparison with the BigBrain histological atlas reveals that electromagnetic-haemodynamic coupling is driven by laminar differentiation and neuron density, suggesting that the mapping between the 2 modalities can be explained by cytoarchitectural variation. Importantly, haemodynamic connectivity cannot be explained by electromagnetic activity in a single frequency band, but rather arises from the mixing of multiple neurophysiological rhythms. Correspondence between the two is largely driven by MEG functional connectivity at the beta (15 to 29 Hz) frequency band. Collectively, these findings demonstrate highly organized but only partly overlapping patterns of connectivity in MEG and fMRI functional networks, opening fundamentally new avenues for studying the relationship between cortical microarchitecture and multimodal connectivity patterns.

Cerebral cortex layer segmentation using diffusion magnetic resonance imaging in vivo with applications to laminar connections and working memory analysis

Understanding the laminar brain structure is of great help in further developing our knowledge of the functions of the brain. However, since most layer segmentation methods are invasive, it is difficult to apply them to the human brain in vivo. To systematically explore the human brain's laminar structure noninvasively, the K-means clustering algorithm was used to automatically segment the left hemisphere into two layers, the superficial and deep layers, using a 7 Tesla (T) diffusion magnetic resonance imaging (dMRI)open dataset. The obtained layer thickness was then compared with the layer thickness of the BigBrain reference dataset, which segmented the neocortex into six layers based on the von Economo atlas. The results show a significant correlation not only between our automatically segmented superficial layer thickness and the thickness of layers 1-3 from the reference histological data, but also between our automatically segmented deep layer thickness and the thickness of layers 4-6 from the reference histological data. Second, we constructed the laminar connections between two pairs of unidirectional connected regions, which is consistent with prior research. Finally, we conducted the laminar analysis of the working memory, which was challenging to do in the past, and explained the conclusions of the functional analysis. Our work successfully demonstrates that it is possible to segment the human cortex noninvasively into layers using dMRI data and further explores the mechanisms of the human brain.

From Analog to Digital Computing: Is Homo sapiens’ Brain on Its Way to Become a Turing Machine?

The abstract basis of modern computation is the formal description of a finite state machine, the Universal Turing Machine, based on manipulation of integers and logic symbols. In this contribution to the discourse on the computer-brain analogy, we discuss the extent to which analog computing, as performed by the mammalian brain, is like and unlike the digital computing of Universal Turing Machines. We begin with ordinary reality being a permanent dialog between continuous and discontinuous worlds. So it is with computing, which can be analog or digital, and is often mixed. The theory behind computers is essentially digital, but efficient simulations of phenomena can be performed by analog devices; indeed, any physical calculation requires implementation in the physical world and is therefore analog to some extent, despite being based on abstract logic and arithmetic. The mammalian brain, comprised of neuronal networks, functions as an analog device and has given rise to artificial neural networks that are implemented as digital algorithms but function as analog models would. Analog constructs compute with the implementation of a variety of feedback and feedforward loops. In contrast, digital algorithms allow the implementation of recursive processes that enable them to generate unparalleled emergent properties. We briefly illustrate how the cortical organization of neurons can integrate signals and make predictions analogically. While we conclude that brains are not digital computers, we speculate on the recent implementation of human writing in the brain as a possible digital path that slowly evolves the brain into a genuine (slow) Turing machine.

GluN3 Subunit Expression Correlates with Increased Vulnerability of Hippocampus and Entorhinal Cortex to Neurodegeneration in a Model of Temporal Lobe Epilepsy

Temporal lobe epilepsy (TLE) is the most common type of epilepsy in adults that is often refractory to anti-epileptic medication therapy. Neither the pathology nor the etiology of TLE are fully characterized, although recent studies have established that the two are causally related. TLE pathology entails a stereotypic pattern of neuron loss in hippocampal and parahippocampal regions, predominantly in CA1 subfield of the hippocampus and layer 3 of the medial entorhinal area (MEA), deemed hallmark pathological features of the disease. Through this work, we address the contribution of glutamatergic N-methyl-D-aspartate receptors (NMDARs) to the pathology (vulnerability and pattern of neuronal loss), and by extension to the pathophysiology (Ca²⁺ induced excitotoxicity), by assaying the spatial expression of their subunit proteins (GluN1, GluN2A, GluN2B and GluN3A) in these regions using ASTA (area specific tissue analysis), a novel methodology for harvesting brain chads from hard-to-reach regions within brain slices for Western blotting. Our data suggest gradient expression of the GluN3A subunit along the mid-lateral extent of layer 3 MEA and along the CA1-subicular axis in the hippocampus, unlike GluN1 or GluN2 subunits which are uniformly distributed. Incorporation of GluN3A in the subunit composition of conventional diheteromeric (d-) NMDARs yield triheteromeric (t-) NMDARs which by virtue of their increased selectivity for Ca²⁺ render neurons vulnerable to excitotoxic damage. Thus, the expression profile of this subunit sheds light on the spatial extent of the pathology observed in these regions and implicates the GluN3 subunit of NMDARs in hippocampal and entorhinal cortical pathology underlying TLE.

Subtly altered topological asymmetry of brain structural covariance networks in autism spectrum disorder across 43 datasets from the ENIGMA consortium

Small average differences in the left-right asymmetry of cerebral cortical thickness have been reported in individuals with autism spectrum disorder (ASD) compared to typically developing controls, affecting widespread cortical regions. The possible impacts of these regional alterations in terms of structural network effects have not previously been characterized. Inter-regional morphological covariance analysis can capture network connectivity between different cortical areas at the macroscale level. Here, we used cortical thickness data from 1455 individuals with ASD and 1560 controls, across 43 independent datasets of the ENIGMA consortium's ASD Working Group, to assess hemispheric asymmetries of intra-individual structural covariance networks, using graph theory-based topological metrics. Compared with typical features of small-world architecture in controls, the ASD sample showed significantly altered average asymmetry of networks involving the fusiform, rostral middle frontal, and medial orbitofrontal cortex, involving higher randomization of the corresponding right-hemispheric networks in ASD. A network involving the superior frontal cortex showed decreased right-hemisphere randomization. Based on comparisons with meta-analyzed functional neuroimaging data, the altered connectivity asymmetry particularly affected networks that subserve executive functions, language-related and sensorimotor processes. These findings provide a network-level characterization of altered left-right brain asymmetry in ASD, based on a large combined sample. Altered asymmetrical brain development in ASD may be partly propagated among spatially distant regions through structural connectivity.

Cytoarchitectonic Maps of the Human Metathalamus in 3D Space

The human metathalamus plays an important role in processing visual and auditory information. Understanding its layers and subdivisions is important to gain insights in its function as a subcortical relay station and involvement in various pathologies. Yet, detailed histological references of the microanatomy in 3D space are still missing. We therefore aim at providing cytoarchitectonic maps of the medial geniculate body (MGB) and its subdivisions in the BigBrain – a high-resolution 3D-reconstructed histological model of the human brain, as well as probabilistic cytoarchitectonic maps of the MGB and lateral geniculate body (LGB). Therefore, histological sections of ten postmortem brains were studied. Three MGB subdivisions (MGBv, MGBd, MGBm) were identified on every 5th BigBrain section, and a deep-learning based tool was applied to map them on every remaining section. The maps were 3D-reconstructed to show the shape and extent of the MGB and its subdivisions with cellular precision. The LGB and MGB were additionally identified in nine other postmortem brains. Probabilistic cytoarchitectonic maps in the MNI “Colin27” and MNI ICBM152 reference spaces were computed which reveal an overall low interindividual variability in topography and extent. The probabilistic maps were included into the Julich-Brain atlas, and are freely available. They can be linked to other 3D data of human brain organization and serve as an anatomical reference for diagnostic, prognostic and therapeutic neuroimaging studies of healthy brains and patients. Furthermore, the high-resolution MGB BigBrain maps provide a basis for data integration, brain modeling and simulation to bridge the larger scale involvement of thalamocortical and local subcortical circuits.

Linking Brain Structure, Activity, and Cognitive Function through Computation

Understanding the human brain is a "Grand Challenge" for 21st century research. Computational approaches enable large and complex datasets to be addressed efficiently, supported by artificial neural networks, modeling and simulation. Dynamic generative multiscale models, which enable the investigation of causation across scales and are guided by principles and theories of brain function, are instrumental for linking brain structure and function. An example of a resource enabling such an integrated approach to neuroscientific discovery is the BigBrain, which spatially anchors tissue models and data across different scales and ensures that multiscale models are supported by the data, making the bridge to both basic neuroscience and medicine. Research at the intersection of neuroscience, computing and robotics has the potential to advance neuro-inspired technologies by taking advantage of a growing body of insights into perception, plasticity and learning. To render data, tools and methods, theories, basic principles and concepts interoperable, the Human Brain Project (HBP) has launched EBRAINS, a digital neuroscience research infrastructure, which brings together a transdisciplinary community of researchers united by the quest to understand the brain, with fascinating insights and perspectives for societal benefits.

Diffusion spectrum imaging predicts hippocampal sclerosis in mesial temporal lobe epilepsy patients

Objectives

Epileptic patients suffer from seizure recurrence after surgery due to the challenging localization. Improvement of the noninvasive imaging‐based approach for a better definition of the abnormalities would be helpful for a better outcome.

Methods

The quantitative anisotropy (QA) of diffusion spectrum imaging (DSI) is a quantitative scalar of evaluating the water diffusivity. Herein, we investigated the association between neuronal diameters or density acquired in literature and QA of DSI as well as the seizure localization in temporal lobe epilepsy. Thirty healthy controls (HCs) and 30 patients with hippocampal sclerosis (HS) were retrospectively analyzed. QA values were calculated and interactively compared between the areas with different neuronal diameter/density acquired from literature in the HCP‐1021 template. Diagnostic tests were performed on Z‐transformed asymmetry indices (AIs) of QA (which exclude physical asymmetry) among HS patients to evaluate its clinical value.

Results

The QA values in HCs conformed with different pyramidal cell distributions ranged from giant to small; corresponding groups were the motor‐sensory, associative, and limbic groups, respectively. Additionally, the QA value was correlated with the neuronal diameter/density in cortical layer IIIc (correlation coefficient with diameter: 0.529, p = 0.035; density: −0.678, p = 0.011). Decreases in cingulum hippocampal segments (Chs) were consistently observed on the sclerosed side in patients. The area under the curve of the Z‐transformed AI in Chs to the lateralization of HS was 0.957 (sensitivity: 0.909, specificity: 0.895).

Interpretation

QA based on DSI is likely to be useful to provide information to reflect the neuronal diameter/density and further facilitate localization of epileptic tissues.

Atlas of lesion locations and postsurgical seizure freedom in focal cortical dysplasia: A MELD study

OBJECTIVE

Drug-resistant focal epilepsy is often caused by focal cortical dysplasias (FCDs). The distribution of these lesions across the cerebral cortex and the impact of lesion location on clinical presentation and surgical outcome are largely unknown. We created a neuroimaging cohort of patients with individually mapped FCDs to determine factors associated with lesion location and predictors of postsurgical outcome.

METHODS

The MELD (Multi-centre Epilepsy Lesion Detection) project collated a retrospective cohort of 580 patients with epilepsy attributed to FCD from 20 epilepsy centers worldwide. Magnetic resonance imaging-based maps of individual FCDs with accompanying demographic, clinical, and surgical information were collected. We mapped the distribution of FCDs, examined for associations between clinical factors and lesion location, and developed a predictive model of postsurgical seizure freedom.

RESULTS

FCDs were nonuniformly distributed, concentrating in the superior frontal sulcus, frontal pole, and temporal pole. Epilepsy onset was typically before the age of 10 years. Earlier epilepsy onset was associated with lesions in primary sensory areas, whereas later epilepsy onset was associated with lesions in association cortices. Lesions in temporal and occipital lobes tended to be larger than frontal lobe lesions. Seizure freedom rates varied with FCD location, from around 30% in visual, motor, and premotor areas to 75% in superior temporal and frontal gyri. The predictive model of postsurgical seizure freedom had a positive predictive value of 70% and negative predictive value of 61%.

SIGNIFICANCE

FCD location is an important determinant of its size, the age at epilepsy onset, and the likelihood of seizure freedom postsurgery. Our atlas of lesion locations can be used to guide the radiological search for subtle lesions in individual patients. Our atlas of regional seizure freedom rates and associated predictive model can be used to estimate individual likelihoods of postsurgical seizure freedom. Data-driven atlases and predictive models are essential for evidence-based, precision medicine and risk counseling in epilepsy.

Learning Cortical Parcellations Using Graph Neural Networks

Deep learning has been applied to magnetic resonance imaging (MRI) for a variety of purposes, ranging from the acceleration of image acquisition and image denoising to tissue segmentation and disease diagnosis. Convolutional neural networks have been particularly useful for analyzing MRI data due to the regularly sampled spatial and temporal nature of the data. However, advances in the field of brain imaging have led to network- and surface-based analyses that are often better represented in the graph domain. In this analysis, we propose a general purpose cortical segmentation method that, given resting-state connectivity features readily computed during conventional MRI pre-processing and a set of corresponding training labels, can generate cortical parcellations for new MRI data. We applied recent advances in the field of graph neural networks to the problem of cortical surface segmentation, using resting-state connectivity to learn discrete maps of the human neocortex. We found that graph neural networks accurately learn low-dimensional representations of functional brain connectivity that can be naturally extended to map the cortices of new datasets. After optimizing over algorithm type, network architecture, and training features, our approach yielded mean classification accuracies of 79.91% relative to a previously published parcellation. We describe how some hyperparameter choices including training and testing data duration, network architecture, and algorithm choice affect model performance.

Computational models of cortical folding: A review of common approaches

The process of gyrification, by which the brain develops the intricate pattern of gyral hills and sulcal valleys, is the result of interactions between biological and mechanical processes during brain development. Researchers have developed a vast array of computational models in order to investigate cortical folding. This review aims to summarize these studies, focusing on five essential elements of the brain that affect development and gyrification and how they are represented in computational models: (i) the constraints of skull, meninges, and cerebrospinal fluid; (ii) heterogeneity of cortical layers and regions; (iii) anisotropic behavior of subcortical fiber tracts; (iv) material properties of brain tissue; and (v) the complex geometry of the brain. Finally, we highlight areas of need for future simulations of brain development.

The glial framework reveals white matter fiber architecture in human and primate brains

[Figure: see text].

Simulating Local Deformations in the Human Cortex Due to Blood Flow-Induced Changes in Mechanical Tissue Properties: Impact on Functional Magnetic Resonance Imaging

Investigating human brain tissue is challenging due to the complexity and the manifold interactions between structures across different scales. Increasing evidence suggests that brain function and microstructural features including biomechanical features are related. More importantly, the relationship between tissue mechanics and its influence on brain imaging results remains poorly understood. As an important example, the study of the brain tissue response to blood flow could have important theoretical and experimental consequences for functional magnetic resonance imaging (fMRI) at high spatial resolutions. Computational simulations, using realistic mechanical models can predict and characterize the brain tissue behavior and give us insights into the consequent potential biases or limitations of in vivo, high-resolution fMRI. In this manuscript, we used a two dimensional biomechanical simulation of an exemplary human gyrus to investigate the relationship between mechanical tissue properties and the respective changes induced by focal blood flow changes. The model is based on the changes in the brain’s stiffness and volume due to the vasodilation evoked by neural activity. Modeling an exemplary gyrus from a brain atlas we assessed the influence of different potential mechanisms: (i) a local increase in tissue stiffness (at the level of a single anatomical layer), (ii) an increase in local volume, and (iii) a combination of both effects. Our simulation results showed considerable tissue displacement because of these temporary changes in mechanical properties. We found that the local volume increase causes more deformation and consequently higher displacement of the gyrus. These displacements introduced considerable artifacts in our simulated fMRI measurements. Our results underline the necessity to consider and characterize the tissue displacement which could be responsible for fMRI artifacts.

Neurodevelopment of the association cortices: Patterns, mechanisms, and implications for psychopathology

The human brain undergoes a prolonged period of cortical development that spans multiple decades. During childhood and adolescence, cortical development progresses from lower-order, primary and unimodal cortices with sensory and motor functions to higher-order, transmodal association cortices subserving executive, socioemotional, and mentalizing functions. The spatiotemporal patterning of cortical maturation thus proceeds in a hierarchical manner, conforming to an evolutionarily rooted, sensorimotor-to-association axis of cortical organization. This developmental program has been characterized by data derived from multimodal human neuroimaging and is linked to the hierarchical unfolding of plasticity-related neurobiological events. Critically, this developmental program serves to enhance feature variation between lower-order and higher-order regions, thus endowing the brain's association cortices with unique functional properties. However, accumulating evidence suggests that protracted plasticity within late-maturing association cortices, which represents a defining feature of the human developmental program, also confers risk for diverse developmental psychopathologies.

The BigBrainWarp toolbox for integration of BigBrain 3D histology with multimodal neuroimaging

Neuroimaging stands to benefit from emerging ultrahigh-resolution 3D histological atlases of the human brain; the first of which is 'BigBrain'. Here, we review recent methodological advances for the integration of BigBrain with multi-modal neuroimaging and introduce a toolbox, 'BigBrainWarp', that combines these developments. The aim of BigBrainWarp is to simplify workflows and support the adoption of best practices. This is accomplished with a simple wrapper function that allows users to easily map data between BigBrain and standard MRI spaces. The function automatically pulls specialised transformation procedures, based on ongoing research from a wide collaborative network of researchers. Additionally, the toolbox improves accessibility of histological information through dissemination of ready-to-use cytoarchitectural features. Finally, we demonstrate the utility of BigBrainWarp with three tutorials and discuss the potential of the toolbox to support multi-scale investigations of brain organisation.

Cortical patterning of morphometric similarity gradient reveals diverged hierarchical organization in sensory-motor cortices

The topological organization of the cerebral cortex provides hierarchical axes, namely gradients, which reveal systematic variations of brain structure and function. However, the hierarchical organization of macroscopic brain morphology and how it constrains cortical function along the organizing axes remains unclear. We map the gradient of cortical morphometric similarity (MS) connectome, combining multiple features conceptualized as a "fingerprint" of an individual's brain. The principal MS gradient is anchored by motor and sensory cortices at two extreme ends, which are reliable and reproducible. Notably, divergences between motor and sensory hierarchies are consistent with the laminar histological thickness gradient but contrary to the canonical functional connectivity (FC) gradient. Moreover, the MS dissociates with FC gradients in the higher-order association cortices. The MS gradient recapitulates fundamental properties of cortical organization, from gene expression and cyto- and myeloarchitecture to evolutionary expansion. Collectively, our findings provide a heuristic hierarchical organization of cortical morphological neuromarkers.

Profiling Transcription Initiation in Peripheral Leukocytes Reveals Severity-Associated Cis-Regulatory Elements in Critical COVID-19

The contribution of transcription factors (TFs) and gene regulatory programs in the immune response to COVID-19 and their relationship to disease outcome is not fully understood. Analysis of genome-wide changes in transcription at both promoter-proximal and distal cis-regulatory DNA elements, collectively termed the 'active cistrome,' offers an unbiased assessment of TF activity identifying key pathways regulated in homeostasis or disease. Here, we profiled the active cistrome from peripheral leukocytes of critically ill COVID-19 patients to identify major regulatory programs and their dynamics during SARS-CoV-2 associated acute respiratory distress syndrome (ARDS). We identified TF motifs that track the severity of COVID- 19 lung injury, disease resolution, and outcome. We used unbiased clustering to reveal distinct cistrome subsets delineating the regulation of pathways, cell types, and the combinatorial activity of TFs. We found critical roles for regulatory networks driven by stimulus and lineage determining TFs, showing that STAT and E2F/MYB regulatory programs targeting myeloid cells are activated in patients with poor disease outcomes and associated with single nucleotide genetic variants implicated in COVID-19 susceptibility. Integration with single-cell RNA-seq found that STAT and E2F/MYB activation converged in specific neutrophils subset found in patients with severe disease. Collectively we demonstrate that cistrome analysis facilitates insight into disease mechanisms and provides an unbiased approach to evaluate global changes in transcription factor activity and stratify patient disease severity.

LayNii: A software suite for layer-fMRI

High-resolution fMRI in the sub-millimeter regime allows researchers to resolve brain activity across cortical layers and columns non-invasively. While these high-resolution data make it possible to address novel questions of directional information flow within and across brain circuits, the corresponding data analyses are challenged by MRI artifacts, including image blurring, image distortions, low SNR, and restricted coverage. These challenges often result in insufficient spatial accuracy of conventional analysis pipelines. Here we introduce a new software suite that is specifically designed for layer-specific functional MRI: LayNii. This toolbox is a collection of command-line executable programs written in C/C++ and is distributed opensource and as pre-compiled binaries for Linux, Windows, and macOS. LayNii is designed for layer-fMRI data that suffer from SNR and coverage constraints and thus cannot be straightforwardly analyzed in alternative software packages. Some of the most popular programs of LayNii contain 'layerification' and columnarization in the native voxel space of functional data as well as many other layer-fMRI specific analysis tasks: layer-specific smoothing, model-based vein mitigation of GE-BOLD data, quality assessment of artifact dominated sub-millimeter fMRI, as well as analyses of VASO data.

Linking cortical circuit models to human cognition with laminar fMRI

Laboratory animal research has provided significant knowledge into the function of cortical circuits at the laminar level, which has yet to be fully leveraged towards insights about human brain function on a similar spatiotemporal scale. The use of functional magnetic resonance imaging (fMRI) in conjunction with neural models provides new opportunities to gain important insights from current knowledge. During the last five years, human studies have demonstrated the value of high-resolution fMRI to study laminar-specific activity in the human brain. This is mostly performed at ultra-high-field strengths (≥ 7 T) and is known as laminar fMRI. Advancements in laminar fMRI are beginning to open new possibilities for studying questions in basic cognitive neuroscience. In this paper, we first review recent methodological advances in laminar fMRI and describe recent human laminar fMRI studies. Then, we discuss how the use of laminar fMRI can help bridge the gap between cortical circuit models and human cognition.

Convolutional neural networks for cytoarchitectonic brain mapping at large scale

Human brain atlases provide spatial reference systems for data characterizing brain organization at different levels, coming from different brains. Cytoarchitecture is a basic principle of the microstructural organization of the brain, as regional differences in the arrangement and composition of neuronal cells are indicators of changes in connectivity and function. Automated scanning procedures and observer-independent methods are prerequisites to reliably identify cytoarchitectonic areas, and to achieve reproducible models of brain segregation. Time becomes a key factor when moving from the analysis of single regions of interest towards high-throughput scanning of large series of whole-brain sections. Here we present a new workflow for mapping cytoarchitectonic areas in large series of cell-body stained histological sections of human postmortem brains. It is based on a Deep Convolutional Neural Network (CNN), which is trained on a pair of section images with annotations, with a large number of un-annotated sections in between. The model learns to create all missing annotations in between with high accuracy, and faster than our previous workflow based on observer-independent mapping. The new workflow does not require preceding 3D-reconstruction of sections, and is robust against histological artefacts. It processes large data sets with sizes in the order of multiple Terabytes efficiently. The workflow was integrated into a web interface, to allow access without expertise in deep learning and batch computing. Applying deep neural networks for cytoarchitectonic mapping opens new perspectives to enable high-resolution models of brain areas, introducing CNNs to identify borders of brain areas.

EEG Phase-Amplitude Coupling Strength and Phase Preference: Association with Age over the First Three Years after Birth

Phase-amplitude coupling (PAC), the coupling of the phase of slower electrophysiological oscillations with the amplitude of faster oscillations, is thought to facilitate dynamic integration of neural activity in the brain. Although the brain undergoes dramatic change and development during the first few years of life, how PAC changes through this developmental period has not been extensively studied. Here, we examined PAC through electroencephalography (EEG) data collected during an awake, eyes-open EEG collection paradigm in 98 children between the ages of three months and three years. We employed non-parametric clustering methods to identify areas of significant PAC across a range of frequency pairs and electrode locations, and examined how PAC strength and phase preference develops in these areas. We found that PAC, primarily between the α-β and γ frequencies, was positively correlated with age from early infancy to early childhood (p = 2.035 × 10-6). Additionally, we found γ over anterior electrodes coupled with the rising phase of the α-β waveform, while γ over posterior electrodes coupled with the falling phase of the α-β waveform; this regionalized phase preference became more prominent with age. This opposing trend may reflect each region's specialization toward feedback or feedforward processing, respectively, suggesting opportunities for back translation in future studies.

Layer-specific parameters of intracortical microstimulation of the somatosensory cortex

Objective. Intracortical microstimulation of the primary somatosensory cortex (S1) has shown great progress in restoring touch sensations to patients with paralysis. Stimulation parameters such as amplitude, phase duration, and frequency can influence the quality of the evoked percept as well as the amount of charge necessary to elicit a response. Previous studies in V1 and auditory cortices have shown that the behavioral responses to stimulation amplitude and phase duration change across cortical depth. However, this depth-dependent response has yet to be investigated in S1. Similarly, to our knowledge, the response to microstimulation frequency across cortical depth remains unexplored.Approach. To assess these questions, we implanted rats in S1 with a microelectrode with electrode-sites spanning all layers of the cortex. A conditioned avoidance behavioral paradigm was used to measure detection thresholds and responses to phase duration and frequency across cortical depth.Main results. Analogous to other cortical areas, the sensitivity to charge and strength-duration chronaxies in S1 varied across cortical layers. Likewise, the sensitivity to microstimulation frequency was layer dependent.Significance. These findings suggest that cortical depth can play an important role in the fine-tuning of stimulation parameters and in the design of intracortical neuroprostheses for clinical applications.

CIVET-Macaque: An automated pipeline for MRI-based cortical surface generation and cortical thickness in macaques

The MNI CIVET pipeline for automated extraction of cortical surfaces and evaluation of cortical thickness from in-vivo human MRI has been extended for processing macaque brains. Processing is performed based on the NIMH Macaque Template (NMT), as the reference template, with the anatomical parcellation of the surface following the D99 and CHARM atlases. The modifications needed to adapt CIVET to the macaque brain are detailed. Results have been obtained using CIVET-macaque to process the anatomical scans of the 31 macaques used to generate the NMT and another 95 macaques from the PRIME-DE initiative. It is anticipated that the open usage of CIVET-macaque will promote collaborative efforts in data collection and processing, sharing, and automated analyses from which the non-human primate brain imaging field will advance.

Systematically disrupted functional gradient of the cortical connectome in generalized epilepsy: Initial discovery and independent sample replication

Genetic generalized epilepsy is a network disorder typically involving distributed areas identified by classical neuroanatomy. However, the finer topological relationships in terms of continuous spatial arrangement between these systems are still ambiguous. Connectome gradients provide the topological representations of human macroscale hierarchy in an abstract low-dimensional space by embedding the functional connectome into a set of axes. Leveraging connectome gradients, we systematically scrutinized abnormalities of functional connectome gradient in patients with genetic generalized epilepsy with tonic-clonic seizure (GGE-GTCS, n = 78) compared to healthy controls (HC, n = 85), and further examined the reproducibility across multiple processing configurations and in an independent validation sample (patients with GGE-GTCS, n = 28; HC, n = 31). Our findings demonstrated an extended principal gradient at different spatial scales, network-level and vertex-level, in patients with GGE-GTCS. We found consistent results across processing parameters and in validation sample. The extended principal gradient revealed the excessive functional segregation between unimodal and transmodal systems associated with duration of epilepsy and age at seizure onset in patients. Furthermore, the connectivity profile of regions with abnormal principal gradients verified the disrupted functional hierarchy revealed by gradients. Together, our findings provided a novel view of functional system hierarchy alterations, which facilitated a continuous spatial arrangement of macroscale networks, to increase our understanding of the functional connectome hierarchy in generalized epilepsy.

The natural axis of transmitter receptor distribution in the human cerebral cortex

Transmitter receptors constitute a key component of the molecular machinery for intercellular communication in the brain. Recent efforts have mapped the density of diverse transmitter receptors across the human cerebral cortex with an unprecedented level of detail. Here, we distill these observations into key organizational principles. We demonstrate that receptor densities form a natural axis in the human cerebral cortex, reflecting decreases in differentiation at the level of laminar organization and a sensory-to-association axis at the functional level. Along this natural axis, key organizational principles are discerned: progressive molecular diversity (increase of the diversity of receptor density); excitation/inhibition (increase of the ratio of excitatory-to-inhibitory receptor density); and mirrored, orderly changes of the density of ionotropic and metabotropic receptors. The uncovered natural axis formed by the distribution of receptors aligns with the axis that is formed by other dimensions of cortical organization, such as the myelo- and cytoarchitectonic levels. Therefore, the uncovered natural axis constitutes a unifying organizational feature linking multiple dimensions of the cerebral cortex, thus bringing order to the heterogeneity of cortical organization.

Topographic gradients of intrinsic dynamics across neocortex

The intrinsic dynamics of neuronal populations are shaped by both microscale attributes and macroscale connectome architecture. Here we comprehensively characterize the rich temporal patterns of neural activity throughout the human brain. Applying massive temporal feature extraction to regional haemodynamic activity, we systematically estimate over 6000 statistical properties of individual brain regions' time-series across the neocortex. We identify two robust spatial gradients of intrinsic dynamics, one spanning a ventromedial-dorsolateral axis and dominated by measures of signal autocorrelation, and the other spanning a unimodal-transmodal axis and dominated by measures of dynamic range. These gradients reflect spatial patterns of gene expression, intracortical myelin and cortical thickness, as well as structural and functional network embedding. Importantly, these gradients are correlated with patterns of meta-analytic functional activation, differentiating cognitive versus affective processing and sensory versus higher-order cognitive processing. Altogether, these findings demonstrate a link between microscale and macroscale architecture, intrinsic dynamics, and cognition.

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.

Histology-Based Average Template of the Marmoset Cortex With Probabilistic Localization of Cytoarchitectural Areas

The rapid adoption of marmosets in neuroscience has created a demand for three dimensional (3D) atlases of the brain of this species to facilitate data integration in a common reference space. We report on a new open access template of the marmoset cortex (the Nencki-Monash, or NM template), representing a morphological average of 20 brains of young adult individuals, obtained by 3D reconstructions generated from Nissl-stained serial sections. The method used to generate the template takes into account morphological features of the individual brains, as well as the borders of clearly defined cytoarchitectural areas. This has resulted in a resource which allows direct estimates of the most likely coordinates of each cortical area, as well as quantification of the margins of error involved in assigning voxels to areas, and preserves quantitative information about the laminar structure of the cortex. We provide spatial transformations between the NM and other available marmoset brain templates, thus enabling integration with magnetic resonance imaging (MRI) and tracer-based connectivity data. The NM template combines some of the main advantages of histology-based atlases (e.g. information about the cytoarchitectural structure) with features more commonly associated with MRI-based templates (isotropic nature of the dataset, and probabilistic analyses). The underlying workflow may be found useful in the future development of 3D brain atlases that incorporate information about the variability of areas in species for which it may be impractical to ensure homogeneity of the sample in terms of age, sex and genetic background.

A Protocol for Cortical Type Analysis of the Human Neocortex Applied on Histological Samples, the Atlas of Von Economo and Koskinas, and Magnetic Resonance Imaging

The human cerebral cortex is parcellated in hundreds of areas using neuroanatomy and imaging methods. Alternatively, cortical areas can be classified into few cortical types according to their degree of laminar differentiation. Cortical type analysis is based on the gradual and systematic variation of laminar features observed across the entire cerebral cortex in Nissl stained sections and has profound implications for understanding fundamental aspects of evolution, development, connections, function, and pathology of the cerebral cortex. In this protocol paper, we explain the general principles of cortical type analysis and provide tables with the fundamental features of laminar structure that are studied for this analysis. We apply cortical type analysis to the micrographs of the Atlas of the human cerebral cortex of von Economo and Koskinas and provide tables and maps with the areas of this Atlas and their corresponding cortical type. Finally, we correlate the cortical type maps with the T1w/T2w ratio from widely used reference magnetic resonance imaging scans. The analysis, tables and maps of the human cerebral cortex shown in this protocol paper can be used to predict patterns of connections between areas according to the principles of the Structural Model and determine their level in cortical hierarchies. Cortical types can also predict the spreading of abnormal proteins in neurodegenerative diseases to the level of cortical layers. In summary, cortical type analysis provides a theoretical and practical framework for directed studies of connectivity, synaptic plasticity, and selective vulnerability to neurologic and psychiatric diseases in the human neocortex.

Cortical morphology at birth reflects spatiotemporal patterns of gene expression in the fetal human brain

Interruption to gestation through preterm birth can significantly impact cortical development and have long-lasting adverse effects on neurodevelopmental outcome. We compared cortical morphology captured by high-resolution, multimodal magnetic resonance imaging (MRI) in n = 292 healthy newborn infants (mean age at birth = 39.9 weeks) with regional patterns of gene expression in the fetal cortex across gestation (n = 156 samples from 16 brains, aged 12 to 37 postconceptional weeks [pcw]). We tested the hypothesis that noninvasive measures of cortical structure at birth mirror areal differences in cortical gene expression across gestation, and in a cohort of n = 64 preterm infants (mean age at birth = 32.0 weeks), we tested whether cortical alterations observed after preterm birth were associated with altered gene expression in specific developmental cell populations. Neonatal cortical structure was aligned to differential patterns of cell-specific gene expression in the fetal cortex. Principal component analysis (PCA) of 6 measures of cortical morphology and microstructure showed that cortical regions were ordered along a principal axis, with primary cortex clearly separated from heteromodal cortex. This axis was correlated with estimated tissue maturity, indexed by differential expression of genes expressed by progenitor cells and neurons, and engaged in stem cell differentiation, neuron migration, and forebrain development. Preterm birth was associated with altered regional MRI metrics and patterns of differential gene expression in glial cell populations. The spatial patterning of gene expression in the developing cortex was thus mirrored by regional variation in cortical morphology and microstructure at term, and this was disrupted by preterm birth. This work provides a framework to link molecular mechanisms to noninvasive measures of cortical development in early life and highlights novel pathways to injury in neonatal populations at increased risk of neurodevelopmental disorder.

Interruption to gestation through preterm birth can significantly impact cortical development and have long-lasting adverse effects on neurodevelopmental outcome. A large neuroimaging study of newborn infants reveals how their cortical structure at birth is associated with patterns of gene expression in the fetal cortex and how this relationship is affected by preterm birth.