Search

Abstract The brain is a complex non-equilibrium system capable of expressing many different dynamics as well as the transitions between them. We hypothesized that the level of non-equilibrium can serve as a signature of a given brain state, which was quantified using the arrow of time (the level of irreversibility). Using this thermodynamic framework, the irreversibility of emergent cortical activity was quantified from local field potential recordings in male Lister-hooded rats at different anesthesia levels and during the sleep-wake cycle. This measure was carried out on five distinct brain states: slow-wave sleep, awake, deep anesthesia–slow waves, light anesthesia–slow waves, and microarousals. Low levels of irreversibility were associated with synchronous activity found both in deep anesthesia and slow-wave sleep states, suggesting that slow waves were the state closest to the thermodynamic equilibrium (maximum symmetry), thus requiring minimum energy. Higher levels of irreversibility were found when brain dynamics became more asynchronous, for example, in wakefulness. These changes were also reflected in the hierarchy of cortical dynamics across different cortical areas. The neural dynamics associated with different brain states were characterized by different degrees of irreversibility and hierarchy, also acting as markers of brain state transitions. This could open new routes to monitoring, controlling, and even changing brain states in health and disease.

Abstract By connecting old and recent notions, different spatial scales, and research domains, we introduce a novel framework on the consequences of brain injury focusing on a key role of slow waves. We argue that the long-standing finding of EEG slow waves after brain injury reflects the intrusion of sleep-like cortical dynamics during wakefulness; we illustrate how these dynamics are generated and how they can lead to functional network disruption and behavioral impairment. Finally, we outline a scenario whereby post-injury slow waves can be modulated to reawaken parts of the brain that have fallen asleep to optimize rehabilitation strategies and promote recovery.

Abstract The primate including the human hippocampus implicated in episodic memory and navigation represents a spatial view, very different from the place representations in rodents. To understand this system in humans, and the computations performed, the pathway for this spatial view information to reach the hippocampus was analysed in humans. Whole-brain effective connectivity was measured with magnetoencephalography between 30 visual cortical regions and 150 other cortical regions using the HCP-MMP1 atlas in 21 participants while performing a 0-back scene memory task. In a ventromedial visual stream, V1–V4 connect to the ProStriate region where the retrosplenial scene area is located. The ProStriate region has connectivity to ventromedial visual regions VMV1–3 and VVC. These ventromedial regions connect to the medial parahippocampal region PHA1–3, which, with the VMV regions, include the parahippocampal scene area. The medial parahippocampal regions have effective connectivity to the entorhinal cortex, perirhinal cortex, and hippocampus. In contrast, when viewing faces, the effective connectivity was more through a ventrolateral visual cortical stream via the fusiform face cortex to the inferior temporal visual cortex regions TE2p and TE2a. A ventromedial visual cortical ‘Where’ stream to the hippocampus for spatial scenes was supported by diffusion topography in 171 HCP participants at 7 T.

Focal brain injuries, such as stroke, cause local structural damage as well as alteration of neuronal activity in distant brain regions. Experimental evidence suggests that one of these changes is the appearance of sleep-like slow waves in the otherwise awake individual. This pattern is prominent in areas surrounding the damaged region and can extend to connected brain regions in a way consistent with the individual's specific long-range connectivity patterns. In this paper we present a generative whole-brain model based on (f)MRI data that, in combination with the disconnection mask associated with a given patient, explains the effects of the sleep-like slow waves originated in the vicinity of the lesion area on the distant brain activity. Our model reveals new aspects of their interaction, being able to reproduce functional connectivity patterns of stroke patients and offering a detailed, causal understanding of how stroke-related effects, in particular slow waves, spread throughout the brain. The presented findings demonstrate that the model effectively captures the links between stroke occurrences, sleep-like slow waves, and their subsequent spread across the human brain.

Abstract Background The hypothesis of decreased neural inhibition in dementia has been sparsely studied in functional magnetic resonance imaging (fMRI) data across patients with different dementia subtypes, and the role of social and demographic heterogeneities on this hypothesis remains to be addressed. Methods We inferred regional inhibition by fitting a biophysical whole-brain model (dynamic mean field model with realistic inter-areal connectivity) to fMRI data from 414 participants, including patients with Alzheimer’s disease, behavioral variant frontotemporal dementia, and controls. We then investigated the effect of disease condition, and demographic and clinical variables on the local inhibitory feedback, a variable related to the maintenance of balanced neural excitation/inhibition. Results Decreased local inhibitory feedback was inferred from the biophysical modeling results in dementia patients, specific to brain areas presenting neurodegeneration. This loss of local inhibition correlated positively with years with disease, and showed differences regarding the gender and geographical origin of the patients. The model correctly reproduced known disease-related changes in functional connectivity. Conclusions Results suggest a critical link between abnormal neural and circuit-level excitability levels, the loss of grey matter observed in dementia, and the reorganization of functional connectivity, while highlighting the sensitivity of the underlying biophysical mechanism to demographic and clinical heterogeneities in the patient population.

Summary: The brain is organized hierarchically to process sensory signals. But, how do functional connections within and across areas contribute to this hierarchical order? We addressed this problem in the thalamocortical network, while monkeys detected vibrotactile stimulus. During this task, we quantified neural variability and directed functional connectivity in simultaneously recorded neurons sharing the cutaneous receptive field within and across VPL and areas 3b and 1. Before stimulus onset, VPL and area 3b exhibited similar fast dynamics while area 1 showed slower timescales. During the stimulus presence, inter-trial neural variability increased along the network VPL-3b-1 while VPL established two main feedforward pathways with areas 3b and 1 to process the stimulus. This lower variability of VPL and area 3b was found to regulate feedforward thalamocortical pathways. Instead, intra-cortical interactions were only anticipated by higher intrinsic timescales in area 1. Overall, our results provide evidence of hierarchical functional roles along the thalamocortical network.

Abstract Adolescence is a timed process with an onset, tempo, and duration. Nevertheless, the temporal dimension, especially the pace of maturation, remains an insufficiently studied aspect of developmental progression. The primary objective is to estimate the precise influence of pubertal maturational tempo on the configuration of associative brain regions. To this end, the connection between maturational stages and the level of hierarchical organization of large‐scale brain networks in 12‐13‐year‐old females is analyzed. Skeletal maturity is used as a proxy for pubertal progress. The degree of maturity is defined by the difference between bone age and chronological age. To assess the level of hierarchical organization in the brain, the temporal dynamic of closed eye resting state high‐density electroencephalography (EEG) in the alpha frequency range is analyzed. Different levels of hierarchical order are captured by the measured asymmetry in the directionality of information flow between different regions. The calculated EEG‐based entropy production of participant groups is then compared with accelerated, average, and decelerated maturity. Results indicate that an average maturational trajectory optimally aligns with cerebral hierarchical order, and both accelerated and decelerated timelines result in diminished cortical organization. This suggests that a “Goldilocks rule” of brain development is favoring a particular maturational tempo.

Abstract Brain dynamic functional connectivity characterises transient connections between brain regions. Features of brain dynamics have been linked to emotion and cognition in adult individuals, and atypical patterns have been associated with neurodevelopmental conditions such as autism. Although reliable functional brain networks have been consistently identified in neonates, little is known about the early development of dynamic functional connectivity. In this study we characterise dynamic functional connectivity with functional magnetic resonance imaging (fMRI) in the first few weeks of postnatal life in term-born (n = 324) and preterm-born (n = 66) individuals. We show that a dynamic landscape of brain connectivity is already established by the time of birth in the human brain, characterised by six transient states of neonatal functional connectivity with changing dynamics through the neonatal period. The pattern of dynamic connectivity is atypical in preterm-born infants, and associated with atypical social, sensory, and repetitive behaviours measured by the Quantitative Checklist for Autism in Toddlers (Q-CHAT) scores at 18 months of age.

Abstract Whole-brain models have proven to be useful to understand the emergence of collective activity among neural populations or brain regions. These models combine connectivity matrices, or connectomes, with local node dynamics, noise, and, eventually, transmission delays. Multiple choices for the local dynamics have been proposed. Among them, nonlinear oscillators corresponding to a supercritical Hopf bifurcation have been used to link brain connectivity and collective phase and amplitude dynamics in different brain states. Here, we studied the linear fluctuations of this model to estimate its stationary statistics, i.e., the instantaneous and lagged covariances and the power spectral densities. This linear approximation—that holds in the case of heterogeneous parameters and time-delays—allows analytical estimation of the statistics and it can be used for fast parameter explorations to study changes in brain state, changes in brain activity due to alterations in structural connectivity, and modulations of parameter due to non-equilibrium dynamics.

Abstract Background Alzheimer’s disease is a neurodegenerative condition associated with the accumulation of two misfolded proteins, amyloid-beta (A $$\beta$$ β ) and tau. We study their effect on neuronal activity, with the aim of assessing their individual and combined impact. Methods We use a whole-brain dynamic model to find the optimal parameters that best describe the effects of A $$\beta$$ β and tau on the excitation-inhibition balance of the local nodes. Results We found a clear dominance of A $$\beta$$ β over tau in the early disease stages (MCI), while tau dominates over A $$\beta$$ β in the latest stages (AD). We identify crucial roles for A $$\beta$$ β and tau in complex neuronal dynamics and demonstrate the viability of using regional distributions to define models of large-scale brain function in AD. Conclusions Our study provides further insight into the dynamics and complex interplay between these two proteins, opening the path for further investigations on biomarkers and candidate therapeutic targets in-silico.

A fundamental challenge in neuroscience is accurately defining brain states and predicting how and where to perturb the brain to force a transition. Here, we investigated resting-state fMRI data of patients suffering from disorders of consciousness (DoC) after coma (minimally conscious and unresponsive wakefulness states) and healthy controls. We applied model-free and model-based approaches to help elucidate the underlying brain mechanisms of patients with DoC. The model-free approach allowed us to characterize brain states in DoC and healthy controls as a probabilistic metastable substate (PMS) space. The PMS of each group was defined by a repertoire of unique patterns (i.e., metastable substates) with different probabilities of occurrence. In the model-based approach, we adjusted the PMS of each DoC group to a causal whole-brain model. This allowed us to explore optimal strategies for promoting transitions by applying off-line in silico probing. Furthermore, this approach enabled us to evaluate the impact of local perturbations in terms of their global effects and sensitivity to stimulation, which is a model-based biomarker providing a deeper understanding of the mechanisms underlying DoC. Our results show that transitions were obtained in a synchronous protocol, in which the somatomotor network, thalamus, precuneus and insula were the most sensitive areas to perturbation. This motivates further work to continue understanding brain function and treatments of disorders of consciousness.

Traumatic Brain Injury (TBI) is a prevalent disorder mostly characterized by persistent impairments in cognitive function that poses a substantial burden on caregivers and the healthcare system worldwide. Crucially, severity classification is primarily based on clinical evaluations, which are non-specific and poorly predictive of long-term disability. In this Mini Review, we first provide a description of our model-free and model-based approaches within the turbulent dynamics framework as well as our vision on how they can potentially contribute to provide new neuroimaging biomarkers for TBI. In addition, we report the main findings of our recent study examining longitudinal changes in moderate-severe TBI (msTBI) patients during a one year spontaneous recovery by applying the turbulent dynamics framework (model-free approach) and the Hopf whole-brain computational model (model-based approach) combined with in silico perturbations. Given the neuroinflammatory response and heightened risk for neurodegeneration after TBI, we also offer future directions to explore the association with genomic information. Moreover, we discuss how whole-brain computational modeling may advance our understanding of the impact of structural disconnection on whole-brain dynamics after msTBI in light of our recent findings. Lastly, we suggest future avenues whereby whole-brain computational modeling may assist the identification of optimal brain targets for deep brain stimulation to promote TBI recovery.

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

Abstract Scale invariance is a characteristic of neural activity. How this property emerges from neural interactions remains a fundamental question. Here, we studied the relation between scale-invariant brain dynamics and structural connectivity by analyzing human resting-state (rs-) fMRI signals, together with diffusion MRI (dMRI) connectivity and its approximation as an exponentially decaying function of the distance between brain regions. We analyzed the rs-fMRI dynamics using functional connectivity and a recently proposed phenomenological renormalization group (PRG) method that tracks the change of collective activity after successive coarse-graining at different scales. We found that brain dynamics display power-law correlations and power-law scaling as a function of PRG coarse-graining based on functional or structural connectivity. Moreover, we modeled the brain activity using a network of spins interacting through large-scale connectivity and presenting a phase transition between ordered and disordered phases. Within this simple model, we found that the observed scaling features were likely to emerge from critical dynamics and connections exponentially decaying with distance. In conclusion, our study tests the PRG method using large-scale brain activity and theoretical models and suggests that scaling of rs-fMRI activity relates to criticality.

Brain signal irreversibility has been shown to be a promising approach to study neural dynamics. Nevertheless, the relation with cortical hierarchy and the influence of different electrophysiological features is not completely understood. In this study, we recorded local field potentials (LFPs) during spontaneous behavior, including awake and sleep periods, using custom micro-electrocorticographic (μECoG) arrays implanted in ferrets. In contrast to humans, ferrets remain less time in each state across the sleep-wake cycle. We deployed a diverse set of metrics in order to measure the levels of complexity of the different behavioral states. In particular, brain irreversibility, which is a signature of non-equilibrium dynamics, captured by the arrow of time of the signal, revealed the hierarchical organization of the ferret's cortex. We found different signatures of irreversibility and functional hierarchy of large-scale dynamics in three different brain states (active awake, quiet awake, and deep sleep), showing a lower level of irreversibility in the deep sleep stage, compared to the other. Irreversibility also allowed us to disentangle the influence of different cortical areas and frequency bands in this process, showing a predominance of the parietal cortex and the theta band. Furthermore, when inspecting the embedded dynamic through a Hidden Markov Model, the deep sleep stage was revealed to have a lower switching rate and lower entropy production. These results suggest functional hierarchies in organization that can be revealed through thermodynamic features and information theory metrics.

Abstract To better understand how network structure shapes intelligent behavior, we developed a learning algorithm that we used to build personalized brain network models for 650 Human Connectome Project participants. We found that participants with higher intelligence scores took more time to solve difficult problems, and that slower solvers had higher average functional connectivity. With simulations we identified a mechanistic link between functional connectivity, intelligence, processing speed and brain synchrony for trading accuracy with speed in dependence of excitation-inhibition balance. Reduced synchrony led decision-making circuits to quickly jump to conclusions, while higher synchrony allowed for better integration of evidence and more robust working memory. Strict tests were applied to ensure reproducibility and generality of the obtained results. Here, we identify links between brain structure and function that enable to learn connectome topology from noninvasive recordings and map it to inter-individual differences in behavior, suggesting broad utility for research and clinical applications.

Abstract Brain circuits display modular architecture at different scales of organization. Such neural assemblies are typically associated to functional specialization but the mechanisms leading to their emergence and consolidation still remain elusive. In this paper we investigate the role of inhibition in structuring new neural assemblies driven by the entrainment to various inputs. In particular, we focus on the role of partially synchronized dynamics for the creation and maintenance of structural modules in neural circuits by considering a network of excitatory and inhibitory $$\theta$$ θ -neurons with plastic Hebbian synapses. The learning process consists of an entrainment to temporally alternating stimuli that are applied to separate regions of the network. This entrainment leads to the emergence of modular structures. Contrary to common practice in artificial neural networks—where the acquired weights are typically frozen after the learning session—we allow for synaptic adaptation even after the learning phase. We find that the presence of inhibitory neurons in the network is crucial for the emergence and the post-learning consolidation of the modular structures. Indeed networks made of purely excitatory neurons or of neurons not respecting Dale’s principle are unable to form or to maintain the modular architecture induced by the stimuli. We also demonstrate that the number of inhibitory neurons in the network is directly related to the maximal number of neural assemblies that can be consolidated, supporting the idea that inhibition has a direct impact on the memory capacity of the neural network.

Abstract Fast, efficient information transfer is essential for the brain to ensure survival. As recently shown in functional magnetic resonance imaging with high spatial resolution, turbulence appears to offer a fundamental way to facilitate energy and information transfer across spatiotemporal scales in brain dynamics. However, given that this imaging modality is comparably slow and not directly linked with neuronal activity, here we investigated the existence of turbulence in fast whole-brain neural dynamics measured with magnetoencephalography (MEG). The coarse spatial observations in MEG necessitated that we created and validated a empirical measure of turbulence. We found that the measure of edge-centric metastability perfectly detected turbulence in a ring of non-local coupled oscillators where the ground-truth was analytically known, even at a coarse spatial scale of observations. This allowed us to use this measure in the spatially coarse, empirical large-scale MEG data from 89 human participants. We demonstrated turbulence in fast neuronal dynamics and used this to quantify information transfer in the brain. The results demonstrate that the necessary efficiency of brain function is dependent on an underlying turbulent regime.

AbstractLarge variability exists across brain regions in health and disease, considering their cellular and molecular composition, connectivity, and function. Large-scale whole-brain models comprising coupled brain regions provide insights into the underlying dynamics that shape complex patterns of spontaneous brain activity. In particular, biophysically grounded mean-field whole-brain models in the asynchronous regime were used to demonstrate the dynamical consequences of including regional variability. Nevertheless, the role of heterogeneities when brain dynamics are supported by synchronous oscillating state, which is a ubiquitous phenomenon in brain, remains poorly understood. Here, we implemented two models capable of presenting oscillatory behavior with different levels of abstraction: a phenomenological Stuart–Landau model and an exact mean-field model. The fit of these models informed by structural- to functional-weighted MRI signal (T1w/T2w) allowed us to explore the implication of the inclusion of heterogeneities for modeling resting-state fMRI recordings from healthy participants. We found that disease-specific regional functional heterogeneity imposed dynamical consequences within the oscillatory regime in fMRI recordings from neurodegeneration with specific impacts on brain atrophy/structure (Alzheimer’s patients). Overall, we found that models with oscillations perform better when structural and functional regional heterogeneities are considered, showing that phenomenological and biophysical models behave similarly at the brink of the Hopf bifurcation.

AbstractBeyond the established effects of subthalamic nucleus deep brain stimulation (STN-DBS) in reducing motor symptoms in Parkinson’s disease, recent evidence has highlighted the effect on non-motor symptoms. However, the impact of STN-DBS on disseminated networks remains unclear. This study aimed to perform a quantitative evaluation of network-specific modulation induced by STN-DBS using Leading Eigenvector Dynamics Analysis (LEiDA). We calculated the occupancy of resting-state networks (RSNs) in functional MRI data from 10 patients with Parkinson’s disease implanted with STN-DBS and statistically compared between ON and OFF conditions. STN-DBS was found to specifically modulate the occupancy of networks overlapping with limbic RSNs. STN-DBS significantly increased the occupancy of an orbitofrontal limbic subsystem with respect to both DBS OFF (p = 0.0057) and 49 age-matched healthy controls (p = 0.0033). Occupancy of a diffuse limbic RSN was increased with STN-DBS OFF when compared with healthy controls (p = 0.021), but not when STN-DBS was ON, which indicates rebalancing of this network. These results highlight the modulatory effect of STN-DBS on components of the limbic system, particularly within the orbitofrontal cortex, a structure associated with reward processing. These results reinforce the value of quantitative biomarkers of RSN activity in evaluating the disseminated impact of brain stimulation techniques and the personalization of therapeutic strategies.

Traditionally, in neuroimaging, model-free analyses are used to find significant differences between brain states via signal detection theory. Depending on the a priori assumptions about the underlying data, different spatio-temporal features can be analysed. Alternatively, model-based techniques infer features from the data and compare significance from model parameters. However, to assess transitions from one brain state to another remains a challenge in current paradigms. Here, we introduce a “Dynamic Sensitivity Analysis” framework that quantifies transitions between brain states in terms of stimulation ability to rebalance spatio-temporal brain activity towards a target state such as healthy brain dynamics. In practice, it means building a whole-brain model fitted to the spatio-temporal description of brain dynamics, and applying systematic stimulations in-silico to assess the optimal strategy to drive brain dynamics towards a target state. Further, we show how Dynamic Sensitivity Analysis extends to various brain stimulation paradigms, ultimately contributing to improving the efficacy of personalised clinical interventions.

Mapping brain connections is critical for decoding brain functions. Here, the authors present an integrated resource of awake resting-state fMRI and neuronal tracing data of marmosets to understand structural-functional relationships of brain connections.

There are several models of how serotonergic psychedelic drugs affect brain activity. Here the authors use network control theory and functional MRI data to provide evidence that serotonin receptor agonists LSD and psilocybin flatten the brain’s dynamic landscape, allowing for facile state transitions and more temporally diverse brain activity.

Existing whole-brain models are generally tailored to the modelling of a particular data modality (e.g., fMRI or MEG/EEG). We propose that despite the differing aspects of neural activity each modality captures, they originate from shared network dynamics. Building on the universal principles of self-organising delay-coupled nonlinear systems, we aim to link distinct features of brain activity - captured across modalities - to the dynamics unfolding on a macroscopic structural connectome.To jointly predict connectivity, spatiotemporal and transient features of distinct signal modalities, we consider two large-scale models - the Stuart Landau and Wilson and Cowan models - which generate short-lived 40 Hz oscillations with varying levels of realism. To this end, we measure features of functional connectivity and metastable oscillatory modes (MOMs) in fMRI and MEG signals - and compare them against simulated data.We show that both models can represent MEG functional connectivity (FC), functional connectivity dynamics (FCD) and generate MOMs to a comparable degree. This is achieved by adjusting the global coupling and mean conduction time delay and, in the WC model, through the inclusion of balance between excitation and inhibition. For both models, the omission of delays dramatically decreased the performance. For fMRI, the SL model performed worse for FCD and MOMs, highlighting the importance of balanced dynamics for the emergence of spatiotemporal and transient patterns of ultra-slow dynamics. Notably, optimal working points varied across modalities and no model was able to achieve a correlation with empirical FC higher than 0.4 across modalities for the same set of parameters. Nonetheless, both displayed the emergence of FC patterns that extended beyond the constraints of the anatomical structure. Finally, we show that both models can generate MOMs with empirical-like properties such as size (number of brain regions engaging in a mode) and duration (continuous time interval during which a mode appears).Our results demonstrate the emergence of static and dynamic properties of neural activity at different timescales from networks of delay-coupled oscillators at 40 Hz. Given the higher dependence of simulated FC on the underlying structural connectivity, we suggest that mesoscale heterogeneities in neural circuitry may be critical for the emergence of parallel cross-modal functional networks and should be accounted for in future modelling endeavours.

Characterising brain states during tasks is common practice for many neuroscientific experiments using electrophysiological modalities such as electroencephalography (EEG) and magnetoencephalography (MEG). Brain states are often described in terms of oscillatory power and correlated brain activity, i.e. functional connectivity. It is, however, not unusual to observe weak task induced functional connectivity alterations in the presence of strong task induced power modulations using classical time-frequency representation of the data. Here, we propose that non-reversibility, or the temporal asymmetry in functional interactions, may be more sensitive to characterise task induced brain states than functional connectivity. As a second step, we explore causal mechanisms of non-reversibility in MEG data using whole brain computational models. We include working memory, motor, language tasks and resting-state data from participants of the Human Connectome Project (HCP). Non-reversibility is derived from the lagged amplitude envelope correlation (LAEC), and is based on asymmetry of the forward and reversed cross-correlations of the amplitude envelopes. Using random forests, we find that non-reversibility outperforms functional connectivity in the identification of task induced brain states. Non-reversibility shows especially better sensitivity to capture bottom-up gamma induced brain states across all tasks, but also alpha band associated brain states. Using whole brain computational models we find that asymmetry in the effective connectivity and axonal conduction delays play a major role in shaping non-reversibility across the brain. Our work paves the way for better sensitivity in characterising brain states during both bottom-up as well as top-down modulation in future neuroscientific experiments.

Disorders of consciousness are complex conditions characterised by persistent loss of responsiveness due to brain injury. They present diagnostic challenges and limited options for treatment, and highlight the urgent need for a more thorough understanding of how human consciousness arises from coordinated neural activity. The increasing availability of multimodal neuroimaging data has given rise to a wide range of clinically- and scientifically-motivated modelling efforts, seeking to improve data-driven stratification of patients, to identify causal mechanisms for patient pathophysiology and loss of consciousness more broadly, and to develop simulations as a means of testing in silico potential treatment avenues to restore consciousness. As a dedicated Working Group of clinicians and neuroscientists of the international Curing Coma Campaign, here we provide our framework and vision to understand the diverse statistical and generative computational modelling approaches that are being employed in this fast-growing field. We identify the gaps that exist between the current state-of-the-art in statistical and biophysical computational modelling in human neuroscience, and the aspirational goal of a mature field of modelling disorders of consciousness; which might drive improved treatments and outcomes in the clinic. Finally, we make several recommendations for how the field as a whole can work together to address these challenges.

Favaretto et al. show that the brain rapidly alternates between transient connectivity patterns, with cortical regions flexibly synchronizing with two groups of subcortical regions, and that this dynamic is abnormal in stroke patients.

The mechanisms underlying transient brain rhythms and weekly stable synchronization of distant brain areas and their link with neural activity is still a matter of debate. Here, the authors use a brain network model to study spatio-temporal synchronization dynamics of brain regions and find that there is an optimal regime where spatially and spectrally resolved metastable oscillatory modes, similar to human magnetoencephalography data, emerge.

A unifying turbulent dynamics framework using both model-free and modelbased measures of whole-brain information provides insights into brain states.

A thermodynamic-inspired framework enables researchers to quantify the balance between intrinsic and extrinsic dynamics in brain signals, providing further insight into how the brain behaves during sleep and when under anesthesia.

Summary: Brain states are frequently represented using a unidimensional scale measuring the richness of subjective experience (level of consciousness). This description assumes a mapping between the high-dimensional space of whole-brain configurations and the trajectories of brain states associated with changes in consciousness, yet this mapping and its properties remain unclear. We combine whole-brain modeling, data augmentation, and deep learning for dimensionality reduction to determine a mapping representing states of consciousness in a low-dimensional space, where distances parallel similarities between states. An orderly trajectory from wakefulness to patients with brain injury is revealed in a latent space whose coordinates represent metrics related to functional modularity and structure-function coupling, increasing alongside loss of consciousness. Finally, we investigate the effects of model perturbations, providing geometrical interpretation for the stability and reversibility of states. We conclude that conscious awareness depends on functional patterns encoded as a low-dimensional trajectory within the vast space of brain configurations.

Spatiotemporal oscillations underlie all cognitive brain functions. Large-scale brain models, constrained by neuroimaging data, aim to trace the principles underlying such macroscopic neural activity from the intricate and multi-scale structure of the brain. Despite substantial progress in the field, many aspects about the mechanisms behind the onset of spatiotemporal neural dynamics are still unknown. In this work we establish a simple framework for the emergence of complex brain dynamics, including high-dimensional chaos and travelling waves. The model consists of a complex network of 90 brain regions, whose structural connectivity is obtained from tractography data. The activity of each brain area is governed by a Jansen neural mass model and we normalize the total input received by each node so it amounts the same across all brain areas. This assumption allows for the existence of an homogeneous invariant manifold, i.e., a set of different stationary and oscillatory states in which all nodes behave identically. Stability analysis of these homogeneous solutions unveils a transverse instability of the synchronized state, which gives rise to different types of spatiotemporal dynamics, such as chaotic alpha activity. Additionally, we illustrate the ubiquity of this route towards complex spatiotemporal activity in a network of next generation neural mass models. Altogehter, our results unveil the bifurcation landscape that underlies the emergence of function from structure in the brain.

Cortical function emerges from the interactions of multi-scale networks that may be studied at a high level using neural mass models (NMM) that represent the mean activity of large numbers of neurons. Here, we provide first a new framework called laminar NMM, or LaNMM for short, where we combine conduction physics with NMMs to simulate electrophysiological measurements. Then, we employ this framework to infer the location of oscillatory generators from laminar-resolved data collected from the prefrontal cortex in the macaque monkey. We define a minimal model capable of generating coupled slow and fast oscillations, and we optimize LaNMM-specific parameters to fit multi-contact recordings. We rank the candidate models using an optimization function that evaluates the match between the functional connectivity (FC) of the model and data, where FC is defined by the covariance between bipolar voltage measurements at different cortical depths. The family of best solutions reproduces the FC of the observed electrophysiology by selecting locations of pyramidal cells and their synapses that result in the generation of fast activity at superficial layers and slow activity across most depths, in line with recent literature proposals. In closing, we discuss how this hybrid modeling framework can be more generally used to infer cortical circuitry.

The treatment of neurodegenerative diseases is hindered by lack of interventions capable of steering multimodal whole-brain dynamics towards patterns indicative of preserved brain health. To address this problem, we combined deep learning with a model capable of reproducing whole-brain functional connectivity in patients diagnosed with Alzheimer’s disease (AD) and behavioral variant frontotemporal dementia (bvFTD). These models included disease-specific atrophy maps as priors to modulate local parameters, revealing increased stability of hippocampal and insular dynamics as signatures of brain atrophy in AD and bvFTD, respectively. Using variational autoencoders, we visualized different pathologies and their severity as the evolution of trajectories in a low-dimensional latent space. Finally, we perturbed the model to reveal key AD- and bvFTD-specific regions to induce transitions from pathological to healthy brain states. Overall, we obtained novel insights on disease progression and control by means of external stimulation, while identifying dynamical mechanisms that underlie functional alterations in neurodegeneration.

Abstract The stock market is a bellwether of socio-economic changes that may directly affect individual well-being. Using large-scale UK-biobank data generated over 14 years, we applied specification curve analysis to rigorously identify significant associations between the local stock market index (FTSE100) and 479,791 UK residents’ mood, as well as their alcohol intake and blood pressure adjusting the results for a large number of potential confounders, including age, sex, linear and non-linear effects of time, research site, other stock market indexes. Furthermore, we found similar associations between FTSE100 and volumetric measures of affective brain regions in a subsample (n = 39,755; measurements performed over 5.5 years), which were particularly strong around phase transitions characterized by maximum volatility in the market. The main findings did not depend on applied effect-size estimation criteria (linear methods or mutual information criterion) and were replicated in two independent US-based studies (Parkinson’s Progression Markers Initiative; n = 424; performed over 2.5 years and MyConnectome; n = 1; 81 measurements over 1.5 years). Our results suggest that phase transitions in the society, indexed by stock market, exhibit close relationships with human mood, health and the affective brain from an individual to population level.

We investigate how brain activity can be supported by a turbulent regime based on the deviations of a self-similar scaling of high-order structure functions within the phenomenological Kolmogorov's theory. By analyzing a large neuroimaging data set, we establish the relationship between scaling exponents and their order, showing that brain activity has more than one invariant scale, and thus orders higher than 2 are needed to accurately describe its underlying statistical properties. Furthermore, we build whole-brain models of coupled oscillators to show that high-order information allows for a better description of the brain's empirical information transmission and reactivity.

AbstractPsychedelic drugs show promise as safe and effective treatments for neuropsychiatric disorders, yet their mechanisms of action are not fully understood. A fundamental hypothesis is that psychedelics work by dose-dependently changing the functional hierarchy of brain dynamics, but it is unclear whether different psychedelics act similarly. Here, we investigated the changes in the brain’s functional hierarchy associated with two different psychedelics (LSD and psilocybin). Using a novel turbulence framework, we were able to determine the vorticity, that is, the local level of synchronization, that allowed us to extend the standard global time-based measure of metastability to become a local-based measure of both space and time. This framework produced detailed signatures of turbulence-based hierarchical change for each psychedelic drug, revealing consistent and discriminate effects on a higher level network, that is, the default mode network. Overall, our findings directly support a prior hypothesis that psychedelics modulate (i.e., “compress”) the functional hierarchy and provide a quantification of these changes for two different psychedelics. Implications for therapeutic applications of psychedelics are discussed.

AbstractSpontaneous brain activity changes across states of consciousness. A particular consciousness-mediated configuration is the anticorrelations between the default mode network and other brain regions. What this antagonistic organization implies about consciousness to date remains inconclusive. In this Perspective Article, we propose that anticorrelations are the physiological expression of the concept of segregation, namely the brain’s capacity to show selectivity in the way areas will be functionally connected. We postulate that this effect is mediated by the process of neural inhibition, by regulating global and local inhibitory activity. While recognizing that this effect can also result from other mechanisms, neural inhibition helps the understanding of how network metastability is affected after disrupting local and global neural balance. In combination with relevant theories of consciousness, we suggest that anticorrelations are a physiological prior that can work as a marker of preserved consciousness. We predict that if the brain is not in a state to host anticorrelations, then most likely the individual does not entertain subjective experience. We believe that this link between anticorrelations and the underlying physiology will help not only to comprehend how consciousness happens, but also conceptualize effective interventions for treating consciousness disorders in which anticorrelations seem particularly affected.

A topic of growing interest in computational neuroscience is the discovery of fundamental principles underlying global dynamics and the self-organization of the brain. In particular, the notion that the brain operates near criticality has gained considerable support, and recent work has shown that the dynamics of different brain states may be modeled by pairwise maximum entropy Ising models at various distances from a phase transition, i.e., from criticality. Here we aim to characterize two brain states (psychedelics-induced and placebo) as captured by functional magnetic resonance imaging (fMRI), with features derived from the Ising spin model formalism (system temperature, critical point, susceptibility) and from algorithmic complexity. We hypothesized, along the lines of the entropic brain hypothesis, that psychedelics drive brain dynamics into a more disordered state at a higher Ising temperature and increased complexity. We analyze resting state blood-oxygen-level-dependent (BOLD) fMRI data collected in an earlier study from fifteen subjects in a control condition (placebo) and during ingestion of lysergic acid diethylamide (LSD). Working with the automated anatomical labeling (AAL) brain parcellation, we first create "archetype" Ising models representative of the entire dataset (global) and of the data in each condition. Remarkably, we find that such archetypes exhibit a strong correlation with an average structural connectome template obtained from dMRI (r = 0.6). We compare the archetypes from the two conditions and find that the Ising connectivity in the LSD condition is lower than in the placebo one, especially in homotopic links (interhemispheric connectivity), reflecting a significant decrease of homotopic functional connectivity in the LSD condition. The global archetype is then personalized for each individual and condition by adjusting the system temperature. The resulting temperatures are all near but above the critical point of the model in the paramagnetic (disordered) phase. The individualized Ising temperatures are higher in the LSD condition than in the placebo condition (p = 9 × 10-5). Next, we estimate the Lempel-Ziv-Welch (LZW) complexity of the binarized BOLD data and the synthetic data generated with the individualized model using the Metropolis algorithm for each participant and condition. The LZW complexity computed from experimental data reveals a weak statistical relationship with condition (p = 0.04 one-tailed Wilcoxon test) and none with Ising temperature (r(13) = 0.13, p = 0.65), presumably because of the limited length of the BOLD time series. Similarly, we explore complexity using the block decomposition method (BDM), a more advanced method for estimating algorithmic complexity. The BDM complexity of the experimental data displays a significant correlation with Ising temperature (r(13) = 0.56, p = 0.03) and a weak but significant correlation with condition (p = 0.04, one-tailed Wilcoxon test). This study suggests that the effects of LSD increase the complexity of brain dynamics by loosening interhemispheric connectivity-especially homotopic links. In agreement with earlier work using the Ising formalism with BOLD data, we find the brain state in the placebo condition is already above the critical point, with LSD resulting in a shift further away from criticality into a more disordered state.

Spontaneous fluctuations in brain activity exhibit complex spatiotemporal patterns across animal species. Here the authors show that sensory-motor regions and spatial heterogeneity in excitation-inhibition balance might shape multi-stability in brain dynamics.

Autism spectrum disorder (ASD) is characterised by synaptic surplus and atypical functional connectivity. Here, the authors show that synaptic pathology in Tsc2 haploinsufficient mice is associated with autism-like behavior and cortico-striatal hyperconnectivity, and that analogous functional hyperconnectivity signatures can be linked to mTOR-pathway dysfunction in subgroups of children with idiopathic ASD.

López-González et al study the fMRI brain dynamics and their underlying mechanism from patients that suffered brain injuries leading to a disorder of consciousness as well as from healthy subjects undergoing propofol-induced sedation. They show that pathological and pharmacological low-level states of consciousness display disrupted synchronization patterns, higher constraint to the anatomy and a loss of heterogeneity and stability in the structural hubs compared to conscious states.

Despite decades of research, there is still a lack of understanding of the role and generating mechanisms of the ubiquitous fluctuations and oscillations found in recordings of brain dynamics. Here, we used whole-brain computational models capable of presenting different dynamical regimes to reproduce empirical data's turbulence level. We showed that the model's fluctuations regime fitted to turbulence more faithfully reproduces the empirical functional connectivity compared to oscillatory and noise regimes. By applying global and local strength-dependent perturbations and subsequently measuring the responsiveness of the model, we revealed each regime's computational capacity demonstrating that brain dynamics is shifted towards fluctuations to provide much-needed flexibility. Importantly, fluctuation regime stimulation in a brain region within a given resting state network modulates that network, aligned with previous empirical and computational studies. Furthermore, this framework generates specific, testable empirical predictions for human stimulation studies using strength-dependent rather than constant perturbation. Overall, the whole-brain models fitted to the level of empirical turbulence together with functional connectivity unveil that the fluctuation regime best captures empirical data, and the strength-dependent perturbative framework demonstrates how this regime provides maximal flexibility to the human brain.

Managing cognitive load depends on adequate resource allocation by the human brain through the engagement of metastable substates, which are large-scale functional networks that change over time. We employed a novel analysis method, deep autoencoder dynamical analysis (DADA), with 100 healthy adults selected from the Human Connectome Project (HCP) data set in rest and six cognitive tasks. The deep autoencoder of DADA described seven recurrent stochastic metastable substates from the functional connectome of BOLD phase coherence matrices. These substates were significantly differentiated in terms of their probability of appearance, time duration, and spatial attributes. We found that during different cognitive tasks, there was a higher probability of having more connected substates dominated by a high degree of connectivity in the thalamus. In addition, compared with those during tasks, resting brain dynamics have a lower level of predictability, indicating a more uniform distribution of metastability between substates, quantified by higher entropy. These novel findings provide empirical evidence for the philosophically motivated cognitive theory, suggesting on-line and off-line as two fundamentally distinct modes of cognition. On-line cognition refers to task-dependent engagement with the sensory input, while off-line cognition is a slower, environmentally detached mode engaged with decision and planning. Overall, the DADA framework provides a bridge between neuroscience and cognitive theory that can be further explored in the future.