Stability of infraslow correlation structure in time-shifted intracranial EEG signals

It is increasingly understood that the epilepsies are characterized by network pathology that can span multiple spatial and temporal scales. Recent work indicates that infraslow (<0.2 Hz) envelope correlations may form a basis for distant spatial coupling in the brain. We speculated that infraslow correlation structure may be preserved even with some time lag between signals. To this end, we studied intracranial EEG (icEEG) data collected from 22 medically refractory epilepsy patients. For each patient, we selected hour-long background, awake icEEG epochs before and after antiseizure medication (ASM) taper. For each epoch, we selected 5,000 random electrode contact pairs and estimated magnitude-squared coherence (MSC) below 0.15 Hz of band power time-series in the traditional EEG frequency bands. Using these same contact pairs, we shifted one signal of the pair by random durations in 15-s increments between 0 and 300 s. We aggregated these data across all patients to determine how infraslow MSC varies with duration of lag. We further examined the effect of ASM taper on infraslow correlation structure. We also used surrogate data to empirically characterize MSC estimator and to set optimal parameters for estimation specifically for the study of infraslow activity. Our empirical analysis of the MSC estimator showed that hour-long segments with MSC computed using 3-min windows with 50% overlap was sufficient to capture infraslow envelope correlations while minimizing estimator bias and variance. The mean MSC decreased monotonically with increasing time lag until 105 s of lag, then plateaued between 106 and 300 s. Significantly nonzero infraslow envelope MSC was preserved in all frequency bands until about 1 min of time lag, both pre- and post-ASM taper. We also saw a slight, but significant increase in infraslow MSC post-ASM taper, consistent with prior work. These results provide evidence for the feasibility of examining infraslow activity via its modulation of higher-frequency activity in the absence of DC-coupled recordings. The use of surrogate data also provides a general methodology for benchmarking measures used in network neuroscience studies. Finally, our study points to the clinical relevance of infraslow activity in assessing seizure risk.