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.

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.

Macroscale intrinsic dynamics are associated with microcircuit function in focal and generalized epilepsies.

Epilepsies are a group of neurological disorders characterized by abnormal spontaneous brain activity, involving multiscale changes in brain functional organizations. However, it is not clear to what extent the epilepsy-related perturbations of spontaneous brain activity affect macroscale intrinsic dynamics and microcircuit organizations, that supports their pathological relevance. We collect a sample of patients with temporal lobe epilepsy (TLE) and genetic generalized epilepsy with tonic-clonic seizure (GTCS), as well as healthy controls. We extract massive temporal features of fMRI BOLD time-series to characterize macroscale intrinsic dynamics, and simulate microcircuit neuronal dynamics used a large-scale biological model. Here we show whether macroscale intrinsic dynamics and microcircuit dysfunction are differed in epilepsies, and how these changes are linked. Differences in macroscale gradient of time-series features are prominent in the primary network and default mode network in TLE and GTCS. Biophysical simulations indicate reduced recurrent connection within somatomotor microcircuits in both subtypes, and even more reduced in GTCS. We further demonstrate strong spatial correlations between differences in the gradient of macroscale intrinsic dynamics and microcircuit dysfunction in epilepsies. These results emphasize the impact of abnormal neuronal activity on primary network and high-order networks, suggesting a systematic abnormality of brain hierarchical organization.

Dynamic primitives of brain network interaction.

What dynamic processes underly functional brain networks? Functional connectivity (FC) and functional connectivity dynamics (FCD) are used to represent the patterns and dynamics of functional brain networks. FC(D) is related to the synchrony of brain activity: when brain areas oscillate in a coordinated manner this yields a high correlation between their signal time series. To explain the processes underlying FC(D) we review how synchronized oscillations emerge from coupled neural populations in brain network models (BNMs). From detailed spiking networks to more abstract population models, there is strong support for the idea that the brain operates near critical instabilities that give rise to multistable or metastable dynamics that in turn lead to the intermittently synchronized slow oscillations underlying FC(D). We explore further consequences from these fundamental mechanisms and how they fit with reality. We conclude by highlighting the need for integrative brain models that connect separate mechanisms across levels of description and spatiotemporal scales and link them with cognitive function.

Sensory-motor cortices shape functional connectivity dynamics in the human brain.

Large-scale biophysical circuit models provide mechanistic insights into the micro-scale and macro-scale properties of brain organization that shape complex patterns of spontaneous brain activity. We developed a spatially heterogeneous large-scale dynamical circuit model that allowed for variation in local synaptic properties across the human cortex. Here we show that parameterizing local circuit properties with both anatomical and functional gradients generates more realistic static and dynamic resting-state functional connectivity (FC). Furthermore, empirical and simulated FC dynamics demonstrates remarkably similar sharp transitions in FC patterns, suggesting the existence of multiple attractors. Time-varying regional fMRI amplitude may track multi-stability in FC dynamics. Causal manipulation of the large-scale circuit model suggests that sensory-motor regions are a driver of FC dynamics. Finally, the spatial distribution of sensory-motor drivers matches the principal gradient of gene expression that encompasses certain interneuron classes, suggesting that heterogeneity in excitation-inhibition balance might shape multi-stability in FC dynamics.

Dynamic mode decomposition of resting-state and task fMRI.

Component analysis is a powerful tool to identify dominant patterns of interactions in multivariate datasets. In the context of fMRI data, methods such as principal component analysis or independent component analysis have been used to identify the brain networks shaping functional connectivity (FC). Importantly, these approaches are static in the sense that they ignore the temporal information contained in fMRI time series. Therefore, the corresponding components provide a static characterization of FC. Building upon recent findings suggesting that FC dynamics encode richer information about brain functional organization, we use a dynamic extension of component analysis to identify dynamic modes (DMs) of fMRI time series. We demonstrate the feasibility and relevance of this approach using resting-state and motor-task fMRI data of 730 healthy subjects of the Human Connectome Project (HCP). In resting-state, dominant DMs have strong resemblance with classical resting-state networks, with an additional temporal characterization of the networks in terms of oscillatory periods and damping times. In motor-task conditions, dominant DMs reveal interactions between several brain areas, including but not limited to the posterior parietal cortex and primary motor areas, that are not found with classical activation maps. Finally, we identify two canonical components linking the temporal properties of the resting-state DMs with 158 behavioral and demographic HCP measures. Altogether, these findings illustrate the benefits of the proposed dynamic component analysis framework, making it a promising tool to characterize the spatio-temporal organization of brain activity.

Inversion of a large-scale circuit model reveals a cortical hierarchy in the dynamic resting human brain.

We considered a large-scale dynamical circuit model of human cerebral cortex with region-specific microscale properties. The model was inverted using a stochastic optimization approach, yielding markedly better fit to new, out-of-sample resting functional magnetic resonance imaging (fMRI) data. Without assuming the existence of a hierarchy, the estimated model parameters revealed a large-scale cortical gradient. At one end, sensorimotor regions had strong recurrent connections and excitatory subcortical inputs, consistent with localized processing of external stimuli. At the opposing end, default network regions had weak recurrent connections and excitatory subcortical inputs, consistent with their role in internal thought. Furthermore, recurrent connection strength and subcortical inputs provided complementary information for differentiating the levels of the hierarchy, with only the former showing strong associations with other macroscale and microscale proxies of cortical hierarchies (meta-analysis of cognitive functions, principal resting fMRI gradient, myelin, and laminar-specific neuronal density). Overall, this study provides microscale insights into a macroscale cortical hierarchy in the dynamic resting brain.