Repertoires of Spike Avalanches Are Modulated by Behavior and Novelty

Parabolic avalanche scaling in the synchronization of cortical cell assemblies

Neurons in the cerebral cortex fire coincident action potentials during ongoing activity and in response to sensory inputs. These synchronized cell assemblies are fundamental to cortex function, yet basic dynamical aspects of their size and duration are largely unknown. Using 2-photon imaging of neurons in the superficial cortex of awake mice, we show that synchronized cell assemblies organize as scale-invariant avalanches that quadratically grow with duration. The quadratic avalanche scaling was only found for correlated neurons, required temporal coarse-graining to compensate for spatial subsampling of the imaged cortex, and suggested cortical dynamics to be critical as demonstrated in simulations of balanced E/I-networks. The corresponding time course of an inverted parabola with exponent of χ = 2 described cortical avalanches of coincident firing for up to 5 s duration over an area of 1 mm². These parabolic avalanches maximized temporal complexity in the ongoing activity of prefrontal and somatosensory cortex and in visual responses of primary visual cortex. Our results identify a scale-invariant temporal order in the synchronization of highly diverse cortical cell assemblies in the form of parabolic avalanches.

Microscale dynamics of electrophysiological markers of epilepsy

OBJECTIVE

Interictal discharges (IIDs) and high frequency oscillations (HFOs) are established neurophysiologic biomarkers of epilepsy, while microseizures are less well studied. We used custom poly(3,4-ethylenedioxythiophene) polystyrene sulfonate (PEDOT:PSS) microelectrodes to better understand these markers' microscale spatial dynamics.

METHODS

Electrodes with spatial resolution down to 50 µm were used to record intraoperatively in 30 subjects. IIDs' degree of spread and spatiotemporal paths were generated by peak-tracking followed by clustering. Repeating HFO patterns were delineated by clustering similar time windows. Multi-unit activity (MUA) was analyzed in relation to IID and HFO timing.

RESULTS

We detected IIDs encompassing the entire array in 93% of subjects, while localized IIDs, observed across < 50% of channels, were seen in 53%. IIDs traveled along specific paths. HFOs appeared in small, repeated spatiotemporal patterns. Finally, we identified microseizure events that spanned 50-100 µm. HFOs covaried with MUA, but not with IIDs.

CONCLUSIONS

Overall, these data suggest that irritable cortex micro-domains may form part of an underlying pathologic architecture which could contribute to the seizure network.

SIGNIFICANCE

These results, supporting the possibility that epileptogenic cortex comprises a mosaic of irritable domains, suggests that microscale approaches might be an important perspective in devising novel seizure control therapies.

Selective Participation of Single Cortical Neurons in Neuronal Avalanches

Neuronal avalanches are scale-invariant neuronal population activity patterns in the cortex that emerge in vivo in the awake state and in vitro during balanced excitation and inhibition. Theory and experiments suggest that avalanches indicate a state of cortex that improves numerous aspects of information processing by allowing for the transient and selective formation of local as well as system-wide spanning neuronal groups. If avalanches are indeed involved with information processing, one might expect that single neurons would participate in avalanche patterns selectively. Alternatively, all neurons could participate proportionally to their own activity in each avalanche as would be expected for a population rate code. Distinguishing these hypotheses, however, has been difficult as robust avalanche analysis requires technically challenging measures of their intricate organization in space and time at the population level, while also recording sub- or suprathreshold activity from individual neurons with high temporal resolution. Here, we identify repeated avalanches in the ongoing local field potential (LFP) measured with high-density microelectrode arrays in the cortex of awake nonhuman primates and in acute cortex slices from young and adult rats. We studied extracellular unit firing in vivo and intracellular responses of pyramidal neurons in vitro. We found that single neurons participate selectively in specific LFP-based avalanche patterns. Furthermore, we show in vitro that manipulating the balance of excitation and inhibition abolishes this selectivity. Our results support the view that avalanches represent the selective, scale-invariant formation of neuronal groups in line with the idea of Hebbian cell assemblies underlying cortical information processing.

Dynamical phase transition in spike neuronal firing patterns of hippocampal cells

There is increasing evidence that the brain resides in a state of criticality. The purpose of the present work is to characterize the dynamics of individual hippocampal CA1 pyramidal cells and to investigate how it is influenced by changes in Kv7.2/7.3 (M-channel) ion channel modulation, which is known to be key in determining the neuronal excitability. We show that the resting activity of CA1 neurons exhibit random dynamics with low information content, while changes in M-channel modulation move the neuronal activity near a phase transition to richer non-trivial dynamics. We interpret these results as the basis upon which the state of self-organized criticality is built.

Brain songs framework used for discovering the relevant timescale of the human brain

A key unresolved problem in neuroscience is to determine the relevant timescale for understanding spatiotemporal dynamics across the whole brain. While resting state fMRI reveals networks at an ultraslow timescale (below 0.1 Hz), other neuroimaging modalities such as MEG and EEG suggest that much faster timescales may be equally or more relevant for discovering spatiotemporal structure. Here, we introduce a novel way to generate whole-brain neural dynamical activity at the millisecond scale from fMRI signals. This method allows us to study the different timescales through binning the output of the model. These timescales can then be investigated using a method (poetically named brain songs) to extract the spacetime motifs at a given timescale. Using independent measures of entropy and hierarchy to characterize the richness of the dynamical repertoire, we show that both methods find a similar optimum at a timescale of around 200 ms in resting state and in task data.

Oscillatory Patterns of Phase Cone Formations near to Epileptic Spikes Derived from 256-Channel Scalp EEG Data

Our objective was to determine if there are any distinguishable phase cone clustering patterns present near to epileptic spikes. These phase cones arise from episodic phase shifts due to the coordinated activity of cortical neurons at or near to state transitions and can be extracted from the high-density scalp EEG recordings. The phase cone clustering activities in the low gamma band (30-50 Hz) and in the ripple band (80-150 Hz) were extracted from the analytic phase after taking Hilbert transform of the 256-channel high density (dEEG) data of adult patients. We used three subjects in this study. Spatiotemporal contour plots of the unwrapped analytic phase with 1.0 ms intervals were constructed using a montage layout of 256 electrode positions. Stable phase cone patterns were selected based on the criteria that the sign of the spatial gradient did not change for at least three consecutive time samples and the frame velocity was within the range of propagation velocities of cortical axons. These plots exhibited dynamical formation of phase cones which were higher in the seizure area as compared with the nearby surrounding brain areas. Spatiotemporal oscillatory patterns were also visible during ±5 sec period from the location of the spike. These results suggest that the phase cone activity might be useful for noninvasive localization of epileptic sites and also for examining the cortical neurodynamics near to epileptic spikes.

Hysteresis, neural avalanches, and critical behavior near a first-order transition of a spiking neural network

Many experimental results, both in vivo and in vitro, support the idea that the brain cortex operates near a critical point and at the same time works as a reservoir of precise spatiotemporal patterns. However, the mechanism at the basis of these observations is still not clear. In this paper we introduce a model which combines both these features, showing that scale-free avalanches are the signature of a system posed near the spinodal line of a first-order transition, with many spatiotemporal patterns stored as dynamical metastable attractors. Specifically, we studied a network of leaky integrate-and-fire neurons whose connections are the result of the learning of multiple spatiotemporal dynamical patterns, each with a randomly chosen ordering of the neurons. We found that the network shows a first-order transition between a low-spiking-rate disordered state (down), and a high-rate state characterized by the emergence of collective activity and the replay of one of the stored patterns (up). The transition is characterized by hysteresis, or alternation of up and down states, depending on the lifetime of the metastable states. In both cases, critical features and neural avalanches are observed. Notably, critical phenomena occur at the edge of a discontinuous phase transition, as recently observed in a network of glow lamps.

The characteristic patterns of neuronal avalanches in mice under anesthesia and at rest: An investigation using constrained artificial neural networks

Local perturbations within complex dynamical systems can trigger cascade-like events that spread across significant portions of the system. Cascades of this type have been observed across a broad range of scales in the brain. Studies of these cascades, known as neuronal avalanches, usually report the statistics of large numbers of avalanches, without probing the characteristic patterns produced by the avalanches themselves. This is partly due to limitations in the extent or spatiotemporal resolution of commonly used neuroimaging techniques. In this study, we overcome these limitations by using optical voltage (genetically encoded voltage indicators) imaging. This allows us to record cortical activity in vivo across an entire cortical hemisphere, at both high spatial (~30um) and temporal (~20ms) resolution in mice that are either in an anesthetized or awake state. We then use artificial neural networks to identify the characteristic patterns created by neuronal avalanches in our data. The avalanches in the anesthetized cortex are most accurately classified by an artificial neural network architecture that simultaneously connects spatial and temporal information. This is in contrast with the awake cortex, in which avalanches are most accurately classified by an architecture that treats spatial and temporal information separately, due to the increased levels of spatiotemporal complexity. This is in keeping with reports of higher levels of spatiotemporal complexity in the awake brain coinciding with features of a dynamical system operating close to criticality.

Neocortical activity is stimulus- and scale-invariant

Mounting evidence supports the hypothesis that the cortex operates near a critical state, defined as the transition point between order (large-scale activity) and disorder (small-scale activity). This criticality is manifested by power law distribution of the size and duration of spontaneous cascades of activity, which are referred as neuronal avalanches. The existence of such neuronal avalanches has been confirmed by several studies both in vitro and in vivo, among different species and across multiple spatial scales. However, despite the prevalence of scale free activity, still very little is known concerning whether and how the scale-free nature of cortical activity is altered during external stimulation. To address this question, we performed in vivo two-photon population calcium imaging of layer 2/3 neurons in primary visual cortex of behaving mice during visual stimulation and conducted statistical analyses on the inferred spike trains. Our investigation for each mouse and condition revealed power law distributed neuronal avalanches, and irregular spiking individual neurons. Importantly, both the avalanche and the spike train properties remained largely unchanged for different stimuli, while the cross-correlation structure varied with stimuli. Our results establish that microcircuits in the visual cortex operate near the critical regime, while rearranging functional connectivity in response to varying sensory inputs.