Neuronal Avalanches in Neocortical Circuits

Neuronal avalanches in acute cortex slices. A, Light microscopic picture of the acute coronal brain slice and position of the microelectrode array. Cx, Cortex. B, Overplot of three different periods in which spontaneous, synchronized LFPs are visible. Three different periods (1-3) are shown (1 sec each) that reveal the occurrence of neuronal avalanches in three different, partially overlapping locations. For spatial location in the slice, see A. C, Raster plot of LFP activity in response to bath application of the NMDA receptor agonist NMDA and the D₁ receptor agonist SKF-38393. Note that synchronized events are visible across the array for ∼2000 sec after which the activity subsides. D, Individual event size distribution at Δt = IEI_avg from each acute slice experiment over plotted (n = 9 acute slices). Black, Number of electrodes; gray, summed LFPs; broken line, power law with α = -3/2). E, Average event size distribution from data shown in D. Inset, Lifetime distributions of avalanches display a power law in initial portion with characteristic slope of -2 and exponential cutoff. Broken line, Power law with exponent of α = -2.

Characteristic exponent for neuronal avalanche sizes is -3/2. A, IED versus IEI_avg for original and rescaled grid sizes. Red, Average; black, individual cultures. B, Power laws at Δt = IEI_avg for each culture have characteristic exponent α ≈ -1.5. Black, Number of electrodes; blue, LFP; average for all cultures. C, Average slopes for cultures (left) and acute slices (right). D, At Δt = IEI_avg and corresponding IED, the slope α is independent of array size. Icons indicate resampled arrays at IED = 200, 400, and 600 μm. E, Resampled power laws for summed LFP values (same arrays as in D). F, Cutoff point of the power law is determined by the number of electrodes in the array (n = 15, 30, 60; IED = 200 μm). G, Reduction in inhibition in the presence of the GABA_A receptor antagonist picrotoxin destroys the power law and renders the event size distribution bimodal. Note the presence of a large hump at higher values, indicating epileptic discharge. H, The initial slope of the event size distribution is significantly steeper (p < 0.05) in the presence of picrotoxin. Same color code as in G. I, Average event size distribution for refractory period set to 0 msec at Δt = 4 msec (three cultures). Broken line in red indicates slope of -3/2.

Size distributions for avalanches follow power laws independently of bin width Δt. A, Probability distribution of avalanche sizes (number of electrodes activated) in log-log coordinates at different Δt (average for n = 7 cultures). The linear part of each function indicates power law. Cutoff given by maximal number of electrodes (n = ∼60). Inset, Dependence of slope α on Δt: α(Δt) ∼ Δt ^{-0.16 ± 0.01} (R² = 0.99 ± 0.01; averages for all cultures). Circles, Electrodes; squares, LFP. B, Probability distribution of avalanche size distributions based on summed LFPs as a function of bin width Δt. Inset, Single culture; overplot of power laws for all seven cultures at Δt = 1 msec expressed in multiples of average LFP size.

Activity within synchronized periods is composed of avalanches. A, Raster of spontaneous activity (top) shows correlated periods containing spatiotemporal patterns (middle) and an avalanche of three frames in the original coordinates of the multielectrode array (bottom). Avalanches were defined as sequences of continuous activity that were preceded and terminated by a bin width of Δt with no activity. Dots, LFP times (dot sizes proportional to LFP amplitudes); 1 and 8, columns and rows on the grid. B, Cultured cortical networks produce thousands of avalanches of different durations per hour.

Lifetime distributions of avalanches display a power law in initial portion with characteristic slope of -2 and exponential cutoff. A, Lifetime distributions as a function of bin width Δt. B, Normalized time t/Δt collapsed from inset giving scale-invariant lifetime distributions (average for Δt = 1, 2, 4, 8, and 16 msec for n = 7 cultures). Broken line indicates slope of -2.

Network dynamics in cultured networks are characterized by a critical branching parameter of σ = 1, suggesting a state of optimal information transmission. A, Estimate of branching parameter σ from individual avalanches. σ = the ratio of descendant electrodes to ancestor electrodes. B, Sketch depicting the critical behavior of a branching process over time. If σ > 1, the size of the avalanche will grow over time, taking over the network (epilepsy), whereas at σ < 1, the avalanche will diminish quickly in size. Only at σ = 1 (critical) can avalanches persist at all scales. C, Values of σ for individual cultures (circles) based on single ancestor (left) or multiple ancestor (right) calculations. Boxes, Means ± SD. D, Phase plot of (σ, α) as a function of Δt. Note that the trajectory passes through critical point (1,-1.5) at the average population IEI of 4.2 msec. σ calculated from single ancestor avalanches. Circles, LFPs; squares, number of electrodes; ± SE.