The criticality hypothesis: how local cortical networks might optimize information processing

The dynamics of the stress neuromatrix

Stressful stimuli in healthy subjects trigger activation of a consistent and reproducible set of brain regions; yet, the notion that there is a single and constant stress neuromatrix is not sustainable. Indeed, after chronic stress exposure there is activation of many brain regions outside that network. This suggests that there is a distinction between the acute and the chronic stress neuromatrix. Herein, a new working model is proposed to understand the shift between these networks. The understanding of the factors that modulate these networks and their interplay will allow for a more comprehensive and holistic perspective of how the brain shifts 'back and forth' from a healthy to a stressed pattern and, ultimately, how the latter can be a trigger for several neurological and psychiatric conditions.

Optimal Information Representation and Criticality in an Adaptive Sensory Recurrent Neuronal Network

Recurrent connections play an important role in cortical function, yet their exact contribution to the network computation remains unknown. The principles guiding the long-term evolution of these connections are poorly understood as well. Therefore, gaining insight into their computational role and into the mechanism shaping their pattern would be of great importance. To that end, we studied the learning dynamics and emergent recurrent connectivity in a sensory network model based on a first-principle information theoretic approach. As a test case, we applied this framework to a model of a hypercolumn in the visual cortex and found that the evolved connections between orientation columns have a "Mexican hat" profile, consistent with empirical data and previous modeling work. Furthermore, we found that optimal information representation is achieved when the network operates near a critical point in its dynamics. Neuronal networks working near such a phase transition are most sensitive to their inputs and are thus optimal in terms of information representation. Nevertheless, a mild change in the pattern of interactions may cause such networks to undergo a transition into a different regime of behavior in which the network activity is dominated by its internal recurrent dynamics and does not reflect the objective input. We discuss several mechanisms by which the pattern of interactions can be driven into this supercritical regime and relate them to various neurological and neuropsychiatric phenomena.

Near-Critical Dynamics in Stimulus-Evoked Activity of the Human Brain and Its Relation to Spontaneous Resting-State Activity

In recent years, numerous studies have found that the brain at resting state displays many features characteristic of a critical state. Here we examine whether stimulus-evoked activity can also be regarded as critical. Additionally, we investigate the relation between resting-state activity and stimulus-evoked activity from the perspective of criticality. We found that cortical activity measured by magnetoencephalography (MEG) is near critical and organizes as neuronal avalanches at both resting-state and stimulus-evoked activities. Moreover, a significantly high intrasubject similarity between avalanche size and duration distributions at both cognitive states was found, suggesting that the distributions capture specific features of the individual brain dynamics. When comparing different subjects, a higher intersubject consistency was found for stimulus-evoked activity than for resting state. This was expressed by the distance between avalanche size and duration distributions of different participants and was supported by the spatial spreading of the avalanches involved. During the course of stimulus-evoked activity, time locked to the stimulus onset, we demonstrate fluctuations in the gain of the neuronal system and thus short timescale deviations from the critical state. Nonetheless, the overall near-critical state in stimulus-evoked activity is retained over longer timescales, in close proximity and with a high correlation to spontaneous (not time-locked) resting-state activity. Spatially, the observed fluctuations in gain manifest through anticorrelative activations of brain sites involved, suggesting a switch between task-negative (default mode) and task-positive networks and assigning the changes in excitation-inhibition balance to nodes within these networks. Overall, this study offers a novel outlook on evoked activity through the framework of criticality.

SIGNIFICANCE STATEMENT

The organization of stimulus-evoked activity and ongoing cortical activity is a topic of high importance. The article addresses several general questions. What is the spatiotemporal organization of stimulus-evoked cortical activity in healthy human subjects? Are there deviations from excitation-inhibition balance during stimulus-evoked activity? What is the relationship between stimulus-evoked activity and ongoing resting-state activity? Using magnetoencephalography (MEG), we demonstrate that stimulus-evoked activity in humans follows a critical branching process that produces neuronal avalanches. Additionally, we investigate the spatiotemporal relationship between resting-state activity and stimulus-evoked activity from the perspective of critical dynamics. These analyses reveal new aspects of this complex relationship and offer novel insights into the interplay between excitation and inhibition that were not observed previously using conventional approaches.

Burst predicting neurons survive an in vitro glutamate injury model of cerebral ischemia

Neuronal activity in vitro exhibits network bursts characterized by brief periods of increased spike rates. Recent work shows that a subpopulation of neurons reliably predicts the occurrence of network bursts. Here, we examined the role of burst predictors in cultures undergoing an in vitro model of cerebral ischemia. Dissociated primary cortical neurons were plated on multielectrode arrays and spontaneous activity was recorded at 17 days in vitro (DIV). This activity was characterized by neuronal avalanches where burst statistics followed a power law. We identified burst predictors as channels that consistently fired immediately prior to network bursts. The timing of these predictors relative to bursts followed a skewed distribution that differed sharply from a null model based on branching ratio. A portion of cultures were subjected to an excitotoxic insult (DIV 18). Propidium iodine and fluorescence imaging confirmed cell death in these cultures. While the insult did not alter the distribution of avalanches, it resulted in alterations in overall spike rates. Burst predictors, however, maintained baseline levels of activity. The resilience of burst predictors following excitotoxic insult suggests a key role of these units in maintaining network activity following injury, with implications for the selective effects of ischemia in the brain.

Structural versus dynamical origins of mean-field behavior in a self-organized critical model of neuronal avalanches

Critical dynamics of cortical neurons have been intensively studied over the past decade. Neuronal avalanches provide the main experimental as well as theoretical tools to consider criticality in such systems. Experimental studies show that critical neuronal avalanches show mean-field behavior. There are structural as well as recently proposed [Phys. Rev. E 89, 052139 (2014)] dynamical mechanisms that can lead to mean-field behavior. In this work we consider a simple model of neuronal dynamics based on threshold self-organized critical models with synaptic noise. We investigate the role of high-average connectivity, random long-range connections, as well as synaptic noise in achieving mean-field behavior. We employ finite-size scaling in order to extract critical exponents with good accuracy. We conclude that relevant structural mechanisms responsible for mean-field behavior cannot be justified in realistic models of the cortex. However, strong dynamical noise, which can have realistic justifications, always leads to mean-field behavior regardless of the underlying structure. Our work provides a different (dynamical) origin than the conventionally accepted (structural) mechanisms for mean-field behavior in neuronal avalanches.

Local inhibitory plasticity tunes macroscopic brain dynamics and allows the emergence of functional brain networks

Rich, spontaneous brain activity has been observed across a range of different temporal and spatial scales. These dynamics are thought to be important for efficient neural functioning. A range of experimental evidence suggests that these neural dynamics are maintained across a variety of different cognitive states, in response to alterations of the environment and to changes in brain configuration (e.g., across individuals, development and in many neurological disorders). This suggests that the brain has evolved mechanisms to maintain rich dynamics across a broad range of situations. Several mechanisms based around homeostatic plasticity have been proposed to explain how these dynamics emerge from networks of neurons at the microscopic scale. Here we explore how a homeostatic mechanism may operate at the macroscopic scale: in particular, focusing on how it interacts with the underlying structural network topology and how it gives rise to well-described functional connectivity networks. We use a simple mean-field model of the brain, constrained by empirical white matter structural connectivity where each region of the brain is simulated using a pool of excitatory and inhibitory neurons. We show, as with the microscopic work, that homeostatic plasticity regulates network activity and allows for the emergence of rich, spontaneous dynamics across a range of brain configurations, which otherwise show a very limited range of dynamic regimes. In addition, the simulated functional connectivity of the homeostatic model better resembles empirical functional connectivity network. To accomplish this, we show how the inhibitory weights adapt over time to capture important graph theoretic properties of the underlying structural network. Therefore, this work presents suggests how inhibitory homeostatic mechanisms facilitate stable macroscopic dynamics to emerge in the brain, aiding the formation of functional connectivity networks.

Cognitive Flexibility through Metastable Neural Dynamics Is Disrupted by Damage to the Structural Connectome

Current theory proposes that healthy neural dynamics operate in a metastable regime, where brain regions interact to simultaneously maximize integration and segregation. Metastability may confer important behavioral properties, such as cognitive flexibility. It is increasingly recognized that neural dynamics are constrained by the underlying structural connections between brain regions. An important challenge is, therefore, to relate structural connectivity, neural dynamics, and behavior. Traumatic brain injury (TBI) is a pre-eminent structural disconnection disorder whereby traumatic axonal injury damages large-scale connectivity, producing characteristic cognitive impairments, including slowed information processing speed and reduced cognitive flexibility, that may be a result of disrupted metastable dynamics. Therefore, TBI provides an experimental and theoretical model to examine how metastable dynamics relate to structural connectivity and cognition. Here, we use complementary empirical and computational approaches to investigate how metastability arises from the healthy structural connectome and relates to cognitive performance. We found reduced metastability in large-scale neural dynamics after TBI, measured with resting-state functional MRI. This reduction in metastability was associated with damage to the connectome, measured using diffusion MRI. Furthermore, decreased metastability was associated with reduced cognitive flexibility and information processing. A computational model, defined by empirically derived connectivity data, demonstrates how behaviorally relevant changes in neural dynamics result from structural disconnection. Our findings suggest how metastable dynamics are important for normal brain function and contingent on the structure of the human connectome.

Importance of being Nernst: Synaptic activity and functional relevance in stem cell-derived neurons

Functional synaptogenesis and network emergence are signature endpoints of neurogenesis. These behaviors provide higher-order confirmation that biochemical and cellular processes necessary for neurotransmitter release, post-synaptic detection and network propagation of neuronal activity have been properly expressed and coordinated among cells. The development of synaptic neurotransmission can therefore be considered a defining property of neurons. Although dissociated primary neuron cultures readily form functioning synapses and network behaviors in vitro, continuously cultured neurogenic cell lines have historically failed to meet these criteria. Therefore, in vitro-derived neuron models that develop synaptic transmission are critically needed for a wide array of studies, including molecular neuroscience, developmental neurogenesis, disease research and neurotoxicology. Over the last decade, neurons derived from various stem cell lines have shown varying ability to develop into functionally mature neurons. In this review, we will discuss the neurogenic potential of various stem cells populations, addressing strengths and weaknesses of each, with particular attention to the emergence of functional behaviors. We will propose methods to functionally characterize new stem cell-derived neuron (SCN) platforms to improve their reliability as physiological relevant models. Finally, we will review how synaptically active SCNs can be applied to accelerate research in a variety of areas. Ultimately, emphasizing the critical importance of synaptic activity and network responses as a marker of neuronal maturation is anticipated to result in in vitro findings that better translate to efficacious clinical treatments.

Degree Correlations Optimize Neuronal Network Sensitivity to Sub-Threshold Stimuli

Information processing in the brain crucially depends on the topology of the neuronal connections. We investigate how the topology influences the response of a population of leaky integrate-and-fire neurons to a stimulus. We devise a method to calculate firing rates from a self-consistent system of equations taking into account the degree distribution and degree correlations in the network. We show that assortative degree correlations strongly improve the sensitivity for weak stimuli and propose that such networks possess an advantage in signal processing. We moreover find that there exists an optimum in assortativity at an intermediate level leading to a maximum in input/output mutual information.

Relationship of fast- and slow-timescale neuronal dynamics in human MEG and SEEG

A growing body of evidence suggests that the neuronal dynamics are poised at criticality. Neuronal avalanches and long-range temporal correlations (LRTCs) are hallmarks of such critical dynamics in neuronal activity and occur at fast (subsecond) and slow (seconds to hours) timescales, respectively. The critical dynamics at different timescales can be characterized by their power-law scaling exponents. However, insight into the avalanche dynamics and LRTCs in the human brain has been largely obtained with sensor-level MEG and EEG recordings, which yield only limited anatomical insight and results confounded by signal mixing. We investigated here the relationship between the human neuronal dynamics at fast and slow timescales using both source-reconstructed MEG and intracranial stereotactical electroencephalography (SEEG). Both MEG and SEEG revealed avalanche dynamics that were characterized parameter-dependently by power-law or truncated-power-law size distributions. Both methods also revealed robust LRTCs throughout the neocortex with distinct scaling exponents in different functional brain systems and frequency bands. The exponents of power-law regimen neuronal avalanches and LRTCs were strongly correlated across subjects. Qualitatively similar power-law correlations were also observed in surrogate data without spatial correlations but with scaling exponents distinct from those of original data. Furthermore, we found that LRTCs in the autonomous nervous system, as indexed by heart-rate variability, were correlated in a complex manner with cortical neuronal avalanches and LRTCs in MEG but not SEEG. These scalp and intracranial data hence show that power-law scaling behavior is a pervasive but neuroanatomically inhomogeneous property of neuronal dynamics in central and autonomous nervous systems.

Fractal analysis of MRI data for the characterization of patients with schizophrenia and bipolar disorder

Fractal geometry can be used to analyze shape and patterns in brain images. With this study we use fractals to analyze T1 data of patients affected by schizophrenia or bipolar disorder, with the aim of distinguishing between healthy and pathological brains using the complexity of brain structure, in particular of grey matter, as a marker of disease. 39 healthy volunteers, 25 subjects affected by schizophrenia and 11 patients affected by bipolar disorder underwent an MRI session. We evaluated fractal dimension of the brain cortex and its substructures, calculated with an algorithm based on the box-count algorithm. We modified this algorithm, with the aim of avoiding the segmentation processing step and using all the information stored in the image grey levels. Moreover, to increase sensitivity to local structural changes, we computed a value of fractal dimension for each slice of the brain or of the particular structure. To have reference values in comparing healthy subjects with patients, we built a template by averaging fractal dimension values of the healthy volunteers data. Standard deviation was evaluated and used to create a confidence interval. We also performed a slice by slice t-test to assess the difference at slice level between the three groups. Consistent average fractal dimension values were found across all the structures in healthy controls, while in the pathological groups we found consistent differences, indicating a change in brain and structures complexity induced by these disorders.

Spin-glass model predicts metastable brain states that diminish in anesthesia

Patterns of resting state connectivity change dynamically and may represent modes of cognitive information processing. The diversity of connectivity patterns (global brain states) reflects the information capacity of the brain and determines the state of consciousness. In this work, computer simulation was used to explore the repertoire of global brain states as a function of cortical activation level. We implemented a modified spin glass model to describe UP/DOWN state transitions of neuronal populations at a mesoscopic scale based on resting state BOLD fMRI data. Resting state fMRI was recorded in 20 participants and mapped to 10,000 cortical regions (sites) defined on a group-aligned cortical surface map. Each site represented the population activity of a ~20 mm(2) area of the cortex. Cross-correlation matrices of the mapped BOLD time courses of the set of sites were calculated and averaged across subjects. In the model, each cortical site was allowed to interact with the 16 other sites that had the highest pair-wise correlation values. All sites stochastically transitioned between UP and DOWN states under the net influence of their 16 pairs. The probability of local state transitions was controlled by a single parameter T corresponding to the level of global cortical activation. To estimate the number of distinct global states, first we ran 10,000 simulations at T = 0. Simulations were started from random configurations that converged to one of several distinct patterns. Using hierarchical clustering, at 99% similarity, close to 300 distinct states were found. At intermediate T, metastable state configurations were formed suggesting critical behavior with a sharp increase in the number of metastable states at an optimal T. Both reduced activation (anesthesia, sleep) and increased activation (hyper-activation) moved the system away from equilibrium, presumably incompatible with conscious mentation. During equilibrium, the diversity of large-scale brain states was maximum, compatible with maximum information capacity-a presumed condition of consciousness.

Quasicritical brain dynamics on a nonequilibrium Widom line

Is the brain really operating at a critical point? We study the nonequilibrium properties of a neural network which models the dynamics of the neocortex and argue for optimal quasicritical dynamics on the Widom line where the correlation length and information transmission are optimized. We simulate the network and introduce an analytical mean-field approximation, characterize the nonequilibrium phase transitions, and present a nonequilibrium phase diagram, which shows that in addition to an ordered and disordered phase, the system exhibits a "quasiperiodic" phase corresponding to synchronous activity in simulations, which may be related to the pathological synchronization associated with epilepsy.

A reaction-diffusion model of cholinergic retinal waves

Prior to receiving visual stimuli, spontaneous, correlated activity in the retina, called retinal waves, drives activity-dependent developmental programs. Early-stage waves mediated by acetylcholine (ACh) manifest as slow, spreading bursts of action potentials. They are believed to be initiated by the spontaneous firing of Starburst Amacrine Cells (SACs), whose dense, recurrent connectivity then propagates this activity laterally. Their inter-wave interval and shifting wave boundaries are the result of the slow after-hyperpolarization of the SACs creating an evolving mosaic of recruitable and refractory cells, which can and cannot participate in waves, respectively. Recent evidence suggests that cholinergic waves may be modulated by the extracellular concentration of ACh. Here, we construct a simplified, biophysically consistent, reaction-diffusion model of cholinergic retinal waves capable of recapitulating wave dynamics observed in mice retina recordings. The dense, recurrent connectivity of SACs is modeled through local, excitatory coupling occurring via the volume release and diffusion of ACh. In addition to simulation, we are thus able to use non-linear wave theory to connect wave features to underlying physiological parameters, making the model useful in determining appropriate pharmacological manipulations to experimentally produce waves of a prescribed spatiotemporal character. The model is used to determine how ACh mediated connectivity may modulate wave activity, and how parameters such as the spontaneous activation rate and sAHP refractory period contribute to critical wave size variability.

Pathologies in functional connectivity, feedback control and robustness: a global workspace perspective on autism spectrum disorders

We study the background to problems of functional connectivity in autism spectrum disorders within the neurocognitive framework of the global workspace model. This we proceed to do by observing network irregularities detracting from that of a well-formed small world network architecture. This is discussed in terms of pathologies in functional connectivity and lack of central coherence disrupting inter-network communication thus impairing effective cognitive action. A typical coherence-connectivity measure as a by-product of various neuroimaging results is considered. This is related to a model of feedback control in which a coherence function in the frequency domain is modified by an environmentally determined interaction parameter. With respect to the latter, we discuss the stability question that in theory may counterbalance inessential metabolic costs and incoherence of processing. We suggest that factors such as local overconnectivity and global underconnectivity, along with acute over-expenditure of metabolic costs give rise to instability within the connective core of the workspace.

Self-organized criticality as a fundamental property of neural systems

The neural criticality hypothesis states that the brain may be poised in a critical state at a boundary between different types of dynamics. Theoretical and experimental studies show that critical systems often exhibit optimal computational properties, suggesting the possibility that criticality has been evolutionarily selected as a useful trait for our nervous system. Evidence for criticality has been found in cell cultures, brain slices, and anesthetized animals. Yet, inconsistent results were reported for recordings in awake animals and humans, and current results point to open questions about the exact nature and mechanism of criticality, as well as its functional role. Therefore, the criticality hypothesis has remained a controversial proposition. Here, we provide an account of the mathematical and physical foundations of criticality. In the light of this conceptual framework, we then review and discuss recent experimental studies with the aim of identifying important next steps to be taken and connections to other fields that should be explored.

Marginally subcritical dynamics explain enhanced stimulus discriminability under attention

Recent experimental and theoretical work has established the hypothesis that cortical neurons operate close to a critical state which describes a phase transition from chaotic to ordered dynamics. Critical dynamics are suggested to optimize several aspects of neuronal information processing. However, although critical dynamics have been demonstrated in recordings of spontaneously active cortical neurons, little is known about how these dynamics are affected by task-dependent changes in neuronal activity when the cortex is engaged in stimulus processing. Here we explore this question in the context of cortical information processing modulated by selective visual attention. In particular, we focus on recent findings that local field potentials (LFPs) in macaque area V4 demonstrate an increase in γ-band synchrony and a simultaneous enhancement of object representation with attention. We reproduce these results using a model of integrate-and-fire neurons where attention increases synchrony by enhancing the efficacy of recurrent interactions. In the phase space spanned by excitatory and inhibitory coupling strengths, we identify critical points and regions of enhanced discriminability. Furthermore, we quantify encoding capacity using information entropy. We find a rapid enhancement of stimulus discriminability with the emergence of synchrony in the network. Strikingly, only a narrow region in the phase space, at the transition from subcritical to supercritical dynamics, supports the experimentally observed discriminability increase. At the supercritical border of this transition region, information entropy decreases drastically as synchrony sets in. At the subcritical border, entropy is maximized under the assumption of a coarse observation scale. Our results suggest that cortical networks operate at such near-critical states, allowing minimal attentional modulations of network excitability to substantially augment stimulus representation in the LFPs.

Zipf's law and criticality in multivariate data without fine-tuning

The joint probability distribution of states of many degrees of freedom in biological systems, such as firing patterns in neural networks or antibody sequence compositions, often follows Zipf's law, where a power law is observed on a rank-frequency plot. This behavior has been shown to imply that these systems reside near a unique critical point where the extensive parts of the entropy and energy are exactly equal. Here, we show analytically, and via numerical simulations, that Zipf-like probability distributions arise naturally if there is a fluctuating unobserved variable (or variables) that affects the system, such as a common input stimulus that causes individual neurons to fire at time-varying rates. In statistics and machine learning, these are called latent-variable or mixture models. We show that Zipf's law arises generically for large systems, without fine-tuning parameters to a point. Our work gives insight into the ubiquity of Zipf's law in a wide range of systems.

Frustrated hierarchical synchronization and emergent complexity in the human connectome network

The spontaneous emergence of coherent behavior through synchronization plays a key role in neural function, and its anomalies often lie at the basis of pathologies. Here we employ a parsimonious (mesoscopic) approach to study analytically and computationally the synchronization (Kuramoto) dynamics on the actual human-brain connectome network. We elucidate the existence of a so-far-uncovered intermediate phase, placed between the standard synchronous and asynchronous phases, i.e. between order and disorder. This novel phase stems from the hierarchical modular organization of the connectome. Where one would expect a hierarchical synchronization process, we show that the interplay between structural bottlenecks and quenched intrinsic frequency heterogeneities at many different scales, gives rise to frustrated synchronization, metastability, and chimera-like states, resulting in a very rich and complex phenomenology. We uncover the origin of the dynamic freezing behind these features by using spectral graph theory and discuss how the emerging complex synchronization patterns relate to the need for the brain to access –in a robust though flexible way– a large variety of functional attractors and dynamical repertoires without ad hoc fine-tuning to a critical point.

Perampanel, an antagonist of α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid receptors, for the treatment of epilepsy: studies in human epileptic brain and nonepileptic brain and in rodent models

Perampanel [Fycompa, 2-(2-oxo-1-phenyl-5-pyridin-2-yl-1,2-dihydropyridin-3-yl)benzonitrile hydrate 4:3; Eisai Inc., Woodcliff Lake, NJ] is an AMPA (α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid) receptor antagonist used as an adjunctive treatment of partial-onset seizures. We asked whether perampanel has AMPA receptor antagonist activity in both the cerebral cortex and hippocampus associated with antiepileptic efficacy and also in the cerebellum associated with motor side effects in rodent and human brains. We also asked whether epileptic or nonepileptic human cortex is similarly responsive to AMPA receptor antagonism by perampanel. In rodent models, perampanel decreased epileptic-like activity in multiple seizure models. However, doses of perampanel that had anticonvulsant effects were within the same range as those engendering motor side effects. Perampanel inhibited native rat and human AMPA receptors from the hippocampus as well as the cerebellum that were reconstituted into Xenopus oocytes. In addition, with the same technique, we found that perampanel inhibited AMPA receptors from hippocampal tissue that had been removed from a patient who underwent surgical resection for refractory epilepsy. Perampanel inhibited AMPA receptor-mediated ion currents from all the tissues investigated with similar potency (IC50 values ranging from 2.6 to 7.0 μM). Cortical slices from the left temporal lobe derived from the same patient were studied in a 60-microelectrode array. Large field potentials were evoked on at least 45 channels of the array, and 10 μM perampanel decreased their amplitude and firing rate. Perampanel also produced a 33% reduction in the branching parameter, demonstrating the effects of perampanel at the network level. These data suggest that perampanel blocks AMPA receptors globally across the brain to account for both its antiepileptic and side-effect profile in rodents and epileptic patients.

Information-based fitness and the emergence of criticality in living systems

Empirical evidence suggesting that living systems might operate in the vicinity of critical points, at the borderline between order and disorder, has proliferated in recent years, with examples ranging from spontaneous brain activity to flock dynamics. However, a well-founded theory for understanding how and why interacting living systems could dynamically tune themselves to be poised in the vicinity of a critical point is lacking. Here we use tools from statistical mechanics and information theory to show that complex adaptive or evolutionary systems can be much more efficient in coping with diverse heterogeneous environmental conditions when operating at criticality. Analytical as well as computational evolutionary and adaptive models vividly illustrate that a community of such systems dynamically self-tunes close to a critical state as the complexity of the environment increases while they remain noncritical for simple and predictable environments. A more robust convergence to criticality emerges in coevolutionary and coadaptive setups in which individuals aim to represent other agents in the community with fidelity, thereby creating a collective critical ensemble and providing the best possible tradeoff between accuracy and flexibility. Our approach provides a parsimonious and general mechanism for the emergence of critical-like behavior in living systems needing to cope with complex environments or trying to efficiently coordinate themselves as an ensemble.

Spike avalanches in vivo suggest a driven, slightly subcritical brain state

In self-organized critical (SOC) systems avalanche size distributions follow power-laws. Power-laws have also been observed for neural activity, and so it has been proposed that SOC underlies brain organization as well. Surprisingly, for spiking activity in vivo, evidence for SOC is still lacking. Therefore, we analyzed highly parallel spike recordings from awake rats and monkeys, anesthetized cats, and also local field potentials from humans. We compared these to spiking activity from two established critical models: the Bak-Tang-Wiesenfeld model, and a stochastic branching model. We found fundamental differences between the neural and the model activity. These differences could be overcome for both models through a combination of three modifications: (1) subsampling, (2) increasing the input to the model (this way eliminating the separation of time scales, which is fundamental to SOC and its avalanche definition), and (3) making the model slightly sub-critical. The match between the neural activity and the modified models held not only for the classical avalanche size distributions and estimated branching parameters, but also for two novel measures (mean avalanche size, and frequency of single spikes), and for the dependence of all these measures on the temporal bin size. Our results suggest that neural activity in vivo shows a mélange of avalanches, and not temporally separated ones, and that their global activity propagation can be approximated by the principle that one spike on average triggers a little less than one spike in the next step. This implies that neural activity does not reflect a SOC state but a slightly sub-critical regime without a separation of time scales. Potential advantages of this regime may be faster information processing, and a safety margin from super-criticality, which has been linked to epilepsy.

Optimized Gamma Synchronization Enhances Functional Binding of Fronto-Parietal Cortices in Mathematically Gifted Adolescents during Deductive Reasoning

As enhanced fronto-parietal network has been suggested to support reasoning ability of math-gifted adolescents, the main goal of this EEG source analysis is to investigate the temporal binding of the gamma-band (30–60 Hz) synchronization between frontal and parietal cortices in adolescents with exceptional mathematical ability, including the functional connectivity of gamma neurocognitive network, the temporal dynamics of fronto-parietal network (phase-locking durations and network lability in time domain), and the self-organized criticality of synchronizing oscillation. Compared with the average-ability subjects, the math-gifted adolescents show a highly integrated fronto-parietal network due to distant gamma phase-locking oscillations, which is indicated by lower modularity of the global network topology, more “connector bridges” between the frontal and parietal cortices and less “connector hubs” in the sensorimotor cortex. The time domain analysis finds that, while maintaining more stable phase dynamics of the fronto-parietal coupling, the math-gifted adolescents are characterized by more extensive fronto-parietal connection reconfiguration. The results from sample fitting in the power-law model further find that the phase-locking durations in the math-gifted brain abides by a wider interval of the power-law distribution. This phase-lock distribution mechanism could represent a relatively optimized pattern for the functional binding of frontal–parietal network, which underlies stable fronto-parietal connectivity and increases flexibility of timely network reconfiguration.

Power law scaling in synchronization of brain signals depends on cognitive load

As it has several features that optimize information processing, it has been proposed that criticality governs the dynamics of nervous system activity. Indications of such dynamics have been reported for a variety of in vitro and in vivo recordings, ranging from in vitro slice electrophysiology to human functional magnetic resonance imaging. However, there still remains considerable debate as to whether the brain actually operates close to criticality or in another governing state such as stochastic or oscillatory dynamics. A tool used to investigate the criticality of nervous system data is the inspection of power-law distributions. Although the findings are controversial, such power-law scaling has been found in different types of recordings. Here, we studied whether there is a power law scaling in the distribution of the phase synchronization derived from magnetoencephalographic recordings during executive function tasks performed by children with and without autism. Characterizing the brain dynamics that is different between autistic and non-autistic individuals is important in order to find differences that could either aid diagnosis or provide insights as to possible therapeutic interventions in autism. We report in this study that power law scaling in the distributions of a phase synchrony index is not very common and its frequency of occurrence is similar in the control and the autism group. In addition, power law scaling tends to diminish with increased cognitive load (difficulty or engagement in the task). There were indications of changes in the probability distribution functions for the phase synchrony that were associated with a transition from power law scaling to lack of power law (or vice versa), which suggests the presence of phenomenological bifurcations in brain dynamics associated with cognitive load. Hence, brain dynamics may fluctuate between criticality and other regimes depending upon context and behaviors.

Scale-free functional connectivity of the brain is maintained in anesthetized healthy participants but not in patients with unresponsive wakefulness syndrome

Loss of consciousness in anesthetized healthy participants and in patients with unresponsive wakefulness syndrome (UWS) is associated with substantial alterations of functional connectivity across large-scale brain networks. Yet, a prominent distinction between the two cases is that after anesthesia, brain connectivity and consciousness are spontaneously restored, whereas in patients with UWS this restoration fails to occur, but why? A possible explanation is that the self-organizing capability of the brain is compromised in patients with UWS but not in healthy participants undergoing anesthesia. According to the theory of self-organized criticality, many natural complex systems, including the brain, evolve spontaneously to a critical state wherein system behaviors display spatial and/or temporal scale-invariant characteristics. Here we tested the hypothesis that the scale-free property of brain network organization is in fact fundamentally different between anesthetized healthy participants and UWS patients. We introduced a novel, computationally efficient approach to determine anatomical-functional parcellation of the whole-brain network at increasingly finer spatial scales. We found that in healthy participants, scale-free distributions of node size and node degree were present across wakefulness, propofol sedation, and recovery, despite significant propofol-induced functional connectivity changes. In patients with UWS, the scale-free distribution of node degree was absent, reflecting a fundamental difference between the two groups in adaptive reconfiguration of functional interaction between network components. The maintenance of scale-invariance across propofol sedation in healthy participants suggests the presence of persistent, on-going self-organizing processes to a critical state – a capacity that is compromised in patients with UWS.

The control of global brain dynamics: opposing actions of frontoparietal control and default mode networks on attention

Understanding how dynamic changes in brain activity control behavior is a major challenge of cognitive neuroscience. Here, we consider the brain as a complex dynamic system and define two measures of brain dynamics: the synchrony of brain activity, measured by the spatial coherence of the BOLD signal across regions of the brain; and metastability, which we define as the extent to which synchrony varies over time. We investigate the relationship among brain network activity, metastability, and cognitive state in humans, testing the hypothesis that global metastability is "tuned" by network interactions. We study the following two conditions: (1) an attentionally demanding choice reaction time task (CRT); and (2) an unconstrained "rest" state. Functional MRI demonstrated increased synchrony, and decreased metastability was associated with increased activity within the frontoparietal control/dorsal attention network (FPCN/DAN) activity and decreased default mode network (DMN) activity during the CRT compared with rest. Using a computational model of neural dynamics that is constrained by white matter structure to test whether simulated changes in FPCN/DAN and DMN activity produce similar effects, we demonstate that activation of the FPCN/DAN increases global synchrony and decreases metastability. DMN activation had the opposite effects. These results suggest that the balance of activity in the FPCN/DAN and DMN might control global metastability, providing a mechanistic explanation of how attentional state is shifted between an unfocused/exploratory mode characterized by high metastability, and a focused/constrained mode characterized by low metastability.

Optimal system size for complex dynamics in random neural networks near criticality

In this article, we consider a model of dynamical agents coupled through a random connectivity matrix, as introduced by Sompolinsky et al. [Phys. Rev. Lett. 61(3), 259-262 (1988)] in the context of random neural networks. When system size is infinite, it is known that increasing the disorder parameter induces a phase transition leading to chaotic dynamics. We observe and investigate here a novel phenomenon in the sub-critical regime for finite size systems: the probability of observing complex dynamics is maximal for an intermediate system size when the disorder is close enough to criticality. We give a more general explanation of this type of system size resonance in the framework of extreme values theory for eigenvalues of random matrices.

Dynamical correlation patterns and corresponding community structure in neural spontaneous activity at criticality

It has been considered that the state in the vicinity of a critical point, which is the point between ordered and disordered states, can underlie and facilitate information processing of the brain in various aspects. In this research, we numerically study the influence of criticality on one aspect of brain information processing, i.e., the community structure, which is an important characteristic of complex networks. We examine community structure of the functional connectivity in simulated brain spontaneous activity, which is based on dynamical correlations between neural activity patterns at different positions. The brain spontaneous activity is simulated by a neural field model whose parameter covers subcritical, critical, and supercritical regions. Then, the corresponding dynamical correlation patterns and community structure are compared. In the critical region, we found some distinctive properties, namely high correlation and correlation switching, high modularity and a low number of modules, high stability of the dynamical functional connectivity, and moderate flexibility of the community structure across temporal scales. We also discuss how these characteristics might improve information processing of the brain.

Short-term synaptic plasticity in the deterministic Tsodyks-Markram model leads to unpredictable network dynamics

Short-term synaptic plasticity strongly affects the neural dynamics of cortical networks. The Tsodyks and Markram (TM) model for short-term synaptic plasticity accurately accounts for a wide range of physiological responses at different types of cortical synapses. Here, we report a route to chaotic behavior via a Shilnikov homoclinic bifurcation that dynamically organizes some of the responses in the TM model. In particular, the presence of such a homoclinic bifurcation strongly affects the shape of the trajectories in the phase space and induces highly irregular transient dynamics; indeed, in the vicinity of the Shilnikov homoclinic bifurcation, the number of population spikes and their precise timing are unpredictable and highly sensitive to the initial conditions. Such an irregular deterministic dynamics has its counterpart in stochastic/network versions of the TM model: The existence of the Shilnikov homoclinic bifurcation generates complex and irregular spiking patterns and--acting as a sort of springboard--facilitates transitions between the down-state and unstable periodic orbits. The interplay between the (deterministic) homoclinic bifurcation and stochastic effects may give rise to some of the complex dynamics observed in neural systems.

Cerebral energy metabolism and the brain's functional network architecture: an integrative review

Recent functional magnetic resonance imaging (fMRI) studies have emphasized the contributions of synchronized activity in distributed brain networks to cognitive processes in both health and disease. The brain's 'functional connectivity' is typically estimated from correlations in the activity time series of anatomically remote areas, and postulated to reflect information flow between neuronal populations. Although the topological properties of functional brain networks have been studied extensively, considerably less is known regarding the neurophysiological and biochemical factors underlying the temporal coordination of large neuronal ensembles. In this review, we highlight the critical contributions of high-frequency electrical oscillations in the γ-band (30 to 100 Hz) to the emergence of functional brain networks. After describing the neurobiological substrates of γ-band dynamics, we specifically discuss the elevated energy requirements of high-frequency neural oscillations, which represent a mechanistic link between the functional connectivity of brain regions and their respective metabolic demands. Experimental evidence is presented for the high oxygen and glucose consumption, and strong mitochondrial performance required to support rhythmic cortical activity in the γ-band. Finally, the implications of mitochondrial impairments and deficits in glucose metabolism for cognition and behavior are discussed in the context of neuropsychiatric and neurodegenerative syndromes characterized by large-scale changes in the organization of functional brain networks.

The evolution of intelligent developmental systems

This chapter aims to understand the relations between the evolution and development of complex cognitive functions by emphasizing the context of complex, changeable environments. What evolves and develops in such contexts cannot be achieved by linear deterministic processes based on stable "codes". Rather, what is needed, even in the molecular ensembles of single-cell organisms, are "intelligent" systems with nonlinear dynamic processing, sensitive to informational structures, not just elements, in environments. This is the view emerging in recent molecular biology. The research is also constructing a new "biologic" of both evolution and development, providing a clearer rationale for transitions into more complex forms, including epigenetic, physiological, nervous, cognitive, and human sociocognitive forms. This chapter explains how these transitions form a nested hierarchical system in which the dynamics within and between levels creates emergent abilities so often underestimated or even demeaned in previous accounts, especially regarding human cognition. The implications of the view for human development in modern societies are also briefly considered.

Optimal channel efficiency in a sensory network

Spontaneous neural activity has been increasingly recognized as a subject of key relevance in neuroscience. It exhibits nontrivial spatiotemporal structure reflecting the organization of the underlying neural network and has proved to be closely intertwined with stimulus-induced activity patterns. As an additional contribution in this regard, we report computational studies that strongly suggest that a stimulus-free feature rules the behavior of an important psychophysical measure of the sensibility of a sensory system to a stimulus, the so-called dynamic range. Indeed in this paper we show that the entropy of the distribution of avalanche lifetimes (information efficiency, since it can be interpreted as the efficiency of the network seen as a communication channel) always accompanies the dynamic range in the benchmark model for sensory systems. Specifically, by simulating the Kinouchi-Copelli (KC) model on two broad families of model networks, we generically observed that both quantities always increase or decrease together as functions of the average branching ratio (the control parameter of the KC model) and that the information efficiency typically exhibits critical optimization jointly with the dynamic range (i.e., both quantities are optimized at the same value of that control parameter, that turns out to be the critical point of a nonequilibrium phase transition). In contrast with the practice of taking power laws to identify critical points in most studies describing measured neuronal avalanches, we rely on data collapses as more robust signatures of criticality to claim that critical optimization may happen even when the distribution of avalanche lifetimes is not a power law, as suggested by a recent experiment. Finally, we note that the entropy of the size distribution of avalanches (information capacity) does not always follow the dynamic range and the information efficiency when they are critically optimized, despite being more widely used than the latter to describe the computational capabilities of a neural network. This strongly suggests that dynamical rules allowing a proper temporal matching of the states of the interacting neurons is the key for achieving good performance in information processing, rather than increasing the number of available units.

Development of antipsychotic medications with novel mechanisms of action based on computational modeling of hippocampal neuropathology

A large number of cellular level abnormalities have been identified in the hippocampus of schizophrenic subjects. Nonetheless, it remains uncertain how these pathologies interact at a system level to create clinical symptoms, and this has hindered the development of more effective antipsychotic medications. Using a 72-processor supercomputer, we created a tissue level hippocampal simulation, featuring multicompartmental neuron models with multiple ion channel subtypes and synaptic channels with realistic temporal dynamics. As an index of the schizophrenic phenotype, we used the specific inability of the model to attune to 40 Hz (gamma band) stimulation, a well-characterized abnormality in schizophrenia. We examined several possible combinations of putatively schizophrenogenic cellular lesions by systematically varying model parameters representing NMDA channel function, dendritic spine density, and GABA system integrity, conducting 910 trials in total. Two discrete “clusters” of neuropathological changes were identified. The most robust was characterized by co-occurring modest reductions in NMDA system function (-30%) and dendritic spine density (-30%). Another set of lesions had greater NMDA hypofunction along with low level GABA system dysregulation. To the schizophrenic model, we applied the effects of 1,500 virtual medications, which were implemented by varying five model parameters, independently, in a graded manner; the effects of known drugs were also applied. The simulation accurately distinguished agents that are known to lack clinical efficacy, and identified novel mechanisms (e.g., decrease in AMPA conductance decay time constant, increase in projection strength of calretinin-positive interneurons) and combinations of mechanisms that could re-equilibrate model behavior. These findings shed light on the mechanistic links between schizophrenic neuropathology and the gamma band oscillatory abnormalities observed in the illness. As such, they generate specific falsifiable hypotheses, which can guide postmortem and other laboratory research. Significantly, this work also suggests specific non-obvious targets for potential pharmacologic agents.

Topological and dynamical complexity of random neural networks

Random neural networks are dynamical descriptions of randomly interconnected neural units. These show a phase transition to chaos as a disorder parameter is increased. The microscopic mechanisms underlying this phase transition are unknown and, similar to spin glasses, shall be fundamentally related to the behavior of the system. In this Letter, we investigate the explosion of complexity arising near that phase transition. We show that the mean number of equilibria undergoes a sharp transition from one equilibrium to a very large number scaling exponentially with the dimension on the system. Near criticality, we compute the exponential rate of divergence, called topological complexity. Strikingly, we show that it behaves exactly as the maximal Lyapunov exponent, a classical measure of dynamical complexity. This relationship unravels a microscopic mechanism leading to chaos which we further demonstrate on a simpler disordered system, suggesting a deep and underexplored link between topological and dynamical complexity.

Hierarchical networks, power laws, and neuronal avalanches

We show that in networks with a hierarchical architecture, critical dynamical behaviors can emerge even when the underlying dynamical processes are not critical. This finding provides explicit insight into current studies of the brain's neuronal network showing power-law avalanches in neural recordings, and provides a theoretical justification of recent numerical findings. Our analysis shows how the hierarchical organization of a network can itself lead to power-law distributions of avalanche sizes and durations, scaling laws between anomalous exponents, and universal functions-even in the absence of self-organized criticality or critical points. This hierarchy-induced phenomenon is independent of, though can potentially operate in conjunction with, standard dynamical mechanisms for generating power laws.

Developing neuronal networks: self-organized criticality predicts the future

Self-organized criticality emerged in neural activity is one of the key concepts to describe the formation and the function of developing neuronal networks. The relationship between critical dynamics and neural development is both theoretically and experimentally appealing. However, whereas it is well-known that cortical networks exhibit a rich repertoire of activity patterns at different stages during in vitro maturation, dynamical activity patterns through the entire neural development still remains unclear. Here we show that a series of metastable network states emerged in the developing and “aging” process of hippocampal networks cultured from dissociated rat neurons. The unidirectional sequence of state transitions could be only observed in networks showing power-law scaling of distributed neuronal avalanches. Our data suggest that self-organized criticality may guide spontaneous activity into a sequential succession of homeostatically-regulated transient patterns during development, which may help to predict the tendency of neural development at early ages in the future.

Dissipation of 'dark energy' by cortex in knowledge retrieval

We have devised a thermodynamic model of cortical neurodynamics expressed at the classical level by neural networks and at the quantum level by dissipative quantum field theory. Our model is based on features in the spatial images of cortical activity newly revealed by high-density electrode arrays. We have incorporated the mechanism and necessity for so-called dark energy in knowledge retrieval. We have extended the model first using the Carnot cycle to define our measures for energy, entropy and temperature, and then using the Rankine cycle to incorporate criticality and phase transitions. We describe the dynamics of two interactive fields of neural activity that express knowledge, one at high and the other at low energy density, and the two operators that create and annihilate the fields. We postulate that the extremely high density of energy sequestered briefly in cortical activity patterns can account for the vividness, richness of associations, and emotional intensity of memories recalled by stimuli.

Neuromodulation: selected approaches and challenges

The brain operates through complex interactions in the flow of information and signal processing within neural networks. The 'wiring' of such networks, being neuronal or glial, can physically and/or functionally go rogue in various pathological states. Neuromodulation, as a multidisciplinary venture, attempts to correct such faulty nets. In this review, selected approaches and challenges in neuromodulation are discussed. The use of water-dispersible carbon nanotubes has been proven effective in the modulation of neurite outgrowth in culture and in aiding regeneration after spinal cord injury in vivo. Studying neural circuits using computational biology and analytical engineering approaches brings to light geometrical mapping of dynamics within neural networks, much needed information for stimulation interventions in medical practice. Indeed, sophisticated desynchronization approaches used for brain stimulation have been successful in coaxing 'misfiring' neuronal circuits to resume productive firing patterns in various human disorders. Devices have been developed for the real-time measurement of various neurotransmitters as well as electrical activity in the human brain during electrical deep brain stimulation. Such devices can establish the dynamics of electrochemical changes in the brain during stimulation. With increasing application of nanomaterials in devices for electrical and chemical recording and stimulating in the brain, the era of cellular, and even intracellular, precision neuromodulation will soon be upon us.

Neuronal avalanches, epileptic quakes and other transient forms of neurodynamics

Power-law behaviors in brain activity in healthy animals, in the form of neuronal avalanches, potentially benefit the computational activities of the brain, including information storage, transmission and processing. In contrast, power-law behaviors associated with seizures, in the form of epileptic quakes, potentially interfere with the brain's computational activities. This review draws attention to the potential roles played by homeostatic mechanisms and multistable time-delayed recurrent inhibitory loops in the generation of power-law phenomena. Moreover, it is suggested that distinctions between health and disease are scale-dependent. In other words, what is abnormal and defines disease it is not the propagation of neural activity but the propagation of activity in a neural population that is large enough to interfere with the normal activities of the brain. From this point of view, epilepsy is a disease that results from a failure of mechanisms, possibly located in part in the cortex itself or in the deep brain nuclei and brainstem, which truncate or otherwise confine the spatiotemporal scales of these power-law phenomena.

Critical-state dynamics of avalanches and oscillations jointly emerge from balanced excitation/inhibition in neuronal networks

Criticality has gained widespread interest in neuroscience as an attractive framework for understanding the character and functional implications of variability in brain activity. The metastability of critical systems maximizes their dynamic range, storage capacity, and computational power. Power-law scaling-a hallmark of criticality-has been observed on different levels, e.g., in the distribution of neuronal avalanches in vitro and in vivo, but also in the decay of temporal correlations in behavioral performance and ongoing oscillations in humans. An unresolved issue is whether power-law scaling on different organizational levels in the brain-and possibly in other hierarchically organized systems-can be related. Here, we show that critical-state dynamics of avalanches and oscillations jointly emerge in a neuronal network model when excitation and inhibition is balanced. The oscillatory activity of the model was qualitatively similar to what is typically observed in recordings of human resting-state MEG. We propose that homeostatic plasticity mechanisms tune this balance in healthy brain networks, and that it is essential for critical behavior on multiple levels of neuronal organization with ensuing functional benefits. Based on our network model, we introduce a concept of multi-level criticality in which power-law scaling can emerge on multiple time scales in oscillating networks.

Analytical investigation of self-organized criticality in neural networks

Dynamical criticality has been shown to enhance information processing in dynamical systems, and there is evidence for self-organized criticality in neural networks. A plausible mechanism for such self-organization is activity-dependent synaptic plasticity. Here, we model neurons as discrete-state nodes on an adaptive network following stochastic dynamics. At a threshold connectivity, this system undergoes a dynamical phase transition at which persistent activity sets in. In a low-dimensional representation of the macroscopic dynamics, this corresponds to a transcritical bifurcation. We show analytically that adding activity-dependent rewiring rules, inspired by homeostatic plasticity, leads to the emergence of an attractive steady state at criticality and present numerical evidence for the system's evolution to such a state.

Synchronisation effects on the behavioural performance and information dynamics of a simulated minimally cognitive robotic agent

Oscillatory activity is ubiquitous in nervous systems, with solid evidence that synchronisation mechanisms underpin cognitive processes. Nevertheless, its informational content and relationship with behaviour are still to be fully understood. In addition, cognitive systems cannot be properly appreciated without taking into account brain-body- environment interactions. In this paper, we developed a model based on the Kuramoto Model of coupled phase oscillators to explore the role of neural synchronisation in the performance of a simulated robotic agent in two different minimally cognitive tasks. We show that there is a statistically significant difference in performance and evolvability depending on the synchronisation regime of the network. In both tasks, a combination of information flow and dynamical analyses show that networks with a definite, but not too strong, propensity for synchronisation are more able to reconfigure, to organise themselves functionally and to adapt to different behavioural conditions. The results highlight the asymmetry of information flow and its behavioural correspondence. Importantly, it also shows that neural synchronisation dynamics, when suitably flexible and reconfigurable, can generate minimally cognitive embodied behaviour.