Structural Attributes and Principles of the Neocortical Connectome in the Marmoset Monkey

From comparative connectomics to large-scale working memory modeling in macaque and marmoset

Although macaques and marmosets are both primates of choice for studying the brain mechanisms of cognition, they differ in key aspects of anatomy and behavior. Interestingly, recent connectomic analysis revealed that strong top-down projections from the prefrontal cortex to the posterior parietal cortex, present in macaques and important for executive function, are absent in marmosets. Here, we propose a consensus mapping that bridges the two species' cortical atlases and allows for direct area-to-area comparison of their connectomes, which are then used to build comparative computational large-scale modeling of the frontoparietal circuit for working memory. We found that the macaque model exhibits resilience against distractors, a prerequisite for normal working memory function. By contrast, the marmoset model is sensitive to distractibility commonly observed behaviorally in this species. Surprisingly, this contrasting trend can be swapped by scaling intrafrontal and frontoparietal connections. Finally, the relevance to primate ethology and evolution is discussed.

Complex harmonics reveal low-dimensional manifolds of critical brain dynamics

The brain needs to perform time-critical computations to ensure survival. A potential solution lies in the nonlocal, distributed computation at the whole-brain level made possible by criticality and amplified by the rare long-range connections found in the brain's unique anatomical structure. This nonlocality can be captured by the mathematical structure of Schrödinger's wave equation, which is at the heart of the complex harmonics decomposition (CHARM) framework that performs the necessary dimensional manifold reduction able to extract nonlocality in critical spacetime brain dynamics. Using a large neuroimaging dataset of over 1000 people, CHARM captured the critical, nonlocal and long-range nature of brain dynamics and the underlying mechanisms were established using a precise whole-brain model. Equally, CHARM revealed the significantly different critical dynamics of wakefulness and sleep. Overall, CHARM is a promising theoretical framework for capturing the low-dimensionality of the complex network dynamics observed in neuroscience and provides evidence that networks of brain regions rather than individual brain regions are the key computational engines of critical brain dynamics.

Multi-scale spiking network model of human cerebral cortex

Although the structure of cortical networks provides the necessary substrate for their neuronal activity, the structure alone does not suffice to understand the activity. Leveraging the increasing availability of human data, we developed a multi-scale, spiking network model of human cortex to investigate the relationship between structure and dynamics. In this model, each area in one hemisphere of the Desikan-Killiany parcellation is represented by a $1\,\mathrm{mm^{2}}$ column with a layered structure. The model aggregates data across multiple modalities, including electron microscopy, electrophysiology, morphological reconstructions, and diffusion tensor imaging, into a coherent framework. It predicts activity on all scales from the single-neuron spiking activity to the area-level functional connectivity. We compared the model activity with human electrophysiological data and human resting-state functional magnetic resonance imaging (fMRI) data. This comparison reveals that the model can reproduce aspects of both spiking statistics and fMRI correlations if the inter-areal connections are sufficiently strong. Furthermore, we study the propagation of a single-spike perturbation and macroscopic fluctuations through the network. The open-source model serves as an integrative platform for further refinements and future in silico studies of human cortical structure, dynamics, and function.

Horizontal cortical connections shape intrinsic traveling waves into feature-selective motifs that regulate perceptual sensitivity

Intrinsic cortical activity forms traveling waves that modulate sensory-evoked responses and perceptual sensitivity. These intrinsic traveling waves (iTWs) may arise from the coordination of synaptic activity through long-range feature-dependent horizontal connectivity within cortical areas. In a spiking network model that incorporates feature-selective patchy connections, we observe iTW motifs that result from shifts in excitatory/inhibitory balance as action potentials traverse these patchy connections. To test whether feature-selective motifs occur in vivo, we examined data recorded in the middle temporal visual area (Area MT) of marmosets performing a visual detection task. We find that some iTWs form motifs that are feature selective, exhibiting direction-selective modulations in spiking activity. Further, motifs modulate the gain of target-evoked responses and perceptual sensitivity if the target matches the preference of the motif. These results suggest that iTWs are shaped by the patchy horizontal fiber projections in the cortex and can regulate neural and perceptual sensitivity in a feature-selective manner.

Bilateral connections from the amygdala to extrastriate visual cortex in the marmoset monkey

It is known that the primate amygdala forms projections to many areas of the ipsilateral cortex, but the extent to which it forms connections with the contralateral visual cortex remains less understood. Based on retrograde tracer injections in marmoset monkeys, we report that the amygdala forms widespread projections to the ipsilateral extrastriate cortex, including V1 and areas in both the dorsal (MT, V4T, V3a, 19M, and PG/PFG) and the ventral (VLP and TEO) streams. In addition, contralateral projections were found to target each of the extrastriate areas, but not V1. In both hemispheres, the tracer-labeled neurons were exclusively located in the basolateral nuclear complex. The number of labeled neurons in the contralateral amygdala was small relative to the ipsilateral connection (1.2% to 5.8%). The percentage of contralateral connections increased progressively with hierarchical level. An injection in the corpus callosum demonstrated that at least some of the amygdalo-cortical connections cross through this fiber tract, in addition to the previously documented path through the anterior commissure. Our results expand knowledge of the amygdalofugal projections to the extrastriate cortex, while also revealing pathways through which visual stimuli conveying affective content can directly influence early stages of neural processing in the contralateral visual field.

Distribution of calbindin-positive neurons across areas and layers of the marmoset cerebral cortex

The diversity of the mammalian cerebral cortex demands technical approaches to map the spatial distribution of neurons with different biochemical identities. This issue is magnified in the case of the primate cortex, characterized by a large number of areas with distinctive cytoarchitectures. To date, no full map of the distribution of cells expressing a specific protein has been reported for the cortex of any primate. Here we have charted the 3-dimensional distribution of neurons expressing the calcium-binding protein calbindin (CB+ neurons) across the entire marmoset cortex, using a combination of immunohistochemistry, automated cell identification, computerized reconstruction, and cytoarchitecture-aware registration. CB+ neurons formed a heterogeneous population, which together corresponded to 10-20% of the cortical neurons. They occurred in higher proportions in areas corresponding to low hierarchical levels of processing, such as sensory cortices. Although CB+ neurons were concentrated in the supragranular and granular layers, there were clear global trends in their laminar distribution. For example, their relative density in infragranular layers increased with hierarchical level along sensorimotor processing streams, and their density in layer 4 was lower in areas involved in sensorimotor integration, action planning and motor control. These results reveal new quantitative aspects of the cytoarchitectural organization of the primate cortex, and demonstrate an approach to mapping the full distribution of neurochemically distinct cells throughout the brain which is readily applicable to most other mammalian species.

Individualized epidemic spreading models predict epilepsy surgery outcomes: A pseudo-prospective study

Epilepsy surgery is the treatment of choice for drug-resistant epilepsy patients, but up to 50% of patients continue to have seizures one year after the resection. In order to aid presurgical planning and predict postsurgical outcome on a patient-by-patient basis, we developed a framework of individualized computational models that combines epidemic spreading with patient-specific connectivity and epileptogeneity maps: the Epidemic Spreading Seizure and Epilepsy Surgery framework (ESSES). ESSES parameters were fitted in a retrospective study (N = 15) to reproduce invasive electroencephalography (iEEG)-recorded seizures. ESSES reproduced the iEEG-recorded seizures, and significantly better so for patients with good (seizure-free, SF) than bad (nonseizure-free, NSF) outcome. We illustrate here the clinical applicability of ESSES with a pseudo-prospective study (N = 34) with a blind setting (to the resection strategy and surgical outcome) that emulated presurgical conditions. By setting the model parameters in the retrospective study, ESSES could be applied also to patients without iEEG data. ESSES could predict the chances of good outcome after any resection by finding patient-specific model-based optimal resection strategies, which we found to be smaller for SF than NSF patients, suggesting an intrinsic difference in the network organization or presurgical evaluation results of NSF patients. The actual surgical plan overlapped more with the model-based optimal resection, and had a larger effect in decreasing modeled seizure propagation, for SF patients than for NSF patients. Overall, ESSES could correctly predict 75% of NSF and 80.8% of SF cases pseudo-prospectively. Our results show that individualised computational models may inform surgical planning by suggesting alternative resections and providing information on the likelihood of a good outcome after a proposed resection. This is the first time that such a model is validated with a fully independent cohort and without the need for iEEG recordings.

Network analysis of marmoset cortical connections reveals pFC and sensory clusters

A new analysis is presented of the retrograde tracer measurements of connections between anatomical areas of the marmoset cortex. The original normalisation of raw data yields the fractional link weight measure, FLNe. That is re-examined to consider other possible measures that reveal the underlying in link weights. Predictions arising from both are used to examine network modules and hubs. With inclusion of the in weights the InfoMap algorithm identifies eight structural modules in marmoset cortex. In and out hubs and major connector nodes are identified using module assignment and participation coefficients. Time evolving network tracing around the major hubs reveals medium sized clusters in pFC, temporal, auditory and visual areas; the most tightly coupled and significant of which is in the pFC. A complementary viewpoint is provided by examining the highest traffic links in the cortical network, and reveals parallel sensory flows to pFC and via association areas to frontal areas.

The meso-connectomes of mouse, marmoset, and macaque: network organization and the emergence of higher cognition

The recent publications of the inter-areal connectomes for mouse, marmoset, and macaque cortex have allowed deeper comparisons across rodent vs. primate cortical organization. In general, these show that the mouse has very widespread, "all-to-all" inter-areal connectivity (i.e. a "highly dense" connectome in a graph theoretical framework), while primates have a more modular organization. In this review, we highlight the relevance of these differences to function, including the example of primary visual cortex (V1) which, in the mouse, is interconnected with all other areas, therefore including other primary sensory and frontal areas. We argue that this dense inter-areal connectivity benefits multimodal associations, at the cost of reduced functional segregation. Conversely, primates have expanded cortices with a modular connectivity structure, where V1 is almost exclusively interconnected with other visual cortices, themselves organized in relatively segregated streams, and hierarchically higher cortical areas such as prefrontal cortex provide top-down regulation for specifying precise information for working memory storage and manipulation. Increased complexity in cytoarchitecture, connectivity, dendritic spine density, and receptor expression additionally reveal a sharper hierarchical organization in primate cortex. Together, we argue that these primate specializations permit separable deconstruction and selective reconstruction of representations, which is essential to higher cognition.

Decomposing cortical activity through neuronal tracing connectome-eigenmodes in marmosets

Deciphering the complex relationship between neuroanatomical connections and functional activity in primate brains remains a daunting task, especially regarding the influence of monosynaptic connectivity on cortical activity. Here, we investigate the anatomical-functional relationship and decompose the neuronal-tracing connectome of marmoset brains into a series of eigenmodes using graph signal processing. These cellular connectome eigenmodes effectively constrain the cortical activity derived from resting-state functional MRI, and uncover a patterned cellular-functional decoupling. This pattern reveals a spatial gradient from coupled dorsal-posterior to decoupled ventral-anterior cortices, and recapitulates micro-structural profiles and macro-scale hierarchical cortical organization. Notably, these marmoset-derived eigenmodes may facilitate the inference of spontaneous cortical activity and functional connectivity of homologous areas in humans, highlighting the potential generalizing of the connectomic constraints across species. Collectively, our findings illuminate how neuronal-tracing connectome eigenmodes constrain cortical activity and improve our understanding of the brain's anatomical-functional relationship.

Horizontal cortical connections shape intrinsic traveling waves into feature-selective motifs that regulate perceptual sensitivity

Intrinsic, ongoing fluctuations of cortical activity form traveling waves that modulate the gain of sensory-evoked responses and perceptual sensitivity. Several lines of evidence suggest that intrinsic traveling waves (iTWs) may arise, in part, from the coordination of synaptic activity through the recurrent horizontal connectivity within cortical areas, which include long range patchy connections that link neurons with shared feature preferences. In a spiking network model with anatomical topology that incorporates feature-selective patchy connections, which we call the Balanced Patchy Network (BPN), we observe repeated iTWs, which we refer to as motifs. In the model, motifs stem from fluctuations in the excitability of like-tuned neurons that result from shifts in E/I balance as action potentials traverse these patchy connections. To test if feature-selective motifs occur in vivo, we examined data previously recorded using multielectrode arrays in Area MT of marmosets trained to perform a threshold visual detection task. Using a newly developed method for comparing the similarity of wave patterns we found that some iTWs can be grouped into motifs. As predicted by the BPN, many of these motifs are feature selective, exhibiting direction-selective modulations in ongoing spiking activity. Further, motifs modulate the gain of the response evoked by a target and perceptual sensitivity to the target if the target matches the preference of the motif. These results provide evidence that iTWs are shaped by the patterns of horizontal fiber projections in the cortex and that patchy connections enable iTWs to regulate neural and perceptual sensitivity in a feature selective manner.

Towards a biologically annotated brain connectome

The brain is a network of interleaved neural circuits. In modern connectomics, brain connectivity is typically encoded as a network of nodes and edges, abstracting away the rich biological detail of local neuronal populations. Yet biological annotations for network nodes - such as gene expression, cytoarchitecture, neurotransmitter receptors or intrinsic dynamics - can be readily measured and overlaid on network models. Here we review how connectomes can be represented and analysed as annotated networks. Annotated connectomes allow us to reconceptualize architectural features of networks and to relate the connection patterns of brain regions to their underlying biology. Emerging work demonstrates that annotated connectomes help to make more veridical models of brain network formation, neural dynamics and disease propagation. Finally, annotations can be used to infer entirely new inter-regional relationships and to construct new types of network that complement existing connectome representations. In summary, biologically annotated connectomes offer a compelling way to study neural wiring in concert with local biological features.

Pattern of ventral temporal lobe interconnections in rhesus macaques

The entorhinal cortex (EC, A28) is linked through reciprocal pathways with nearby perirhinal and visual, auditory, and multimodal association cortices in the temporal lobe, in pathways associated with the flow of information for memory processing. The density and laminar organization of these pathways is not well understood in primates. We studied interconnections within the ventral temporal lobe in young adult rhesus monkeys of both sexes with the aid of neural tracers injected in temporal areas (Ts1, Ts2, TE1, area 36, temporal polar area TPro, and area 28) to determine the density and laminar distribution of projection neurons within the temporal lobe. These temporal areas can be categorized into three different cortical types based on their laminar architecture: the sensory association areas Ts1, Ts2, and TE1 have six layers (eulaminate); the perirhinal limbic areas TPro and area 36 have an incipient layer IV (dysgranular); and area 28 lacks layer IV (agranular). We found that (1) temporal areas that are similar in laminar architecture by cortical type are strongly interconnected, and (2) the laminar pattern of connections is dependent on the difference in cortical laminar structure between linked areas. Thus, agranular A28 is more strongly connected with other agranular/dysgranular areas than with eulaminate cortices. Further, A28 predominantly projected via feedback-like pathways that originated in the deep layers, and received feedforward-like projections from areas of greater laminar differentiation, which emanated from the upper layers. Our results are consistent with the Structural Model, which relates the density and laminar distribution of connections to the relationship of the laminar structure between the linked areas. These connections were viewed in the context of the inhibitory microenvironment of A28, which is the key recipient of pathways from the cortex and of the output of hippocampus. Our findings revealed a higher population of calretinin (CR)-expressing neurons in EC, with a significantly higher density in its lateral division. Medial EC had a higher density of CR neurons in the deep layers, particularly in layer Va. In contrast, parvalbumin (PV) neurons were more densely distributed in the deep layers of the lateral subdivisions of rostral EC, especially in layer Va, whereas the densities of calbindin (CB) neurons in the medial and lateral EC were comparable in all layers, except for layer IIIa, in which medial EC had a higher CB population than the lateral. The pattern of connections in the inhibitory microenvironment of EC, which sends and receives input from the hippocampus, may shed light on signal propagation in this network associated with diverse aspects of memory, and disruptions in neurologic and psychiatric diseases that affect this region.

The regional variation of laminar thickness in the human isocortex is related to cortical hierarchy and interregional connectivity

The human isocortex consists of tangentially organized layers with unique cytoarchitectural properties. These layers show spatial variations in thickness and cytoarchitecture across the neocortex, which is thought to support function through enabling targeted corticocortical connections. Here, leveraging maps of the 6 cortical layers based on 3D human brain histology, we aimed to quantitatively characterize the systematic covariation of laminar structure in the cortex and its functional consequences. After correcting for the effect of cortical curvature, we identified a spatial pattern of changes in laminar thickness covariance from lateral frontal to posterior occipital regions, which differentiated the dominance of infra- versus supragranular layer thickness. Corresponding to the laminar regularities of cortical connections along cortical hierarchy, the infragranular-dominant pattern of laminar thickness was associated with higher hierarchical positions of regions, mapped based on resting-state effective connectivity in humans and tract-tracing of structural connections in macaques. Moreover, we show that regions with similar laminar thickness patterns have a higher likelihood of structural connections and strength of functional connections. In sum, here we characterize the organization of laminar thickness in the human isocortex and its association with cortico-cortical connectivity, illustrating how laminar organization may provide a foundational principle of cortical function.

A marmoset brain cell census reveals regional specialization of cellular identities

The mammalian brain is composed of many brain structures, each with its own ontogenetic and developmental history. We used single-nucleus RNA sequencing to sample over 2.4 million brain cells across 18 locations in the common marmoset, a New World monkey primed for genetic engineering, and examined gene expression patterns of cell types within and across brain structures. The adult transcriptomic identity of most neuronal types is shaped more by developmental origin than by neurotransmitter signaling repertoire. Quantitative mapping of GABAergic types with single-molecule FISH (smFISH) reveals that interneurons in the striatum and neocortex follow distinct spatial principles, and that lateral prefrontal and other higher-order cortical association areas are distinguished by high proportions of VIP+ neurons. We use cell type-specific enhancers to drive AAV-GFP and reconstruct the morphologies of molecularly resolved interneuron types in neocortex and striatum. Our analyses highlight how lineage, local context, and functional class contribute to the transcriptional identity and biodistribution of primate brain cell types.

Winning the Lottery With Neural Connectivity Constraints: Faster Learning Across Cognitive Tasks With Spatially Constrained Sparse RNNs

Recurrent neural networks (RNNs) are often used to model circuits in the brain and can solve a variety of difficult computational problems requiring memory, error correction, or selection (Hopfield, 1982; Maass et al., 2002; Maass, 2011). However, fully connected RNNs contrast structurally with their biological counterparts, which are extremely sparse (about 0.1%). Motivated by the neocortex, where neural connectivity is constrained by physical distance along cortical sheets and other synaptic wiring costs, we introduce locality masked RNNs (LM-RNNs) that use task-agnostic predetermined graphs with sparsity as low as 4%. We study LM-RNNs in a multitask learning setting relevant to cognitive systems neuroscience with a commonly used set of tasks, 20-Cog-tasks (Yang et al., 2019). We show through reductio ad absurdum that 20-Cog-tasks can be solved by a small pool of separated autapses that we can mechanistically analyze and understand. Thus, these tasks fall short of the goal of inducing complex recurrent dynamics and modular structure in RNNs. We next contribute a new cognitive multitask battery, Mod-Cog, consisting of up to 132 tasks that expands by about seven-fold the number of tasks and task complexity of 20-Cog-tasks. Importantly, while autapses can solve the simple 20-Cog-tasks, the expanded task set requires richer neural architectures and continuous attractor dynamics. On these tasks, we show that LM-RNNs with an optimal sparsity result in faster training and better data efficiency than fully connected networks.

Tau protein spreads through functionally connected neurons in Alzheimer's disease: a combined MEG/PET study

Recent studies on Alzheimer's disease (AD) suggest that tau proteins spread through the brain following neuronal connections. Several mechanisms could be involved in this process: spreading between brain regions that interact strongly (functional connectivity); through the pattern of anatomical connections (structural connectivity); or simple diffusion. Using magnetoencephalography (MEG), we investigated which spreading pathways influence tau protein spreading by modelling the tau propagation process using an epidemic spreading model. We compared the modelled tau depositions with 18F-flortaucipir PET binding potentials at several stages of the AD continuum. In this cross-sectional study, we analysed source-reconstructed MEG data and dynamic 100-min 18F-flortaucipir PET from 57 subjects positive for amyloid-β pathology [preclinical AD (n = 16), mild cognitive impairment (MCI) due to AD (n = 16) and AD dementia (n = 25)]. Cognitively healthy subjects without amyloid-β pathology were included as controls (n = 25). Tau propagation was modelled as an epidemic process (susceptible-infected model) on MEG-based functional networks [in alpha (8-13 Hz) and beta (13-30 Hz) bands], a structural or diffusion network, starting from the middle and inferior temporal lobe. The group-level network of the control group was used as input for the model to predict tau deposition in three stages of the AD continuum. To assess performance, model output was compared to the group-specific tau deposition patterns as measured with 18F-flortaucipir PET. We repeated the analysis by using networks of the preceding disease stage and/or using regions with most observed tau deposition during the preceding stage as seeds. In the preclinical AD stage, the functional networks predicted most of the modelled tau-PET binding potential, with best correlations between model and tau-PET [corrected amplitude envelope correlation (AEC-c) alpha C = 0.584; AEC-c beta C = 0.569], followed by the structural network (C = 0.451) and simple diffusion (C = 0.451). Prediction accuracy declined for the MCI and AD dementia stages, although the correlation between modelled tau and tau-PET binding remained highest for the functional networks (C = 0.384; C = 0.376). Replacing the control-network with the network from the preceding disease stage and/or alternative seeds improved prediction accuracy in MCI but not in the dementia stage. These results suggest that in addition to structural connections, functional connections play an important role in tau spread, and highlight that neuronal dynamics play a key role in promoting this pathological process. Aberrant neuronal communication patterns should be taken into account when identifying targets for future therapy. Our results also suggest that this process is more important in earlier disease stages (preclinical AD/MCI); possibly, in later stages, other processes may be influential.

Functional geometry of auditory cortical resting state networks derived from intracranial electrophysiology

Understanding central auditory processing critically depends on defining underlying auditory cortical networks and their relationship to the rest of the brain. We addressed these questions using resting state functional connectivity derived from human intracranial electroencephalography. Mapping recording sites into a low-dimensional space where proximity represents functional similarity revealed a hierarchical organization. At a fine scale, a group of auditory cortical regions excluded several higher-order auditory areas and segregated maximally from the prefrontal cortex. On mesoscale, the proximity of limbic structures to the auditory cortex suggested a limbic stream that parallels the classically described ventral and dorsal auditory processing streams. Identities of global hubs in anterior temporal and cingulate cortex depended on frequency band, consistent with diverse roles in semantic and cognitive processing. On a macroscale, observed hemispheric asymmetries were not specific for speech and language networks. This approach can be applied to multivariate brain data with respect to development, behavior, and disorders.

Local and long-distance organization of prefrontal cortex circuits in the marmoset brain

The prefrontal cortex (PFC) has dramatically expanded in primates, but its organization and interactions with other brain regions are only partially understood. We performed high-resolution connectomic mapping of the marmoset PFC and found two contrasting corticocortical and corticostriatal projection patterns: "patchy" projections that formed many columns of submillimeter scale in nearby and distant regions and "diffuse" projections that spread widely across the cortex and striatum. Parcellation-free analyses revealed representations of PFC gradients in these projections' local and global distribution patterns. We also demonstrated column-scale precision of reciprocal corticocortical connectivity, suggesting that PFC contains a mosaic of discrete columns. Diffuse projections showed considerable diversity in the laminar patterns of axonal spread. Altogether, these fine-grained analyses reveal important principles of local and long-distance PFC circuits in marmosets and provide insights into the functional organization of the primate brain.

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.

The Brain/MINDS Marmoset Connectivity Resource: An open-access platform for cellular-level tracing and tractography in the primate brain

The primate brain has unique anatomical characteristics, which translate into advanced cognitive, sensory, and motor abilities. Thus, it is important that we gain insight on its structure to provide a solid basis for models that will clarify function. Here, we report on the implementation and features of the Brain/MINDS Marmoset Connectivity Resource (BMCR), a new open-access platform that provides access to high-resolution anterograde neuronal tracer data in the marmoset brain, integrated to retrograde tracer and tractography data. Unlike other existing image explorers, the BMCR allows visualization of data from different individuals and modalities in a common reference space. This feature, allied to an unprecedented high resolution, enables analyses of features such as reciprocity, directionality, and spatial segregation of connections. The present release of the BMCR focuses on the prefrontal cortex (PFC), a uniquely developed region of the primate brain that is linked to advanced cognition, including the results of 52 anterograde and 164 retrograde tracer injections in the cortex of the marmoset. Moreover, the inclusion of tractography data from diffusion MRI allows systematic analyses of this noninvasive modality against gold-standard cellular connectivity data, enabling detection of false positives and negatives, which provide a basis for future development of tractography. This paper introduces the BMCR image preprocessing pipeline and resources, which include new tools for exploring and reviewing the data.

Critical scaling of whole-brain resting-state dynamics

Scale invariance is a characteristic of neural activity. How this property emerges from neural interactions remains a fundamental question. Here, we studied the relation between scale-invariant brain dynamics and structural connectivity by analyzing human resting-state (rs-) fMRI signals, together with diffusion MRI (dMRI) connectivity and its approximation as an exponentially decaying function of the distance between brain regions. We analyzed the rs-fMRI dynamics using functional connectivity and a recently proposed phenomenological renormalization group (PRG) method that tracks the change of collective activity after successive coarse-graining at different scales. We found that brain dynamics display power-law correlations and power-law scaling as a function of PRG coarse-graining based on functional or structural connectivity. Moreover, we modeled the brain activity using a network of spins interacting through large-scale connectivity and presenting a phase transition between ordered and disordered phases. Within this simple model, we found that the observed scaling features were likely to emerge from critical dynamics and connections exponentially decaying with distance. In conclusion, our study tests the PRG method using large-scale brain activity and theoretical models and suggests that scaling of rs-fMRI activity relates to criticality.

Geometric constraints on human brain function

The anatomy of the brain necessarily constrains its function, but precisely how remains unclear. The classical and dominant paradigm in neuroscience is that neuronal dynamics are driven by interactions between discrete, functionally specialized cell populations connected by a complex array of axonal fibres1-3. However, predictions from neural field theory, an established mathematical framework for modelling large-scale brain activity4-6, suggest that the geometry of the brain may represent a more fundamental constraint on dynamics than complex interregional connectivity7,8. Here, we confirm these theoretical predictions by analysing human magnetic resonance imaging data acquired under spontaneous and diverse task-evoked conditions. Specifically, we show that cortical and subcortical activity can be parsimoniously understood as resulting from excitations of fundamental, resonant modes of the brain's geometry (that is, its shape) rather than from modes of complex interregional connectivity, as classically assumed. We then use these geometric modes to show that task-evoked activations across over 10,000 brain maps are not confined to focal areas, as widely believed, but instead excite brain-wide modes with wavelengths spanning over 60 mm. Finally, we confirm predictions that the close link between geometry and function is explained by a dominant role for wave-like activity, showing that wave dynamics can reproduce numerous canonical spatiotemporal properties of spontaneous and evoked recordings. Our findings challenge prevailing views and identify a previously underappreciated role of geometry in shaping function, as predicted by a unifying and physically principled model of brain-wide dynamics.

An anatomical and connectivity atlas of the marmoset cerebellum

The cerebellum is essential for motor control and cognitive functioning, engaging in bidirectional communication with the cerebral cortex. The common marmoset, a small non-human primate, offers unique advantages for studying cerebello-cerebral circuits. However, the marmoset cerebellum is not well described in published resources. In this study, we present a comprehensive atlas of the marmoset cerebellum comprising (1) fine-detailed anatomical atlases and surface-analysis tools of the cerebellar cortex based on ultra-high-resolution ex vivo MRI, (2) functional connectivity and gradient patterns of the cerebellar cortex revealed by awake resting-state fMRI, and (3) structural-connectivity mapping of cerebellar nuclei using high-resolution diffusion MRI tractography. The atlas elucidates the anatomical details of the marmoset cerebellum, reveals distinct gradient patterns of intra-cerebellar and cerebello-cerebral functional connectivity, and maps the topological relationship of cerebellar nuclei in cerebello-cerebral circuits. As version 5 of the Marmoset Brain Mapping project, this atlas is publicly available at https://marmosetbrainmapping.org/MBMv5.html.

Toward Best Practices for Imaging Transcriptomics of the Human Brain

Modern brainwide transcriptional atlases provide unprecedented opportunities for investigating the molecular correlates of brain organization, as quantified using noninvasive neuroimaging. However, integrating neuroimaging data with transcriptomic measures is not straightforward, and careful consideration is required to make valid inferences. In this article, we review recent work exploring how various methodological choices affect 3 main phases of imaging transcriptomic analyses, including 1) processing of transcriptional atlas data; 2) relating transcriptional measures to independently derived neuroimaging phenotypes; and 3) evaluating the functional implications of identified associations through gene enrichment analyses. Our aim is to facilitate the development of standardized and reproducible approaches for this rapidly growing field. We identify sources of methodological variability, key choices that can affect findings, and considerations for mitigating false positive and/or spurious results. Finally, we provide an overview of freely available open-source toolboxes implementing current best-practice procedures across all 3 analysis phases.

Neuronal vulnerability to brain aging and neurodegeneration in cognitively impaired marmoset monkeys (Callithrix jacchus)

The investigation of neurobiological and neuropathological changes that affect synaptic integrity and function with aging is key to understanding why the aging brain is vulnerable to Alzheimer's disease. We investigated the cellular characteristics in the cerebral cortex of behaviorally characterized marmosets, based on their trajectories of cognitive learning as they transitioned to old age. We found increased astrogliosis, increased phagocytic activity of microglial cells and differences in resting and reactive microglial cell phenotypes in cognitively impaired compared to nonimpaired marmosets. Differences in amyloid beta deposition were not related to cognitive trajectory. However, we found age-related changes in density and morphology of dendritic spines in pyramidal neurons of layer 3 in the dorsolateral prefrontal cortex and the CA1 field of the hippocampus between cohorts. Overall, our data suggest that an accelerated aging process, accompanied by neurodegeneration, that takes place in cognitively impaired aged marmosets and affects the plasticity of dendritic spines in cortical areas involved in cognition and points to mechanisms of neuronal vulnerability to aging.

Marmoset core visual object recognition behavior is comparable to that of macaques and humans

Among the smallest simian primates, the common marmoset offers promise as an experimentally tractable primate model for neuroscience with translational potential to humans. However, given its exceedingly small brain and body, the gap in perceptual and cognitive abilities between marmosets and humans requires study. Here, we performed a comparison of marmoset behavior to that of three other species in the domain of high-level vision. We first found that marmosets outperformed rats - a marmoset-sized rodent - on a simple recognition task, with marmosets robustly recognizing objects across views. On a more challenging invariant object recognition task used previously in humans, marmosets also achieved high performance. Notably, across hundreds of images, marmosets' image-by-image behavior was highly similar to that of humans - nearly as human-like as macaque behavior. Thus, core aspects of visual perception are conserved across monkeys and humans, and marmosets present salient behavioral advantages over other small model organisms for visual neuroscience.

Connectome topology of mammalian brains and its relationship to taxonomy and phylogeny

Network models of anatomical connections allow for the extraction of quantitative features describing brain organization, and their comparison across brains from different species. Such comparisons can inform our understanding of between-species differences in brain architecture and can be compared to existing taxonomies and phylogenies. Here we performed a quantitative comparative analysis using the MaMI database (Tel Aviv University), a collection of brain networks reconstructed from ex vivo diffusion MRI spanning 125 species and 12 taxonomic orders or superorders. We used a broad range of metrics to measure between-mammal distances and compare these estimates to the separation of species as derived from taxonomy and phylogeny. We found that within-taxonomy order network distances are significantly closer than between-taxonomy network distances, and this relation holds for several measures of network distance. Furthermore, to estimate the evolutionary divergence between species, we obtained phylogenetic distances across 10,000 plausible phylogenetic trees. The anatomical network distances were rank-correlated with phylogenetic distances 10,000 times, creating a distribution of coefficients that demonstrate significantly positive correlations between network and phylogenetic distances. Collectively, these analyses demonstrate species-level organization across scales and informational sources: we relate brain networks distances, derived from MRI, with evolutionary distances, derived from genotyping data.

Neuronal density and expression of calcium-binding proteins across the layers of the superior colliculus in the common marmoset (Callithrix jacchus)

The superior colliculus (SC) is a layered midbrain structure with functions that include polysensory and sensorimotor integration. Here, we describe the distribution of different immunohistochemically identified classes of neurons in the SC of adult marmoset monkeys (Callithrix jacchus). Neuronal nuclei (NeuN) staining was used to determine the overall neuronal density in the different SC layers. In addition, we studied the distribution of neurons expressing different calcium-binding proteins (calbindin [CB], parvalbumin [PV] and calretinin [CR]). Our results indicate that neuronal density in the SC decreases from superficial to deep layers. Although the neuronal density within the same layer varies little across the mediolateral axis, it tends to be lower at rostral levels, compared to caudal levels. Cells expressing different calcium-binding proteins display differential gradients of density according to depth. Both CB- and CR-expressing neurons show markedly higher densities in the stratum griseum superficiale (SGS), compared to the stratum opticum and intermediate and deep layers. However, CR-expressing neurons are twice as common as CB-expressing neurons outside the SGS. The distribution of PV-expressing cells follows a shallow density gradient from superficial to deep layers. When normalized relative to total neuronal density, the proportion of CR-expressing neurons increases between the superficial and intermediate layers, whereas that of CB-expressing neurons declines toward the deep layers. The proportion of PV-expressing neurons remains constant across layers. Our data provide layer-specific and accurate estimates of neuronal density, which may be important for the generation of biophysical models of how the primate SC transforms sensory inputs into motor signals.

An integrated resource for functional and structural connectivity of the marmoset brain

Comprehensive integration of structural and functional connectivity data is required to model brain functions accurately. While resources for studying the structural connectivity of non-human primate brains already exist, their integration with functional connectivity data has remained unavailable. Here we present a comprehensive resource that integrates the most extensive awake marmoset resting-state fMRI data available to date (39 marmoset monkeys, 710 runs, 12117 mins) with previously published cellular-level neuronal tracing data (52 marmoset monkeys, 143 injections) and multi-resolution diffusion MRI datasets. The combination of these data allowed us to (1) map the fine-detailed functional brain networks and cortical parcellations, (2) develop a deep-learning-based parcellation generator that preserves the topographical organization of functional connectivity and reflects individual variabilities, and (3) investigate the structural basis underlying functional connectivity by computational modeling. This resource will enable modeling structure-function relationships and facilitate future comparative and translational studies of primate brains.

Theory of the Multiregional Neocortex: Large-Scale Neural Dynamics and Distributed Cognition

The neocortex is a complex neurobiological system with many interacting regions. How these regions work together to subserve flexible behavior and cognition has become increasingly amenable to rigorous research. Here, I review recent experimental and theoretical work on the modus operandi of a multiregional cortex. These studies revealed several general principles for the neocortical interareal connectivity, low-dimensional macroscopic gradients of biological properties across cortical areas, and a hierarchy of timescales for information processing. Theoretical work suggests testable predictions regarding differential excitation and inhibition along feedforward and feedback pathways in the cortical hierarchy. Furthermore, modeling of distributed working memory and simple decision-making has given rise to a novel mathematical concept, dubbed bifurcation in space, that potentially explains how different cortical areas, with a canonical circuit organization but gradients of biological heterogeneities, are able to subserve their respective (e.g., sensory coding versus executive control) functions in a modularly organized brain.

Epidemic models characterize seizure propagation and the effects of epilepsy surgery in individualized brain networks based on MEG and invasive EEG recordings

Epilepsy surgery is the treatment of choice for drug-resistant epilepsy patients. However, seizure-freedom is currently achieved in only 2/3 of the patients after surgery. In this study we have developed an individualized computational model based on MEG brain networks to explore seizure propagation and the efficacy of different virtual resections. Eventually, the goal is to obtain individualized models to optimize resection strategy and outcome. We have modelled seizure propagation as an epidemic process using the susceptible-infected (SI) model on individual brain networks derived from presurgical MEG. We included 10 patients who had received epilepsy surgery and for whom the surgery outcome at least one year after surgery was known. The model parameters were tuned in in order to reproduce the patient-specific seizure propagation patterns as recorded with invasive EEG. We defined a personalized search algorithm that combined structural and dynamical information to find resections that maximally decreased seizure propagation for a given resection size. The optimal resection for each patient was defined as the smallest resection leading to at least a 90% reduction in seizure propagation. The individualized model reproduced the basic aspects of seizure propagation for 9 out of 10 patients when using the resection area as the origin of epidemic spreading, and for 10 out of 10 patients with an alternative definition of the seed region. We found that, for 7 patients, the optimal resection was smaller than the resection area, and for 4 patients we also found that a resection smaller than the resection area could lead to a 100% decrease in propagation. Moreover, for two cases these alternative resections included nodes outside the resection area. Epidemic spreading models fitted with patient specific data can capture the fundamental aspects of clinically observed seizure propagation, and can be used to test virtual resections in silico. Combined with optimization algorithms, smaller or alternative resection strategies, that are individually targeted for each patient, can be determined with the ultimate goal to improve surgery outcome. MEG-based networks can provide a good approximation of structural connectivity for computational models of seizure propagation, and facilitate their clinical use.

Directed functional and structural connectivity in a large-scale model for the mouse cortex

Inferring the structural connectivity from electrophysiological measurements is a fundamental challenge in systems neuroscience. Directed functional connectivity measures, such as the generalized partial directed coherence (GPDC), provide estimates of the causal influence between areas. However, the relation between causality estimates and structural connectivity is still not clear. We analyzed this problem by evaluating the effectiveness of GPDC to estimate the connectivity of a ground-truth, data-constrained computational model of a large-scale network model of the mouse cortex. The model contains 19 cortical areas composed of spiking neurons, with areas connected by long-range projections with weights obtained from a tract-tracing cortical connectome. We show that GPDC values provide a reasonable estimate of structural connectivity, with an average Pearson correlation over simulations of 0.74. Moreover, even in a typical electrophysiological recording scenario containing five areas, the mean correlation was above 0.6. These results suggest that it may be possible to empirically estimate structural connectivity from functional connectivity even when detailed whole-brain recordings are not achievable.

We analyzed the relationship between structural and directed functional connectivity by evaluating the effectiveness of generalized partial directed coherence (GPDC) to estimate the connectivity of a ground-truth, data-constrained computational model of a large-scale network model of the mouse cortex. We show that GPDC values provide a reasonable estimate of structural connectivity even in a typical electrophysiological recording scenario containing few areas. These results suggest that it may be possible to empirically estimate structural connectivity from functional connectivity even when detailed whole-brain recordings are not achievable.

Brain rhythms define distinct interaction networks with differential dependence on anatomy

Cognitive functions are subserved by rhythmic neuronal synchronization across widely distributed brain areas. In 105 area pairs, we investigated functional connectivity (FC) through coherence, power correlation, and Granger causality (GC) in the theta, beta, high-beta, and gamma rhythms. Between rhythms, spatial FC patterns were largely independent. Thus, the rhythms defined distinct interaction networks. Importantly, networks of coherence and GC were not explained by the spatial distributions of the strengths of the rhythms. Those networks, particularly the GC networks, contained clear modules, with typically one dominant rhythm per module. To understand how this distinctiveness and modularity arises on a common anatomical backbone, we correlated, across 91 area pairs, the metrics of functional interaction with those of anatomical projection strength. Anatomy was primarily related to coherence and GC, with the largest effect sizes for GC. The correlation differed markedly between rhythms, being less pronounced for the beta and strongest for the gamma rhythm.

Anatomical and functional connectomes underlying hierarchical visual processing in mouse visual system

Over the last 10 years, there has been a surge in interest in the rodent visual system resulting from the discovery of visual processing functions shared with primates V1, and of a complex anatomical structure in the extrastriate visual cortex. This surprisingly intricate visual system was elucidated by recent investigations using rapidly growing genetic tools primarily available in the mouse. Here, we examine the structural and functional connections of visual areas that have been identified in mice mostly during the past decade, and the impact of these findings on our understanding of brain functions associated with vision. Special attention is paid to structure-function relationships arising from the hierarchical organization, which is a prominent feature of the primate visual system. Recent evidence supports the existence of a hierarchical organization in rodents that contains levels that are poorly resolved relative to those observed in primates. This shallowness of the hierarchy indicates that the mouse visual system incorporates abundant non-hierarchical processing. Thus, the mouse visual system provides a unique opportunity to study non-hierarchical processing and its relation to hierarchical processing.

Co-registration of Imaging Modalities (MRI, CT and PET) to Perform Frameless Stereotaxic Robotic Injections in the Common Marmoset

The common marmoset has emerged as a popular model in neuroscience research, in part due to its reproductive efficiency, genetic and neuroanatomical similarities to humans and the successful generation of transgenic lines. Stereotaxic procedures in marmosets are guided by 2D stereotaxic atlases, which are constructed with a limited number of animals and fail to account for inter-individual variability in skull and brain size. Here, we developed a frameless imaging-guided stereotaxic system that improves upon traditional approaches by using subject-specific registration of computed tomography (CT), magnetic resonance imaging (MRI) and positron emission tomography (PET) data to identify a surgical target, namely the putamen, in two marmosets. The skull surface was laser-scanned to create a point cloud that was registered to the 3D reconstruction of the skull from CT. Reconstruction of the skull, as well as of the brain from MR images, was crucial for surgical planning. Localisation and injection into the putamen was done using a 6-axis robotic arm controlled by a surgical navigation software (Brainsight™). Integration of subject-specific registration and frameless stereotaxic navigation allowed target localisation specific to each animal. Injection of alpha-synuclein fibrils into the putamen triggered progressive neurodegeneration of the nigro-striatal system, a key feature of Parkinson's disease. Four months post-surgery, a PET scan found evidence of nigro-striatal denervation, supporting accurate targeting of the putamen during co-registration and subsequent surgery. Our results suggest that this approach, coupled with frameless stereotaxic neuronavigation, is accurate in localising surgical targets and can be used to assess endpoints for longitudinal studies.

Interpreting the Entire Connectivity of Individual Neurons in Micropatterned Neural Culture With an Integrated Connectome Analyzer of a Neuronal Network (iCANN)

The function of a neural circuit can be determined by the following: (1) characteristics of individual neurons composing the circuit, (2) their distinct connection structure, and (3) their neural circuit activity. However, prior research on correlations between these three factors revealed many limitations. In particular, profiling and modeling of the connectivity of complex neural circuits at the cellular level are highly challenging. To reduce the burden of the analysis, we suggest a new approach with simplification of the neural connection in an array of honeycomb patterns on 2D, using a microcontact printing technique. Through a series of guided neuronal growths in defined honeycomb patterns, a simplified neuronal circuit was achieved. Our approach allowed us to obtain the whole network connectivity at cellular resolution using a combination of stochastic multicolor labeling via viral transfection. Therefore, we were able to identify several types of hub neurons with distinct connectivity features. We also compared the structural differences between different circuits using three-node motif analysis. This new model system, iCANN, is the first experimental model of neural computation at the cellular level, providing neuronal circuit structures for the study of the relationship between anatomical structure and function of the neuronal network.

Morphometric analysis of astrocytes in vocal production circuits of common marmoset (Callithrix jacchus)

Astrocytes, the star-shaped glial cells, are the most abundant non-neuronal cell population in the central nervous system. They play a key role in modulating activities of neural networks, including those involved in complex motor behaviors. Common marmosets (Callithrix jacchus), the most vocal non-human primate (NHP), have been used to study the physiology of vocalization and social vocal production. However, the neural circuitry involved in vocal production is not fully understood. In addition, even less is known about the involvement of astrocytes in this circuit. To understand the role, that astrocytes may play in the complex behavior of vocalization, the initial step may be to study their structural properties in the cortical and subcortical regions that are known to be involved in vocalization. Here, in the common marmoset, we identify all astrocytic subtypes seen in other primate's brains, including intralaminar astrocytes. In addition, we reveal detailed structural characteristics of astrocytes and perform morphometric analysis of astrocytes residing in the cortex and midbrain regions that are associated with vocal production. We found that cortical astrocytes in these regions illustrate a higher level of complexity when compared to those in the midbrain. We hypothesize that this complexity that is expressed in cortical astrocytes may reflect their functions to meet the metabolic/structural needs of these regions.