D.M. Groppe 1,2,3 , S. Bickel 4 , C. Keller 1,2,3,4 , L. Entz 5,6 , A.D. Mehta 1,2,3 1 Departments of Neurology and Neurosurgery, North Shore Long Island Jewish Health System; 2 Feinstein Institute for Medical Research, Hofstra University School of Medicine; 3 Harvey Cushing Institutes of Neuroscience; 4 Albert Einstein College of Medicine; 5 Institute for Psychology of the Hungarian Academy of Sciences; 6 Hungarian National Institute of Neuroscience Questions •What are the dominant frequencies of oscillatory resting activity of different human neocortical areas as measured by the electrocorticogram (ECoG)? •How reliable are such dominant frequencies within an individual? Participants & Procedure •11 subjects with pharmacologically intractable epilepsy (5 female) •Mean subject age: 32 (SD=14) years •Two clips of ECoG recordings were taken from 24 hour clinical recordings of each patient based on video footage of the subjects •During each clip, patients were awake and resting, passively watching TV, or potentially listening to the conversations of others •Mean clip duration: 6.2 (SD=1.7) minutes •Median time between clips: 33.8 (SIQR=25.3) hours • ~7 Hz oscillations, are generally the dominant mode of sub-13 Hz activity of the cortical areas sampled. This contradicts previous estimates from scalp EEG the scalp which generally found predominant ~10 Hz activity 1,2 and muddies the conventional distinction between theta and alpha bands. • Theta peaks most reliably found in fusiform, inferior-middle temporal, parahippocampal, lateral orbitofrontal cortex. • Beta peaks found most reliably in middle temporal, parsopercularis, parstriangularis, pre & post central cortex. It is not clear to what extent these peaks are harmonics of lower frequencies • Gamma peaks found in entorhinal cortex. May be related to pathology • Peaks reproduce across days and are evident in .5 min of data. But reproducibility might be improved by better control of subject state and environment. Dominant frequencies of resting human neocortical oscillations as measured by the electrocorticogram Introduction Methods Results Conclusions Background Neocortical electromagnetic activity often exhibits oscillations at different frequencies. These oscillations are apparent in the activity of the resting brain and are modulated by sensory stimulation, behavior, and cognitive state. Previous research on humans has found that the most dominant frequencies of these oscillations tend varies across neocortical areas. For example, using EEG recordings of the resting brain, Nunez et al. have found alpha activity (8-12 Hz) to be widespread and theta activity (4-7 Hz) to be somewhat frontally distributed 1,2 . Such systematic differences in neural dynamics likely reflect differences in connectivity between brain areas and resting function and may be a useful biomarker of neural disorders 3 . However characterizing the dominant frequencies of neocortical areas in humans has been difficult due the poor spatial resolution of noninvasive measures of brain electromagnetic activity. Our Study We measured magnitude of neocortical oscillatory activity in the resting brain using ECoG recordings obtained via intracranial electrodes implanted on the surface of the brain in patients with pharmacologically intractable epilepsy undergoing evaluation for surgery. ECoG Recording & Preprocessing •On average 109 (SD=10) subdural electrodes were implanted in each subject •ECoG recorded at 500-2000 Hz using an XLtek clinical system •ECoG referenced to electrodes over interhemispheric fissure and re-referenced to average reference after electrodes over ictal onset regions or with poor contacts were rejected (leaving a mean of 88, SD=10, usable electrodes per subject) •Time periods in which the patient moved or that contained interictal spikes were manually identified and rejected, leaving a mean of 5.6 (SD=1.1) minutes of data for analysis •ECoG whitened by taking its temporal derivative to heighten spectral power peaks 1. P. L. Nunez, B. M. Wingeier, R. B. Silberstein, Spatial-temporal structures of human alpha rhythms: theory, microcurrent sources, multiscale measurements, and global binding of local networks, Hum. Brain Mapp. 13, 125 (2001). 2. P. L. Nunez, L. Reid, R. G. Bickford, The relationship of head size to alpha frequency with implications to a brain wave model, Electroencephalogr. Clin. Neurophysiol. 44, 344 (1978). 3. D. Zhang, M. E. Raichle, Disease and the brain's dark energy, Nat. Rev. Neurol. 6, 15 (2010). References A B C D E Spectral Analysis •Segments of data divided into 1 sec epochs with 0.5 sec overlap •Power spectrum density (PSD) for each epoch estimated with the mean of two Slepian tapers and discrete Fourier transform •Mean PSD for each electrode derived by taking the 5% trimmed mean of all epochs to ignore outliers •PSD normalized to unit area ignoring line noise frequencies and frequencies above 164 Hz MRI Coregistration and Cortical Parcellation •(A) Postimplant CT scan, (B) Manual identification of electrodes, (C) Coregistration of CT scan to preimplant MRI, (D) Coregistration of electrode locations to FreeSurfer pial surface, (E) Electrode locations snapped to pial surface •Cortical areas automatically defined by FreeSurfer using the Desikan-Killiany atlas 1-12 Hz Peak Frequency Mean PSD Across Subjects per Cortical Area: Shaded regions represent 95% confidence intervals. n indicates the number of subjects contributing to each average. Dashed lines indicate peaks. Max Power in Bands: Mean frequency of maximal power value in select bands. Error bars=95% CI. Solid lines indicate frequency range boundaries. Bars that differ from boundaries indicate PSD peaks. 13-30 Hz Peak Frequency 31-127 Hz Peak Frequency Peak frequencies across all electrodes: Histogram of PSD peaks across all electrodes. Blank regions indicate frequencies ignored due to line noise. Uniform predominance of peaks above ~50 Hz reflects residual 1/f PSD distribution. 10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 160 0 0.05 0.09 0.14 0.19 0.24 0.28 Hz Relative Frequency 7 Hz 14 Hz 19 Hz 23 Hz ! " #$ %& $! #&" &’$ ’#& ( ()# ()& ()% ()! ()’ ()$ ()* ()" +,-./012 14 50-6/ 7809 :6;<12=6> 7-166!7809 71--;8./012 #!#& ?@ #%!%( ?@ Contact: [email protected] Test-retest peak reliability: Mean correlation between maximal PSD frequencies in theta-alpha and beta bands across the two clips of data for each patient. Various amounts of the data in the first clip were used to estimate the amount data necessary to get reliable estimates. 0 2 4 6 8 10 12 14 superior frontal rostral middle frontal medial orbitofrontal lateral orbitofrontal parsorbitalis parstriangularis parsopercularis caudal middle frontal precentral postcentral middle temporal inferior temporal entorhinal parahippocampal fusiform precuneus inferior parietal lingual lateral occipital Hz Alpha/Theta Boundary 8 10 12 14 16 18 20 22 24 superior frontal rostral middle frontal medial orbitofrontal lateral orbitofrontal parsorbitalis parstriangularis parsopercularis caudal middle frontal precentral postcentral middle temporal inferior temporal entorhinal parahippocampal fusiform precuneus inferior parietal lingual lateral occipital Hz 20 30 40 50 60 70 80 90 superior frontal rostral middle frontal medial orbitofrontal lateral orbitofrontal parsorbitalis parstriangularis parsopercularis caudal middle frontal precentral postcentral middle temporal inferior temporal entorhinal parahippocampal fusiform precuneus inferior parietal lingual lateral occipital Hz