Modeling neuronal avalanches and long-range temporal correlations at the emergence of collective oscillations: Continuously varying exponents mimic M/EEG results

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.

Optimal input reverberation and homeostatic self-organization toward the edge of synchronization

Transient or partial synchronization can be used to do computations, although a fully synchronized network is sometimes related to the onset of epileptic seizures. Here, we propose a homeostatic mechanism that is capable of maintaining a neuronal network at the edge of a synchronization transition, thereby avoiding the harmful consequences of a fully synchronized network. We model neurons by maps since they are dynamically richer than integrate-and-fire models and more computationally efficient than conductance-based approaches. We first describe the synchronization phase transition of a dense network of neurons with different tonic spiking frequencies coupled by gap junctions. We show that at the transition critical point, inputs optimally reverberate through the network activity through transient synchronization. Then, we introduce a local homeostatic dynamic in the synaptic coupling and show that it produces a robust self-organization toward the edge of this phase transition. We discuss the potential biological consequences of this self-organization process, such as its relation to the Brain Criticality hypothesis, its input processing capacity, and how its malfunction could lead to pathological synchronization and the onset of seizure-like activity.

Low-dimensional criticality embedded in high-dimensional awake brain dynamics

Whether cortical neurons operate in a strongly or weakly correlated dynamical regime determines fundamental information processing capabilities and has fueled decades of debate. We offer a resolution of this debate; we show that two important dynamical regimes, typically considered incompatible, can coexist in the same local cortical circuit by separating them into two different subspaces. In awake mouse motor cortex, we find a low-dimensional subspace with large fluctuations consistent with criticality-a dynamical regime with moderate correlations and multi-scale information capacity and transmission. Orthogonal to this critical subspace, we find a high-dimensional subspace containing a desynchronized dynamical regime, which may optimize input discrimination. The critical subspace is apparent only at long timescales, which explains discrepancies among some previous studies. Using a computational model, we show that the emergence of a low-dimensional critical subspace at large timescales agrees with established theory of critical dynamics. Our results suggest that the cortex leverages its high dimensionality to multiplex dynamical regimes across different subspaces.

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 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.

Bistability and criticality in the stochastic Wilson-Cowan model

We study a stochastic version of the Wilson-Cowan model of neural dynamics, where the response function of neurons grows faster than linearly above the threshold. The model shows a region of parameters where two attractive fixed points of the dynamics exist simultaneously. One fixed point is characterized by lower activity and scale-free critical behavior, while the second fixed point corresponds to a higher (supercritical) persistent activity, with small fluctuations around a mean value. When the number of neurons is not too large, the system can switch between these two different states with a probability depending on the parameters of the network. Along with alternation of states, the model displays a bimodal distribution of the avalanches of activity, with a power-law behavior corresponding to the critical state, and a bump of very large avalanches due to the high-activity supercritical state. The bistability is due to the presence of a first-order (discontinuous) transition in the phase diagram, and the observed critical behavior is connected with the line where the low-activity state becomes unstable (spinodal line).

Scale-free behavioral dynamics directly linked with scale-free cortical dynamics

Naturally occurring body movements and collective neural activity both exhibit complex dynamics, often with scale-free, fractal spatiotemporal structure. Scale-free dynamics of both brain and behavior are important because each is associated with functional benefits to the organism. Despite their similarities, scale-free brain activity and scale-free behavior have been studied separately, without a unified explanation. Here, we show that scale-free dynamics of mouse behavior and neurons in the visual cortex are strongly related. Surprisingly, the scale-free neural activity is limited to specific subsets of neurons, and these scale-free subsets exhibit stochastic winner-take-all competition with other neural subsets. This observation is inconsistent with prevailing theories of scale-free dynamics in neural systems, which stem from the criticality hypothesis. We develop a computational model which incorporates known cell-type-specific circuit structure, explaining our findings with a new type of critical dynamics. Our results establish neural underpinnings of scale-free behavior and clear behavioral relevance of scale-free neural activity.

Power spectrum and critical exponents in the 2D stochastic Wilson-Cowan model

The power spectrum of brain activity is composed by peaks at characteristic frequencies superimposed to a background that decays as a power law of the frequency, [Formula: see text], with an exponent [Formula: see text] close to 1 (pink noise). This exponent is predicted to be connected with the exponent [Formula: see text] related to the scaling of the average size with the duration of avalanches of activity. "Mean field" models of neural dynamics predict exponents [Formula: see text] and [Formula: see text] equal or near 2 at criticality (brown noise), including the simple branching model and the fully-connected stochastic Wilson-Cowan model. We here show that a 2D version of the stochastic Wilson-Cowan model, where neuron connections decay exponentially with the distance, is characterized by exponents [Formula: see text] and [Formula: see text] markedly different from those of mean field, respectively around 1 and 1.3. The exponents [Formula: see text] and [Formula: see text] of avalanche size and duration distributions, equal to 1.5 and 2 in mean field, decrease respectively to [Formula: see text] and [Formula: see text]. This seems to suggest the possibility of a different universality class for the model in finite dimension.

Sampling effects and measurement overlap can bias the inference of neuronal avalanches

To date, it is still impossible to sample the entire mammalian brain with single-neuron precision. This forces one to either use spikes (focusing on few neurons) or to use coarse-sampled activity (averaging over many neurons, e.g. LFP). Naturally, the sampling technique impacts inference about collective properties. Here, we emulate both sampling techniques on a simple spiking model to quantify how they alter observed correlations and signatures of criticality. We describe a general effect: when the inter-electrode distance is small, electrodes sample overlapping regions in space, which increases the correlation between the signals. For coarse-sampled activity, this can produce power-law distributions even for non-critical systems. In contrast, spike recordings do not suffer this particular bias and underlying dynamics can be identified. This may resolve why coarse measures and spikes have produced contradicting results in the past.

Theoretical foundations of studying criticality in the brain

Criticality is hypothesized as a physical mechanism underlying efficient transitions between cortical states and remarkable information-processing capacities in the brain. While considerable evidence generally supports this hypothesis, nonnegligible controversies persist regarding the ubiquity of criticality in neural dynamics and its role in information processing. Validity issues frequently arise during identifying potential brain criticality from empirical data. Moreover, the functional benefits implied by brain criticality are frequently misconceived or unduly generalized. These problems stem from the nontriviality and immaturity of the physical theories that analytically derive brain criticality and the statistic techniques that estimate brain criticality from empirical data. To help solve these problems, we present a systematic review and reformulate the foundations of studying brain criticality, that is, ordinary criticality (OC), quasi-criticality (qC), self-organized criticality (SOC), and self-organized quasi-criticality (SOqC), using the terminology of neuroscience. We offer accessible explanations of the physical theories and statistical techniques of brain criticality, providing step-by-step derivations to characterize neural dynamics as a physical system with avalanches. We summarize error-prone details and existing limitations in brain criticality analysis and suggest possible solutions. Moreover, we present a forward-looking perspective on how optimizing the foundations of studying brain criticality can deepen our understanding of various neuroscience questions.

How critical is brain criticality?

Criticality is the singular state of complex systems poised at the brink of a phase transition between order and randomness. Such systems display remarkable information-processing capabilities, evoking the compelling hypothesis that the brain may itself be critical. This foundational idea is now drawing renewed interest thanks to high-density data and converging cross-disciplinary knowledge. Together, these lines of inquiry have shed light on the intimate link between criticality, computation, and cognition. Here, we review these emerging trends in criticality neuroscience, highlighting new data pertaining to the edge of chaos and near-criticality, and making a case for the distance to criticality as a useful metric for probing cognitive states and mental illness. This unfolding progress in the field contributes to establishing criticality theory as a powerful mechanistic framework for studying emergent function and its efficiency in both biological and artificial neural networks.

Disentangling the critical signatures of neural activity

The critical brain hypothesis has emerged as an attractive framework to understand neuronal activity, but it is still widely debated. In this work, we analyze data from a multi-electrodes array in the rat's cortex and we find that power-law neuronal avalanches satisfying the crackling-noise relation coexist with spatial correlations that display typical features of critical systems. In order to shed a light on the underlying mechanisms at the origin of these signatures of criticality, we introduce a paradigmatic framework with a common stochastic modulation and pairwise linear interactions inferred from our data. We show that in such models power-law avalanches that satisfy the crackling-noise relation emerge as a consequence of the extrinsic modulation, whereas scale-free correlations are solely determined by internal interactions. Moreover, this disentangling is fully captured by the mutual information in the system. Finally, we show that analogous power-law avalanches are found in more realistic models of neural activity as well, suggesting that extrinsic modulation might be a broad mechanism for their generation.

Optimal reinforcement learning near the edge of a synchronization transition

Recent experimental and theoretical studies have indicated that the putative criticality of cortical dynamics may correspond to a synchronization phase transition. The critical dynamics near such a critical point needs further investigation specifically when compared to the critical behavior near the standard absorbing state phase transition. Since the phenomena of learning and self-organized criticality (SOC) at the edge of synchronization transition can emerge jointly in spiking neural networks due to the presence of spike-timing dependent plasticity (STDP), it is tempting to ask the following: what is the relationship between synchronization and learning in neural networks? Further, does learning benefit from SOC at the edge of synchronization transition? In this paper, we intend to address these important issues. Accordingly, we construct a biologically inspired model of a cognitive system which learns to perform stimulus-response tasks. We train this system using a reinforcement learning rule implemented through dopamine-modulated STDP. We find that the system exhibits a continuous transition from synchronous to asynchronous neural oscillations upon increasing the average axonal time delay. We characterize the learning performance of the system and observe that it is optimized near the synchronization transition. We also study neuronal avalanches in the system and provide evidence that optimized learning is achieved in a slightly supercritical state.

Universality, criticality and complexity of information propagation in social media

Statistical laws of information avalanches in social media appear, at least according to existing empirical studies, not robust across systems. As a consequence, radically different processes may represent plausible driving mechanisms for information propagation. Here, we analyze almost one billion time-stamped events collected from several online platforms – including Telegram, Twitter and Weibo – over observation windows longer than ten years, and show that the propagation of information in social media is a universal and critical process. Universality arises from the observation of identical macroscopic patterns across platforms, irrespective of the details of the specific system at hand. Critical behavior is deduced from the power-law distributions, and corresponding hyperscaling relations, characterizing size and duration of avalanches of information. Statistical testing on our data indicates that a mixture of simple and complex contagion characterizes the propagation of information in social media. Data suggest that the complexity of the process is correlated with the semantic content of the information that is propagated.

The authors identify characteristic patterns that describe the propagation of information in online social media platforms. They show that, depending on the topic, the information flows can spread as simple or complex contagion processes, operating at a critical regime.

Criticality enhances the multilevel reliability of stimulus responses in cortical neural networks

Cortical neural networks exhibit high internal variability in spontaneous dynamic activities and they can robustly and reliably respond to external stimuli with multilevel features–from microscopic irregular spiking of neurons to macroscopic oscillatory local field potential. A comprehensive study integrating these multilevel features in spontaneous and stimulus–evoked dynamics with seemingly distinct mechanisms is still lacking. Here, we study the stimulus–response dynamics of biologically plausible excitation–inhibition (E–I) balanced networks. We confirm that networks around critical synchronous transition states can maintain strong internal variability but are sensitive to external stimuli. In this dynamical region, applying a stimulus to the network can reduce the trial-to-trial variability and shift the network oscillatory frequency while preserving the dynamical criticality. These multilevel features widely observed in different experiments cannot simultaneously occur in non-critical dynamical states. Furthermore, the dynamical mechanisms underlying these multilevel features are revealed using a semi-analytical mean-field theory that derives the macroscopic network field equations from the microscopic neuronal networks, enabling the analysis by nonlinear dynamics theory and linear noise approximation. The generic dynamical principle revealed here contributes to a more integrative understanding of neural systems and brain functions and incorporates multimodal and multilevel experimental observations. The E–I balanced neural network in combination with the effective mean-field theory can serve as a mechanistic modeling framework to study the multilevel neural dynamics underlying neural information and cognitive processes.

The complexity and variability of brain dynamical activity range from neuronal spiking and neural avalanches to oscillatory local field potentials of local neural circuits in both spontaneous and stimulus-evoked states. Such multilevel variable brain dynamics are functionally and behaviorally relevant and are principal components of the underlying circuit organization. To more comprehensively clarify their neural mechanisms, we use a bottom-up approach to study the stimulus–response dynamics of neural circuits. Our model assumes the following key biologically plausible components: excitation–inhibition (E–I) neuronal interaction and chemical synaptic coupling. We show that the circuits with E–I balance have a special dynamic sub-region, the critical region. Circuits around this region could account for the emergence of multilevel brain response patterns, both ongoing and stimulus-induced, observed in different experiments, including the reduction of trial-to-trial variability, effective modulation of gamma frequency, and preservation of criticality in the presence of a stimulus. We further analyze the corresponding nonlinear dynamical principles using a novel and highly generalizable semi-analytical mean-field theory. Our computational and theoretical studies explain the cross-level brain dynamical organization of spontaneous and evoked states in a more integrative manner.

Neuronal Avalanches Across the Rat Somatosensory Barrel Cortex and the Effect of Single Whisker Stimulation

Since its first experimental signatures, the so called "critical brain hypothesis" has been extensively studied. Yet, its actual foundations remain elusive. According to a widely accepted teleological reasoning, the brain would be poised to a critical state to optimize the mapping of the noisy and ever changing real-world inputs, thus suggesting that primary sensory cortical areas should be critical. We investigated whether a single barrel column of the somatosensory cortex of the anesthetized rat displays a critical behavior. Neuronal avalanches were recorded across all cortical layers in terms of both multi-unit activities and population local field potentials, and their behavior during spontaneous activity compared to the one evoked by a controlled single whisker deflection. By applying a maximum likelihood statistical method based on timeseries undersampling to fit the avalanches distributions, we show that neuronal avalanches are power law distributed for both multi-unit activities and local field potentials during spontaneous activity, with exponents that are spread along a scaling line. Instead, after the tactile stimulus, activity switches to a transient across-layers synchronization mode that appears to dominate the cortical representation of the single sensory input.

Less is more: Wiring-economical modular networks support self-sustained firing-economical neural avalanches for efficient processing

The brain network is notably cost-efficient, while the fundamental physical and dynamic mechanisms underlying its economical optimization in network structure and activity have not been determined. In this study, we investigate the intricate cost-efficient interplay between structure and dynamics in biologically plausible spatial modular neuronal network models. We observe that critical avalanche states from excitation-inhibition balance under modular network topology with less wiring cost can also achieve lower costs in firing but with strongly enhanced response sensitivity to stimuli. We derive mean-field equations that govern the macroscopic network dynamics through a novel approximate theory. The mechanism of low firing cost and stronger response in the form of critical avalanches is explained as a proximity to a Hopf bifurcation of the modules when increasing their connection density. Our work reveals the generic mechanism underlying the cost-efficient modular organization and critical dynamics widely observed in neural systems, providing insights into brain-inspired efficient computational designs.

Impact of GABAA and GABAB Inhibition on Cortical Dynamics and Perturbational Complexity during Synchronous and Desynchronized States

Quantitative estimations of spatiotemporal complexity of cortical activity patterns are used in the clinic as a measure of consciousness levels, but the cortical mechanisms involved are not fully understood. We used a version of the perturbational complexity index (PCI) adapted to multisite recordings from the ferret (either sex) cerebral cortex in vitro (sPCI) to investigate the role of GABAergic inhibition in cortical complexity. We studied two dynamical states: slow-wave activity (synchronous state) and desynchronized activity, that express low and high causal complexity respectively. Progressive blockade of GABAergic inhibition during both regimes revealed its impact on the emergent cortical activity and on sPCI. Gradual GABAA receptor blockade resulted in higher synchronization, being able to drive the network from a desynchronized to a synchronous state, with a progressive decrease of complexity (sPCI). Blocking GABAB receptors also resulted in a reduced sPCI, in particular when in a synchronous, slow wave state. Our findings demonstrate that physiological levels of inhibition contribute to the generation of dynamical richness and spatiotemporal complexity. However, if inhibition is diminished or enhanced, cortical complexity decreases. Using a computational model, we explored a larger parameter space in this relationship and demonstrate a link between excitatory/inhibitory balance and the complexity expressed by the cortical network.SIGNIFICANCE STATEMENT The spatiotemporal complexity of the activity expressed by the cerebral cortex is a highly revealing feature of the underlying network's state. Complexity varies with physiological brain states: it is higher during awake than during sleep states. But it also informs about pathologic states: in disorders of consciousness, complexity is lower in an unresponsive wakefulness syndrome than in a minimally conscious state. What are the network parameters that modulate complexity? Here we investigate how inhibition, mediated by either GABAA or GABAA receptors, influences cortical complexity. And we do this departing from two extreme functional states: a highly synchronous, slow-wave state, and a desynchronized one that mimics wakefulness. We find that there is an optimal level of inhibition in which complexity is highest.

Role of anaxonic local neurons in the crossover to continuously varying exponents for avalanche activity

Local anaxonic neurons with graded potential release are important ingredients of nervous systems, present in the olfactory bulb system of mammalians and in the human visual system, as well as in arthropods and nematodes. We develop a neuronal network model including both axonic and anaxonic neurons and monitor the activity tuned by the following parameters: the decay length of the graded potential in local neurons, the fraction of local neurons, the largest eigenvalue of the adjacency matrix, and the range of connections of the local neurons. Tuning the fraction of local neurons, we derive the phase diagram including two transition lines: a critical line separating subcritical and supercritical regions, characterized by power-law distributions of avalanche sizes and durations, and a bifurcation line. We find that the overall behavior of the system is controlled by a parameter tuning the relevance of local neuron transmission with respect to the axonal one. The statistical properties of spontaneous activity are affected by local neurons at large fractions and on the condition that the graded potential transmission dominates the axonal one. In this case the scaling properties of spontaneous activity exhibit continuously varying exponents, rather than the mean-field branching model universality class.

Percolation theory of self-exciting temporal processes

We investigate how the properties of inhomogeneous patterns of activity, appearing in many natural and social phenomena, depend on the temporal resolution used to define individual bursts of activity. To this end, we consider time series of microscopic events produced by a self-exciting Hawkes process, and leverage a percolation framework to study the formation of macroscopic bursts of activity as a function of the resolution parameter. We find that the very same process may result in different distributions of avalanche size and duration, which are understood in terms of the competition between the 1D percolation and the branching process universality class. Pure regimes for the individual classes are observed at specific values of the resolution parameter corresponding to the critical points of the percolation diagram. A regime of crossover characterized by a mixture of the two universal behaviors is observed in a wide region of the diagram. The hybrid scaling appears to be a likely outcome for an analysis of the time series based on a reasonably chosen, but not precisely adjusted, value of the resolution parameter.

Subsampled Directed-Percolation Models Explain Scaling Relations Experimentally Observed in the Brain

Recent experimental results on spike avalanches measured in the urethane-anesthetized rat cortex have revealed scaling relations that indicate a phase transition at a specific level of cortical firing rate variability. The scaling relations point to critical exponents whose values differ from those of a branching process, which has been the canonical model employed to understand brain criticality. This suggested that a different model, with a different phase transition, might be required to explain the data. Here we show that this is not necessarily the case. By employing two different models belonging to the same universality class as the branching process (mean-field directed percolation) and treating the simulation data exactly like experimental data, we reproduce most of the experimental results. We find that subsampling the model and adjusting the time bin used to define avalanches (as done with experimental data) are sufficient ingredients to change the apparent exponents of the critical point. Moreover, experimental data is only reproduced within a very narrow range in parameter space around the phase transition.

Hopf Bifurcation in Mean Field Explains Critical Avalanches in Excitation-Inhibition Balanced Neuronal Networks: A Mechanism for Multiscale Variability

Cortical neural circuits display highly irregular spiking in individual neurons but variably sized collective firing, oscillations and critical avalanches at the population level, all of which have functional importance for information processing. Theoretically, the balance of excitation and inhibition inputs is thought to account for spiking irregularity and critical avalanches may originate from an underlying phase transition. However, the theoretical reconciliation of these multilevel dynamic aspects in neural circuits remains an open question. Herein, we study excitation-inhibition (E-I) balanced neuronal network with biologically realistic synaptic kinetics. It can maintain irregular spiking dynamics with different levels of synchrony and critical avalanches emerge near the synchronous transition point. We propose a novel semi-analytical mean-field theory to derive the field equations governing the network macroscopic dynamics. It reveals that the E-I balanced state of the network manifesting irregular individual spiking is characterized by a macroscopic stable state, which can be either a fixed point or a periodic motion and the transition is predicted by a Hopf bifurcation in the macroscopic field. Furthermore, by analyzing public data, we find the coexistence of irregular spiking and critical avalanches in the spontaneous spiking activities of mouse cortical slice in vitro, indicating the universality of the observed phenomena. Our theory unveils the mechanism that permits complex neural activities in different spatiotemporal scales to coexist and elucidates a possible origin of the criticality of neural systems. It also provides a novel tool for analyzing the macroscopic dynamics of E-I balanced networks and its relationship to the microscopic counterparts, which can be useful for large-scale modeling and computation of cortical dynamics.

Spike-phase coupling as an order parameter in a leaky integrate-and-fire model

It is known that the leaky integrate-and-fire neural model shows a transition from irregular to synchronous firing by increasing the coupling between the neurons. However, a quantitative characterization of this order-disorder transition, that is, the determination of the order of transition and also the critical exponents in the case of continuous transition, is not entirely known. In this work, we consider a network of N excitatory neurons with local connections, residing on a square lattice with periodic boundary conditions. The cooperation between neurons K plays the role of the control parameter that generates criticality when the critical cooperation strength K_{c} is adopted. We introduce the population-averaged voltage (PAV) as a representative value of the network's cooperative activity. Then, we show that the coupling between the timing of spikes and the phase of temporal fluctuations of PAV defined as m resorts to identify a Kuramoto order parameter. By increasing K, we find a continuous transition from irregular spiking to a phase-locked state at the critical point K_{c}. We deploy the finite-size scaling analysis to calculate the critical exponents of this transition. To explore the formal indicator of criticality, we study the neuronal avalanches profile at this critical point and find a scaling behavior with the exponents in a fair agreement with the experimental values both in vivo and in vitro.

Measurement of excitation-inhibition ratio in autism spectrum disorder using critical brain dynamics

Balance between excitation (E) and inhibition (I) is a key principle for neuronal network organization and information processing. Consistent with this notion, excitation-inhibition imbalances are considered a pathophysiological mechanism in many brain disorders including autism spectrum disorder (ASD). However, methods to measure E/I ratios in human brain networks are lacking. Here, we present a method to quantify a functional E/I ratio (fE/I) from neuronal oscillations, and validate it in healthy subjects and children with ASD. We define structural E/I ratio in an in silico neuronal network, investigate how it relates to power and long-range temporal correlations (LRTC) of the network's activity, and use these relationships to design the fE/I algorithm. Application of this algorithm to the EEGs of healthy adults showed that fE/I is balanced at the population level and is decreased through GABAergic enforcement. In children with ASD, we observed larger fE/I variability and stronger LRTC compared to typically developing children (TDC). Interestingly, visual grading for EEG abnormalities that are thought to reflect E/I imbalances revealed elevated fE/I and LRTC in ASD children with normal EEG compared to TDC or ASD with abnormal EEG. We speculate that our approach will help understand physiological heterogeneity also in other brain disorders.

Pre-stimulus phase and amplitude regulation of phase-locked responses are maximized in the critical state

Understanding why identical stimuli give differing neuronal responses and percepts is a central challenge in research on attention and consciousness. Ongoing oscillations reflect functional states that bias processing of incoming signals through amplitude and phase. It is not known, however, whether the effect of phase or amplitude on stimulus processing depends on the long-term global dynamics of the networks generating the oscillations. Here, we show, using a computational model, that the ability of networks to regulate stimulus response based on pre-stimulus activity requires near-critical dynamics-a dynamical state that emerges from networks with balanced excitation and inhibition, and that is characterized by scale-free fluctuations. We also find that networks exhibiting critical oscillations produce differing responses to the largest range of stimulus intensities. Thus, the brain may bring its dynamics close to the critical state whenever such network versatility is required.

Spike-Timing-Dependent Plasticity With Axonal Delay Tunes Networks of Izhikevich Neurons to the Edge of Synchronization Transition With Scale-Free Avalanches

Critical brain hypothesis has been intensively studied both in experimental and theoretical neuroscience over the past two decades. However, some important questions still remain: (i) What is the critical point the brain operates at? (ii) What is the regulatory mechanism that brings about and maintains such a critical state? (iii) The critical state is characterized by scale-invariant behavior which is seemingly at odds with definitive brain oscillations? In this work we consider a biologically motivated model of Izhikevich neuronal network with chemical synapses interacting via spike-timing-dependent plasticity (STDP) as well as axonal time delay. Under generic and physiologically relevant conditions we show that the system is organized and maintained around a synchronization transition point as opposed to an activity transition point associated with an absorbing state phase transition. However, such a state exhibits experimentally relevant signs of critical dynamics including scale-free avalanches with finite-size scaling as well as critical branching ratios. While the system displays stochastic oscillations with highly correlated fluctuations, it also displays dominant frequency modes seen as sharp peaks in the power spectrum. The role of STDP as well as time delay is crucial in achieving and maintaining such critical dynamics, while the role of inhibition is not as crucial. In this way we provide possible answers to all three questions posed above. We also show that one can achieve supercritical or subcritical dynamics if one changes the average time delay associated with axonal conduction.

Oscillations and collective excitability in a model of stochastic neurons under excitatory and inhibitory coupling

We study a model with excitable neurons modeled as stochastic units with three states, representing quiescence, firing, and refractoriness. The transition rates between quiescence and firing depend exponentially on the number of firing neighbors, whereas all other rates are kept constant. This model class was shown to exhibit collective oscillations (synchronization) if neurons are spiking autonomously, but not if neurons are in the excitable regime. In both cases, neurons were restricted to interact through excitatory coupling. Here we show that a plethora of collective phenomena appear if inhibitory coupling is added. Besides the usual transition between an absorbing and an active phase, the model with excitatory and inhibitory neurons can also undergo reentrant transitions to an oscillatory phase. In the mean-field description, oscillations can emerge through supercritical or subcritical Hopf bifurcations, as well as through infinite period bifurcations. The model has bistability between active and oscillating behavior, as well as collective excitability, a regime where the system can display a peak of global activity when subject to a sufficiently strong perturbation. We employ a variant of the Shinomoto-Kuramoto order parameter to characterize the phase transitions and their system-size dependence.

The scale-invariant, temporal profile of neuronal avalanches in relation to cortical γ-oscillations

Activity cascades are found in many complex systems. In the cortex, they arise in the form of neuronal avalanches that capture ongoing and evoked neuronal activities at many spatial and temporal scales. The scale-invariant nature of avalanches suggests that the brain is in a critical state, yet predictions from critical theory on the temporal unfolding of avalanches have yet to be confirmed in vivo. Here we show in awake nonhuman primates that the temporal profile of avalanches follows a symmetrical, inverted parabola spanning up to hundreds of milliseconds. This parabola constrains how avalanches initiate locally, extend spatially and shrink as they evolve in time. Importantly, parabolas of different durations can be collapsed with a scaling exponent close to 2 supporting critical generational models of neuronal avalanches. Spontaneously emerging, transient γ-oscillations coexist with and modulate these avalanche parabolas thereby providing a temporal segmentation to inherently scale-invariant, critical dynamics. Our results identify avalanches and oscillations as dual principles in the temporal organization of brain activity.

Neuronal avalanches and time-frequency representations in stimulus-evoked activity

Neuronal avalanches are a hallmark feature of critical dynamics in the brain. While the theoretical framework of a critical branching processes is generally accepted for describing avalanches during ongoing brain activity, there is a current debate about the corresponding dynamical description during stimulus-evoked activity. As the brain activity evoked by external stimuli considerably varies in magnitude across time, it is not clear whether the parameters that govern the neuronal avalanche analysis (a threshold or a temporal scale) should be adaptively altered to accommodate these changes. Here, the relationship between neuronal avalanches and time-frequency representations of stimulus-evoked activity is explored. We show that neuronal avalanche metrics, calculated under a fixed threshold and temporal scale, reflect genuine changes in the underlying dynamics. In particular, event-related synchronization and de-synchronization are shown to align with variations in the power-law exponents of avalanche size distributions and the branching parameter (neural gain), as well as in the spatio-temporal spreading of avalanches. Nonetheless, the scale-invariant behavior associated with avalanches is shown to be a robust feature of healthy brain dynamics, preserved across various periods of stimulus-evoked activity and frequency bands. Taken together, the combined results suggest that throughout stimulus-evoked responses the operating point of the dynamics may drift within an extended-critical-like region.

Time-series thresholding and the definition of avalanche size

Avalanches whose sizes and durations are distributed as power laws appear in many contexts, from physics to geophysics and biology. Here we show that there is a hidden peril in thresholding continuous times series-from either empirical or synthetic data-for the identification of avalanches. In particular, we consider two possible alternative definitions of avalanche size used, e.g., in the empirical determination of avalanche exponents in the analysis of neural-activity data. By performing analytical and computational studies of an Ornstein-Uhlenbeck process (taken as a guiding example) we show that (1) if relatively large threshold values are employed to determine the beginning and ending of avalanches and (2) if-as sometimes done in the literature-avalanche sizes are defined as the total area (above zero) of the avalanche, then true asymptotic scaling behavior is not seen, instead the observations are dominated by transient effects. This problem-that we have detected in some recent works-leads to misinterpretations of the resulting scaling regimes.

Criticality between Cortical States

Since the first measurements of neuronal avalanches, the critical brain hypothesis has gained traction. However, if the brain is critical, what is the phase transition? For several decades, it has been known that the cerebral cortex operates in a diversity of regimes, ranging from highly synchronous states (with higher spiking variability) to desynchronized states (with lower spiking variability). Here, using both new and publicly available data, we test independent signatures of criticality and show that a phase transition occurs in an intermediate value of spiking variability, in both anesthetized and freely moving animals. The critical exponents point to a universality class different from mean-field directed percolation. Importantly, as the cortex hovers around this critical point, the avalanche exponents follow a linear relation that encompasses previous experimental results from different setups and is reproduced by a model.

Single-Cell Membrane Potential Fluctuations Evince Network Scale-Freeness and Quasicriticality

What information single neurons receive about general neural circuit activity is a fundamental question for neuroscience. Somatic membrane potential (V _m) fluctuations are driven by the convergence of synaptic inputs from a diverse cross-section of upstream neurons. Furthermore, neural activity is often scale-free, implying that some measurements should be the same, whether taken at large or small scales. Together, convergence and scale-freeness support the hypothesis that single V _m recordings carry useful information about high-dimensional cortical activity. Conveniently, the theory of "critical branching networks" (one purported explanation for scale-freeness) provides testable predictions about scale-free measurements that are readily applied to V _m fluctuations. To investigate, we obtained whole-cell current-clamp recordings of pyramidal neurons in visual cortex of turtles with unknown genders. We isolated fluctuations in V _m below the firing threshold and analyzed them by adapting the definition of "neuronal avalanches" (i.e., spurts of population spiking). The V _m fluctuations which we analyzed were scale-free and consistent with critical branching. These findings recapitulated results from large-scale cortical population data obtained separately in complementary experiments using microelectrode arrays described previously (Shew et al., 2015). Simultaneously recorded single-unit local field potential did not provide a good match, demonstrating the specific utility of V _m Modeling shows that estimation of dynamical network properties from neuronal inputs is most accurate when networks are structured as critical branching networks. In conclusion, these findings extend evidence of critical phenomena while also establishing subthreshold pyramidal neuron V _m fluctuations as an informative gauge of high-dimensional cortical population activity.SIGNIFICANCE STATEMENT The relationship between membrane potential (V _m) dynamics of single neurons and population dynamics is indispensable to understanding cortical circuits. Just as important to the biophysics of computation are emergent properties such as scale-freeness, where critical branching networks offer insight. This report makes progress on both fronts by comparing statistics from single-neuron whole-cell recordings with population statistics obtained with microelectrode arrays. Not only are fluctuations of somatic V _m scale-free, they match fluctuations of population activity. Thus, our results demonstrate appropriation of the brain's own subsampling method (convergence of synaptic inputs) while extending the range of fundamental evidence for critical phenomena in neural systems from the previously observed mesoscale (fMRI, LFP, population spiking) to the microscale, namely, V _m fluctuations.