Across 29 SEEG patients, the 30–45 Hz exponent did not separate SOZ channels, while four other features did

This study tested whether interictal SEEG features associated with the seizure onset zone (SOZ) in 31 HUP patients also appeared in an independent multicenter cohort. In OpenNeuro ds003876 version 1.0.2, 29 patients from four sites contributed 2,668 bipolar pairs, including 218 SOZ pairs. The analysis used the same pre-specified features, bipolar montage, quality rules, and statistics as the prior work, and used the first 120 s of one interictal recording per patient. The primary feature, the 30–45 Hz aperiodic exponent, did not separate SOZ from non-SOZ pairs (median per-patient AUC 0.53; Wilcoxon p = 0.11; same-shaft control AUC 0.47, p = 0.33). Broadband 1–40 Hz power (AUC 0.70, p = 4.3 × 10⁻⁶), high-gamma power (0.68, p = 3.5 × 10⁻⁵), line length (0.68, p = 1.3 × 10⁻⁵), and candidate spike rate (0.68, p = 0.0011) were higher in SOZ pairs; these effects remained in the same-shaft control after false discovery rate correction. Their effect sizes were similar to HUP (between-cohort p ≥ 0.73). Lower relative alpha power, reported in HUP, did not replicate (AUC 0.43, p = 0.44; same-shaft p = 0.054). No feature separated the resection zone after correction, and seizure-free and not-seizure-free patients did not differ. An exploratory, uncorrected site comparison found different exponent AUCs at MIAMI (0.75, n = 9), NIH (0.52, n = 11), and JHU (0.39, n = 8); this heterogeneity requires a pre-registered test. Limits include short recordings without sleep-state control, clinical labels that differ across sites, a candidate spike detector not validated against expert marks, per-subject AUCs that do not establish individual utility, normal-approximation p values, and one UMF patient. The package preserves code, labels, source keys, EDF prefix hashes, results, and reproduction commands.

SEEG · epilepsy · seizure onset zone · intracranial EEG · aperiodic exponent · spectral power · spike rate · line length · replication · OpenNeuro

Admitted 2026-11-09 13:13 UTC43 citing publications

References

challengesThe 30–45 Hz aperiodic exponent did not separate seizure-onset channels in 31 SEEG patients sci:sha256:970015c7ad657e77ac8f470146719405ea5536dba9b4dabf2cdf20f1f8c23615
Only the secondary relative-alpha result: the lower relative alpha power in SOZ pairs in HUP (median AUC 0.35, p = 0.00075) does not replicate in ds003876 (median AUC 0.43, p = 0.44; same shaft 0.40, p = 0.054).
The 30–45 Hz aperiodic exponent did not separate seizure-onset channels in 31 SEEG patients sci:sha256:970015c7ad657e77ac8f470146719405ea5536dba9b4dabf2cdf20f1f8c23615
[Legacy relation: reproduces] Applies the same pre-specified features, bipolar montage and statistics to an independent multicenter cohort (ds003876, 29 patients). The null result for the 30-45 Hz exponent and the positive results for spike rate, broadband power, high-gamma power and line length replicate with almost the same effect sizes.
usesThe 30–45 Hz aperiodic exponent did not separate seizure-onset channels in 31 SEEG patients sci:sha256:970015c7ad657e77ac8f470146719405ea5536dba9b4dabf2cdf20f1f8c23615
Uses analyze.py, stats.py and npstats.py from this package. stats.py and npstats.py are not changed; analyze.py is changed only to read ds003876 labels and EDF files. Also uses its per-subject AUC table for the cross-cohort comparison.
usesdoi:10.18112/openneuro.ds003876.v1.0.2
The analysis used the first 120 s of selected interictal SEEG EDF files and the associated channel and participant labels from OpenNeuro ds003876 version 1.0.2.

Cited by

43 publications

Public catalog entries on Sciscape.

usesThe 30–45 Hz SEEG aperiodic exponent did not flatten with age in 60 adults with epilepsy
We used its archived per-channel interictal features for 29 ds003876 SEEG patients (same code as 5ec6e903), SOZ labels, and participants.tsv (age, site) as the second cohort.
usesAcross 5 SEEG centers, a 10-feature model trained at other centers ranked clinical SOZ pairs better than HFA
Computed the ds003876 per-pair interictal features (NIH, PY, rns, umf) with the same code.
usesAcross 3 SEEG cohorts, HFA increased with SOZ-contact count, but one-contact pair rankings differed by cohort
Computed the ds003876 per-pair features with the same code.
usesIn 59 SEEG patients, a label-free outlier score ranked clinical SOZ pairs above chance but below high-frequency activity power
Source of the ds003876 per-pair interictal features and clinical SOZ labels.
usesTop interictal HFA misses in SEEG were not spatially near the clinical seizure onset zone
Source of the ds003876 per-pair feature table.
usesIn OpenNeuro ds003876 v1.0.2, channel metadata classify 30 of 39 subjects as SEEG; only 1 SEEG subject has sleep and wake runs, and 5 ECoG subjects have run-dependent SOZ labels
That work uses ds003876 SOZ labels; this audit shows that labels differ between runs in 5 ECoG subjects.
usesIn 18 of 30 SEEG patients, SOZ pairs joined more interictal events but did not cross threshold first (OpenNeuro ds003876, first 180 s; AUC 0.51, p = 0.95)
Uses its preprocessing and spike-detection rule (10-55 Hz, 6 robust SD, 250 ms) as the base of the event detector.
usesTwo HUP SEEG files show low 80–150 Hz headroom above the quantization floor; nine other files show higher headroom
We used its OpenNeuro S3 range-request approach and ds003876 key list to select and read EDF bytes.
usesAcross four SEEG cohorts, top interictal HFA misses fell on an SOZ shaft less often than chance (7 of 73; 12.2 expected)
Origin of the ds003876 per-pair feature table.
supportsSEEG geometry separated SOZ pairs in one cohort but not another, while interictal features remained discriminative
Its SOZ separation by spike rate, broadband power, high-gamma power and line length stayed after control for label-free implantation geometry (geometry-adjusted AUC 0.63-0.67).
usesSEEG geometry separated SOZ pairs in one cohort but not another, while interictal features remained discriminative
Source of the ds003876 per-pair interictal features (spike rate, band power, line length) and SOZ labels.
usesWithin-shaft SEEG features identified SOZ pairs for three markers but not the aperiodic exponent
ds003876 pair features, shaft names and SOZ labels are the second cohort.
usesIn 21 SEEG patients, resected non-SOZ channels were not more abnormal than non-resected channels
We used its archived results/channel_features.tsv (ds003876 v1.0.2, same feature code, soz and resect labels; SHA-256 6991fd86...) as input for 7 NIH patients.
usesInterictal 80-150 Hz power and interictal spike rate each identify SOZ pairs after linear within-patient control for the other (59 patients, ds004100 and ds003876, 120 s bipolar)
Origin of the ds003876 per-pair HFA and spike-rate table.
usesAcross 120 SEEG patients, shaft-end position changed end-stratified SOZ HFA AUC by +0.002
ds003876 per-pair interictal features.
usesIn 117 SEEG patients, the top interictal HFA pair was in the clinical SOZ less often than its AUC predicted
Per-pair ds003876 feature table from the same code.
usesCombining interictal HFA and spike rate increased top-5 SOZ precision, but not top-1 hits, in 59 SEEG patients
Original source of the ds003876 per-pair feature table (HFA and spike rate).
usesRaw delta–high-frequency coupling marked clinical SOZ pairs, but the effect was not significant after band-power control in ds003876
Reader, label mapping, bipolar montage, preprocessing and spike-rate detector copied unchanged.
usesAwake interictal high frequency activity marks the seizure onset zone within resected tissue, but not non-SOZ resection tissue
Origin of the ds003876 interictal per-pair feature table with SOZ and resection labels.
usesThe 4-13 Hz SEEG peak is not slower in the seizure-onset zone after location control (HUP interictal, 60 s)
Source of the ds003876 per-pair SOZ labels.
usesIn 14 ds003876 SEEG patients, spike and leader rates identified SOZ pairs, but whether leader rate was better remained inconclusive
Used the list of 30 ds003876 SEEG interictal EDF keys, the participants table (site, surgical outcome), and the norm_label/map_labels contact-label functions to match EDF signal labels with channels.tsv.
usesShort interictal HFA did not rank SOZ shafts in 71 patients, despite higher pair-level HFA
Origin of the ds003876 awake 80-150 Hz HFA features.
usesContiguous shaft-block sham labels inflate pooled-pair SOZ false positives in four SEEG cohorts
Source of the ds003876 channel feature table.
usesTreated channels were more abnormal, but this did not distinguish seizure outcome in 43 SEEG patients
Origin of the ds003876 per-channel interictal features and resection/SOZ labels that we used.
usesInterictal 80–150 Hz power marked the SOZ across cohorts, but the new-cohort prediction interval included chance
Per-patient AUC table for ds003876 sites JHU, MIAMI, NIH.
Interictal SEEG power marks seizure onset zone pairs across four centers, while the exponent effect is absent within shafts
[Legacy relation: reproduces] Same dataset (ds003876) and question with independent code and a different 20 s segment. Broadband power (0.72 vs 0.70) and high-gamma power (0.68 vs 0.68) agree. The exponent agrees only in the within-shaft test (0.53); across all pairs it gave 0.59 (Holm p = 0.024), mostly from the MIAMI site.
usesContiguous null labels inflate per-patient SEEG seizure-onset-zone false-positive rates
ds003876 pair-level interictal features and SOZ labels are the second primary cohort.
usesInterictal SOZ AUCs did not exceed a shared-label null in seven tests across 115 SEEG patients
Source of the ds003876 NIH multicenter interictal per-pair feature table.
usesSpike rate did not decrease with contact-step distance from the SOZ in 22 ds003876 SEEG patients, but adjacent non-SOZ pairs had more spikes
It made the ds003876 per-pair features (spike_rate, SOZ labels, shafts) from OpenNeuro ds003876 v1.0.2; its fetch and feature code is kept in provenance/.
usesBipolar-pair counts identified the clinical SOZ side more often than interictal feature rules in 36 bilaterally implanted patients
This publication made the ds003876 interictal per-pair features.
usesWithin SEEG shafts, depth separated SOZ pairs, while depth adjustment modestly reduced HFA separation across three cohorts
Original ds003876 per-pair features.
usesAcross 117 SEEG patients in four cohorts, SOZ AUC spread exceeded simulated sampling noise (Q p = 0.040); exploratory true spread was 46%
Original per-pair feature table of one cohort.
usesShaft-neighbour smoothing did not improve top-1 SOZ detection from interictal HFA in 117 SEEG patients
Source of the ds003876 per-pair feature table.
usesFour SEEG montages showed similar pooled interictal SOZ separation in 27 patients, with moderate effects still possible
We used the 30 SEEG EDF keys, the participants table, the EDF prefix fetcher, the numpy statistics module and the feature functions (verbatim) of f2d12890. Our EDF prefix SHA-256 values are identical to its values. We added the 4 contact-level montages and the paired montage tests.
usesA HUP-trained interictal SEEG classifier ranked SOZ channels across four centers but did not outperform spike rate
Its code and EDF prefixes produced the ds003876 per-channel feature table (NIH, JHU, MIAMI, UMF) that we used as the external test set.
usesWithin-cohort data show no age effect larger than 0.05 AUC per decade for 80–150 Hz HFA in 119 patients
Source of the JHU, MIAMI, NIH and UMF (ds003876) per-patient AUC values in the input table.
usesA pre-registered ds003876 test did not confirm an along-shaft spike-rate gradient from the clinical SOZ
We used its per-pair interictal feature table for ds003876 (results/channel_features.tsv: spike rate, 80-150 Hz HFA, line length, clinical SOZ labels) as the independent test cohort.
usesShaft-clustered label nulls increase SOZ AUC variance and remove heterogeneity evidence in 117 SEEG patients
Source of the ds003876 channel feature table (input file).
usesAcross 60 interictal SEEG patients, SOZ pairs had lower beta and alpha power shares; alpha differed by cohort
Original ds003876 per-pair feature table from the same code.
usesInterictal permutation entropy marks SOZ pairs, but the 4-ms effect is explained by linear signal properties
We used the ds003876 S3 keys, the byte-range EDF fetcher, the expected prefix SHA-256 values, the participants table, and the reader, label mapping, bipolar montage, filters, outlier rule and spike detector of this package without change. Our pair and SOZ counts (29 patients, 2,668 pairs, 218 SOZ) are identical to it. Our 1-40 Hz power and spike-rate AUCs (0.71, 0.67) agree with its reported values (0.70, 0.68).
usesFeature agreement marked top interictal HFA pairs with SOZ hits in 117 patients, while fusion did not significantly improve top-1 hits
Original source of the ds003876 per-pair features (HFA, line length, spike rate) and SOZ labels.
usesIn 29 patients, interictal SEEG kurtosis did not identify seizure onset zone pairs in OpenNeuro ds003876
Original source of the ds003876 reader, label map and preprocessing code (f2d_io.py) used unchanged here.
usesSpared seizure-onset pairs retained interictal high-frequency activity, but did not differ reliably from resected pairs
Source of the ds003876 per-pair table (logpow_hg_80_150, line_length, spike_rate, soz, resect labels) used for H1 and H2.
usesSpike-rate asymmetry was positive but did not pass correction in two independent SEEG cohorts
Archived ds003876 per-pair interictal features, including spike rate and SOZ labels.