1
|
Zhang C, Hu W, Wu Y, Li G, Yang C, Wu T. Altered Directed-Connectivity Network in Temporal Lobe Epilepsy: A MEG Study. SENSORS (BASEL, SWITZERLAND) 2025; 25:1356. [PMID: 40096174 PMCID: PMC11902853 DOI: 10.3390/s25051356] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/08/2025] [Revised: 02/10/2025] [Accepted: 02/17/2025] [Indexed: 03/19/2025]
Abstract
Temporal lobe epilepsy (TLE) is considered a network disorder rather than a localized lesion, making it essential to study the network mechanisms underlying TLE. In this study, we constructed directed brain networks based on clinical MEG data using the Granger Causality Analysis (GCA) method, aiming to provide new insights into the network mechanisms of TLE. MEG data from 13 lTLE and 21 rTLE patients and 14 healthy controls (HCs) were analyzed. The preprocessed MEG data were used to construct directed brain networks using the GCA method and undirected brain networks using the Pearson Correlation Coefficient (PCC) method. Graph theoretical analysis extracted global and local topologies from the binary matrix, and SVM classified topologies with significant differences (p < 0.05). Comparative studies were performed on connectivity strengths, graph theory metrics, and SVM classifications between GCA and PCC, with an additional analysis of GCA-weighted network connectivity. The results show that TLE patients showed significantly increased functional connectivity based on GCA compared to the control group; similarities of the hub brain regions between lTLE and rTLE patients and the cortical-limbic-thalamic-cortical loop were identified; TLE patients exhibited a significant increase in GCA-based Global Clustering Coefficient (GCC) and Global Local Efficiency (GLE); most brain regions with abnormal local topological properties in TLE patients overlapped with their hub regions. The directionality of brain connectivity has played a significantly more pivotal role in research on TLE. GCA may be a potential tool in MEG analysis to distinguish TLE patients and HC effectively.
Collapse
Affiliation(s)
- Chen Zhang
- College of Chemistry and Life Science, Beijing University of Technology, Beijing 100124, China; (C.Z.); (Y.W.); (G.L.)
| | - Wenhan Hu
- Department of Neurosurgery, Tiantan Hospital, Beijing 100070, China;
| | - Yutong Wu
- College of Chemistry and Life Science, Beijing University of Technology, Beijing 100124, China; (C.Z.); (Y.W.); (G.L.)
| | - Guangfei Li
- College of Chemistry and Life Science, Beijing University of Technology, Beijing 100124, China; (C.Z.); (Y.W.); (G.L.)
| | - Chunlan Yang
- College of Chemistry and Life Science, Beijing University of Technology, Beijing 100124, China; (C.Z.); (Y.W.); (G.L.)
| | - Ting Wu
- Department of Radiology, Jiangsu Province Hospital of Chinese Medicine, Affiliated Hospital of Nanjing University of Chinese Medicine, Nanjing 210000, China
| |
Collapse
|
2
|
Liu X, Han J, Zhang X, Zhou Q, Huang Z, Wang Y, Zhang J, Lin Y. Dynamic evolution of frontal-temporal network connectivity in temporal lobe epilepsy: A magnetoencephalography study. Hum Brain Mapp 2024; 45:e70033. [PMID: 39319686 PMCID: PMC11423264 DOI: 10.1002/hbm.70033] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/03/2024] [Revised: 09/06/2024] [Accepted: 09/11/2024] [Indexed: 09/26/2024] Open
Abstract
Temporal lobe epilepsy (TLE) frequently involves an intricate, extensive epileptic frontal-temporal network. This study aimed to investigate the interactions between temporal and frontal regions and the dynamic patterns of the frontal-temporal network in TLE patients with different disease durations. The magnetoencephalography data of 36 postoperative seizure-free patients with long-term follow-up of at least 1 year, and 21 age- and sex-matched healthy subjects were included in this study. Patients were initially divided into LONG-TERM (n = 18, DURATION >10 years) and SHORT-TERM (n = 18, DURATION ≤10 years) groups based on 10-year disease duration. For reliability, supplementary analyses were conducted with alternative cutoffs, creating three groups: 0 < DURATION ≤7 years (n = 11), 7 < DURATION ≤14 years (n = 11), and DURATION >14 years (n = 14). This study examined the intraregional phase-amplitude coupling (PAC) between theta phase and alpha amplitude across the whole brain. The interregional directed phase transfer entropy (dPTE) between frontal and temporal regions in the alpha and theta bands, and the interregional cross-frequency directionality (CFD) between temporal and frontal regions from the theta phase to the alpha amplitude were further computed and compared among groups. Partial correlation analysis was conducted to investigate correlations between intraregional PAC, interregional dPTE connectivity, interregional CFD, and disease duration. Whole-brain intraregional PAC analyses revealed enhanced theta phase-alpha amplitude coupling within the ipsilateral temporal and frontal regions in TLE patients, and the ipsilateral temporal PAC was positively correlated with disease duration (r = 0.38, p <.05). Interregional dPTE analyses demonstrated a gradual increase in frontal-to-temporal connectivity within the alpha band, while the direction of theta-band connectivity reversed from frontal-to-temporal to temporal-to-frontal as the disease duration increased. Interregional CFD analyses revealed that the inhibitory effect of frontal regions on temporal regions gradually increased with prolonged disease duration (r = -0.36, p <.05). This study clarified the intrinsic reciprocal connectivity between temporal and frontal regions with TLE duration. We propose a dynamically reorganized triple-stage network that transitions from balanced networks to constrained networks and further develops into imbalanced networks as the disease duration increases.
Collapse
Affiliation(s)
- Xinyan Liu
- School of Biological Science and Medical EngineeringBeihang UniversityBeijingChina
- Beijing Advanced Innovation Center for Big Data‐Based Precision MedicineBeihang UniversityBeijingChina
- Beijing Advanced Innovation Centre for Biomedical EngineeringBeihang UniversityBeijingChina
| | - Jiaqi Han
- Department of Neurology, Xuanwu HospitalCapital Medical UniversityBeijingChina
| | - Xiating Zhang
- Department of Neurology, Xuanwu HospitalCapital Medical UniversityBeijingChina
- Department of Neurologythe First Hospital of Hebei Medical UniversityShijiazhuangHebeiChina
| | - Qilin Zhou
- Department of Neurology, Xuanwu HospitalCapital Medical UniversityBeijingChina
| | - Zhaoyang Huang
- Department of Neurology, Xuanwu HospitalCapital Medical UniversityBeijingChina
- Department of Neurologythe First Hospital of Hebei Medical UniversityShijiazhuangHebeiChina
| | - Yuping Wang
- Department of Neurology, Xuanwu HospitalCapital Medical UniversityBeijingChina
- Department of Neurologythe First Hospital of Hebei Medical UniversityShijiazhuangHebeiChina
- Beijing Key Laboratory of NeuromodulationXuanwu Hospital, Capital Medical UniversityBeijingChina
| | - Jicong Zhang
- School of Biological Science and Medical EngineeringBeihang UniversityBeijingChina
- Beijing Advanced Innovation Center for Big Data‐Based Precision MedicineBeihang UniversityBeijingChina
- Beijing Advanced Innovation Centre for Biomedical EngineeringBeihang UniversityBeijingChina
- Hefei Innovation Research InstituteBeihang UniversityBeijingChina
| | - Yicong Lin
- Department of Neurology, Xuanwu HospitalCapital Medical UniversityBeijingChina
- Department of Neurologythe First Hospital of Hebei Medical UniversityShijiazhuangHebeiChina
- Beijing Key Laboratory of NeuromodulationXuanwu Hospital, Capital Medical UniversityBeijingChina
| |
Collapse
|
3
|
Krishnan B, Tousseyn S, Taylor K, Wu G, Serletis D, Najm I, Bulacio J, Alexopoulos AV. Measurable transitions during seizures in intracranial EEG: A stereoelectroencephalography and SPECT study. Clin Neurophysiol 2024; 161:80-92. [PMID: 38452427 DOI: 10.1016/j.clinph.2024.02.022] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2023] [Revised: 02/05/2024] [Accepted: 02/17/2024] [Indexed: 03/09/2024]
Abstract
OBJECTIVE Ictal Single Photon Emission Computed Tomography (SPECT) and stereo-electroencephalography (SEEG) are diagnostic techniques used for the management of patients with drug-resistant focal epilepsies. While hyperperfusion patterns in ictal SPECT studies reveal seizure onset and propagation pathways, the role of ictal hypoperfusion remains poorly understood. The goal of this study was to systematically characterize the spatio-temporal information flow dynamics between differently perfused brain regions using stereo-EEG recordings. METHODS We identified seizure-free patients after resective epilepsy surgery who had prior ictal SPECT and SEEG investigations. We estimated directional connectivity between the epileptogenic-zone (EZ), non-resected areas of hyperperfusion, hypoperfusion, and baseline perfusion during the interictal, preictal, ictal, and postictal periods. RESULTS Compared to the background, we noted significant information flow (1) during the preictal period from the EZ to the baseline and hyperperfused regions, (2) during the ictal onset from the EZ to all three regions, and (3) during the period of seizure evolution from the area of hypoperfusion to all three regions. CONCLUSIONS Hypoperfused brain regions were found to indirectly interact with the EZ during the ictal period. SIGNIFICANCE Our unique study, combining intracranial electrophysiology and perfusion imaging, presents compelling evidence of dynamic changes in directional connectivity between brain regions during the transition from interictal to ictal states.
Collapse
Affiliation(s)
- Balu Krishnan
- Neurological Institute, Epilepsy Center, Cleveland Clinic, Cleveland, OH, USA.
| | - Simon Tousseyn
- Academic Center for Epileptology, Kempenhaeghe and Maastricht UMC+, Heeze, The Netherlands; School for Mental Health and Neuroscience (MHeNs), University Maastricht (UM), Maastricht, The Netherlands
| | - Kenneth Taylor
- Neurological Institute, Epilepsy Center, Cleveland Clinic, Cleveland, OH, USA
| | - Guiyun Wu
- Neurological Institute, Epilepsy Center, Cleveland Clinic, Cleveland, OH, USA
| | - Demitre Serletis
- Neurological Institute, Epilepsy Center, Cleveland Clinic, Cleveland, OH, USA
| | - Imad Najm
- Neurological Institute, Epilepsy Center, Cleveland Clinic, Cleveland, OH, USA
| | - Juan Bulacio
- Neurological Institute, Epilepsy Center, Cleveland Clinic, Cleveland, OH, USA
| | | |
Collapse
|
4
|
Fernández-Martín R, Feys O, Juvené E, Aeby A, Urbain C, De Tiège X, Wens V. Towards the automated detection of interictal epileptiform discharges with magnetoencephalography. J Neurosci Methods 2024; 403:110052. [PMID: 38151188 DOI: 10.1016/j.jneumeth.2023.110052] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2023] [Revised: 12/08/2023] [Accepted: 12/21/2023] [Indexed: 12/29/2023]
Abstract
BACKGROUND The analysis of clinical magnetoencephalography (MEG) in patients with epilepsy traditionally relies on visual identification of interictal epileptiform discharges (IEDs), which is time consuming and dependent on subjective criteria. NEW METHOD Here, we explore the ability of Independent Components Analysis (ICA) and Hidden Markov Modeling (HMM) to automatically detect and localize IEDs. We tested our pipelines on resting-state MEG recordings from 10 school-aged children with (multi)focal epilepsy. RESULTS In focal epilepsy patients, both pipelines successfully detected visually identified IEDs, but also revealed unidentified low-amplitude IEDs. Success was more mitigated in patients with multifocal epilepsy, as our automated pipeline missed IED activity associated with some foci-an issue that could be alleviated by post-hoc manual selection of epileptiform ICs or HMM states. COMPARISON WITH EXISTING METHODS We compared our results with visual IED detection by an experienced clinical magnetoencephalographer, getting heightened sensitivity and requiring minimal input from clinical practitioners. CONCLUSIONS IED detection based on ICA or HMM represents an efficient way to identify IED localization and timing. The development of these automatic IED detection algorithms provide a step forward in clinical MEG practice by decreasing the duration of MEG analysis and enhancing its sensitivity.
Collapse
Affiliation(s)
- Raquel Fernández-Martín
- Université libre de Bruxelles (ULB), ULB Neuroscience Institute (UNI), Laboratoire de Neuroanatomie et de Neuroimagerie translationnelles (LNbT), Brussels, Belgium.
| | - Odile Feys
- Université libre de Bruxelles (ULB), ULB Neuroscience Institute (UNI), Laboratoire de Neuroanatomie et de Neuroimagerie translationnelles (LNbT), Brussels, Belgium; Université libre de Bruxelles (ULB), Hôpital Universitaire de Bruxelles (H.U.B.), Hôpital Erasme, Department of Neurology, Brussels, Belgium
| | - Elodie Juvené
- Université libre de Bruxelles (ULB), Hôpital Universitaire de Bruxelles (H.U.B.), Department of Pediatric Neurology, Brussels, Belgium
| | - Alec Aeby
- Université libre de Bruxelles (ULB), Hôpital Universitaire de Bruxelles (H.U.B.), Department of Pediatric Neurology, Brussels, Belgium
| | - Charline Urbain
- Université libre de Bruxelles (ULB), ULB Neuroscience Institute (UNI), Laboratoire de Neuroanatomie et de Neuroimagerie translationnelles (LNbT), Brussels, Belgium; Université libre de Bruxelles (ULB), ULB Neuroscience Institute (UNI), Centre for Research in Cognition and Neurosciences (CRCN), Neuropsychology and Functional Neuroimaging Research Unit (UR2NF), Brussels, Belgium
| | - Xavier De Tiège
- Université libre de Bruxelles (ULB), ULB Neuroscience Institute (UNI), Laboratoire de Neuroanatomie et de Neuroimagerie translationnelles (LNbT), Brussels, Belgium; Université libre de Bruxelles (ULB), Hôpital Universitaire de Bruxelles (H.U.B.), Hôpital Erasme, Service of translational Neuroimaging, Brussels, Belgium
| | - Vincent Wens
- Université libre de Bruxelles (ULB), ULB Neuroscience Institute (UNI), Laboratoire de Neuroanatomie et de Neuroimagerie translationnelles (LNbT), Brussels, Belgium; Université libre de Bruxelles (ULB), Hôpital Universitaire de Bruxelles (H.U.B.), Hôpital Erasme, Service of translational Neuroimaging, Brussels, Belgium
| |
Collapse
|
5
|
Ye S, Bagić A, He B. Disentanglement of Resting State Brain Networks for Localizing Epileptogenic Zone in Focal Epilepsy. Brain Topogr 2024; 37:152-168. [PMID: 38112884 PMCID: PMC10771380 DOI: 10.1007/s10548-023-01025-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2023] [Accepted: 11/20/2023] [Indexed: 12/21/2023]
Abstract
The objective of this study is to extract pathological brain networks from interictal period of E/MEG recordings to localize epileptic foci for presurgical evaluation. We proposed here a resting state E/MEG analysis framework, to disentangle brain functional networks represented by neural oscillations. By using an Embedded Hidden Markov Model, we constructed a state space for resting state recordings consisting of brain states with different spatiotemporal patterns. Functional connectivity analysis along with graph theory was applied on the extracted brain states to quantify the network features of the extracted brain states, based on which the source location of pathological states is determined. The method is evaluated by computer simulations and our simulation results revealed the proposed framework can extract brain states with high accuracy regarding both spatial and temporal profiles. We further evaluated the framework as compared with intracranial EEG defined seizure onset zone in 10 patients with drug-resistant focal epilepsy who underwent MEG recordings and were seizure free after surgical resection. The real patient data analysis showed very good localization results using the extracted pathological brain states in 6/10 patients, with localization error of about 15 mm as compared to the seizure onset zone. We show that the pathological brain networks can be disentangled from the resting-state electromagnetic recording and could be identified based on the connectivity features. The framework can serve as a useful tool in extracting brain functional networks from noninvasive resting state electromagnetic recordings, and promises to offer an alternative to aid presurgical evaluation guiding intracranial EEG electrodes implantation.
Collapse
Affiliation(s)
- Shuai Ye
- Department of Biomedical Engineering, Carnegie Mellon University, 5000 Forbes Avenue, Pittsburgh, PA, 15213, USA
| | - Anto Bagić
- Department of Neurology, University of Pittsburgh Comprehensive Epilepsy Center (UPCEC), University of Pittsburgh Medical School, Pittsburgh, PA, USA
| | - Bin He
- Department of Biomedical Engineering, Carnegie Mellon University, 5000 Forbes Avenue, Pittsburgh, PA, 15213, USA.
| |
Collapse
|
6
|
Krishnan B, Tousseyn S, Wang ZI, Murakami H, Wu G, Burgess R, Iasemidis L, Najm I, Alexopoulos AV. Novel noninvasive identification of patient-specific epileptic networks in focal epilepsies: Linking single-photon emission computed tomography perfusion during seizures with resting-state magnetoencephalography dynamics. Hum Brain Mapp 2023; 44:1695-1710. [PMID: 36480260 PMCID: PMC9921232 DOI: 10.1002/hbm.26168] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/24/2022] [Revised: 08/31/2022] [Accepted: 11/18/2022] [Indexed: 12/13/2022] Open
Abstract
Single-photon emission computed tomography (SPECT) during seizures and magnetoencephalography (MEG) during the interictal state are noninvasive modalities employed in the localization of the epileptogenic zone in patients with drug-resistant focal epilepsy (DRFE). The present study aims to investigate whether there exists a preferentially high MEG functional connectivity (FC) among those regions of the brain that exhibit hyperperfusion or hypoperfusion during seizures. We studied MEG and SPECT data in 30 consecutive DRFE patients who had resective epilepsy surgery. We parcellated each ictal perfusion map into 200 regions of interest (ROIs) and generated ROI time series using source modeling of MEG data. FC between ROIs was quantified using coherence and phase-locking value. We defined a generalized linear model to relate the connectivity of each ROI, ictal perfusion z score, and distance between ROIs. We compared the coefficients relating perfusion z score to FC of each ROI and estimated the connectivity within and between resected and unresected ROIs. We found that perfusion z scores were strongly correlated with the FC of hyper-, and separately, hypoperfused ROIs across patients. High interictal connectivity was observed between hyperperfused brain regions inside and outside the resected area. High connectivity was also observed between regions of ictal hypoperfusion. Importantly, the ictally hypoperfused regions had a low interictal connectivity to regions that became hyperperfused during seizures. We conclude that brain regions exhibiting hyperperfusion during seizures highlight a preferentially connected interictal network, whereas regions of ictal hypoperfusion highlight a separate, discrete and interconnected, interictal network.
Collapse
Affiliation(s)
- Balu Krishnan
- Neurological InstituteEpilepsy Center, Cleveland ClinicClevelandOhioUSA
| | - Simon Tousseyn
- Academic Center for EpileptologyKempenhaeghe and Maastricht UMC+HeezeThe Netherlands
| | - Zhong Irene Wang
- Neurological InstituteEpilepsy Center, Cleveland ClinicClevelandOhioUSA
| | - Hiroatsu Murakami
- Neurological InstituteEpilepsy Center, Cleveland ClinicClevelandOhioUSA
| | - Guiyun Wu
- Neurological InstituteEpilepsy Center, Cleveland ClinicClevelandOhioUSA
| | - Richard Burgess
- Neurological InstituteEpilepsy Center, Cleveland ClinicClevelandOhioUSA
| | - Leonidas Iasemidis
- Department of Translational NeuroscienceBarrow Neurological InstituteScottsdaleArizonaUSA
- Department of NeurologyBarrow Neurological InstituteScottsdaleArizonaUSA
| | - Imad Najm
- Neurological InstituteEpilepsy Center, Cleveland ClinicClevelandOhioUSA
| | | |
Collapse
|
7
|
Rampp S, Kaltenhäuser M, Müller-Voggel N, Doerfler A, Kasper BS, Hamer HM, Brandner S, Buchfelder M. MEG Node Degree for Focus Localization: Comparison with Invasive EEG. Biomedicines 2023; 11:biomedicines11020438. [PMID: 36830974 PMCID: PMC9953213 DOI: 10.3390/biomedicines11020438] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2023] [Revised: 01/23/2023] [Accepted: 01/30/2023] [Indexed: 02/05/2023] Open
Abstract
Epilepsy surgery is a viable therapy option for patients with pharmacoresistant focal epilepsies. A prerequisite for postoperative seizure freedom is the localization of the epileptogenic zone, e.g., using electro- and magnetoencephalography (EEG/MEG). Evidence shows that resting state MEG contains subtle alterations, which may add information to the workup of epilepsy surgery. Here, we investigate node degree (ND), a graph-theoretical parameter of functional connectivity, in relation to the seizure onset zone (SOZ) determined by invasive EEG (iEEG) in a consecutive series of 50 adult patients. Resting state data were subjected to whole brain, all-to-all connectivity analysis using the imaginary part of coherence. Graphs were described using parcellated ND. SOZ localization was investigated on a lobar and sublobar level. On a lobar level, all frequency bands except alpha showed significantly higher maximal ND (mND) values inside the SOZ compared to outside (ratios 1.11-1.20, alpha 1.02). Area-under-the-curve (AUC) was 0.67-0.78 for all expected alpha (0.44, ns). On a sublobar level, mND inside the SOZ was higher for all frequency bands (1.13-1.38, AUC 0.58-0.78) except gamma (1.02). MEG ND is significantly related to SOZ in delta, theta and beta bands. ND may provide new localization tools for presurgical evaluation of epilepsy surgery.
Collapse
Affiliation(s)
- Stefan Rampp
- Department of Neurosurgery, University Hospital Erlangen, 91054 Erlangen, Germany
- Department of Neurosurgery, University Hospital Halle (Saale), 06120 Halle (Saale), Germany
- Correspondence: ; Tel.: +49-9131-85-46921; Fax: +49-9131-85-34476
| | - Martin Kaltenhäuser
- Department of Neurosurgery, University Hospital Erlangen, 91054 Erlangen, Germany
| | - Nadia Müller-Voggel
- Department of Neurosurgery, University Hospital Erlangen, 91054 Erlangen, Germany
| | - Arnd Doerfler
- Department of Neuroradiology, University Hospital Erlangen, 91054 Erlangen, Germany
| | - Burkhard S. Kasper
- Epilepsy Center, Department of Neurology, University Hospital Erlangen, 91054 Erlangen, Germany
| | - Hajo M. Hamer
- Epilepsy Center, Department of Neurology, University Hospital Erlangen, 91054 Erlangen, Germany
| | - Sebastian Brandner
- Department of Neurosurgery, University Hospital Erlangen, 91054 Erlangen, Germany
| | - Michael Buchfelder
- Department of Neurosurgery, University Hospital Erlangen, 91054 Erlangen, Germany
| |
Collapse
|
8
|
Alamoudi OA, Ilyas A, Pati S, Iasemidis L. Interictal localization of the epileptogenic zone: Utilizing the observed resonance behavior in the spectral band of surrounding inhibition. Front Neurosci 2022; 16:993678. [PMID: 36578827 PMCID: PMC9791262 DOI: 10.3389/fnins.2022.993678] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2022] [Accepted: 11/16/2022] [Indexed: 12/14/2022] Open
Abstract
Introduction The gold standard for identification of the epileptogenic zone (EZ) continues to be the visual inspection of electrographic changes around seizures' onset by experienced electroencephalography (EEG) readers. Development of an epileptogenic focus localization tool that can delineate the EZ from analysis of interictal (seizure-free) periods is still an open question of great significance for improved diagnosis (e.g., presurgical evaluation) and treatment of epilepsy (e.g., surgical outcome). Methods We developed an EZ interictal localization algorithm (EZILA) based on novel analysis of intracranial EEG (iEEG) using a univariate periodogram-type power measure, a straight-forward ranking approach, a robust dimensional reduction method and a clustering technique. Ten patients with temporal and extra temporal lobe epilepsies, and matching the inclusion criteria of having iEEG recordings at the epilepsy monitoring unit (EMU) and being Engel Class I ≥12 months post-surgery, were recruited in this study. Results In a nested k-fold cross validation statistical framework, EZILA assigned the highest score to iEEG channels within the EZ in all patients (10/10) during the first hour of the iEEG recordings and up to their first typical clinical seizure in the EMU (i.e., early interictal period). To further validate EZILA's performance, data from two new (Engel Class I) patients were analyzed in a double-blinded fashion; the EZILA successfully localized iEEG channels within the EZ from interictal iEEG in both patients. Discussion Out of the sampled brain regions, iEEG channels in the EZ were most frequently and maximally active in seizure-free (interictal) periods across patients in specific narrow gamma frequency band (∼60-80 Hz), which we have termed focal frequency band (FFB). These findings are consistent with the hypothesis that the EZ may interictally be regulated (controlled) by surrounding inhibitory neurons with resonance characteristics within this narrow gamma band.
Collapse
Affiliation(s)
- Omar A. Alamoudi
- Biomedical Engineering Program, Faculty of Engineering, King Abdulaziz University, Jeddah, Saudi Arabia,Neurology Department, Texas Institute for Restorative Neurotechnologies (TIRN), University of Texas Medical School, Houston, TX, United States,*Correspondence: Omar A. Alamoudi,
| | - Adeel Ilyas
- Neurology Department, Texas Institute for Restorative Neurotechnologies (TIRN), University of Texas Medical School, Houston, TX, United States,Department of Neurological Surgery, University of Alabama at Birmingham, Birmingham, AL, United States,Vivian L. Smith Department of Neurosurgery, McGovern Medical School at University of Texas (UT) Health Houston, Houston, TX, United States
| | - Sandipan Pati
- Neurology Department, Texas Institute for Restorative Neurotechnologies (TIRN), University of Texas Medical School, Houston, TX, United States
| | - Leon Iasemidis
- Biomedical Engineering Department, Arizona State University, Tempe, AZ, United States,Department of Translational Neuroscience, Barrow Neurological Institute, Phoenix, AZ, United States
| |
Collapse
|
9
|
Ren Z, Han X, Wang B. The performance evaluation of the state-of-the-art EEG-based seizure prediction models. Front Neurol 2022; 13:1016224. [PMID: 36504642 PMCID: PMC9732735 DOI: 10.3389/fneur.2022.1016224] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2022] [Accepted: 11/09/2022] [Indexed: 11/26/2022] Open
Abstract
The recurrent and unpredictable nature of seizures can lead to unintentional injuries and even death. The rapid development of electroencephalogram (EEG) and Artificial Intelligence (AI) technologies has made it possible to predict seizures in real-time through brain-machine interfaces (BCI), allowing advanced intervention. To date, there is still much room for improvement in predictive seizure models constructed by EEG using machine learning (ML) and deep learning (DL). But, the most critical issue is how to improve the performance and generalization of the model, which involves some confusing conceptual and methodological issues. This review focuses on analyzing several factors affecting the performance of seizure prediction models, focusing on the aspects of post-processing, seizure occurrence period (SOP), seizure prediction horizon (SPH), and algorithms. Furthermore, this study presents some new directions and suggestions for building high-performance prediction models in the future. We aimed to clarify the concept for future research in related fields and improve the performance of prediction models to provide a theoretical basis for future applications of wearable seizure detection devices.
Collapse
Affiliation(s)
- Zhe Ren
- Department of Neurology, Zhengzhou University People's Hospital, Zhengzhou, China,Department of Neurology, Henan Provincial People's Hospital, Zhengzhou, China
| | - Xiong Han
- Department of Neurology, Zhengzhou University People's Hospital, Zhengzhou, China,Department of Neurology, Henan Provincial People's Hospital, Zhengzhou, China,*Correspondence: Xiong Han
| | - Bin Wang
- Department of Neurology, Zhengzhou University People's Hospital, Zhengzhou, China,Department of Neurology, Henan Provincial People's Hospital, Zhengzhou, China
| |
Collapse
|
10
|
Cao M, Vogrin SJ, Peterson ADH, Woods W, Cook MJ, Plummer C. Dynamical Network Models From EEG and MEG for Epilepsy Surgery—A Quantitative Approach. Front Neurol 2022; 13:837893. [PMID: 35422755 PMCID: PMC9001937 DOI: 10.3389/fneur.2022.837893] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2021] [Accepted: 03/01/2022] [Indexed: 11/16/2022] Open
Abstract
There is an urgent need for more informative quantitative techniques that non-invasively and objectively assess strategies for epilepsy surgery. Invasive intracranial electroencephalography (iEEG) remains the clinical gold standard to investigate the nature of the epileptogenic zone (EZ) before surgical resection. However, there are major limitations of iEEG, such as the limited spatial sampling and the degree of subjectivity inherent in the analysis and clinical interpretation of iEEG data. Recent advances in network analysis and dynamical network modeling provide a novel aspect toward a more objective assessment of the EZ. The advantage of such approaches is that they are data-driven and require less or no human input. Multiple studies have demonstrated success using these approaches when applied to iEEG data in characterizing the EZ and predicting surgical outcomes. However, the limitations of iEEG recordings equally apply to these studies—limited spatial sampling and the implicit assumption that iEEG electrodes, whether strip, grid, depth or stereo EEG (sEEG) arrays, are placed in the correct location. Therefore, it is of interest to clinicians and scientists to see whether the same analysis and modeling techniques can be applied to whole-brain, non-invasive neuroimaging data (from MRI-based techniques) and neurophysiological data (from MEG and scalp EEG recordings), thus removing the limitation of spatial sampling, while safely and objectively characterizing the EZ. This review aims to summarize current state of the art non-invasive methods that inform epilepsy surgery using network analysis and dynamical network models. We also present perspectives on future directions and clinical applications of these promising approaches.
Collapse
Affiliation(s)
- Miao Cao
- Center for MRI Research, Peking University, Beijing, China
- Department of Medicine, The University of Melbourne, Melbourne, VIC, Australia
- Centre for Clinical Neurosciences and Neurological Research, St Vincent's Hospital Melbourne, Melbourne, VIC, Australia
| | - Simon J. Vogrin
- Department of Medicine, The University of Melbourne, Melbourne, VIC, Australia
- Centre for Clinical Neurosciences and Neurological Research, St Vincent's Hospital Melbourne, Melbourne, VIC, Australia
- School of Health Sciences, Swinburne University of Technology, Melbourne, VIC, Australia
| | - Andre D. H. Peterson
- Department of Medicine, The University of Melbourne, Melbourne, VIC, Australia
- Centre for Clinical Neurosciences and Neurological Research, St Vincent's Hospital Melbourne, Melbourne, VIC, Australia
| | - William Woods
- School of Health Sciences, Swinburne University of Technology, Melbourne, VIC, Australia
| | - Mark J. Cook
- Department of Medicine, The University of Melbourne, Melbourne, VIC, Australia
- Centre for Clinical Neurosciences and Neurological Research, St Vincent's Hospital Melbourne, Melbourne, VIC, Australia
| | - Chris Plummer
- Department of Medicine, The University of Melbourne, Melbourne, VIC, Australia
- Centre for Clinical Neurosciences and Neurological Research, St Vincent's Hospital Melbourne, Melbourne, VIC, Australia
- School of Health Sciences, Swinburne University of Technology, Melbourne, VIC, Australia
- *Correspondence: Chris Plummer
| |
Collapse
|
11
|
Sorelli M, Hutson TN, Iasemidis L, Bocchi L. Linear and Nonlinear Directed Connectivity Analysis of the Cardio-Respiratory System in Type 1 Diabetes. FRONTIERS IN NETWORK PHYSIOLOGY 2022; 2:840829. [PMID: 36926087 PMCID: PMC10013013 DOI: 10.3389/fnetp.2022.840829] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/21/2021] [Accepted: 02/14/2022] [Indexed: 12/31/2022]
Abstract
In this study, we explored the possibility of developing non-invasive biomarkers for patients with type 1 diabetes (T1D) by quantifying the directional couplings between the cardiac, vascular, and respiratory systems, treating them as interconnected nodes in a network configuration. Towards this goal, we employed a linear directional connectivity measure, the directed transfer function (DTF), estimated by a linear multivariate autoregressive modelling of ECG, respiratory and skin perfusion signals, and a nonlinear method, the dynamical Bayesian inference (DBI) analysis of bivariate phase interactions. The physiological data were recorded concurrently for a relatively short time period (5 min) from 10 healthy control subjects and 10 T1D patients. We found that, in both control and T1D subjects, breathing had greater influence on the heart and perfusion with respect to the opposite coupling direction and that, by both employed methods of analysis, the causal influence of breathing on the heart was significantly decreased (p < 0.05) in T1D patients compared to the control group. These preliminary results, although obtained from a limited number of subjects, provide a strong indication for the usefulness of a network-based multi-modal analysis for the development of biomarkers of T1D-related complications from short-duration data, as well as their potential in the exploration of the pathophysiological mechanisms that underlie this devastating and very widespread disease.
Collapse
Affiliation(s)
- Michele Sorelli
- European Laboratory for Non-Linear Spectroscopy, University of Florence, Florence, Italy
- Department of Physics and Astronomy, University of Florence, Florence, Italy
| | - T. Noah Hutson
- Department of Translational Neuroscience, Barrow Neurological Institute, Phoenix, AZ, United States
| | - Leonidas Iasemidis
- Department of Translational Neuroscience, Barrow Neurological Institute, Phoenix, AZ, United States
| | - Leonardo Bocchi
- European Laboratory for Non-Linear Spectroscopy, University of Florence, Florence, Italy
- Department of Information Engineering, University of Florence, Florence, Italy
| |
Collapse
|
12
|
Vogel S, Kaltenhäuser M, Kim C, Müller-Voggel N, Rössler K, Dörfler A, Schwab S, Hamer H, Buchfelder M, Rampp S. MEG Node Degree Differences in Patients with Focal Epilepsy vs. Controls-Influence of Experimental Conditions. Brain Sci 2021; 11:1590. [PMID: 34942895 PMCID: PMC8699109 DOI: 10.3390/brainsci11121590] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2021] [Revised: 11/25/2021] [Accepted: 11/27/2021] [Indexed: 11/16/2022] Open
Abstract
Drug-resistant epilepsy can be most limiting for patients, and surgery represents a viable therapy option. With the growing research on the human connectome and the evidence of epilepsy being a network disorder, connectivity analysis may be able to contribute to our understanding of epilepsy and may be potentially developed into clinical applications. In this magnetoencephalographic study, we determined the whole-brain node degree of connectivity levels in patients and controls. Resting-state activity was measured at five frequency bands in 15 healthy controls and 15 patients with focal epilepsy of different etiologies. The whole-brain all-to-all imaginary part of coherence in source space was then calculated. Node degree was determined and parcellated and was used for further statistical evaluation. In comparison to controls, we found a significantly higher overall node degree in patients with lesional and non-lesional epilepsy. Furthermore, we examined the conditions of high/reduced vigilance and open/closed eyes in controls, to analyze whether patient node degree levels can be achieved. We evaluated intraclass-correlation statistics (ICC) to evaluate the reproducibility. Connectivity and specifically node degree analysis could present new tools for one of the most common neurological diseases, with potential applications in epilepsy diagnostics.
Collapse
Affiliation(s)
- Stephan Vogel
- Department of Neurosurgery, University Hospital Erlangen, 91054 Erlangen, Germany; (M.K.); (C.K.); (N.M.-V.); (M.B.); (S.R.)
- Friedrich Alexander University Erlangen Nürnberg (FAU), 91054 Erlangen, Germany
| | - Martin Kaltenhäuser
- Department of Neurosurgery, University Hospital Erlangen, 91054 Erlangen, Germany; (M.K.); (C.K.); (N.M.-V.); (M.B.); (S.R.)
| | - Cora Kim
- Department of Neurosurgery, University Hospital Erlangen, 91054 Erlangen, Germany; (M.K.); (C.K.); (N.M.-V.); (M.B.); (S.R.)
| | - Nadia Müller-Voggel
- Department of Neurosurgery, University Hospital Erlangen, 91054 Erlangen, Germany; (M.K.); (C.K.); (N.M.-V.); (M.B.); (S.R.)
| | - Karl Rössler
- Department of Neurosurgery, Medical University Vienna, 1090 Vienna, Austria;
| | - Arnd Dörfler
- Department of Neuroradiology, University Hospital Erlangen, 91054 Erlangen, Germany;
| | - Stefan Schwab
- Department of Neurology, University Hospital Erlangen, 91054 Erlangen, Germany; (S.S.); (H.H.)
| | - Hajo Hamer
- Department of Neurology, University Hospital Erlangen, 91054 Erlangen, Germany; (S.S.); (H.H.)
| | - Michael Buchfelder
- Department of Neurosurgery, University Hospital Erlangen, 91054 Erlangen, Germany; (M.K.); (C.K.); (N.M.-V.); (M.B.); (S.R.)
| | - Stefan Rampp
- Department of Neurosurgery, University Hospital Erlangen, 91054 Erlangen, Germany; (M.K.); (C.K.); (N.M.-V.); (M.B.); (S.R.)
- Department of Neurosurgery, University Hospital Halle (Saale), 06120 Halle (Saale), Germany
| |
Collapse
|
13
|
Xu N, Shan W, Qi J, Wu J, Wang Q. Presurgical Evaluation of Epilepsy Using Resting-State MEG Functional Connectivity. Front Hum Neurosci 2021; 15:649074. [PMID: 34276321 PMCID: PMC8283278 DOI: 10.3389/fnhum.2021.649074] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/03/2021] [Accepted: 06/07/2021] [Indexed: 11/21/2022] Open
Abstract
Epilepsy is caused by abnormal electrical discharges (clinically identified by electrophysiological recording) in a specific part of the brain [originating in only one part of the brain, namely, the epileptogenic zone (EZ)]. Epilepsy is now defined as an archetypical hyperexcited neural network disorder. It can be investigated through the network analysis of interictal discharges, ictal discharges, and resting-state functional connectivity. Currently, there is an increasing interest in embedding resting-state connectivity analysis into the preoperative evaluation of epilepsy. Among the various neuroimaging technologies employed to achieve brain functional networks, magnetoencephalography (MEG) with the excellent temporal resolution is an ideal tool for estimating the resting-state connectivity between brain regions, which can reveal network abnormalities in epilepsy. What value does MEG resting-state functional connectivity offer for epileptic presurgical evaluation? Regarding this topic, this paper introduced the origin of MEG and the workflow of constructing source-space functional connectivity based on MEG signals. Resting-state functional connectivity abnormalities correlate with epileptogenic networks, which are defined by the brain regions involved in the production and propagation of epileptic activities. This paper reviewed the evidence of altered epileptic connectivity based on low- or high-frequency oscillations (HFOs) and the evidence of the advantage of using simultaneous MEG and intracranial electroencephalography (iEEG) recordings. More importantly, this review highlighted that MEG-based resting-state functional connectivity has the potential to predict postsurgical outcomes. In conclusion, resting-state MEG functional connectivity has made a substantial progress toward serving as a candidate biomarker included in epileptic presurgical evaluations.
Collapse
Affiliation(s)
- Na Xu
- Department of Neurology, Beijing Tiantan Hospital, Capital Medical University, Beijing, China
| | - Wei Shan
- Department of Neurology, Beijing Tiantan Hospital, Capital Medical University, Beijing, China
| | - Jing Qi
- Department of Neurology, Beijing Tiantan Hospital, Capital Medical University, Beijing, China
| | - Jianping Wu
- Department of Neurology, Beijing Tiantan Hospital, Capital Medical University, Beijing, China
- National Clinical Research Center for Neurological Diseases, Beijing, China
- Advanced Innovation Center for Human Brain Protection, Capital Medical University, Beijing, China
| | - Qun Wang
- Department of Neurology, Beijing Tiantan Hospital, Capital Medical University, Beijing, China
- National Clinical Research Center for Neurological Diseases, Beijing, China
- Beijing Institute of Brain Disorders, Collaborative Innovation Center for Brain Disorders, Capital Medical University, Beijing, China
- Beijing Key Laboratory of Neuromodulation, Beijing, China
| |
Collapse
|
14
|
Lin Y, Du P, Sun H, Liang Y, Wang Z, Cui Y, Chen K, Xia Y, Yao D, Yu L, Guo D. Identifying Refractory Epilepsy Without Structural Abnormalities by Fusing the Common Spatial Patterns of Functional and Effective EEG Networks. IEEE Trans Neural Syst Rehabil Eng 2021; 29:708-717. [PMID: 33830925 DOI: 10.1109/tnsre.2021.3071785] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
Drug refractory epilepsy (RE) is believed to be associated with structural lesions, but some RE patients show no significant structural abnormalities (RE-no-SA) on conventional magnetic resonance imaging scans. Since most of the medically controlled epilepsy (MCE) patients also do not exhibit structural abnormalities, a reliable assessment needs to be developed to differentiate RE-no-SA patients and MCE patients to avoid misdiagnosis and inappropriate treatment. Using resting-state scalp electroencephalogram (EEG) datasets, we extracted the spatial pattern of network (SPN) features from the functional and effective EEG networks of both RE-no-SA patients and MCE patients. Compared to the performance of traditional resting-state EEG network properties, the SPN features exhibited remarkable superiority in classifying these two groups of epilepsy patients, and accuracy values of 90.00% and 80.00% were obtained for the SPN features of the functional and effective EEG networks, respectively. By further fusing the SPN features of functional and effective networks, we demonstrated that the highest accuracy value of 96.67% could be reached, with a sensitivity of 100% and specificity of 92.86%. Overall, these findings not only indicate that the fused functional and effective SPN features are promising as reliable measurements for distinguishing RE-no-SA patients and MCE patients but also may provide a new perspective to explore the complex neurophysiology of refractory epilepsy.
Collapse
|
15
|
Hutson TN, Rezaei F, Gautier NM, Indumathy J, Glasscock E, Iasemidis L. Directed Connectivity Analysis of the Neuro-Cardio- and Respiratory Systems Reveals Novel Biomarkers of Susceptibility to SUDEP. IEEE OPEN JOURNAL OF ENGINEERING IN MEDICINE AND BIOLOGY 2020; 1:301-311. [PMID: 34223181 PMCID: PMC8249082 DOI: 10.1109/ojemb.2020.3036544] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2020] [Revised: 10/16/2020] [Accepted: 11/02/2020] [Indexed: 01/11/2023] Open
Abstract
Goal: Sudden unexpected death in epilepsy (SUDEP) is the leading cause of epilepsy-related mortality and its pathophysiological mechanisms remain unknown. We set to record and analyze for the first time concurrent electroencephalographic (EEG), electrocardiographic (ECG), and unrestrained whole-body plethysmographic (Pleth) signals from control (WT - wild type) and SUDEP-prone mice (KO- knockout Kcna1 animal model). Employing multivariate autoregressive models (MVAR) we measured all tri-organ effective directional interactions by the generalized partial directed coherence (GPDC) in the frequency domain over time (hours). When compared to the control (WT) animals, the SUDEP-prone (KO) animals exhibited (p < 0.001) reduced afferent and efferent interactions between the heart and the brain over the full frequency spectrum (0-200Hz), enhanced efferent interactions from the brain to the lungs and from the heart to the lungs at high (>90 Hz) frequencies (especially during periods with seizure activity), and decreased feedback from the lungs to the brain at low (<40 Hz) frequencies. These results show that impairment in the afferent and efferent pathways in the holistic neuro-cardio-respiratory network could lead to SUDEP, and effective connectivity measures and their dynamics could serve as novel biomarkers of susceptibility to SUDEP and seizures respectively.
Collapse
Affiliation(s)
- T. Noah Hutson
- Department of Biomedical EngineeringLouisiana Tech UniversityRustonLA71272USA
| | - Farnaz Rezaei
- Department of Mathematics and StatisticsLouisiana Tech UniversityRustonLA71272USA
| | - Nicole M. Gautier
- Department of Cellular Biology and AnatomyLouisiana State University Health Sciences CenterShreveportLA71130USA
| | - Jagadeeswaran Indumathy
- Department of PhysiologyJawaharlal Institute of Postgraduate Medical Education and ResearchPuducherryIndia
| | - Edward Glasscock
- Department of Biological SciencesSouthern Methodist UniversityDallasTX75275USA
| | - Leonidas Iasemidis
- Department of Biomedical EngineeringLouisiana Tech UniversityRustonLA71272USA
- Center for Biomedical Engineering and Rehabilitation ScienceLouisiana Tech UniversityRustonLA71272USA
| |
Collapse
|
16
|
Narasimhan S, Kundassery KB, Gupta K, Johnson GW, Wills KE, Goodale SE, Haas K, Rolston JD, Naftel RP, Morgan VL, Dawant BM, González HFJ, Englot DJ. Seizure-onset regions demonstrate high inward directed connectivity during resting-state: An SEEG study in focal epilepsy. Epilepsia 2020; 61:2534-2544. [PMID: 32944945 PMCID: PMC7899016 DOI: 10.1111/epi.16686] [Citation(s) in RCA: 50] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2020] [Revised: 08/15/2020] [Accepted: 08/17/2020] [Indexed: 01/06/2023]
Abstract
OBJECTIVE In patients with medically refractory focal epilepsy, stereotactic-electroencephalography (SEEG) can aid in localizing epileptogenic regions for surgical treatment. SEEG, however, requires long hospitalizations to record seizures, and ictal interpretation can be incomplete or inaccurate. Our recent work showed that non-directed resting-state analyses may identify brain regions as epileptogenic or uninvolved. Our present objective is to map epileptogenic networks in greater detail and more accurately identify seizure-onset regions using directed resting-state SEEG connectivity. METHODS In 25 patients with focal epilepsy who underwent SEEG, 2 minutes of resting-state, artifact-free, SEEG data were selected and functional connectivity was estimated. Using standard clinical interpretation, brain regions were classified into four categories: ictogenic, early propagation, irritative, or uninvolved. Three non-directed connectivity measures (mutual information [MI] strength, and imaginary coherence between and within regions) and four directed measures (partial directed coherence [PDC] and directed transfer function [DTF], inward and outward strength) were calculated. Logistic regression was used to generate a predictive model of ictogenicity. RESULTS Ictogenic regions had the highest and uninvolved regions had the lowest MI strength. Although both PDC and DTF inward strengths were highest in ictogenic regions, outward strengths did not differ among categories. A model incorporating directed and nondirected connectivity measures demonstrated an area under the receiver-operating characteristic (ROC) curve (AUC) of 0.88 in predicting ictogenicity of individual regions. The AUC of this model was 0.93 when restricted to patients with favorable postsurgical seizure outcomes. SIGNIFICANCE Directed connectivity measures may help identify epileptogenic networks without requiring ictal recordings. Greater inward but not outward connectivity in ictogenic regions at rest may represent broad inhibitory input to prevent seizure generation.
Collapse
Affiliation(s)
- Saramati Narasimhan
- Department of Neurological Surgery, Vanderbilt University Medical Center, Nashville, Tennessee
- Department of Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, Tennessee
- Department of Biomedical Engineering, Vanderbilt University, Nashville, Tennessee
| | - Keshav B. Kundassery
- Department of Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, Tennessee
| | - Kanupriya Gupta
- Department of Neurological Surgery, Vanderbilt University Medical Center, Nashville, Tennessee
- Department of Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, Tennessee
| | - Graham W. Johnson
- Department of Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, Tennessee
- Department of Biomedical Engineering, Vanderbilt University, Nashville, Tennessee
| | - Kristin E. Wills
- Department of Neurological Surgery, Vanderbilt University Medical Center, Nashville, Tennessee
- Department of Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, Tennessee
| | - Sarah E. Goodale
- Department of Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, Tennessee
- Department of Biomedical Engineering, Vanderbilt University, Nashville, Tennessee
| | - Kevin Haas
- Department of Neurology, Vanderbilt University Medical Center, Nashville, Tennessee
| | - John D. Rolston
- Department of Neurosurgery, University of Utah, Salt Lake City, Utah
| | - Robert P. Naftel
- Department of Neurological Surgery, Vanderbilt University Medical Center, Nashville, Tennessee
| | - Victoria L. Morgan
- Department of Neurological Surgery, Vanderbilt University Medical Center, Nashville, Tennessee
- Department of Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, Tennessee
- Department of Biomedical Engineering, Vanderbilt University, Nashville, Tennessee
- Department of Radiology and Radiological Sciences, Vanderbilt University Medical Center, Nashville, Tennessee
| | - Benoit M. Dawant
- Department of Biomedical Engineering, Vanderbilt University, Nashville, Tennessee
- Department of Electrical Engineering and Computer Science, Vanderbilt University, Nashville, Tennessee
| | - Hernán F. J. González
- Department of Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, Tennessee
- Department of Biomedical Engineering, Vanderbilt University, Nashville, Tennessee
| | - Dario J. Englot
- Department of Neurological Surgery, Vanderbilt University Medical Center, Nashville, Tennessee
- Department of Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, Tennessee
- Department of Biomedical Engineering, Vanderbilt University, Nashville, Tennessee
- Department of Radiology and Radiological Sciences, Vanderbilt University Medical Center, Nashville, Tennessee
- Department of Electrical Engineering and Computer Science, Vanderbilt University, Nashville, Tennessee
| |
Collapse
|
17
|
Abstract
Concise history of fascinating magnetoencephalography (MEG) technology and catalog of very selected milestone preclinical and clinical MEG studies are provided as the background. The focus is the societal context defining a journey of MEG to and through clinical practice and formation of the American Clinical MEG Society (ACMEGS). We aspired to provide an objective historic perspective and document contributions of many professionals while focusing on the role of ACMEGS in the growth and maturation of clinical MEG field. The ACMEGS was born (2006) out of inevitability to address two vital issues-fair reimbursement and proper clinical acceptance. A beacon of accountable MEG practice and utilization is now an expanding professional organization with the highest level of competence in practice of clinical MEG and clinical credibility. The ACMEGS facilitated a favorable disposition of insurances toward MEG in the United States by combining the national replication of the grassroots efforts and teaming up with the strategic partners-particularly the American Academy of Neurology (AAN), published two Position Statements (2009 and 2017), the world's only set of MEG Clinical Practice Guidelines (CPGs; 2011) and surveys of clinical MEG practice (2011 and 2020) and use (2020). In addition to the annual ACMEGS Course (2012), we directly engaged MEG practitioners through an Invitational Summit (2019). The Society remains focused on the improvements and expansion of clinical practice, education, clinical training, and constructive engagement of vendors in these issues and pivotal studies toward additional MEG indications. The ACMEGS not only had the critical role in the progress of Clinical MEG in the United States and beyond since 2006 but positioned itself as the field leader in the future.
Collapse
|
18
|
Ramaraju S, Wang Y, Sinha N, McEvoy AW, Miserocchi A, de Tisi J, Duncan JS, Rugg-Gunn F, Taylor PN. Removal of Interictal MEG-Derived Network Hubs Is Associated With Postoperative Seizure Freedom. Front Neurol 2020; 11:563847. [PMID: 33071948 PMCID: PMC7543719 DOI: 10.3389/fneur.2020.563847] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2020] [Accepted: 08/20/2020] [Indexed: 01/21/2023] Open
Abstract
Objective: To investigate whether MEG network connectivity was associated with epilepsy duration, to identify functional brain network hubs in patients with refractory focal epilepsy, and assess if their surgical removal was associated with post-operative seizure freedom. Methods: We studied 31 patients with drug refractory focal epilepsy who underwent resting state magnetoencephalography (MEG), and structural magnetic resonance imaging (MRI) as part of pre-surgical evaluation. Using the structural MRI, we generated 114 cortical regions of interest, performed surface reconstruction and MEG source localization. Representative source localized signals for each region were correlated with each other to generate a functional brain network. We repeated this procedure across three randomly chosen one-minute epochs. Network hubs were defined as those with the highest intra-hemispheric mean correlations. Post-operative MRI identified regions that were surgically removed. Results: Greater mean MEG network connectivity was associated with a longer duration of epilepsy. Patients who were seizure free after surgery had more hubs surgically removed than patients who were not seizure free (AUC = 0.76, p = 0.01) consistently across three randomly chosen time segments. Conclusion: Our results support a growing literature implicating network hub involvement in focal epilepsy, the removal of which by surgery is associated with greater chance of post-operative seizure freedom.
Collapse
Affiliation(s)
- Sriharsha Ramaraju
- Interdisciplinary Computing and Complex BioSystems Group, CNNP Lab, School of Computing, Newcastle University, Newcastle upon Tyne, United Kingdom
| | - Yujiang Wang
- Interdisciplinary Computing and Complex BioSystems Group, CNNP Lab, School of Computing, Newcastle University, Newcastle upon Tyne, United Kingdom.,Department of Clinical and Experimental Epilepsy, UCL Queen Square Institute of Neurology, London, United Kingdom.,Faculty of Medical Science, Newcastle University, Newcastle upon Tyne, United Kingdom
| | - Nishant Sinha
- Interdisciplinary Computing and Complex BioSystems Group, CNNP Lab, School of Computing, Newcastle University, Newcastle upon Tyne, United Kingdom.,Faculty of Medical Science, Newcastle University, Newcastle upon Tyne, United Kingdom
| | - Andrew W McEvoy
- Department of Clinical and Experimental Epilepsy, UCL Queen Square Institute of Neurology, London, United Kingdom
| | - Anna Miserocchi
- Department of Clinical and Experimental Epilepsy, UCL Queen Square Institute of Neurology, London, United Kingdom
| | - Jane de Tisi
- Department of Clinical and Experimental Epilepsy, UCL Queen Square Institute of Neurology, London, United Kingdom
| | - John S Duncan
- Department of Clinical and Experimental Epilepsy, UCL Queen Square Institute of Neurology, London, United Kingdom
| | - Fergus Rugg-Gunn
- Department of Clinical and Experimental Epilepsy, UCL Queen Square Institute of Neurology, London, United Kingdom
| | - Peter N Taylor
- Interdisciplinary Computing and Complex BioSystems Group, CNNP Lab, School of Computing, Newcastle University, Newcastle upon Tyne, United Kingdom.,Department of Clinical and Experimental Epilepsy, UCL Queen Square Institute of Neurology, London, United Kingdom.,Faculty of Medical Science, Newcastle University, Newcastle upon Tyne, United Kingdom
| |
Collapse
|
19
|
Abstract
Magnetoencephalography is the noninvasive measurement of miniscule magnetic fields produced by brain electrical currents, and is used most fruitfully to evaluate epilepsy patients. While other modalities infer brain function indirectly by measuring changes in blood flow, metabolism, and oxygenation, magnetoencephalography measures neuronal and synaptic function directly with submillisecond temporal resolution. The brain's magnetic field is recorded by neuromagnetometers surrounding the head in a helmet-shaped sensor array. Because magnetic signals are not distorted by anatomy, magnetoencephalography allows for a more accurate measurement and localization of brain activities than electroencephalography. Magnetoencephalography has become an indispensable part of the armamentarium at epilepsy centers.
Collapse
Affiliation(s)
- Richard C Burgess
- Epilepsy Center, Neurological Institute, The Cleveland Clinic Foundation, 9500 Euclid Avenue, Cleveland, OH 44195, USA.
| |
Collapse
|
20
|
He B, Astolfi L, Valdés-Sosa PA, Marinazzo D, Palva SO, Bénar CG, Michel CM, Koenig T. Electrophysiological Brain Connectivity: Theory and Implementation. IEEE Trans Biomed Eng 2019; 66:10.1109/TBME.2019.2913928. [PMID: 31071012 PMCID: PMC6834897 DOI: 10.1109/tbme.2019.2913928] [Citation(s) in RCA: 102] [Impact Index Per Article: 17.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
We review the theory and algorithms of electrophysiological brain connectivity analysis. This tutorial is aimed at providing an introduction to brain functional connectivity from electrophysiological signals, including electroencephalography (EEG), magnetoencephalography (MEG), electrocorticography (ECoG), stereoelectroencephalography (SEEG). Various connectivity estimators are discussed, and algorithms introduced. Important issues for estimating and mapping brain functional connectivity with electrophysiology are discussed.
Collapse
Affiliation(s)
- Bin He
- Department of Biomedical Engineering, Carnegie Mellon University, Pittsburgh, USA
| | - Laura Astolfi
- Department of Computer, Control and Management Engineering, University of Rome Sapienza, and with IRCCS Fondazione Santa Lucia, Rome, Italy
| | | | | | | | | | | | | |
Collapse
|
21
|
Burgess RC. Magnetoencephalography for localizing and characterizing the epileptic focus. HANDBOOK OF CLINICAL NEUROLOGY 2019; 160:203-214. [PMID: 31277848 DOI: 10.1016/b978-0-444-64032-1.00013-8] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/04/2022]
Abstract
Magnetoencephalography (MEG) is the noninvasive measurement of the miniscule magnetic fields produced by electrical currents flowing in the brain-the same neuroelectric activity that produces the EEG. MEG is one of several diagnostic tests employed in the evaluation of patients with epilepsy, but without the need to expose the patient to any potentially harmful agents. MEG is especially important in those being considered for epilepsy surgery, in whom accurate localization of the epileptic focus is paramount. While other modalities infer brain function indirectly by measuring changes in blood flow, metabolism, oxygenation, etc., MEG, as well as EEG, measures neuronal and synaptic function directly and, like EEG, MEG enjoys submillisecond temporal resolution. The measurement of magnetic fields provides information not only about the amplitude of the current but also its orientation. MEG picks up the magnetic field from neuromagnetometers surrounding the head in a helmet-shaped array of sensors. Clinical whole-head systems currently have 200-300 magnetic sensors, thereby offering very high resolution. The magnetic signals are not distorted by anatomy, because magnetic susceptibility is the same for all tissues, including the skull. Hence, MEG allows for a more accurate measurement and localization of brain activities than does EEG. Because one of its primary strengths is the ability to precisely localize electromagnetic activity within brain areas, MEG results are always coregistered to the patient's MRI. When combined in this way with structural imaging, it has been called magnetic source imaging (MSI), but MEG is properly understood as a clinical neurophysiologic diagnostic test. Signal processing and clinical interpretation in magnetoencephalography require sophisticated noise reduction and computerized mathematical modeling. Technological advances in these areas have brought MEG to the point where it is now part of routine clinical practice. MEG has become an indispensable part of the armamentarium at epilepsy centers where MEG laboratories are located, especially when patients are MRI-negative or where results of other structural and functional tests are not entirely concordant.
Collapse
Affiliation(s)
- Richard C Burgess
- Department of Neurology, Cleveland Clinic Foundation, Cleveland, OH, United States.
| |
Collapse
|
22
|
Hutson T, Pizarro D, Pati S, Iasemidis LD. Predictability and Resetting in a Case of Convulsive Status Epilepticus. Front Neurol 2018; 9:172. [PMID: 29623064 PMCID: PMC5874309 DOI: 10.3389/fneur.2018.00172] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2017] [Accepted: 03/06/2018] [Indexed: 11/13/2022] Open
Abstract
In this case study, we present evidence of resetting of brain dynamics following convulsive status epilepticus (SE) that was treated successfully with antiepileptic medications (AEDs). The measure of effective inflow (EI), a novel measure of network connectivity, was applied to the continuously recorded multichannel intracranial stereoelectroencephalographic (SEEG) signals before, during and after SE. Results from this analysis indicate trends of progressive reduction of EI over hours up to the onset of SE, mainly at sites of the epileptogenic focus with reversal of those trends upon successful treatment of SE by AEDs. The proposed analytical framework is promising for elucidation of the pathology of neuronal network dynamics that could lead to SE, evaluation of the efficacy of SE treatment strategies, as well as the development of biomarkers for susceptibility to SE.
Collapse
Affiliation(s)
- Timothy Hutson
- Department of Biomedical Engineering, Louisiana Tech University, Ruston, LA, United States
| | - Diana Pizarro
- Department of Neurology, University of Alabama at Birmingham, Birmingham, AL, United States
| | - Sandipan Pati
- Department of Neurology, University of Alabama at Birmingham, Birmingham, AL, United States
| | - Leon D Iasemidis
- Department of Biomedical Engineering, Louisiana Tech University, Ruston, LA, United States
| |
Collapse
|
23
|
Wang Q, Teng P, Luan G. Magnetoencephalography in Preoperative Epileptic Foci Localization: Enlightenment from Cognitive Studies. Front Comput Neurosci 2017; 11:58. [PMID: 28701945 PMCID: PMC5487414 DOI: 10.3389/fncom.2017.00058] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2017] [Accepted: 06/12/2017] [Indexed: 02/02/2023] Open
Abstract
Over 30% epileptic patients are refractory to medication, who are amenable to neurosurgical treatment. Non-invasive brain imaging technologies including video-electroencephalogram (EEG), magnetic resonance imaging (MRI), and magnetoencephalography (MEG) are widely used in presurgical assessment of epileptic patients. This review mainly discussed the current development of clinical MEG imaging as a diagnose approach, and its correlations with the golden standard intracranial electroencephalogram (iEEG). More importantly, this review discussed the possible applications of functional networks in preoperative epileptic foci localization in future studies.
Collapse
Affiliation(s)
- Qian Wang
- Beijing Key Laboratory of Epilepsy, Sanbo Brain Hospital, Capital Medical UniversityBeijing, China.,Department of Neurosurgery, Epilepsy Center, Sanbo Brain Hospital, Capital Medical UniversityBeijing, China
| | - Pengfei Teng
- Department of Neurosurgery, Epilepsy Center, Sanbo Brain Hospital, Capital Medical UniversityBeijing, China
| | - Guoming Luan
- Beijing Key Laboratory of Epilepsy, Sanbo Brain Hospital, Capital Medical UniversityBeijing, China.,Department of Neurosurgery, Epilepsy Center, Sanbo Brain Hospital, Capital Medical UniversityBeijing, China.,Beijing Institute for Brain Disorders, Capital Medical UniversityBeijing, China
| |
Collapse
|
24
|
Perescis MFJ, de Bruin N, Heijink L, Kruse C, Vinogradova L, Lüttjohann A, van Luijtelaar G, van Rijn CM. Cannabinoid antagonist SLV326 induces convulsive seizures and changes in the interictal EEG in rats. PLoS One 2017; 12:e0165363. [PMID: 28151935 PMCID: PMC5289424 DOI: 10.1371/journal.pone.0165363] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2016] [Accepted: 10/11/2016] [Indexed: 11/18/2022] Open
Abstract
Cannabinoid CB1 antagonists have been investigated for possible treatment of e.g. obesity-related disorders. However, clinical application was halted due to their symptoms of anxiety and depression. In addition to these adverse effects, we have shown earlier that chronic treatment with the CB1 antagonist rimonabant may induce EEG-confirmed convulsive seizures. In a regulatory repeat-dose toxicity study violent episodes of “muscle spasms” were observed in Wistar rats, daily dosed with the CB1 receptor antagonist SLV326 during 5 months. The aim of the present follow-up study was to investigate whether these violent movements were of an epileptic origin. In selected SLV326-treated and control animals, EEG and behavior were monitored for 24 hours. 25% of SLV326 treated animals showed 1 to 21 EEG-confirmed generalized convulsive seizures, whereas controls were seizure-free. The behavioral seizures were typical for a limbic origin. Moreover, interictal spikes were found in 38% of treated animals. The frequency spectrum of the interictal EEG of the treated rats showed a lower theta peak frequency, as well as lower gamma power compared to the controls. These frequency changes were state-dependent: they were only found during high locomotor activity. It is concluded that long term blockade of the endogenous cannabinoid system can provoke limbic seizures in otherwise healthy rats. Additionally, SLV326 alters the frequency spectrum of the EEG when rats are highly active, suggesting effects on complex behavior and cognition.
Collapse
Affiliation(s)
- Martin F. J. Perescis
- Donders Institute for Brain, Cognition and Behaviour, Radboud University, Nijmegen, The Netherlands
- HAS University of Applied Sciences, ‘s-Hertogenbosch, The Netherlands
- * E-mail:
| | - Natasja de Bruin
- Abbott Healthcare Products BV (formerly Solvay Pharmaceuticals), Weesp, The Netherlands
| | - Liesbeth Heijink
- Abbott Healthcare Products BV (formerly Solvay Pharmaceuticals), Weesp, The Netherlands
| | - Chris Kruse
- Abbott Healthcare Products BV (formerly Solvay Pharmaceuticals), Weesp, The Netherlands
| | - Lyudmila Vinogradova
- Institute of Higher Nervous Activity and Neurophysiology, Russian Academy of Sciences, Moscow, Russia
| | - Annika Lüttjohann
- Institut für Physiologie I, Westfälische Wilhelms Universität Münster, Münster, Germany
| | - Gilles van Luijtelaar
- Donders Institute for Brain, Cognition and Behaviour, Radboud University, Nijmegen, The Netherlands
| | - Clementina M. van Rijn
- Donders Institute for Brain, Cognition and Behaviour, Radboud University, Nijmegen, The Netherlands
| |
Collapse
|
25
|
Nissen IA, Stam CJ, Reijneveld JC, van Straaten IECW, Hendriks EJ, Baayen JC, De Witt Hamer PC, Idema S, Hillebrand A. Identifying the epileptogenic zone in interictal resting-state MEG source-space networks. Epilepsia 2016; 58:137-148. [PMID: 27888520 DOI: 10.1111/epi.13622] [Citation(s) in RCA: 59] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 11/01/2016] [Indexed: 01/07/2023]
Abstract
OBJECTIVE In one third of patients, seizures remain after epilepsy surgery, meaning that improved preoperative evaluation methods are needed to identify the epileptogenic zone. A potential framework for such a method is network theory, as it can be applied to noninvasive recordings, even in the absence of epileptiform activity. Our aim was to identify the epileptogenic zone on the basis of hub status of local brain areas in interictal magnetoencephalography (MEG) networks. METHODS Preoperative eyes-closed resting-state MEG recordings were retrospectively analyzed in 22 patients with refractory epilepsy, of whom 14 were seizure-free 1 year after surgery. Beamformer-based time series were reconstructed for 90 cortical and subcortical automated anatomic labeling (AAL) regions of interest (ROIs). Broadband functional connectivity was estimated using the phase lag index in artifact-free epochs without interictal epileptiform abnormalities. A minimum spanning tree was generated to represent the network, and the hub status of each ROI was calculated using betweenness centrality, which indicates the centrality of a node in a network. The correspondence of resection cavity to hub values was evaluated on four levels: resection cavity, lobar, hemisphere, and temporal versus extratemporal areas. RESULTS Hubs were localized within the resection cavity in 8 of 14 seizure-free patients and in zero of 8 patients who were not seizure-free (57% sensitivity, 100% specificity, 73% accuracy). Hubs were localized in the lobe of resection in 9 of 14 seizure-free patients and in zero of 8 patients who were not seizure-free (64% sensitivity, 100% specificity, 77% accuracy). For the other two levels, the true negatives are unknown; hence, only sensitivity could be determined: hubs coincided with both the resection hemisphere and the resection location (temporal versus extratemporal) in 11 of 14 seizure-free patients (79% sensitivity). SIGNIFICANCE Identifying hubs noninvasively before surgery is a valuable approach with the potential of indicating the epileptogenic zone in patients without interictal abnormalities.
Collapse
Affiliation(s)
- Ida A Nissen
- Department of Clinical Neurophysiology and MEG Center, VU University Medical Center, Amsterdam, The Netherlands
| | - Cornelis J Stam
- Department of Clinical Neurophysiology and MEG Center, VU University Medical Center, Amsterdam, The Netherlands
| | - Jaap C Reijneveld
- Brain Tumor Center Amsterdam & Department of Neurology, VU University Medical Center, Amsterdam, The Netherlands
| | - Ilse E C W van Straaten
- Department of Clinical Neurophysiology and MEG Center, VU University Medical Center, Amsterdam, The Netherlands
| | - Eef J Hendriks
- Department of Radiology & Nuclear Medicine, VU University Medical Center, Amsterdam, The Netherlands
| | - Johannes C Baayen
- Neurosurgical Center Amsterdam, VU University Medical Center, Amsterdam, The Netherlands
| | - Philip C De Witt Hamer
- Neurosurgical Center Amsterdam, VU University Medical Center, Amsterdam, The Netherlands
| | - Sander Idema
- Neurosurgical Center Amsterdam, VU University Medical Center, Amsterdam, The Netherlands
| | - Arjan Hillebrand
- Department of Clinical Neurophysiology and MEG Center, VU University Medical Center, Amsterdam, The Netherlands
| |
Collapse
|