Homeostatic plasticity and emergence of functional networks in a whole-brain model at criticality

Aperiodic and Hurst EEG exponents across early human brain development: A systematic review

In electroencephalographic (EEG) data, power-frequency slope exponents (1/f_β) can provide non-invasive markers of in vivo neural activity excitation-inhibition (E:I) balance. E:I balance may be altered in neurodevelopmental conditions; hence, understanding how 1/fβ evolves across infancy/childhood has implications for developing early assessments/interventions. This systematic review (PROSPERO-ID: CRD42023363294) explored the early maturation (0-26 yrs) of resting-state EEG 1/f measures (aperiodic [AE], power law [PLE] and Hurst [HE] exponents), including studies containing ≥1 1/f measures and ≥10 typically developing participants. Five databases (including Embase and Scopus) were searched during March 2023. Forty-two studies were identified (Nparticipants=3478). Risk of bias was assessed using the Quality Assessment with Diverse Studies tool. Narrative synthesis of HE data suggests non-stationary EEG activity occurs throughout development. Age-related trends were complex, with rapid decreases in AEs during infancy and heterogenous changes thereafter. Regionally, AE maxima shifted developmentally, potentially reflecting spatial trends in maturing brain connectivity. This work highlights the importance of further characterising the development of 1/f measures to better understand how E:I balance shapes brain and cognitive development.

Mean-field solution of the neural dynamics in a Greenberg-Hastings model with excitatory and inhibitory units

We present a mean-field solution of the dynamics of a Greenberg-Hastings neural network with both excitatory and inhibitory units. We analyze the dynamical phase transitions that appear in the stationary state as the model parameters are varied. Analytical solutions are compared with numerical simulations of the microscopic model defined on a fully connected network. We found that the stationary state of this system exhibits a first-order dynamical phase transition (with the associated hysteresis) when the fraction of inhibitory units f is smaller than some critical value f_{t}≲1/2, even for a finite system. Moreover, any solution for f<1/2 can be mapped to a solution for purely excitatory systems (f=0). In finite systems, when the system is dominated by inhibition (f>f_{t}), the first-order transition is replaced by a pseudocritical one, namely a continuous crossover between regions of low and high activity that resembles the finite size behavior of a continuous phase transition order parameter. However, in the thermodynamic limit (i.e., infinite-system-size limit), we found that f_{t}→1/2 and the activity for the inhibition dominated case (f≥f_{t}) becomes negligible for any value of the parameters, while the first-order transition between low- and high-activity phases for f Collapse

Key Words Collapse

MESH Headings Collapse

Grants Collapse

Affiliation(s)

Joaquin Almeira Instituto de Física Enrique Gaviola (IFEG-CONICET), Facultad de Matemática Astronomía Física y Computación, Universidad Nacional de Córdoba, 5000 Córdoba, Argentina Consejo Nacional de Investigaciones Científcas y Tecnológicas (CONICET), 1425 Buenos Aires, Argentina Facultad de Matemática Astronomía Física y Computación, Universidad Nacional de Córdoba, 5000 Córdoba, Argentina
Tomas S Grigera Consejo Nacional de Investigaciones Científcas y Tecnológicas (CONICET), 1425 Buenos Aires, Argentina Instituto de Física de Líquidos y Sistemas Biológicos (IFLYSIB), CONICET and Universidad Nacional de La Plata, 1900 La Plata, Argentina Departamento de Física, Facultad de Ciencias Exactas, Universidad Nacional de La Plata, 1900 La Plata, Argentina Istituto dei Sistemi Complessi, Consiglio Nazionale delle Ricerche, 00185 Rome, Italy
Daniel A Martin Consejo Nacional de Investigaciones Científcas y Tecnológicas (CONICET), 1425 Buenos Aires, Argentina Instituto de Ciencias Físicas (ICIFI-CONICET), Center for Complex Systems and Brain Sciences (CEMSC3), Escuela de Ciencia y Tecnología, Universidad Nacional de Gral. San Martín, Campus Miguelete, San Martín, 1650 Buenos Aires, Argentina
Dante R Chialvo Consejo Nacional de Investigaciones Científcas y Tecnológicas (CONICET), 1425 Buenos Aires, Argentina Instituto de Ciencias Físicas (ICIFI-CONICET), Center for Complex Systems and Brain Sciences (CEMSC3), Escuela de Ciencia y Tecnología, Universidad Nacional de Gral. San Martín, Campus Miguelete, San Martín, 1650 Buenos Aires, Argentina
Sergio A Cannas Instituto de Física Enrique Gaviola (IFEG-CONICET), Facultad de Matemática Astronomía Física y Computación, Universidad Nacional de Córdoba, 5000 Córdoba, Argentina Consejo Nacional de Investigaciones Científcas y Tecnológicas (CONICET), 1425 Buenos Aires, Argentina Facultad de Matemática Astronomía Física y Computación, Universidad Nacional de Córdoba, 5000 Córdoba, Argentina

Collapse

Multi-modal and multi-model interrogation of large-scale functional brain networks

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.

Extensive topographic remapping and functional sharpening in the adult rat visual pathway upon first visual experience

Understanding the dynamics of stability/plasticity balances during adulthood is pivotal for learning, disease, and recovery from injury. However, the brain-wide topography of sensory remapping remains unknown. Here, using a first-of-its-kind setup for delivering patterned visual stimuli in a rodent magnetic resonance imaging (MRI) scanner, coupled with biologically inspired computational models, we noninvasively mapped brain-wide properties-receptive fields (RFs) and spatial frequency (SF) tuning curves-that were insofar only available from invasive electrophysiology or optical imaging. We then tracked the RF dynamics in the chronic visual deprivation model (VDM) of plasticity and found that light exposure progressively promoted a large-scale topographic remapping in adult rats. Upon light exposure, the initially unspecialized visual pathway progressively evidenced sharpened RFs (smaller and more spatially selective) and enhanced SF tuning curves. Our findings reveal that visual experience following VDM reshapes both structure and function of the visual system and shifts the stability/plasticity balance in adults.

Investigating structural and functional aspects of the brain's criticality in stroke

This paper addresses the question of the brain's critical dynamics after an injury such as a stroke. It is hypothesized that the healthy brain operates near a phase transition (critical point), which provides optimal conditions for information transmission and responses to inputs. If structural damage could cause the critical point to disappear and thus make self-organized criticality unachievable, it would offer the theoretical explanation for the post-stroke impairment of brain function. In our contribution, however, we demonstrate using network models of the brain, that the dynamics remain critical even after a stroke. In cases where the average size of the second-largest cluster of active nodes, which is one of the commonly used indicators of criticality, shows an anomalous behavior, it results from the loss of integrity of the network, quantifiable within graph theory, and not from genuine non-critical dynamics. We propose a new simple model of an artificial stroke that explains this anomaly. The proposed interpretation of the results is confirmed by an analysis of real connectomes acquired from post-stroke patients and a control group. The results presented refer to neurobiological data; however, the conclusions reached apply to a broad class of complex systems that admit a critical state.

Anomalous Self-Organization in Active Piles

Inspired by recent observations on active self-organized critical (SOC) systems, we designed an active pile (or ant pile) model with two ingredients: beyond-threshold toppling and under-threshold active motions. By including the latter component, we were able to replace the typical power-law distribution for geometric observables with a stretched exponential fat-tailed distribution, where the exponent and decay rate are dependent on the activity's strength (ζ). This observation helped us to uncover a hidden connection between active SOC systems and α-stable Levy systems. We demonstrate that one can partially sweep α-stable Levy distributions by changing ζ. The system undergoes a crossover towards Bak-Tang-Weisenfeld (BTW) sandpiles with a power-law behavior (SOC fixed point) below a crossover point ζ<ζ*≈0.1.

Recovery of neural dynamics criticality in personalized whole-brain models of stroke

The critical brain hypothesis states that biological neuronal networks, because of their structural and functional architecture, work near phase transitions for optimal response to internal and external inputs. Criticality thus provides optimal function and behavioral capabilities. We test this hypothesis by examining the influence of brain injury (strokes) on the criticality of neural dynamics estimated at the level of single participants using directly measured individual structural connectomes and whole-brain models. Lesions engender a sub-critical state that recovers over time in parallel with behavior. The improvement of criticality is associated with the re-modeling of specific white-matter connections. We show that personalized whole-brain dynamical models poised at criticality track neural dynamics, alteration post-stroke, and behavior at the level of single participants.

Similar local neuronal dynamics may lead to different collective behavior

This report is concerned with the relevance of the microscopic rules that implement individual neuronal activation, in determining the collective dynamics, under variations of the network topology. To fix ideas we study the dynamics of two cellular automaton models, commonly used, rather in-distinctively, as the building blocks of large-scale neuronal networks. One model, due to Greenberg and Hastings (GH), can be described by evolution equations mimicking an integrate-and-fire process, while the other model, due to Kinouchi and Copelli (KC), represents an abstract branching process, where a single active neuron activates a given number of postsynaptic neurons according to a prescribed "activity" branching ratio. Despite the apparent similarity between the local neuronal dynamics of the two models, it is shown that they exhibit very different collective dynamics as a function of the network topology. The GH model shows qualitatively different dynamical regimes as the network topology is varied, including transients to a ground (inactive) state, continuous and discontinuous dynamical phase transitions. In contrast, the KC model only exhibits a continuous phase transition, independently of the network topology. These results highlight the importance of paying attention to the microscopic rules chosen to model the interneuronal interactions in large-scale numerical simulations, in particular when the network topology is far from a mean-field description. One such case is the extensive work being done in the context of the Human Connectome, where a wide variety of types of models are being used to understand the brain collective dynamics.

A heuristic perspective on non-variational free energy modulation at the sleep-like edge

BACKGROUND

The variational Free Energy Principle (FEP) establishes that a neural system minimizes a free energy function of their internal state through environmental sensing entailing beliefs about hidden states in their environment.

PROBLEM

Because sensations are drastically reduced during sleep, it is still unclear how a self-organizing neural network can modulate free energy during sleep transitions.

GOAL

To address this issue, we study how network's state-dependent changes in energy, entropy and free energy connect with changes at the synaptic level in the absence of sensing during a sleep-like transition.

APPROACH

We use simulations of a physically plausible, environmentally isolated neuronal network that self-organize after inducing a thalamic input to show that the reduction of non-variational free energy depends sensitively upon thalamic input at a slow, rhythmic Poisson (delta) frequency due to spike timing dependent plasticity.

METHODS

We define a non-variational free energy in terms of the relative difference between the energy and entropy of the network from the initial distribution (prior to activity dependent plasticity) to the nonequilibrium steady-state distribution (after plasticity). We repeated the analysis under different levels of thalamic drive - as defined by the number of cortical neurons in receipt of thalamic input.

RESULTS

Entraining slow activity with thalamic input induces a transition from a gamma (awake-like state) to a delta (sleep-like state) mode of activity, which can be characterized through a modulation of network's energy and entropy (non-variational free energy) of the ensuing dynamics. The self-organizing response to low and high thalamic drive also showed characteristic differences in the spectrum of frequency content due to spike timing dependent plasticity.

CONCLUSIONS

The modulation of this non-variational free energy in a network that self-organizes, seems to be an organizational network principle. This could open a window to new empirically testable hypotheses about state changes in a neural network.

Prognostic Structural Neural Markers of MRI in Response to Mechanical Thrombectomy for Basilar Artery Occlusion

Objective: Mechanical thrombectomy (MT) has been an effective first-line therapeutic strategy for ischemic stroke. With impairment characteristics separating it from anterior circulation stroke, we aimed to explore prognostic structural neural markers for basilar artery occlusion (BAO) after MT. Methods: Fifty-four BAO patients with multi-modal magnetic resonance imaging at admission from the multicenter real-world designed BASILAR research were enrolled in this study. Features including volumes for cortical structures and subcortical regions, locations and volumes of infarctions, and white matter hyperintensity (WMH) volumes were recorded from all individuals. The impact features were identified using ANCOVA and logistic analysis. Another cohort (n = 21) was further recruited to verify the prognostic roles of screened prognostic structures. Results: For the primary clinical outcome, decreased brainstem volume and total infarction volumes from mesencephalon and midbrain were significantly related to reduced 90-day modified Rankin score (mRS) after MT treatment. WMH volume, WMH grade, average cortex thickness, white matter volume, and gray matter volume did not exhibit a remarkable relationship with the prognosis of BAO. The increased left caudate volume was obviously associated with early symptomatic recovery after MT. The prognostic role of the ratio of pons and midbrain infarct volume in brainstem was further confirmed in another cohort with area under the curve (AUC) = 0.77. Conclusions: This study was the first to provide comprehensive structural markers for the prognostic evaluation of BAO. The fully automatic and semiautomatic segmentation approaches in our study supported that the proportion of mesencephalon and midbrain infarct volume in brainstem was a crucial prognostic structural neural marker for BAO.

Default mode network activity in depression subtypes

Depression continues to carry a major disease burden worldwide, with limitations on the success of traditional pharmacological or psychological treatments. Recent approaches have therefore focused upon the neurobiological underpinnings of depression, and on the "individualization" of depression symptom profiles. One such model of depression has divided the standard diagnostic criteria into four "depression subtypes", with neurological and behavioral pathways. At the same time, attention has been focused upon the region of the brain known as the "default mode network" (DMN) and its role in attention and problem-solving. However, to date, no review has been published of the links between the DMN and the four subtypes of depression. By searching the literature studies from the last 20 years, 62 relevant papers were identified, and their findings are described for the association they demonstrate between aspects of the DMN and the four depression subtypes. It is apparent from this review that there are potential positive clinical and therapeutic outcomes from focusing upon DMN activation and connectivity, via psychological therapies, transcranial magnetic stimulation, and some emerging pharmacological models.

Cholinergic neuromodulation of inhibitory interneurons facilitates functional integration in whole-brain models

Segregation and integration are two fundamental principles of brain structural and functional organization. Neuroimaging studies have shown that the brain transits between different functionally segregated and integrated states, and neuromodulatory systems have been proposed as key to facilitate these transitions. Although whole-brain computational models have reproduced this neuromodulatory effect, the role of local inhibitory circuits and their cholinergic modulation has not been studied. In this article, we consider a Jansen & Rit whole-brain model in a network interconnected using a human connectome, and study the influence of the cholinergic and noradrenergic neuromodulatory systems on the segregation/integration balance. In our model, we introduce a local inhibitory feedback as a plausible biophysical mechanism that enables the integration of whole-brain activity, and that interacts with the other neuromodulatory influences to facilitate the transition between different functional segregation/integration regimes in the brain.

The thermodynamic brain and the evolution of intellect: the role of mental energy

The living state is low entropy, highly complex organization, yet it is part of the energy cycle of the environment. Due to the recurring presence of the resting state, stimulus and its response form a thermodynamic cycle of perception that can be modeled by the Carnot engine. The endothermic reversed Carnot engine relies on energy from the environment to increase entropy (i.e., the synaptic complexity of the resting state). High entropy relies on mental energy, which represents intrinsic motivation and focuses on the future. It increases freedom of action. The Carnot engine can model exothermic, negative emotional states, which direct the focus on the past. The organism dumps entropy and energy to its environment, in the form of aggravation, anxiety, criticism, and physical violence. The loss of mental energy curtails freedom of action, forming apathy, depression, mental diseases, and immune problems. Our improving intuition about the brain's intelligent computations will allow the development of new treatments for mental disease and novel find applications in robotics and artificial intelligence.

Critical synchronization dynamics of the Kuramoto model on connectome and small world graphs

The hypothesis, that cortical dynamics operates near criticality also suggests, that it exhibits universal critical exponents which marks the Kuramoto equation, a fundamental model for synchronization, as a prime candidate for an underlying universal model. Here, we determined the synchronization behavior of this model by solving it numerically on a large, weighted human connectome network, containing 836733 nodes, in an assumed homeostatic state. Since this graph has a topological dimension d < 4, a real synchronization phase transition is not possible in the thermodynamic limit, still we could locate a transition between partially synchronized and desynchronized states. At this crossover point we observe power-law–tailed synchronization durations, with τ_t ≃ 1.2(1), away from experimental values for the brain. For comparison, on a large two-dimensional lattice, having additional random, long-range links, we obtain a mean-field value: τ_t ≃ 1.6(1). However, below the transition of the connectome we found global coupling control-parameter dependent exponents 1 < τ_t ≤ 2, overlapping with the range of human brain experiments. We also studied the effects of random flipping of a small portion of link weights, mimicking a network with inhibitory interactions, and found similar results. The control-parameter dependent exponent suggests extended dynamical criticality below the transition point.

The road ahead in clinical network neuroscience

Clinical network neuroscience, the study of brain network topology in neurological and psychiatric diseases, has become a mainstay field within clinical neuroscience. Being a multidisciplinary group of clinical network neuroscience experts based in The Netherlands, we often discuss the current state of the art and possible avenues for future investigations. These discussions revolve around questions like "How do dynamic processes alter the underlying structural network?" and "Can we use network neuroscience for disease classification?" This opinion paper is an incomplete overview of these discussions and expands on ten questions that may potentially advance the field. By no means intended as a review of the current state of the field, it is instead meant as a conversation starter and source of inspiration to others.

Robustness of Griffiths effects in homeostatic connectome models

I provide numerical evidence for the robustness of the Griffiths phase (GP) reported previously in dynamical threshold model simulations on a large human brain network with N=836733 connected nodes. The model, with equalized network sensitivity, is extended in two ways: introduction of refractory states or by randomized time-dependent thresholds. The nonuniversal power-law dynamics in an extended control parameter region survives these modifications for a short refractory state and weak disorder. In case of temporal disorder the GP shrinks and for stronger heterogeneity disappears, leaving behind a mean-field type of critical transition. Activity avalanche size distributions below the critical point decay faster than in the original model, but the addition of inhibitory interactions sets it back to the range of experimental values.