Utilization of independent component analysis for accurate pathological ripple detection in intracranial EEG recordings recorded extra- and intra-operatively
Introduction
Approximately 30% of patients with epilepsy continue to have disabling seizures despite treatment with multiple antiepileptic drugs (Kwan and Brodie, 2000). Of those patients with focal epilepsy that are resistant to multiple medications, resective surgery is an intervention that has been proven to reduce the seizure burden, improve the patients’ quality of life, and reduce mortality rate (Wiebe et al., 2001, Sperling et al., 2016). The goal of resective epilepsy surgery is to identify and remove epileptogenic brain regions while minimizing residual neurological deficits. High frequency oscillations (HFOs), which consist of brief bursts of energy with spectral content ranging between 80 and 600 Hz, have shown significant promise as a potential biomarker of epileptogenic tissue (Engel et al., 2009, Gotman, 2010, Jacobs et al., 2012). HFOs with a spectral content in the 80–250 Hz band are commonly referred to as ripples, while those in the 250–600 Hz band are termed fast ripples (Staba et al., 2002, Bragin et al., 2002).
HFOs can be identified by visual inspection of intracranial EEG (iEEG), or using automated and unsupervised detection software (Zelmann et al., 2012, Burnos et al., 2014, Gliske et al., 2016). One barrier to utilizing HFOs for clinical decision making is that inter-reader agreement on what constitutes an event is often poor (Spring et al., 2017). Therefore, it is difficult to validate automated HFO detectors. Furthermore, HFO detectors may generate clinically informative results (Weiss et al., 2013, Weiss et al., 2015), in the absence of a gold standard comparison biomarker that can be confirmed visually. This paradox can perhaps be solved by developing a procedure to allow experts to code HFOs from iEEG signals using classes (Jacobs et al., 2008) that generate an agreed upon gold standard for evaluating automated procedures.
HFO detection is performed using recordings in bipolar montage in order to reduce artifact, originating from muscle or the reference electrode, that can mimic HFO events. There are no previously published studies that utilize an automated HFO detector to define events in macroelectrode recordings recorded in referential montage. In theory, performing HFO detection in referential montage provides two advantages. First, referential montage increases the spatial resolution of the iEEG as compared with bipolar montage by increasing the total number of recordings. Second, referential montage could improve HFO detection by increasing the signal to noise ratio, since the dipole generators of the HFO could be distributed across multiple macro-electrode sites, which could obscure the HFO due to in phase cancellation (de la Prida et al., 2016).
In order to identify HFOs in referential montage using an automated detector, it is essential to develop a strategy for reducing or eliminating artifact, that would otherwise be eliminated by a bipolar montage. Independent component analysis (ICA) is a signal processing approach that can separate signal sources based on minimizing mutual information, and assuring that each component has a non-Gaussian distribution (Bell and Sejnowski, 1995, Cardoso, 1997). ICA has been utilized to reduce artifact in scalp EEG (Delorme et al., 2007, McMenamin et al., 2010), and has also been used to remove the scalp reference signal from iEEG recordings (Hu et al., 2007). In this paper, we tested the hypothesis that applying ICA to band-pass filtered iEEG recordings in referential montage could be utilized to accurately detect ripple events even when the recordings are contaminated by artifact.
Section snippets
Patients
Recordings were selected from 11 patients who underwent intracranial monitoring with depth electrodes between 2014 and 2016 at University of California Los Angeles (UCLA) and from five patients at Thomas Jefferson University (TJU) in 2016–2017 for the purpose of localization of the seizure onset zone. The inclusion criteria were at least one night and day of intracranial recording with 2000 Hz sampling rate and at least 4 h of interictal EEG uninterrupted by seizures for the UCLA patients. For
Patient description
The recordings were selected from eleven patients who underwent intracranial monitoring with depth electrodes during sleep, and five patients in which the intracranial monitoring was performed in the operating room using depth electrodes. Among the patients in whom the sleep recordings were selected, three had unilateral mesial temporal lobe epilepsy (MTLE), three patients had bilateral MTLE, three patients had MTLE plus the involvement of neocortical regions, and two patients had neocortical
Discussion
In summary, herein we report that (1) by processing band-pass filtered (80–600 Hz) iEEG signals recorded in referential montage with Infomax ICA we could reduce or eliminate muscle artifact, and demarcate artefactual HFOs, (2) manual validation demonstrated that when a Hilbert detector was used to detect the ripple events in the post-ICA processed referential recordings, it was moderately sensitive and very precise, (3) both the true and false ripple on spike events defined using the detector
Conclusions
We report a novel approach to identify and classify ripple events in iEEG recordings in the referential montage that eliminates or dramatically reduces artifact using Infomax ICA. The detector could also separate true ripple on spike events from false ripple on spike events that result from filter ringing, and accurately define the magnitude and spectral content of the true events. The true and false ripple on spike rates defined by the detector accurately classified the SOZ in a diverse group
Acknowledgements
The authors would like to thank Mr. Dale Wyeth and Mr. Edmund Wyeth at Thomas Jefferson University, and Mr. Kirk Shattuck at University of California Los Angeles for their technical assistance with the experiments.
Conflict of interest statement
The author(s) declare(s) that there is no conflict of interest regarding the publication of this article.
Funding
This work was supported by NIH/NINDS K23NS094633 (SAW).
References (41)
- et al.
Pitfalls of high-pass filtering for detecting epileptic oscillations: a technical note on “false” ripples
Clin Neurophysiol
(2010) - et al.
The morphology of high frequency oscillations does not improve delineating the epileptogenic zone
Clin Neurophysiol
(2016) - et al.
Enhanced detection of artifacts in EEG data using higher-order statistics and independent component analysis
Neuroimage
(2007) - et al.
Universal automated high frequency oscillation detector for real-time, long term EEG
Clin Neurophysiol
(2016) - et al.
High-frequency oscillations (HFOs) in clinical epilepsy
Prog Neurobiol
(2012) - et al.
Validation of ICA-based myogenic artifact correction for scalp and source-localized EEG
Neuroimage
(2010) - et al.
Interrater reliability of visually evaluated high frequency oscillations
Clin Neurophysiol
(2017) - et al.
In vivo laminar electrophysiology co-registered with histology in the hippocampus of patients with temporal lobe epilepsy
Exp Neurol
(2004) - et al.
A comparison between detectors of high frequency oscillations
Clin Neurophysiol
(2012) - et al.
Ictal clinical and scalp-EEG findings differentiating temporal lobe epilepsies from temporal 'plus' epilepsies
Brain
(2007)
Temporal plus epilepsy is a major determinant of temporal lobe surgery failures
Brain
Information-maximization approach to blind separation and blind deconvolution
Technology
Local generation of fast ripples in epileptic brain
J Neurosci
Human intracranial high frequency oscillations detected by automatic time-frequency analysis
PLoS ONE
Infomax and maximum likelihood for blind source separation
IEEE Signal Process Lett
Brain oscillations and the importance of waveform shape
Trends Cognitive Sci
High-frequency oscillations: what is normal and what is not?
Epilepsia
High frequency oscillations: the new EEG frontier?
Epilepsia
Cited by (29)
Delta oscillation coupled propagating fast ripples precede epileptiform discharges in patients with focal epilepsy
2022, Neurobiology of DiseaseMEG detection of high frequency oscillations and intracranial-EEG validation in pediatric epilepsy surgery
2021, Clinical NeurophysiologyEffects of depth electrode montage and single-pulse electrical stimulation sites on neuronal responses and effective connectivity
2020, Clinical NeurophysiologyCitation Excerpt :Such an original reference is often referred to as a system reference (Lesser et al., 2010; Nariai et al., 2011). We clinically assessed iEEG signals using average montage (Crone et al., 2001) as well as bipolar montage (Zijlmans et al., 2011; Isnard et al., 2018; Shimamoto et al., 2018). We excluded SPES-induced neuronal responses recorded at SOZ, irritative zone, or artifactual channels from further analyses of the present study.
Patient phenotypes and clinical outcomes in invasive monitoring for epilepsy: An individual patient data meta-analysis
2020, Epilepsy and BehaviorCitation Excerpt :However, there was no significant difference (p = 0.0565) in the rate of seizure freedom between these groups regardless of resection at 53.6% (CI: 50.3–56.8%) of SEEG patients and 62.2% (CI: 59.0–65.3%) of SDE patients. Additionally, 29 [19–47] studies reported SEEG- or SDE-guided resection status for 548 patients (92.1%) within the total cohort (Table 2). Of these patients, 300 (54.7%) underwent SEEG and 248 (45.3%) underwent SDE.
Graph index complexity as a novel surrogate marker of high frequency oscillations in delineating the seizure onset zone
2020, Clinical NeurophysiologyDoes data cleaning improve brain state classification?
2019, Journal of Neuroscience MethodsCitation Excerpt :The resulting AUCs/t-statistics may be found in Fig. 4. Recently, a semi-automatic method has been evaluated to identify epochs with high frequency oscillatory (HFO) activity, a characteristic of epileptic EEG and iEEG (Waldman et al., 2018; Shimamoto et al., 2018; Weiss et al., 2018). The algorithm, run only on depth electrodes for optimal signal-to-noise, annotates an epoch with one of the following labels that are then visually validated: (0) No HFO event detected, (1) Ripple superimposed in an interictal discharge, (2) Sharply contoured epileptiform spike, (3) Ripple on oscillation, (4) Fast ripple superimposed in an interictal discharge, (5) Fast ripple on oscillation.