1
|
Lahtinen J, Moura F, Samavaki M, Siltanen S, Pursiainen S. In silicostudy of the effects of cerebral circulation on source localization using a dynamical anatomical atlas of the human head. J Neural Eng 2023; 20. [PMID: 36808911 DOI: 10.1088/1741-2552/acbdc1] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/03/2022] [Accepted: 02/21/2023] [Indexed: 02/23/2023]
Abstract
Objective.This study focuses on the effects of dynamical vascular modeling on source localization errors in electroencephalography (EEG). Our aim of thisin silicostudy is to (a) find out the effects of cerebral circulation on the accuracy of EEG source localization estimates, and (b) evaluate its relevance with respect to measurement noise and interpatient variation.Approach.We employ a four-dimensional (3D + T) statistical atlas of the electrical properties of the human head with a cerebral circulation model to generate virtual patients with different cerebral circulatory conditions for EEG source localization analysis. As source reconstruction techniques, we use the linearly constraint minimum variance (LCMV) beamformer, standardized low-resolution brain electromagnetic tomography (sLORETA), and the dipole scan (DS).Main results.Results indicate that arterial blood flow affects source localization at different depths and with varying significance. The average flow rate plays an important role in source localization performance, while the pulsatility effects are very small. In cases where a personalized model of the head is available, blood circulation mismodeling causes localization errors, especially in the deep structures of the brain where the main cerebral arteries are located. When interpatient variations are considered, the results show differences up to 15 mm for sLORETA and LCMV beamformer and 10 mm for DS in the brainstem and entorhinal cortices regions. In regions far from the main arteries vessels, the discrepancies are smaller than 3 mm. When measurement noise is added and interpatient differences are considered in a deep dipolar source, the results indicate that the effects of conductivity mismatch are detectable even for moderate measurement noise. The signal-to-noise ratio limit for sLORETA and LCMV beamformer is 15 dB, while the limit is under 30 dB for DS.Significance.Localization of the brain activity via EEG constitutes an ill-posed inverse problem, where any modeling uncertainty, e.g. a slight amount of noise in the data or material parameter discrepancies, can lead to a significant deviation of the estimated activity, especially in the deep structures of the brain. Proper modeling of the conductivity distribution is necessary in order to obtain an appropriate source localization. In this study, we show that the conductivity of the deep brain structures is particularly impacted by blood flow-induced changes in conductivity because large arteries and veins access the brain through that region.
Collapse
Affiliation(s)
- Joonas Lahtinen
- Computing Sciences Unit, Faculty of Information Technology and Communication Sciences, Tampere University, Tampere, Finland
| | - Fernando Moura
- Department of Mathematics and Statistics, University of Helsinki, Helsinki, Finland.,Engineering, Modelling and Applied Social Sciences Center, Federal University of ABC, São Bernardo do Campo, São Paulo, Brazil
| | - Maryam Samavaki
- Computing Sciences Unit, Faculty of Information Technology and Communication Sciences, Tampere University, Tampere, Finland
| | - Samuli Siltanen
- Department of Mathematics and Statistics, University of Helsinki, Helsinki, Finland
| | - Sampsa Pursiainen
- Computing Sciences Unit, Faculty of Information Technology and Communication Sciences, Tampere University, Tampere, Finland
| |
Collapse
|
2
|
Ricci L, Matarrese M, Peters JM, Tamilia E, Madsen JR, Pearl PL, Papadelis C. Virtual implantation using conventional scalp EEG delineates seizure onset and predicts surgical outcome in children with epilepsy. Clin Neurophysiol 2022; 139:49-57. [PMID: 35526353 PMCID: PMC10026594 DOI: 10.1016/j.clinph.2022.04.009] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/09/2022] [Revised: 04/04/2022] [Accepted: 04/08/2022] [Indexed: 11/16/2022]
Abstract
OBJECTIVE Delineation of the seizure onset zone (SOZ) is required in children with drug resistant epilepsy (DRE) undergoing neurosurgery. Intracranial EEG (icEEG) serves as gold standard but has limitations. Here, we examine the utility of virtual implantation with electrical source imaging (ESI) on ictal scalp EEG for mapping the SOZ and predict surgical outcome. METHODS We retrospectively analyzed EEG data from 35 children with DRE who underwent surgery and dichotomized into seizure-free (SF) and non-seizure-free (NSF). We estimated virtual sensors (VSs) at brain locations that matched icEEG implantation and compared ictal patterns at VSs vs icEEG. We calculated the agreement between VSs SOZ and clinically defined SOZ and built receiver operating characteristic (ROC) curves to test whether it predicted outcome. RESULTS Twenty-one patients were SF after surgery. Moderate agreement between virtual and icEEG patterns was observed (kappa = 0.45, p < 0.001). Virtual SOZ agreement with clinically defined SOZ was higher in SF vs NSF patients (66.6% vs 41.6%, p = 0.01). Anatomical concordance of virtual SOZ with clinically defined SOZ predicted outcome (AUC = 0.73; 95% CI: 0.57-0.89; sensitivity = 66.7%; specificity = 78.6%; accuracy = 71.4%). CONCLUSIONS Virtual implantation on ictal scalp EEG can approximate the SOZ and predict outcome. SIGNIFICANCE SOZ mapping with VSs may contribute to tailoring icEEG implantation and predict outcome.
Collapse
Affiliation(s)
- Lorenzo Ricci
- Unit of Neurology, Neurophysiology, Neurobiology, Department of Medicine, University Campus Bio-Medico of Rome, Rome, Italy
| | - Margherita Matarrese
- Unit of Non-Linear Physics and Mathematical Modelling, Engineering Department, University Campus Bio-Medico of Rome, Rome, Italy; Jane and John Justin Neurosciences Center, Cook Children's Health Care System, Fort Worth, TX, USA; Department of Bioengineering, University of Texas at Arlington, Arlington, TX, USA
| | - Jurriaan M Peters
- Division of Epilepsy and Clinical Neurophysiology, Department of Neurology, Boston Children's Hospital, Harvard Medical School, Boston, MA, USA
| | - Eleonora Tamilia
- Division of Newborn Medicine, Department of Medicine, Boston Children's Hospital, Harvard Medical School, Boston, MA, USA
| | - Joseph R Madsen
- Division of Epilepsy Surgery, Department of Neurosurgery, Boston Children's Hospital, Harvard Medical School, Boston, MA, USA
| | - Phillip L Pearl
- Division of Epilepsy and Clinical Neurophysiology, Department of Neurology, Boston Children's Hospital, Harvard Medical School, Boston, MA, USA
| | - Christos Papadelis
- Jane and John Justin Neurosciences Center, Cook Children's Health Care System, Fort Worth, TX, USA; Department of Bioengineering, University of Texas at Arlington, Arlington, TX, USA; School of Medicine, Texas Christian University, Fort Worth, TX, USA.
| |
Collapse
|
3
|
Validating EEG, MEG and Combined MEG and EEG Beamforming for an Estimation of the Epileptogenic Zone in Focal Cortical Dysplasia. Brain Sci 2022; 12:brainsci12010114. [PMID: 35053857 PMCID: PMC8796031 DOI: 10.3390/brainsci12010114] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2021] [Revised: 01/04/2022] [Accepted: 01/06/2022] [Indexed: 02/04/2023] Open
Abstract
MEG and EEG source analysis is frequently used for the presurgical evaluation of pharmacoresistant epilepsy patients. The source localization of the epileptogenic zone depends, among other aspects, on the selected inverse and forward approaches and their respective parameter choices. In this validation study, we compare the standard dipole scanning method with two beamformer approaches for the inverse problem, and we investigate the influence of the covariance estimation method and the strength of regularization on the localization performance for EEG, MEG, and combined EEG and MEG. For forward modelling, we investigate the difference between calibrated six-compartment and standard three-compartment head modelling. In a retrospective study, two patients with focal epilepsy due to focal cortical dysplasia type IIb and seizure freedom following lesionectomy or radiofrequency-guided thermocoagulation (RFTC) used the distance of the localization of interictal epileptic spikes to the resection cavity resp. RFTC lesion as reference for good localization. We found that beamformer localization can be sensitive to the choice of the regularization parameter, which has to be individually optimized. Estimation of the covariance matrix with averaged spike data yielded more robust results across the modalities. MEG was the dominant modality and provided a good localization in one case, while it was EEG for the other. When combining the modalities, the good results of the dominant modality were mostly not spoiled by the weaker modality. For appropriate regularization parameter choices, the beamformer localized better than the standard dipole scan. Compared to the importance of an appropriate regularization, the sensitivity of the localization to the head modelling was smaller, due to similar skull conductivity modelling and the fixed source space without orientation constraint.
Collapse
|
4
|
The Value of Source Localization for Clinical Magnetoencephalography: Beyond the Equivalent Current Dipole. J Clin Neurophysiol 2020; 37:537-544. [DOI: 10.1097/wnp.0000000000000487] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022] Open
|
5
|
Comparison of beamformer implementations for MEG source localization. Neuroimage 2020; 216:116797. [PMID: 32278091 PMCID: PMC7322560 DOI: 10.1016/j.neuroimage.2020.116797] [Citation(s) in RCA: 31] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2019] [Revised: 02/18/2020] [Accepted: 03/31/2020] [Indexed: 01/09/2023] Open
Abstract
Beamformers are applied for estimating spatiotemporal characteristics of neuronal sources underlying measured MEG/EEG signals. Several MEG analysis toolboxes include an implementation of a linearly constrained minimum-variance (LCMV) beamformer. However, differences in implementations and in their results complicate the selection and application of beamformers and may hinder their wider adoption in research and clinical use. Additionally, combinations of different MEG sensor types (such as magnetometers and planar gradiometers) and application of preprocessing methods for interference suppression, such as signal space separation (SSS), can affect the results in different ways for different implementations. So far, a systematic evaluation of the different implementations has not been performed. Here, we compared the localization performance of the LCMV beamformer pipelines in four widely used open-source toolboxes (MNE-Python, FieldTrip, DAiSS (SPM12), and Brainstorm) using datasets both with and without SSS interference suppression. We analyzed MEG data that were i) simulated, ii) recorded from a static and moving phantom, and iii) recorded from a healthy volunteer receiving auditory, visual, and somatosensory stimulation. We also investigated the effects of SSS and the combination of the magnetometer and gradiometer signals. We quantified how localization error and point-spread volume vary with the signal-to-noise ratio (SNR) in all four toolboxes. When applied carefully to MEG data with a typical SNR (3–15 dB), all four toolboxes localized the sources reliably; however, they differed in their sensitivity to preprocessing parameters. As expected, localizations were highly unreliable at very low SNR, but we found high localization error also at very high SNRs for the first three toolboxes while Brainstorm showed greater robustness but with lower spatial resolution. We also found that the SNR improvement offered by SSS led to more accurate localization. Different beamformer implementations are reported to sometimes yield differing source estimates for the same MEG data. We compared beamformers in four major open-source MEG analysis toolboxes. All toolboxes provide consistent and accurate results with 3–15-dB input SNR. However, localization errors are high at very high input SNR for the tested scalar beamformers. We discuss the critical differences between the implementations.
Collapse
|
6
|
Gofshteyn JS, Le T, Kessler S, Kamens R, Carr C, Gaetz W, Bloy L, Roberts TPL, Schwartz ES, Marsh ED. Synthetic aperture magnetometry and excess kurtosis mapping of Magnetoencephalography (MEG) is predictive of epilepsy surgical outcome in a large pediatric cohort. Epilepsy Res 2019; 155:106151. [PMID: 31247475 PMCID: PMC6699633 DOI: 10.1016/j.eplepsyres.2019.106151] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2019] [Revised: 05/23/2019] [Accepted: 06/09/2019] [Indexed: 12/23/2022]
Abstract
OBJECTIVE Resective surgery is the most effective treatment option for patients with refractory epilepsy; however identification of patients who will benefit from epilepsy surgery remains challenging. Synthetic aperture magnetometry and excess kurtosis mapping (SAM(g2)) of magnetoencephalography (MEG) is a non-invasive tool that warrants further examination in the pediatric epilepsy population. Here, we examined the utility of MEG with SAM(g2) to determine if MEG epileptiform foci correlates with surgical outcome and to develop a predictive model incorporating MEG information to best assess likelihood of seizure improvement/freedom from resective surgery. METHODS 564 subjects who had MEG at the Children's Hospital of Philadelphia between 2010-2015 were screened. Clinical epilepsy history and prior electrographic records were extracted and reviewed and correlated with MEG findings. MEG assessments were made by both a neurologist and neuroradiologist. Predictive models were developed to assess the utility of MEG in determining Engel class at one year and five years after resective epilepsy surgery. RESULTS The number of MEG spike foci was highly associated with Engel class outcome at both one year and five years; however, using MEG data in isolation was not significantly predictive of 5 year surgical outcome. When combined with clinical factors; scalp EEG (single ictal onset zone), MRI (lesional or not), age and sex in a logistic regression model MEG foci was significant for Engel class outcome at both 1 year (p = 0.03) and 5 years (0.02). The percent correctly classified for Engel class at one year was 78.43% and the positive predictive value was 71.43. SIGNIFICANCE MEG using SAM(g2) analysis in an important non-invasive tool in the identification of those patients who will benefit most from surgery. Integrating MEG data analysis into pre-surgical evaluation can help to predict epilepsy outcome after resective surgery in the pediatric population if utilized with skilled interpretation.
Collapse
Affiliation(s)
- J S Gofshteyn
- Division of Pediatric Neurology, Department of Pediatrics, Weill Cornell Medicine, New York, NY, United States; New-York Presbyterian Hospital, New York, NY, United States
| | - T Le
- Division of Pediatric Neurology, The Children's Hospital of Philadelphia, United States
| | - S Kessler
- Division of Pediatric Neurology, The Children's Hospital of Philadelphia, United States; Departments of Neurology and Pediatrics, Perelman School of Medicine at the University of Pennsylvania, United States
| | - R Kamens
- Division of Pediatric Neurology, The Children's Hospital of Philadelphia, United States
| | - C Carr
- Division of Pediatric Neurology, The Children's Hospital of Philadelphia, United States
| | - W Gaetz
- Division of Neuroradiology, Department of Radiology, The Children's Hospital of Philadelphia, United States; Department of Radiology, Perelman School of Medicine at the University of Pennsylvania, United States
| | - L Bloy
- Division of Neuroradiology, Department of Radiology, The Children's Hospital of Philadelphia, United States; Department of Radiology, Perelman School of Medicine at the University of Pennsylvania, United States
| | - T P L Roberts
- Division of Neuroradiology, Department of Radiology, The Children's Hospital of Philadelphia, United States; Department of Radiology, Perelman School of Medicine at the University of Pennsylvania, United States
| | - E S Schwartz
- Division of Neuroradiology, Department of Radiology, The Children's Hospital of Philadelphia, United States; Department of Radiology, Perelman School of Medicine at the University of Pennsylvania, United States
| | - E D Marsh
- Division of Pediatric Neurology, The Children's Hospital of Philadelphia, United States; Departments of Neurology and Pediatrics, Perelman School of Medicine at the University of Pennsylvania, United States.
| |
Collapse
|
7
|
Jobst C, Ferrari P, Isabella S, Cheyne D. BrainWave: A Matlab Toolbox for Beamformer Source Analysis of MEG Data. Front Neurosci 2018; 12:587. [PMID: 30186107 PMCID: PMC6113377 DOI: 10.3389/fnins.2018.00587] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2018] [Accepted: 08/06/2018] [Indexed: 01/13/2023] Open
Abstract
BrainWave is an easy-to-use Matlab toolbox for the analysis of magnetoencephalography data. It provides a graphical user interface for performing minimum-variance beamforming analysis with rapid and interactive visualization of evoked and induced brain activity. This article provides an overview of the main features of BrainWave with a step-by-step demonstration of how to proceed from raw experimental data to group source images and time series analyses. This includes data selection and pre-processing, magnetic resonance image co-registration and normalization procedures, and the generation of volumetric (whole-brain) or cortical surface based source images, and corresponding source time series as virtual sensor waveforms and their time-frequency representations. We illustrate these steps using example data from a recently published study on response inhibition (Isabella et al., 2015) using the sustained attention to response task paradigm in 12 healthy adult participants. In this task participants were required to press a button with their right index finger to a rapidly presented series of numerical digits and withhold their response to an infrequently presented target digit. This paradigm elicited movement-locked brain responses, as well as task-related modulation of brain rhythmic activity in different frequency bands (e.g., theta, beta, and gamma), and is used to illustrate two different types of source reconstruction implemented in the BrainWave toolbox: (1) event-related beamforming of averaged brain responses and (2) beamformer analysis of modulation of rhythmic brain activity using the synthetic aperture magnetometry algorithm. We also demonstrate the ability to generate group contrast images between different response types, using the example of frontal theta activation patterns during error responses (failure to withhold on target trials). BrainWave is free academic software available for download at http://cheynelab.utoronto.ca/brainwave along with supporting software and documentation. The development of the BrainWave toolbox was supported by grants from the Canadian Institutes of Health Research, the National Research and Engineering Research Council of Canada, and the Ontario Brain Institute.
Collapse
Affiliation(s)
- Cecilia Jobst
- Program in Neurosciences and Mental Health, The Hospital for Sick Children, Toronto, ON, Canada
| | - Paul Ferrari
- MEG Laboratory, Dell Children's Medical Centre of Central Texas, Austin, TX, United States
| | - Silvia Isabella
- Program in Neurosciences and Mental Health, The Hospital for Sick Children, Toronto, ON, Canada
| | - Douglas Cheyne
- Program in Neurosciences and Mental Health, The Hospital for Sick Children, Toronto, ON, Canada.,Department of Medical Imaging, University of Toronto, Toronto, ON, Canada
| |
Collapse
|
8
|
Magnetoencephalography: Clinical and Research Practices. Brain Sci 2018; 8:brainsci8080157. [PMID: 30126121 PMCID: PMC6120049 DOI: 10.3390/brainsci8080157] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2018] [Revised: 08/07/2018] [Accepted: 08/11/2018] [Indexed: 11/25/2022] Open
Abstract
Magnetoencephalography (MEG) is a neurophysiological technique that detects the magnetic fields associated with brain activity. Synthetic aperture magnetometry (SAM), a MEG magnetic source imaging technique, can be used to construct both detailed maps of global brain activity as well as virtual electrode signals, which provide information that is similar to invasive electrode recordings. This innovative approach has demonstrated utility in both clinical and research settings. For individuals with epilepsy, MEG provides valuable, nonredundant information. MEG accurately localizes the irritative zone associated with interictal spikes, often detecting epileptiform activity other methods cannot, and may give localizing information when other methods fail. These capabilities potentially greatly increase the population eligible for epilepsy surgery and improve planning for those undergoing surgery. MEG methods can be readily adapted to research settings, allowing noninvasive assessment of whole brain neurophysiological activity, with a theoretical spatial range down to submillimeter voxels, and in both humans and nonhuman primates. The combination of clinical and research activities with MEG offers a unique opportunity to advance translational research from bench to bedside and back.
Collapse
|
9
|
Juárez-Martinez EL, Nissen IA, Idema S, Velis DN, Hillebrand A, Stam CJ, van Straaten ECW. Virtual localization of the seizure onset zone: Using non-invasive MEG virtual electrodes at stereo-EEG electrode locations in refractory epilepsy patients. Neuroimage Clin 2018; 19:758-766. [PMID: 30009129 PMCID: PMC6041424 DOI: 10.1016/j.nicl.2018.06.001] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2018] [Revised: 05/22/2018] [Accepted: 06/01/2018] [Indexed: 12/11/2022]
Abstract
In some patients with medically refractory epilepsy, EEG with intracerebrally placed electrodes (stereo-electroencephalography, SEEG) is needed to locate the seizure onset zone (SOZ) for successful epilepsy surgery. SEEG has limitations and entails risk of complications because of its invasive character. Non-invasive magnetoencephalography virtual electrodes (MEG-VEs) may overcome SEEG limitations and optimize electrode placement making SOZ localization safer. Our purpose was to assess whether interictal activity measured by MEG-VEs and SEEG at identical anatomical locations were comparable, and whether MEG-VEs activity properties could determine the location of a later resected brain area (RA) as an approximation of the SOZ. We analyzed data from nine patients who underwent MEG and SEEG evaluation, and surgery for medically refractory epilepsy. MEG activity was retrospectively reconstructed using beamforming to obtain VEs at the anatomical locations corresponding to those of SEEG electrodes. Spectral, functional connectivity and functional network properties were obtained for both, MEG-VEs and SEEG time series, and their correlation and reliability were established. Based on these properties, the approximation of the SOZ was characterized by the differences between RA and non-RA (NRA). We found significant positive correlation and reliability between MEG-VEs and SEEG spectral measures (particularly in delta [0.5-4 Hz], alpha2 [10-13 Hz], and beta [13-30 Hz] bands) and broadband functional connectivity. Both modalities showed significantly slower activity and a tendency towards increased broadband functional connectivity in the RA compared to the NRA. Our findings show that spectral and functional connectivity properties of non-invasively obtained MEG-VEs match those of invasive SEEG recordings, and can characterize the SOZ. This suggests that MEG-VEs might be used for optimal SEEG planning and fewer depth electrode implantations, making the localization of the SOZ safer and more successful.
Collapse
Affiliation(s)
| | - Ida A Nissen
- Department of Neurology and Clinical Neurophysiology, Amsterdam, the Netherlands
| | - Sander Idema
- Department of Neurosurgery, VU University Medical Center, Amsterdam, the Netherlands
| | - Demetrios N Velis
- Department of Neurology and Clinical Neurophysiology, Amsterdam, the Netherlands
| | - Arjan Hillebrand
- Department of Neurology and Clinical Neurophysiology, Amsterdam, the Netherlands
| | - Cornelis J Stam
- Department of Neurology and Clinical Neurophysiology, Amsterdam, the Netherlands
| | | |
Collapse
|
10
|
Magnetoencephalographic Spike Analysis in Patients With Focal Cortical Dysplasia: What Defines a "Dipole Cluster"? Pediatr Neurol 2018; 83:25-31. [PMID: 29685607 PMCID: PMC5988951 DOI: 10.1016/j.pediatrneurol.2018.03.004] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 12/29/2017] [Accepted: 03/09/2018] [Indexed: 11/22/2022]
Abstract
BACKGROUND The purpose of this study is to clarify the source distribution patterns of magnetoencephalographic spikes correlated with postsurgical seizure-free outcome in pediatric patients with focal cortical dysplasia. PATIENTS AND METHODS Thirty-two patients with pathologically confirmed focal cortical dysplasia were divided into seizure-free and seizure-persistent groups according to their surgical outcomes based on Engel classification. In each patient, presurgical magnetoencephalography was reviewed. Dipole sources of magnetoencephalographic spikes were calculated according to a single dipole model. We obtained the following quantitative indices for evaluating dipole distribution: maximum distance over all pairs of dipoles, standard deviation of the distances between each dipole and the mean coordinate of all dipoles, average nearest neighbor distance, the rate of dipoles located within 10, 20, and 30 mm from the mean coordinate, and the rate of dipoles included in the resection. These indices were compared between the two patient groups. RESULTS Average nearest neighbor distance was significantly smaller in the seizure-free group than in the seizure-persistent group (P = 0.008). The rates of dipoles located within 10, 20, and 30 mm from the mean coordinate were significantly higher in the seizure-free group (P = 0.001, 0.001, 0.005, respectively). The maximum distance, standard deviation, and resection rate of dipoles did not show a significant difference between the two groups. CONCLUSIONS A spatially restricted dipole distribution of magnetoencephalographic spikes is correlated with postsurgical seizure-free outcomes in patients with focal cortical dysplasia. The distribution can be assessed by quantitative indices that are clinically useful in the presurgical evaluation of these patients.
Collapse
|
11
|
Neugebauer F, Möddel G, Rampp S, Burger M, Wolters CH. The Effect of Head Model Simplification on Beamformer Source Localization. Front Neurosci 2017; 11:625. [PMID: 29209157 PMCID: PMC5701642 DOI: 10.3389/fnins.2017.00625] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2017] [Accepted: 10/26/2017] [Indexed: 11/13/2022] Open
Abstract
Beamformers are a widely-used tool in brain analysis with magnetoencephalography (MEG) and electroencephalography (EEG). For the construction of the beamformer filters realistic head volume conductor modeling is necessary for accurately computing the EEG and MEG leadfields, i.e., for solving the EEG and MEG forward problem. In this work, we investigate the influence of including realistic head tissue compartments into a finite element method (FEM) model on the beamformer's localization ability. Specifically, we investigate the effect of including cerebrospinal fluid, gray matter, and white matter distinction, as well as segmenting the skull bone into compacta and spongiosa, and modeling white matter anisotropy. We simulate an interictal epileptic measurement with white sensor noise. Beamformer filters are constructed with unit gain, unit array gain, and unit noise gain constraint. Beamformer source positions are determined by evaluating power and excess sample kurtosis (g2) of the source-waveforms at all source space nodes. For both modalities, we see a strong effect of modeling the cerebrospinal fluid and white and gray matter. Depending on the source position, both effects can each be in the magnitude of centimeters, rendering their modeling necessary for successful localization. Precise skull modeling mainly effected the EEG up to a few millimeters, while both modalities could profit from modeling white matter anisotropy to a smaller extent of 5-10 mm. The unit noise gain or neural activity index beamformer behaves similarly to the array gain beamformer when noise strength is sufficiently high. Variance localization seems more robust against modeling errors than kurtosis.
Collapse
Affiliation(s)
- Frank Neugebauer
- Institute for Biomagnetism und Biosignalanalysis, University of Münster, Münster, Germany
| | - Gabriel Möddel
- Department of Sleep Medicine and Neuromuscular Disorders, Epilepsy Center Münster-Osnabrück, University of Münster, Münster, Germany
| | - Stefan Rampp
- Department of Neurosurgery, University Hospital Erlangen, Erlangen, Germany
| | - Martin Burger
- Institute for Computational and Applied Mathematics, University of Münster, Münster, Germany
| | - Carsten H. Wolters
- Institute for Biomagnetism und Biosignalanalysis, University of Münster, Münster, Germany
| |
Collapse
|
12
|
Hillebrand A, Nissen IA, Ris-Hilgersom I, Sijsma NCG, Ronner HE, van Dijk BW, Stam CJ. Detecting epileptiform activity from deeper brain regions in spatially filtered MEG data. Clin Neurophysiol 2016; 127:2766-2769. [PMID: 27417050 DOI: 10.1016/j.clinph.2016.05.272] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/22/2014] [Revised: 05/19/2016] [Accepted: 05/19/2016] [Indexed: 10/21/2022]
Affiliation(s)
- A Hillebrand
- Department of Clinical Neurophysiology and Magnetoencephalography Center, VU University Medical Center, Amsterdam, The Netherlands.
| | - I A Nissen
- Department of Clinical Neurophysiology and Magnetoencephalography Center, VU University Medical Center, Amsterdam, The Netherlands
| | - I Ris-Hilgersom
- Department of Clinical Neurophysiology and Magnetoencephalography Center, VU University Medical Center, Amsterdam, The Netherlands
| | - N C G Sijsma
- Department of Clinical Neurophysiology and Magnetoencephalography Center, VU University Medical Center, Amsterdam, The Netherlands
| | - H E Ronner
- Department of Clinical Neurophysiology and Magnetoencephalography Center, VU University Medical Center, Amsterdam, The Netherlands
| | - B W van Dijk
- Department of Clinical Neurophysiology and Magnetoencephalography Center, VU University Medical Center, Amsterdam, The Netherlands; Department of Physics and Medical Technology, VU University Medical Center, Amsterdam, The Netherlands
| | - C J Stam
- Department of Clinical Neurophysiology and Magnetoencephalography Center, VU University Medical Center, Amsterdam, The Netherlands
| |
Collapse
|
13
|
Wennberg R, Cheyne D. Reliability of MEG source imaging of anterior temporal spikes: analysis of an intracranially characterized spike focus. Clin Neurophysiol 2013; 125:903-18. [PMID: 24210513 DOI: 10.1016/j.clinph.2013.08.032] [Citation(s) in RCA: 40] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2013] [Revised: 07/28/2013] [Accepted: 08/21/2013] [Indexed: 11/19/2022]
Abstract
OBJECTIVE To assess the reliability of MEG source imaging (MSI) of anterior temporal spikes through detailed analysis of the localization and orientation of source solutions obtained for a large number of spikes that were separately confirmed by intracranial EEG to be focally generated within a single, well-characterized spike focus. METHODS MSI was performed on 64 identical right anterior temporal spikes from an anterolateral temporal neocortical spike focus. The effects of different volume conductors (sphere and realistic head model), removal of noise with low frequency filters (LFFs) and averaging multiple spikes were assessed in terms of the reliability of the source solutions. RESULTS MSI of single spikes resulted in scattered dipole source solutions that showed reasonable reliability for localization at the lobar level, but only for solutions with a goodness-of-fit exceeding 80% using a LFF of 3 Hz. Reliability at a finer level of intralobar localization was limited. Spike averaging significantly improved the reliability of source solutions and averaging 8 or more spikes reduced dependency on goodness-of-fit and data filtering. CONCLUSIONS MSI performed on topographically identical individual spikes from an intracranially defined classical anterior temporal lobe spike focus was limited by low reliability (i.e., scattered source solutions) in terms of fine, sublobar localization within the ipsilateral temporal lobe. Spike averaging significantly improved reliability. SIGNIFICANCE MSI performed on individual anterior temporal spikes is limited by low reliability. Reduction of background noise through spike averaging significantly improves the reliability of MSI solutions.
Collapse
Affiliation(s)
- Richard Wennberg
- Krembil Neuroscience Centre, Division of Neurology, Toronto Western Hospital, University of Toronto, 399 Bathurst Street, Toronto, ON M5T 2S8, Canada.
| | - Douglas Cheyne
- Program in Neurosciences and Mental Health, Research Institute, The Hospital for Sick Children, University of Toronto, Toronto, ON M5G 1X8, Canada
| |
Collapse
|