Multiscale Aspects of Generation of High-Gamma Activity during Seizures in Human Neocortex

Methodology. A, Acute slice recording setup with concurrent extracellular and intracellular recordings. SLA and PDSs were induced with the addition of bicuculline and NMDA (black arrow). B, Local seizure activity from patients was recorded with MEAs placed concurrently with cortical grids that measured the ECoG. The power of the HG activity was determined in both individual microelectrode signals (magenta) and averaged MEA activity (blue).

Jitter analysis and synchrony. Cartoon of the jitter analysis in which a randomized delay (shift) was applied to each signal prior to computing the signal average. The jitter simulates levels of asynchrony between the signals, while the average represents the effect of volume conduction (Eq. 1). A, Flowchart of jitter methodology. Each signal within the ensemble is jittered by applying a randomized delay (τ_i) ranging from 0 to τ_max. The jittered signals are then averaged to represent the compound activity of the signals due to volume conduction. We perform power spectral analysis on the averaged activity and determine the amount of HG power in the signal. The amount of HG power is then plotted against the maximum delay (τ_max; i.e., the degree of asynchrony). This procedure is repeated for a range of increasing maximum delays (in this example, a range from 0 to 8 ms in steps of 2 ms). B, Set of randomly jittered signals (blue dots mark the center of each signal) drawn from a uniform distribution characterized by a maximum delay, τ_max. The red trace is the time average of the jittered dataset. C, The power spectra of the individual signals in the inset of A, and the average of these power spectra (in red). D, The relationship between the power in a time-averaged signal and the degree of asynchrony determined by the maximum delay (τ_max) of the applied jitter. The two gray panels depict example datasets: one in which the signals are perfectly synchronized (No jitter, i.e., Max Delay = 0) and one with significant jitter (maximum delay = 4 ms). To quantify the amount of synchrony in a set of signals, we computed the synchrony ratio by dividing the power from the time-averaged power spectrum by the power from the averaged power spectra (i.e., the values marked as 2 in D divided by those marked by 1 in C).

Seizure-like activity produces significant amounts of HG power. A, Example trace of intracellular activities recorded in an acute slice of human neocortical tissue. The arrow marks bath application of bicuculline and NMDA. Insets show examples of PDS and baseline (Base) signals. Intracellular PDSs are characterized by large depolarizations of ∼30 mV above the resting potential that includes depolarization block. This resulted in reduced interburst spike rates, compared with physiological bursts, and decreased spike amplitudes during the burst. Power spectra of intracellular signals show increased power in the HG band (80-150 Hz, gray box) during PDS activity (black trace) compared with baseline (gray trace). B, Same as in A for the extracellular signal recorded simultaneously with the intracellular signal. Insets depict SLA and network baseline activity (Base). Power spectra show an increase in the HG band power during SLA compared with baseline. C, Intracellular, nonsaturating bursting activity from acute human neocortical slice. Baseline (Base) and bursting (Burst) activities are shown in the insets. Power spectra show increased HG power during bursting over baseline. D, Extracellular bursting activity simultaneously recorded in the same slice as C. Power spectra do not show a significant increase of HG power during the network burst. E, Bar graph of intracellular HG power across states, including the mean HG power and SEM, shows a large increase in HG power during both PDSs and nonsaturating bursts. Note the discontinuity in the vertical axis. Dark blue, PDS; light blue, corresponding baseline activity; red, nonsaturating bursts; pink, corresponding baseline activity. F, Bar graph of extracellular HG power across states shows that extracellular recordings from networks with SLA (that included single neurons with PDSs) show significantly more HG power than network recordings from nonsaturating cellular bursts. Dark blue, SLA; light blue, corresponding baseline activity; red, network bursting of networks with nonsaturating cellular bursts; pink, corresponding baseline activity. **p < 0.02.

Power of average HG activity depends on synchrony levels. A, B, Sequential burst events from intracellular (IC) PDS activity (A) and extracellular (EC) SLA slice recordings (B) were detected and used to represent signals generated across multiple locations. Detected burst events (synchronous case, black traces) or randomly selected epochs (asynchronous case, green traces) were then superimposed to mimic volume conduction. These superimposed signals were then averaged, and power spectra were computed (bottom plots; black, synchronous case; green, asynchronous case; gray boxes, HG band 80-150 Hz). C, D, Bar graphs of the mean HG power and SEM of the averaged signals (gray, baseline; green, asynchronous average; black, synchronous average). Result shows a large difference in the amount of HG power (note the discontinuity in the vertical axes) between the synchronous and asynchronous scenarios (**p < 0.02). E, F, Jitter plots showing the amount of HG power retained in a series of averages in which the superimposed signals undergo randomized delays prior to averaging. For each average, these randomized delays are drawn from a uniform distribution ranging from zero to Max Delay (abscissa). Insets detail the drop in HG power retained with maximum delays of <50 ms. Both IC and EC recordings show a dramatic drop in HG power with maximum delays of <10 ms.

Microelectrodes in the core record higher levels of HG power compared with those in the penumbra. Seizure activity from MEA recordings was separated post hoc into the following two distinct territories: (1) the ictal core, defined as the cortical area where low-frequency activity is correlated with high levels of multiunit activity (two patients, four seizures); and (2) the ictal penumbra, where the low-frequency fluctuations are accompanied by uncorrelated and small changes in firing (two patients, four seizures). A, Top, Example microelectrode recording from core territories during a seizure. Bottom, Heat map of HG power detected in 100 ms overlapping epochs across the seizure (50 ms overlap). Channel # corresponds to the single microelectrode channels of the MEA. B, Same as in A for a recording taken from the penumbral territory. C, Examples of power spectra from microelectrodes in the core (orange) and penumbra (yellow) during fully developed seizure. Note that the core shows more power in the HG band. D, Bar graph of the mean and SEM of the HG power from the entire seizure for microelectrode recordings in the ictal core compared with those in the penumbra shows significantly more HG power in the core territories; note the discontinuity in the vertical axis (**p < 0.02).

Averaged MEA activity shows HG power predominantly in the core. A, Example ECoG trace recorded from a cortical grid electrode 1 cm from the MEA location (black), Pseudo-ECoG (orange) produced by averaging the microelectrode channels for a single seizure in the ictal core. Bottom, A heat map of HG power of the Pseudo-ECoG in the ictal core territory. To compute the heat map, we subdivided signals into 100 ms overlapping epochs (50 ms overlap). B, Same as in A for a recording in the penumbra. C, Bar graph of mean and SEM of the HG power of the Pseudo-ECoG activity across the entire seizure in core and penumbra shows more HG power in the averaged activity of the core; note the discontinuity in the vertical axis. D, Since computing the averaged temporal activity takes into account differences in phase between HG activities at each microelectrode, synchrony was measured as a percentage of the HG power retained after temporal averaging compared with the averaged spectra of the individual microelectrodes (Fig. 2; Eq. 4). More power was retained in the core, suggesting a higher level of HG synchrony between these microelectrodes (*p < 0.05).

Synchrony is increased during seizure activity in the core. P, Patient number; S, seizure number. A, Pseudo-ECoG recorded from ictal core territories. In the bar labeled Analysis 1, HG power is determined from the power spectrum of the time average. In contrast, in the Analysis 2 case, HG power is determined in the average spectrum of the individual power spectra of the MEA signals. The bar labeled as Synchrony Ratio represents the ratio between the two analyses (Fig. 2; Eq. 4). Heat maps and synchrony were determined in 100 ms overlapping epochs (overlap, 50 ms). The heat maps represent the averaged spectral activity, and the averaged temporal activity and their synchrony ratio. Right-hand graphs show the result of our jitter analysis, depicting the amount of HG power (ordinate) as a function of the maximum delay (abscissa). Jitter plots from core seizures show a slow decline in HG power. B, Same as in A for seizures recorded from penumbral territories. Note different scales used for the ordinate across the jitter plots.