Detecting Spontaneous Neural Oscillation Events in Primate Auditory Cortex

Stereotyped ERPs in NHP A1: supragranular, granular, infragranular layers (left to right; 50-dB clicks in NHP). Top, Average ERP waveforms (click at 0 ms). Middle, Wavelet transform spectrograms (values are in units of Power; note that these spectrograms are not normalized). Bottom, Apparent oscillation peaks (average of spectrogram over time).

Validation of event detection algorithm. A, Example of 1-cycle, 11-cycle simulated 10-Hz α signals (green bounding boxes) added to supragranular CSD; amplitude 1.5 mV/mm². B, 1-cycle signal from A was detected as two cycles. C, 11-cycle detected as 11.6 cycles. B, C, Top, Wavelet transform spectrograms with detected event in boxes; bottom: raw (red), filtered (blue) signal. D, Detected number of cycles as a function of the actual burst duration (dotted black line) for BOSC (blue) and OEvent (red). E, Peak frequency detected was generally close to the frequency of the simulated signal (horizontal gray line) using both methods in most frequency bands but varied across different burst durations (SEM: dashed lines; insets in D, E show the results for low-γ frequency band). F, G, Number of cycles and peak frequency detected in NHP A1 recordings using BOSC (blue) and OEvent (red).

Oscillation events in at least one of the bands occupy the majority of recording time. Example from NHP supragranular A1. Oscillation events (red bounding boxes) occur in one or more frequency bands during this 3-s period. A long δ event is detected from t = 0 to near the end of this 3-s period (red box appears as line across bottom). Red dots show peak frequency; height of box indicates frequency spread. Normalized wavelet transform spectrogram (top) of signal shown at bottom (spectrogram values normalized by median power at a given frequency). This example shows events in δ, θ, α, β, low γ, and γ bands.

Examples of oscillation events from NHP A1 supragranular layer. Normalized Morlet wavelet spectrograms demonstrating individual events (red bounding-box) with raw (red) and filtered (blue) waveforms below (black trace: period outside of detected oscillation). Note that x- and y-scales differ for different bands; spectrogram power (color) is in median normalized units. White text in the spectrograms specifies the events’ power relative to the median, the peak frequency of the event (and frequency range), the number of cycles, and the correlation value between the raw and filtered waveforms (filter-match).

Examples of oscillation events from human STG iEEG. Normalized Morlet wavelet spectrograms demonstrating individual events (red bounding-box) with raw (red) and filtered (blue) waveforms below [x-y-scales differ for different bands; spectrogram power (color) is in median normalized units; iEEG time-series values are in units of Volts]. Time of 0 ms corresponds to the wavelet phase of 0 radians (local maxima) closest to the event’s peak power at threshold detection. White text in the spectrograms specifies the events’ power relative to the median, the peak frequency of the event (and frequency range), the number of cycles, and the correlation value between the raw and filtered waveforms (filter-match).

Individual oscillation events from NHP A1. δ (A), θ (B), α (C), β (D), γ (E), high γ (F). In each case, vertical axis arranges waveforms with decreasing numbers of cycles from top to bottom; each row organizes waveforms with decreasing filter-match between raw and filtered signal from left to right. Horizontal scale bars 1 s except: E, 200 ms; F, 100 ms. Examples in green bounding boxes in A, F are described in the text. Note that waveforms are normalized to allow easier visual comparison.

Oscillation event rates and active time varied by frequency band. A, Higher frequency events are more frequent. B, ATR: lower and higher frequency events fill much of the recording duration (HGamma: high-γ bands; mean ± SEM in both A, B; see Extended Data Table 8-1, Table 8-2).

Event features. A, Number of cycles. B, Number of local peaks in the time domain of the filtered waveform. C, Intraevent Fspan = log(maxF/minF); 0.7 is freq doubling. D, Filter-match r value (see Extended Data Tables 9-1, 9-2, 9-3, 9-4 for comparisons).

Parameters of OEvent. Influence of different parameter combinations on the detected number of cycles for each frequency band. The results are averaged across all recordings and cortical layers. Error bars indicate SEM. Each row corresponds to a different wavelet width (mspecwidth = 5, 7, 9), and each column corresponds to a different threshold for merging overlapping events (overlap = 0.35, 0.5, 0.65). Colored lines indicate different detection thresholds. Higher detection thresholds resulted in lower detected number of cycles for all bands (t test, p < 0.01 for all bands, FDR corrected). No significant differences were observed between different wavelet widths or overlap thresholds within any of the bands.

IEIs suggest rhythmic recurrence in all bands. A, CV2 values (mean and SEM) in all bands < 1 (p < 0.05, one-sided Wilcoxon signed-rank test). B, FF (mean and SEM) in all bands except for LGamma < 1 (p < 0.05, one-sided Wilcoxon signed-rank test).

Individual oscillation events are categorized by the frequency of maximum amplitude (peak frequency; x-axis) but may spread to lower or higher frequency bands. A, B, % of oscillation events with minimum and maximum frequency in a specific band. C, % of oscillation events that are band-limited, defined as having both minF and maxF in the same band as the event’s peak. In some frequency bands (i.e., δ, γ, and high γ) individual events tend to be more band limited (C), while in other bands (i.e., θ, α, and β) as many as 60−80% of the events spread to adjacent bands (x-axis location indicates categorization of an oscillation event based on its peak frequency; color indicates which frequency band an oscillation event spreads to, e.g., in A, ∼75%, 25% of θ events have a minimum frequency in the θ, δ bands, respectively).