Cyto- and receptor architectonic mapping of the human brain

Multimodal mapping of macaque monkey somatosensory cortex

The somatosensory cortex is a brain region responsible for receiving and processing sensory information from across the body and is structurally and functionally heterogeneous. Since the chemoarchitectonic segregation of the cerebral cortex can be revealed by transmitter receptor distribution patterns, by using a quantitative multireceptor architectonical analysis, we determined the number and extent of distinct areas of the macaque somatosensory cortex. We identified three architectonically distinct cortical entities within the primary somatosensory cortex (i.e., 3bm, 3bli, 3ble), four within the anterior parietal cortex (i.e., 3am, 3al, 1 and 2) and six subdivisions (i.e., S2l, S2m, PVl, PVm, PRl and PRm) within the lateral fissure. We provide an ultra-high resolution 3D atlas of macaque somatosensory areas in stereotaxic space, which integrates cyto- and receptor architectonic features of identified areas. Multivariate analyses of the receptor fingerprints revealed four clusters of identified areas based on the degree of (dis)similarity of their receptor architecture. Each of these clusters can be associated with distinct levels of somatosensory processing, further demonstrating that the functional segregation of cortical areas is underpinned by differences in their molecular organization.

Patient-specific models link neurotransmitter receptor mechanisms with motor and visuospatial axes of Parkinson's disease

Parkinson's disease involves multiple neurotransmitter systems beyond the classical dopaminergic circuit, but their influence on structural and functional alterations is not well understood. Here, we use patient-specific causal brain modeling to identify latent neurotransmitter receptor-mediated mechanisms contributing to Parkinson's disease progression. Combining the spatial distribution of 15 receptors from post-mortem autoradiography with 6 neuroimaging-derived pathological factors, we detect a diverse set of receptors influencing gray matter atrophy, functional activity dysregulation, microstructural degeneration, and dendrite and dopaminergic transporter loss. Inter-individual variability in receptor mechanisms correlates with symptom severity along two distinct axes, representing motor and psychomotor symptoms with large GABAergic and glutamatergic contributions, and cholinergically-dominant visuospatial, psychiatric and memory dysfunction. Our work demonstrates that receptor architecture helps explain multi-factorial brain re-organization, and suggests that distinct, co-existing receptor-mediated processes underlie Parkinson's disease.

Cytoarchitectonic, receptor distribution and functional connectivity analyses of the macaque frontal lobe

Based on quantitative cyto- and receptor architectonic analyses, we identified 35 prefrontal areas, including novel subdivisions of Walker's areas 10, 9, 8B, and 46. Statistical analysis of receptor densities revealed regional differences in lateral and ventrolateral prefrontal cortex. Indeed, structural and functional organization of subdivisions encompassing areas 46 and 12 demonstrated significant differences in the interareal levels of α2 receptors. Furthermore, multivariate analysis included receptor fingerprints of previously identified 16 motor areas in the same macaque brains and revealed 5 clusters encompassing frontal lobe areas. We used the MRI datasets from the non-human primate data sharing consortium PRIME-DE to perform functional connectivity analyses using the resulting frontal maps as seed regions. In general, rostrally located frontal areas were characterized by bigger fingerprints, that is, higher receptor densities, and stronger regional interconnections. Whereas more caudal areas had smaller fingerprints, but showed a widespread connectivity pattern with distant cortical regions. Taken together, this study provides a comprehensive insight into the molecular structure underlying the functional organization of the cortex and, thus, reconcile the discrepancies between the structural and functional hierarchical organization of the primate frontal lobe. Finally, our data are publicly available via the EBRAINS and BALSA repositories for the entire scientific community.

Gradients of neurotransmitter receptor expression in the macaque cortex

Dynamics and functions of neural circuits depend on interactions mediated by receptors. Therefore, a comprehensive map of receptor organization across cortical regions is needed. In this study, we used in vitro receptor autoradiography to measure the density of 14 neurotransmitter receptor types in 109 areas of macaque cortex. We integrated the receptor data with anatomical, genetic and functional connectivity data into a common cortical space. We uncovered a principal gradient of receptor expression per neuron. This aligns with the cortical hierarchy from sensory cortex to higher cognitive areas. A second gradient, driven by serotonin 5-HT1A receptors, peaks in the anterior cingulate, default mode and salience networks. We found a similar pattern of 5-HT1A expression in the human brain. Thus, the macaque may be a promising translational model of serotonergic processing and disorders. The receptor gradients may enable rapid, reliable information processing in sensory cortical areas and slow, flexible integration in higher cognitive areas.

Neuromodulatory control of complex adaptive dynamics in the brain

How is the massive dimensionality and complexity of the microscopic constituents of the nervous system brought under sufficiently tight control so as to coordinate adaptive behaviour? A powerful means for striking this balance is to poise neurons close to the critical point of a phase transition, at which a small change in neuronal excitability can manifest a nonlinear augmentation in neuronal activity. How the brain could mediate this critical transition is a key open question in neuroscience. Here, I propose that the different arms of the ascending arousal system provide the brain with a diverse set of heterogeneous control parameters that can be used to modulate the excitability and receptivity of target neurons-in other words, to act as control parameters for mediating critical neuronal order. Through a series of worked examples, I demonstrate how the neuromodulatory arousal system can interact with the inherent topological complexity of neuronal subsystems in the brain to mediate complex adaptive behaviour.

Relationships between neurotransmitter receptor densities and expression levels of their corresponding genes in the human hippocampus

Neurotransmitter receptors are key molecules in signal transmission, their alterations are associated with brain dysfunction. Relationships between receptors and their corresponding genes are poorly understood, especially in humans. We combined in vitro receptor autoradiography and RNA sequencing to quantify, in the same tissue samples (7 subjects), the densities of 14 receptors and expression levels of their corresponding 43 genes in the Cornu Ammonis (CA) and dentate gyrus (DG) of human hippocampus. Significant differences in receptor densities between both structures were found only for metabotropic receptors, whereas significant differences in RNA expression levels mostly pertained ionotropic receptors. Receptor fingerprints of CA and DG differ in shapes but have similar sizes; the opposite holds true for their "RNA fingerprints", which represent the expression levels of multiple genes in a single area. In addition, the correlation coefficients between receptor densities and corresponding gene expression levels vary widely and the mean correlation strength was weak-to-moderate. Our results suggest that receptor densities are not only controlled by corresponding RNA expression levels, but also by multiple regionally specific post-translational factors.

Mapping the primate thalamus: systematic approach to analyze the distribution of subcortical neuromodulatory afferents

Neuromodulatory afferents to thalamic nuclei are key for information transmission and thus play critical roles in sensory, motor, and limbic processes. Over the course of the last decades, diverse attempts have been made to map and describe subcortical neuromodulatory afferents to the primate thalamus, including axons using acetylcholine, serotonin, dopamine, noradrenaline, adrenaline, and histamine. Our group has been actively involved in this endeavor. The published descriptions on neuromodulatory afferents to the primate thalamus have been made in different laboratories and are not fully comparable due to methodological divergences (for example, fixation procedures, planes of cutting, techniques used to detect the afferents, different criteria for identification of thalamic nuclei…). Such variation affects the results obtained. Therefore, systematic methodological and analytical approaches are much needed. The present article proposes reproducible methodological and terminological frameworks for primate thalamic mapping. We suggest the use of standard stereotaxic planes to produce and present maps of the primate thalamus, as well as the use of the Anglo-American school terminology (vs. the German school terminology) for identification of thalamic nuclei. Finally, a public repository of the data collected under agreed-on frameworks would be a useful tool for looking up and comparing data on the structure and connections of primate thalamic nuclei. Important and agreed-on efforts are required to create, manage, and fund a unified and homogeneous resource of data on the primate thalamus. Likewise, a firm commitment of the institutions to preserve experimental brain material is much needed because neuroscience work with non-human primates is becoming increasingly rare, making earlier material still more valuable.

Correspondence between gene expression and neurotransmitter receptor and transporter density in the human brain

Neurotransmitter receptors modulate signaling between neurons. Thus, neurotransmitter receptors and transporters play a key role in shaping brain function. Due to the lack of comprehensive neurotransmitter receptor/transporter density datasets, microarray gene expression measuring mRNA transcripts is often used as a proxy for receptor densities. In the present report, we comprehensively test the spatial correlation between gene expression and protein density for a total of 27 neurotransmitter receptors, receptor binding-sites, and transporters across 9 different neurotransmitter systems, using both PET and autoradiography radioligand-based imaging modalities. We find poor spatial correspondences between gene expression and density for all neurotransmitter receptors and transporters except four single-protein metabotropic receptors (5-HT_1A, CB₁, D₂, and MOR). These expression-density associations are related to gene differential stability and can vary between cortical and subcortical structures. Altogether, we recommend using direct measures of receptor and transporter density when relating neurotransmitter systems to brain structure and function.

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.

Closing the mechanistic gap: the value of microarchitecture in understanding cognitive networks

Cognitive neuroscience aims to provide biologically relevant accounts of cognition. Contemporary research linking spatial patterns of neural activity to psychological constructs describes 'where' hypothesised functions occur, but not 'how' these regions contribute to cognition. Technological, empirical, and conceptual advances allow this mechanistic gap to be closed by embedding patterns of functional activity in macro- and microscale descriptions of brain organisation. Recent work on the default mode network (DMN) and the multiple demand network (MDN), for example, highlights a microarchitectural landscape that may explain how activity in these networks integrates varied information, thus providing an anatomical foundation that will help to explain how these networks contribute to many different cognitive states. This perspective highlights emerging insights into how microarchitecture can constrain network accounts of human cognition.

Receptor architecture of macaque and human early visual areas: not equal, but comparable

Existing cytoarchitectonic maps of the human and macaque posterior occipital cortex differ in the number of areas they display, thus hampering identification of homolog structures. We applied quantitative in vitro receptor autoradiography to characterize the receptor architecture of the primary visual and early extrastriate cortex in macaque and human brains, using previously published cytoarchitectonic criteria as starting point of our analysis. We identified 8 receptor architectonically distinct areas in the macaque brain (mV1d, mV1v, mV2d, mV2v, mV3d, mV3v, mV3A, mV4v), and their respective counterpart areas in the human brain (hV1d, hV1v, hV2d, hV2v, hV3d, hV3v, hV3A, hV4v). Mean densities of 14 neurotransmitter receptors were quantified in each area, and ensuing receptor fingerprints used for multivariate analyses. The 1st principal component segregated macaque and human early visual areas differ. However, the 2nd principal component showed that within each species, area-specific differences in receptor fingerprints were associated with the hierarchical processing level of each area. Subdivisions of V2 and V3 were found to cluster together in both species and were segregated from subdivisions of V1 and from V4v. Thus, comparative studies like this provide valuable architectonic insights into how differences in underlying microstructure impact evolutionary changes in functional processing of the primate brain and, at the same time, provide strong arguments for use of macaque monkey brain as a suitable animal model for translational studies.

A dopamine gradient controls access to distributed working memory in the large-scale monkey cortex

Dopamine is required for working memory, but how it modulates the large-scale cortex is unknown. Here, we report that dopamine receptor density per neuron, measured by autoradiography, displays a macroscopic gradient along the macaque cortical hierarchy. This gradient is incorporated in a connectome-based large-scale cortex model endowed with multiple neuron types. The model captures an inverted U-shaped dependence of working memory on dopamine and spatial patterns of persistent activity observed in over 90 experimental studies. Moreover, we show that dopamine is crucial for filtering out irrelevant stimuli by enhancing inhibition from dendrite-targeting interneurons. Our model revealed that an activity-silent memory trace can be realized by facilitation of inter-areal connections and that adjusting cortical dopamine induces a switch from this internal memory state to distributed persistent activity. Our work represents a cross-level understanding from molecules and cell types to recurrent circuit dynamics underlying a core cognitive function distributed across the primate cortex.

Personalized brain models identify neurotransmitter receptor changes in Alzheimer's disease

Alzheimer's disease (AD) involves many neurobiological alterations from molecular to macroscopic spatial scales, but we currently lack integrative, mechanistic brain models characterizing how factors across different biological scales interact to cause clinical deterioration in a way that is subject-specific or personalized. Neurotransmitter receptors, as important signaling molecules and potential drug targets, are key mediators of interactions between many neurobiological processes altered in AD. We present a neurotransmitter receptor-enriched multifactorial brain model, which integrates spatial distribution patterns of 15 neurotransmitter receptors from post-mortem autoradiography with multiple in-vivo neuroimaging modalities (tau, amyloid-β and glucose PET, and structural, functional and arterial spin labeling MRI) in a personalized, generative, whole-brain formulation. Applying this data-driven model to a heterogeneous aged population (N = 423, ADNI data), we observed that personalized receptor-neuroimaging interactions explained about 70% (± 20%) of the across-population variance in longitudinal changes to the six neuroimaging modalities, and up to 39.7% (P < 0.003, FWE-corrected) of inter-individual variability in AD cognitive deterioration via an axis primarily affecting executive function. Notably, based on their contribution to the clinical severity in AD, we found significant functional alterations to glutamatergic interactions affecting tau accumulation and neural activity dysfunction, and GABAergic interactions concurrently affecting neural activity dysfunction, amyloid and tau distributions, as well as significant cholinergic receptor effects on tau accumulation. Overall, GABAergic alterations had the largest effect on cognitive impairment (particularly executive function) in our AD cohort (N = 25). Furthermore, we demonstrate the clinical applicability of this approach by characterizing subjects based on individualized 'fingerprints' of receptor alterations. This study introduces the first robust, data-driven framework for integrating several neurotransmitter receptors, multi-modal neuroimaging and clinical data in a flexible and interpretable brain model. It enables further understanding of the mechanistic neuropathological basis of neurodegenerative progression and heterogeneity, and constitutes a promising step towards implementing personalized, neurotransmitter-based treatments.

In vivo synaptic density relates to glucose metabolism at rest in healthy subjects, but is strongly modulated by regional differences

Preclinical and postmortem studies have suggested that regional synaptic density and glucose consumption (CMRGlc) are strongly related. However, the relation between synaptic density and cerebral glucose metabolism in the human brain has not directly been assessed in vivo. Using [¹¹C]UCB-J binding to synaptic vesicle glycoprotein 2 A (SV2A) as indicator for synaptic density and [¹⁸F]FDG for measuring cerebral glucose consumption, we studied twenty healthy female subjects (age 29.6 ± 9.9 yrs) who underwent a single-day dual-tracer protocol (GE Signa PET-MR). Global measures of absolute and relative CMRGlc and specific binding of [¹¹C]UCB-J were indeed highly significantly correlated (r > 0.47, p < 0.001). However, regional differences in relative [¹⁸F]FDG and [¹¹C]UCB-J uptake were observed, with up to 19% higher [¹¹C]UCB-J uptake in the medial temporal lobe (MTL) and up to 17% higher glucose metabolism in frontal and motor-related areas and thalamus. This pattern has a considerable overlap with the brain regions showing different levels of aerobic glycolysis. Regionally varying energy demands of inhibitory and excitatory synapses at rest may also contribute to this difference. Being unaffected by astroglial and/or microglial energy demands, changes in synaptic density in the MTL may therefore be more sensitive to early detection of pathological conditions compared to changes in glucose metabolism.

Distribution of the Noradrenaline Innervation and Adrenoceptors in the Macaque Monkey Thalamus

Noradrenaline (NA) in the thalamus has important roles in physiological, pharmacological, and pathological neuromodulation. In this work, a complete characterization of NA axons and Alpha adrenoceptors distributions is provided.

NA axons, revealed by immunohistochemistry against the synthesizing enzyme and the NA transporter, are present in all thalamic nuclei. The most densely innervated ones are the midline nuclei, intralaminar nuclei (paracentral and parafascicular), and the medial sector of the mediodorsal nucleus (MDm). The ventral motor nuclei and most somatosensory relay nuclei receive a moderate NA innervation. The pulvinar complex receives a heterogeneous innervation. The lateral geniculate nucleus (GL) has the lowest NA innervation.

Alpha adrenoceptors were analyzed by in vitro quantitative autoradiography. Alpha-1 receptor densities are higher than Alpha-2 densities. Overall, axonal densities and Alpha adrenoceptor densities coincide; although some mismatches were identified. The nuclei with the highest Alpha-1 values are MDm, the parvocellular part of the ventral posterior medial nucleus, medial pulvinar, and midline nuclei. The nucleus with the lowest Alpha-1 receptor density is GL. Alpha-2 receptor densities are highest in the lateral dorsal, centromedian, medial and inferior pulvinar, and midline nuclei.

These results suggest a role for NA in modulating thalamic involvement in consciousness, limbic, cognitive, and executive functions.

Computational models link cellular mechanisms of neuromodulation to large-scale neural dynamics

Decades of neurobiological research have disclosed the diverse manners in which the response properties of neurons are dynamically modulated to support adaptive cognitive functions. This neuromodulation is achieved through alterations in the biophysical properties of the neuron. However, changes in cognitive function do not arise directly from the modulation of individual neurons, but are mediated by population dynamics in mesoscopic neural ensembles. Understanding this multiscale mapping is an important but nontrivial issue. Here, we bridge these different levels of description by showing how computational models parametrically map classic neuromodulatory processes onto systems-level models of neural activity. The ensuing critical balance of systems-level activity supports perception and action, although our knowledge of this mapping remains incomplete. In this way, quantitative models that link microscale neuronal neuromodulation to systems-level brain function highlight gaps in knowledge and suggest new directions for integrating theoretical and experimental work.

Organization of the macaque monkey inferior parietal lobule based on multimodal receptor architectonics

The macaque monkey inferior parietal lobe (IPL) is a structurally heterogeneous brain region, although the number of areas it contains and the anatomical/functional relationship of identified subdivisions remains controversial. Neurotransmitter receptor distribution patterns not only reveal the position of the cortical borders, but also segregate areas associated to different functional systems. Thus we carried out a multimodal quantitative analysis of the cyto- and receptor architecture of the macaque IPL to determine the number and extent of distinct areas it encompasses. We identified four areas on the IPL convexity arranged in a caudo-rostral sequence, as well as two areas in the parietal operculum, which we projected onto the Yerkes19 surface. We found rostral areas to have relatively smaller receptor fingerprints than the caudal ones, which is in an agreement with the functional gradient along the caudo-rostral axis described in previous studies. The hierarchical analysis segregated IPL areas into two clusters: the caudal one, contains areas involved in multisensory integration and visual-motor functions, and rostral cluster, encompasses areas active during motor planning and action-related functions. The results of the present study provide novel insights into clarifying the homologies between human and macaque IPL areas. The ensuing 3D map of the macaque IPL, and the receptor fingerprints are made publicly available to the neuroscientific community via the Human Brain Project and BALSA repositories for future cyto- and/or receptor architectonically driven analyses of functional imaging studies in non-human primates.

The Neurotransmitter Receptor Architecture of the Mouse Olfactory System

The uptake, transmission and processing of sensory olfactory information is modulated by inhibitory and excitatory receptors in the olfactory system. Previous studies have focused on the function of individual receptors in distinct brain areas, but the receptor architecture of the whole system remains unclear. Here, we analyzed the receptor profiles of the whole olfactory system of adult male mice. We examined the distribution patterns of glutamatergic (AMPA, kainate, mGlu_2/3, and NMDA), GABAergic (GABA_A, GABA_A(BZ), and GABA_B), dopaminergic (D_1/5) and noradrenergic (α₁ and α₂) neurotransmitter receptors by quantitative in vitro receptor autoradiography combined with an analysis of the cyto- and myelo-architecture. We observed that each subarea of the olfactory system is characterized by individual densities of distinct neurotransmitter receptor types, leading to a region- and layer-specific receptor profile. Thereby, the investigated receptors in the respective areas and strata showed a heterogeneous expression. Generally, we detected high densities of mGlu_2/3Rs, GABA_A(BZ)Rs and GABA_BRs. Noradrenergic receptors revealed a highly heterogenic distribution, while the dopaminergic receptor D_1/5 displayed low concentrations, except in the olfactory tubercle and the dorsal endopiriform nucleus. The similarities and dissimilarities of the area-specific multireceptor profiles were analyzed by a hierarchical cluster analysis. A three-cluster solution was found that divided the areas into the (1) olfactory relay stations (main and accessory olfactory bulb), (2) the olfactory cortex (anterior olfactory cortex, dorsal peduncular cortex, taenia tecta, piriform cortex, endopiriform nucleus, entorhinal cortex, orbitofrontal cortex) and the (3) olfactory tubercle, constituting its own cluster. The multimodal receptor-architectonic analysis of each component of the olfactory system provides new insights into its neurochemical organization and future possibilities for pharmaceutic targeting.

Bridging Scales in Alzheimer's Disease: Biological Framework for Brain Simulation With The Virtual Brain

Despite the acceleration of knowledge and data accumulation in neuroscience over the last years, the highly prevalent neurodegenerative disease of AD remains a growing problem. Alzheimer's Disease (AD) is the most common cause of dementia and represents the most prevalent neurodegenerative disease. For AD, disease-modifying treatments are presently lacking, and the understanding of disease mechanisms continues to be incomplete. In the present review, we discuss candidate contributing factors leading to AD, and evaluate novel computational brain simulation methods to further disentangle their potential roles. We first present an overview of existing computational models for AD that aim to provide a mechanistic understanding of the disease. Next, we outline the potential to link molecular aspects of neurodegeneration in AD with large-scale brain network modeling using The Virtual Brain (www.thevirtualbrain.org), an open-source, multiscale, whole-brain simulation neuroinformatics platform. Finally, we discuss how this methodological approach may contribute to the understanding, improved diagnostics, and treatment optimization of AD.

A Connectomic Hypothesis for the Hominization of the Brain

Cognitive abilities of the human brain, including language, have expanded dramatically in the course of our recent evolution from nonhuman primates, despite only minor apparent changes at the gene level. The hypothesis we propose for this paradox relies upon fundamental features of human brain connectivity, which contribute to a characteristic anatomical, functional, and computational neural phenotype, offering a parsimonious framework for connectomic changes taking place upon the human-specific evolution of the genome. Many human connectomic features might be accounted for by substantially increased brain size within the global neural architecture of the primate brain, resulting in a larger number of neurons and areas and the sparsification, increased modularity, and laminar differentiation of cortical connections. The combination of these features with the developmental expansion of upper cortical layers, prolonged postnatal brain development, and multiplied nongenetic interactions with the physical, social, and cultural environment gives rise to categorically human-specific cognitive abilities including the recursivity of language. Thus, a small set of genetic regulatory events affecting quantitative gene expression may plausibly account for the origins of human brain connectivity and cognition.

Multimodal 3D atlas of the macaque monkey motor and premotor cortex

In the present study we reevaluated the parcellation scheme of the macaque frontal agranular cortex by implementing quantitative cytoarchitectonic and multireceptor analyses, with the purpose to integrate and reconcile the discrepancies between previously published maps of this region. We applied an observer-independent and statistically testable approach to determine the position of cytoarchitectonic borders. Analysis of the regional and laminar distribution patterns of 13 different transmitter receptors confirmed the position of cytoarchitectonically identified borders. Receptor densities were extracted from each area and visualized as its "receptor fingerprint". Hierarchical and principal components analyses were conducted to detect clusters of areas according to the degree of (dis)similarity of their fingerprints. Finally, functional connectivity pattern of each identified area was analyzed with areas of prefrontal, cingulate, somatosensory and lateral parietal cortex and the results were depicted as "connectivity fingerprints" and seed-to-vertex connectivity maps. We identified 16 cyto- and receptor architectonically distinct areas, including novel subdivisions of the primary motor area 4 (i.e. 4a, 4p, 4m) and of premotor areas F4 (i.e. F4s, F4d, F4v), F5 (i.e. F5s, F5d, F5v) and F7 (i.e. F7d, F7i, F7s). Multivariate analyses of receptor fingerprints revealed three clusters, which first segregated the subdivisions of area 4 with F4d and F4s from the remaining premotor areas, then separated ventrolateral from dorsolateral and medial premotor areas. The functional connectivity analysis revealed that medial and dorsolateral premotor and motor areas show stronger functional connectivity with areas involved in visual processing, whereas 4p and ventrolateral premotor areas presented a stronger functional connectivity with areas involved in somatomotor responses. For the first time, we provide a 3D atlas integrating cyto- and multi-receptor architectonic features of the macaque motor and premotor cortex. This atlas constitutes a valuable resource for the analysis of functional experiments carried out with non-human primates, for modeling approaches with realistic synaptic dynamics, as well as to provide insights into how brain functions have developed by changes in the underlying microstructure and encoding strategies during evolution.

Anatomy of nerve fiber bundles at micrometer-resolution in the vervet monkey visual system

Although the primate visual system has been extensively studied, detailed spatial organization of white matter fiber tracts carrying visual information between areas has not been fully established. This is mainly due to the large gap between tracer studies and diffusion-weighted MRI studies, which focus on specific axonal connections and macroscale organization of fiber tracts, respectively. Here we used 3D polarization light imaging (3D-PLI), which enables direct visualization of fiber tracts at micrometer resolution, to identify and visualize fiber tracts of the visual system, such as stratum sagittale, inferior longitudinal fascicle, vertical occipital fascicle, tapetum and dorsal occipital bundle in vervet monkey brains. Moreover, 3D-PLI data provide detailed information on cortical projections of these tracts, distinction between neighboring tracts, and novel short-range pathways. This work provides essential information for interpretation of functional and diffusion-weighted MRI data, as well as revision of wiring diagrams based upon observations in the vervet visual system.

Receptor-driven, multimodal mapping of cortical areas in the macaque monkey intraparietal sulcus

The intraparietal sulcus (IPS) is structurally and functionally heterogeneous. We performed a quantitative cyto-/myelo- and receptor architectonical analysis to provide a multimodal map of the macaque IPS. We identified 17 cortical areas, including novel areas PEipe, PEipi (external and internal subdivisions of PEip), and MIPd. Multivariate analyses of receptor densities resulted in a grouping of areas based on the degree of (dis)similarity of their receptor architecture: a cluster encompassing areas located in the posterior portion of the IPS and associated mainly with the processing of visual information, a cluster including areas found in the anterior portion of the IPS and involved in sensorimotor processing, and an ‘intermediate’ cluster of multimodal association areas. Thus, differences in cyto-/myelo- and receptor architecture segregate the cortical ribbon within the IPS, and receptor fingerprints provide novel insights into the relationship between the structural and functional segregation of this brain region in the macaque monkey.

Spatiotemporal Differences in the Regional Cortical Plate and Subplate Volume Growth during Fetal Development

The regional specification of the cerebral cortex can be described by protomap and protocortex hypotheses. The protomap hypothesis suggests that the regional destiny of cortical neurons and the relative size of the cortical area are genetically determined early during embryonic development. The protocortex hypothesis suggests that the regional growth rate is predominantly shaped by external influences. In order to determine regional volumes of cortical compartments (cortical plate (CP) or subplate (SP)) and estimate their growth rates, we acquired T2-weighted in utero MRIs of 40 healthy fetuses and grouped them into early (<25.5 GW), mid- (25.5-31.6 GW), and late (>31.6 GW) prenatal periods. MRIs were segmented into CP and SP and further parcellated into 22 gyral regions. No significant difference was found between periods in regional volume fractions of the CP or SP. However, during the early and mid-prenatal periods, we found significant differences in relative growth rates (% increase per GW) between regions of cortical compartments. Thus, the relative size of these regions are most likely conserved and determined early during development whereas more subtle growth differences between regions are fine-tuned later, during periods of peak thalamocortical growth. This is in agreement with both the protomap and protocortex hypothesis.

Developmental Changes of Glutamate and GABA Receptor Densities in Wistar Rats

Neurotransmitters and their receptors are key molecules of signal transduction and subject to various changes during pre- and postnatal development. Previous studies addressed ontogeny at the level of neurotransmitters and expression of neurotransmitter receptor subunits. However, developmental changes in receptor densities to this day are not well understood. Here, we analyzed developmental changes in excitatory glutamate and inhibitory γ-aminobutyric acid (GABA) receptors in adjacent sections of the rat brain by means of quantitative in vitro receptor autoradiography. Receptor densities of the ionotropic glutamatergic receptors α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA), kainate and N-methyl-D-aspartate (NMDA) as well as of the ionotropic GABA_A and metabotropic GABA_B receptors were investigated using specific high-affinity ligands. For each receptor binding site, significant density differences were demonstrated in the investigated regions of interest [olfactory bulb, striatum, hippocampus, and cerebellum] and developmental stages [postnatal day (P) 0, 10, 20, 30 and 90]. In particular, we showed that the glutamatergic and GABAergic receptor densities were already present between P0 and P10 in all regions of interest, which may indicate the early relevance of these receptors for brain development. A transient increase of glutamatergic receptor densities in the hippocampus was found, indicating their possible involvement in synaptic plasticity. We demonstrated a decline of NMDA receptor densities in the striatum and hippocampus from P30 to P90, which could be due to synapse elimination, a process that redefines neuronal networks in postnatal brains. Furthermore, the highest increase in GABA_A receptor densities from P10 to P20 coincides with the developmental shift from excitatory to inhibitory GABA transmission. Moreover, the increase from P10 to P20 in GABA_A receptor densities in the cerebellum corresponds to a point in time when functional GABAergic synapses are formed. Taken together, the present data reveal differential changes in glutamate and GABA receptor densities during postnatal rat brain development, which may contribute to their specific functions during ontogenesis, thus providing a deeper understanding of brain ontogenesis and receptor function.

Receptor density pattern confirms and enhances the anatomic-functional features of the macaque superior parietal lobule areas

The macaque monkey superior parietal lobule (SPL) is part of a neuronal network involved in the integration of information from visual and somatosensory cortical areas for execution of reaching and grasping movements. We applied quantitative in vitro receptor autoradiography to analyse the distribution patterns of 15 different receptors for glutamate, GABA, acetylcholine, serotonin, dopamine, and adenosine in the SPL of three adult male Macaca fascicularis monkeys. For each area, mean (averaged over all cortical layers) receptor densities were visualized as a receptor fingerprint of that area. Multivariate analyses were conducted to detect clusters of areas according to the degree of (dis)similarity of their receptor organization. Differences in regional and laminar receptor distributions confirm the location and extent of areas V6, V6Av, V6Ad, PEc, PEci, and PGm as found in cytoarchitectonic and functional studies, but also enable the definition of three subdivisions within area PE. Receptor densities are higher in supra- than in infragranular layers, with the exception of kainate, M₂, and adenosine receptors. Glutamate and GABAergic receptors are the most expressed in all areas analysed. Hierarchical cluster analyses demonstrate that SPL areas are organized in two groups, an organization that corresponds to the visual or sensory-motor characteristics of those areas. Finally, based on present results and in the framework of our current understanding of the structural and functional organization of the primate SPL, we propose a novel pattern of homologies between human and macaque SPL areas.

Human visual cortex is organized along two genetically opposed hierarchical gradients with unique developmental and evolutionary origins

Human visual cortex is organized with striking consistency across individuals. While recent findings demonstrate an unexpected coupling between functional and cytoarchitectonic regions relative to the folding of human visual cortex, a unifying principle linking these anatomical and functional features of the cortex remains elusive. To fill this gap in knowledge, we combined independent and ground truth measurements of cytoarchitectonic regions and genetic tissue characterization within human occipitotemporal cortex. Using a data-driven approach, we examined whether differential gene expression among cytoarchitectonic areas could contribute to the arealization of occipitotemporal cortex into a hierarchy based on transcriptomics. This approach revealed two opposing gene expression gradients: one that contains a series of genes with expression magnitudes that ascend from posterior (e.g., areas human occipital [hOc]1, hOc2, hOc3, etc.) to anterior cytoarchitectonic areas (e.g., areas fusiform gyrus [FG]1–FG4) and another that contains a separate series of genes that show a descending gradient from posterior to anterior areas. Using data from the living human brain, we show that each of these gradients correlates strongly with variations in measures related to either thickness or myelination of cortex, respectively. We further reveal that these genetic gradients emerge along unique trajectories in human development: the ascending gradient is present at 10–12 gestational weeks, while the descending gradient emerges later (19–24 gestational weeks). Interestingly, it is not until early childhood (before 5 years of age) that the two expression gradients achieve their adult-like mean expression values. Additional analyses in nonhuman primates (NHPs) reveal that homologous genes do not generate the same ascending and descending expression gradients as in humans. We discuss these findings relative to previously proposed hierarchies based on functional and cytoarchitectonic features of visual cortex. Altogether, these findings bridge macroscopic features of human cytoarchitectonic areas in visual cortex with microscopic features of cellular organization and genetic expression, which, despite the complexity of this multiscale correspondence, can be described by a sparse subset (approximately 200) of genes. These findings help pinpoint the genes contributing to healthy cortical development and explicate the cortical biology distinguishing humans from other primates, as well as establishing essential groundwork for understanding future work linking genetic mutations with the function and development of the human brain.

The expression of a sparse subset of human genes forms two opposed gradients that capture the processing hierarchy of visual cortex; these transcription gradients emerge at different points during human development and distinguish human from nonhuman primates.

Myelin and Modeling: Bootstrapping Cortical Microcircuits

Histological studies of myelin-stained sectioned cadaver brain and in vivo myelin-weighted magnetic resonance imaging (MRI) show that the cerebral cortex is organized into cortical areas with generally well-defined boundaries, which have consistent internal patterns of myelination. The process of myelination is largely driven by neural experience, in which the axonal passage of action potentials stimulates neighboring oligodendrocytes to perform their task. This bootstrapping process, such that the traffic of action potentials facilitates increased traffic, suggests the hypothesis that the specific pattern of myelination (myeloarchitecture) in each cortical area reveals the principal cortical microcircuits required for the function of that area. If this idea is correct, the observable sequential maturation of specific brain areas can provide evidence for models of the stages of cognitive development.

Rostral Anterior Cingulate Glutamine/Glutamate Disbalance in Major Depressive Disorder Depends on Symptom Severity

BACKGROUND

Patients with major depressive disorder (MDD) show glutamatergic deficits in the ventral anterior cingulate cortex. The glutamine/glutamate (Gln/Glu) ratio was proposed to be connected to glutamatergic cycling, which is hypothesized to be dysregulated in MDD. As an indicator of regional metabolite status, this ratio might be a robust state marker sensitive to clinical heterogeneity.

METHODS

Thirty-two MDD patients (mean age 40.88 ± 13.66 years, 19 women) and control subjects (mean age 33.09 ± 8.24 years, 19 women) were compared for pregenual anterior cingulate cortex levels of Gln/Glu, Gln/total creatine (tCr), Glu/tCr, and gamma-aminobutyric acid/tCr as determined by high-field magnetic resonance spectroscopy. We tested if symptom severity (Hamilton Depression Rating Scale) and anhedonia (Snaith-Hamilton Pleasure Scale) influence the relation of metabolites to clinical symptoms.

RESULTS

MDD patients showed higher Gln/Glu. This was driven by marginally higher Gln/tCr and nonsignificantly lower Glu/tCr. Groups defined by severity moderated relationship between Gln/Glu and the Hamilton Depression Rating Scale. Moreover, severe cases differed from both control subjects and moderate cases. Groups defined by the Snaith-Hamilton Pleasure Scale also displayed differential relationship between Gln/Glu and levels of anhedonia, predominantly driven by Gln/tCr.

CONCLUSIONS

We elaborate previous accounts of metabolite deficits in the anterior cingulate cortex toward increased Gln/Glu. There is a moderated relationship between severity and the ratio, which suggests consideration of different mechanisms or disease state for the respective subgroups in future studies.

The Structural Model: a theory linking connections, plasticity, pathology, development and evolution of the cerebral cortex

The classical theory of cortical systematic variation has been independently described in reptiles, monotremes, marsupials and placental mammals, including primates, suggesting a common bauplan in the evolution of the cortex. The Structural Model is based on the systematic variation of the cortex and is a platform for advancing testable hypotheses about cortical organization and function across species, including humans. The Structural Model captures the overall laminar structure of areas by dividing the cortical architectonic continuum into discrete categories (cortical types), which can be used to test hypotheses about cortical organization. By type, the phylogenetically ancient limbic cortices-which form a ring at the base of the cerebral hemisphere-are agranular if they lack layer IV, or dysgranular if they have an incipient granular layer IV. Beyond the dysgranular areas, eulaminate type cortices have six layers. The number and laminar elaboration of eulaminate areas differ depending on species or cortical system within a species. The construct of cortical type retains the topology of the systematic variation of the cortex and forms the basis for a predictive Structural Model, which has successfully linked cortical variation to the laminar pattern and strength of cortical connections, the continuum of plasticity and stability of areas, the regularities in the distribution of classical and novel markers, and the preferential vulnerability of limbic areas to neurodegenerative and psychiatric diseases. The origin of cortical types has been recently traced to cortical development, and helps explain the variability of diseases with an onset in ontogeny.