Self-organized criticality as a fundamental property of neural systems

Common neural signatures of psychedelics: Frequency-specific energy changes and repertoire expansion revealed using connectome-harmonic decomposition

The search for the universal laws of human brain function is still on-going but progress is being made. Here we describe the novel concepts of connectome harmonics and connectome-harmonic decomposition, which can be used to characterize the brain activity associated with any mental state. We use this new frequency-specific language to describe the brain activity elicited by psilocybin and LSD and find remarkably similar effects in terms of increases in total energy and power, as well as frequency-specific energy changes and repertoire expansion. In addition, we find enhanced signatures of criticality suggesting that the brain dynamics tune toward criticality in both psychedelic elicited states. Overall, our findings provide new evidence for the remarkable ability of psychedelics to change the spatiotemporal dynamics of the human brain.

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

Understanding the relationship between large-scale structural and functional brain networks remains a crucial issue in modern neuroscience. Recently, there has been growing interest in investigating the role of homeostatic plasticity mechanisms, across different spatiotemporal scales, in regulating network activity and brain functioning against a wide range of environmental conditions and brain states (e.g., during learning, development, ageing, neurological diseases). In the present study, we investigate how the inclusion of homeostatic plasticity in a stochastic whole-brain model, implemented as a normalization of the incoming node's excitatory input, affects the macroscopic activity during rest and the formation of functional networks. Importantly, we address the structure-function relationship both at the group and individual-based levels. In this work, we show that normalization of the node's excitatory input improves the correspondence between simulated neural patterns of the model and various brain functional data. Indeed, we find that the best match is achieved when the model control parameter is in its critical value and that normalization minimizes both the variability of the critical points and neuronal activity patterns among subjects. Therefore, our results suggest that the inclusion of homeostatic principles lead to more realistic brain activity consistent with the hallmarks of criticality. Our theoretical framework open new perspectives in personalized brain modeling with potential applications to investigate the deviation from criticality due to structural lesions (e.g. stroke) or brain disorders.

Synaptic noise facilitates the emergence of self-organized criticality in the Caenorhabditis elegans neuronal network

Avalanches with power-law distributed size parameters have been observed in neuronal networks. This observation might be a manifestation of self-organized criticality (SOC). Yet, the physiological mechanisms of this behaviour are currently unknown. Describing synaptic noise as transmission failures mainly originating from the probabilistic nature of neurotransmitter release, this study investigates the potential of this noise as a mechanism for driving the functional architecture of the neuronal networks towards SOC. To this end, a simple finite state neuron model, with activity dependent and synapse specific failure probabilities, was built based on the known anatomical connectivity data of the nematode Ceanorhabditis elegans. Beginning from random values, it was observed that synaptic noise levels picked out a set of synapses and consequently an active subnetwork that generates power-law distributed neuronal avalanches. The findings of this study bring up the possibility that synaptic failures might be a component of physiological processes underlying SOC in neuronal networks.

Coexistence of scale-invariant and rhythmic behavior in self-organized criticality

Scale-free behavior as well as oscillations are frequently observed in the activity of many natural systems. One important example is the cortical tissues of mammalian brain where both phenomena are simultaneously observed. Rhythmic oscillations as well as critical (scale-free) dynamics are thought to be important, but theoretically incompatible, features of a healthy brain. Motivated by the above, we study the possibility of the coexistence of scale-free avalanches along with rhythmic behavior within the framework of self-organized criticality. In particular, we add an oscillatory perturbation to local threshold condition of the continuous Zhang model and characterize the subsequent activity of the system. We observe regular oscillations embedded in well-defined avalanches which exhibit scale-free size and duration in line with observed neuronal avalanches. The average amplitude of such oscillations are shown to decrease with increasing frequency consistent with real brain oscillations. Furthermore, it is shown that optimal amplification of oscillations occur at the critical point, further providing evidence for functional advantages of criticality.

Characteristic Fluctuations Around Stable Attractor Dynamics Extracted from Highly Nonstationary Electroencephalographic Recordings

Since the discovery of electrical activity of the brain, electroencephalographic (EEG) recordings constitute one of the most popular techniques of brain research. However, EEG signals are highly nonstationary and one should expect that averages of the cross-correlation coefficient, which may take positive and negative values with equal probability, (almost) vanish when estimated over long data segments. Instead, we found that the average zero-lag cross-correlation matrix estimated with a running window over the whole night of sleep EEGs, or of resting state during eyes-open and eyes-closed conditions of healthy subjects shows a characteristic correlation pattern containing pronounced nonzero values. A similar correlation structure has already been encountered in scalp EEG signals containing focal onset seizures. Therefore, we conclude that this structure is independent of the physiological state. Because of its pronounced similarity across subjects, we believe that it depicts a generic feature of the brain dynamics. Namely, we interpret this pattern as a manifestation of a dynamical ground state of the brain activity, necessary to preserve an efficient operational mode, or, expressed in terms of dynamical system theory, we interpret it as a "shadow" of the evolution on (or close to) an attractor in phase space. Nonstationary dynamical aspects of higher cerebral processes should manifest in deviations from this stable pattern. We confirm this hypothesis through a correlation analysis of EEG recordings of 10 healthy subjects during night sleep, 20 recordings of 9 epilepsy patients, and 42 recordings of 21 healthy subjects in resting state during eyes-open and eyes-closed conditions. In particular, we show that the estimation of deviations from the stationary correlation structures provides a more significant differentiation of physiological states and more homogeneous results across subjects.

Beta-Rhythm Oscillations and Synchronization Transition in Network Models of Izhikevich Neurons: Effect of Topology and Synaptic Type

Despite their significant functional roles, beta-band oscillations are least understood. Synchronization in neuronal networks have attracted much attention in recent years with the main focus on transition type. Whether one obtains explosive transition or a continuous transition is an important feature of the neuronal network which can depend on network structure as well as synaptic types. In this study we consider the effect of synaptic interaction (electrical and chemical) as well as structural connectivity on synchronization transition in network models of Izhikevich neurons which spike regularly with beta rhythms. We find a wide range of behavior including continuous transition, explosive transition, as well as lack of global order. The stronger electrical synapses are more conducive to synchronization and can even lead to explosive synchronization. The key network element which determines the order of transition is found to be the clustering coefficient and not the small world effect, or the existence of hubs in a network. These results are in contrast to previous results which use phase oscillator models such as the Kuramoto model. Furthermore, we show that the patterns of synchronization changes when one goes to the gamma band. We attribute such a change to the change in the refractory period of Izhikevich neurons which changes significantly with frequency.

Robustness of functional networks at criticality against structural defects

The robustness of dynamical properties of neuronal networks against structural damages is a central problem in computational and experimental neuroscience. Research has shown that the cortical network of a healthy brain works near a critical state and, moreover, that functional neuronal networks often have scale-free and small-world properties. In this work, we study how the robustness of simple functional networks at criticality is affected by structural defects. In particular, we consider a two-dimensional Ising model at the critical temperature and investigate how its functional network changes with the increasing degree of structural defects. We show that the scale-free and small-world properties of the functional network at criticality are robust against large degrees of structural lesions while the system remains below the percolation limit. Although the Ising model is only a conceptual description of a two-state neuron, our research reveals fundamental robustness properties of functional networks derived from classical statistical mechanics models.

Network Intervention, a Method to Address Complex Therapeutic Strategies

Objective: Network-based approaches emerged as powerful tools for studying complex diseases. Our intention in this article was to raise awareness of the benefits of new therapeutic strategy in biological networks context and provide an introduction to this topic.

Methods: This article will discuss the rational for network intervention, and outline some of the important aspects of deciphering targets activities in the network and future embodiments of network intervention. We also present examples of network intervention based on the strategies these approaches use.

Results: Network intervention seeks for target combinations to perturb a specific subset of nodes in disease networks to inhibit the bypass mechanisms at systems level. Experimental results derived from our studies are discussed, with conclusions that lead to future research directions. A simple diagram is designed to give a way to find the minimum number of external input required for a network intervention based on the graph theory and get the analytical value of the least input.

Conclusion: Creating network intervention that addresses blindness and unthinking action in this way could, therefore, provide more benefit than multi-target therapy. We hope that this article will give readers an appreciation for a new therapeutic strategy that has been proposed for improving clinical benefit by adopting network-based approaches as well as insight into their properties.

Scale invariance in natural and artificial collective systems: a review

Self-organized collective coordinated behaviour is an impressive phenomenon, observed in a variety of natural and artificial systems, in which coherent global structures or dynamics emerge from local interactions between individual parts. If the degree of collective integration of a system does not depend on size, its level of robustness and adaptivity is typically increased and we refer to it as scale-invariant. In this review, we first identify three main types of self-organized scale-invariant systems: scale-invariant spatial structures, scale-invariant topologies and scale-invariant dynamics. We then provide examples of scale invariance from different domains in science, describe their origins and main features and discuss potential challenges and approaches for designing and engineering artificial systems with scale-invariant properties.

Neurofeedback Tunes Scale-Free Dynamics in Spontaneous Brain Activity

Brain oscillations exhibit long-range temporal correlations (LRTCs), which reflect the regularity of their fluctuations: low values representing more random (decorrelated) while high values more persistent (correlated) dynamics. LRTCs constitute supporting evidence that the brain operates near criticality, a state where neuronal activities are balanced between order and randomness. Here, healthy adults used closed-loop brain training (neurofeedback, NFB) to reduce the amplitude of alpha oscillations, producing a significant increase in spontaneous LRTCs post-training. This effect was reproduced in patients with post-traumatic stress disorder, where abnormally random dynamics were reversed by NFB, correlating with significant improvements in hyperarousal. Notably, regions manifesting abnormally low LRTCs (i.e., excessive randomness) normalized toward healthy population levels, consistent with theoretical predictions about self-organized criticality. Hence, when exposed to appropriate training, spontaneous cortical activity reveals a residual capacity for "self-tuning" its own temporal complexity, despite manifesting the abnormal dynamics seen in individuals with psychiatric disorder. Lastly, we observed an inverse-U relationship between strength of LRTC and oscillation amplitude, suggesting a breakdown of long-range dependence at high/low synchronization extremes, in line with recent computational models. Together, our findings offer a broader mechanistic framework for motivating research and clinical applications of NFB, encompassing disorders with perturbed LRTCs.

A Basic Phase Diagram of Neuronal Dynamics

The extreme complexity of the brain has attracted the attention of neuroscientists and other researchers for a long time. More recently, the neuromorphic hardware has matured to provide a new powerful tool to study neuronal dynamics. Here, we study neuronal dynamics using different settings on a neuromorphic chip built with flexible parameters of neuron models. Our unique setting in the network of leaky integrate-and-fire (LIF) neurons is to introduce a weak noise environment. We observed three different types of collective neuronal activities, or phases, separated by sharp boundaries, or phase transitions. From this, we construct a rudimentary phase diagram of neuronal dynamics and demonstrate that a noise-induced chaotic phase (N-phase), which is dominated by neuronal avalanche activity (intermittent aperiodic neuron firing), emerges in the presence of noise and its width grows with the noise intensity. The dynamics can be manipulated in this N-phase. Our results and comparison with clinical data is consistent with the literature and our previous work showing that healthy brain must reside in the N-phase. We argue that the brain phase diagram with further refinement may be used for the diagnosis and treatment of mental disease and also suggest that the dynamics may be manipulated to serve as a means of new information processing (e.g., for optimization). Neuromorphic chips, similar to the one we used but with a variety of neuron models, may be used to further enhance the understanding of human brain function and accelerate the development of neuroscience research.

Efficient Coding and Energy Efficiency Are Promoted by Balanced Excitatory and Inhibitory Synaptic Currents in Neuronal Network

Selective pressure may drive neural systems to process as much information as possible with the lowest energy cost. Recent experiment evidence revealed that the ratio between synaptic excitation and inhibition (E/I) in local cortex is generally maintained at a certain value which may influence the efficiency of energy consumption and information transmission of neural networks. To understand this issue deeply, we constructed a typical recurrent Hodgkin-Huxley network model and studied the general principles that governs the relationship among the E/I synaptic current ratio, the energy cost and total amount of information transmission. We observed in such a network that there exists an optimal E/I synaptic current ratio in the network by which the information transmission achieves the maximum with relatively low energy cost. The coding energy efficiency which is defined as the mutual information divided by the energy cost, achieved the maximum with the balanced synaptic current. Although background noise degrades information transmission and imposes an additional energy cost, we find an optimal noise intensity that yields the largest information transmission and energy efficiency at this optimal E/I synaptic transmission ratio. The maximization of energy efficiency also requires a certain part of energy cost associated with spontaneous spiking and synaptic activities. We further proved this finding with analytical solution based on the response function of bistable neurons, and demonstrated that optimal net synaptic currents are capable of maximizing both the mutual information and energy efficiency. These results revealed that the development of E/I synaptic current balance could lead a cortical network to operate at a highly efficient information transmission rate at a relatively low energy cost. The generality of neuronal models and the recurrent network configuration used here suggest that the existence of an optimal E/I cell ratio for highly efficient energy costs and information maximization is a potential principle for cortical circuit networks.

Supraspinal modulation of neuronal synchronization by nociceptive stimulation induces an enduring reorganization of dorsal horn neuronal connectivity

KEY POINTS

The state of central sensitization induced by the intradermic injection of capsaicin leads to structured (non-random) changes in functional connectivity between dorsal horn neuronal populations distributed along the spinal lumbar segments in anaesthetized cats. The capsaicin-induced changes in neuronal connectivity and the concurrent increase in secondary hyperalgesia are transiently reversed by the systemic administration of small doses of lidocaine, a clinically effective procedure to treat neuropathic pain. The effects of both capsaicin and lidocaine are greatly attenuated in spinalized preparations, showing that supraspinal influences play a significant role in the shaping of nociceptive-induced changes in dorsal horn functional neuronal connectivity. We conclude that changes in functional connectivity between segmental populations of dorsal horn neurones induced by capsaicin and lidocaine result from a cooperative adaptive interaction between supraspinal and spinal neuronal networks, a process that may have a relevant role in the pathogenesis of chronic pain and analgesia.

ABSTRACT

Despite a profusion of information on the molecular and cellular mechanisms involved in the central sensitization produced by intense nociceptive stimulation, the changes in the patterns of functional connectivity between spinal neurones associated with the development of secondary hyperalgesia and allodynia remain largely unknown. Here we show that the state of central sensitization produced by the intradermal injection of capsaicin is associated with structured transformations in neuronal synchronization that lead to an enduring reorganization of the functional connectivity within a segmentally distributed ensemble of dorsal horn neurones. These changes are transiently reversed by the systemic administration of small doses of lidocaine, a clinically effective procedure to treat neuropathic pain. Lidocaine also reduces the capsaicin-induced facilitation of the spinal responses evoked by weak mechanical stimulation of the skin in the region of secondary but not primary hyperalgesia. The effects of both intradermic capsaicin and systemic lidocaine on the segmental correlation and coherence between ongoing cord dorsum potentials and on the responses evoked by tactile stimulation in the region of secondary hyperalgesia are greatly attenuated in spinalized preparations, showing that supraspinal influences are involved in the reorganization of the nociceptive-induced structured patterns of dorsal horn neuronal connectivity. We conclude that the structured reorganization of the functional connectivity between the dorsal horn neurones induced by capsaicin nociceptive stimulation results from cooperative interactions between supraspinal and spinal networks, a process that may have a relevant role in the shaping of the spinal state in the pathogenesis of chronic pain and analgesia.

Determination of the Edge of Criticality in Echo State Networks Through Fisher Information Maximization

It is a widely accepted fact that the computational capability of recurrent neural networks (RNNs) is maximized on the so-called "edge of criticality." Once the network operates in this configuration, it performs efficiently on a specific application both in terms of: 1) low prediction error and 2) high short-term memory capacity. Since the behavior of recurrent networks is strongly influenced by the particular input signal driving the dynamics, a universal, application-independent method for determining the edge of criticality is still missing. In this paper, we aim at addressing this issue by proposing a theoretically motivated, unsupervised method based on Fisher information for determining the edge of criticality in RNNs. It is proved that Fisher information is maximized for (finite-size) systems operating in such critical regions. However, Fisher information is notoriously difficult to compute and requires the analytic form of the probability density function ruling the system behavior. This paper takes advantage of a recently developed nonparametric estimator of the Fisher information matrix and provides a method to determine the critical region of echo state networks (ESNs), a particular class of recurrent networks. The considered control parameters, which indirectly affect the ESN performance, are explored to identify those configurations lying on the edge of criticality and, as such, maximizing Fisher information and computational performance. Experimental results on benchmarks and real-world data demonstrate the effectiveness of the proposed method.

Answering Schrödinger's question: A free-energy formulation

The free-energy principle (FEP) is a formal model of neuronal processes that is widely recognised in neuroscience as a unifying theory of the brain and biobehaviour. More recently, however, it has been extended beyond the brain to explain the dynamics of living systems, and their unique capacity to avoid decay. The aim of this review is to synthesise these advances with a meta-theoretical ontology of biological systems called variational neuroethology, which integrates the FEP with Tinbergen's four research questions to explain biological systems across spatial and temporal scales. We exemplify this framework by applying it to Homo sapiens, before translating variational neuroethology into a systematic research heuristic that supplies the biological, cognitive, and social sciences with a computationally tractable guide to discovery.

Quantum coherence-driven self-organized criticality and nonequilibrium light localization

Self-organized criticality emerges in dynamical complex systems driven out of equilibrium and characterizes a wide range of classical phenomena in physics, geology, and biology. We report on a quantum coherence-controlled self-organized critical transition observed in the light localization behavior of a coherence-driven nanophotonic configuration. Our system is composed of a gain-enhanced plasmonic heterostructure controlled by a coherent drive, in which photons close to the stopped-light regime interact in the presence of the active nonlinearities, eventually synchronizing their dynamics. In this system, on the basis of analytical and corroborating full-wave Maxwell-Bloch computations, we observe quantum coherence-controlled self-organized criticality in the emergence of light localization arising from the synchronization of the photons. It is associated with two first-order phase transitions: one pertaining to the synchronization of the dynamics of the photons and the second pertaining to an inversionless lasing transition by the coherent drive. The so-attained light localization, which is robust to dissipation, fluctuations, and many-body interactions, exhibits scale-invariant power laws and absence of finely tuned control parameters. We also found that, in this nonequilibrium dynamical system, the effective critical "temperature" of the system drops to zero, whereupon one enters the quantum self-organized critical regime.

Optimization by Self-Organized Criticality

Self-organized criticality (SOC) is a phenomenon observed in certain complex systems of multiple interacting components, e.g., neural networks, forest fires, and power grids, that produce power-law distributed avalanche sizes. Here, we report the surprising result that the avalanches from an SOC process can be used to solve non-convex optimization problems. To generate avalanches, we use the Abelian sandpile model on a graph that mirrors the graph of the optimization problem. For optimization, we map the avalanche areas onto search patterns for optimization, while the SOC process receives no feedback from the optimization itself. The resulting method can be applied without parameter tuning to a wide range of optimization problems, as demonstrated on three problems: finding the ground-state of an Ising spin glass, graph coloring, and image segmentation. We find that SOC search is more efficient compared to other random search methods, including simulated annealing, and unlike annealing, it is parameter free, thereby eliminating the time-consuming requirement to tune an annealing temperature schedule.

Catecholamines alter the intrinsic variability of cortical population activity and perception

The ascending modulatory systems of the brain stem are powerful regulators of global brain state. Disturbances of these systems are implicated in several major neuropsychiatric disorders. Yet, how these systems interact with specific neural computations in the cerebral cortex to shape perception, cognition, and behavior remains poorly understood. Here, we probed into the effect of two such systems, the catecholaminergic (dopaminergic and noradrenergic) and cholinergic systems, on an important aspect of cortical computation: its intrinsic variability. To this end, we combined placebo-controlled pharmacological intervention in humans, recordings of cortical population activity using magnetoencephalography (MEG), and psychophysical measurements of the perception of ambiguous visual input. A low-dose catecholaminergic, but not cholinergic, manipulation altered the rate of spontaneous perceptual fluctuations as well as the temporal structure of "scale-free" population activity of large swaths of the visual and parietal cortices. Computational analyses indicate that both effects were consistent with an increase in excitatory relative to inhibitory activity in the cortical areas underlying visual perceptual inference. We propose that catecholamines regulate the variability of perception and cognition through dynamically changing the cortical excitation-inhibition ratio. The combined readout of fluctuations in perception and cortical activity we established here may prove useful as an efficient and easily accessible marker of altered cortical computation in neuropsychiatric disorders.

Locating Order-Disorder Phase Transition in a Cardiac System

To prevent sudden cardiac death, predicting where in the cardiac system an order-disorder phase transition into ventricular fibrillation begins is as important as when it begins. We present a computationally efficient, information-theoretic approach to predicting the locations of the wavebreaks. Such wavebreaks initiate fibrillation in a cardiac system where the order-disorder behavior is controlled by a single driving component, mimicking electrical misfiring from the pulmonary veins or from the Purkinje fibers. Communication analysis between the driving component and each component of the system reveals that channel capacity, mutual information and transfer entropy can locate the wavebreaks. This approach is applicable to interventional therapies to prevent sudden death, and to a wide range of systems to mitigate or prevent imminent phase transitions.

Inferring multi-scale neural mechanisms with brain network modelling

The neurophysiological processes underlying non-invasive brain activity measurements are incompletely understood. Here, we developed a connectome-based brain network model that integrates individual structural and functional data with neural population dynamics to support multi-scale neurophysiological inference. Simulated populations were linked by structural connectivity and, as a novelty, driven by electroencephalography (EEG) source activity. Simulations not only predicted subjects' individual resting-state functional magnetic resonance imaging (fMRI) time series and spatial network topologies over 20 minutes of activity, but more importantly, they also revealed precise neurophysiological mechanisms that underlie and link six empirical observations from different scales and modalities: (1) resting-state fMRI oscillations, (2) functional connectivity networks, (3) excitation-inhibition balance, (4, 5) inverse relationships between α-rhythms, spike-firing and fMRI on short and long time scales, and (6) fMRI power-law scaling. These findings underscore the potential of this new modelling framework for general inference and integration of neurophysiological knowledge to complement empirical studies.

Neuroscientists can use various techniques to measure activity within the brain without opening up the skull. One of the most common is electroencephalography, or EEG for short. A net of electrodes is attached to the scalp and reveals the patterns of electrical activity occurring in brain tissue. But while EEG is good at revealing electrical activity across the surface of the scalp, it is less effective at linking the observed activity to specific locations in the brain.

Another widely used technique is functional magnetic resonance imaging, or fMRI. A patient, or healthy volunteer, lies inside a scanner containing a large magnet. The scanner tracks changes in the level of oxygen at different regions of the brain to provide a measure of how the activity of these regions changes over time. In contrast to EEG, fMRI is good at pinpointing the location of brain activity, but it is an indirect measure of brain activity as it depends on blood flow and several other factors. In terms of understanding how the brain works, EEG and fMRI thus provide different pieces of the puzzle. But there is no easy way to fit these pieces together.

Other areas of science have used computer models to merge different sources of data to obtain new insights into complex processes. Schirner et al. now adopt this approach to reveal the workings of the brain that underly signals like EEG and fMRI.

After recording structural MRI data from healthy volunteers, Schirner et al. built a computer model of each person’s brain. They then ran simulations with each individual model stimulating it with the person’s EEG to predict the fMRI activity of the same individual. Comparing these predictions with real fMRI data collected at the same time as the EEG confirmed that the predictions were accurate. Importantly, the brain models also displayed many features of neural activity that previously could only be measured by implanting electrodes into the brain.

This new approach provides a way of combining experimental data with theories about how the nervous system works. The resulting models can help generate and test ideas about the mechanisms underlying brain activity. Building models of different brains based on data from individual people could also help reveal the biological basis of differences between individuals. This could in turn provide insights into why some individuals are more vulnerable to certain brain diseases and open up new ways to treat these diseases.

Critical and Supercritical Spatiotemporal Calcium Dynamics in Beta Cells

A coordinated functioning of beta cells within pancreatic islets is mediated by oscillatory membrane depolarization and subsequent changes in cytoplasmic calcium concentration. While gap junctions allow for intraislet information exchange, beta cells within islets form complex syncytia that are intrinsically nonlinear and highly heterogeneous. To study spatiotemporal calcium dynamics within these syncytia, we make use of computational modeling and confocal high-speed functional multicellular imaging. We show that model predictions are in good agreement with experimental data, especially if a high degree of heterogeneity in the intercellular coupling term is assumed. In particular, during the first few minutes after stimulation, the probability distribution of calcium wave sizes is characterized by a power law, thus indicating critical behavior. After this period, the dynamics changes qualitatively such that the number of global intercellular calcium events increases to the point where the behavior becomes supercritical. To better mimic normal in vivo conditions, we compare the described behavior during supraphysiological non-oscillatory stimulation with the behavior during exposure to a slightly lower and oscillatory glucose challenge. In the case of this protocol, we observe only critical behavior in both experiment and model. Our results indicate that the loss of oscillatory changes, along with the rise in plasma glucose observed in diabetes, could be associated with a switch to supercritical calcium dynamics and loss of beta cell functionality.

Connectome-harmonic decomposition of human brain activity reveals dynamical repertoire re-organization under LSD

Recent studies have started to elucidate the effects of lysergic acid diethylamide (LSD) on the human brain but the underlying dynamics are not yet fully understood. Here we used 'connectome-harmonic decomposition', a novel method to investigate the dynamical changes in brain states. We found that LSD alters the energy and the power of individual harmonic brain states in a frequency-selective manner. Remarkably, this leads to an expansion of the repertoire of active brain states, suggestive of a general re-organization of brain dynamics given the non-random increase in co-activation across frequencies. Interestingly, the frequency distribution of the active repertoire of brain states under LSD closely follows power-laws indicating a re-organization of the dynamics at the edge of criticality. Beyond the present findings, these methods open up for a better understanding of the complex brain dynamics in health and disease.

A model of individualized canonical microcircuits supporting cognitive operations

Major cognitive functions such as language, memory, and decision-making are thought to rely on distributed networks of a large number of basic elements, called canonical microcircuits. In this theoretical study we propose a novel canonical microcircuit model and find that it supports two basic computational operations: a gating mechanism and working memory. By means of bifurcation analysis we systematically investigate the dynamical behavior of the canonical microcircuit with respect to parameters that govern the local network balance, that is, the relationship between excitation and inhibition, and key intrinsic feedback architectures of canonical microcircuits. We relate the local behavior of the canonical microcircuit to cognitive processing and demonstrate how a network of interacting canonical microcircuits enables the establishment of spatiotemporal sequences in the context of syntax parsing during sentence comprehension. This study provides a framework for using individualized canonical microcircuits for the construction of biologically realistic networks supporting cognitive operations.

Comparative approximations of criticality in a neural and quantum regime

Under a variety of conditions, stochastic and non-linear systems with many degrees of freedom tend to evolve towards complexity and criticality. Over the last decades, a steady proliferation of models re: far-from-equilibrium thermodynamics of metastable, many-valued systems arose, serving as attributes of a 'critical' attractor landscape. Building off recent data citing trademark aspects of criticality in the brain-including: power-laws, scale-free (1/f) behavior (scale invariance, or scale independence), critical slowing, and avalanches-it has been conjectured that operating at criticality entails functional advantages such as: optimized neural computation and information processing; boosted memory; large dynamical ranges; long-range communication; and an increased ability to react to highly diverse stimuli. In short, critical dynamics provide a necessary condition for neurobiologically significant elements of brain dynamics. Theoretical predictions have been verified in specific models such as Boolean networks, liquid state machines, and neural networks. These findings inspired the neural criticality hypothesis, proposing that the brain operates in a critical state because the associated optimal computational capabilities provide an evolutionarily advantage. This paper develops in three parts: after developing the critical landscape, we will then shift gears to rediscover another inroad to criticality via stochastic quantum field theory and dissipative dynamics. The existence of these two approaches deserves some consideration, given both neural and quantum criticality hypotheses propose specific mechanisms that leverage the same phenomena. This suggests that understanding the quantum approach could help to shed light on brain-based modeling. In the third part, we will turn to Whitehead's actual entities and modes of perception in order to demonstrate a concomitant logic underwriting both models. In the discussion, I briefly motivate a reading of criticality and its properties as responsive to the characterization of tenets from Eastern wisdom traditions.

Temporal profiles of avalanches on networks

An avalanche or cascade occurs when one event causes one or more subsequent events, which in turn may cause further events in a chain reaction. Avalanching dynamics are studied in many disciplines, with a recent focus on average avalanche shapes, i.e., the temporal profiles of avalanches of fixed duration. At the critical point of the dynamics, the rescaled average avalanche shapes for different durations collapse onto a single universal curve. We apply Markov branching process theory to derive an equation governing the average avalanche shape for cascade dynamics on networks. Analysis of the equation at criticality demonstrates that nonsymmetric average avalanche shapes (as observed in some experiments) occur for certain combinations of dynamics and network topology. We give examples using numerical simulations of models for information spreading, neural dynamics, and behavior adoption and we propose simple experimental tests to quantify whether cascading systems are in the critical state.

Mechanisms of self-organized criticality in social processes of knowledge creation

In online social dynamics, a robust scale invariance appears as a key feature of collaborative efforts that lead to new social value. The underlying empirical data thus offers a unique opportunity to study the origin of self-organized criticality (SOC) in social systems. In contrast to physical systems in the laboratory, various human attributes of the actors play an essential role in the process along with the contents (cognitive, emotional) of the communicated artifacts. As a prototypical example, we consider the social endeavor of knowledge creation via Questions and Answers (Q&A). Using a large empirical data set from one of such Q&A sites and theoretical modeling, we reveal fundamental characteristics of SOC by investigating the temporal correlations at all scales and the role of cognitive contents to the avalanches of the knowledge-creation process. Our analysis shows that the universal social dynamics with power-law inhomogeneities of the actions and delay times provides the primary mechanism for self-tuning towards the critical state; it leads to the long-range correlations and the event clustering in response to the external driving by the arrival of new users. In addition, the involved cognitive contents (systematically annotated in the data and observed in the model) exert important constraints that identify unique classes of the knowledge-creation avalanches. Specifically, besides determining a fine structure of the developing knowledge networks, they affect the values of scaling exponents and the geometry of large avalanches and shape the multifractal spectrum. Furthermore, we find that the level of the activity of the communities that share the knowledge correlates with the fluctuations of the innovation rate, implying that the increase of innovation may serve as the active principle of self-organization. To identify relevant parameters and unravel the role of the network evolution underlying the process in the social system under consideration, we compare the social avalanches to the avalanche sequences occurring in the field-driven physical model of disordered solids, where the factors contributing to the collective dynamics are better understood.

Criticality predicts maximum irregularity in recurrent networks of excitatory nodes

A rigorous understanding of brain dynamics and function requires a conceptual bridge between multiple levels of organization, including neural spiking and network-level population activity. Mounting evidence suggests that neural networks of cerebral cortex operate at a critical regime, which is defined as a transition point between two phases of short lasting and chaotic activity. However, despite the fact that criticality brings about certain functional advantages for information processing, its supporting evidence is still far from conclusive, as it has been mostly based on power law scaling of size and durations of cascades of activity. Moreover, to what degree such hypothesis could explain some fundamental features of neural activity is still largely unknown. One of the most prevalent features of cortical activity in vivo is known to be spike irregularity of spike trains, which is measured in terms of the coefficient of variation (CV) larger than one. Here, using a minimal computational model of excitatory nodes, we show that irregular spiking (CV > 1) naturally emerges in a recurrent network operating at criticality. More importantly, we show that even at the presence of other sources of spike irregularity, being at criticality maximizes the mean coefficient of variation of neurons, thereby maximizing their spike irregularity. Furthermore, we also show that such a maximized irregularity results in maximum correlation between neuronal firing rates and their corresponding spike irregularity (measured in terms of CV). On the one hand, using a model in the universality class of directed percolation, we propose new hallmarks of criticality at single-unit level, which could be applicable to any network of excitable nodes. On the other hand, given the controversy of the neural criticality hypothesis, we discuss the limitation of this approach to neural systems and to what degree they support the criticality hypothesis in real neural networks. Finally, we discuss the limitations of applying our results to real networks and to what degree they support the criticality hypothesis.

Self-organizing Complex Networks: individual versus global rules

We introduce a form of Self-Organized Criticality (SOC) inspired by the new generation of evolutionary game theory, which ranges from physiology to sociology. The single individuals are the nodes of a composite network, equivalent to two interacting subnetworks, one leading to strategy choices made by the individuals under the influence of the choices of their nearest neighbors and the other measuring the Prisoner's Dilemma Game payoffs of these choices. The interaction between the two networks is established by making the imitation strength K increase or decrease according to whether the last two payoffs increase or decrease upon increasing or decreasing K. Although each of these imitation strengths is selected selfishly, and independently of the others as well, the social system spontaneously evolves toward the state of cooperation. Criticality is signaled by temporal complexity, namely the occurrence of non-Poisson renewal events, the time intervals between two consecutive crucial events being given by an inverse power law index μ = 1.3 rather than by avalanches with an inverse power law distribution as in the original form of SOC. This new phenomenon is herein labeled self-organized temporal criticality (SOTC). We compare this bottom-up self-organization process to the adoption of a global choice rule based on assigning to all the units the same value K, with the time evolution of common K being determined by consciousness of the social benefit, a top-down process implying the action of a leader. In this case self-organization is impeded by large intensity fluctuations and the global social benefit turns out to be much weaker. We conclude that the SOTC model fits the requests of a manifesto recently proposed by a number of European social scientists.

Biological modelling of a computational spiking neural network with neuronal avalanches

In recent years, an increasing number of studies have demonstrated that networks in the brain can self-organize into a critical state where dynamics exhibit a mixture of ordered and disordered patterns. This critical branching phenomenon is termed neuronal avalanches. It has been hypothesized that the homeostatic level balanced between stability and plasticity of this critical state may be the optimal state for performing diverse neural computational tasks. However, the critical region for high performance is narrow and sensitive for spiking neural networks (SNNs). In this paper, we investigated the role of the critical state in neural computations based on liquid-state machines, a biologically plausible computational neural network model for real-time computing. The computational performance of an SNN when operating at the critical state and, in particular, with spike-timing-dependent plasticity for updating synaptic weights is investigated. The network is found to show the best computational performance when it is subjected to critical dynamic states. Moreover, the active-neuron-dominant structure refined from synaptic learning can remarkably enhance the robustness of the critical state and further improve computational accuracy. These results may have important implications in the modelling of spiking neural networks with optimal computational performance.This article is part of the themed issue 'Mathematical methods in medicine: neuroscience, cardiology and pathology'.

Avalanche and edge-of-chaos criticality do not necessarily co-occur in neural networks

There are indications that for optimizing neural computation, neural networks may operate at criticality. Previous approaches have used distinct fingerprints of criticality, leaving open the question whether the different notions would necessarily reflect different aspects of one and the same instance of criticality, or whether they could potentially refer to distinct instances of criticality. In this work, we choose avalanche criticality and edge-of-chaos criticality and demonstrate for a recurrent spiking neural network that avalanche criticality does not necessarily entrain dynamical edge-of-chaos criticality. This suggests that the different fingerprints may pertain to distinct phenomena.

Chaos and Correlated Avalanches in Excitatory Neural Networks with Synaptic Plasticity

A collective chaotic phase with power law scaling of activity events is observed in a disordered mean field network of purely excitatory leaky integrate-and-fire neurons with short-term synaptic plasticity. The dynamical phase diagram exhibits two transitions from quasisynchronous and asynchronous regimes to the nontrivial, collective, bursty regime with avalanches. In the homogeneous case without disorder, the system synchronizes and the bursty behavior is reflected into a period doubling transition to chaos for a two dimensional discrete map. Numerical simulations show that the bursty chaotic phase with avalanches exhibits a spontaneous emergence of persistent time correlations and enhanced Kolmogorov complexity. Our analysis reveals a mechanism for the generation of irregular avalanches that emerges from the combination of disorder and deterministic underlying chaotic dynamics.

Adaptive Local Information Transfer in Random Boolean Networks

Living systems such as gene regulatory networks and neuronal networks have been supposed to work close to dynamical criticality, where their information-processing ability is optimal at the whole-system level. We investigate how this global information-processing optimality is related to the local information transfer at each individual-unit level. In particular, we introduce an internal adjustment process of the local information transfer and examine whether the former can emerge from the latter. We propose an adaptive random Boolean network model in which each unit rewires its incoming arcs from other units to balance stability of its information processing based on the measurement of the local information transfer pattern. First, we show numerically that random Boolean networks can self-organize toward near dynamical criticality in our model. Second, the proposed model is analyzed by a mean-field theory. We recognize that the rewiring rule has a bootstrapping feature. The stationary indegree distribution is calculated semi-analytically and is shown to be close to dynamical criticality in a broad range of model parameter values.

Power-law statistics and universal scaling in the absence of criticality

Critical states are sometimes identified experimentally through power-law statistics or universal scaling functions. We show here that such features naturally emerge from networks in self-sustained irregular regimes away from criticality. In these regimes, statistical physics theory of large interacting systems predict a regime where the nodes have independent and identically distributed dynamics. We thus investigated the statistics of a system in which units are replaced by independent stochastic surrogates and found the same power-law statistics, indicating that these are not sufficient to establish criticality. We rather suggest that these are universal features of large-scale networks when considered macroscopically. These results put caution on the interpretation of scaling laws found in nature.

Dissipative stochastic sandpile model on small-world networks: Properties of nondissipative and dissipative avalanches

A dissipative stochastic sandpile model is constructed and studied on small-world networks in one and two dimensions with different shortcut densities ϕ, where ϕ=0 represents regular lattice and ϕ=1 represents random network. The effect of dimension, network topology, and specific dissipation mode (bulk or boundary) on the the steady-state critical properties of nondissipative and dissipative avalanches along with all avalanches are analyzed. Though the distributions of all avalanches and nondissipative avalanches display stochastic scaling at ϕ=0 and mean-field scaling at ϕ=1, the dissipative avalanches display nontrivial critical properties at ϕ=0 and 1 in both one and two dimensions. In the small-world regime (2^{-12}≤ϕ≤0.1), the size distributions of different types of avalanches are found to exhibit more than one power-law scaling with different scaling exponents around a crossover toppling size s_{c}. Stochastic scaling is found to occur for s<s_{c} and the mean-field scaling is found to occur for s>s_{c}. As different scaling forms are found to coexist in a single probability distribution, a coexistence scaling theory on small world network is developed and numerically verified.

Could a Neuroscientist Understand a Microprocessor?

There is a popular belief in neuroscience that we are primarily data limited, and that producing large, multimodal, and complex datasets will, with the help of advanced data analysis algorithms, lead to fundamental insights into the way the brain processes information. These datasets do not yet exist, and if they did we would have no way of evaluating whether or not the algorithmically-generated insights were sufficient or even correct. To address this, here we take a classical microprocessor as a model organism, and use our ability to perform arbitrary experiments on it to see if popular data analysis methods from neuroscience can elucidate the way it processes information. Microprocessors are among those artificial information processing systems that are both complex and that we understand at all levels, from the overall logical flow, via logical gates, to the dynamics of transistors. We show that the approaches reveal interesting structure in the data but do not meaningfully describe the hierarchy of information processing in the microprocessor. This suggests current analytic approaches in neuroscience may fall short of producing meaningful understanding of neural systems, regardless of the amount of data. Additionally, we argue for scientists using complex non-linear dynamical systems with known ground truth, such as the microprocessor as a validation platform for time-series and structure discovery methods.

More Severe Insomnia Complaints in People with Stronger Long-Range Temporal Correlations in Wake Resting-State EEG

The complaints of people suffering from Insomnia Disorder (ID) concern both sleep and daytime functioning. However, little is known about wake brain temporal dynamics in people with ID. We therefore assessed possible alterations in Long-Range Temporal Correlations (LRTC) in the amplitude fluctuations of band-filtered oscillations in electroencephalography (EEG) recordings. We investigated whether LRTC differ between cases with ID and matched controls. Within both groups, we moreover investigated whether individual differences in subjective insomnia complaints are associated with LRTC. Resting-state high-density EEG (256-channel) was recorded in 52 participants with ID and 43 age- and sex-matched controls, during Eyes Open (EO) and Eyes Closed (EC). Detrended fluctuation analysis was applied to the amplitude envelope of band-filtered EEG oscillations (theta, alpha, sigma, beta-1, beta-2) to obtain the Hurst exponents (H), as measures of LRTC. Participants rated their subjective insomnia complaints using the Insomnia Severity Index (ISI). Through general linear models, we evaluated whether H, aggregated across electrodes and frequencies, differed between cases and controls, or showed within-group associations with individual differences in ISI. Additionally, we characterized the spatio-spectral profiles of group differences and associations using non-parametric statistics. H did not differ between cases with ID and controls in any of the frequency bands, neither during EO nor EC. During EO, however, within-group associations between H and ISI indicated that individuals who experienced worse sleep quality had stronger LRTC. Spatio-spectral profiles indicated that the associations held most prominently for the amplitude fluctuations of parietal theta oscillations within the ID group, and of centro-frontal beta-1 oscillations in controls. While people suffering from insomnia experience substantially worse sleep quality than controls, their brain dynamics express similar strength of LRTC. In each group, however, individuals experiencing worse sleep quality tend to have stronger LRTC during eyes open wakefulness, in a spatio-spectral range specific for each group. Taken together, the findings indicate that subjective insomnia complaints involve distinct dynamical processes in people with ID and controls. The findings are in agreement with recent reports on decreasing LRTC with sleep depth, and with the hypothesis that sleep balances brain excitability.

Phase transitions and self-organized criticality in networks of stochastic spiking neurons

Phase transitions and critical behavior are crucial issues both in theoretical and experimental neuroscience. We report analytic and computational results about phase transitions and self-organized criticality (SOC) in networks with general stochastic neurons. The stochastic neuron has a firing probability given by a smooth monotonic function Φ(V) of the membrane potential V, rather than a sharp firing threshold. We find that such networks can operate in several dynamic regimes (phases) depending on the average synaptic weight and the shape of the firing function Φ. In particular, we encounter both continuous and discontinuous phase transitions to absorbing states. At the continuous transition critical boundary, neuronal avalanches occur whose distributions of size and duration are given by power laws, as observed in biological neural networks. We also propose and test a new mechanism to produce SOC: the use of dynamic neuronal gains - a form of short-term plasticity probably located at the axon initial segment (AIS) - instead of depressing synapses at the dendrites (as previously studied in the literature). The new self-organization mechanism produces a slightly supercritical state, that we called SOSC, in accord to some intuitions of Alan Turing.

Emergence and maintenance of excitability: kinetics over structure

The capacity to generate action potentials in neurons and other excitable cells requires tuning of both ionic channel expression and kinetics in a large parameter space. Alongside studies that extend traditional focus on control-based regulation of structural parameters (channel densities), there is a budding interest in self-organization of kinetic parameters. In this picture, ionic channels are continually forced by activity in-and-out of a pool of states not available for the mechanism of excitability. The process, acting on expressed structure, provides a bed for generation of a spectrum of excitability modes. Driven by microscopic fluctuations over a broad range of temporal scales, self-organization of kinetic parameters enriches the concepts and tools used in the study of development of excitability.

Critical phenomena at a first-order phase transition in a lattice of glow lamps: Experimental findings and analogy to neural activity

Networks of non-linear electronic oscillators have shown potential as physical models of neural dynamics. However, two properties of brain activity, namely, criticality and metastability, remain under-investigated with this approach. Here, we present a simple circuit that exhibits both phenomena. The apparatus consists of a two-dimensional square lattice of capacitively coupled glow (neon) lamps. The dynamics of lamp breakdown (flash) events are controlled by a DC voltage globally connected to all nodes via fixed resistors. Depending on this parameter, two phases having distinct event rate and degree of spatiotemporal order are observed. The transition between them is hysteretic, thus a first-order one, and it is possible to enter a metastability region, wherein, approaching a spinodal point, critical phenomena emerge. Avalanches of events occur according to power-law distributions having exponents ≈3/2 for size and ≈2 for duration, and fractal structure is evident as power-law scaling of the Fano factor. These critical exponents overlap observations in biological neural networks; hence, this circuit may have value as building block to realize corresponding physical models.

The topology of large Open Connectome networks for the human brain

The structural human connectome (i.e. the network of fiber connections in the brain) can be analyzed at ever finer spatial resolution thanks to advances in neuroimaging. Here we analyze several large data sets for the human brain network made available by the Open Connectome Project. We apply statistical model selection to characterize the degree distributions of graphs containing up to nodes and edges. A three-parameter generalized Weibull (also known as a stretched exponential) distribution is a good fit to most of the observed degree distributions. For almost all networks, simple power laws cannot fit the data, but in some cases there is statistical support for power laws with an exponential cutoff. We also calculate the topological (graph) dimension D and the small-world coefficient σ of these networks. While σ suggests a small-world topology, we found that D < 4 showing that long-distance connections provide only a small correction to the topology of the embedding three-dimensional space.

Repertoires of Spike Avalanches Are Modulated by Behavior and Novelty

Neuronal avalanches measured as consecutive bouts of thresholded field potentials represent a statistical signature that the brain operates near a critical point. In theory, criticality optimizes stimulus sensitivity, information transmission, computational capability and mnemonic repertoires size. Field potential avalanches recorded via multielectrode arrays from cortical slice cultures are repeatable spatiotemporal activity patterns. It remains unclear whether avalanches of action potentials observed in forebrain regions of freely-behaving rats also form recursive repertoires, and whether these have any behavioral relevance. Here, we show that spike avalanches, recorded from hippocampus (HP) and sensory neocortex of freely-behaving rats, constitute distinct families of recursive spatiotemporal patterns. A significant number of those patterns were specific to a behavioral state. Although avalanches produced during sleep were mostly similar to others that occurred during waking, the repertoire of patterns recruited during sleep differed significantly from that of waking. More importantly, exposure to novel objects increased the rate at which new patterns arose, also leading to changes in post-exposure repertoires, which were significantly different from those before the exposure. A significant number of families occurred exclusively during periods of whisker contact with objects, but few were associated with specific objects. Altogether, the results provide original evidence linking behavior and criticality at the spike level: spike avalanches form repertoires that emerge in waking, recur during sleep, are diversified by novelty and contribute to object representation.

Stochastic Oscillation in Self-Organized Critical States of Small Systems: Sensitive Resting State in Neural Systems

Self-organized critical states (SOCs) and stochastic oscillations (SOs) are simultaneously observed in neural systems, which appears to be theoretically contradictory since SOCs are characterized by scale-free avalanche sizes but oscillations indicate typical scales. Here, we show that SOs can emerge in SOCs of small size systems due to temporal correlation between large avalanches at the finite-size cutoff, resulting from the accumulation-release process in SOCs. In contrast, the critical branching process without accumulation-release dynamics cannot exhibit oscillations. The reconciliation of SOCs and SOs is demonstrated both in the sandpile model and robustly in biologically plausible neuronal networks. The oscillations can be suppressed if external inputs eliminate the prominent slow accumulation process, providing a potential explanation of the widely studied Berger effect or event-related desynchronization in neural response. The features of neural oscillations and suppression are confirmed during task processing in monkey eye-movement experiments. Our results suggest that finite-size, columnar neural circuits may play an important role in generating neural oscillations around the critical states, potentially enabling functional advantages of both SOCs and oscillations for sensitive response to transient stimuli.

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 connectivity in in vitro neuronal assemblies

Complex network topologies represent the necessary substrate to support complex brain functions. In this work, we reviewed in vitro neuronal networks coupled to Micro-Electrode Arrays (MEAs) as biological substrate. Networks of dissociated neurons developing in vitro and coupled to MEAs, represent a valid experimental model for studying the mechanisms governing the formation, organization and conservation of neuronal cell assemblies. In this review, we present some examples of the use of statistical Cluster Coefficients and Small World indices to infer topological rules underlying the dynamics exhibited by homogeneous and engineered neuronal networks.

Structural determinants of criticality in biological networks

Many adaptive evolutionary systems display spatial and temporal features, such as long-range correlations, typically associated with the critical point of a phase transition in statistical physics. Empirical and theoretical studies suggest that operating near criticality enhances the functionality of biological networks, such as brain and gene networks, in terms for instance of information processing, robustness, and evolvability. While previous studies have explained criticality with specific system features, we still lack a general theory of critical behavior in biological systems. Here we look at this problem from the complex systems perspective, since in principle all critical biological circuits have in common the fact that their internal organization can be described as a complex network. An important question is how self-similar structure influences self-similar dynamics. Modularity and heterogeneity, for instance, affect the location of critical points and can be used to tune the system toward criticality. We review and discuss recent studies on the criticality of neuronal and genetic networks, and discuss the implications of network theory when assessing the evolutionary features of criticality.