Hierarchical Heterogeneity across Human Cortex Shapes Large-Scale Neural Dynamics

Exploring dynamical whole-brain models in high-dimensional parameter spaces

Personalized modeling of the resting-state brain activity implies the usage of dynamical whole-brain models with high-dimensional model parameter spaces. However, the practical benefits and mathematical challenges originating from such approaches have not been thoroughly documented, leaving the question of the value and utility of high-dimensional approaches unanswered. Studying a whole-brain model of coupled phase oscillators, we proceeded from low-dimensional scenarios featuring 2-3 global model parameters only to high-dimensional cases, where we additionally equipped every brain region with a specific local model parameter. To enable the parameter optimizations for the high-dimensional model fitting to empirical data, we applied two dedicated mathematical optimization algorithms (Bayesian Optimization, Covariance Matrix Adaptation Evolution Strategy). We thereby optimized up to 103 parameters simultaneously with the aim to maximize the correlation between simulated and empirical functional connectivity separately for 272 subjects. The obtained model parameters demonstrated increased variability within subjects and reduced reliability across repeated optimization runs in high-dimensional spaces. Nevertheless, the quality of the model validation (goodness-of-fit, GoF) improved considerably and remained very stable and reliable together with the simulated functional connectivity. Applying the modeling results to phenotypical data, we found significantly higher prediction accuracies for sex classification when the GoF or coupling parameter values optimized in the high-dimensional spaces were considered as features. Our results elucidate the model fitting in high-dimensional parameter spaces and can contribute to an improved dynamical brain modeling as well as its application to the frameworks of inter-individual variability and brain-behavior relationships.

Selective control of prefrontal neural timescales by parietal cortex

Intrinsic neural timescales quantify how long spontaneous neuronal activity patterns persist, reflecting dynamics of endogenous fluctuations. We measured intrinsic timescales of frontal eye field (FEF) neurons and examined their changes during posterior parietal cortex (PPC) inactivation. We observed two distinct classes of FEF neurons based on their intrinsic timescales: short-timescale neurons (∼25 ms) or long-timescale neurons (∼100 ms). Short-timescale neurons showed stronger transient visual responses, suggesting their role in rapid visual processing, whereas long-timescale neurons exhibited pronounced sustained attentional modulation, suggesting their role in maintaining stimulus-driven attention. During PPC inactivation, intrinsic timescales increased in both neuron types, with a significantly larger effect in short-timescale neurons. In addition, PPC inactivation reduced attentional modulation, particularly in long-timescale neurons. Our findings provide the first causal evidence linking intrinsic local neural timescales to long-range inter-area communications. These findings also suggest the presence of at least two distinct network motifs that support different neuronal dynamics and functional computations within the FEF.

Mapping effective connectivity by virtually perturbing a surrogate brain

Effective connectivity (EC), which reflects the causal interactions between brain regions, is fundamental to understanding information processing in the brain; however, traditional methods for obtaining EC, which rely on neural responses to stimulation, are often invasive or limited in spatial coverage, making them unsuitable for whole-brain EC mapping in humans. Here, to address this gap, we introduce Neural Perturbational Inference (NPI), a data-driven framework for mapping whole-brain EC. NPI employs an artificial neural network trained to model large-scale neural dynamics, serving as a computational surrogate of the brain. By systematically perturbing all regions in the surrogate brain and analyzing the resulting responses in other regions, NPI maps the directionality, strength and excitatory/inhibitory properties of brain-wide EC. Validation of NPI on generative models with known ground-truth EC demonstrates its superiority over existing methods such as Granger causality and dynamic causal modeling. When applied to resting-state functional magnetic resonance imaging data across diverse datasets, NPI reveals consistent, structurally supported EC patterns. Furthermore, comparisons with cortico-cortical evoked potential data show a strong resemblance between NPI-inferred EC and real stimulation propagation patterns. By transitioning from correlational to causal understandings of brain functionality, NPI marks a stride in decoding the brain's functional architecture and facilitating both neuroscience studies and clinical applications.

Uncertainty mapping and probabilistic tractography using Simulation-based Inference in diffusion MRI: A comparison with classical Bayes

Simulation-Based Inference (SBI) has recently emerged as a powerful framework for Bayesian inference: Neural networks are trained on simulations from a forward model, and learn to rapidly estimate posterior distributions. We here present an SBI framework for parametric spherical deconvolution of diffusion MRI data of the brain. We demonstrate its utility for estimating white matter fibre orientations, mapping uncertainty of voxel-based estimates and performing probabilistic tractography by spatially propagating fibre orientation uncertainty. We conduct an extensive comparison against established Bayesian methods based on Markov-Chain Monte-Carlo (MCMC) and find that: a) in-silico training can lead to calibrated SBI networks with accurate parameter estimates and uncertainty mapping for both single- and multi-shell diffusion MRI, b) SBI allows amortised inference of the posterior distribution of model parameters given unseen observations, which is orders of magnitude faster than MCMC, c) SBI-based tractography yields reconstructions that have a high level of agreement with their MCMC-based counterparts, equal to or higher than scan-rescan reproducibility of estimates. We further demonstrate how SBI design considerations (such as dealing with noise, defining priors and handling model selection) can affect performance, allowing us to identify optimal practices. Taken together, our results show that SBI provides a powerful alternative to classical Bayesian inference approaches for fast and accurate model estimation and uncertainty mapping in MRI.

Optimizing Biophysical Large-Scale Brain Circuit Models With Deep Neural Networks

Biophysical modeling provides mechanistic insights into brain function, from single-neuron dynamics to large-scale circuit models bridging macro-scale brain activity with microscale measurements. Biophysical models are governed by biologically meaningful parameters, many of which can be experimentally measured. Some parameters are unknown, and optimizing their values can dramatically improve adherence to experimental data, significantly enhancing biological plausibility. Previous optimization methods - such as exhaustive search, gradient descent, evolutionary strategies and Bayesian optimization - require repeated, computationally expensive numerical integration of biophysical differential equations, limiting scalability to population-level datasets. Here, we introduce DELSSOME (DEep Learning for Surrogate Statistics Optimization in MEan field modeling), a framework that bypasses numerical integration by directly predicting whether model parameters produce realistic brain dynamics. When applied to the widely used feedback inhibition control (FIC) mean field model, DELSSOME achieves a 2000× speedup over Euler integration. By embedding DELSSOME within an evolutionary optimization strategy, trained models generalize to new datasets without additional tuning, enabling a 50× speedup in FIC model estimation while preserving neurobiological insights. The massive acceleration facilitates large-scale mechanistic modeling in population-level neuroscience, unlocking new opportunities for understanding brain function.

Adaptive Whole-Brain Dynamics Predictive Method: Relevancy to Mental Disorders

The Hopf whole-brain model, based on structural connectivity, overcomes limitations of traditional structural or functional connectivity-focused methods by incorporating heterogeneity parameters, quantifying dynamic brain characteristics in healthy and diseased states. Traditional parameter fitting techniques lack precision, restricting broader use. To address this, we validated parameter fitting methods using simulated networks and synthetic models, introducing improvements such as individual-specific initialization and optimized gradient descent, which reduced individual data loss. We also developed an approximate loss function and gradient adjustment mechanism, enhancing parameter fitting accuracy and stability. Applying this refined method to datasets for major depressive disorder (MDD) and autism spectrum disorder (ASD), we identified differences in brain regions between patients and healthy controls, explaining related anomalies. This rigorous validation is crucial for clinical application, paving the way for precise neuropathological identification and novel treatments in neuropsychiatric research, demonstrating substantial potential in clinical neurology.

Modeling the interplay between regional heterogeneity and critical dynamics underlying brain functional networks

The human brain exhibits heterogeneity across regions and network connectivity patterns; However, how these heterogeneities contribute to whole-brain network functions and cognitive capacities remains unclear. In this study, we focus on the regional heterogeneity reflected in local dynamics and study how it contributes to the emergence of functional connectivity patterns, network ignition dynamics of the empirical brains. We find that the level of synchrony among voxelwise neural activities measured from the fMRI data is significantly correlated with the transcriptional variations in excitatory and inhibitory receptor gene expression. Consequently, we construct heterogeneous whole-brain network models with nodal excitability calibrated by the synchronization measure of regional dynamics. We demonstrate that as the extent of heterogeneity increases, the models operating around the critical point between order and disorder generate simulated functional connectivity networks increasingly similar to empirical resting-state or working memory task-evoked function connectivity networks. Furthermore, the heterogeneous models can predict individual differences in resting-state and task-evoked reconfiguration of the functional connectivity, as well as the comparative causal effect of empirical brain networks-that is, how the dynamics of one brain region affect whole-brain synchronization. Overall, this work demonstrates the viability of using regional heterogeneous functional signals to improve the performance of the whole-brain models, and illustrates how regional heterogeneity in human brains interplays with structural connections and critical dynamics to contribute to the emergence of functional connectivity networks.

Network Mechanisms Underlying the Regional Diversity of Variance and Time Scales of the Brain's Spontaneous Activity Fluctuations

The brain's activity fluctuations have different temporal scales across the brain regions, with associative regions displaying slower timescales than sensory areas. This hierarchy of timescales has been shown to correlate with both structural brain connectivity and intrinsic regional properties. Here, using publicly available human resting-state fMRI and dMRI data, it was found that, while more structurally connected brain regions presented activity fluctuations with longer timescales, their activity fluctuations presented lower variance. The opposite relationships between the structural connectivity and the variance and temporal scales of resting-state fluctuations, respectively, were not trivially explained by simple network propagation principles. To understand these structure-function relationships, two commonly used whole-brain models were studied, namely, the Hopf and Wilson-Cowan models. These models use the brain's connectome to couple local nodes (representing brain regions) displaying noise-driven oscillations. The models show that the variance and temporal scales of activity fluctuations can oppositely relate to connectivity within specific parameter regions, even when all nodes have the same intrinsic dynamics-but also when intrinsic dynamics are constrained by the myelinization-related macroscopic gradient. These results show that, setting aside intrinsic regional differences, connectivity and network state are sufficient to explain the regional differences in fluctuations' scales. State dependence supports the vision that structure-function relationships can serve as biomarkers of altered brain states. Finally, the results indicate that the hierarchies of timescales and variances reflect a balance between stability and responsivity, with greater and faster responsiveness at the network periphery, while the network core ensures overall system robustness.

Understanding cognitive processes across spatial scales of the brain

Cognition arises from neural operations at multiple spatial scales, from individual neurons to large-scale networks. Despite extensive research on coding principles and emergent cognitive processes across brain areas, investigation across scales has been limited. Here, we propose ways to test the idea that different cognitive processes emerge from distinct information coding principles at various scales, which collectively give rise to complex behavior. This approach involves comparing brain-behavior associations and the underlying neural geometry across scales, alongside an investigation of global and local scale interactions. Bridging findings across species and techniques through open science and collaborations is essential to comprehensively understand the multiscale brain and its functions.

Test-retest reliability of resting-state EEG intrinsic neural timescales

Intrinsic neural timescales, which reflect the duration of neural information storage within local brain regions and capacity for information integration, are typically measured using autocorrelation windows (ACWs). Extraction of intrinsic neural timescales from resting-state brain activity has been extensively applied in psychiatric disease research. Given the potential of intrinsic neural timescales as a neural marker for psychiatric disorders, investigating their reliability is crucial. This study, using an open-source database, aimed to evaluate the test-retest reliability of ACW-0 and ACW-50 under both eyes-open and eyes-closed conditions across three sessions. The intraclass correlation coefficients (ICCs) were employed to quantify the reliability of the intrinsic neural timescales. Our results showed that intrinsic neural timescales exhibited good reliability (ICC > 0.6) at the whole-brain level across different index types and eye states. Spatially, except for the right temporal region in the eyes-open condition, all other regions showed moderate-to-high ICCs. Over 60% of the electrodes demonstrated moderate-to-high intrinsic neural timescale ICCs under both eyes-open and eyes-closed conditions, with ACW-0 being more stable than ACW-50. Moreover, in the new dataset, the above results were consistently reproduced. The present study comprehensively assessed the reliability of intrinsic neural timescale under various conditions, providing robust evidence for their stability in neuroscience and psychiatry.

A multimodal characterization of low-dimensional thalamocortical structural connectivity patterns

The human thalamus is a heterogeneous subcortical structure coordinating whole-brain activity. Investigations of its internal organization reveal differentiable subnuclei, however, a consensus on subnuclei boundaries remains absent. Recent work suggests that thalamic organization additionally reflects continuous axes transcending nuclear boundaries. Here, we study how low-dimensional axes of thalamocortical structural connectivity relate to intrathalamic microstructural features, functional connectivity, and structural covariance. Using diffusion MRI, we compute a thalamocortical structural connectome and derive two main axes of thalamic organization. The principal axis, extending from medial to lateral, relates to intrathalamic myelin, and functional connectivity organization. The secondary axis corresponds to the core-matrix cell distribution. Lastly, exploring multimodal associations globally, we observe the principal axis consistently differentiating limbic, frontoparietal, and default mode network nodes from dorsal and ventral attention networks across modalities. However, the link with sensory modalities varies. In sum, we show the coherence between lower dimensional patterns of thalamocortical structural connectivity and various modalities, shedding light on multiscale thalamic organization.

Excitation-Inhibition Balance, Neural Criticality, and Activities in Neuronal Circuits

Neural activities in local circuits exhibit complex and multilevel dynamic features. Individual neurons spike irregularly, which is believed to originate from receiving balanced amounts of excitatory and inhibitory inputs, known as the excitation-inhibition balance. The spatial-temporal cascades of clustered neuronal spikes occur in variable sizes and durations, manifested as neural avalanches with scale-free features. These may be explained by the neural criticality hypothesis, which posits that neural systems operate around the transition between distinct dynamic states. Here, we summarize the experimental evidence for and the underlying theory of excitation-inhibition balance and neural criticality. Furthermore, we review recent studies of excitatory-inhibitory networks with synaptic kinetics as a simple solution to reconcile these two apparently distinct theories in a single circuit model. This provides a more unified understanding of multilevel neural activities in local circuits, from spontaneous to stimulus-response dynamics.

Connectional axis of individual functional variability: Patterns, structural correlates, and relevance for development and cognition

The human cerebral cortex exhibits intricate interareal functional synchronization at the macroscale, with substantial individual variability in these functional connections. However, the spatial organization of functional connectivity (FC) variability across the human connectome edges and its significance in cognitive development remain unclear. Here, we identified a connectional axis in the edge-level FC variability. The variability declined continuously along this axis from within-network to between-network connections, and from the edges linking association networks to those linking the sensorimotor and association networks. This connectional axis of functional variability is associated with spatial pattern of structural connectivity variability. Moreover, the connectional variability axis evolves in youth with an increasing flatter axis slope. We also observed that the slope of connectional variability axis was positively related to the performance in the higher-order cognition. Together, our results reveal a connectional axis in functional variability that is linked with structural connectome variability, refines during development, and is relevant to cognition.

Precision data-driven modeling of cortical dynamics reveals person-specific mechanisms underpinning brain electrophysiology

Task-free brain activity affords unique insight into the functional structure of brain network dynamics and has been used to identify neural markers of individual differences. In this work, we present an algorithmic optimization framework that directly inverts and parameterizes brain-wide dynamical-systems models involving hundreds of interacting neural populations, from single-subject M/EEG time-series recordings. This technique provides a powerful neurocomputational tool for interrogating mechanisms underlying individual brain dynamics ("precision brain models") and making quantitative predictions. We extensively validate the models' performance in forecasting future brain activity and predicting individual variability in key M/EEG metrics. Last, we demonstrate the power of our technique in resolving individual differences in the generation of alpha and beta-frequency oscillations. We characterize subjects based upon model attractor topology and a dynamical-systems mechanism by which these topologies generate individual variation in the expression of alpha vs. beta rhythms. We trace these phenomena back to global variation in excitatory-inhibitory balance, highlighting the explanatory power of our framework to generate mechanistic insights.

Functional connectome gradient predicts clinical symptoms of chronic insomnia disorder

Insomnia is the second most prevalent psychiatric disorder worldwide, but the understanding of the pathophysiology of insomnia remains fragmented. In this study, we calculated the connectome gradient in 50 chronic insomnia disorder (CID) patients and 38 healthy controls (HC) to assess changes due to insomnia and utilized these gradients in a connectome-based predictive modeling (CPM) to predict clinical symptoms associated with insomnia. The results suggested that insomnia led to significant alterations in the functional gradients of some brain areas. Specifically, the gradient scores in the middle frontal gyrus, superior anterior cingulate gyrus, and right nucleus accumbens were significantly higher in the CID patients than in the HC group, whereas the scores in the middle occipital gyrus, right fusiform gyrus, and right postcentral gyrus were significantly lower than in the HC group. Further correlation analysis revealed that the right middle frontal gyrus is positively correlated with the self-rating anxiety scale (r=0.3702). Additionally, the prediction model built with functional gradients could well predict the sleep quality (r=0.5858), anxiety (r=0.6150), and depression (r=0.4022) levels of insomnia patients. This offers an objective depiction of the clinical diagnosis of insomnia, yielding a beneficial impact on the identification of effective biomarkers and the comprehension of insomnia.

Hierarchical gradients of multiple timescales in the mammalian forebrain

Many anatomical and physiological features of cortical circuits, ranging from the biophysical properties of synapses to the connectivity patterns among different neuron types, exhibit consistent variation along the hierarchical axis from sensory to association areas. Notably, the temporal correlation of neural activity at rest, known as the intrinsic timescale, increases systematically along this hierarchy in both primates and rodents, analogous to the increasing scale and complexity of spatial receptive fields. However, how the timescales for task-related activity vary across brain regions and whether their hierarchical organization appears consistently across different mammalian species remain unexplored. Here, we show that both the intrinsic timescale and those of task-related activity follow a similar hierarchical gradient in the cortices of monkeys, rats, and mice. We also found that these timescales covary similarly in both the cortex and basal ganglia, whereas the timescales of thalamic activity are shorter than cortical timescales and do not conform to the hierarchical order predicted by their cortical projections. These results suggest that the hierarchical gradient of cortical timescales might represent a universal feature of intracortical circuits in the mammalian brain.

Cognitive network interactions through communication subspaces in large-scale models of the neocortex

Neocortex-wide neural activity is organized into distinct networks of areas engaged in different cognitive processes. To elucidate the underlying mechanism of flexible network reconfiguration, we developed connectivity-constrained macaque and human whole-cortex models. In our model, within-area connectivity consists of a mixture of symmetric, asymmetric, and random motifs that give rise to stable (attractor) or transient (sequential) heterogeneous dynamics. Assuming sparse low-rank plus random inter-areal connectivity constrained by cognitive networks' activation maps, we show that our model captures key aspects of the cognitive networks' dynamics and interactions observed experimentally. In particular, the anti-correlation between the default mode network and the dorsal attention network. Communication between networks is shaped by the alignment of long-range communication subspaces with local connectivity motifs and is switchable in a bottom-up salience-dependent routing mechanism. Furthermore, the frontoparietal multiple-demand network displays a coexistence of stable and dynamic coding, suitable for top-down cognitive control. Our work provides a theoretical framework for understanding the dynamic routing in the cortical networks during cognition.

Neural mass modeling for the masses: Democratizing access to whole-brain biophysical modeling with FastDMF

Different whole-brain computational models have been recently developed to investigate hypotheses related to brain mechanisms. Among these, the Dynamic Mean Field (DMF) model is particularly attractive, combining a biophysically realistic model that is scaled up via a mean-field approach and multimodal imaging data. However, an important barrier to the widespread usage of the DMF model is that current implementations are computationally expensive, supporting only simulations on brain parcellations that consider less than 100 brain regions. Here, we introduce an efficient and accessible implementation of the DMF model: the FastDMF. By leveraging analytical and numerical advances-including a novel estimation of the feedback inhibition control parameter and a Bayesian optimization algorithm-the FastDMF circumvents various computational bottlenecks of previous implementations, improving interpretability, performance, and memory use. Furthermore, these advances allow the FastDMF to increase the number of simulated regions by one order of magnitude, as confirmed by the good fit to fMRI data parcellated at 90 and 1,000 regions. These advances open the way to the widespread use of biophysically grounded whole-brain models for investigating the interplay between anatomy, function, and brain dynamics and to identify mechanistic explanations of recent results obtained from fine-grained neuroimaging recordings.

Decoding cortical chronotopy-Comparing the influence of different cortical organizational schemes

The brain's diverse intrinsic timescales enable us to perceive stimuli with varying temporal persistency. This study aimed to uncover the cortical organizational schemes underlying these variations, revealing the neural architecture for processing a wide range of sensory experiences. We collected resting-state fMRI, task-fMRI, and diffusion-weighted imaging data from 47 individuals. Based on this data, we extracted six organizational schemes: (1) the structural Rich Club (RC) architecture, shown to synchronize the connectome; (2) the structural Diverse Club architecture, as an alternative to the RC based on the network's module structure; (3) the functional uni-to-multimodal gradient, reflected in a wide range of structural and functional features; and (4) the spatial posterior/lateral-to-anterior/medial gradient, established for hierarchical levels of cognitive control. Also, we explored the effects of (5) structural graph theoretical measures of centrality and (6) cytoarchitectural differences. Using Bayesian model comparison, we contrasted the impact of these organizational schemes on (1) intrinsic resting-state timescales and (2) inter-subject correlation (ISC) from a task involving hierarchically nested digit sequences. As expected, resting-state timescales were slower in structural network hubs, hierarchically higher areas defined by the functional and spatial gradients, and thicker cortical regions. ISC analysis demonstrated hints for the engagement of higher cortical areas with more temporally persistent stimuli. Finally, the model comparison identified the uni-to-multimodal gradient as the best organizational scheme for explaining the chronotopy in both task and rest. Future research should explore the microarchitectural features that shape this gradient, elucidating how our brain adapts and evolves across different modes of processing.

Whole brain functional connectivity: Insights from next generation neural mass modelling incorporating electrical synapses

The ready availability of brain connectome data has both inspired and facilitated the modelling of whole brain activity using networks of phenomenological neural mass models that can incorporate both interaction strength and tract length between brain regions. Recently, a new class of neural mass model has been developed from an exact mean field reduction of a network of spiking cortical cell models with a biophysically realistic model of the chemical synapse. Moreover, this new population dynamics model can naturally incorporate electrical synapses. Here we demonstrate the ability of this new modelling framework, when combined with data from the Human Connectome Project, to generate patterns of functional connectivity (FC) of the type observed in both magnetoencephalography and functional magnetic resonance neuroimaging. Some limited explanatory power is obtained via an eigenmode description of frequency-specific FC patterns, obtained via a linear stability analysis of the network steady state in the neigbourhood of a Hopf bifurcation. However, direct numerical simulations show that empirical data is more faithfully recapitulated in the nonlinear regime, and exposes a key role of gap junction coupling strength in generating empirically-observed neural activity, and associated FC patterns and their evolution. Thereby, we emphasise the importance of maintaining known links with biological reality when developing multi-scale models of brain dynamics. As a tool for the study of dynamic whole brain models of the type presented here we further provide a suite of C++ codes for the efficient, and user friendly, simulation of neural mass networks with multiple delayed interactions.

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.

Hierarchical individual variation and socioeconomic impact on personalized functional network topography in children

BACKGROUND

The spatial layout of large-scale functional brain networks exhibits considerable inter-individual variability, especially in the association cortex. Research has demonstrated a link between early socioeconomic status (SES) and variations in both brain structure and function, which are further associated with cognitive and mental health outcomes. However, the extent to which SES is associated with individual differences in personalized functional network topography during childhood remains largely unexplored.

METHODS

We used a machine learning approach-spatially regularized non-negative matrix factorization (NMF)-to delineate 17 personalized functional networks in children aged 9-10 years, utilizing high-quality functional MRI data from 6001 participants in the Adolescent Brain Cognitive Development study. Partial least square regression approach with repeated random twofold cross-validation was used to evaluate the association between the multivariate pattern of functional network topography and three SES factors, including family income-to-needs ratio, parental education, and neighborhood disadvantage.

RESULTS

We found that individual variations in personalized functional network topography aligned with the hierarchical sensorimotor-association axis across the cortex. Furthermore, we observed that functional network topography significantly predicted the three SES factors from unseen individuals. The associations between functional topography and SES factors were also hierarchically organized along the sensorimotor-association cortical axis, exhibiting stronger positive associations in the higher-order association cortex. Additionally, we have made the personalized functional networks publicly accessible.

CONCLUSIONS

These results offer insights into how SES influences neurodevelopment through personalized functional neuroanatomy in childhood, highlighting the cortex-wide, hierarchically organized plasticity of the functional networks in response to diverse SES backgrounds.

Comparing structure-function relationships in brain networks using EEG and fNIRS

Identifying relationships between structural and functional networks is crucial for understanding the large-scale organization of the human brain. The potential contribution of emerging techniques like functional near-infrared spectroscopy to investigate the structure-functional relationship has yet to be explored. In our study, using simultaneous Electroencephalography (EEG) and Functional near-infrared spectroscopy (fNIRS) recordings from 18 subjects, we characterize global and local structure-function coupling using source-reconstructed EEG and fNIRS signals in both resting state and motor imagery tasks, as this relationship during task periods remains underexplored. Employing the mathematical framework of graph signal processing, we investigate how this relationship varies across electrical and hemodynamic networks and different brain states. Results show that fNIRS structure-function coupling resembles slower-frequency EEG coupling at rest, with variations across brain states and oscillations. Locally, the relationship is heterogeneous, with greater coupling in the sensory cortex and increased decoupling in the association cortex, following the unimodal to transmodal gradient. Discrepancies between EEG and fNIRS are noted, particularly in the frontoparietal network. Cross-band representations of neural activity revealed lower correspondence between electrical and hemodynamic activity in the transmodal cortex, irrespective of brain state while showing specificity for the somatomotor network during a motor imagery task. Overall, these findings initiate a multimodal comprehension of structure-function relationship and brain organization when using affordable functional brain imaging.

The localization of coma

Coma and disorders of consciousness (DoC) are common manifestations of acute severe brain injuries. Research into their neuroanatomical basis can be traced from Hippocrates to the present day. Lesions causing DoC have traditionally been conceptualized as decreasing "alertness" from damage to the ascending arousal system, and/or, reducing level of "awareness" due to structural or functional impairment of large-scale brain networks. Within this framework, pharmacological and neuromodulatory interventions to promote recovery from DoC have hitherto met with limited success. This is partly due to inter-individual heterogeneity of brain injury patterns, and an incomplete understanding of brain network properties that characterize consciousness. Advances in multiscale computational modelling of brain dynamics have opened a unique opportunity to explore the causal mechanisms of brain activity at the biophysical level. These models can provide a novel approach for selection and optimization of potential interventions by simulation of brain network dynamics individualized for each patient.

Comparison of whole-brain task-modulated functional connectivity methods for fMRI task connectomics

Higher brain functions require flexible integration of information across widely distributed brain regions depending on the task context. Resting-state functional magnetic resonance imaging (fMRI) has provided substantial insight into large-scale intrinsic brain network organisation, yet the principles of rapid context-dependent reconfiguration of that intrinsic network organisation are much less understood. A major challenge for task connectome mapping is the absence of a gold standard for deriving whole-brain task-modulated functional connectivity matrices. Here, we perform biophysically realistic simulations to control the ground-truth task-modulated functional connectivity over a wide range of experimental settings. We reveal the best-performing methods for different types of task designs and their fundamental limitations. Importantly, we demonstrate that rapid (100 ms) modulations of oscillatory neuronal synchronisation can be recovered from sluggish haemodynamic fluctuations even at typically low fMRI temporal resolution (2 s). Finally, we provide practical recommendations on task design and statistical analysis to foster task connectome mapping.

Competitive interactions shape brain dynamics and computation across species

Adaptive cognition relies on cooperation across anatomically distributed brain circuits. However, specialised neural systems are also in constant competition for limited processing resources. How does the brain's network architecture enable it to balance these cooperative and competitive tendencies? Here we use computational whole-brain modelling to examine the dynamical and computational relevance of cooperative and competitive interactions in the mammalian connectome. Across human, macaque, and mouse we show that the architecture of the models that most faithfully reproduce brain activity, consistently combines modular cooperative interactions with diffuse, long-range competitive interactions. The model with competitive interactions consistently outperforms the cooperative-only model, with excellent fit to both spatial and dynamical properties of the living brain, which were not explicitly optimised but rather emerge spontaneously. Competitive interactions in the effective connectivity produce greater levels of synergistic information and local-global hierarchy, and lead to superior computational capacity when used for neuromorphic computing. Altogether, this work provides a mechanistic link between network architecture, dynamical properties, and computation in the mammalian brain.

Comparison of whole-brain task-modulated functional connectivity methods for fMRI task connectomics

Higher brain functions require flexible integration of information across widely distributed brain regions depending on the task context. Resting-state functional magnetic resonance imaging (fMRI) has provided substantial insight into large-scale intrinsic brain network organisation, yet the principles of rapid context-dependent reconfiguration of that intrinsic network organisation are much less understood. A major challenge for task connectome mapping is the absence of a gold standard for deriving whole-brain task-modulated functional connectivity matrices. Here, we perform biophysically realistic simulations to control the ground-truth task-modulated functional connectivity over a wide range of experimental settings. We reveal the best-performing methods for different types of task designs and their fundamental limitations. Importantly, we demonstrate that rapid (100 ms) modulations of oscillatory neuronal synchronisation can be recovered from sluggish haemodynamic fluctuations even at typically low fMRI temporal resolution (2 s). Finally, we provide practical recommendations on task design and statistical analysis to foster task connectome mapping.

Macroscale connectome topographical structure reveals the biomechanisms of brain dysfunction in Alzheimer's disease

The intricate spatial configurations of brain networks offer essential insights into understanding the specific patterns of brain abnormalities and the underlying biological mechanisms associated with Alzheimer's disease (AD), normal aging, and other neurodegenerative disorders. This study investigated alterations in the topographical structure of the brain related to aging and neurodegenerative diseases by analyzing brain gradients derived from structural MRI data across multiple cohorts (n = 7323). The analysis identified distinct gradient patterns in AD, aging, and other neurodegenerative conditions. Gene enrichment analysis indicated that inorganic ion transmembrane transport was the most significant term in normal aging, while chemical synaptic transmission is a common enrichment term across various neurodegenerative diseases. Moreover, the findings show that each disorder exhibits unique dysfunctional neurophysiological characteristics. These insights are pivotal for elucidating the distinct biological mechanisms underlying AD, thereby enhancing our understanding of its unique clinical phenotypes in contrast to normal aging and other neurodegenerative disorders.

Large-Scale Mechanistic Models of Brain Circuits with Biophysically and Morphologically Detailed Neurons

Understanding the brain requires studying its multiscale interactions from molecules to networks. The increasing availability of large-scale datasets detailing brain circuit composition, connectivity, and activity is transforming neuroscience. However, integrating and interpreting this data remains challenging. Concurrently, advances in supercomputing and sophisticated modeling tools now enable the development of highly detailed, large-scale biophysical circuit models. These mechanistic multiscale models offer a method to systematically integrate experimental data, facilitating investigations into brain structure, function, and disease. This review, based on a Society for Neuroscience 2024 MiniSymposium, aims to disseminate recent advances in large-scale mechanistic modeling to the broader community. It highlights (1) examples of current models for various brain regions developed through experimental data integration; (2) their predictive capabilities regarding cellular and circuit mechanisms underlying experimental recordings (e.g., membrane voltage, spikes, local-field potential, electroencephalography/magnetoencephalography) and brain function; and (3) their use in simulating biomarkers for brain diseases like epilepsy, depression, schizophrenia, and Parkinson's, aiding in understanding their biophysical underpinnings and developing novel treatments. The review showcases state-of-the-art models covering hippocampus, somatosensory, visual, motor, auditory cortical, and thalamic circuits across species. These models predict neural activity at multiple scales and provide insights into the biophysical mechanisms underlying sensation, motor behavior, brain signals, neural coding, disease, pharmacological interventions, and neural stimulation. Collaboration with experimental neuroscientists and clinicians is essential for the development and validation of these models, particularly as datasets grow. Hence, this review aims to foster interest in detailed brain circuit models, leading to cross-disciplinary collaborations that accelerate brain research.

Temporal neural dynamics of understanding communicative intentions from speech prosody

Understanding the correct intention of a speaker is critical for social interaction. Speech prosody is an important source for understanding speakers' intentions during verbal communication. However, the neural dynamics by which the human brain translates the prosodic cues into a mental representation of communicative intentions in real time remains unclear. Here, we recorded EEG (electroencephalograph) while participants listened to dialogues. The prosodic features of the critical words at the end of sentences were manipulated to signal either suggestion, warning, or neutral intentions. The results showed that suggestion and warning intentions evoked enhanced late positive event-related potentials (ERPs) compared to the neutral condition. Linear mixed-effects model (LMEM) regression and representational similarity analysis (RSA) analyses revealed that these ERP effects were distinctively correlated with prosodic acoustic analysis, emotional valence evaluation, and intention interpretation in different time windows; The onset latency significantly increased as the processing level of abstractness and communicative intentionality increased. Neural representations of intention and emotional information emerged and parallelly persisted over a long time window, guiding the correct identification of communicative intention. These results provide new insights into understanding the structural components of intention processing and their temporal neural dynamics underlying communicative intention comprehension from speech prosody in online social interactions.

Analyzing asymmetry in brain hierarchies with a linear state-space model of resting-state fMRI data

This study challenges the traditional focus on zero-lag statistics in resting-state functional magnetic resonance imaging (rsfMRI) research. Instead, it advocates for considering time-lag interactions to unveil the directionality and asymmetries of the brain hierarchy. Effective connectivity (EC), the state matrix in dynamical causal modeling (DCM), is a commonly used metric for studying dynamical properties and causal interactions within a linear state-space system description. Here, we focused on how time-lag statistics are incorporated within the framework of DCM resulting in an asymmetric EC matrix. Our approach involves decomposing the EC matrix, revealing a steady-state differential cross-covariance matrix that is responsible for modeling information flow and introducing time-irreversibility. Specifically, the system's dynamics, influenced by the off-diagonal part of the differential covariance, exhibit a curl steady-state flow component that breaks detailed balance and diverges the dynamics from equilibrium. Our empirical findings indicate that the EC matrix's outgoing strengths correlate with the flow described by the differential cross covariance, while incoming strengths are primarily driven by zero-lag covariance, emphasizing conditional independence over directionality.

Group-common and individual-specific effects of structure-function coupling in human brain networks with graph neural networks

The human cerebral cortex is organized into functionally segregated but synchronized regions bridged by the structural connectivity of white matter pathways. While structure-function coupling has been implicated in cognitive development and neuropsychiatric disorders, it remains unclear to what extent the structure-function coupling reflects a group-common characteristic or varies across individuals, at both the global and regional brain levels. By leveraging two independent, high-quality datasets, we found that the graph neural network accurately predicted unseen individuals' functional connectivity from structural connectivity, reflecting a strong structure-function coupling. This coupling was primarily driven by network topology and was substantially stronger than that of the correlation approaches. Moreover, we observed that structure-function coupling was dominated by group-common effects, with subtle yet significant individual-specific effects. The regional group and individual effects of coupling were hierarchically organized across the cortex along a sensorimotor-association axis, with lower group and higher individual effects in association cortices. These findings emphasize the importance of considering both group and individual effects in understanding cortical structure-function coupling, suggesting insights into interpreting individual differences of the coupling and informing connectivity-guided therapeutics.

Dendrites contribute to the gradient of intrinsic timescales encompassing cortical and subcortical brain networks

Introduction

Cytoarchitectonic studies have uncovered a correlation between higher levels of cortical hierarchy and reduced dendritic size. This hierarchical organization extends to the brain's timescales, revealing longer intrinsic timescales at higher hierarchical levels. However, estimating the contribution of single-neuron dendritic morphology to the hierarchy of timescales, which is typically characterized at a macroscopic level, remains challenging.

Method

Here we mapped the intrinsic timescales of six functional networks using functional magnetic resonance imaging (fMRI) data, and characterized the influence of neuronal dendritic size on intrinsic timescales of brain regions, utilizing a multicompartmental neuronal modeling approach based on digitally reconstructed neurons.

Results

The fMRI results revealed a hierarchy of intrinsic timescales encompassing both cortical and subcortical brain regions. The neuronal modeling indicated that neurons with larger dendritic structures exhibit shorter intrinsic timescales. Together these findings highlight the contribution of dendrites at the neuronal level to the hierarchy of intrinsic timescales at the whole-brain level.

Discussion

This study sheds light on the intricate relationship between neuronal structure, cytoarchitectonic maps, and the hierarchy of timescales in the brain.

Hierarchical gradients of multiple timescales in the mammalian forebrain

Many anatomical and physiological features of cortical circuits, ranging from the biophysical properties of synapses to the connectivity patterns among different neuron types, exhibit consistent variation along the hierarchical axis from sensory to association areas. Notably, the scale of temporal correlation of neural activity at rest, known as the intrinsic timescale, increases systematically along this hierarchy in both primates and rodents, analogous to the growing scale and complexity of spatial receptive fields. However, how the timescales for task-related activity vary across brain regions and whether their hierarchical organization appears consistently across different mammalian species remain unexplored. Here, we show that both the intrinsic timescale and the timescales of task-related activity follow a similar hierarchical gradient in the cortices of monkeys, rats, and mice. We also found that these timescales covary similarly in both the cortex and basal ganglia, whereas the timescales of thalamic activity are shorter than cortical timescales and do not conform to the hierarchical order predicted by their cortical projections. These results suggest that the hierarchical gradient of cortical timescales might be a universal feature of intra-cortical circuits in the mammalian brain.

A synergistic workspace for human consciousness revealed by Integrated Information Decomposition

How is the information-processing architecture of the human brain organised, and how does its organisation support consciousness? Here, we combine network science and a rigorous information-theoretic notion of synergy to delineate a 'synergistic global workspace', comprising gateway regions that gather synergistic information from specialised modules across the human brain. This information is then integrated within the workspace and widely distributed via broadcaster regions. Through functional MRI analysis, we show that gateway regions of the synergistic workspace correspond to the human brain's default mode network, whereas broadcasters coincide with the executive control network. We find that loss of consciousness due to general anaesthesia or disorders of consciousness corresponds to diminished ability of the synergistic workspace to integrate information, which is restored upon recovery. Thus, loss of consciousness coincides with a breakdown of information integration within the synergistic workspace of the human brain. This work contributes to conceptual and empirical reconciliation between two prominent scientific theories of consciousness, the Global Neuronal Workspace and Integrated Information Theory, while also advancing our understanding of how the human brain supports consciousness through the synergistic integration of information.

Unravelling consciousness and brain function through the lens of time, space, and information

Disentangling how cognitive functions emerge from the interplay of brain dynamics and network architecture is among the major challenges that neuroscientists face. Pharmacological and pathological perturbations of consciousness provide a lens to investigate these complex challenges. Here, we review how recent advances about consciousness and the brain's functional organisation have been driven by a common denominator: decomposing brain function into fundamental constituents of time, space, and information. Whereas unconsciousness increases structure-function coupling across scales, psychedelics may decouple brain function from structure. Convergent effects also emerge: anaesthetics, psychedelics, and disorders of consciousness can exhibit similar reconfigurations of the brain's unimodal-transmodal functional axis. Decomposition approaches reveal the potential to translate discoveries across species, with computational modelling providing a path towards mechanistic integration.

Interacting ventral temporal gradients of timescales and functional connectivity and their relationships to visual behavior

Cortical gradients in endogenous and stimulus-evoked neurodynamic timescales, and long-range cortical interactions, provide organizational constraints to the brain and influence neural populations' roles in cognition. It is unclear how these functional gradients interrelate and which influence behavior. Here, intracranial recordings from 4,090 electrode contacts in 35 individuals map gradients of neural timescales and functional connectivity to assess their interactions along category-selective ventral temporal cortex. Endogenous and stimulus-evoked information processing timescales were not significantly correlated with one another suggesting that local neural timescales are context dependent and may arise through distinct neurophysiological mechanisms. Endogenous neural timescales correlated with functional connectivity even after removing the effects of shared anatomical gradients. Neural timescales and functional connectivity correlated with how strongly a population's activity predicted behavior in a simple visual task. These results suggest both interrelated and distinct neurophysiological processes give rise to different functional connectivity and neural timescale gradients, which together influence behavior.

Spatial and temporal pattern of structure-function coupling of human brain connectome with development

Brain structural circuitry shapes a richly patterned functional synchronization, supporting for complex cognitive and behavioural abilities. However, how coupling of structural connectome (SC) and functional connectome (FC) develops and its relationships with cognitive functions and transcriptomic architecture remain unclear. We used multimodal magnetic resonance imaging data from 439 participants aged 5.7-21.9 years to predict functional connectivity by incorporating intracortical and extracortical structural connectivity, characterizing SC-FC coupling. Our findings revealed that SC-FC coupling was strongest in the visual and somatomotor networks, consistent with evolutionary expansion, myelin content, and functional principal gradient. As development progressed, SC-FC coupling exhibited heterogeneous alterations dominated by an increase in cortical regions, broadly distributed across the somatomotor, frontoparietal, dorsal attention, and default mode networks. Moreover, we discovered that SC-FC coupling significantly predicted individual variability in general intelligence, mainly influencing frontoparietal and default mode networks. Finally, our results demonstrated that the heterogeneous development of SC-FC coupling is positively associated with genes in oligodendrocyte-related pathways and negatively associated with astrocyte-related genes. This study offers insight into the maturational principles of SC-FC coupling in typical development.

Adolescent maturation of cortical excitation-inhibition balance based on individualized biophysical network modeling

The balance of excitation and inhibition is a key functional property of cortical microcircuits which changes through the lifespan. Adolescence is considered a crucial period for the maturation of excitation-inhibition balance. This has been primarily observed in animal studies, yet human in vivo evidence on adolescent maturation of the excitation-inhibition balance at the individual level is limited. Here, we developed an individualized in vivo marker of regional excitation-inhibition balance in human adolescents, estimated using large-scale simulations of biophysical network models fitted to resting-state functional magnetic resonance imaging data from two independent cross-sectional (N = 752) and longitudinal (N = 149) cohorts. We found a widespread relative increase of inhibition in association cortices paralleled by a relative age-related increase of excitation, or lack of change, in sensorimotor areas across both datasets. This developmental pattern co-aligned with multiscale markers of sensorimotor-association differentiation. The spatial pattern of excitation-inhibition development in adolescence was robust to inter-individual variability of structural connectomes and modeling configurations. Notably, we found that alternative simulation-based markers of excitation-inhibition balance show a variable sensitivity to maturational change. Taken together, our study highlights an increase of inhibition during adolescence in association areas using cross sectional and longitudinal data, and provides a robust computational framework to estimate microcircuit maturation in vivo at the individual level.

In vivo whole-cortex marker of excitation-inhibition ratio indexes cortical maturation and cognitive ability in youth

A balanced excitation-inhibition ratio (E/I ratio) is critical for healthy brain function. Normative development of cortex-wide E/I ratio remains unknown. Here, we noninvasively estimate a putative marker of whole-cortex E/I ratio by fitting a large-scale biophysically plausible circuit model to resting-state functional MRI (fMRI) data. We first confirm that our model generates realistic brain dynamics in the Human Connectome Project. Next, we show that the estimated E/I ratio marker is sensitive to the gamma-aminobutyric acid (GABA) agonist benzodiazepine alprazolam during fMRI. Alprazolam-induced E/I changes are spatially consistent with positron emission tomography measurement of benzodiazepine receptor density. We then investigate the relationship between the E/I ratio marker and neurodevelopment. We find that the E/I ratio marker declines heterogeneously across the cerebral cortex during youth, with the greatest reduction occurring in sensorimotor systems relative to association systems. Importantly, among children with the same chronological age, a lower E/I ratio marker (especially in the association cortex) is linked to better cognitive performance. This result is replicated across North American (8.2 to 23.0 y old) and Asian (7.2 to 7.9 y old) cohorts, suggesting that a more mature E/I ratio indexes improved cognition during normative development. Overall, our findings open the door to studying how disrupted E/I trajectories may lead to cognitive dysfunction in psychopathology that emerges during youth.

Sustained attention operates via dissociable neural mechanisms across different eccentric locations

In primates, foveal and peripheral vision have distinct neural architectures and functions. However, it has been debated if selective attention operates via the same or different neural mechanisms across eccentricities. We tested these alternative accounts by examining the effects of selective attention on the steady-state visually evoked potential (SSVEP) and the fronto-parietal signal measured via EEG from human subjects performing a sustained visuospatial attention task. With a negligible level of eye movements, both SSVEP and SND exhibited the heterogeneous patterns of attentional modulations across eccentricities. Specifically, the attentional modulations of these signals peaked at the parafoveal locations and such modulations wore off as visual stimuli appeared closer to the fovea or further away towards the periphery. However, with a relatively higher level of eye movements, the heterogeneous patterns of attentional modulations of these neural signals were less robust. These data demonstrate that the top-down influence of covert visuospatial attention on early sensory processing in human cortex depends on eccentricity and the level of saccadic responses. Taken together, the results suggest that sustained visuospatial attention operates differently across different eccentric locations, providing new understanding of how attention augments sensory representations regardless of where the attended stimulus appears.

Body mass index-dependent shifts along large-scale gradients in human cortical organization explain dietary regulatory success

Making healthy dietary choices is essential for keeping weight within a normal range. Yet many people struggle with dietary self-control despite good intentions. What distinguishes neural processing in those who succeed or fail to implement healthy eating goals? Does this vary by weight status? To examine these questions, we utilized an analytical framework of gradients that characterize systematic spatial patterns of large-scale neural activity, which have the advantage of considering the entire suite of processes subserving self-control and potential regulatory tactics at the whole-brain level. Using an established laboratory food task capturing brain responses in natural and regulatory conditions (N = 123), we demonstrate that regulatory changes of dietary brain states in the gradient space predict individual differences in dietary success. Better regulators required smaller shifts in brain states to achieve larger goal-consistent changes in dietary behaviors, pointing toward efficient network organization. This pattern was most pronounced in individuals with lower weight status (low-BMI, body mass index) but absent in high-BMI individuals. Consistent with prior work, regulatory goals increased activity in frontoparietal brain circuits. However, this shift in brain states alone did not predict variance in dietary success. Instead, regulatory success emerged from combined changes along multiple gradients, showcasing the interplay of different large-scale brain networks subserving dietary control and possible regulatory strategies. Our results provide insights into how the brain might solve the problem of dietary control: Dietary success may be easier for people who adopt modes of large-scale brain activation that do not require significant reconfigurations across contexts and goals.

Cell-type-specific optogenetic fMRI on basal forebrain reveals functional network basis of behavioral preference

The basal forebrain (BF) is a complex structure that plays key roles in regulating various brain functions. However, it remains unclear how cholinergic and non-cholinergic BF neurons modulate large-scale functional networks and their relevance in intrinsic and extrinsic behaviors. With an optimized awake mouse optogenetic fMRI approach, we revealed that optogenetic stimulation of four BF neuron types evoked distinct cell-type-specific whole-brain BOLD activations, which could be attributed to BF-originated low-dimensional structural networks. Additionally, optogenetic activation of VGLUT2, ChAT, and PV neurons in the BF modulated the preference for locomotion, exploration, and grooming, respectively. Furthermore, we uncovered the functional network basis of the above BF-modulated behavioral preference through a decoding model linking the BF-modulated BOLD activation, low-dimensional structural networks, and behavioral preference. To summarize, we decoded the functional network basis of differential behavioral preferences with cell-type-specific optogenetic fMRI on the BF and provided an avenue for investigating mouse behaviors from a whole-brain view.

Information decomposition and the informational architecture of the brain

To explain how the brain orchestrates information-processing for cognition, we must understand information itself. Importantly, information is not a monolithic entity. Information decomposition techniques provide a way to split information into its constituent elements: unique, redundant, and synergistic information. We review how disentangling synergistic and redundant interactions is redefining our understanding of integrative brain function and its neural organisation. To explain how the brain navigates the trade-offs between redundancy and synergy, we review converging evidence integrating the structural, molecular, and functional underpinnings of synergy and redundancy; their roles in cognition and computation; and how they might arise over evolution and development. Overall, disentangling synergistic and redundant information provides a guiding principle for understanding the informational architecture of the brain and cognition.

In-vivo whole-cortex marker of excitation-inhibition ratio indexes cortical maturation and cognitive ability in youth

A balanced excitation-inhibition ratio (E/I ratio) is critical for healthy brain function. Normative development of cortex-wide E/I ratio remains unknown. Here we non-invasively estimate a putative marker of whole-cortex E/I ratio by fitting a large-scale biophysically-plausible circuit model to resting-state functional MRI (fMRI) data. We first confirm that our model generates realistic brain dynamics in the Human Connectome Project. Next, we show that the estimated E/I ratio marker is sensitive to the GABA-agonist benzodiazepine alprazolam during fMRI. Alprazolam-induced E/I changes are spatially consistent with positron emission tomography measurement of benzodiazepine receptor density. We then investigate the relationship between the E/I ratio marker and neurodevelopment. We find that the E/I ratio marker declines heterogeneously across the cerebral cortex during youth, with the greatest reduction occurring in sensorimotor systems relative to association systems. Importantly, among children with the same chronological age, a lower E/I ratio marker (especially in association cortex) is linked to better cognitive performance. This result is replicated across North American (8.2 to 23.0 years old) and Asian (7.2 to 7.9 years old) cohorts, suggesting that a more mature E/I ratio indexes improved cognition during normative development. Overall, our findings open the door to studying how disrupted E/I trajectories may lead to cognitive dysfunction in psychopathology that emerges during youth.

Spatial and temporal pattern of structure-function coupling of human brain connectome with development

Brain structural circuitry shapes a richly patterned functional synchronization, supporting for complex cognitive and behavioural abilities. However, how coupling of structural connectome (SC) and functional connectome (FC) develops and its relationships with cognitive functions and transcriptomic architecture remain unclear. We used multimodal magnetic resonance imaging data from 439 participants aged 5.7 to 21.9 years to predict functional connectivity by incorporating intracortical and extracortical structural connectivity, characterizing SC-FC coupling. Our findings revealed that SC-FC coupling was strongest in the visual and somatomotor networks, consistent with evolutionary expansion, myelin content, and functional principal gradient. As development progressed, SC-FC coupling exhibited heterogeneous alterations dominated by an increase in cortical regions, broadly distributed across the somatomotor, frontoparietal, dorsal attention, and default mode networks. Moreover, we discovered that SC-FC coupling significantly predicted individual variability in general intelligence, mainly influencing frontoparietal and default mode networks. Finally, our results demonstrated that the heterogeneous development of SC-FC coupling is positively associated with genes in oligodendrocyte-related pathways and negatively associated with astrocyte-related genes. This study offers insight into the maturational principles of SC-FC coupling in typical development.

Where do the symptoms come from in depression? Topography and dynamics matter

This scientific commentary refers to 'Brain dynamics predictive of response to psilocybin for treatment-resistant depression', by Vohryzek et al. (https://doi.org/10.1093/braincomms/fcae049).

Decreased intrinsic neural timescale in treatment-naïve adolescent depression

BACKGROUND

Major depressive disorder (MDD) is mainly characterized by its core dysfunction in higher-order brain cortices involved in emotional and cognitive processes, whose neurobiological basis remains unclear. In this study, we applied a relatively new developed resting-state functional magnetic resonance imaging (rs-fMRI) method of intrinsic neural timescale (INT), which reflects how long neural information is stored in a local brain area and reflects an ability of information integration, to investigate the local intrinsic neural dynamics using univariate and multivariate analyses in adolescent depression.

METHOD

Based on the rs-fMRI data of sixty-six treatment-naïve adolescents with MDD and fifty-two well-matched healthy controls (HCs), we calculated an INT by assessing the magnitude of autocorrelation of the resting-state brain activity, and then compared the difference of INT between the two groups. Correlation between abnormal INT and clinical features was performed. We also utilized multivariate pattern analysis to determine whether INT could differentiate MDD patients from HCs at the individual level.

RESULT

Compared with HCs, patients with MDD showed shorter INT widely distributed in cortical and partial subcortical regions. Interestingly, the decreased INT in the left hippocampus was related to disease severity of MDD. Furthermore, INT can distinguish MDD patients from HCs with the most discriminative regions located in the dorsolateral prefrontal cortex, angular, middle occipital gyrus, and cerebellar posterior lobe.

CONCLUSION

Our research aids in advancing understanding the brain abnormalities of treatment-naïve adolescents with MDD from the perspective of the local neural dynamics, highlighting the significant role of INT in understanding neurophysiological mechanisms. This study shows that the altered intrinsic timescales of local neural signals widely distributed in higher-order brain cortices regions may be the neurodynamic basis of cognitive and emotional disturbances in MDD patients, and provides preliminary support for the suggestion that these could be used to aid the identification of MDD patients in clinical practice.

Dynamic nonreversibility view of intrinsic brain organization and brain dynamic analysis of repetitive transcranial magnitude stimulation

Intrinsic neural activities are characterized as endless spontaneous fluctuation over multiple time scales. However, how the intrinsic brain organization changes over time under local perturbation remains an open question. By means of statistical physics, we proposed an approach to capture whole-brain dynamics based on estimating time-varying nonreversibility and k-means clustering of dynamic varying nonreversibility patterns. We first used synthetic fMRI to investigate the effects of window parameters on the temporal variability of varying nonreversibility. Second, using real test-retest fMRI data, we examined the reproducibility, reliability, biological, and physiological correlation of the varying nonreversibility substates. Finally, using repetitive transcranial magnetic stimulation-fMRI data, we investigated the modulation effects of repetitive transcranial magnetic stimulation on varying nonreversibility substate dynamics. The results show that: (i) as window length increased, the varying nonreversibility variance decreased, while the sliding step almost did not alter it; (ii) the global high varying nonreversibility states and low varying nonreversibility states were reproducible across multiple datasets and different window lengths; and (iii) there were increased low varying nonreversibility states and decreased high varying nonreversibility states when the left frontal lobe was stimulated, but not the occipital lobe. Taken together, these results provide a thermodynamic equilibrium perspective of intrinsic brain organization and reorganization under local perturbation.

Exploring global and local processes underlying alterations in resting-state functional connectivity and dynamics in schizophrenia

Introduction

We examined changes in large-scale functional connectivity and temporal dynamics and their underlying mechanisms in schizophrenia (ScZ) through measurements of resting-state functional magnetic resonance imaging (rs-fMRI) data and computational modelling.

Methods

The rs-fMRI measurements from patients with chronic ScZ (n=38) and matched healthy controls (n=43), were obtained through the public schizConnect repository. Computational models were constructed based on diffusion-weighted MRI scans and fit to the experimental rs-fMRI data.

Results

We found decreased large-scale functional connectivity across sensory and association areas and for all functional subnetworks for the ScZ group. Additionally global synchrony was reduced in patients while metastability was unaltered. Perturbations of the computational model revealed that decreased global coupling and increased background noise levels both explained the experimentally found deficits better than local changes to the GABAergic or glutamatergic system.

Discussion

The current study suggests that large-scale alterations in ScZ are more likely the result of global rather than local network changes.