The highways and byways of the brain.

Brain functions rely on the communication network formed by axonal fibers. However, the number of axons connecting different brain regions is unknown. A study in PLoS Biology addresses this question and finds that most areas of the human cerebral cortex are linked by an astoundingly small number of fibers.

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

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

Structural basis of envelope and phase intrinsic coupling modes in the cerebral cortex.

Intrinsic coupling modes (ICMs) can be observed in ongoing brain activity at multiple spatial and temporal scales. Two families of ICMs can be distinguished: phase and envelope ICMs. The principles that shape these ICMs remain partly elusive, in particular their relation to the underlying brain structure. Here we explored structure-function relationships in the ferret brain between ICMs quantified from ongoing brain activity recorded with chronically implanted micro-ECoG arrays and structural connectivity (SC) obtained from high-resolution diffusion MRI tractography. Large-scale computational models were used to explore the ability to predict both types of ICMs. Importantly, all investigations were conducted with ICM measures that are sensitive or insensitive to volume conduction effects. The results show that both types of ICMs are significantly related to SC, except for phase ICMs when using measures removing zero-lag coupling. The correlation between SC and ICMs increases with increasing frequency which is accompanied by reduced delays. Computational models produced results that were highly dependent on the specific parameter settings. The most consistent predictions were derived from measures solely based on SC. Overall, the results demonstrate that patterns of cortical functional coupling as reflected in both phase and envelope ICMs are both related, albeit to different degrees, to the underlying structural connectivity in the cerebral cortex.

Systematic perturbation of an artificial neural network: A step towards quantifying causal contributions in the brain.

Lesion inference analysis is a fundamental approach for characterizing the causal contributions of neural elements to brain function. This approach has gained new prominence through the arrival of modern perturbation techniques with unprecedented levels of spatiotemporal precision. While inferences drawn from brain perturbations are conceptually powerful, they face methodological difficulties. Particularly, they are challenged to disentangle the true causal contributions of the involved elements, since often functions arise from coalitions of distributed, interacting elements, and localized perturbations have unknown global consequences. To elucidate these limitations, we systematically and exhaustively lesioned a small artificial neural network (ANN) playing a classic arcade game. We determined the functional contributions of all nodes and links, contrasting results from sequential single-element perturbations with simultaneous perturbations of multiple elements. We found that lesioning individual elements, one at a time, produced biased results. By contrast, multi-site lesion analysis captured crucial details that were missed by single-site lesions. We conclude that even small and seemingly simple ANNs show surprising complexity that needs to be addressed by multi-lesioning for a coherent causal characterization.

Brain connectivity meets reservoir computing.

The connectivity of Artificial Neural Networks (ANNs) is different from the one observed in Biological Neural Networks (BNNs). Can the wiring of actual brains help improve ANNs architectures? Can we learn from ANNs about what network features support computation in the brain when solving a task? At a meso/macro-scale level of the connectivity, ANNs' architectures are carefully engineered and such those design decisions have crucial importance in many recent performance improvements. On the other hand, BNNs exhibit complex emergent connectivity patterns at all scales. At the individual level, BNNs connectivity results from brain development and plasticity processes, while at the species level, adaptive reconfigurations during evolution also play a major role shaping connectivity. Ubiquitous features of brain connectivity have been identified in recent years, but their role in the brain's ability to perform concrete computations remains poorly understood. Computational neuroscience studies reveal the influence of specific brain connectivity features only on abstract dynamical properties, although the implications of real brain networks topologies on machine learning or cognitive tasks have been barely explored. Here we present a cross-species study with a hybrid approach integrating real brain connectomes and Bio-Echo State Networks, which we use to solve concrete memory tasks, allowing us to probe the potential computational implications of real brain connectivity patterns on task solving. We find results consistent across species and tasks, showing that biologically inspired networks perform as well as classical echo state networks, provided a minimum level of randomness and diversity of connections is allowed. We also present a framework, bio2art, to map and scale up real connectomes that can be integrated into recurrent ANNs. This approach also allows us to show the crucial importance of the diversity of interareal connectivity patterns, stressing the importance of stochastic processes determining neural networks connectivity in general.

Reaction-diffusion models in weighted and directed connectomes.

Connectomes represent comprehensive descriptions of neural connections in a nervous system to better understand and model central brain function and peripheral processing of afferent and efferent neural signals. Connectomes can be considered as a distinctive and necessary structural component alongside glial, vascular, neurochemical, and metabolic networks of the nervous systems of higher organisms that are required for the control of body functions and interaction with the environment. They are carriers of functional phenomena such as planning behavior and cognition, which are based on the processing of highly dynamic neural signaling patterns. In this study, we examine more detailed connectomes with edge weighting and orientation properties, in which reciprocal neuronal connections are also considered. Diffusion processes are a further necessary condition for generating dynamic bioelectric patterns in connectomes. Based on our precise connectome data, we investigate different diffusion-reaction models to study the propagation of dynamic concentration patterns in control and lesioned connectomes. Therefore, differential equations for modeling diffusion were combined with well-known reaction terms to allow the use of connection weights, connectivity orientation and spatial distances. Three reaction-diffusion systems Gray-Scott, Gierer-Meinhardt and Mimura-Murray were investigated. For this purpose, implicit solvers were implemented in a numerically stable reaction-diffusion system within the framework of neuroVIISAS. The implemented reaction-diffusion systems were applied to a subconnectome which shapes the mechanosensitive pathway that is strongly affected in the multiple sclerosis demyelination disease. It was found that demyelination modeling by connectivity weight modulation changes the oscillations of the target region, i.e. the primary somatosensory cortex, of the mechanosensitive pathway. In conclusion, a new application of reaction-diffusion systems to weighted and directed connectomes has been realized. Because the implementation was realized in the neuroVIISAS framework many possibilities for the study of dynamic reaction-diffusion processes in empirical connectomes as well as specific randomized network models are available now.

Brain simulation as a cloud service: The Virtual Brain on EBRAINS.

The Virtual Brain (TVB) is now available as open-source services on the cloud research platform EBRAINS (ebrains.eu). It offers software for constructing, simulating and analysing brain network models including the TVB simulator; magnetic resonance imaging (MRI) processing pipelines to extract structural and functional brain networks; combined simulation of large-scale brain networks with small-scale spiking networks; automatic conversion of user-specified model equations into fast simulation code; simulation-ready brain models of patients and healthy volunteers; Bayesian parameter optimization in epilepsy patient models; data and software for mouse brain simulation; and extensive educational material. TVB cloud services facilitate reproducible online collaboration and discovery of data assets, models, and software embedded in scalable and secure workflows, a precondition for research on large cohort data sets, better generalizability, and clinical translation.

A blueprint of mammalian cortical connectomes.

The cerebral cortex of mammals exhibits intricate interareal wiring. Moreover, mammalian cortices differ vastly in size, cytological composition, and phylogenetic distance. Given such complexity and pronounced species differences, it is a considerable challenge to decipher organizational principles of mammalian connectomes. Here, we demonstrate species-specific and species-general unifying principles linking the physical, cytological, and connectional dimensions of architecture in the mouse, cat, marmoset, and macaque monkey. The existence of connections is related to the cytology of cortical areas, in addition to the role of physical distance, but this relation is attenuated in mice and marmoset monkeys. The cytoarchitectonic cortical gradients, and not the rostrocaudal axis of the cortex, are closely linked to the laminar origin of connections, a principle that allows the extrapolation of this connectional feature to humans. Lastly, a network core, with a central role under different modes of network communication, characterizes all cortical connectomes. We observe a displacement of the network core in mammals, with a shift of the core of cats and macaque monkeys toward the less neuronally dense areas of the cerebral cortex. This displacement has functional ramifications but also entails a potential increased degree of vulnerability to pathology. In sum, our results sketch out a blueprint of mammalian connectomes consisting of species-specific and species-general links between the connectional, physical, and cytological dimensions of the cerebral cortex, possibly reflecting variations and persistence of evolutionarily conserved mechanisms and cellular phenomena. Our framework elucidates organizational principles that encompass but also extend beyond the wiring economy principle imposed by the physical embedding of the cerebral cortex.

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.

Systematic modelling of the development of laminar projection origins in the cerebral cortex: Interactions of spatio-temporal patterns of neurogenesis and cellular heterogeneity.

The architectonic type principle conceptualizes structural connections between brain areas in terms of the relative architectonic differentiation of connected areas. It has previously been shown that spatio-temporal interactions between the time and place of neurogenesis could underlie multiple features of empirical mammalian connectomes, such as projection existence and the distribution of projection strengths. However, so far no mechanistic explanation for the emergence of typically observed laminar patterns of projection origins and terminations has been tested. Here, we expand an in silico model of the developing cortical sheet to explore which factors could potentially constrain the development of laminar projection patterns. We show that manipulations which rely solely on spatio-temporal interactions, namely the relative density of laminar compartments, a delay in the neurogenesis of infragranular layers relative to layer 1, and a delay in the neurogenesis of supragranular layers relative to infragranular layers, do not result in the striking correlation between supragranular contribution to projections and the relative differentiation of areas that is typically observed in the mammalian cortex. In contrast, we find that if we introduce systematic variation in cell-intrinsic properties, coupling them with architectonic differentiation, the resulting laminar projection patterns closely mirror the empirically observed patterns. We also find that the spatio-temporal interactions posited to occur during neurogenesis are necessary for the formation of the characteristic laminar patterns. Hence, our results indicate that the specification of the laminar patterns of projection origins may result from systematic variation in a number of cell-intrinsic properties, superimposed on the previously identified spatio-temporal interactions which are sufficient for the emergence of the architectonic type principle on the level of inter-areal connectivity in silico.

Revisiting 'brain modes' in a new computational era: approaches for the characterization of brain-behavioural associations.

The study of brain-function relationships is undergoing a conceptual and methodological transformation due to the emergence of network neuroscience and the development of multivariate methods for lesion-deficit inferences. Anticipating this process, in 1998 Godefroy and co-workers conceptualized the potential of four elementary typologies of brain-behaviour relationships named 'brain modes' (unicity, equivalence, association, summation) as building blocks able to describe the association between intact or lesioned brain regions and cognitive processes or neurological deficits. In the light of new multivariate lesion inference and network approaches, we critically revisit and update the original theoretical notion of brain modes, and provide real-life clinical examples that support their existence. To improve the characterization of elementary units of brain-behavioural relationships further, we extend such conceptualization with a fifth brain mode (mutual inhibition/masking summation). We critically assess the ability of these five brain modes to account for any type of brain-function relationship, and discuss past versus future contributions in redefining the anatomical basis of human cognition. We also address the potential of brain modes for predicting the behavioural consequences of lesions and their future role in the design of cognitive neurorehabilitation therapies.

Game theoretical mapping of white matter contributions to visuospatial attention in stroke patients with hemineglect.

White matter bundles linking gray matter nodes are key anatomical players to fully characterize associations between brain systems and cognitive functions. Here we used a multivariate lesion inference approach grounded in coalitional game theory (multiperturbation Shapley value analysis, MSA) to infer causal contributions of white matter bundles to visuospatial orienting of attention. Our work is based on the characterization of the lesion patterns of 25 right hemisphere stroke patients and the causal analysis of their impact on three neuropsychological tasks: line bisection, letter cancellation, and bells cancellation. We report that, out of the 11 white matter bundles included in our MSA coalitions, the optic radiations, the inferior fronto-occipital fasciculus and the anterior cingulum were the only tracts to display task-invariant contributions (positive, positive, and negative, respectively) to the tasks. We also report task-dependent influences for the branches of the superior longitudinal fasciculus and the posterior cingulum. By extending prior findings to white matter tracts linking key gray matter nodes, we further characterize from a network perspective the anatomical basis of visual and attentional orienting processes. The knowledge about interactions patterns mediated by white matter tracts linking cortical nodes of attention orienting networks, consolidated by further studies, may help develop and customize brain stimulation approaches for the rehabilitation of visuospatial neglect.

Neuron density fundamentally relates to architecture and connectivity of the primate cerebral cortex.

Studies of structural brain connectivity have revealed many intriguing features of complex cortical networks. To advance integrative theories of cortical organization, an understanding is required of how connectivity interrelates with other aspects of brain structure. Recent studies have suggested that interareal connectivity may be related to a variety of macroscopic as well as microscopic architectonic features of cortical areas. However, it is unclear how these features are inter-dependent and which of them most strongly and fundamentally relate to structural corticocortical connectivity. Here, we systematically investigated the relation of a range of microscopic and macroscopic architectonic features of cortical organization, namely layer III pyramidal cell soma cross section, dendritic synapse count, dendritic synapse density and dendritic tree size as well as area neuron density, to multiple properties of cortical connectivity, using a comprehensive, up-to-date structural connectome of the primate brain. Importantly, relationships were investigated by multi-variate analyses to account for the interrelations of features. Of all considered factors, the classical architectonic parameter of neuron density most strongly and consistently related to essential features of cortical connectivity (existence and laminar patterns of projections, area degree), and in conjoint analyses largely abolished effects of cellular morphological features. These results confirm neuron density as a central architectonic indicator of the primate cerebral cortex that is closely related to essential aspects of brain connectivity and is also highly indicative of further features of the architectonic organization of cortical areas, such as the considered cellular morphological measures. Our findings integrate several aspects of cortical micro- and macroscopic organization, with implications for cortical development and function.

Comprehensive computational modelling of the development of mammalian cortical connectivity underlying an architectonic type principle.

The architectonic type principle relates patterns of cortico-cortical connectivity to the relative architectonic differentiation of cortical regions. One mechanism through which the observed close relation between cortical architecture and connectivity may be established is the joint development of cortical areas and their connections in developmental time windows. Here, we describe a theoretical exploration of the possible mechanistic underpinnings of the architectonic type principle, by performing systematic computational simulations of cortical development. The main component of our in silico model was a developing two-dimensional cortical sheet, which was gradually populated by neurons that formed cortico-cortical connections. To assess different explanatory mechanisms, we varied the spatiotemporal trajectory of the simulated neurogenesis. By keeping the rules governing axon outgrowth and connection formation constant across all variants of simulated development, we were able to create model variants which differed exclusively by the specifics of when and where neurons were generated. Thus, all differences in the resulting connectivity were due to the variations in spatiotemporal growth trajectories. Our results demonstrated that a prescribed targeting of interareal connection sites was not necessary for obtaining a realistic replication of the experimentally observed relation between connection patterns and architectonic differentiation. Instead, we found that spatiotemporal interactions within the forming cortical sheet were sufficient if a small number of empirically well-grounded assumptions were met, namely planar, expansive growth of the cortical sheet around two points of origin as neurogenesis progressed, stronger architectonic differentiation of cortical areas for later neurogenetic time windows, and stochastic connection formation. Thus, our study highlights a potential mechanism of how relative architectonic differentiation and cortical connectivity become linked during development. We successfully predicted connectivity in two species, cat and macaque, from simulated cortico-cortical connection networks, which further underscored the general applicability of mechanisms through which the architectonic type principle can explain cortical connectivity in terms of the relative architectonic differentiation of cortical regions.

Toward a theory of coactivation patterns in excitable neural networks.

The relationship between the structural connectivity (SC) and functional connectivity (FC) of neural systems is of central importance in brain network science. It is an open question, however, how the SC-FC relationship depends on specific topological features of brain networks or the models used for describing neural dynamics. Using a basic but general model of discrete excitable units that follow a susceptible-excited-refractory activity cycle (SER model), we here analyze how the network activity patterns underlying functional connectivity are shaped by the characteristic topological features of the network. We develop an analytical framework for describing the contribution of essential topological elements, such as common inputs and pacemakers, to the coactivation of nodes, and demonstrate the validity of the approach by comparison of the analytical predictions with numerical simulations of various exemplar networks. The present analytic framework may serve as an initial step for the mechanistic understanding of the contributions of brain network topology to brain dynamics.

Intrinsic Functional Connectivity Resembles Cortical Architecture at Various Levels of Isoflurane Anesthesia.

Cortical single neuron activity and local field potential patterns change at different depths of general anesthesia. Here, we investigate the associated network level changes of functional connectivity. We recorded ongoing electrocorticographic (ECoG) activity from temporo-parieto-occipital cortex of 6 ferrets at various levels of isoflurane/nitrous oxide anesthesia and determined functional connectivity by computing amplitude envelope correlations. Through hierarchical clustering, we derived typical connectivity patterns corresponding to light, intermediate and deep anesthesia. Generally, amplitude correlation strength increased strongly with depth of anesthesia across all cortical areas and frequency bands. This was accompanied, at the deepest level, by the emergence of burst-suppression activity in the ECoG signal and a change of the spectrum of the amplitude envelope. Normalization of functional connectivity to the distribution of correlation coefficients showed that the topographical patterns remained similar across depths of anesthesia, reflecting the functional association of the underlying cortical areas. Thus, while strength and temporal properties of amplitude co-modulation vary depending on the activity of local neural circuits, their network-level interaction pattern is presumably most strongly determined by the underlying structural connectivity.

Features of spatial and functional segregation and integration of the primate connectome revealed by trade-off between wiring cost and efficiency.

The primate connectome, possessing a characteristic global topology and specific regional connectivity profiles, is well organized to support both segregated and integrated brain function. However, the organization mechanisms shaping the characteristic connectivity and its relationship to functional requirements remain unclear. The primate brain connectome is shaped by metabolic economy as well as functional values. Here, we explored the influence of two competing factors and additional advanced functional requirements on the primate connectome employing an optimal trade-off model between neural wiring cost and the representative functional requirement of processing efficiency. Moreover, we compared this model with a generative model combining spatial distance and topological similarity, with the objective of statistically reproducing multiple topological features of the network. The primate connectome indeed displays a cost-efficiency trade-off and that up to 67% of the connections were recovered by optimal combination of the two basic factors of wiring economy and processing efficiency, clearly higher than the proportion of connections (56%) explained by the generative model. While not explicitly aimed for, the trade-off model captured several key topological features of the real connectome as the generative model, yet better explained the connectivity of most regions. The majority of the remaining 33% of connections unexplained by the best trade-off model were long-distance links, which are concentrated on few cortical areas, termed long-distance connectors (LDCs). The LDCs are mainly non-hubs, but form a densely connected group overlapping on spatially segregated functional modalities. LDCs are crucial for both functional segregation and integration across different scales. These organization features revealed by the optimization analysis provide evidence that the demands of advanced functional segregation and integration among spatially distributed regions may play a significant role in shaping the cortical connectome, in addition to the basic cost-efficiency trade-off. These findings also shed light on inherent vulnerabilities of brain networks in diseases.

Game theoretical mapping of causal interactions underlying visuo-spatial attention in the human brain based on stroke lesions.

Anatomical studies conducted in neurological conditions have developed our understanding of the causal relationships between brain lesions and their clinical consequences. The analysis of lesion patterns extended across brain networks has been particularly useful in offering new insights on brain-behavior relationships. Here we applied multiperturbation Shapley value Analysis (MSA), a multivariate method based on coalitional game theory inferring causal regional contributions to specific behavioral outcomes from the characteristic functional deficits after stroke lesions. We established the causal patterns of contributions and interactions of nodes of the attentional orienting network on the basis of lesion and behavioral data from 25 right hemisphere stroke patients tested in visuo-spatial attention tasks. We calculated the percentage of damaged voxels for five right hemisphere cortical regions contributing to attentional orienting, involving seven specific Brodmann Areas (BA): Frontal Eye Fields, (FEF-BA6), Intraparietal Sulcus (IPS-BA7), Inferior Frontal Gyrus (IFG-BA44/BA45), Temporo-Parietal Junction (TPJ-BA39/BA40), and Inferior Occipital Gyrus (IOG-BA19). We computed the MSA contributions of these seven BAs to three behavioral clinical tests (line bisection, bells cancellation, and letter cancelation). Our analyses indicated IPS as the main contributor to the attentional orienting and also revealed synergistic influences among IPS, TPJ, and IOG (for bells cancellation and line bisection) and between TPJ and IFG (for bells and letter cancellation tasks). The findings demonstrate the ability of the MSA approach to infer plausible causal contributions of relevant right hemisphere sites in poststroke visuo-spatial attention and awareness disorders. Hum Brain Mapp 38:3454-3471, 2017. © 2017 Wiley Periodicals, Inc.

Modular topology emerges from plasticity in a minimalistic excitable network model.

Topological features play a major role in the emergence of complex brain network dynamics underlying brain function. Specific topological properties of brain networks, such as their modular organization, have been widely studied in recent years and shown to be ubiquitous across spatial scales and species. However, the mechanisms underlying the generation and maintenance of such features are still unclear. Using a minimalistic network model with excitable nodes and discrete deterministic dynamics, we studied the effects of a local Hebbian plasticity rule on global network topology. We found that, despite the simple model set-up, the plasticity rule was able to reorganize the global network topology into a modular structure. The structural reorganization was accompanied by enhanced correlations between structural and functional connectivity, and the final network organization reflected features of the dynamical model. These findings demonstrate the potential of simple plasticity rules for structuring the topology of brain connectivity.

Bridging Cytoarchitectonics and Connectomics in Human Cerebral Cortex.

The rich variation in cytoarchitectonics of the human cortex is well known to play an important role in the differentiation of cortical information processing, with functional multimodal areas noted to display more branched, more spinous, and an overall more complex cytoarchitecture. In parallel, connectome studies have suggested that also the macroscale wiring profile of brain areas may have an important contribution in shaping neural processes; for example, multimodal areas have been noted to display an elaborate macroscale connectivity profile. However, how these two scales of brain connectivity are related-and perhaps interact-remains poorly understood. In this communication, we combined data from the detailed mappings of early twentieth century cytoarchitectonic pioneers Von Economo and Koskinas (1925) on the microscale cellular structure of the human cortex with data on macroscale connectome wiring as derived from high-resolution diffusion imaging data from the Human Connectome Project. In a cross-scale examination, we show evidence of a significant association between cytoarchitectonic features of human cortical organization-in particular the size of layer 3 neurons-and whole-brain corticocortical connectivity. Our findings suggest that aspects of microscale cytoarchitectonics and macroscale connectomics are related. SIGNIFICANCE STATEMENT: One of the most widely known and perhaps most fundamental properties of the human cortex is its rich variation in cytoarchitectonics. At the same time, neuroimaging studies have also revealed cortical areas to vary in their level of macroscale connectivity. Here, we provide evidence that aspects of local cytoarchitecture are associated with aspects of global macroscale connectivity, providing insight into the question of how the scales of micro-organization and macro-organization of the human cortex are related.