Structural-functional brain network coupling predicts human cognitive ability

Individual differences in general cognitive ability (GCA) have a biological basis within the structure and function of the human brain. Network neuroscience investigations revealed neural correlates of GCA in structural as well as in functional brain networks. However, whether the relationship between structural and functional networks, the structural-functional brain network coupling (SC-FC coupling), is related to individual differences in GCA remains an open question. We used data from 1030 adults of the Human Connectome Project, derived structural connectivity from diffusion weighted imaging, functional connectivity from resting-state fMRI, and assessed GCA as a latent g-factor from 12 cognitive tasks. Two similarity measures and six communication measures were used to model possible functional interactions arising from structural brain networks. SC-FC coupling was estimated as the degree to which these measures align with the actual functional connectivity, providing insights into different neural communication strategies. At the whole-brain level, higher GCA was associated with higher SC-FC coupling, but only when considering path transitivity as neural communication strategy. Taking region-specific variations in the SC-FC coupling strategy into account and differentiating between positive and negative associations with GCA, allows for prediction of individual cognitive ability scores in a cross-validated prediction framework (correlation between predicted and observed scores: r = 0.25, p < .001). The same model also predicts GCA scores in a completely independent sample (N = 567, r = 0.19, p < .001). Our results propose structural-functional brain network coupling as a neurobiological correlate of GCA and suggest brain region-specific coupling strategies as neural basis of efficient information processing predictive of cognitive ability.

Spatial transcriptomic patterns underlying amyloid-β and tau pathology are associated with cognitive dysfunction in Alzheimer's disease

Amyloid-β (Aβ) and tau proteins accumulate within distinct neuronal systems in Alzheimer's disease (AD). Although it is not clear why certain brain regions are more vulnerable to Aβ and tau pathologies than others, gene expression may play a role. We study the association between brain-wide gene expression profiles and regional vulnerability to Aβ (gene-to-Aβ associations) and tau (gene-to-tau associations) pathologies by leveraging two large independent AD cohorts. We identify AD susceptibility genes and gene modules in a gene co-expression network with expression profiles specifically related to regional vulnerability to Aβ and tau pathologies in AD. In addition, we identify distinct biochemical pathways associated with the gene-to-Aβ and the gene-to-tau associations. These findings may explain the discordance between regional Aβ and tau pathologies. Finally, we propose an analytic framework, linking the identified gene-to-pathology associations to cognitive dysfunction in AD at the individual level, suggesting potential clinical implication of the gene-to-pathology associations.

Relation of connectome topology to brain volume across 103 mammalian species

The brain connectome is an embedded network of anatomically interconnected brain regions, and the study of its topological organization in mammals has become of paramount importance due to its role in scaffolding brain function and behavior. Unlike many other observable networks, brain connections incur material and energetic cost, and their length and density are volumetrically constrained by the skull. Thus, an open question is how differences in brain volume impact connectome topology. We address this issue using the MaMI database, a diverse set of mammalian connectomes reconstructed from 201 animals, covering 103 species and 12 taxonomy orders, whose brain size varies over more than 4 orders of magnitude. Our analyses focus on relationships between volume and modular organization. After having identified modules through a multiresolution approach, we observed how connectivity features relate to the modular structure and how these relations vary across brain volume. We found that as the brain volume increases, modules become more spatially compact and dense, comprising more costly connections. Furthermore, we investigated how spatial embedding shapes network communication, finding that as brain volume increases, nodes' distance progressively impacts communication efficiency. We identified modes of variation in network communication policies, as smaller and bigger brains show higher efficiency in routing- and diffusion-based signaling, respectively. Finally, bridging network modularity and communication, we found that in larger brains, modular structure imposes stronger constraints on network signaling. Altogether, our results show that brain volume is systematically related to mammalian connectome topology and that spatial embedding imposes tighter restrictions on larger brains.

Edge time series components of functional connectivity and cognitive function in Alzheimer's disease

Understanding the interrelationships of brain function as measured by resting-state magnetic resonance imaging and neuropsychological/behavioral measures in Alzheimer's disease is key for advancement of neuroimaging analysis methods in clinical research. The edge time-series framework recently developed in the field of network neuroscience, in combination with other network science methods, allows for investigations of brain-behavior relationships that are not possible with conventional functional connectivity methods. Data from the Indiana Alzheimer's Disease Research Center sample (53 cognitively normal control, 47 subjective cognitive decline, 32 mild cognitive impairment, and 20 Alzheimer's disease participants) were used to investigate relationships between functional connectivity components, each derived from a subset of time points based on co-fluctuation of regional signals, and measures of domain-specific neuropsychological functions. Multiple relationships were identified with the component approach that were not found with conventional functional connectivity. These involved attentional, limbic, frontoparietal, and default mode systems and their interactions, which were shown to couple with cognitive, executive, language, and attention neuropsychological domains. Additionally, overlapping results were obtained with two different statistical strategies (network contingency correlation analysis and network-based statistics correlation). Results demonstrate that connectivity components derived from edge time-series based on co-fluctuation reveal disease-relevant relationships not observed with conventional static functional connectivity.

Network architecture of intrinsic connectivity in a mammalian spinal cord (the central nervous system's caudal sector)

The vertebrate spinal cord (SP) is the long, thin extension of the brain forming the central nervous system's caudal sector. Functionally, the SP directly mediates motor and somatic sensory interactions with most parts of the body except the face, and it is the preferred model for analyzing relatively simple reflex behaviors. Here, we analyze the organization of axonal connections between the 50 gray matter regions forming the bilaterally symmetric rat SP. The assembled dataset suggests that there are about 385 of a possible 2,450 connections between the 50 regions for a connection density of 15.7%. Multiresolution consensus cluster analysis reveals a hierarchy of structure-function subsystems in this neural network, with 4 subsystems at the top level and 12 at the bottom-level. The top-level subsystems include a) a bilateral subsystem related most clearly to somatic and autonomic motor functions and centered in the ventral horn and intermediate zone; b) a bilateral subsystem associated with general somatosensory functions and centered in the base, neck, and head of the dorsal horn; and c) a pair of unilateral, bilaterally symmetric subsystems associated with nociceptive information processing and occupying the apex of the dorsal horn. The intrinsic SP network displayed no hubs, rich club, or small-world attributes, which are common measures of global functionality. Advantages and limitations of our methodology are discussed in some detail. The present work is part of a comprehensive project to assemble and analyze the neurome of a mammalian nervous system and its interactions with the body.

Intrinsic circuitry of the rhombicbrain (central nervous system's intermediate sector) in a mammal

The rhombicbrain (rhombencephalon or intermediate sector) is the vertebrate central nervous system part between the forebrain-midbrain (rostral sector) and spinal cord (caudal sector), and it has three main divisions: pons, cerebellum, and medulla. Using a data-driven approach, here we examine intrinsic rhombicbrain (intrarhombicbrain) network architecture that in rat consists of 52,670 possible axonal connections between 230 gray matter regions (115 bilaterally symmetrical pairs). Our analysis indicates that only 8,089 (15.4%) of these connections exist. Multiresolution consensus cluster analysis yields a nested hierarchy model of rhombicbrain subsystems that at the top level are associated with 1) the cerebellum and vestibular nuclei, 2) orofacial-pharyngeal-visceral integration, and 3) auditory connections; the bottom level has 68 clusters, ranging in size from 2 to 11 regions. The model provides a basis for functional hypothesis development and interrogation. More granular network analyses performed on the intrinsic connectivity of individual and combined main rhombicbrain divisions (pons, cerebellum, medulla, pons + cerebellum, and pons + medulla) demonstrate the mutability of network architecture in response to the addition or subtraction of connections. Clear differences between the structure-function network architecture of the rhombicbrain and forebrain-midbrain are discussed, with a stark comparison provided by the subsystem and small-world organization of the cerebellar cortex and cerebral cortex. Future analysis of the connections within and between the forebrain-midbrain and rhombicbrain will provide a model of brain neural network architecture in a mammal.

Edge Time Series Components of Functional Connectivity and Cognitive Function in Alzheimer's Disease

Understanding the interrelationships of brain function as measured by resting-state magnetic resonance imaging and neuropsychological/behavioral measures in Alzheimer's disease is key for advancement of neuroimaging analysis methods in clinical research. The edge time-series framework recently developed in the field of network neuroscience, in combination with other network science methods, allows for investigations of brain-behavior relationships that are not possible with conventional functional connectivity methods. Data from the Indiana Alzheimer's Disease Research Center sample (53 cognitively normal control, 47 subjective cognitive decline, 32 mild cognitive impairment, and 20 Alzheimer's disease participants) were used to investigate relationships between functional connectivity components, each derived from a subset of time points based on co-fluctuation of regional signals, and measures of domain-specific neuropsychological functions. Multiple relationships were identified with the component approach that were not found with conventional functional connectivity. These involved attentional, limbic, frontoparietal, and default mode systems and their interactions, which were shown to couple with cognitive, executive, language, and attention neuropsychological domains. Additionally, overlapping results were obtained with two different statistical strategies (network contingency correlation analysis and network-based statistics correlation). Results demonstrate that connectivity components derived from edge time-series based on co-fluctuation reveal disease-relevant relationships not observed with conventional static functional connectivity.

Living on the edge: network neuroscience beyond nodes

Network neuroscience has emphasized the connectional properties of neural elements - cells, populations, and regions. This has come at the expense of the anatomical and functional connections that link these elements to one another. A new perspective - namely one that emphasizes 'edges' - may prove fruitful in addressing outstanding questions in network neuroscience. We highlight one recently proposed 'edge-centric' method and review its current applications, merits, and limitations. We also seek to establish conceptual and mathematical links between this method and previously proposed approaches in the network science and neuroimaging literature. We conclude by presenting several avenues for future work to extend and refine existing edge-centric analysis.

Hierarchical organization of spontaneous co-fluctuations in densely sampled individuals using fMRI

Edge time series decompose functional connectivity into its framewise contributions. Previous studies have focused on characterizing the properties of high-amplitude frames (time points when the global co-fluctuation amplitude takes on its largest value), including their cluster structure. Less is known about middle- and low-amplitude co-fluctuations (peaks in co-fluctuation time series but of lower amplitude). Here, we directly address those questions, using data from two dense-sampling studies: the MyConnectome project and Midnight Scan Club. We develop a hierarchical clustering algorithm to group peak co-fluctuations of all magnitudes into nested and multiscale clusters based on their pairwise concordance. At a coarse scale, we find evidence of three large clusters that, collectively, engage virtually all canonical brain systems. At finer scales, however, each cluster is dissolved, giving way to increasingly refined patterns of co-fluctuations involving specific sets of brain systems. We also find an increase in global co-fluctuation magnitude with hierarchical scale. Finally, we comment on the amount of data needed to estimate co-fluctuation pattern clusters and implications for brain-behavior studies. Collectively, the findings reported here fill several gaps in current knowledge concerning the heterogeneity and richness of co-fluctuation patterns as estimated with edge time series while providing some practical guidance for future studies.

A low dimensional embedding of brain dynamics enhances diagnostic accuracy and behavioral prediction in stroke

Large-scale brain networks reveal structural connections as well as functional synchronization between distinct regions of the brain. The latter, referred to as functional connectivity (FC), can be derived from neuroimaging techniques such as functional magnetic resonance imaging (fMRI). FC studies have shown that brain networks are severely disrupted by stroke. However, since FC data are usually large and high-dimensional, extracting clinically useful information from this vast amount of data is still a great challenge, and our understanding of the functional consequences of stroke remains limited. Here, we propose a dimensionality reduction approach to simplify the analysis of this complex neural data. By using autoencoders, we find a low-dimensional representation encoding the fMRI data which preserves the typical FC anomalies known to be present in stroke patients. By employing the latent representations emerging from the autoencoders, we enhanced patients' diagnostics and severity classification. Furthermore, we showed how low-dimensional representation increased the accuracy of recovery prediction.

Brain network communication: concepts, models and applications

Understanding communication and information processing in nervous systems is a central goal of neuroscience. Over the past two decades, advances in connectomics and network neuroscience have opened new avenues for investigating polysynaptic communication in complex brain networks. Recent work has brought into question the mainstay assumption that connectome signalling occurs exclusively via shortest paths, resulting in a sprawling constellation of alternative network communication models. This Review surveys the latest developments in models of brain network communication. We begin by drawing a conceptual link between the mathematics of graph theory and biological aspects of neural signalling such as transmission delays and metabolic cost. We organize key network communication models and measures into a taxonomy, aimed at helping researchers navigate the growing number of concepts and methods in the literature. The taxonomy highlights the pros, cons and interpretations of different conceptualizations of connectome signalling. We showcase the utility of network communication models as a flexible, interpretable and tractable framework to study brain function by reviewing prominent applications in basic, cognitive and clinical neurosciences. Finally, we provide recommendations to guide the future development, application and validation of network communication models.

Neuroscience Needs Network Science

The brain is a complex system comprising a myriad of interacting neurons, posing significant challenges in understanding its structure, function, and dynamics. Network science has emerged as a powerful tool for studying such interconnected systems, offering a framework for integrating multiscale data and complexity. To date, network methods have significantly advanced functional imaging studies of the human brain and have facilitated the development of control theory-based applications for directing brain activity. Here, we discuss emerging frontiers for network neuroscience in the brain atlas era, addressing the challenges and opportunities in integrating multiple data streams for understanding the neural transitions from development to healthy function to disease. We underscore the importance of fostering interdisciplinary opportunities through workshops, conferences, and funding initiatives, such as supporting students and postdoctoral fellows with interests in both disciplines. By bringing together the network science and neuroscience communities, we can develop novel network-based methods tailored to neural circuits, paving the way toward a deeper understanding of the brain and its functions, as well as offering new challenges for network science.

Spatial transcriptomic patterns underlying regional vulnerability to amyloid-β and tau pathologies and their relationships to cognitive dysfunction in Alzheimer's disease

Amyloid-β (Aβ) and tau proteins accumulate within distinct neuronal systems in Alzheimer's disease (AD). Although it is not clear why certain brain regions are more vulnerable to Aβ and tau pathologies than others, gene expression may play a role. We studied the association between brain-wide gene expression profiles and regional vulnerability to Aβ (gene-to-Aβ associations) and tau (gene-to-tau associations) pathologies leveraging two large independent cohorts (n = 715) of participants along the AD continuum. We identified several AD susceptibility genes and gene modules in a gene co-expression network with expression profiles related to regional vulnerability to Aβ and tau pathologies in AD. In particular, we found that the positive APOE -to-tau association was only seen in the AD cohort, whereas patients with AD and frontotemporal dementia shared similar positive MAPT -to-tau association. Some AD candidate genes showed sex-dependent negative gene-to-Aβ and gene-to-tau associations. In addition, we identified distinct biochemical pathways associated with the gene-to-Aβ and the gene-to-tau associations. Finally, we proposed a novel analytic framework, linking the identified gene-to-pathology associations to cognitive dysfunction in AD at the individual level, suggesting potential clinical implication of the gene-to-pathology associations. Taken together, our study identified distinct gene expression profiles and biochemical pathways that may explain the discordance between regional Aβ and tau pathologies, and filled the gap between gene-to-pathology associations and cognitive dysfunction in individual AD patients that may ultimately help identify novel personalized pathogenetic biomarkers and therapeutic targets.

One Sentence Summary

We identified replicable cognition-related associations between regional gene expression profiles and selectively regional vulnerability to amyloid-β and tau pathologies in AD.

Co-evolving dynamics and topology in a coupled oscillator model of resting brain function

Dynamic models of ongoing BOLD fMRI brain dynamics and models of communication strategies have been two important approaches to understanding how brain network structure constrains function. However, dynamic models have yet to widely incorporate one of the most important insights from communication models: the brain may not use all of its connections in the same way or at the same time. Here we present a variation of a phase delayed Kuramoto coupled oscillator model that dynamically limits communication between nodes on each time step. An active subgraph of the empirically derived anatomical brain network is chosen in accordance with the local dynamic state on every time step, thus coupling dynamics and network structure in a novel way. We analyze this model with respect to its fit to empirical time-averaged functional connectivity, finding that, with the addition of only one parameter, it significantly outperforms standard Kuramoto models with phase delays. We also perform analyses on the novel time series of active edges it produces, demonstrating a slowly evolving topology moving through intermittent episodes of integration and segregation. We hope to demonstrate that the exploration of novel modeling mechanisms and the investigation of dynamics of networks in addition to dynamics on networks may advance our understanding of the relationship between brain structure and function.

Partial entropy decomposition reveals higher-order information structures in human brain activity

The standard approach to modeling the human brain as a complex system is with a network, where the basic unit of interaction is a pairwise link between two brain regions. While powerful, this approach is limited by the inability to assess higher-order interactions involving three or more elements directly. In this work, we explore a method for capturing higher-order dependencies in multivariate data: the partial entropy decomposition (PED). Our approach decomposes the joint entropy of the whole system into a set of nonnegative atoms that describe the redundant, unique, and synergistic interactions that compose the system's structure. PED gives insight into the mathematics of functional connectivity and its limitation. When applied to resting-state fMRI data, we find robust evidence of higher-order synergies that are largely invisible to standard functional connectivity analyses. Our approach can also be localized in time, allowing a frame-by-frame analysis of how the distributions of redundancies and synergies change over the course of a recording. We find that different ensembles of regions can transiently change from being redundancy-dominated to synergy-dominated and that the temporal pattern is structured in time. These results provide strong evidence that there exists a large space of unexplored structures in human brain data that have been largely missed by a focus on bivariate network connectivity models. This synergistic structure is dynamic in time and likely will illuminate interesting links between brain and behavior. Beyond brain-specific application, the PED provides a very general approach for understanding higher-order structures in a variety of complex systems.

Fine-scale dynamics of functional connectivity in the face-processing network during movie watching

Mapping the human face-processing network is typically done during rest or using isolated, static face images, overlooking widespread cortical interactions obtained in response to naturalistic face dynamics and context. To determine how inter-subject functional correlation (ISFC) relates to face recognition scores, we measure cortical connectivity patterns in response to a dynamic movie in typical adults (N = 517). We find a positive correlation with recognition scores in edges connecting the occipital visual and anterior temporal regions and a negative correlation in edges connecting the attentional dorsal, frontal default, and occipital visual regions. We measure the inter-subject stimulus-evoked response at a single TR resolution and demonstrate that co-fluctuations in face-selective edges are related to activity in core face-selective regions and that the ISFC patterns peak during boundaries between movie segments rather than during the presence of faces. Our approach demonstrates how face processing is linked to fine-scale dynamics in attentional, memory, and perceptual neural circuitry.

Few temporally distributed brain connectivity states predict human cognitive abilities

Human functional brain connectivity can be temporally decomposed into states of high and low cofluctuation, defined as coactivation of brain regions over time. Rare states of particularly high cofluctuation have been shown to reflect fundamentals of intrinsic functional network architecture and to be highly subject-specific. However, it is unclear whether such network-defining states also contribute to individual variations in cognitive abilities - which strongly rely on the interactions among distributed brain regions. By introducing CMEP, a new eigenvector-based prediction framework, we show that as few as 16 temporally separated time frames (< 1.5% of 10min resting-state fMRI) can significantly predict individual differences in intelligence (N = 263, p < .001). Against previous expectations, individual's network-defining time frames of particularly high cofluctuation do not predict intelligence. Multiple functional brain networks contribute to the prediction, and all results replicate in an independent sample (N = 831). Our results suggest that although fundamentals of person-specific functional connectomes can be derived from few time frames of highest connectivity, temporally distributed information is necessary to extract information about cognitive abilities. This information is not restricted to specific connectivity states, like network-defining high-cofluctuation states, but rather reflected across the entire length of the brain connectivity time series.

Neuroscience needs Network Science

The brain is a complex system comprising a myriad of interacting elements, posing significant challenges in understanding its structure, function, and dynamics. Network science has emerged as a powerful tool for studying such intricate systems, offering a framework for integrating multiscale data and complexity. Here, we discuss the application of network science in the study of the brain, addressing topics such as network models and metrics, the connectome, and the role of dynamics in neural networks. We explore the challenges and opportunities in integrating multiple data streams for understanding the neural transitions from development to healthy function to disease, and discuss the potential for collaboration between network science and neuroscience communities. We underscore the importance of fostering interdisciplinary opportunities through funding initiatives, workshops, and conferences, as well as supporting students and postdoctoral fellows with interests in both disciplines. By uniting the network science and neuroscience communities, we can develop novel network-based methods tailored to neural circuits, paving the way towards a deeper understanding of the brain and its functions.

Multivariate information theory uncovers synergistic subsystems of the human cerebral cortex

One of the most well-established tools for modeling the brain is the functional connectivity network, which is constructed from pairs of interacting brain regions. While powerful, the network model is limited by the restriction that only pairwise dependencies are considered and potentially higher-order structures are missed. Here, we explore how multivariate information theory reveals higher-order dependencies in the human brain. We begin with a mathematical analysis of the O-information, showing analytically and numerically how it is related to previously established information theoretic measures of complexity. We then apply the O-information to brain data, showing that synergistic subsystems are widespread in the human brain. Highly synergistic subsystems typically sit between canonical functional networks, and may serve an integrative role. We then use simulated annealing to find maximally synergistic subsystems, finding that such systems typically comprise ≈10 brain regions, recruited from multiple canonical brain systems. Though ubiquitous, highly synergistic subsystems are invisible when considering pairwise functional connectivity, suggesting that higher-order dependencies form a kind of shadow structure that has been unrecognized by established network-based analyses. We assert that higher-order interactions in the brain represent an under-explored space that, accessible with tools of multivariate information theory, may offer novel scientific insights.

Communication dynamics in the human connectome shape the cortex-wide propagation of direct electrical stimulation

Communication between gray matter regions underpins all facets of brain function. We study inter-areal communication in the human brain using intracranial EEG recordings, acquired following 29,055 single-pulse direct electrical stimulations in a total of 550 individuals across 20 medical centers (average of 87 ± 37 electrode contacts per subject). We found that network communication models-computed on structural connectivity inferred from diffusion MRI-can explain the causal propagation of focal stimuli, measured at millisecond timescales. Building on this finding, we show that a parsimonious statistical model comprising structural, functional, and spatial factors can accurately and robustly predict cortex-wide effects of brain stimulation (R²=46% in data from held-out medical centers). Our work contributes toward the biological validation of concepts in network neuroscience and provides insight into how connectome topology shapes polysynaptic inter-areal signaling. We anticipate that our findings will have implications for research on neural communication and the design of brain stimulation paradigms.

Uncovering individual differences in fine-scale dynamics of functional connectivity

Functional connectivity (FC) profiles contain subject-specific features that are conserved across time and have potential to capture brain-behavior relationships. Most prior work has focused on spatial features (nodes and systems) of these FC fingerprints, computed over entire imaging sessions. We propose a method for temporally filtering FC, which allows selecting specific moments in time while also maintaining the spatial pattern of node-based activity. To this end, we leverage a recently proposed decomposition of FC into edge time series (eTS). We systematically analyze functional magnetic resonance imaging frames to define features that enhance identifiability across multiple fingerprinting metrics, similarity metrics, and data sets. Results show that these metrics characteristically vary with eTS cofluctuation amplitude, similarity of frames within a run, transition velocity, and expression of functional systems. We further show that data-driven optimization of features that maximize fingerprinting metrics isolates multiple spatial patterns of system expression at specific moments in time. Selecting just 10% of the data can yield stronger fingerprints than are obtained from the full data set. Our findings support the idea that FC fingerprints are differentially expressed across time and suggest that multiple distinct fingerprints can be identified when spatial and temporal characteristics are considered simultaneously.

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.

Information-processing dynamics in neural networks of macaque cerebral cortex reflect cognitive state and behavior

One of the essential functions of biological neural networks is the processing of information. This includes everything from processing sensory information to perceive the environment, up to processing motor information to interact with the environment. Due to methodological limitations, it has been historically unclear how information processing changes during different cognitive or behavioral states and to what extent information is processed within or between the network of neurons in different brain areas. In this study, we leverage recent advances in the calculation of information dynamics to explore neural-level processing within and between the frontoparietal areas AIP, F5, and M1 during a delayed grasping task performed by three macaque monkeys. While information processing was high within all areas during all cognitive and behavioral states of the task, interareal processing varied widely: During visuomotor transformation, AIP and F5 formed a reciprocally connected processing unit, while no processing was present between areas during the memory period. Movement execution was processed globally across all areas with predominance of processing in the feedback direction. Furthermore, the fine-scale network structure reconfigured at the neuron level in response to different grasping conditions, despite no differences in the overall amount of information present. These results suggest that areas dynamically form higher-order processing units according to the cognitive or behavioral demand and that the information-processing network is hierarchically organized at the neuron level, with the coarse network structure determining the behavioral state and finer changes reflecting different conditions.

The complex brain: connectivity, dynamics, information

Most would agree, the brain is complex. But, beyond metaphor, does the brain's complexity demand a paradigm shift in how we study its structure and function? I argue that complexity manifests in three domains - connectivity, dynamics, and information - and that unlocking their interactions will greatly advance our understanding of brain and cognition.

Multi-modal and multi-subject modular organization of human brain networks

The human brain is a complex network of anatomically interconnected brain areas. Spontaneous neural activity is constrained by this architecture, giving rise to patterns of statistical dependencies between the activity of remote neural elements. The non-trivial relationship between structural and functional connectivity poses many unsolved challenges about cognition, disease, development, learning and aging. While numerous studies have focused on statistical relationships between edge weights in anatomical and functional networks, less is known about dependencies between their modules and communities. In this work, we investigate and characterize the relationship between anatomical and functional modular organization of the human brain, developing a novel multi-layer framework that expands the classical concept of multi-layer modularity. By simultaneously mapping anatomical and functional networks estimated from different subjects into communities, this approach allows us to carry out a multi-subject and multi-modal analysis of the brain's modular organization. Here, we investigate the relationship between anatomical and functional modules during resting state, finding unique and shared structures. The proposed framework constitutes a methodological advance in the context of multi-layer network analysis and paves the way to further investigate the relationship between structural and functional network organization in clinical cohorts, during cognitively demanding tasks, and in developmental or lifespan studies.

A connectomics-based taxonomy of mammals

Mammalian taxonomies are conventionally defined by morphological traits and genetics. How species differ in terms of neural circuits and whether inter-species differences in neural circuit organization conform to these taxonomies is unknown. The main obstacle for the comparison of neural architectures have been differences in network reconstruction techniques, yielding species-specific connectomes that are not directly comparable to one another. Here we comprehensively chart connectome organization across the mammalian phylogenetic spectrum using a common reconstruction protocol. We analyze the mammalian MRI (MaMI) data set, a database that encompasses high-resolution ex vivo structural and diffusion magnetic resonance imaging (MRI) scans of 124 species across 12 taxonomic orders and 5 superorders, collected using a unified MRI protocol. We assess similarity between species connectomes using two methods: similarity of Laplacian eigenspectra and similarity of multiscale topological features. We find greater inter-species similarities among species within the same taxonomic order, suggesting that connectome organization reflects established taxonomic relationships defined by morphology and genetics. While all connectomes retain hallmark global features and relative proportions of connection classes, inter-species variation is driven by local regional connectivity profiles. By encoding connectomes into a common frame of reference, these findings establish a foundation for investigating how neural circuits change over phylogeny, forging a link from genes to circuits to behaviour.

Editorial: Focus feature on consciousness and cognition

Consciousness and cognition are an increasing focus of theoretical and experimental research in neuroscience, leveraging the methods and tools of brain dynamics and connectivity. This Focus Feature brings together a collection of articles that examine the various roles of brain networks in computational and dynamic models, and in studies of physiological and neuroimaging processes that underpin and enable behavioral and cognitive function.

Edge-centric analysis of time-varying functional brain networks with applications in autism spectrum disorder

The interaction between brain regions changes over time, which can be characterized using time-varying functional connectivity (tvFC). The common approach to estimate tvFC uses sliding windows and offers limited temporal resolution. An alternative method is to use the recently proposed edge-centric approach, which enables the tracking of moment-to-moment changes in co-fluctuation patterns between pairs of brain regions. Here, we first examined the dynamic features of edge time series and compared them to those in the sliding window tvFC (sw-tvFC). Then, we used edge time series to compare subjects with autism spectrum disorder (ASD) and healthy controls (CN). Our results indicate that relative to sw-tvFC, edge time series captured rapid and bursty network-level fluctuations that synchronize across subjects during movie-watching. The results from the second part of the study suggested that the magnitude of peak amplitude in the collective co-fluctuations of brain regions (estimated as root sum square (RSS) of edge time series) is similar in CN and ASD. However, the trough-to-trough duration in RSS signal is greater in ASD, compared to CN. Furthermore, an edge-wise comparison of high-amplitude co-fluctuations showed that the within-network edges exhibited greater magnitude fluctuations in CN. Our findings suggest that high-amplitude co-fluctuations captured by edge time series provide details about the disruption of functional brain dynamics that could potentially be used in developing new biomarkers of mental disorders.

Self-organization of in vitro neuronal assemblies drives to complex network topology

Activity-dependent self-organization plays an important role in the formation of specific and stereotyped connectivity patterns in neural circuits. By combining neuronal cultures, and tools with approaches from network neuroscience and information theory, we can study how complex network topology emerges from local neuronal interactions. We constructed effective connectivity networks using a transfer entropy analysis of spike trains recorded from rat embryo dissociated hippocampal neuron cultures between 6 and 35 days in vitro to investigate how the topology evolves during maturation. The methodology for constructing the networks considered the synapse delay and addressed the influence of firing rate and population bursts as well as spurious effects on the inference of connections. We found that the number of links in the networks grew over the course of development, shifting from a segregated to a more integrated architecture. As part of this progression, three significant aspects of complex network topology emerged. In agreement with previous in silico and in vitro studies, a small-world architecture was detected, largely due to strong clustering among neurons. Additionally, the networks developed in a modular topology, with most modules comprising nearby neurons. Finally, highly active neurons acquired topological characteristics that made them important nodes to the network and integrators of modules. These findings leverage new insights into how neuronal effective network topology relates to neuronal assembly self-organization mechanisms.

Edge-centric analysis of stroke patients: An alternative approach for biomarkers of lesion recovery

Most neuroimaging studies of post-stroke recovery rely on analyses derived from standard node-centric functional connectivity to map the distributed effects in stroke patients. Here, given the importance of nonlocal and diffuse damage, we use an edge-centric approach to functional connectivity in order to provide an alternative description of the effects of this disorder. These techniques allow for the rendering of metrics such as normalized entropy, which describes the diversity of edge communities at each node. Moreover, the approach enables the identification of high amplitude co-fluctuations in fMRI time series. We found that normalized entropy is associated with stroke lesion severity and continually increases across the time of patients' recovery. Furthermore, high amplitude co-fluctuations not only relate to the lesion severity but are also associated with patients' level of recovery. The current study is the first edge-centric application for a clinical population in a longitudinal dataset and demonstrates how a different perspective for functional data analysis can further characterize topographic modulations of brain dynamics.

Network communication models narrow the gap between the modular organization of structural and functional brain networks

Structural and functional brain networks are modular. Canonical functional systems, such as the default mode network, are well-known modules of the human brain and have been implicated in a large number of cognitive, behavioral and clinical processes. However, modules delineated in structural brain networks inferred from tractography generally do not recapitulate canonical functional systems. Neuroimaging evidence suggests that functional connectivity between regions in the same systems is not always underpinned by anatomical connections. As such, direct structural connectivity alone would be insufficient to characterize the functional modular organization of the brain. Here, we demonstrate that augmenting structural brain networks with models of indirect (polysynaptic) communication unveils a modular network architecture that more closely resembles the brain's established functional systems. We find that diffusion models of polysynaptic connectivity, particularly communicability, narrow the gap between the modular organization of structural and functional brain networks by 20-60%, whereas routing models based on single efficient paths do not improve mesoscopic structure-function correspondence. This suggests that functional modules emerge from the constraints imposed by local network structure that facilitates diffusive neural communication. Our work establishes the importance of modeling polysynaptic communication to understand the structural basis of functional systems.

Network Analysis of Time Series: Novel Approaches to Network Neuroscience

In the last two decades, there has been an explosion of interest in modeling the brain as a network, where nodes correspond variously to brain regions or neurons, and edges correspond to structural or statistical dependencies between them. This kind of network construction, which preserves spatial, or structural, information while collapsing across time, has become broadly known as "network neuroscience." In this work, we provide an alternative application of network science to neural data: network-based analysis of non-linear time series and review applications of these methods to neural data. Instead of preserving spatial information and collapsing across time, network analysis of time series does the reverse: it collapses spatial information, instead preserving temporally extended dynamics, typically corresponding to evolution through some kind of phase/state-space. This allows researchers to infer a, possibly low-dimensional, "intrinsic manifold" from empirical brain data. We will discuss three methods of constructing networks from nonlinear time series, and how to interpret them in the context of neural data: recurrence networks, visibility networks, and ordinal partition networks. By capturing typically continuous, non-linear dynamics in the form of discrete networks, we show how techniques from network science, non-linear dynamics, and information theory can extract meaningful information distinct from what is normally accessible in standard network neuroscience approaches.

Multitask brain network reconfiguration is inversely associated with human intelligence

Intelligence describes the general cognitive ability level of a person. It is one of the most fundamental concepts in psychological science and is crucial for the effective adaption of behavior to varying environmental demands. Changing external task demands have been shown to induce reconfiguration of functional brain networks. However, whether neural reconfiguration between different tasks is associated with intelligence has not yet been investigated. We used functional magnetic resonance imaging data from 812 subjects to show that higher scores of general intelligence are related to less brain network reconfiguration between resting state and seven different task states as well as to network reconfiguration between tasks. This association holds for all functional brain networks except the motor system and replicates in two independent samples (n = 138 and n = 184). Our findings suggest that the intrinsic network architecture of individuals with higher intelligence scores is closer to the network architecture as required by various cognitive demands. Multitask brain network reconfiguration may, therefore, represent a neural reflection of the behavioral positive manifold - the essence of the concept of general intelligence. Finally, our results support neural efficiency theories of cognitive ability and reveal insights into human intelligence as an emergent property from a distributed multitask brain network.

Edges in brain networks: Contributions to models of structure and function

Network models describe the brain as sets of nodes and edges that represent its distributed organization. So far, most discoveries in network neuroscience have prioritized insights that highlight distinct groupings and specialized functional contributions of network nodes. Importantly, these functional contributions are determined and expressed by the web of their interrelationships, formed by network edges. Here, we underscore the important contributions made by brain network edges for understanding distributed brain organization. Different types of edges represent different types of relationships, including connectivity and similarity among nodes. Adopting a specific definition of edges can fundamentally alter how we analyze and interpret a brain network. Furthermore, edges can associate into collectives and higher order arrangements, describe time series, and form edge communities that provide insights into brain network topology complementary to the traditional node-centric perspective. Focusing on the edges, and the higher order or dynamic information they can provide, discloses previously underappreciated aspects of structural and functional network organization.

Modular community structure of the face network supports face recognition

Face recognition is dependent on computations conducted in specialized brain regions and the communication among them, giving rise to the face-processing network. We examined whether modularity of this network may underlie the vast individual differences found in human face recognition abilities. Modular networks, characterized by strong within and weaker between-network connectivity, were previously suggested to promote efficacy and reduce interference among cognitive systems and also correlated with better cognitive abilities. The study was conducted in a large sample (n = 409) with diffusion-weighted imaging, resting-state fMRI, and a behavioral face recognition measure. We defined a network of face-selective regions and derived a novel measure of communication along with structural and functional connectivity among them. The modularity of this network was positively correlated with recognition abilities even when controlled for age. Furthermore, the results were specific to the face network when compared with the place network or to spatially permuted null networks. The relation to behavior was also preserved at the individual-edge level such that a larger correlation to behavior was found within hemispheres and particularly within the right hemisphere. This study provides the first evidence of modularity-behavior relationships in the domain of face processing and more generally in visual perception.

Additive and interaction effects of working memory and motor sequence training on brain functional connectivity

Although shared behavioral and neural mechanisms between working memory (WM) and motor sequence learning (MSL) have been suggested, the additive and interactive effects of training have not been studied. This study aimed at investigating changes in brain functional connectivity (FC) induced by sequential (WM + MSL and MSL + WM) and combined (WM × MSL) training programs. 54 healthy subjects (27 women; mean age: 30.2 ± 8.6 years) allocated to three training groups underwent twenty-four 40-min training sessions over 6 weeks and four cognitive assessments including functional MRI. A double-baseline approach was applied to account for practice effects. Test performances were compared using linear mixed-effects models and t-tests. Resting state fMRI data were analysed using FSL. Processing speed, verbal WM and manual dexterity increased following training in all groups. MSL + WM training led to additive effects in processing speed and verbal WM. Increased FC was found after training in a network including the right angular gyrus, left superior temporal sulcus, right superior parietal gyrus, bilateral middle temporal gyri and left precentral gyrus. No difference in FC was found between double baselines. Results indicate distinct patterns of resting state FC modulation related to sequential and combined WM and MSL training suggesting a relevance of the order of training performance. These observations could provide new insight for the planning of effective training/rehabilitation.

Modular origins of high-amplitude cofluctuations in fine-scale functional connectivity dynamics

The topology of structural brain networks shapes brain dynamics, including the correlation structure of brain activity (functional connectivity) as estimated from functional neuroimaging data. Empirical studies have shown that functional connectivity fluctuates over time, exhibiting patterns that vary in the spatial arrangement of correlations among segregated functional systems. Recently, an exact decomposition of functional connectivity into frame-wise contributions has revealed fine-scale dynamics that are punctuated by brief and intermittent episodes (events) of high-amplitude cofluctuations involving large sets of brain regions. Their origin is currently unclear. Here, we demonstrate that similar episodes readily appear in silico using computational simulations of whole-brain dynamics. As in empirical data, simulated events contribute disproportionately to long-time functional connectivity, involve recurrence of patterned cofluctuations, and can be clustered into distinct families. Importantly, comparison of event-related patterns of cofluctuations to underlying patterns of structural connectivity reveals that modular organization present in the coupling matrix shapes patterns of event-related cofluctuations. Our work suggests that brief, intermittent events in functional dynamics are partly shaped by modular organization of structural connectivity.

The human connectome in Alzheimer disease - relationship to biomarkers and genetics

The pathology of Alzheimer disease (AD) damages structural and functional brain networks, resulting in cognitive impairment. The results of recent connectomics studies have now linked changes in structural and functional network organization in AD to the patterns of amyloid-β and tau accumulation and spread, providing insights into the neurobiological mechanisms of the disease. In addition, the detection of gene-related connectome changes might aid in the early diagnosis of AD and facilitate the development of personalized therapeutic strategies that are effective at earlier stages of the disease spectrum. In this article, we review studies of the associations between connectome changes and amyloid-β and tau pathologies as well as molecular genetics in different subtypes and stages of AD. We also highlight the utility of connectome-derived computational models for replicating empirical findings and for tracking and predicting the progression of biomarker-indicated AD pathophysiology.

Mapping individual differences across brain network structure to function and behavior with connectome embedding

The connectome, a comprehensive map of the brain’s anatomical connections, is often summarized as a matrix comprising all dyadic connections among pairs of brain regions. This representation cannot capture higher-order relations within the brain graph. Connectome embedding (CE) addresses this limitation by creating compact vectorized representations of brain nodes capturing their context in the global network topology. Here, nodes “context” is defined as random walks on the brain graph and as such, represents a generative model of diffusive communication around nodes. Applied to group-averaged structural connectivity, CE was previously shown to capture relations between inter-hemispheric homologous brain regions and uncover putative missing edges from the network reconstruction. Here we extend this framework to explore individual differences with a novel embedding alignment approach. We test this approach in two lifespan datasets (NKI: n = 542; Cam-CAN: n = 601) that include diffusion-weighted imaging, resting-state fMRI, demographics and behavioral measures. We demonstrate that modeling functional connectivity with CE substantially improves structural to functional connectivity mapping both at the group and subject level. Furthermore, age-related differences in this structure-function mapping, are preserved and enhanced. Importantly, CE captures individual differences by out-of-sample prediction of age and intelligence. The resulting predictive accuracy was higher compared to using structural connectivity and functional connectivity. We attribute these findings to the capacity of the CE to incorporate aspects of both anatomy (the structural graph) and function (diffusive communication). Our novel approach allows mapping individual differences in the connectome through structure to function and behavior.

Topological analysis of differential effects of ketamine and propofol anaesthesia on brain dynamics

Research has found that the vividness of conscious experience is related to brain dynamics. Despite both being anaesthetics, propofol and ketamine produce different subjective states: we explore the different effects of these two anaesthetics on the structure of dynamic attractors reconstructed from electrophysiological activity recorded from cerebral cortex of two macaques. We used two methods: the first embeds the recordings in a continuous high-dimensional manifold on which we use topological data analysis to infer the presence of higher-order dynamics. The second reconstruction, an ordinal partition network embedding, allows us to create a discrete state-transition network, which is amenable to information-theoretic analysis and contains rich information about state-transition dynamics. We find that the awake condition generally had the 'richest' structure, visiting the most states, the presence of pronounced higher-order structures, and the least deterministic dynamics. By contrast, the propofol condition had the most dissimilar dynamics, transitioning to a more impoverished, constrained, low-structure regime. The ketamine condition, interestingly, seemed to combine aspects of both: while it was generally less complex than the awake condition, it remained well above propofol in almost all measures. These results provide deeper and more comprehensive insights than what is typically gained by using point-measures of complexity.

Temporal stability of the ventral attention network and general cognition along the Alzheimer's disease spectrum

Understanding the interrelationships of clinical manifestations of Alzheimer's disease (AD) and functional connectivity (FC) as the disease progresses is necessary for use of FC as a potential neuroimaging biomarker. Degradation of resting-state networks in AD has been observed when FC is estimated over the entire scan, however, the temporal dynamics of these networks are less studied. We implemented a novel approach to investigate the modular structure of static (sFC) and time-varying (tvFC) connectivity along the AD spectrum in a two-sample Discovery/Validation design (n = 80 and 81, respectively). Cortical FC networks were estimated across 4 diagnostic groups (cognitively normal, subjective cognitive decline, mild cognitive impairment, and AD) for whole scan (sFC) and with sliding window correlation (tvFC). Modularity quality (across a range of spatial scales) did not differ in either sFC or tvFC. For tvFC, group differences in temporal stability within and between multiple resting state networks were observed; however, these differences were not consistent between samples. Correlation analyses identified a relationship between global cognition and temporal stability of the ventral attention network, which was reproduced in both samples. While the ventral attention system has been predominantly studied in task-evoked designs, the relationship between its intrinsic dynamics at-rest and general cognition along the AD spectrum highlights its relevance regarding clinical manifestation of the disease.

Subject identification using edge-centric functional connectivity

Group-level studies do not capture individual differences in network organization, an important prerequisite for understanding neural substrates shaping behavior and for developing interventions in clinical conditions. Recent studies have employed 'fingerprinting' analyses on functional connectivity to identify subjects' idiosyncratic features. Here, we develop a complementary approach based on an edge-centric model of functional connectivity, which focuses on the co-fluctuations of edges. We first show whole-brain edge functional connectivity (eFC) to be a robust substrate that improves identifiability over nodal FC (nFC) across different datasets and parcellations. Next, we characterize subjects' identifiability at different spatial scales, from single nodes to the level of functional systems and clusters using k-means clustering. Across spatial scales, we find that heteromodal brain regions exhibit consistently greater identifiability than unimodal, sensorimotor, and limbic regions. Lastly, we show that identifiability can be further improved by reconstructing eFC using specific subsets of its principal components. In summary, our results highlight the utility of the edge-centric network model for capturing meaningful subject-specific features and sets the stage for future investigations into individual differences using edge-centric models.

Dynamic expression of brain functional systems disclosed by fine-scale analysis of edge time series

Functional connectivity (FC) describes the statistical dependence between neuronal populations or brain regions in resting-state fMRI studies and is commonly estimated as the Pearson correlation of time courses. Clustering or community detection reveals densely coupled sets of regions constituting resting-state networks or functional systems. These systems manifest most clearly when FC is sampled over longer epochs but appear to fluctuate on shorter timescales. Here, we propose a new approach to reveal temporal fluctuations in neuronal time series. Unwrapping FC signal correlations yields pairwise co-fluctuation time series, one for each node pair or edge, and allows tracking of fine-scale dynamics across the network. Co-fluctuations partition the network, at each time step, into exactly two communities. Sampled over time, the overlay of these bipartitions, a binary decomposition of the original time series, very closely approximates functional connectivity. Bipartitions exhibit characteristic spatiotemporal patterns that are reproducible across participants and imaging runs, capture individual differences, and disclose fine-scale temporal expression of functional systems. Our findings document that functional systems appear transiently and intermittently, and that FC results from the overlay of many variable instances of system expression. Potential applications of this decomposition of functional connectivity into a set of binary patterns are discussed.

Anesthetics fragment hippocampal network activity, alter spine dynamics, and affect memory consolidation

General anesthesia is characterized by reversible loss of consciousness accompanied by transient amnesia. Yet, long-term memory impairment is an undesirable side effect. How different types of general anesthetics (GAs) affect the hippocampus, a brain region central to memory formation and consolidation, is poorly understood. Using extracellular recordings, chronic 2-photon imaging, and behavioral analysis, we monitor the effects of isoflurane (Iso), medetomidine/midazolam/fentanyl (MMF), and ketamine/xylazine (Keta/Xyl) on network activity and structural spine dynamics in the hippocampal CA1 area of adult mice. GAs robustly reduced spiking activity, decorrelated cellular ensembles, albeit with distinct activity signatures, and altered spine dynamics. CA1 network activity under all 3 anesthetics was different to natural sleep. Iso anesthesia most closely resembled unperturbed activity during wakefulness and sleep, and network alterations recovered more readily than with Keta/Xyl and MMF. Correspondingly, memory consolidation was impaired after exposure to Keta/Xyl and MMF, but not Iso. Thus, different anesthetics distinctly alter hippocampal network dynamics, synaptic connectivity, and memory consolidation, with implications for GA strategy appraisal in animal research and clinical settings.

Differential effects of propofol and ketamine on critical brain dynamics

Whether the brain operates at a critical "tipping" point is a long standing scientific question, with evidence from both cellular and systems-scale studies suggesting that the brain does sit in, or near, a critical regime. Neuroimaging studies of humans in altered states of consciousness have prompted the suggestion that maintenance of critical dynamics is necessary for the emergence of consciousness and complex cognition, and that reduced or disorganized consciousness may be associated with deviations from criticality. Unfortunately, many of the cellular-level studies reporting signs of criticality were performed in non-conscious systems (in vitro neuronal cultures) or unconscious animals (e.g. anaesthetized rats). Here we attempted to address this knowledge gap by exploring critical brain dynamics in invasive ECoG recordings from multiple sessions with a single macaque as the animal transitioned from consciousness to unconsciousness under different anaesthetics (ketamine and propofol). We use a previously-validated test of criticality: avalanche dynamics to assess the differences in brain dynamics between normal consciousness and both drug-states. Propofol and ketamine were selected due to their differential effects on consciousness (ketamine, but not propofol, is known to induce an unusual state known as "dissociative anaesthesia"). Our analyses indicate that propofol dramatically restricted the size and duration of avalanches, while ketamine allowed for more awake-like dynamics to persist. In addition, propofol, but not ketamine, triggered a large reduction in the complexity of brain dynamics. All states, however, showed some signs of persistent criticality when testing for exponent relations and universal shape-collapse. Further, maintenance of critical brain dynamics may be important for regulation and control of conscious awareness.