Network Neuroscience Research Articles

CovSTATIS: A multi-table technique for network neuroscience

Similarity analyses between multiple correlation or covariance tables constitute the cornerstone of network neuroscience. Here, we introduce covSTATIS, a versatile, linear, unsupervised multi-table method designed to identify structured patterns in multi-table data, and allow for the simultaneous extraction and interpretation of both individual and group-level features. With covSTATIS, multiple similarity tables can now be easily integrated, without requiring a priori data simplification, complex black-box implementations, user-dependent specifications, or supervised frameworks. Applications of covSTATIS, a tutorial with Open Data and source code are provided. CovSTATIS offers a promising avenue for advancing the theoretical and analytic landscape of network neuroscience.

Mindfulness Meditation and Network Neuroscience: Review, Synthesis, and Future Directions

Mindfulness Meditation and Network Neuroscience: Review, Synthesis, and Future Directions

Contrasting topologies of synchronous and asynchronous functional brain networks

ABSTRACT We generated asynchronous functional networks (aFNs) using a novel method called optimal causation entropy and compared aFN topology with the correlation-based synchronous functional networks (sFNs), which are commonly used in network neuroscience studies. Functional magnetic resonance imaging (fMRI) time series from 212 participants of the National Consortium on Alcohol and Neurodevelopment in Adolescence study were used to generate aFNs and sFNs. As a demonstration of how aFNs and sFNs can be used in tandem, we used multivariate mixed effects models to determine whether age interacted with node efficiency to influence connection probabilities in the two networks. After adjusting for differences in network density, aFNs had higher global efficiency but lower local efficiency than the sFNs. In the aFNs, nodes with the highest outgoing global efficiency tended to be in the brainstem and orbitofrontal cortex. aFN nodes with the highest incoming global efficiency tended to be members of the default mode network in sFNs. Age interacted with node global efficiency in aFNs and node local efficiency in sFNs to influence connection probability. We conclude that the sFN and aFN both offer information about functional brain connectivity, which the other type of network does not.

Segregation of the regional radiomics similarity network exhibited an increase from late childhood to early adolescence: A developmental investigation

Brain development is characterized by an increase in structural and functional segregation, which supports the specialization of cognitive processes within the context of network neuroscience. In this study, we investigated age-related changes in morphological segregation using individual Regional Radiomics Similarity Networks (R2SNs) constructed with a longitudinal dataset of 494 T1-weighted MR scans from 309 typically developing children aged 6.2 to 13 years at baseline. Segertation indices were defined as the relative difference in connectivity strengths within and between modules and cacluated at the global, system and local levels. Linear mixed-effect models revealed longitudinal increases in both global and system segregation indices, particularly within the limbic and dorsal attention network, and decreases within the ventral attention network. Superior performance in working memory and inhibitory control was associated with higher system-level segregation indices in default, frontoparietal, ventral attention, somatomotor and subcortical systems, and lower local segregation indices in visual network regions, regardless of age. Furthermore, gene enrichment analysis revealed correlations between age-related changes in local segregation indices and regional expression levels of genes related to developmental processes. These findings provide novel insights into typical brain developmental changes using R2SN-derived segregation indices, offering a valuable tool for understanding human brain structural and cognitive maturation.

A Principled Framework to Assess the Information-Theoretic Fitness of Brain Functional Sub-Circuits

In systems and network neuroscience, many common practices in brain connectomic analysis are often not properly scrutinized. One such practice is mapping a predetermined set of sub-circuits, like functional networks (FNs), onto subjects’ functional connectomes (FCs) without adequately assessing the information-theoretic appropriateness of the partition. Another practice that goes unchallenged is thresholding weighted FCs to remove spurious connections without justifying the chosen threshold. This paper leverages recent theoretical advances in Stochastic Block Models (SBMs) to formally define and quantify the information-theoretic fitness (e.g., prominence) of a predetermined set of FNs when mapped to individual FCs under different fMRI task conditions. Our framework allows for evaluating any combination of FC granularity, FN partition, and thresholding strategy, thereby optimizing these choices to preserve the important topological features of the human brain connectomes. By applying to the Human Connectome Project with Schaefer parcellations at multiple levels of granularity, the framework showed that the common thresholding value of 0.25 was indeed information-theoretically valid for group-average FCs, despite its previous lack of justification. Our results pave the way for the proper use of FNs and thresholding methods, and provide insights for future research in individualized parcellations.

Distributed associations among white matter hyperintensities and structural brain networks with fluid cognition in healthy aging

White matter hyperintensities (WMHs) are associated with age-related cognitive impairment and increased risk of Alzheimer’s disease. However, the manner by which WMHs contribute to cognitive impairment is unclear. Using a combination of predictive modeling and network neuroscience, we investigated the relationship between structural white matter connectivity and age, fluid cognition, and WMHs in 68 healthy adults (18–78 years). Consistent with previous work, WMHs were increased in older adults and exhibited a strong negative association with fluid cognition. Extending previous work, using predictive modeling, we demonstrated that age, WMHs, and fluid cognition were jointly associated with widespread alterations in structural connectivity. Subcortical-cortical connections between the thalamus/basal ganglia and frontal and parietal regions of the default mode and frontoparietal networks were most prominent. At the network level, both age and WMHs were negatively associated with network density and communicability, and positively associated with modularity. Spatially, WMHs were most prominent in arterial zones served by the middle cerebral artery and associated lenticulostriate branches that supply subcortical regions. Finally, WMHs overlapped with all major white matter tracts, most prominently in tracts that facilitate subcortical-cortical communication and are implicated in fluid cognition, including the anterior thalamic-radiations and forceps minor. Finally, results of mediation analyses suggest that whole-brain WMH load influences age-related decline in fluid cognition. Thus, across multiple levels of analysis, we showed that WMHs were increased in older adults and associated with altered structural white matter connectivity and network topology involving subcortical-cortical pathways critical for fluid cognition.

The cognitive triad network - oscillation - behaviour links individual differences in EEG theta frequency with task performance and effective connectivity

We reconcile two significant lines of Cognitive Neuroscience research: the relationship between the structural and functional architecture of the brain and behaviour on the one hand and the functional significance of oscillatory brain processes to behavioural performance on the other. Network neuroscience proposes that the three elements, behavioural performance, EEG oscillation frequency, and network connectivity should be tightly connected at the individual level. Young and old healthy adults were recruited as a proxy for performance variation. An auditory inhibitory control task was used to demonstrate that task performance correlates with the individual EEG frontal theta frequency. Older adults had a significantly slower theta frequency, and both theta frequency and task performance correlated with the strengths of two network connections that involve the main areas of inhibitory control and speech processing. The results suggest that both the recruited functional network and the oscillation frequency induced by the task are specific to the task, are inseparable, and mark individual differences that directly link structure and function to behaviour in health and disease.

Hierarchical communities in the larval Drosophila connectome: Links to cellular annotations and network topology

One of the longstanding aims of network neuroscience is to link a connectome's topological properties-i.e., features defined from connectivity alone-with an organism's neurobiology. One approach for doing so is to compare connectome properties with annotational maps. This type of analysis is popular at the meso-/macroscale, but is less common at the nano-scale, owing to a paucity of neuron-level connectome data. However, recent methodological advances have made possible the reconstruction of whole-brain connectomes at single-neuron resolution for a select set of organisms. These include the fruit fly, Drosophila melanogaster, and its developing larvae. In addition to fine-scale descriptions of connectivity, these datasets are accompanied by rich annotations. Here, we use a variant of the stochastic blockmodel to detect multilevel communities in the larval Drosophila connectome. We find that communities partition neurons based on function and cell type and that most interact assortatively, reflecting the principle of functional segregation. However, a small number of communities interact nonassortatively, forming form a "rich-club" of interneurons that receive sensory/ascending inputs and deliver outputs along descending pathways. Next, we investigate the role of community structure in shaping communication patterns. We find that polysynaptic signaling follows specific trajectories across modular hierarchies, with interneurons playing a key role in mediating communication routes between modules and hierarchical scales. Our work suggests a relationship between system-level architecture and the biological function and classification of individual neurons. We envision our study as an important step toward bridging the gap between complex systems and neurobiological lines of investigation in brain sciences.

Stable individual differences from dynamic patterns of function: brain network flexibility predicts openness/intellect, intelligence, and psychoticism.

A growing understanding of the nature of brain function has led to increased interest in interpreting the properties of large-scale brain networks. Methodological advances in network neuroscience provide means to decompose these networks into smaller functional communities and measure how they reconfigure over time as an index of their dynamic and flexible properties. Recent evidence has identified associations between flexibility and a variety of traits pertaining to complex cognition including creativity and working memory. The present study used measures of dynamic resting-state functional connectivity in data from the Human Connectome Project (n = 994) to test associations with Openness/Intellect, general intelligence, and psychoticism, three traits that involve flexible cognition. Using a machine-learning cross-validation approach, we identified reliable associations of intelligence with cohesive flexibility of parcels in large communities across the cortex, of psychoticism with disjoint flexibility, and of Openness/Intellect with overall flexibility among parcels in smaller communities. These findings are reasonably consistent with previous theories of the neural correlates of these traits and help to expand on previous associations of behavior with dynamic functional connectivity, in the context of broad personality dimensions.

Task-specific topology of brain networks supporting working memory and inhibition.

Network neuroscience explores the brain's connectome, demonstrating that dynamic neural networks support cognitive functions. This study investigates how distinct cognitive abilities-working memory and cognitive inhibitory control-are supported by unique brain network configurations constructed by estimating whole-brain networks using mutual information. The study involved 195 participants who completed the Sternberg Item Recognition task and Flanker tasks while undergoing electroencephalography recording. A mixed-effects linear model analyzed the influence of network metrics on cognitive performance, considering individual differences and task-specific dynamics. The findings indicate that working memory and cognitive inhibitory control are associated with different network attributes, with working memory relying on distributed networks and cognitive inhibitory control on more segregated ones. Our analysis suggests that both strong and weak connections contribute to cognitive processes, with weak connections potentially leading to a more stable and support networks of memory and cognitive inhibitory control. The findings indirectly support the network neuroscience theory of intelligence, suggesting different functional topology of networks inherent to various cognitive functions. Nevertheless, we propose that understanding individual variations in cognitive abilities requires recognizing both shared and unique processes within the brain's network dynamics.

Stability of infraslow correlation structure in time-shifted intracranial EEG signals.

It is increasingly understood that the epilepsies are characterized by network pathology that can span multiple spatial and temporal scales. Recent work indicates that infraslow (<0.2Hz) envelope correlations may form a basis for distant spatial coupling in the brain. We speculated that infraslow correlation structure may be preserved even with some time lag between signals. To this end, we studied intracranial EEG (icEEG) data collected from 22 medically refractory epilepsy patients. For each patient, we selected hour-long background, awake icEEG epochs before and after antiseizure medication (ASM) taper. For each epoch, we selected 5,000 random electrode contact pairs and estimated magnitude-squared coherence (MSC) below 0.15Hz of band power time-series in the traditional EEG frequency bands. Using these same contact pairs, we shifted one signal of the pair by random durations in 15-s increments between 0 and 300s. We aggregated these data across all patients to determine how infraslow MSC varies with duration of lag. We further examined the effect of ASM taper on infraslow correlation structure. We also used surrogate data to empirically characterize MSC estimator and to set optimal parameters for estimation specifically for the study of infraslow activity. Our empirical analysis of the MSC estimator showed that hour-long segments with MSC computed using 3-min windows with 50% overlap was sufficient to capture infraslow envelope correlations while minimizing estimator bias and variance. The mean MSC decreased monotonically with increasing time lag until 105s of lag, then plateaued between 106 and 300s. Significantly nonzero infraslow envelope MSC was preserved in all frequency bands until about 1min of time lag, both pre- and post-ASM taper. We also saw a slight, but significant increase in infraslow MSC post-ASM taper, consistent with prior work. These results provide evidence for the feasibility of examining infraslow activity via its modulation of higher-frequency activity in the absence of DC-coupled recordings. The use of surrogate data also provides a general methodology for benchmarking measures used in network neuroscience studies. Finally, our study points to the clinical relevance of infraslow activity in assessing seizure risk.

Regression and alignment for functional data and network topology.

In the brain, functional connections form a network whose topological organization can be described by graph-theoretic network diagnostics. These include characterizations of the community structure, such as modularity and participation coefficient, which have been shown to change over the course of childhood and adolescence. To investigate if such changes in the functional network are associated with changes in cognitive performance during development, network studies often rely on an arbitrary choice of preprocessing parameters, in particular the proportional threshold of network edges. Because the choice of parameter can impact the value of the network diagnostic, and therefore downstream conclusions, we propose to circumvent that choice by conceptualizing the network diagnostic as a function of the parameter. As opposed to a single value, a network diagnostic curve describes the connectome topology at multiple scales-from the sparsest group of the strongest edges to the entire edge set. To relate these curves to executive function and other covariates, we use scalar-on-function regression, which is more flexible than previous functional data-based models used in network neuroscience. We then consider how systematic differences between networks can manifest in misalignment of diagnostic curves, and consequently propose a supervised curve alignment method that incorporates auxiliary information from other variables. Our algorithm performs both functional regression and alignment via an iterative, penalized, and nonlinear likelihood optimization. The illustrated method has the potential to improve the interpretability and generalizability of neuroscience studies where the goal is to study heterogeneity among a mixture of function- and scalar-valued measures.

A network control theory pipeline for studying the dynamics of the structural connectome.

Network control theory (NCT) is a simple and powerful tool for studying how network topology informs and constrains the dynamics of a system. Compared to other structure-function coupling approaches, the strength of NCT lies in its capacity to predict the patterns of external control signals that may alter the dynamics of a system in a desired way. An interesting development for NCT in the neuroscience field is its application to study behavior and mental health symptoms. To date, NCT has been validated to study different aspects of the human structural connectome. NCT outputs can be monitored throughout developmental stages to study the effects of connectome topology on neural dynamics and, separately, to test the coherence of empirical datasets with brain function and stimulation. Here, we provide a comprehensive pipeline for applying NCT to structural connectomes by following two procedures. The main procedure focuses on computing the control energy associated with the transitions between specific neural activity states. The second procedure focuses on computing average controllability, which indexes nodes' general capacity to control the dynamics of the system. We provide recommendations for comparing NCT outputs against null network models, and we further support this approach with a Python-based software package called 'network control theory for python'. The procedures in this protocol are appropriate for users with a background in network neuroscience and experience in dynamical systems theory.

Emergence and reconfiguration of modular structure for artificial neural networks during continual familiarity detection.

Advances in artificial intelligence enable neural networks to learn a wide variety of tasks, yet our understanding of the learning dynamics of these networks remains limited. Here, we study the temporal dynamics during learning of Hebbian feedforward neural networks in tasks of continual familiarity detection. Drawing inspiration from network neuroscience, we examine the network's dynamic reconfiguration, focusing on how network modules evolve throughout learning. Through a comprehensive assessment involving metrics like network accuracy, modular flexibility, and distribution entropy across diverse learning modes, our approach reveals various previously unknown patterns of network reconfiguration. We find that the emergence of network modularity is a salient predictor of performance and that modularization strengthens with increasing flexibility throughout learning. These insights not only elucidate the nuanced interplay of network modularity, accuracy, and learning dynamics but also bridge our understanding of learning in artificial and biological agents.

Contrasting topologies of synchronous and asynchronous functional brain networks.

We generated asynchronous functional networks (aFNs) using a novel method called optimal causation entropy (oCSE) and compared aFN topology to the correlation-based synchronous functional networks (sFNs) which are commonly used in network neuroscience studies. Functional magnetic resonance imaging (fMRI) time series from 212 participants of the National Consortium on Alcohol and NeuroDevelopment in Adolescence (NCANDA) study were used to generate aFNs and sFNs. As a demonstration of how aFNs and sFNs can be used in tandem, we used multivariate mixed effects models to determine whether age interacted with node efficiency to influence connection probabilities in the two networks. After adjusting for differences in network density, aFNs had higher global efficiency but lower local efficiency than the sFNs. In the aFNs, nodes with the highest outgoing global efficiency tended to be in the brainstem and orbitofrontal cortex. aFN nodes with the highest incoming global efficiency tended to be members of the Default Mode Network (DMN) in sFNs. Age interacted with node global efficiency in aFNs and node local efficiency in sFNs to influence connection probability. We conclude that the sFN and aFN both offer information about functional brain connectivity which the other type of network does not.

Brain functional gradient and structure features in adolescent and adult autism spectrum disorders.

Understanding how function and structure are organized and their coupling with clinical traits in individuals with autism spectrum disorder (ASD) is a primary goal in network neuroscience research for ASD. Atypical brain functional networks and structures in individuals with ASD have been reported, but whether these associations show heterogeneous hierarchy modeling in adolescents and adults with ASD remains to be clarified. In this study, 176 adolescent and 74 adult participants with ASD without medication or comorbidities and sex, age matched healthy controls (HCs) from 19 research groups from the openly shared Autism Brain Imaging Data Exchange II database were included. To investigate the relationship between the functional gradient, structural changes, and clinical symptoms of brain networks in adolescents and adults with ASD, functional gradient and voxel-based morphometry (VBM) analyses based on 1000 parcels defined by Schaefer mapped to Yeo's seven-network atlas were performed. Pearson's correlation was calculated between the gradient scores, gray volume and density, and clinical traits. The subsystem-level analysis showed that the second gradient scores of the default mode networks and frontoparietal network in patients with ASD were relatively compressed compared to adolescent HCs. Adult patients with ASD showed an overall compression gradient of 1 in the ventral attention networks. In addition, the gray density and volumes of the subnetworks showed no significant differences between the ASD and HC groups at the adolescent stage. However, adults with ASD showed decreased gray density in the limbic network. Moreover, numerous functional gradient parameters, but not VBM parameters, in adolescents with ASD were considerably correlated with clinical traits in contrast to those in adults with ASD. Our findings proved that the atypical changes in adolescent ASD mainly involve the brain functional network, while in adult ASD, the changes are more related to brain structure, including gray density and volume. These changes in functional gradients or structures are markedly correlated with clinical traits in patients with ASD. Our study provides a novel understanding of the pathophysiology of the structure-function hierarchy in ASD.

Simulation-based inference on virtual brain models of disorders

Connectome-based models, also known as virtual brain models (VBMs), have been well established in network neuroscience to investigate pathophysiological causes underlying a large range of brain diseases. The integration of an individual’s brain imaging data in VBMs has improved patient-specific predictivity, although Bayesian estimation of spatially distributed parameters remains challenging even with state-of-the-art Monte Carlo sampling. VBMs imply latent nonlinear state space models driven by noise and network input, necessitating advanced probabilistic machine learning techniques for widely applicable Bayesian estimation. Here we present simulation-based inference on VBMs (SBI-VBMs), and demonstrate that training deep neural networks on both spatio-temporal and functional features allows for accurate estimation of generative parameters in brain disorders. The systematic use of brain stimulation provides an effective remedy for the non-identifiability issue in estimating the degradation limited to smaller subset of connections. By prioritizing model structure over data, we show that the hierarchical structure in SBI-VBMs renders the inference more effective, precise and biologically plausible. This approach could broadly advance precision medicine by enabling fast and reliable prediction of patient-specific brain disorders.

Effects of packetization on communication dynamics in brain networks.

Computational studies in network neuroscience build models of communication dynamics in the connectome that help us understand the structure-function relationships of the brain. In these models, the dynamics of cortical signal transmission in brain networks are approximated with simple propagation strategies such as random walks and shortest path routing. Furthermore, the signal transmission dynamics in brain networks can be associated with the switching architectures of engineered communication systems (e.g., message switching and packet switching). However, it has been unclear how propagation strategies and switching architectures are related in models of brain network communication. Here, we investigate the effects of the difference between packet switching and message switching (i.e., whether signals are packetized or not) on the transmission completion time of propagation strategies when simulating signal propagation in mammalian brain networks. The results show that packetization in the connectome with hubs increases the time of the random walk strategy and does not change that of the shortest path strategy, but decreases that of more plausible strategies for brain networks that balance between communication speed and information requirements. This finding suggests an advantage of packet-switched communication in the connectome and provides new insights into modeling the communication dynamics in brain networks.

The multiscale topological organization of the functional brain network in adolescent PTSD.

The experience of an extremely aversive event can produce enduring deleterious behavioral, and neural consequences, among which posttraumatic stress disorder (PTSD) is a representative example. Although adolescence is a period of great exposure to potentially traumatic events, the effects of trauma during adolescence remain understudied in clinical neuroscience. In this exploratory work, we aim to study the whole-cortex functional organization of 14 adolescents with PTSD using a data-driven method tailored to our population of interest. To do so, we built on the network neuroscience framework and specifically on multilayer (multisubject) community analysis to study the functional connectivity of the brain. We show, across different topological scales (the number of communities composing the cortex), a hyper-colocalization between regions belonging to occipital and pericentral regions and hypo-colocalization in middle temporal, posterior-anterior medial, and frontal cortices in the adolescent PTSD group compared to a nontrauma exposed group of adolescents. These preliminary results raise the question of an altered large-scale cortical organization in adolescent PTSD, opening an interesting line of research for future investigations.

From Connectivity to Controllability: Unraveling the Brain Biomarkers of Major Depressive Disorder

Major Depressive Disorder (MDD) is a significant neurological condition associated with aberrations in brain functional networks. Traditional studies have predominantly analyzed these from a network topology perspective. However, given the brain’s dynamic and complex nature, exploring its mechanisms from a network control standpoint provides a fresh and insightful framework. This research investigates the integration of network controllability and machine learning to pinpoint essential biomarkers for MDD using functional magnetic resonance imaging (fMRI) data. By employing network controllability methods, we identify crucial brain regions that are instrumental in facilitating transitions between brain states. These regions demonstrate the brain’s ability to navigate various functional states, emphasizing the utility of network controllability metrics as potential biomarkers. Furthermore, these metrics elucidate the complex dynamics of MDD and support the development of precision medicine strategies that incorporate machine learning to improve the precision of diagnostics and the efficacy of treatments. This study underscores the value of merging machine learning with network neuroscience to craft personalized interventions that align with the unique pathological profiles of individuals, ultimately enhancing the management and treatment of MDD.