Universal critical dynamics in high resolution neuronal avalanche data

An electrically resistive sheet of glial cells for amplifying signals of neuronal extracellular recordings

Electrical signals of neuronal cells can be recorded non-invasively and with a high degree of temporal resolution using multielectrode arrays (MEAs). However, signals that are recorded with these devices are small, usually 0.01%-0.1% of intracellular recordings. Here, we show that the amplitude of neuronal signals recorded with MEA devices can be amplified by covering neuronal networks with an electrically resistive sheet. The resistive sheet used in this study is a monolayer of glial cells, supportive cells in the brain. The glial cells were grown on a collagen-gel film that is permeable to oxygen and other nutrients. The impedance of the glial sheet was measured by electrochemical impedance spectroscopy, and equivalent circuit simulations were performed to theoretically investigate the effect of covering the neurons with such a resistive sheet. Finally, the effect of the resistive glial sheet was confirmed experimentally, showing a 6-fold increase in neuronal signals. This technique feasibly amplifies signals of MEA recordings.

How we move is universal: Scaling in the average shape of human activity

Human motor activity is constrained by the rhythmicity of the 24 hours circadian cycle, including the usual 12-15 hours sleep-wake cycle. However, activity fluctuations also appear over a wide range of temporal scales, from days to a few seconds, resulting from the concatenation of a myriad of individual smaller motor events. Furthermore,  individuals present different propensity to wakefulness and thus to motor activity throughout the circadian cycle.  Are activity fluctuations across temporal scales intrinsically different, or is there a universal description encompassing them? Is this description also universal across individuals, considering the aforementioned variability?  Here we establish the presence of universality in motor activity fluctuations based on the empirical study of a month of continuous wristwatch accelerometer recordings. We study the scaling of average fluctuations across temporal scales and determine a universal law characterized by critical exponents $\alpha$, $\tau$ and $1/{ \mu}$. Results are highly reminiscent of the universality described for the average shape of avalanches in systems exhibiting crackling noise. Beyond its theoretical relevance, the present results can be important for developing objective markers of healthy  as well as pathological human motor behavior.Received: 18 October 2015,  Accepted: 10 November 2015; Edited by: E. Mizraji; Reviewed by: J. Lin, Department of Physics, Washington College, Maryland, USA.; DOI: http://dx.doi.org/10.4279/PIP.070017Cite as: D R Chialvo, A M Gonzalez Torrado, E Gudowska-Nowak, J K Ochab, P Montoya, M A Nowak, E Tagliazucchi, Papers in Physics 7, 070017 (2015)This paper, by D R Chialvo, A M Gonzalez Torrado, E Gudowska-Nowak, J K Ochab, P Montoya, M A Nowak, E Tagliazucchi, is licensed under the Creative Commons Attribution License 3.0.

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.

Functional Clusters, Hubs, and Communities in the Cortical Microconnectome

Although relationships between networks of different scales have been observed in macroscopic brain studies, relationships between structures of different scales in networks of neurons are unknown. To address this, we recorded from up to 500 neurons simultaneously from slice cultures of rodent somatosensory cortex. We then measured directed effective networks with transfer entropy, previously validated in simulated cortical networks. These effective networks enabled us to evaluate distinctive nonrandom structures of connectivity at 2 different scales. We have 4 main findings. First, at the scale of 3-6 neurons (clusters), we found that high numbers of connections occurred significantly more often than expected by chance. Second, the distribution of the number of connections per neuron (degree distribution) had a long tail, indicating that the network contained distinctively high-degree neurons, or hubs. Third, at the scale of tens to hundreds of neurons, we typically found 2-3 significantly large communities. Finally, we demonstrated that communities were relatively more robust than clusters against shuffling of connections. We conclude the microconnectome of the cortex has specific organization at different scales, as revealed by differences in robustness. We suggest that this information will help us to understand how the microconnectome is robust against damage.

Thermodynamics and signatures of criticality in a network of neurons

The activity of a neural network is defined by patterns of spiking and silence from the individual neurons. Because spikes are (relatively) sparse, patterns of activity with increasing numbers of spikes are less probable, but, with more spikes, the number of possible patterns increases. This tradeoff between probability and numerosity is mathematically equivalent to the relationship between entropy and energy in statistical physics. We construct this relationship for populations of up to N = 160 neurons in a small patch of the vertebrate retina, using a combination of direct and model-based analyses of experiments on the response of this network to naturalistic movies. We see signs of a thermodynamic limit, where the entropy per neuron approaches a smooth function of the energy per neuron as N increases. The form of this function corresponds to the distribution of activity being poised near an unusual kind of critical point. We suggest further tests of criticality, and give a brief discussion of its functional significance.

Brain Performance versus Phase Transitions

We here illustrate how a well-founded study of the brain may originate in assuming analogies with phase-transition phenomena. Analyzing to what extent a weak signal endures in noisy environments, we identify the underlying mechanisms, and it results a description of how the excitability associated to (non-equilibrium) phase changes and criticality optimizes the processing of the signal. Our setting is a network of integrate-and-fire nodes in which connections are heterogeneous with rapid time-varying intensities mimicking fatigue and potentiation. Emergence then becomes quite robust against wiring topology modification—in fact, we considered from a fully connected network to the Homo sapiens connectome—showing the essential role of synaptic flickering on computations. We also suggest how to experimentally disclose significant changes during actual brain operation.

Irregular spiking of pyramidal neurons organizes as scale-invariant neuronal avalanches in the awake state

Spontaneous fluctuations in neuronal activity emerge at many spatial and temporal scales in cortex. Population measures found these fluctuations to organize as scale-invariant neuronal avalanches, suggesting cortical dynamics to be critical. Macroscopic dynamics, though, depend on physiological states and are ambiguous as to their cellular composition, spatiotemporal origin, and contributions from synaptic input or action potential (AP) output. Here, we study spontaneous firing in pyramidal neurons (PNs) from rat superficial cortical layers in vivo and in vitro using 2-photon imaging. As the animal transitions from the anesthetized to awake state, spontaneous single neuron firing increases in irregularity and assembles into scale-invariant avalanches at the group level. In vitro spike avalanches emerged naturally yet required balanced excitation and inhibition. This demonstrates that neuronal avalanches are linked to the global physiological state of wakefulness and that cortical resting activity organizes as avalanches from firing of local PN groups to global population activity.

DOI:http://dx.doi.org/10.7554/eLife.07224.001

Even when we are not engaged in any specific task, the brain shows coordinated patterns of spontaneous activity that can be monitored using electrodes placed on the scalp. This resting activity shapes the way that the brain responds to subsequent stimuli. Changes in resting activity patterns are seen in various neurological and psychiatric disorders, as well as in healthy individuals following sleep deprivation.

The brain's outer layer is known as the cortex. On a large scale, when monitoring many thousands of neurons, resting activity in the cortex demonstrates propagation in the brain in an organized manner. Specifically, resting activity was found to organize as so-called neuronal avalanches, in which large bursts of neuronal activity are grouped with medium-sized and smaller bursts in a very characteristic order. In fact, the sizes of these bursts—that is, the number of neurons that fire—are found to be scale-invariant, that is, the ratio of large bursts to medium-sized bursts is the same as that of medium-sized to small bursts. Such scale-invariance suggests that neuronal bursts are not independent of one another. However, it is largely unclear how neuronal avalanches arise from individual neurons, which fire simply in a noisy, irregular manner.

Bellay, Klaus et al. have now provided insights into this process by examining patterns of firing of a particular type of neuron—known as a pyramidal cell—in the cortex of rats as they recover from anesthesia. As the animals awaken, the firing of individual pyramidal cells in the cortex becomes even more irregular than under anesthesia. However, by considering the activity of a group of these neurons, Bellay, Klaus et al. realized that it is this more irregular firing that gives rise to neuronal avalanches, and that this occurs only when the animals are awake. Further experiments on individual pyramidal cells grown in the laboratory confirmed that neuronal avalanches emerge spontaneously from the irregular firing of individual neurons. These avalanches depend on there being a balance between two types of activity among the cells: ‘excitatory’ activity that causes other neurons to fire, and ‘inhibitory’ activity that prevents neuronal firing.

Given that resting activity influences the brain's responses to the outside world, the origins of neuronal avalanches are likely to provide clues about the way the brain processes information. Future experiments should also examine the possibility that the emergence of neuronal avalanches marks the transition from unconsciousness to wakefulness within the brain.

DOI:http://dx.doi.org/10.7554/eLife.07224.002

Cascades and cognitive state: focused attention incurs subcritical dynamics

The analysis of neuronal avalanches supports the hypothesis that the human cortex operates with critical neural dynamics. Here, we investigate the relationship between cascades of activity in electroencephalogram data, cognitive state, and reaction time in humans using a multimodal approach. We recruited 18 healthy volunteers for the acquisition of simultaneous electroencephalogram and functional magnetic resonance imaging during both rest and during a visuomotor cognitive task. We compared distributions of electroencephalogram-derived cascades to reference power laws for task and rest conditions. We then explored the large-scale spatial correspondence of these cascades in the simultaneously acquired functional magnetic resonance imaging data. Furthermore, we investigated whether individual variability in reaction times is associated with the amount of deviation from power law form. We found that while resting state cascades are associated with approximate power law form, the task state is associated with subcritical dynamics. Furthermore, we found that electroencephalogram cascades are related to blood oxygen level-dependent activation, predominantly in sensorimotor brain regions. Finally, we found that decreased reaction times during the task condition are associated with increased proximity to power law form of cascade distributions. These findings suggest that the resting state is associated with near-critical dynamics, in which a high dynamic range and a large repertoire of brain states may be advantageous. In contrast, a focused cognitive task induces subcritical dynamics, which is associated with a lower dynamic range, which in turn may reduce elements of interference affecting task performance.

Self-organized criticality in cortical assemblies occurs in concurrent scale-free and small-world networks

The spontaneous activity of cortical networks is characterized by the emergence of different dynamic states. Although several attempts were accomplished to understand the origin of these dynamics, the underlying factors continue to be elusive. In this work, we specifically investigated the interplay between network topology and spontaneous dynamics within the framework of self-organized criticality (SOC). The obtained results support the hypothesis that the emergence of critical states occurs in specific complex network topologies. By combining multi-electrode recordings of spontaneous activity of in vitro cortical assemblies with theoretical models, we demonstrate that different 'connectivity rules' drive the network towards different dynamic states. In particular, scale-free architectures with different degree of small-worldness account better for the variability observed in experimental data, giving rise to different dynamic states. Moreover, in relationship with the balance between excitation and inhibition and percentage of inhibitory hubs, the simulated cortical networks fall in a critical regime.

Cortical burst dynamics predict clinical outcome early in extremely preterm infants

Intermittent bursts of electrical activity are a ubiquitous signature of very early brain activity. Previous studies have largely focused on assessing the amplitudes of these transient cortical bursts or the intervals between them. Recent advances in basic neuroscience have identified the presence of scale-free 'avalanche' processes in bursting patterns of cortical activity in other clinical contexts. Here, we hypothesize that cortical bursts in human preterm infants also exhibit scale-free properties, providing new insights into the nature, temporal evolution, and prognostic value of spontaneous brain activity in the days immediately following preterm birth. We examined electroencephalographic recordings from 43 extremely preterm infants (gestational age 22-28 weeks) and demonstrated that their cortical bursts exhibit scale-free properties as early as 12 h after birth. The scaling relationships of cortical bursts correlate significantly with later mental development-particularly within the first 12 h of life. These findings show that early preterm brain activity is characterized by scale-free dynamics which carry developmental significance, hence offering novel means for rapid and early clinical prediction of neurodevelopmental outcomes.

How can we study reasoning in the brain?

The brain did not develop a dedicated device for reasoning. This fact bears dramatic consequences. While for perceptuo-motor functions neural activity is shaped by the input's statistical properties, and processing is carried out at high speed in hardwired spatially segregated modules, in reasoning, neural activity is driven by internal dynamics and processing times, stages, and functional brain geometry are largely unconstrained a priori. Here, it is shown that the complex properties of spontaneous activity, which can be ignored in a short-lived event-related world, become prominent at the long time scales of certain forms of reasoning. It is argued that the neural correlates of reasoning should in fact be defined in terms of non-trivial generic properties of spontaneous brain activity, and that this implies resorting to concepts, analytical tools, and ways of designing experiments that are as yet non-standard in cognitive neuroscience. The implications in terms of models of brain activity, shape of the neural correlates, methods of data analysis, observability of the phenomenon, and experimental designs are discussed.

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.

Voltage imaging of waking mouse cortex reveals emergence of critical neuronal dynamics

Complex cognitive processes require neuronal activity to be coordinated across multiple scales, ranging from local microcircuits to cortex-wide networks. However, multiscale cortical dynamics are not well understood because few experimental approaches have provided sufficient support for hypotheses involving multiscale interactions. To address these limitations, we used, in experiments involving mice, genetically encoded voltage indicator imaging, which measures cortex-wide electrical activity at high spatiotemporal resolution. Here we show that, as mice recovered from anesthesia, scale-invariant spatiotemporal patterns of neuronal activity gradually emerge. We show for the first time that this scale-invariant activity spans four orders of magnitude in awake mice. In contrast, we found that the cortical dynamics of anesthetized mice were not scale invariant. Our results bridge empirical evidence from disparate scales and support theoretical predictions that the awake cortex operates in a dynamical regime known as criticality. The criticality hypothesis predicts that small-scale cortical dynamics are governed by the same principles as those governing larger-scale dynamics. Importantly, these scale-invariant principles also optimize certain aspects of information processing. Our results suggest that during the emergence from anesthesia, criticality arises as information processing demands increase. We expect that, as measurement tools advance toward larger scales and greater resolution, the multiscale framework offered by criticality will continue to provide quantitative predictions and insight on how neurons, microcircuits, and large-scale networks are dynamically coordinated in the brain.

Dynamical criticality in the collective activity of a population of retinal neurons

Recent experimental results based on multielectrode and imaging techniques have reinvigorated the idea that large neural networks operate near a critical point, between order and disorder. However, evidence for criticality has relied on the definition of arbitrary order parameters, or on models that do not address the dynamical nature of network activity. Here we introduce a novel approach to assess criticality that overcomes these limitations, while encompassing and generalizing previous criteria. We find a simple model to describe the global activity of large populations of ganglion cells in the rat retina, and show that their statistics are poised near a critical point. Taking into account the temporal dynamics of the activity greatly enhances the evidence for criticality, revealing it where previous methods would not. The approach is general and could be used in other biological networks.

Critical and maximally informative encoding between neural populations in the retina

Computation in the brain involves multiple types of neurons, yet the organizing principles for how these neurons work together remain unclear. Information theory has offered explanations for how different types of neurons can maximize the transmitted information by encoding different stimulus features. However, recent experiments indicate that separate neuronal types exist that encode the same filtered version of the stimulus, but then the different cell types signal the presence of that stimulus feature with different thresholds. Here we show that the emergence of these neuronal types can be quantitatively described by the theory of transitions between different phases of matter. The two key parameters that control the separation of neurons into subclasses are the mean and standard deviation (SD) of noise affecting neural responses. The average noise across the neural population plays the role of temperature in the classic theory of phase transitions, whereas the SD is equivalent to pressure or magnetic field, in the case of liquid-gas and magnetic transitions, respectively. Our results account for properties of two recently discovered types of salamander Off retinal ganglion cells, as well as the absence of multiple types of On cells. We further show that, across visual stimulus contrasts, retinal circuits continued to operate near the critical point whose quantitative characteristics matched those expected near a liquid-gas critical point and described by the nearest-neighbor Ising model in three dimensions. By operating near a critical point, neural circuits can maximize information transmission in a given environment while retaining the ability to quickly adapt to a new environment.

Unambiguous reconstruction of network structure using avalanche dynamics

A robust method for inferring the structure of networks is presented based on the one-to-one correspondence between the expected composition of cascades of bursts of activity, called crackling noise or avalanches, and the weight matrix. Using a model of neuronal avalanches as a paradigmatic example, we derive this correspondence exactly by calculating the closed-form expression of the joint probability distribution of avalanche sizes obtained by counting separately the number of elements active in each subnetwork during avalanches.

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.

On the temporal organization of neuronal avalanches

Spontaneous activity of cortex in vitro and in vivo has been shown to organize as neuronal avalanches. Avalanches are cascades of neuronal activity that exhibit a power law in their size and duration distribution, typical features of balanced systems in a critical state. Recently it has been shown that the distribution of quiet times between consecutive avalanches in rat cortex slice cultures displays a non-monotonic behavior with a power law decay at short time scales. This behavior has been attributed to the slow alternation between up and down-states. Here we further characterize the avalanche process and investigate how the functional behavior of the quiet time distribution depends on the fine structure of avalanche sequences. By systematically removing smaller avalanches from the experimental time series we show that size and quiet times are correlated and highlight that avalanche occurrence exhibits the characteristic periodicity of θ and β/γ oscillations, which jointly emerge in most of the analyzed samples. Furthermore, our analysis indicates that smaller avalanches tend to be associated with faster β/γ oscillations, whereas larger ones are associated with slower θ and 1-2 Hz oscillations. In particular, large avalanches corresponding to θ cycles trigger cascades of smaller ones, which occur at β/γ frequency. This temporal structure follows closely the one of nested θ - β/γ oscillations. Finally we demonstrate that, because of the multiple time scales characterizing avalanche dynamics, the distributions of quiet times between avalanches larger than a certain size do not collapse onto a unique function when rescaled by the average occurrence rate. However, when considered separately in the up-state and in the down-state, these distributions are solely controlled by the respective average rate and two different unique function can be identified.

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.

Slow synaptic dynamics in a network: from exponential to power-law forgetting

We investigate a mean-field model of interacting synapses on a directed neural network. Our interest lies in the slow adaptive dynamics of synapses, which are driven by the fast dynamics of the neurons they connect. Cooperation is modeled from the usual Hebbian perspective, while competition is modeled by an original polarity-driven rule. The emergence of a critical manifold culminating in a tricritical point is crucially dependent on the presence of synaptic competition. This leads to a universal 1/t power-law relaxation of the mean synaptic strength along the critical manifold and an equally universal 1/√[t] relaxation at the tricritical point, to be contrasted with the exponential relaxation that is otherwise generic. In turn, this leads to the natural emergence of long- and short-term memory from different parts of parameter space in a synaptic network, which is the most original and important result of our present investigations.

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.

Critical role for resource constraints in neural models

Criticality has emerged as a leading dynamical candidate for healthy and pathological neuronal activity. At the heart of criticality in neural systems is the need for parameters to be tuned to specific values or for the existence of self-organizing mechanisms. Existing models lack precise physiological descriptions for how the brain maintains its tuning near a critical point. In this paper we argue that a key ingredient missing from the field is a formulation of reciprocal coupling between neural activity and metabolic resources. We propose that the constraint of optimizing the balance between energy use and activity plays a major role in tuning brain states to lie near criticality. Important recent findings aligned with our viewpoint have emerged from analyses of disorders that involve severe metabolic disturbances and alter scale-free properties of brain dynamics, including burst suppression. Moreover, we argue that average shapes of neuronal avalanches are a signature of scale-free activity that offers sharper insights into underlying mechanisms than afforded by traditional analyses of avalanche statistics.

Large-scale, high-resolution multielectrode-array recording depicts functional network differences of cortical and hippocampal cultures

Understanding the detailed circuitry of functioning neuronal networks is one of the major goals of neuroscience. Recent improvements in neuronal recording techniques have made it possible to record the spiking activity from hundreds of neurons simultaneously with sub-millisecond temporal resolution. Here we used a 512-channel multielectrode array system to record the activity from hundreds of neurons in organotypic cultures of cortico-hippocampal brain slices from mice. To probe the network structure, we employed a wavelet transform of the cross-correlogram to categorize the functional connectivity in different frequency ranges. With this method we directly compare, for the first time, in any preparation, the neuronal network structures of cortex and hippocampus, on the scale of hundreds of neurons, with sub-millisecond time resolution. Among the three frequency ranges that we investigated, the lower two frequency ranges (gamma (30–80 Hz) and beta (12–30 Hz) range) showed similar network structure between cortex and hippocampus, but there were many significant differences between these structures in the high frequency range (100–1000 Hz). The high frequency networks in cortex showed short tailed degree-distributions, shorter decay length of connectivity density, smaller clustering coefficients, and positive assortativity. Our results suggest that our method can characterize frequency dependent differences of network architecture from different brain regions. Crucially, because these differences between brain regions require millisecond temporal scales to be observed and characterized, these results underscore the importance of high temporal resolution recordings for the understanding of functional networks in neuronal systems.

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.

Scale-free bursting in human cortex following hypoxia at birth

The human brain is fragile in the face of oxygen deprivation. Even a brief interruption of metabolic supply at birth challenges an otherwise healthy neonatal cortex, leading to a cascade of homeostatic responses. During recovery from hypoxia, cortical activity exhibits a period of highly irregular electrical fluctuations known as burst suppression. Here we show that these bursts have fractal properties, with power-law scaling of burst sizes across a remarkable 5 orders of magnitude and a scale-free relationship between burst sizes and durations. Although burst waveforms vary greatly, their average shape converges to a simple form that is asymmetric at long time scales. Using a simple computational model, we argue that this asymmetry reflects activity-dependent changes in the excitatory-inhibitory balance of cortical neurons. Bursts become more symmetric following the resumption of normal activity, with a corresponding reorganization of burst scaling relationships. These findings place burst suppression in the broad class of scale-free physical processes termed crackling noise and suggest that the resumption of healthy activity reflects a fundamental reorganization in the relationship between neuronal activity and its underlying metabolic constraints.

Functional significance of complex fluctuations in brain activity: from resting state to cognitive neuroscience

Behavioral studies have shown that human cognition is characterized by properties such as temporal scale invariance, heavy-tailed non-Gaussian distributions, and long-range correlations at long time scales, suggesting models of how (non observable) components of cognition interact. On the other hand, results from functional neuroimaging studies show that complex scaling and intermittency may be generic spatio-temporal properties of the brain at rest. Somehow surprisingly, though, hardly ever have the neural correlates of cognition been studied at time scales comparable to those at which cognition shows scaling properties. Here, we analyze the meanings of scaling properties and the significance of their task-related modulations for cognitive neuroscience. It is proposed that cognitive processes can be framed in terms of complex generic properties of brain activity at rest and, ultimately, of functional equations, limiting distributions, symmetries, and possibly universality classes characterizing them.

Mean-field behavior as a result of noisy local dynamics in self-organized criticality: neuroscience implications

Motivated by recent experiments in neuroscience which indicate that neuronal avalanches exhibit scale invariant behavior similar to self-organized critical systems, we study the role of noisy (nonconservative) local dynamics on the critical behavior of a sandpile model which can be taken to mimic the dynamics of neuronal avalanches. We find that despite the fact that noise breaks the strict local conservation required to attain criticality, our system exhibits true criticality for a wide range of noise in various dimensions, given that conservation is respected on the average. Although the system remains critical, exhibiting finite-size scaling, the value of critical exponents change depending on the intensity of local noise. Interestingly, for a sufficiently strong noise level, the critical exponents approach and saturate at their mean-field values, consistent with empirical measurements of neuronal avalanches. This is confirmed for both two and three dimensional models. However, the addition of noise does not affect the exponents at the upper critical dimension (D = 4). In addition to an extensive finite-size scaling analysis of our systems, we also employ a useful time-series analysis method to establish true criticality of noisy systems. Finally, we discuss the implications of our work in neuroscience as well as some implications for the general phenomena of criticality in nonequilibrium systems.

Avalanches in self-organized critical neural networks: a minimal model for the neural SOC universality class

The brain keeps its overall dynamics in a corridor of intermediate activity and it has been a long standing question what possible mechanism could achieve this task. Mechanisms from the field of statistical physics have long been suggesting that this homeostasis of brain activity could occur even without a central regulator, via self-organization on the level of neurons and their interactions, alone. Such physical mechanisms from the class of self-organized criticality exhibit characteristic dynamical signatures, similar to seismic activity related to earthquakes. Measurements of cortex rest activity showed first signs of dynamical signatures potentially pointing to self-organized critical dynamics in the brain. Indeed, recent more accurate measurements allowed for a detailed comparison with scaling theory of non-equilibrium critical phenomena, proving the existence of criticality in cortex dynamics. We here compare this new evaluation of cortex activity data to the predictions of the earliest physics spin model of self-organized critical neural networks. We find that the model matches with the recent experimental data and its interpretation in terms of dynamical signatures for criticality in the brain. The combination of signatures for criticality, power law distributions of avalanche sizes and durations, as well as a specific scaling relationship between anomalous exponents, defines a universality class characteristic of the particular critical phenomenon observed in the neural experiments. Thus the model is a candidate for a minimal model of a self-organized critical adaptive network for the universality class of neural criticality. As a prototype model, it provides the background for models that may include more biological details, yet share the same universality class characteristic of the homeostasis of activity in the brain.

Novel features of early burst suppression predict outcome after birth asphyxia

Burst suppression patterns in the electroencephalogram are a reliable marker of recent severe brain insult. Here we analyze statistical properties of bursts occurring in 20 electroencephalographic recordings acquired from hypothermic asphyxic newborns in the hours immediately following birth. We show that the distributions of burst area and duration in these acute data predict later clinical outcome in both structural neuroimaging and neurodevelopment. Our findings indicate the first early electroencephalographic metrics that offer outcome prediction in asphyxic neonates undergoing hypothermia treatment.

Intrinsic noise induces critical behavior in leaky Markovian networks leading to avalanching

The role intrinsic statistical fluctuations play in creating avalanches – patterns of complex bursting activity with scale-free properties – is examined in leaky Markovian networks. Using this broad class of models, we develop a probabilistic approach that employs a potential energy landscape perspective coupled with a macroscopic description based on statistical thermodynamics. We identify six important thermodynamic quantities essential for characterizing system behavior as a function of network size: the internal potential energy, entropy, free potential energy, internal pressure, pressure, and bulk modulus. In agreement with classical phase transitions, these quantities evolve smoothly as a function of the network size until a critical value is reached. At that value, a discontinuity in pressure is observed that leads to a spike in the bulk modulus demarcating loss of thermodynamic robustness. We attribute this novel result to a reallocation of the ground states (global minima) of the system's stationary potential energy landscape caused by a noise-induced deformation of its topographic surface. Further analysis demonstrates that appreciable levels of intrinsic noise can cause avalanching, a complex mode of operation that dominates system dynamics at near-critical or subcritical network sizes. Illustrative examples are provided using an epidemiological model of bacterial infection, where avalanching has not been characterized before, and a previously studied model of computational neuroscience, where avalanching was erroneously attributed to specific neural architectures. The general methods developed here can be used to study the emergence of avalanching (and other complex phenomena) in many biological, physical and man-made interaction networks.

Networks of noisy interacting components arise in diverse scientific disciplines. Here, we develop a mathematical framework to study the underlying causes of a bursting phenomenon in network activity known as avalanching. As prototypical examples, we study a model of disease spreading in a population of individuals and a model of brain activity in a neural network. Although avalanching is well-documented in neural networks, thought to be crucial for learning, information processing, and memory, it has not been studied before in disease spreading. We employ tools originally used to analyze thermodynamic systems to argue that randomness in the actions of individual network components plays a fundamental role in avalanche formation. We show that avalanching is a spontaneous behavior, brought about by a phenomenon reminiscent to a phase transition in statistical mechanics, caused by increasing randomness as the network size decreases. Our work demonstrates that a previously suggested balanced feed-forward network structure is not necessary for neuronal avalanching. Instead, we attribute avalanching to a reallocation of the global minima of the network's stationary potential energy landscape, caused by a noise-induced deformation of its topographic surface.

Critical phenomena and noise-induced phase transitions in neuronal networks

We study numerically and analytically first- and second-order phase transitions in neuronal networks stimulated by shot noise (a flow of random spikes bombarding neurons). Using an exactly solvable cortical model of neuronal networks on classical random networks, we find critical phenomena accompanying the transitions and their dependence on the shot noise intensity. We show that a pattern of spontaneous neuronal activity near a critical point of a phase transition is a characteristic property that can be used to identify the bifurcation mechanism of the transition. We demonstrate that bursts and avalanches are precursors of a first-order phase transition, paroxysmal-like spikes of activity precede a second-order phase transition caused by a saddle-node bifurcation, while irregular spindle oscillations represent spontaneous activity near a second-order phase transition caused by a supercritical Hopf bifurcation. Our most interesting result is the observation of the paroxysmal-like spikes. We show that a paroxysmal-like spike is a single nonlinear event that appears instantly from a low background activity with a rapid onset, reaches a large amplitude, and ends up with an abrupt return to lower activity. These spikes are similar to single paroxysmal spikes and sharp waves observed in electroencephalographic (EEG) measurements. Our analysis shows that above the saddle-node bifurcation, sustained network oscillations appear with a large amplitude but a small frequency in contrast to network oscillations near the Hopf bifurcation that have a small amplitude but a large frequency. We discuss an amazing similarity between excitability of the cortical model stimulated by shot noise and excitability of the Morris-Lecar neuron stimulated by an applied current.

Fading signatures of critical brain dynamics during sustained wakefulness in humans

Sleep encompasses approximately a third of our lifetime, yet its purpose and biological function are not well understood. Without sleep optimal brain functioning such as responsiveness to stimuli, information processing, or learning may be impaired. Such observations suggest that sleep plays a crucial role in organizing or reorganizing neuronal networks of the brain toward states where information processing is optimized. Increasing evidence suggests that cortical neuronal networks operate near a critical state characterized by balanced activity patterns, which supports optimal information processing. However, it remains unknown whether critical dynamics is affected in the course of wake and sleep, which would also impact information processing. Here, we show that signatures of criticality are progressively disturbed during wake and restored by sleep. We demonstrate that the precise power-laws governing the cascading activity of neuronal avalanches and the distribution of phase-lock intervals in human electroencephalographic recordings are increasingly disarranged during sustained wakefulness. These changes are accompanied by a decrease in variability of synchronization. Interpreted in the context of a critical branching process, these seemingly different findings indicate a decline of balanced activity and progressive distance from criticality toward states characterized by an imbalance toward excitation where larger events prevail dynamics. Conversely, sleep restores the critical state resulting in recovered power-law characteristics in activity and variability of synchronization. These findings support the intriguing hypothesis that sleep may be important to reorganize cortical network dynamics to a critical state thereby assuring optimal computational capabilities for the following time awake.

Single-neuron criticality optimizes analog dendritic computation

Active dendritic branchlets enable the propagation of dendritic spikes, whose computational functions remain an open question. Here we propose a concrete function to the active channels in large dendritic trees. Modelling the input-output response of large active dendritic arbors subjected to complex spatio-temporal inputs and exhibiting non-stereotyped dendritic spikes, we find that the dendritic arbor can undergo a continuous phase transition from a quiescent to an active state, thereby exhibiting spontaneous and self-sustained localized activity as suggested by experiments. Analogously to the critical brain hypothesis, which states that neuronal networks self-organize near criticality to take advantage of its specific properties, here we propose that neurons with large dendritic arbors optimize their capacity to distinguish incoming stimuli at the critical state. We suggest that "computation at the edge of a phase transition" is more compatible with the view that dendritic arbors perform an analog rather than a digital dendritic computation.

Controlling avalanche criticality in 2D nano arrays

Many physical systems respond to slowly changing external force through avalanches spanning broad range of sizes. Some systems crackle even without apparent external force, such as bursts of neuronal activity or charge transfer avalanches in 2D molecular layers. Advanced development of theoretical models describing disorder-induced critical phenomena calls for experiments probing the dynamics upon tuneable disorder. Here we show that isomeric structural transitions in 2D organic self-assembled monolayer (SAM) exhibit critical dynamics with experimentally tuneable disorder. The system consists of field effect transistor coupled through SAM to illuminated semiconducting nanocrystals (NCs). Charges photoinduced in NCs are transferred through SAM to the transistor surface and modulate its conductivity. Avalanches of isomeric structural transitions are revealed by measuring the current noise I(t) of the transistor. Accumulated surface traps charges reduce dipole moments of the molecules, decrease their coupling, and thus decrease the critical disorder of the SAM enabling its tuning during experiments.

Universal organization of resting brain activity at the thermodynamic critical point

Thermodynamic criticality describes emergent phenomena in a wide variety of complex systems. In the mammalian cortex, one type of complex dynamics that spontaneously emerges from neuronal interactions has been characterized as neuronal avalanches. Several aspects of neuronal avalanches such as their size and life time distributions are described by power laws with unique exponents, indicating an underlying critical branching process that governs avalanche formation. Here, we show that neuronal avalanches also reflect an organization of brain dynamics close to a thermodynamic critical point. We recorded spontaneous cortical activity in monkeys and humans at rest using high-density intracranial microelectrode arrays and magnetoencephalography, respectively. By numerically changing a control parameter equivalent to thermodynamic temperature, we observed typical critical behavior in cortical activities near the actual physiological condition, including the phase transition of an order parameter, as well as the divergence of susceptibility and specific heat. Finite-size scaling of these quantities allowed us to derive robust critical exponents highly consistent across monkey and humans that uncover a distinct, yet universal organization of brain dynamics. Our results demonstrate that normal brain dynamics at rest resides near or at criticality, which maximizes several aspects of information processing such as input sensitivity and dynamic range.

A place for time: the spatiotemporal structure of neural dynamics during natural audition

We use functional magnetic resonance imaging (fMRI) to analyze neural responses to natural auditory stimuli. We characterize the fMRI time series through the shape of the voxel power spectrum and find that the timescales of neural dynamics vary along a spatial gradient, with faster dynamics in early auditory cortex and slower dynamics in higher order brain regions. The timescale gradient is observed through the unsupervised clustering of the power spectra of individual brains, both in the presence and absence of a stimulus, and is enhanced in the stimulus-locked component that is shared across listeners. Moreover, intrinsically faster dynamics occur in areas that respond preferentially to momentary stimulus features, while the intrinsically slower dynamics occur in areas that integrate stimulus information over longer timescales. These observations connect the timescales of intrinsic neural dynamics to the timescales of information processing, suggesting a temporal organizing principle for neural computation across the cerebral cortex.

Critical avalanches and subsampling in map-based neural networks coupled with noisy synapses

Many different kinds of noise are experimentally observed in the brain. Among them, we study a model of noisy chemical synapse and obtain critical avalanches for the spatiotemporal activity of the neural network. Neurons and synapses are modeled by dynamical maps. We discuss the relevant neuronal and synaptic properties to achieve the critical state. We verify that networks of functionally excitable neurons with fast synapses present power-law avalanches, due to rebound spiking dynamics. We also discuss the measuring of neuronal avalanches by subsampling our data, shedding light on the experimental search for self-organized criticality in neural networks.

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.

Neural avalanches at the critical point between replay and non-replay of spatiotemporal patterns

We model spontaneous cortical activity with a network of coupled spiking units, in which multiple spatio-temporal patterns are stored as dynamical attractors. We introduce an order parameter, which measures the overlap (similarity) between the activity of the network and the stored patterns. We find that, depending on the excitability of the network, different working regimes are possible. For high excitability, the dynamical attractors are stable, and a collective activity that replays one of the stored patterns emerges spontaneously, while for low excitability, no replay is induced. Between these two regimes, there is a critical region in which the dynamical attractors are unstable, and intermittent short replays are induced by noise. At the critical spiking threshold, the order parameter goes from zero to one, and its fluctuations are maximized, as expected for a phase transition (and as observed in recent experimental results in the brain). Notably, in this critical region, the avalanche size and duration distributions follow power laws. Critical exponents are consistent with a scaling relationship observed recently in neural avalanches measurements. In conclusion, our simple model suggests that avalanche power laws in cortical spontaneous activity may be the effect of a network at the critical point between the replay and non-replay of spatio-temporal patterns.

Time scales in cognitive neuroscience

Cognitive neuroscience boils down to describing the ways in which cognitive function results from brain activity. In turn, brain activity shows complex fluctuations, with structure at many spatio-temporal scales. Exactly how cognitive function inherits the physical dimensions of neural activity, though, is highly non-trivial, and so are generally the corresponding dimensions of cognitive phenomena. As for any physical phenomenon, when studying cognitive function, the first conceptual step should be that of establishing its dimensions. Here, we provide a systematic presentation of the temporal aspects of task-related brain activity, from the smallest scale of the brain imaging technique's resolution, to the observation time of a given experiment, through the characteristic time scales of the process under study. We first review some standard assumptions on the temporal scales of cognitive function. In spite of their general use, these assumptions hold true to a high degree of approximation for many cognitive (viz. fast perceptual) processes, but have their limitations for other ones (e.g., thinking or reasoning). We define in a rigorous way the temporal quantifiers of cognition at all scales, and illustrate how they qualitatively vary as a function of the properties of the cognitive process under study. We propose that each phenomenon should be approached with its own set of theoretical, methodological and analytical tools. In particular, we show that when treating cognitive processes such as thinking or reasoning, complex properties of ongoing brain activity, which can be drastically simplified when considering fast (e.g., perceptual) processes, start playing a major role, and not only characterize the temporal properties of task-related brain activity, but also determine the conditions for proper observation of the phenomena. Finally, some implications on the design of experiments, data analyses, and the choice of recording parameters are discussed.

Neuronal avalanches differ from wakefulness to deep sleep--evidence from intracranial depth recordings in humans

Neuronal activity differs between wakefulness and sleep states. In contrast, an attractor state, called self-organized critical (SOC), was proposed to govern brain dynamics because it allows for optimal information coding. But is the human brain SOC for each vigilance state despite the variations in neuronal dynamics? We characterized neuronal avalanches--spatiotemporal waves of enhanced activity--from dense intracranial depth recordings in humans. We showed that avalanche distributions closely follow a power law--the hallmark feature of SOC--for each vigilance state. However, avalanches clearly differ with vigilance states: slow wave sleep (SWS) shows large avalanches, wakefulness intermediate, and rapid eye movement (REM) sleep small ones. Our SOC model, together with the data, suggested first that the differences are mediated by global but tiny changes in synaptic strength, and second, that the changes with vigilance states reflect small deviations from criticality to the subcritical regime, implying that the human brain does not operate at criticality proper but close to SOC. Independent of criticality, the analysis confirms that SWS shows increased correlations between cortical areas, and reveals that REM sleep shows more fragmented cortical dynamics.

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.