Elsevier

Clinical Neurophysiology

Volume 129, Issue 1, January 2018, Pages 296-307
Clinical Neurophysiology

Utilization of independent component analysis for accurate pathological ripple detection in intracranial EEG recordings recorded extra- and intra-operatively

https://doi.org/10.1016/j.clinph.2017.08.036Get rights and content

Highlights

  • Applying independent component analysis (ICA) to intracranial EEG following band-pass filtering (80–600 Hz) reduces artifact.

  • Ripple detection is precise after utilizing ICA to reduce and demarcate artifact.

  • Ripple rates are elevated in the seizure onset zone in recordings performed during sleep and intraoperatively.

Abstract

Objective

To develop and validate a detector that identifies ripple (80–200 Hz) events in intracranial EEG (iEEG) recordings in a referential montage and utilizes independent component analysis (ICA) to eliminate or reduce high-frequency artifact contamination. Also, investigate the correspondence of detected ripples and the seizure onset zone (SOZ).

Methods

iEEG recordings from 16 patients were first band-pass filtered (80–600 Hz) and Infomax ICA was next applied to derive the first independent component (IC1). IC1 was subsequently pruned, and an artifact index was derived to reduce the identification of high-frequency events introduced by the reference electrode signal. A Hilbert detector identified ripple events in the processed iEEG recordings using amplitude and duration criteria. The identified ripple events were further classified and characterized as true or false ripple on spikes, or ripples on oscillations by utilizing a topographical analysis to their time-frequency plot, and confirmed by visual inspection.

Results

The signal to noise ratio was improved by pruning IC1. The precision of the detector for ripple events was 91.27 ± 4.3%, and the sensitivity of the detector was 79.4 ± 3.0% (N = 16 patients, 5842 ripple events). The sensitivity and precision of the detector was equivalent in iEEG recordings obtained during sleep or intra-operatively. Across all the patients, true ripple on spike rates and also the rates of false ripple on spikes, that were generated due to filter ringing, classified the seizure onset zone (SOZ) with an area under the receiver operating curve (AUROC) of >76%. The magnitude and spectral content of true ripple on spikes generated in the SOZ was distinct as compared with the ripples generated in the NSOZ (p < .001).

Conclusions

Utilizing ICA to analyze iEEG recordings in referential montage provides many benefits to the study of high-frequency oscillations. The ripple rates and properties defined using this approach may accurately delineate the seizure onset zone.

Significance

Strategies to improve the spatial resolution of intracranial EEG and reduce artifact can help improve the clinical utility of HFO biomarkers.

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)

  • C. Barba et al.

    Temporal plus epilepsy is a major determinant of temporal lobe surgery failures

    Brain

    (2016)
  • A.J. Bell et al.

    Information-maximization approach to blind separation and blind deconvolution

    Technology

    (1995)
  • A. Bragin et al.

    Local generation of fast ripples in epileptic brain

    J Neurosci

    (2002)
  • S. Burnos et al.

    Human intracranial high frequency oscillations detected by automatic time-frequency analysis

    PLoS ONE

    (2014)
  • J.-F. Cardoso

    Infomax and maximum likelihood for blind source separation

    IEEE Signal Process Lett

    (1997)
  • S.R. Cole et al.

    Brain oscillations and the importance of waveform shape

    Trends Cognitive Sci

    (2016)
  • Davis J, Goadrich M. The relationship between precision-recall and ROC curves. In: ICML '06 proceedings of the 23rd...
  • Eissa TL, Tryba AK, Masrcuccilli CJ, Ben-Mabrouk F, Smith EH, Lew SM, et al. Multiscale aspects of generation of...
  • J. Engel et al.

    High-frequency oscillations: what is normal and what is not?

    Epilepsia

    (2009)
  • J. Gotman

    High frequency oscillations: the new EEG frontier?

    Epilepsia

    (2010)
  • Cited by (29)

    • Effects of depth electrode montage and single-pulse electrical stimulation sites on neuronal responses and effective connectivity

      2020, Clinical Neurophysiology
      Citation 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 Behavior
      Citation 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.

    • Does data cleaning improve brain state classification?

      2019, Journal of Neuroscience Methods
      Citation 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.

    View all citing articles on Scopus
    View full text