Collapse of complexity of brain and body activity due to excessive inhibition and MeCP2 disruption.

Complex body movements require complex dynamics and coordination among neurons in motor cortex. Conversely, a long-standing theoretical notion supposes that if many neurons in motor cortex become excessively synchronized, they may lack the necessary complexity for healthy motor coding. However, direct experimental support for this idea is rare and underlying mechanisms are unclear. Here we recorded three-dimensional body movements and spiking activity of many single neurons in motor cortex of rats with enhanced synaptic inhibition and a transgenic rat model of Rett syndrome (RTT). For both cases, we found a collapse of complexity in the motor system. Reduced complexity was apparent in lower-dimensional, stereotyped brain-body interactions, neural synchrony, and simpler behavior. Our results demonstrate how imbalanced inhibition can cause excessive synchrony among movement-related neurons and, consequently, a stereotyped motor code. Excessive inhibition and synchrony may underlie abnormal motor function in RTT.

Odor modality is transmitted to cortical brain regions from the olfactory bulb.

Odor perception is the impetus for important animal behaviors with two predominate modes of processing: odors pass through the front of the nose (orthonasal) while inhaling and sniffing, or through the rear (retronasal) during exhalation and while eating. Despite the importance of olfaction for an animal's well-being and that ortho and retro naturally occur, it is unknown how the modality (ortho vs. retro) is even transmitted to cortical brain regions, which could significantly affect how odors are processed and perceived. Using multielectrode array recordings in tracheotomized anesthetized rats, which decouples ortho-retro modality from breathing, we show that mitral cells in rat olfactory bulb can reliably and directly transmit orthonasal versus retronasal modality with ethyl butyrate, a common food odor. Drug manipulations affecting synaptic inhibition via GABAA lead to worse decoding of ortho versus retro, independent of whether overall inhibition increases or decreases, suggesting that the olfactory bulb circuit may naturally favor encoding this important aspect of odors. Detailed data analysis paired with a firing rate model that captures population trends in spiking statistics shows how this circuit can encode odor modality. We have not only demonstrated that ortho/retro information is encoded to downstream brain regions but also used modeling to demonstrate a plausible mechanism for this encoding; due to synaptic adaptation, it is the slower time course of the retronasal stimulation that causes retronasal responses to be stronger and less sensitive to inhibitory drug manipulations than orthonasal responses.NEW & NOTEWORTHY Whether ortho (sniffing odors) versus retro (exhalation and eating) is encoded from the olfactory bulb to other brain areas is not completely known. Using multielectrode array recordings in anesthetized rats, we show that the olfactory bulb transmits this information downstream via spikes. Altering inhibition degrades ortho/retro information on average. We use theory and computation to explain our results, which should have implications on cortical processing considering that only food odors occur retronasally.

Differences in olfactory bulb mitral cell spiking with ortho- and retronasal stimulation revealed by data-driven models.

The majority of olfaction studies focus on orthonasal stimulation where odors enter via the front nasal cavity, while retronasal olfaction, where odors enter the rear of the nasal cavity during feeding, is understudied. The coding of retronasal odors via coordinated spiking of neurons in the olfactory bulb (OB) is largely unknown despite evidence that higher level processing is different than orthonasal. To this end, we use multi-electrode array in vivo recordings of rat OB mitral cells (MC) in response to a food odor with both modes of stimulation, and find significant differences in evoked firing rates and spike count covariances (i.e., noise correlations). Differences in spiking activity often have implications for sensory coding, thus we develop a single-compartment biophysical OB model that is able to reproduce key properties of important OB cell types. Prior experiments in olfactory receptor neurons (ORN) showed retro stimulation yields slower and spatially smaller ORN inputs than with ortho, yet whether this is consequential for OB activity remains unknown. Indeed with these specifications for ORN inputs, our OB model captures the salient trends in our OB data. We also analyze how first and second order ORN input statistics dynamically transfer to MC spiking statistics with a phenomenological linear-nonlinear filter model, and find that retro inputs result in larger linear filters than ortho inputs. Finally, our models show that the temporal profile of ORN is crucial for capturing our data and is thus a distinguishing feature between ortho and retro stimulation, even at the OB. Using data-driven modeling, we detail how ORN inputs result in differences in OB dynamics and MC spiking statistics. These differences may ultimately shape how ortho and retro odors are coded.

Tuning network dynamics from criticality to an asynchronous state.

According to many experimental observations, neurons in cerebral cortex tend to operate in an asynchronous regime, firing independently of each other. In contrast, many other experimental observations reveal cortical population firing dynamics that are relatively coordinated and occasionally synchronous. These discrepant observations have naturally led to competing hypotheses. A commonly hypothesized explanation of asynchronous firing is that excitatory and inhibitory synaptic inputs are precisely correlated, nearly canceling each other, sometimes referred to as 'balanced' excitation and inhibition. On the other hand, the 'criticality' hypothesis posits an explanation of the more coordinated state that also requires a certain balance of excitatory and inhibitory interactions. Both hypotheses claim the same qualitative mechanism-properly balanced excitation and inhibition. Thus, a natural question arises: how are asynchronous population dynamics and critical dynamics related, how do they differ? Here we propose an answer to this question based on investigation of a simple, network-level computational model. We show that the strength of inhibitory synapses relative to excitatory synapses can be tuned from weak to strong to generate a family of models that spans a continuum from critical dynamics to asynchronous dynamics. Our results demonstrate that the coordinated dynamics of criticality and asynchronous dynamics can be generated by the same neural system if excitatory and inhibitory synapses are tuned appropriately.

Adaptation towards scale-free dynamics improves cortical stimulus discrimination at the cost of reduced detection.

Fundamental to the function of nervous systems is the ability to reorganize to cope with changing sensory input. Although well-studied in single neurons, how such adaptive versatility manifests in the collective population dynamics and function of cerebral cortex remains unknown. Here we measured population neural activity with microelectrode arrays in turtle visual cortex while visually stimulating the retina. First, we found that, following the onset of stimulation, adaptation tunes the collective population dynamics towards a special regime with scale-free spatiotemporal activity, after an initial large-scale transient response. Concurrently, we observed an adaptive tradeoff between two important aspects of population coding-sensory detection and discrimination. As adaptation tuned the cortex toward scale-free dynamics, stimulus discrimination was enhanced, while stimulus detection was reduced. Finally, we used a network-level computational model to show that short-term synaptic depression was sufficient to mechanistically explain our experimental results. In the model, scale-free dynamics emerge only when the model operates near a special regime called criticality. Together our model and experimental results suggest unanticipated functional benefits and costs of adaptation near criticality in visual cortex.

A theoretical framework for analyzing coupled neuronal networks: Application to the olfactory system.

Determining how synaptic coupling within and between regions is modulated during sensory processing is an important topic in neuroscience. Electrophysiological recordings provide detailed information about neural spiking but have traditionally been confined to a particular region or layer of cortex. Here we develop new theoretical methods to study interactions between and within two brain regions, based on experimental measurements of spiking activity simultaneously recorded from the two regions. By systematically comparing experimentally-obtained spiking statistics to (efficiently computed) model spike rate statistics, we identify regions in model parameter space that are consistent with the experimental data. We apply our new technique to dual micro-electrode array in vivo recordings from two distinct regions: olfactory bulb (OB) and anterior piriform cortex (PC). Our analysis predicts that: i) inhibition within the afferent region (OB) has to be weaker than the inhibition within PC, ii) excitation from PC to OB is generally stronger than excitation from OB to PC, iii) excitation from PC to OB and inhibition within PC have to both be relatively strong compared to presynaptic inputs from OB. These predictions are validated in a spiking neural network model of the OB-PC pathway that satisfies the many constraints from our experimental data. We find when the derived relationships are violated, the spiking statistics no longer satisfy the constraints from the data. In principle this modeling framework can be adapted to other systems and be used to investigate relationships between other neural attributes besides network connection strengths. Thus, this work can serve as a guide to further investigations into the relationships of various neural attributes within and across different regions during sensory processing.

Robust entropy requires strong and balanced excitatory and inhibitory synapses.

It is widely appreciated that balanced excitation and inhibition are necessary for proper function in neural networks. However, in principle, balance could be achieved by many possible configurations of excitatory and inhibitory synaptic strengths and relative numbers of excitatory and inhibitory neurons. For instance, a given level of excitation could be balanced by either numerous inhibitory neurons with weak synapses or a few inhibitory neurons with strong synapses. Among the continuum of different but balanced configurations, why should any particular configuration be favored? Here, we address this question in the context of the entropy of network dynamics by studying an analytically tractable network of binary neurons. We find that entropy is highest at the boundary between excitation-dominant and inhibition-dominant regimes. Entropy also varies along this boundary with a trade-off between high and robust entropy: weak synapse strengths yield high network entropy which is fragile to parameter variations, while strong synapse strengths yield a lower, but more robust, network entropy. In the case where inhibitory and excitatory synapses are constrained to have similar strength, we find that a small, but non-zero fraction of inhibitory neurons, like that seen in mammalian cortex, results in robust and relatively high entropy.

Cortical Entropy, Mutual Information and Scale-Free Dynamics in Waking Mice.

Some neural circuits operate with simple dynamics characterized by one or a few well-defined spatiotemporal scales (e.g. central pattern generators). In contrast, cortical neuronal networks often exhibit richer activity patterns in which all spatiotemporal scales are represented. Such "scale-free" cortical dynamics manifest as cascades of activity with cascade sizes that are distributed according to a power-law. Theory and in vitro experiments suggest that information transmission among cortical circuits is optimized by scale-free dynamics. In vivo tests of this hypothesis have been limited by experimental techniques with insufficient spatial coverage and resolution, i.e., restricted access to a wide range of scales. We overcame these limitations by using genetically encoded voltage imaging to track neural activity in layer 2/3 pyramidal cells across the cortex in mice. As mice recovered from anesthesia, we observed three changes: (a) cortical information capacity increased, (b) information transmission among cortical regions increased and (c) neural activity became scale-free. Our results demonstrate that both information capacity and information transmission are maximized in the awake state in cortical regions with scale-free network dynamics.

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.

Maximizing Sensory Dynamic Range by Tuning the Cortical State to Criticality.

Modulation of interactions among neurons can manifest as dramatic changes in the state of population dynamics in cerebral cortex. How such transitions in cortical state impact the information processing performed by cortical circuits is not clear. Here we performed experiments and computational modeling to determine how somatosensory dynamic range depends on cortical state. We used microelectrode arrays to record ongoing and whisker stimulus-evoked population spiking activity in somatosensory cortex of urethane anesthetized rats. We observed a continuum of different cortical states; at one extreme population activity exhibited small scale variability and was weakly correlated, the other extreme had large scale fluctuations and strong correlations. In experiments, shifts along the continuum often occurred naturally, without direct manipulation. In addition, in both the experiment and the model we directly tuned the cortical state by manipulating inhibitory synaptic interactions. Our principal finding was that somatosensory dynamic range was maximized in a specific cortical state, called criticality, near the tipping point midway between the ends of the continuum. The optimal cortical state was uniquely characterized by scale-free ongoing population dynamics and moderate correlations, in line with theoretical predictions about criticality. However, to reproduce our experimental findings, we found that existing theory required modifications which account for activity-dependent depression. In conclusion, our experiments indicate that in vivo sensory dynamic range is maximized near criticality and our model revealed an unanticipated role for activity-dependent depression in this basic principle of cortical function.

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.

Inhibition causes ceaseless dynamics in networks of excitable nodes.

The collective dynamics of a network of excitable nodes changes dramatically when inhibitory nodes are introduced. We consider inhibitory nodes which may be activated just like excitatory nodes but, upon activating, decrease the probability of activation of network neighbors. We show that, although the direct effect of inhibitory nodes is to decrease activity, the collective dynamics becomes self-sustaining. We explain this counterintuitive result by defining and analyzing a "branching function" which may be thought of as an activity-dependent branching ratio. The shape of the branching function implies that, for a range of global coupling parameters, dynamics are self-sustaining. Within the self-sustaining region of parameter space lies a critical line along which dynamics take the form of avalanches with universal scaling of size and duration, embedded in a ceaseless time series of activity. Our analyses, confirmed by numerical simulation, suggest that inhibition may play a counterintuitive role in excitable networks.

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

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

Sensory input to cortex encoded on low-dimensional periphery-correlated subspaces.

As information about the world is conveyed from the sensory periphery to central neural circuits, it mixes with complex ongoing cortical activity. How do neural populations keep track of sensory signals, separating them from noisy ongoing activity? Here, we show that sensory signals are encoded more reliably in certain low-dimensional subspaces. These coding subspaces are defined by correlations between neural activity in the primary sensory cortex and upstream sensory brain regions; the most correlated dimensions were best for decoding. We analytically show that these correlation-based coding subspaces improve, reaching optimal limits (without an ideal observer), as noise correlations between cortex and upstream regions are reduced. We show that this principle generalizes across diverse sensory stimuli in the olfactory system and the visual system of awake mice. Our results demonstrate an algorithm the cortex may use to multiplex different functions, processing sensory input in low-dimensional subspaces separate from other ongoing functions.

Maximal variability of phase synchrony in cortical networks with neuronal avalanches.

Ongoing interactions among cortical neurons often manifest as network-level synchrony. Understanding the spatiotemporal dynamics of such spontaneous synchrony is important because it may (1) influence network response to input, (2) shape activity-dependent microcircuit structure, and (3) reveal fundamental network properties, such as an imbalance of excitation (E) and inhibition (I). Here we delineate the spatiotemporal character of spontaneous synchrony in rat cortex slice cultures and a computational model over a range of different E-I conditions including disfacilitated (antagonized AMPA, NMDA receptors), unperturbed, and disinhibited (antagonized GABA(A) receptors). Local field potential was recorded with multielectrode arrays during spontaneous burst activity. Synchrony among neuronal groups was quantified based on phase-locking among recording sites. As network excitability was increased from low to high, we discovered three phenomena at an intermediate excitability level: (1) onset of synchrony, (2) maximized variability of synchrony, and (3) neuronal avalanches. Our computational model predicted that these three features occur when the network operates near a unique balanced E-I condition called "criticality." These results were invariant to changes in the measurement spatial extent, spatial resolution, and frequency bands. Our findings indicate that moderate average synchrony, which is required to avoid pathology, occurs over a limited range of E-I conditions and emerges together with maximally variable synchrony. If variable synchrony is detrimental to cortical function, this is a cost paid for moderate average synchrony. However, if variable synchrony is beneficial, then by operating near criticality the cortex may doubly benefit from moderate mean and maximized variability of synchrony.

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

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

As we go about our days, how often do we fidget, compared to how frequently we make larger movements, like walking down the hall? And how rare is a trek across town compared to that same walk down the hall? Animals tend to follow a mathematical law that relates the size of our movements to how often we do them. This law posits that small-to-medium movements and large-to-huge movements are related in the same way, that is, the law is 'scale-free', it holds the same for different scales of movement. Surprisingly, measurements of brain activity also follow this scale-free law: the level of activation of a group of neurons relates to how often they are activated in the same way for different levels of activation. Although body movements and brain activity behave in a mathematically similar way, these two facts had not previously been linked. Jones et al. studied body movements and brain activity in mice, and found that scale-free body movements were linked to scale-free brain activity, but only in certain subsets of neurons. This observation had been hidden because other subsets of neurons compete with scale-free neurons. When the scale-free neurons turn on, the competing groups turn off. When averaged together, these fluctuations cancel out. The findings of Jones et al. provide a new understanding of how brain and body dynamics are orchestrated in healthy organisms. In particular, their results suggest that the complex, multi-scale nature of behavior and body movements may emerge from brain activity operating at a critical tipping point between order and disorder, at the edge of chaos.

Low dimensional criticality embedded in high dimensional awake brain dynamics.

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

Information capacity and transmission are maximized in balanced cortical networks with neuronal avalanches.

The repertoire of neural activity patterns that a cortical network can produce constrains the ability of the network to transfer and process information. Here, we measured activity patterns obtained from multisite local field potential recordings in cortex cultures, urethane-anesthetized rats, and awake macaque monkeys. First, we quantified the information capacity of the pattern repertoire of ongoing and stimulus-evoked activity using Shannon entropy. Next, we quantified the efficacy of information transmission between stimulus and response using mutual information. By systematically changing the ratio of excitation/inhibition (E/I) in vitro and in a network model, we discovered that both information capacity and information transmission are maximized at a particular intermediate E/I, at which ongoing activity emerges as neuronal avalanches. Next, we used our in vitro and model results to correctly predict in vivo information capacity and interactions between neuronal groups during ongoing activity. Close agreement between our experiments and model suggest that neuronal avalanches and peak information capacity arise because of criticality and are general properties of cortical networks with balanced E/I.

Predicting criticality and dynamic range in complex networks: effects of topology.

The collective dynamics of a network of coupled excitable systems in response to an external stimulus depends on the topology of the connections in the network. Here we develop a general theoretical approach to study the effects of network topology on dynamic range, which quantifies the range of stimulus intensities resulting in distinguishable network responses. We find that the largest eigenvalue of the weighted network adjacency matrix governs the network dynamic range. When the largest eigenvalue is exactly one, the system is in a critical state and its dynamic range is maximized. Further, we examine higher order behavior of the steady state system, which predicts that networks with more homogeneous degree distributions should have higher dynamic range. Our analysis, confirmed by numerical simulations, generalizes previous studies in terms of the largest eigenvalue of the adjacency matrix.