Why Brain Criticality Is Clinically Relevant: A Scoping Review

Mapping and modeling age-related changes in intrinsic neural timescales

Intrinsic timescales of brain regions exhibit heterogeneity, escalating with hierarchical levels, and are crucial for the temporal integration of external stimuli. Aging, often associated with cognitive decline, involves progressive neuronal and synaptic loss, reshaping brain structure and dynamics. However, the impact of these structural changes on temporal coding in the aging brain remains unclear. We mapped intrinsic timescales and gray matter volume (GMV) using magnetic resonance imaging (MRI) in young and elderly adults. We found shorter intrinsic timescales across multiple large-scale functional networks in the elderly cohort, and a significant positive association between intrinsic timescales and GMV. Additionally, age-related decline in performance on visual discrimination tasks was linked to a reduction in intrinsic timescales in the cuneus. To explain these age-related shifts, we developed an age-dependent spiking neuron network model. In younger subjects, brain regions were near a critical branching regime, while regions in elderly subjects had fewer neurons and synapses, pushing the dynamics toward a subcritical regime. The model accurately reproduced the empirical results, showing longer intrinsic timescales in young adults due to critical slowing down. Our findings reveal how age-related structural brain changes may drive alterations in brain dynamics, offering testable predictions and informing possible interventions targeting cognitive decline.

Heterogeneous and higher-order cortical connectivity undergirds efficient, robust, and reliable neural codes

We hypothesized that the heterogeneous architecture of biological neural networks provides a substrate to regulate the well-known tradeoff between robustness and efficiency, thereby allowing different subpopulations of the same network to optimize for different objectives. To distinguish between subpopulations, we developed a metric based on the mathematical theory of simplicial complexes that captures the complexity of their connectivity by contrasting its higher-order structure to a random control and confirmed its relevance in several openly available connectomes. Using a biologically detailed cortical model and an electron microscopic dataset, we showed that subpopulations with low simplicial complexity exhibit efficient activity. Conversely, subpopulations of high simplicial complexity play a supporting role in boosting the reliability of the network as a whole, softening the robustness-efficiency tradeoff. Crucially, we found that both types of subpopulations can and do coexist within a single connectome in biological neural networks, due to the heterogeneity of their connectivity.

Criticality is Associated with Future Psychotherapy Response in Patients with Post-Traumatic Stress Disorder-A Pilot Study

Background

Trauma-focused psychotherapy is treatment of choice for post-traumatic stress disorder (PTSD). However, about half of patients do not respond. Recently, there is increased interest in brain criticality, which assesses the phase transition between order and disorder in brain activity. Operating close to this borderline is theorized to facilitate optimal information processing. We studied if brain criticality is related to future response to treatment, hypothesizing that treatment responders' brains function closer to criticality.

Methods

Functional magnetic resonance imaging resting-state scans were acquired from 46 male veterans with PTSD around the start of treatment. Psychotherapy consisted of trauma-focused cognitive behavioral therapy, eye movement desensitization and reprocessing, or a combination thereof. Treatment response was assessed using the Clinician-Administered PTSD Scale, and criticality was assessed using an Ising temperature approach for seven canonical brain networks (ie, the visual, somatomotor, dorsal attention, ventral attention, limbic, frontoparietal and default mode networks) to measure distance to criticality.

Results

The brains of prospective treatment responders were closer to criticality than nonresponders (P = 0.017), while no significant interaction effect between group and brain network was observed (P = 0.486). In addition, average criticality across networks correlated with future treatment response (P = 0.028).

Conclusion

These results show that the brains of prospective PTSD psychotherapy treatment responders operate closer to criticality than nonresponders, and this occurs across the entire brain instead of in separate canonical brain networks. These results suggest that effective psychotherapy is mediated by brains operating closer to criticality.

Role of homeostatic plasticity in critical brain dynamics following focal stroke lesions

Can focal brain lesions, such as those caused by stroke, disrupt critical brain dynamics? What biological mechanisms drive its recovery? In a recent study, we showed that focal lesions generate a sub-critical state that recovers over time in parallel with behavior (Rocha et al., Nat. Commun. 13, 2022). The loss of criticality in a cohort of stroke patients was associated with structural brain disconnections, while its recovery was accompanied by the re-modeling of specific white-matter tracts. These results were challenged by Janarek et al. (Sci. Rep. 13, 2023), who proposed an alternative interpretation for the anomalous monotonic decaying of the second cluster size, which is the neural signature originally used to infer loss of criticality. The present study tackles this controversy and provides evidence that the theoretical framework proposed by Janarek et al. cannot explain the anomalous cluster dynamics observed in our patients. Notably, this invalidates the claim that the brain maintains its critical dynamics regardless of the lesion severity. In addition, we explore biological mechanisms beyond white-matter remodeling that may facilitate the recovery of criticality over time. We considered two distinct scenarios: one where we suppress homeostatic plasticity, and another where we increase the excitability of brain regions. We find that suppressing homeostatic plasticity - specifically, the inhibition-excitation balance - disfavors the emergence of criticality. Conversely, increasing brain excitability can help to restore criticality when the latter is disrupted. Our results suggest that normalizing the excitation-inhibition balance is crucial for supporting recovery of critical brain dynamics.

Homeodynamic feedback inhibition control in whole-brain simulations

Simulations of large-scale brain dynamics are often impacted by overexcitation resulting from heavy-tailed structural network distributions, leading to biologically implausible simulation results. We implement a homeodynamic plasticity mechanism, known from other modeling work, in the widely used Jansen-Rit neural mass model for The Virtual Brain (TVB) simulation framework. We aim at heterogeneously adjusting the inhibitory coupling weights to reach desired dynamic regimes in each brain region. We show that, by using this dynamic approach, we can control the target activity level to obtain biologically plausible brain simulations, including post-synaptic potentials and blood-oxygen-level-dependent functional magnetic resonance imaging (fMRI) activity. We demonstrate that the derived dynamic Feedback Inhibitory Control (dFIC) can be used to enable increased variability of model dynamics. We derive the conditions under which the simulated brain activity converges to a predefined target level analytically and via simulations. We highlight the benefits of dFIC in the context of fitting the TVB model to static and dynamic measures of fMRI empirical data, accounting for global synchronization across the whole brain. The proposed novel method helps computational neuroscientists, especially TVB users, to easily "tune" brain models to desired dynamical regimes depending on the specific requirements of each study. The presented method is a steppingstone towards increased biological realism in brain network models and a valuable tool to better understand their underlying behavior.

Failure in a population: Tauopathy disrupts homeostatic set-points in emergent dynamics despite stability in the constituent neurons

Homeostatic regulation of neuronal activity is essential for robust computation; set-points, such as firing rate, are actively stabilized to compensate for perturbations. The disruption of brain function central to neurodegenerative disease likely arises from impairments of computationally essential set-points. Here, we systematically investigated the effects of tau-mediated neurodegeneration on all known set-points in neuronal activity. We continuously tracked hippocampal neuronal activity across the lifetime of a mouse model of tauopathy. We were unable to detect effects of disease in measures of single-neuron firing activity. By contrast, as tauopathy progressed, there was disruption of network-level neuronal activity, quantified by measuring neuronal pairwise interactions and criticality, a homeostatically controlled, ideal computational regime. Deviations in criticality correlated with symptoms, predicted underlying anatomical pathology, occurred in a sleep-wake-dependent manner, and could be used to reliably classify an animal's genotype. This work illustrates how neurodegeneration may disrupt the computational capacity of neurobiological systems.

Temporal complexity of the BOLD-signal in preterm versus term infants

Preterm birth causes alterations in structural and functional cerebral development that are not fully understood. Here, we investigate whether basic characteristics of BOLD signal itself might differ across preterm, term equivalent, and term infants. Anatomical, fMRI, and diffusion weighted imaging data from 716 neonates born at 23-43 weeks gestational age were obtained from the Developing Human Connectome Project. Hurst exponent (H; a measure of temporal complexity of a time-series) was computed from the power spectral density of the BOLD signal within 13 resting state networks. Using linear mixed effects models to account for scan age and birth age, we found that H increased with age, that earlier birth age contributed to lower H values, and that H increased most in motor and sensory networks. We then tested for a relationship between temporal complexity and structural development using H and DTI-based estimates of myelination and found moderate but significant correlations. These findings suggest that the temporal complexity of BOLD signal in neonates relates to age and tracks with known developmental trajectories in the brain. Elucidating how these signal-based differences might relate to maturing hemodynamics in the preterm brain could yield new information about neurophysiological vulnerabilities during this crucial developmental period.

Oscillatory and nonoscillatory sleep electroencephalographic biomarkers of the epileptic network

OBJECTIVE

In addition to the oscillatory brain activity, the nonoscillatory (scale-free) components of the background electroencephalogram (EEG) may provide further information about the complexity of the underlying neuronal network. As epilepsy is considered a network disease, such scale-free metrics might help to delineate the epileptic network. Here, we performed an analysis of the sleep oscillatory (spindle, slow wave, and rhythmic spectral power) and nonoscillatory (H exponent) intracranial EEG using multiple interictal features to estimate whether and how they deviate from normalcy in 38 adults with drug-resistant epilepsy.

METHODS

To quantify intracranial EEG abnormalities within and outside the seizure onset areas, patients' values were adjusted based on normative maps derived from the open-access Montreal Neurological Institute open iEEG Atlas. In a subset of 29 patients who underwent resective surgery, we estimated the predictive value of these features to identify the epileptogenic zone in those with a good postsurgical outcome.

RESULTS

We found that distinct sleep oscillatory and nonoscillatory metrics behave differently across the epileptic network, with the strongest differences observed for (1) a reduction in spindle activity (spindle rates and rhythmic sigma power in the 10-16 Hz band), (2) a higher rhythmic gamma power (30-80 Hz), and (3) a higher H exponent (steeper 1/f slope). As expected, epileptic spikes were also highest in the seizure onset areas. Furthermore, in surgical patients, the H exponent achieved the highest performance (balanced accuracy of .76) for classifying resected versus nonresected channels in good outcome patients.

SIGNIFICANCE

This work suggests that nonoscillatory components of the intracranial EEG signal could serve as promising interictal sleep candidates of epileptogenicity in patients with drug-resistant epilepsy. Our findings further advance the understanding of epilepsy as a disease, whereby absence or loss of sleep physiology may provide information complementary to pathological epileptic processes.

Role of the Locus Coeruleus Arousal Promoting Neurons in Maintaining Brain Criticality across the Sleep-Wake Cycle

Sleep control depends on a delicate interplay among brain regions. This generates a complex temporal architecture with numerous sleep-stage transitions and intermittent fluctuations to micro-states and brief arousals. These temporal dynamics exhibit hallmarks of criticality, suggesting that tuning to criticality is essential for spontaneous sleep-stage and arousal transitions. However, how the brain maintains criticality remains not understood. Here, we investigate θ- and δ-burst dynamics during the sleep-wake cycle of rats (Sprague-Dawley, adult male) with lesion in the wake-promoting locus coeruleus (LC). We show that, in control rats, θ- and δ-bursts exhibit power-law (θ-bursts, active phase) and exponential-like (δ-bursts, quiescent phase) duration distributions, as well as power-law long-range temporal correlations (LRTCs)-typical of non-equilibrium systems self-organizing at criticality. Furthermore, consecutive θ- and δ-bursts durations are characterized by anti-correlated coupling, indicating a new class of self-organized criticality that emerges from underlying feedback between neuronal populations and brain areas involved in generating arousals and sleep states. In contrast, we uncover that LC lesion leads to alteration of θ- and δ-burst critical features, with change in duration distributions and correlation properties, and increase in θ-δ coupling. Notably, these LC-lesion effects are opposite to those observed for lesions in the sleep-promoting ventrolateral preoptic (VLPO) nucleus. Our findings indicate that critical dynamics of θ- and δ-bursts arise from a balanced interplay of LC and VLPO, which maintains brain tuning to criticality across the sleep-wake cycle-a non-equilibrium behavior in sleep micro-architecture at short timescales that coexists with large-scale sleep-wake homeostasis.

Critical dynamics in spontaneous EEG predict anesthetic-induced loss of consciousness and perturbational complexity

Consciousness has been proposed to be supported by electrophysiological patterns poised at criticality, a dynamical regime which exhibits adaptive computational properties, maximally complex patterns and divergent sensitivity to perturbation. Here, we investigate dynamical properties of the resting-state electroencephalogram (EEG) of healthy subjects undergoing general anesthesia with propofol, xenon or ketamine. Importantly, all participants were unresponsive under anesthesia, while consciousness was retained only during ketamine anesthesia (in the form of vivid dreams), enabling an experimental dissociation between unresponsiveness and unconsciousness. For each condition, we measure (i) avalanche criticality, (ii) chaoticity, and (iii) criticality-related metrics, revealing that states of unconsciousness are characterized by a distancing from both avalanche criticality and the edge of chaos. We then ask whether these same dynamical properties are predictive of the perturbational complexity index (PCI), a TMS-based measure that has shown remarkably high sensitivity in detecting consciousness independently of behavior. We successfully predict individual subjects' PCI values with considerably high accuracy from resting-state EEG dynamical properties alone. Our results establish a firm link between perturbational complexity and criticality, and provide further evidence that criticality is a necessary condition for the emergence of consciousness.

Adolescent-to-adult gains in cognitive flexibility are adaptively supported by reward sensitivity, exploration, and neural variability

Cognitive flexibility exhibits dynamic changes throughout development, with different forms of flexibility showing dissociable developmental trajectories. In this review, we propose that an adolescent-specific mode of flexibility in the face of changing environmental contingencies supports the emergence of adolescent-to-adult gains in cognitive shifting efficiency. We first describe how cognitive shifting abilities monotonically improve from childhood to adulthood, accompanied by increases in brain state flexibility, neural variability, and excitatory/inhibitory balance. We next summarize evidence supporting the existence of a dopamine-driven, adolescent peak in flexible behavior that results in reward seeking, undirected exploration, and environmental sampling. We propose a neurodevelopmental framework that relates these adolescent behaviors to the refinement of neural phenotypes relevant to mature cognitive flexibility, and thus highlight the importance of the adolescent period in fostering healthy neurocognitive trajectories.

Criticality in Alzheimer's and healthy brains: insights from phase-ordering

Criticality, observed during second-order phase transitions, is an emergent phenomenon. The brain operates near criticality where complex systems exhibit high correlations. As a system approaches criticality, it develops "domain"-like regions with competing phases and increased spatio-temporal correlations that diverge. The dynamics of these domains depend on the system's proximity to criticality. This study explores the differences in the proximity to criticality of Alzheimer's-afflicted and cognitively normal brains through the use of a spin-lattice model derived from resting-state fMRI data and investigates the type of criticality found in the human brain - whether it is of the Ising class or something more complex. The temporal correlations in both groups display a stretched exponential nature, indicating closer alignment with the criticality of the spin-glass class rather than the Ising class. Longer relaxation times observed in cognitively normal subjects suggest increased proximity to the phase boundary. The weak distinction observed in the spatial characteristics related to proximity to criticality might once more point to a spin-glass scenario, necessitating nuanced order parameters to distinguish between phase-ordering in Alzheimer's and cognitively normal brains.

Simple model for the prediction of seizure durations

A simple model is used to simulate seizures in a population of spiking excitatory neurons experiencing a uniform effect from inhibitory neurons. A key feature is introduced into the model, i.e., a mechanism that weakens the firing thresholds. This weakening mechanism adds memory to the dynamics. We find a seizure-prone state in a "mode-switching" phase. In this phase, the system can suddenly switch from a "healthy" state with small scale-free avalanches to a "seizure" state with almost periodic large avalanches ("seizures"). Simulations of the model predict statistics for the average time spent in the seizure state (the seizure "duration") that agree with experiments and theoretical examples of similar behavior in neuronal systems. Our study points to. different connections between seizures and fracture and also offers an alternative view on the type of critical point controlling neuronal avalanches.

Neuronal synchrony and critical bistability: Mechanistic biomarkers for localizing the epileptogenic network

OBJECTIVE

Postsurgical seizure freedom in drug-resistant epilepsy (DRE) patients varies from 30% to 80%, implying that in many cases the current approaches fail to fully map the epileptogenic zone (EZ). We aimed to advance a novel approach to better characterize epileptogenicity and investigate whether the EZ encompasses a broader epileptogenic network (EpiNet) beyond the seizure zone (SZ) that exhibits seizure activity.

METHODS

We first used computational modeling to test putative complex systems-driven and systems neuroscience-driven mechanistic biomarkers for epileptogenicity. We then used these biomarkers to extract features from resting-state stereoelectroencephalograms recorded from DRE patients and trained supervised classifiers to localize the SZ against gold standard clinical localization. To further explore the prevalence of pathological features in an extended brain network outside of the clinically identified SZ, we also used unsupervised classification.

RESULTS

Supervised SZ classification trained on individual features achieved accuracies of .6-.7 area under the receiver operating characteristic curve (AUC). Combining all criticality and synchrony features further improved the AUC to .85. Unsupervised classification discovered an EpiNet-like cluster of brain regions, in which 51% of brain regions were outside of the SZ. Brain regions in the EpiNet-like cluster engaged in interareal hypersynchrony and locally exhibited high-amplitude bistability and excessive inhibition, which was strikingly similar to the high seizure risk regime revealed by our computational modeling.

SIGNIFICANCE

The finding that combining biomarkers improves SZ localization accuracy indicates that the novel mechanistic biomarkers for epileptogenicity employed here yield synergistic information. On the other hand, the discovery of SZ-like brain dynamics outside of the clinically defined SZ provides empirical evidence of an extended pathophysiological EpiNet.

Behavioral state-dependent associations between EEG temporal correlations and depressive symptoms

Previous studies have shown abnormal long-range temporal correlations in neuronal oscillations among individuals with Major Depressive Disorders, occurring during both resting states and transitions between resting and task states. However, the understanding of this effect in preclinical individuals with depression remains limited. This study investigated the association between temporal correlations of neuronal oscillations and depressive symptoms during resting and task states in preclinical individuals, specifically focusing on male action video gaming experts. Detrended fluctuation analysis (DFA), Lifetimes, and Waitingtimes were employed to explore temporal correlations across long-range and short-range scales. The results indicated widespread changes from the resting state to the task state across all frequency bands and temporal scales. Rest-task DFA changes in the alpha band exhibited a negative correlation with depressive scores at most electrodes. Significant positive correlations between DFA values and depressive scores were observed in the alpha band during the resting state but not in the task state. Similar patterns of results emerged concerning maladaptive negative emotion regulation strategies. Additionally, short-range temporal correlations in the alpha band echoed the DFA results. These findings underscore the state-dependent relationships between temporal correlations of neuronal oscillations and depressive symptoms, as well as maladaptive emotion regulation strategies, in preclinical individuals.

From pre-processing to advanced dynamic modeling of pupil data

The pupil of the eye provides a rich source of information for cognitive scientists, as it can index a variety of bodily states (e.g., arousal, fatigue) and cognitive processes (e.g., attention, decision-making). As pupillometry becomes a more accessible and popular methodology, researchers have proposed a variety of techniques for analyzing pupil data. Here, we focus on time series-based, signal-to-signal approaches that enable one to relate dynamic changes in pupil size over time with dynamic changes in a stimulus time series, continuous behavioral outcome measures, or other participants' pupil traces. We first introduce pupillometry, its neural underpinnings, and the relation between pupil measurements and other oculomotor behaviors (e.g., blinks, saccades), to stress the importance of understanding what is being measured and what can be inferred from changes in pupillary activity. Next, we discuss possible pre-processing steps, and the contexts in which they may be necessary. Finally, we turn to signal-to-signal analytic techniques, including regression-based approaches, dynamic time-warping, phase clustering, detrended fluctuation analysis, and recurrence quantification analysis. Assumptions of these techniques, and examples of the scientific questions each can address, are outlined, with references to key papers and software packages. Additionally, we provide a detailed code tutorial that steps through the key examples and figures in this paper. Ultimately, we contend that the insights gained from pupillometry are constrained by the analysis techniques used, and that signal-to-signal approaches offer a means to generate novel scientific insights by taking into account understudied spectro-temporal relationships between the pupil signal and other signals of interest.

Does the brain behave like a (complex) network? I. Dynamics

Graph theory is now becoming a standard tool in system-level neuroscience. However, endowing observed brain anatomy and dynamics with a complex network structure does not entail that the brain actually works as a network. Asking whether the brain behaves as a network means asking whether network properties count. From the viewpoint of neurophysiology and, possibly, of brain physics, the most substantial issues a network structure may be instrumental in addressing relate to the influence of network properties on brain dynamics and to whether these properties ultimately explain some aspects of brain function. Here, we address the dynamical implications of complex network, examining which aspects and scales of brain activity may be understood to genuinely behave as a network. To do so, we first define the meaning of networkness, and analyse some of its implications. We then examine ways in which brain anatomy and dynamics can be endowed with a network structure and discuss possible ways in which network structure may be shown to represent a genuine organisational principle of brain activity, rather than just a convenient description of its anatomy and dynamics.

Criticality supports cross-frequency cortical-thalamic information transfer during conscious states

Consciousness is thought to be regulated by bidirectional information transfer between the cortex and thalamus, but the nature of this bidirectional communication - and its possible disruption in unconsciousness - remains poorly understood. Here, we present two main findings elucidating mechanisms of corticothalamic information transfer during conscious states. First, we identify a highly preserved spectral channel of cortical-thalamic communication that is present during conscious states, but which is diminished during the loss of consciousness and enhanced during psychedelic states. Specifically, we show that in humans, mice, and rats, information sent from either the cortex or thalamus via δ/θ/α waves (∼1-13 Hz) is consistently encoded by the other brain region by high γ waves (52-104 Hz); moreover, unconsciousness induced by propofol anesthesia or generalized spike-and-wave seizures diminishes this cross-frequency communication, whereas the psychedelic 5-methoxy-N,N-dimethyltryptamine (5-MeO-DMT) enhances this low-to-high frequency interregional communication. Second, we leverage numerical simulations and neural electrophysiology recordings from the thalamus and cortex of human patients, rats, and mice to show that these changes in cross-frequency cortical-thalamic information transfer may be mediated by excursions of low-frequency thalamocortical electrodynamics toward/away from edge-of-chaos criticality, or the phase transition from stability to chaos. Overall, our findings link thalamic-cortical communication to consciousness, and further offer a novel, mathematically well-defined framework to explain the disruption to thalamic-cortical information transfer during unconscious states.

[Connectome in stroke patients]

Stroke is the main cause of disability among neurological diseases. There are questions of the accuracy of topical diagnosis and rehabilitation prognosis in clinical practice. Answers to these questions may be given by an approach to the study of the nervous system as a dynamic network consisting of a set of brain regions with anatomical and functional connections between them. Active study of the connectome in neurological patients in recent years became possible due to the availability of noninvasive neuroimaging methods. This review covers types of connectome and most accessible methods of obtaining research data for their construction in a neurological hospital. The review also describes resting-state networks that reflect basic brain activity in the absence of tasks. Resting-state connectivity can be used for the diagnosis of patients with severe neurological deficits. Also, changes in resting-state connectivity may indicate recovery after a stroke. The connectome analysis uses graph theory, representing the nervous system as a set of nodes and connections between them, and providing a mathematical framework allowing to study it, and methods of algebraic topology that expand the possibilities of analyzing neuroimaging data beyond graph theory. Attention is paid to the concept of self-organized criticality, which describes the brain as a system located near the critical point, where the transmission of information is most optimized. Also presented are data from studies of self-organized criticality in relation to the dynamics of recovery of patients after stroke.

Fractal Neurodynamics

The neuronal ongoing electrical activity in the brain network, the neurodynamics, reflects the structure and functionality of generating neuronal pools. The activity of neurons due to their excitatory and inhibitory projections is associated with specific brain functions. Here, the purpose was to investigate if the local ongoing electrical activity exhibits its characteristic spectral and fractal features in wakefulness and sleep across and within subjects. Moreover, we aimed to show that measures typical of complex systems catch physiological features missed by linear spectral analyses. For this study, we concentrated on the evaluation of the power spectral density (PSD) and Higuchi fractal dimension (HFD) measures. Relevant clinical impact of the specific features of neurodynamics identification stands primarily in the potential of classifying cortical parcels according to their neurodynamics as well as enhancing the effectiveness of neuromodulation interventions to cure symptoms secondary to neuronal activity unbalances.

Critical-like Brain Dynamics in a Continuum from Second- to First-Order Phase Transition

The classic brain criticality hypothesis postulates that the brain benefits from operating near a continuous second-order phase transition. Slow feedback regulation of neuronal activity could, however, lead to a discontinuous first-order transition and thereby bistable activity. Observations of bistability in awake brain activity have nonetheless remained scarce and its functional significance unclear. Moreover, there is no empirical evidence to support the hypothesis that the human brain could flexibly operate near either a first- or second-order phase transition despite such a continuum being common in models. Here, using computational modeling, we found bistable synchronization dynamics to emerge through elevated positive feedback and occur exclusively in a regimen of critical-like dynamics. We then assessed bistability in vivo with resting-state MEG in healthy adults (7 females, 11 males) and stereo-electroencephalography in epilepsy patients (28 females, 36 males). This analysis revealed that a large fraction of the neocortices exhibited varying degrees of bistability in neuronal oscillations from 3 to 200 Hz. In line with our modeling results, the neuronal bistability was positively correlated with classic assessment of brain criticality across narrow-band frequencies. Excessive bistability was predictive of epileptic pathophysiology in the patients, whereas moderate bistability was positively correlated with task performance in the healthy subjects. These empirical findings thus reveal the human brain as a one-of-a-kind complex system that exhibits critical-like dynamics in a continuum between continuous and discontinuous phase transitions.SIGNIFICANCE STATEMENT In the model, while synchrony per se was controlled by connectivity, increasing positive local feedback led to gradually emerging bistable synchrony with scale-free dynamics, suggesting a continuum between second- and first-order phase transitions in synchrony dynamics inside a critical-like regimen. In resting-state MEG and SEEG, bistability of ongoing neuronal oscillations was pervasive across brain areas and frequency bands and was observed only with concurring critical-like dynamics as the modeling predicted. As evidence for functional relevance, moderate bistability was positively correlated with executive functioning in the healthy subjects, and excessive bistability was associated with epileptic pathophysiology. These findings show that critical-like neuronal dynamics in vivo involves both continuous and discontinuous phase transitions in a frequency-, neuroanatomy-, and state-dependent manner.

Criticality of resting-state EEG predicts perturbational complexity and level of consciousness during anesthesia

Consciousness has been proposed to be supported by electrophysiological patterns poised at criticality, a dynamical regime which exhibits adaptive computational properties, maximally complex patterns and divergent sensitivity to perturbation. Here, we investigated dynamical properties of the resting-state electroencephalogram of healthy subjects undergoing general anesthesia with propofol, xenon or ketamine. We then studied the relation of these dynamic properties with the perturbational complexity index (PCI), which has shown remarkably high sensitivity in detecting consciousness independent of behavior. All participants were unresponsive under anesthesia, while consciousness was retained only during ketamine anesthesia (in the form of vivid dreams)., enabling an experimental dissociation between unresponsiveness and unconsciousness. We estimated (i) avalanche criticality, (ii) chaoticity, and (iii) criticality-related measures, and found that states of unconsciousness were characterized by a distancing from both the edge of activity propagation and the edge of chaos. We were then able to predict individual subjects' PCI (i.e., PCImax) with a mean absolute error below 7%. Our results establish a firm link between the PCI and criticality and provide further evidence for the role of criticality in the emergence of consciousness.

Criticality of neuronal avalanches in human sleep and their relationship with sleep macro- and micro-architecture

Sleep plays a key role in preserving brain function, keeping brain networks in a state that ensures optimal computation. Empirical evidence indicates that this state is consistent with criticality, where scale-free neuronal avalanches emerge. However, the connection between sleep architecture and brain tuning to criticality remains poorly understood. Here, we characterize the critical behavior of avalanches and study their relationship with sleep macro- and micro-architectures, in particular, the cyclic alternating pattern (CAP). We show that avalanches exhibit robust scaling behaviors, with exponents obeying scaling relations consistent with the mean-field directed percolation universality class. We demonstrate that avalanche dynamics is modulated by the NREM-REM cycles and that, within NREM sleep, avalanche occurrence correlates with CAP activation phases-indicating a potential link between CAP and brain tuning to criticality. The results open new perspectives on the collective dynamics underlying CAP function, and on the relationship between sleep architecture, avalanches, and self-organization to criticality.

Improvements in task performance after practice are associated with scale-free dynamics of brain activity

Although practicing a task generally benefits later performance on that same task, there are individual differences in practice effects. One avenue to model such differences comes from research showing that brain networks extract functional advantages from operating in the vicinity of criticality, a state in which brain network activity is more scale-free. We hypothesized that higher scale-free signal from fMRI data, measured with the Hurst exponent (H), indicates closer proximity to critical states. We tested whether individuals with higher H during repeated task performance would show greater practice effects. In Study 1, participants performed a dual-n-back task (DNB) twice during MRI (n = 56). In Study 2, we used two runs of n-back task (NBK) data from the Human Connectome Project sample (n = 599). In Study 3, participants performed a word completion task (CAST) across six runs (n = 44). In all three studies, multivariate analysis was used to test whether higher H was related to greater practice-related performance improvement. Supporting our hypothesis, we found patterns of higher H that reliably correlated with greater performance improvement across participants in all three studies. However, the predictive brain regions were distinct, suggesting that the specific spatial H↑ patterns are not task-general.

Linear thinking does not reflect the newer 21st-century anesthesia concepts. A narrative review

The brain constitutes a good example of a chaotic, nonlinear biological system where large neuronal networks operate chaotically with random connectivity. This critical state is significantly affected by the anesthetic loss of consciousness induced by drugs whose pharmacological behavior has been classically based on linear kinetics and dynamics. Recent developments in pharmacology and brain monitoring during anesthesia suggest a different view that we tried to explore in this article. The concepts of effect-site for hypnotic drugs modeling a maximum effect, electroencephalographic dynamics during induction, maintenance, and recovery from anesthesia are discussed, integrated into this alternative view, and how it may be applied in daily clinical practice.

Tauopathy severely disrupts homeostatic set-points in emergent neural dynamics but not in the activity of individual neurons

The homeostatic regulation of neuronal activity is essential for robust computation; key set-points, such as firing rate, are actively stabilized to compensate for perturbations. From this perspective, the disruption of brain function central to neurodegenerative disease should reflect impairments of computationally essential set-points. Despite connecting neurodegeneration to functional outcomes, the impact of disease on set-points in neuronal activity is unknown. Here we present a comprehensive, theory-driven investigation of the effects of tau-mediated neurodegeneration on homeostatic set-points in neuronal activity. In a mouse model of tauopathy, we examine 27,000 hours of hippocampal recordings during free behavior throughout disease progression. Contrary to our initial hypothesis that tauopathy would impact set-points in spike rate and variance, we found that cell-level set-points are resilient to even the latest stages of disease. Instead, we find that tauopathy disrupts neuronal activity at the network-level, which we quantify using both pairwise measures of neuron interactions as well as measurement of the network's nearness to criticality, an ideal computational regime that is known to be a homeostatic set-point. We find that shifts in network criticality 1) track with symptoms, 2) predict underlying anatomical and molecular pathology, 3) occur in a sleep/wake dependent manner, and 4) can be used to reliably classify an animal's genotype. Our data suggest that the critical set-point is intact, but that homeostatic machinery is progressively incapable of stabilizing hippocampal networks, particularly during waking. This work illustrates how neurodegenerative processes can impact the computational capacity of neurobiological systems, and suggest an important connection between molecular pathology, circuit function, and animal behavior.

Thermal Management in Neuromorphic Materials, Devices, and Networks

Machine learning has experienced unprecedented growth in recent years, often referred to as an "artificial intelligence revolution." Biological systems inspire the fundamental approach for this new computing paradigm: using neural networks to classify large amounts of data into sorting categories. Current machine-learning schemes implement simulated neurons and synapses on standard computers based on a von Neumann architecture. This approach is inefficient in energy consumption, and thermal management, motivating the search for hardware-based systems that imitate the brain. Here, the present state of thermal management of neuromorphic computing technology and the challenges and opportunities of the energy-efficient implementation of neuromorphic devices are considered. The main features of brain-inspired computing and quantum materials for implementing neuromorphic devices are briefly described, the brain criticality and resistive switching-based neuromorphic devices are discussed, the energy and electrical considerations for spiking-based computation are presented, the fundamental features of the brain's thermal regulation are addressed, the physical mechanisms for thermal management and thermoelectric control of materials and neuromorphic devices are analyzed, and challenges and new avenues for implementing energy-efficient computing are described.

Brain criticality predicts individual levels of inter-areal synchronization in human electrophysiological data

Neuronal oscillations and their synchronization between brain areas are fundamental for healthy brain function. Yet, synchronization levels exhibit large inter-individual variability that is associated with behavioral variability. We test whether individual synchronization levels are predicted by individual brain states along an extended regime of critical-like dynamics - the Griffiths phase (GP). We use computational modelling to assess how synchronization is dependent on brain criticality indexed by long-range temporal correlations (LRTCs). We analyze LRTCs and synchronization of oscillations from resting-state magnetoencephalography and stereo-electroencephalography data. Synchronization and LRTCs are both positively linearly and quadratically correlated among healthy subjects, while in epileptogenic areas they are negatively linearly correlated. These results show that variability in synchronization levels is explained by the individual position along the GP with healthy brain areas operating in its subcritical and epileptogenic areas in its supercritical side. We suggest that the GP is fundamental for brain function allowing individual variability while retaining functional advantages of criticality.

Spatial organisation of the mesoscale connectome: A feature influencing synchrony and metastability of network dynamics

Significant research has investigated synchronisation in brain networks, but the bulk of this work has explored the contribution of brain networks at the macroscale. Here we explore the effects of changing network topology on functional dynamics in spatially constrained random networks representing mesoscale neocortex. We use the Kuramoto model to simulate network dynamics and explore synchronisation and critical dynamics of the system as a function of topology in randomly generated networks with a distance-related wiring probability and no preferential attachment term. We show networks which predominantly make short-distance connections smooth out the critical coupling point and show much greater metastability, resulting in a wider range of coupling strengths demonstrating critical dynamics and metastability. We show the emergence of cluster synchronisation in these geometrically-constrained networks with functional organisation occurring along structural connections that minimise the participation coefficient of the cluster. We show that these cohorts of internally synchronised nodes also behave en masse as weakly coupled nodes and show intra-cluster desynchronisation and resynchronisation events related to inter-cluster interaction. While cluster synchronisation appears crucial to healthy brain function, it may also be pathological if it leads to unbreakable local synchronisation which may happen at extreme topologies, with implications for epilepsy research, wider brain function and other domains such as social networks.

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.

The impact of antidepressants on human neurodevelopment: Brain organoids as experimental tools

The use of antidepressants during pregnancy benefits the mother's well-being, but the effects of such substances on neurodevelopment remain poorly understood. Moreover, the consequences of early exposure to antidepressants may not be immediately apparent at birth. In utero exposure to selective serotonin reuptake inhibitors (SSRIs) has been related to developmental abnormalities, including a reduced white matter volume. Several reports have observed an increased incidence of autism spectrum disorder (ASD) and attention deficit hyperactivity disorder (ADHD) after prenatal exposure to SSRIs such as sertraline, the most widely prescribed SSRI. The advent of human-induced pluripotent stem cell (hiPSC) methods and assays now offers appropriate tools to test the consequences of such compounds for neurodevelopment in vitro. In particular, hiPSCs can be used to generate cerebral organoids - self-organized structures that recapitulate the morphology and complex physiology of the developing human brain, overcoming the limitations found in 2D cell culture and experimental animal models for testing drug efficacy and side effects. For example, single-cell RNA sequencing (scRNA-seq) and electrophysiological measurements on organoids can be used to evaluate the impact of antidepressants on the transcriptome and neuronal activity signatures in developing neurons. While the analysis of large-scale transcriptomic data depends on dimensionality reduction methods, electrophysiological recordings rely on temporal data series to discriminate statistical characteristics of neuronal activity, allowing for the rigorous analysis of the effects of antidepressants and other molecules that affect the developing nervous system, especially when applied in combination with relevant human cellular models such as brain organoids.

Intermediately synchronised brain states optimise trade-off between subject specificity and predictive capacity

Functional connectivity (FC) refers to the statistical dependencies between activity of distinct brain areas. To study temporal fluctuations in FC within the duration of a functional magnetic resonance imaging (fMRI) scanning session, researchers have proposed the computation of an edge time series (ETS) and their derivatives. Evidence suggests that FC is driven by a few time points of high-amplitude co-fluctuation (HACF) in the ETS, which may also contribute disproportionately to interindividual differences. However, it remains unclear to what degree different time points actually contribute to brain-behaviour associations. Here, we systematically evaluate this question by assessing the predictive utility of FC estimates at different levels of co-fluctuation using machine learning (ML) approaches. We demonstrate that time points of lower and intermediate co-fluctuation levels provide overall highest subject specificity as well as highest predictive capacity of individual-level phenotypes.

Biases in inverse Ising estimates of near-critical behavior

Inverse Ising inference allows pairwise interactions of complex binary systems to be reconstructed from empirical correlations. Typical estimators used for this inference, such as pseudo-likelihood maximization (PLM), are biased. Using the Sherrington-Kirkpatrick model as a benchmark, we show that these biases are large in critical regimes close to phase boundaries, and they may alter the qualitative interpretation of the inferred model. In particular, we show that the small-sample bias causes models inferred through PLM to appear closer to criticality than one would expect from the data. Data-driven methods to correct this bias are explored and applied to a functional magnetic resonance imaging data set from neuroscience. Our results indicate that additional care should be taken when attributing criticality to real-world data sets.

Brain Criticality EEG analysis for tracking neurodevelopment from Childhood to Adolescence

The brain criticality hypothesis suggests that neural networks and multiple aspects of brain activity self-organize into a critical state, and criticality marks the transition between ordered and disordered states. This hypothesis is appealing from computer science perspective because neural networks at criticality exhibit optimal processing and computing properties while having implications in clinical applications to neurological disorders. In this paper, we introduced brain criticality analysis to track neurodevelopment from childhood to adolescence using the electroencephalogram (EEG) data of 662 subjects aged 5 to 16 years from the Child Mind Institute. We computed brain criticality from long-range temporal correlation (LRTC) using detrended fluctuation analysis (DFA). We also compared the brain criticality analysis with standard EEG power analysis. The results showed a statistically significant increase in brain criticality from childhood to adolescence in the alpha band. A decreasing trend was observed in theta band from EEG power analysis, but a much higher variance was observed compared to the brain criticality analysis. However, the significant results were only observed in some EEG channels, and not observed if the analysis were performed separately with eyes-open and eyes-close condition. Nonetheless, the results suggest that brain criticality may serve as a biomarker of brain development and maturation, but further research is needed to improve brain criticality algorithms and EEG analysis methods.Clinical Relevance- The brain criticality analysis may be used to characterize and predict neurodevelopment in early childhood.

Awake or Sleeping? Maybe Both… A Review of Sleep-Related Dissociative States

Recent studies have begun to understand sleep not only as a whole-brain process but also as a complex local phenomenon controlled by specific neurotransmitters that act in different neural networks, which is called "local sleep". Moreover, the basic states of human consciousness-wakefulness, sleep onset (N1), light sleep (N2), deep sleep (N3), and rapid eye movement (REM) sleep-can concurrently appear, which may result in different sleep-related dissociative states. In this article, we classify these sleep-related dissociative states into physiological, pathological, and altered states of consciousness. Physiological states are daydreaming, lucid dreaming, and false awakenings. Pathological states include sleep paralysis, sleepwalking, and REM sleep behavior disorder. Altered states are hypnosis, anesthesia, and psychedelics. We review the neurophysiology and phenomenology of these sleep-related dissociative states of consciousness and update them with recent studies. We conclude that these sleep-related dissociative states have a significant basic and clinical impact since their study contributes to the understanding of consciousness and the proper treatment of neuropsychiatric diseases.

A systematic approach to brain dynamics: cognitive evolution theory of consciousness

The brain integrates volition, cognition, and consciousness seamlessly over three hierarchical (scale-dependent) levels of neural activity for their emergence: a causal or 'hard' level, a computational (unconscious) or 'soft' level, and a phenomenal (conscious) or 'psyche' level respectively. The cognitive evolution theory (CET) is based on three general prerequisites: physicalism, dynamism, and emergentism, which entail five consequences about the nature of consciousness: discreteness, passivity, uniqueness, integrity, and graduation. CET starts from the assumption that brains should have primarily evolved as volitional subsystems of organisms, not as prediction machines. This emphasizes the dynamical nature of consciousness in terms of critical dynamics to account for metastability, avalanches, and self-organized criticality of brain processes, then coupling it with volition and cognition in a framework unified over the levels. Consciousness emerges near critical points, and unfolds as a discrete stream of momentary states, each volitionally driven from oldest subcortical arousal systems. The stream is the brain's way of making a difference via predictive (Bayesian) processing. Its objective observables could be complexity measures reflecting levels of consciousness and its dynamical coherency to reveal how much knowledge (information gain) the brain acquires over the stream. CET also proposes a quantitative classification of both disorders of consciousness and mental disorders within that unified framework.

A scoping review for building a criticality-based conceptual framework of altered states of consciousness

The healthy conscious brain is thought to operate near a critical state, reflecting optimal information processing and high susceptibility to external stimuli. Conversely, deviations from the critical state are hypothesized to give rise to altered states of consciousness (ASC). Measures of criticality could therefore be an effective way of establishing the conscious state of an individual. Furthermore, characterizing the direction of a deviation from criticality may enable the development of treatment strategies for pathological ASC. The aim of this scoping review is to assess the current evidence supporting the criticality hypothesis, and the use of criticality as a conceptual framework for ASC. Using the PRISMA guidelines, Web of Science and PubMed were searched from inception to February 7th 2022 to find articles relating to measures of criticality across ASC. N = 427 independent papers were initially found on the subject. N = 378 were excluded because they were either: not related to criticality; not related to consciousness; not presenting results from a primary study; presenting model data. N = 49 independent papers were included in the present research, separated in 7 sub-categories of ASC: disorders of consciousness (DOC) (n = 5); sleep (n = 13); anesthesia (n = 18); epilepsy (n = 12); psychedelics and shamanic state of consciousness (n = 4); delirium (n = 1); meditative state (n = 2). Each category included articles suggesting a deviation of the critical state. While most studies were only able to identify a deviation from criticality without being certain of its direction, the preliminary consensus arising from the literature is that non-rapid eye movement (NREM) sleep reflects a subcritical state, epileptic seizures reflect a supercritical state, and psychedelics are closer to the critical state than normal consciousness. This scoping review suggests that, though the literature is limited and methodologically inhomogeneous, ASC are characterized by a deviation from criticality, though its direction is not clearly reported in a majority of studies. Criticality could become, with more extensive research, an effective and objective way to characterize ASC, and help identify therapeutic avenues to improve criticality in pathological brain states. Furthermore, we suggest how anesthesia and psychedelics could potentially be used as neuromodulation techniques to restore criticality in DOC.

Aperiodic brain activity and response to anesthesia vary in disorders of consciousness

In the human electroencephalogram (EEG), oscillatory power peaks co-exist with non-oscillatory, aperiodic activity. Although EEG analysis has traditionally focused exclusively on oscillatory power, recent investigations have shown that the aperiodic EEG component can distinguish conscious wakefulness from sleep and anesthetic-induced unconsciousness. This study investigates the aperiodic EEG component of individuals in a disorder of consciousness (DOC); how it changes in response to exposure to anesthesia; and how it relates to the brain's information richness and criticality. High-density EEG was recorded from 43 individuals in a DOC, with 16 of these individuals undergoing a protocol of propofol anesthesia. The aperiodic component was defined by the spectral slope of the power spectral density. Our results demonstrate that the EEG aperiodic component is more informative about the participants' level of consciousness than the oscillatory component, especially for patients that suffered from a stroke. Importantly, the pharmacologically induced change in the spectral slope from 30-45 Hz positively correlated with individual's pre-anesthetic level of consciousness. The pharmacologically induced loss of information-richness and criticality was associated with individual's pre-anesthetic aperiodic component. During exposure to anesthesia, the aperiodic component was correlated with 3-month recovery status for individuals with DOC. The aperiodic EEG component has been historically neglected; this research highlights the necessity of considering this measure for the assessment of individuals in DOC and future research that seeks to understand the neurophysiological underpinnings of consciousness.

Time-resolved correlation of distributed brain activity tracks E-I balance and accounts for diverse scale-free phenomena

Much of systems neuroscience posits the functional importance of brain activity patterns that lack natural scales of sizes, durations, or frequencies. The field has developed prominent, and sometimes competing, explanations for the nature of this scale-free activity. Here, we reconcile these explanations across species and modalities. First, we link estimates of excitation-inhibition (E-I) balance with time-resolved correlation of distributed brain activity. Second, we develop an unbiased method for sampling time series constrained by this time-resolved correlation. Third, we use this method to show that estimates of E-I balance account for diverse scale-free phenomena without need to attribute additional function or importance to these phenomena. Collectively, our results simplify existing explanations of scale-free brain activity and provide stringent tests on future theories that seek to transcend these explanations.

Structural Modularity Tunes Mesoscale Criticality in Biological Neuronal Networks

Numerous studies suggest that biological neuronal networks self-organize toward a critical state with stable recruitment dynamics. Individual neurons would then statistically activate exactly one further neuron during activity cascades termed neuronal avalanches. Yet, it is unclear if and how this can be reconciled with the explosive recruitment dynamics within neocortical minicolumns in vivo and within neuronal clusters in vitro, which indicates that neurons form supercritical local circuits. Theoretical studies propose that modular networks with a mix of regionally subcritical and supercritical dynamics would create apparently critical dynamics, resolving this inconsistency. Here, we provide experimental support by manipulating the structural self-organization process of networks of cultured rat cortical neurons (either sex). Consistent with the prediction, we show that increasing clustering in neuronal networks developing in vitro strongly correlates with avalanche size distributions transitioning from supercritical to subcritical activity dynamics. Avalanche size distributions approximated a power law in moderately clustered networks, indicating overall critical recruitment. We propose that activity-dependent self-organization can tune inherently supercritical networks toward mesoscale criticality by creating a modular structure in neuronal networks.SIGNIFICANCE STATEMENT Critical recruitment dynamics in neuronal networks are considered optimal for information processing in the brain. However, it remains heavily debated how neuronal networks would self-organize criticality by detailed fine-tuning of connectivity, inhibition, and excitability. We provide experimental support for theoretical considerations that modularity tunes critical recruitment dynamics at the mesoscale level of interacting neuron clusters. This reconciles reports of supercritical recruitment dynamics in local neuron clusters with findings on criticality sampled at mesoscopic network scales. Intriguingly, altered mesoscale organization is a prominent aspect of various neuropathological diseases currently investigated in the framework of criticality. We therefore believe that our findings would also be of interest for clinical scientists searching to link the functional and anatomic signatures of such brain disorders.

Altered spreading of neuronal avalanches in temporal lobe epilepsy relates to cognitive performance: A resting-state hdEEG study

OBJECTIVE

Large aperiodic bursts of activations named neuronal avalanches have been used to characterize whole-brain activity, as their presence typically relates to optimal dynamics. Epilepsy is characterized by alterations in large-scale brain network dynamics. Here we exploited neuronal avalanches to characterize differences in electroencephalography (EEG) basal activity, free from seizures and/or interictal spikes, between patients with temporal lobe epilepsy (TLE) and matched controls.

METHOD

We defined neuronal avalanches as starting when the z-scored source-reconstructed EEG signals crossed a specific threshold in any region and ending when all regions returned to baseline. This technique avoids data manipulation or assumptions of signal stationarity, focusing on the aperiodic, scale-free components of the signals. We computed individual avalanche transition matrices to track the probability of avalanche spreading across any two regions, compared them between patients and controls, and related them to memory performance in patients.

RESULTS

We observed a robust topography of significant edges clustering in regions functionally and structurally relevant for the TLE, such as the entorhinal cortex, the inferior parietal and fusiform area, the inferior temporal gyrus, and the anterior cingulate cortex. We detected a significant correlation between the centrality of the entorhinal cortex in the transition matrix and the long-term memory performance (delay recall Rey-Osterrieth Complex Figure Test).

SIGNIFICANCE

Our results show that the propagation patterns of large-scale neuronal avalanches are altered in TLE during the resting state, suggesting a potential diagnostic application in epilepsy. Furthermore, the relationship between specific patterns of propagation and memory performance support the neurophysiological relevance of neuronal avalanches.

LSD-induced increase of Ising temperature and algorithmic complexity of brain dynamics

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.

Differential associations between scale-free neural dynamics and different levels of cognitive ability

As indicators of cognitive function, scale-free neural dynamics are gaining increasing attention in cognitive neuroscience. Although the functional relevance of scale-free dynamics has been extensively reported, one fundamental question about its association with cognitive ability remains unanswered: is the association universal across a wide spectrum of cognitive abilities or confined to specific domains? Based on dual-process theory, we designed two categories of tasks to analyze two types of cognitive processes-automatic and controlled-and examined their associations with scale-free neural dynamics characterized from resting-state electroencephalography (EEG) recordings obtained from a large sample of human adults (N = 102). Our results showed that resting-state scale-free neural dynamics did not predict individuals' behavioral performance in tasks that primarily engaged the automatic process but did so in tasks that primarily engaged the controlled process. In addition, by fitting the scale-free parameters separately in different frequency bands, we found that the cognitive association of scale-free dynamics was more strongly manifested in higher-band EEG spectrum. Our findings indicate that resting-state scale-free dynamics are not universal neural indicators for all cognitive abilities but are mainly associated with high-level cognition that entails controlled processes. This finding is compatible with the widely claimed role of scale-free dynamics in reflecting properties of complex dynamic systems.

Temporal Structure in Sensorimotor Variability: A Stable Trait, But What For?

Human performance shows substantial endogenous variability over time, and this variability is a robust marker of individual differences. Of growing interest to psychologists is the realisation that variability is not fully random, but often exhibits temporal dependencies. However, their measurement and interpretation come with several controversies. Furthermore, their potential benefit for studying individual differences in healthy and clinical populations remains unclear. Here, we gather new and archival datasets featuring 11 sensorimotor and cognitive tasks across 526 participants, to examine individual differences in temporal structures. We first investigate intra-individual repeatability of the most common measures of temporal structures - to test their potential for capturing stable individual differences. Secondly, we examine inter-individual differences in these measures using: (1) task performance assessed from the same data, (2) meta-cognitive ratings of on-taskness from thought probes occasionally presented throughout the task, and (3) self-assessed attention-deficit related traits. Across all datasets, autocorrelation at lag 1 and Power Spectra Density slope showed high intra-individual repeatability across sessions and correlated with task performance. The Detrended Fluctuation Analysis slope showed the same pattern, but less reliably. The long-term component (d) of the ARFIMA(1,d,1) model showed poor repeatability and no correlation to performance. Overall, these measures failed to show external validity when correlated with either mean subjective attentional state or self-assessed traits between participants. Thus, some measures of serial dependencies may be stable individual traits, but their usefulness in capturing individual differences in other constructs typically associated with variability in performance seems limited. We conclude with comprehensive recommendations for researchers.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42113-022-00162-1.

Evidence of Critical Dynamics in Movements of Bees inside a Hive

Social insects such as honey bees exhibit complex behavioral patterns, and their distributed behavioral coordination enables decision-making at the colony level. It has, therefore, been proposed that a high-level description of their collective behavior might share commonalities with the dynamics of neural processes in brains. Here, we investigated this proposal by focusing on the possibility that brains are poised at the edge of a critical phase transition and that such a state is enabling increased computational power and adaptability. We applied mathematical tools developed in computational neuroscience to a dataset of bee movement trajectories that were recorded within the hive during the course of many days. We found that certain characteristics of the activity of the bee hive system are consistent with the Ising model when it operates at a critical temperature, and that the system's behavioral dynamics share features with the human brain in the resting state.

First few seconds for flow: A comprehensive proposal of the neurobiology and neurodynamics of state onset

Flow is a cognitive state that manifests when there is complete attentional absorption while performing a task. Flow occurs when certain internal as well as external conditions are present, including intense concentration, a sense of control, feedback, and a balance between the challenge of the task and the relevant skillset. Phenomenologically, flow is accompanied by a loss of self-consciousness, seamless integration of action and awareness, and acute changes in time perception. Research has begun to uncover some of the neurophysiological correlates of flow, as well as some of the state's neuromodulatory processes. We comprehensively review this work and consider the neurodynamics of the onset of the state, considering large-scale brain networks, as well as dopaminergic, noradrenergic, and endocannabinoid systems. To accomplish this, we outline an evidence-based hypothetical situation, and consider the flow state in a broader context including other profound alterations in consciousness, such as the psychedelic state and the state of traumatic stress that can induce PTSD. We present a broad theoretical framework which may motivate future testable hypotheses.

Absence of dynamic neural oscillatory response to environmental conditions marks childhood attention deficit hyperactivity disorder

BACKGROUND

Prior research suggests that symptoms of attention deficit hyperactivity disorder (ADHD) and related neurodevelopmental disorders may derive from alterations in the brain's ability to flexibly tune the balance between information integration and segregation and global versus local processing. This balance allows the brain to optimally filter salient stimuli in the environment and can be measured with electroencephalography (EEG) via calculation of the aperiodic spectral slope. A steeper aperiodic slope increases the capacity of global neural networks to process low-salience stimuli, while a flatter aperiodic slope reflects an emphasis on local neural networks that respond preferentially to high-salience input. Although aperiodic slope differences have been reported in ADHD, prior studies have not accounted for differing levels of stimulus input in experimental paradigms. There is evidence to suggest that dynamic shifts in neural oscillation patterns in response to changing environmental conditions could be critical for attention regulation.

METHODS

Using high-density resting EEG, we measured aperiodic spectral slope during low contrast (lights off) and high contrast (lights on) environmental conditions in a sample of 88 7-11-year-old children diagnosed with ADHD and 29 controls (30% female).

RESULTS

While controls showed a flatter aperiodic slope during the high contrast (lights on) as compared to low contrast (lights off) environmental condition, children with ADHD did not show any change in aperiodic slope across conditions.

CONCLUSIONS

This study presents a novel etiological model of biological mechanisms associated with ADHD. Children with ADHD show suboptimal modulation of intrinsic neural activity in response to changing environmental input. The dynamic spectral slope is a promising candidate biomarker for ADHD. The possibility that dynamic spectral slope is associated with cognitive-behavioral regulation more broadly merits further investigation.

Quasicriticality explains variability of human neural dynamics across life span

Aging impacts the brain's structural and functional organization and over time leads to various disorders, such as Alzheimer's disease and cognitive impairment. The process also impacts sensory function, bringing about a general slowing in various perceptual and cognitive functions. Here, we analyze the Cambridge Centre for Ageing and Neuroscience (Cam-CAN) resting-state magnetoencephalography (MEG) dataset-the largest aging cohort available-in light of the quasicriticality framework, a novel organizing principle for brain functionality which relates information processing and scaling properties of brain activity to brain connectivity and stimulus. Examination of the data using this framework reveals interesting correlations with age and gender of test subjects. Using simulated data as verification, our results suggest a link between changes to brain connectivity due to aging and increased dynamical fluctuations of neuronal firing rates. Our findings suggest a platform to develop biomarkers of neurological health.

From the origins to the stream of consciousness and its neural correlates

There are now dozens of very different theories of consciousness, each somehow contributing to our understanding of its nature. The science of consciousness needs therefore not new theories but a general framework integrating insights from those, yet not making it a still-born "Frankenstein" theory. First, the framework must operate explicitly on the stream of consciousness, not on its static description. Second, this dynamical account must also be put on the evolutionary timeline to explain the origins of consciousness. The Cognitive Evolution Theory (CET), outlined here, proposes such a framework. This starts with the assumption that brains have primarily evolved as volitional subsystems of organisms, inherited from primitive (fast and random) reflexes of simplest neural networks, only then resembling error-minimizing prediction machines. CET adopts the tools of critical dynamics to account for metastability, scale-free avalanches, and self-organization which are all intrinsic to brain dynamics. This formalizes the stream of consciousness as a discrete (transitive, irreflexive) chain of momentary states derived from critical brain dynamics at points of phase transitions and mapped then onto a state space as neural correlates of a particular conscious state. The continuous/discrete dichotomy appears naturally between the brain dynamics at the causal level and conscious states at the phenomenal level, each volitionally triggered from arousal centers of the brainstem and cognitively modulated by thalamocortical systems. Their objective observables can be entropy-based complexity measures, reflecting the transient level or quantity of consciousness at that moment.

Structure-function clustering in weighted brain networks

Functional networks, which typically describe patterns of activity taking place across the cerebral cortex, are widely studied in neuroscience. The dynamical features of these networks, and in particular their deviation from the relatively static structural network, are thought to be key to higher brain function. The interactions between such structural networks and emergent function, and the multimodal neuroimaging approaches and common analysis according to frequency band motivate a multilayer network approach. However, many such investigations rely on arbitrary threshold choices that convert dense, weighted networks to sparse, binary structures. Here, we generalise a measure of multiplex clustering to describe weighted multiplexes with arbitrarily-many layers. Moreover, we extend a recently-developed measure of structure-function clustering (that describes the disparity between anatomical connectivity and functional networks) to the weighted case. To demonstrate its utility we combine human connectome data with simulated neural activity and bifurcation analysis. Our results indicate that this new measure can extract neurologically relevant features not readily apparent in analogous single-layer analyses. In particular, we are able to deduce dynamical regimes under which multistable patterns of neural activity emerge. Importantly, these findings suggest a role for brain operation just beyond criticality to promote cognitive flexibility.