Detection of short-lasting and ictal spike-and-wave discharges in around-the-ears EEG recordings in children with absence epilepsy.

PURPOSE: Long-term ambulatory EEG recordings can improve the monitoring of absence epilepsy in children, but signal quality and increased review workload are a concern. We evaluated the feasibility of around-the-ears EEG arrays (cEEGrids) to capture 3-Hz short-lasting and ictal spike-and-wave discharges and assessed the performance of automated detection software in cEEGrids data. We compared patterns of bilateral synchronisation between short-lasting and ictal spike-and-wave discharges. METHODS: We recruited children with suspected generalised epilepsy undergoing routine video-EEG monitoring and performed simultaneous cEEGrids recordings. We used ASSYST software to detect short-lasting 3-Hz spike-and-wave discharges (1-3 s) and ictal spike-and-wave discharges in the cEEGrids data. We assessed data quality and sensitivity of cEEGrids for spike-and-wave discharges in routine EEG. We determined the sensitivity and false detection rate for automated spike-and-wave discharge detection in cEEGrids data. We compared bihemispheric synchrony across the onset of short-lasting and ictal spike-and-wave discharges using the mean phase coherence in the 2-4 Hz frequency band. RESULTS: We included nine children with absence epilepsy (median age = 11 y, range 8-15 y, nine females) and recorded 4 h and 27 min of cEEGrids data. The recordings from seven participants were suitable for quantitative analysis, containing 82 spike-and-wave discharges. The cEEGrids captured 58 % of all spike-and-wave discharges (median individual sensitivity: 100 %, range: 47-100 %). ASSYST detected 82 % of all spike-and-wave discharges (median: 100 %, range: 41-100 %) with a false detection rate of 48/h (median: 6/h, range: 0-154/h). The mean phase coherence significantly increased during short-lasting and ictal spike-and-wave discharges in the 500-ms pre-onset to 1-s post-onset interval. CONCLUSIONS: cEEGrids are of variable quality for monitoring spike-and-wave discharges in children with absence epilepsy. ASSYST could facilitate the detection of short-lasting and ictal spike-and-wave discharges with clear periodic structures but with low specificity. A similar course of bihemispheric synchrony between short-lasting and ictal spike-and-wave discharges indicates that cortico-thalamic driving may be relevant for both types of spike-and-wave discharges.

Sensitivity of magnetoencephalography as a diagnostic tool for epilepsy: a prospective study.

The diagnostic process for epilepsy can be lengthy and stressful, which may delay the start of treatment. The objective of this study was to determine the benefit of routine magnetoencephalography (MEG) with regard to diagnostic gain, compared to routine electroencephalography (EEG), EEG following sleep deprivation (EEGsd), and 24-hour EEG. In this prospective study, patients were included from two centres (Academic Centre for Epileptology Kempenhaeghe, Heeze and Elisabeth-Twee Steden Hospital, Tilburg) and MEG recording took place at a single centre (Amsterdam University Medical Centre, Vrije Universiteit Amsterdam) in The Netherlands. Consecutively referred patients from peripheral hospitals were included between August 2013 and March 2016. Patients were offered routine MEG in addition to EEG examination and MRI for the diagnosis of epilepsy. The final clinical diagnosis was based on all available clinical data and test results at the end of the diagnostic process. Sensitivity, specificity, and positive and negative predictive values were calculated for routine EEG, routine EEG plus additional EEG and MEG. In addition, diagnostic gain associated with MEG, relative to the other modalities, was calculated. Secondary outcome was congruence of localization of epileptiform discharges between MEG and MRI or final clinical diagnosis. Based on a cohort of 138 patients, sensitivity and specificity was shown to be 31.6% and 78.4% for routine MEG, 31.6% and 100% for routine EEG, and 52.6% and 97.3% for routine EEG plus additional EEG, respectively. Routine MEG demonstrated a diagnostic gain of 16.8% compared to routine EEG and 9.5% compared to routine EEG plus additional EEG. In 35.7% of patients with a lesion on MRI that was consistent with the final clinical diagnosis, MEG showed epileptiform discharges in the same area. Routine MEG may provide additional value during the initial diagnosis of epilepsy.

A network approach to investigate the bi-hemispheric synchrony in absence epilepsy.

OBJECTIVE: Our objective was to unravel the dynamics underlying spike-and-wave discharges (SWDs) characteristic for childhood absence epilepsy. METHODS: SWDs were recorded for a cohort of 28 children using magnetoencephalography. Non-linear association analyses and a graph theoretical metric of local connectedness (LoC) were utilized in a sliding window starting one s before till four s after ictal onset. RESULTS: A focal pattern of bilateral frontal and parietal areas with high LoC during the spikes alternated by generalized patterns during the waves was found for all children studied during generalization of the SWDs. In the interval preceding the generalization a focal parietal region was most often (16/28) encountered and less often an occipital (4/28), temporal (5/28) or frontal (3/28) region. 55% of the children with a parietal/occipital focal onset became seizure free after the administration of two anti-epileptic drugs, and only 12.5% with a temporal/frontal focal onset. CONCLUSIONS: The transition from the interictal to the ictal state is for some of the children characterized by dominant LoC at either the parietal/occipital and for others at the frontal/temporal region. SIGNIFICANCE: The focal onset of the SWDs varies in location among the children with a clinical similar profile, who, however, seemingly are differing with regard to seizure control.

Burden of disease of people with epilepsy during an optimized diagnostic trajectory: costs and quality of life.

BACKGROUND: Diagnosing epilepsy can be lengthy and stressful, potentially leading to increased use of healthcare resources and a reduction in quality of life. AIM: This study aims to determine cost and quality of life before and after an optimized diagnostic procedure for people suspected of having epilepsy from a societal perspective with a follow-up of 12 months. In addition, this study aims to differentiate between people diagnosed with epilepsy during the follow-up of the study and the people who are diagnosed as not having epilepsy or for whom diagnosis is still uncertain. METHODS: A questionnaire regarding the use of healthcare resources was used accompanied by the EQ-5D-3 L. Multiple imputations by chained equations with predictive mean matching was used to account for missing data. To investigate the uncertainty of the results, non-parametric bootstrapped (1000 times) was used. RESULTS: In total, 116 people were included in the study. Total average costs per patient made in the previous 3 months had decreased from €4594 before the optimized diagnostic trajectory to €2609 in the 12 months after the optimized diagnostic trajectory. Healthcare costs were the largest expense group (52-66%) and had decreased significantly from baseline measurement to 12 months after baseline (€2395 vs €1581). Productivity costs had decreased from €1367 to €442 per 3 months. Total annual costs were similar between people diagnosed with epilepsy during the follow-up of the study and the people who are diagnosed as not having epilepsy or for whom diagnosis is still uncertain. Quality of Life had significantly increased over the course of 12 months from 0.80 to 0.84 (Dutch tariff). DISCUSSION: This study indicates that an optimized diagnostic trajectory has positively influenced the use of healthcare resources and the quality of life in people with epilepsy. As chronic care patients make diverse costs, future research should identify the long-term costs after an optimized diagnostic trajectory for patients with epilepsy, possibly identifying patients who are at high risk of becoming high-cost users in the future for early intervention.

Modeling of intracerebral interictal epileptic discharges: Evidence for network interactions.

OBJECTIVE: The interictal epileptic discharges (IEDs) occurring in stereotactic EEG (SEEG) recordings are in general abundant compared to ictal discharges, but difficult to interpret due to complex underlying network interactions. A framework is developed to model these network interactions. METHODS: To identify the synchronized neuronal activity underlying the IEDs, the variation in correlation over time of the SEEG signals is related to the occurrence of IEDs using the general linear model. The interdependency is assessed of the brain areas that reflect highly synchronized neural activity by applying independent component analysis, followed by cluster analysis of the spatial distributions of the independent components. The spatiotemporal interactions of the spike clusters reveal the leading or lagging of brain areas. RESULTS: The analysis framework was evaluated for five successfully operated patients, showing that the spike cluster that was related to the MRI-visible brain lesions coincided with the seizure onset zone. The additional value of the framework was demonstrated for two more patients, who were MRI-negative and for whom surgery was not successful. CONCLUSIONS: A network approach is promising in case of complex epilepsies. SIGNIFICANCE: Analysis of IEDs is considered a valuable addition to routine review of SEEG recordings, with the potential to increase the success rate of epilepsy surgery.

Stability metrics for optic radiation tractography: Towards damage prediction after resective surgery.

BACKGROUND: An accurate delineation of the optic radiation (OR) using diffusion MR tractography may reduce the risk of a visual field deficit after temporal lobe resection. However, tractography is prone to generate spurious streamlines, which deviate strongly from neighboring streamlines and hinder a reliable distance measurement between the temporal pole and the Meyer's loop (ML-TP distance). NEW METHOD: Stability metrics are introduced for the automated removal of spurious streamlines near the Meyer's loop. Firstly, fiber-to-bundle coherence (FBC) measures can identify spurious streamlines by estimating their alignment with the surrounding streamline bundle. Secondly, robust threshold selection removes spurious streamlines while preventing an underestimation of the extent of the Meyer's loop. Standardized parameter selection is realized through test-retest evaluation of the variability in ML-TP distance. RESULTS: The variability in ML-TP distance after parameter selection was below 2mm for each of the healthy volunteers studied (N=8). The importance of the stability metrics is illustrated for epilepsy surgery candidates (N=3) for whom the damage to the Meyer's loop was evaluated by comparing the pre- and post-operative OR reconstruction. The difference between predicted and observed damage is in the order of a few millimeters, which is the error in measured ML-TP distance. COMPARISON WITH EXISTING METHOD(S): The stability metrics are a novel method for the robust estimate of the ML-TP distance. CONCLUSIONS: The stability metrics are a promising tool for clinical trial studies, in which the damage to the OR can be related to the visual field deficit that may occur after epilepsy surgery.

Effectiveness of Reference Signal-Based Methods for Removal of EEG Artifacts Due to Subtle Movements During fMRI Scanning.

OBJECTIVE: Subtle motion of an epileptic patient examined with co-registered EEG and functional MRI (EEG-fMRI) may often lead to spurious fMRI activation patterns when true epileptic spikes are contaminated with motion artefacts. In recent years, methods relying on reference signals for correcting these subtle movements in the EEG have emerged. In this study, the performance of two reference-based devices are compared to the template-based method with regard to their ability to remove movement-related artifacts in EEG measured during scanning. METHODS: Measurements were performed with a novel double layer cap consisting of 29 EEG and 29 reference electrodes, and with a current loop cap consisting of 60 electrodes and three current loop wires attached to the cap. EEG was acquired inside the scanner during resting state, as well as when the subject was performing a cued movement task. For the double layer cap recordings, newly developed artifact removal algorithms are introduced and both reference signal-based methods are compared to a template-based correction method. RESULTS: The BCG artifacts occurring at resting state could be removed successfully by both the reference signal-based methods as well as by the template-based method. However, the reference signal-based methods were also capable of removing EEG artifacts induced by subtle movements, whereas the template-based method failed to remove these artifacts. CONCLUSION: Reference signal-based methods enable to correct for artifacts due to subtle movements, which are not removed by commonly used template-based removal algorithms. SIGNIFICANCE: Sensitivity of EEG-fMRI analysis in patients with focal epilepsy is improved by avoiding erroneous detections of subtle movements as epileptic spikes in the EEG.

Structural Brain Network: What is the Effect of LiFE Optimization of Whole Brain Tractography?

Structural brain networks constructed based on diffusion-weighted MRI (dMRI) have provided a systems perspective to explore the organization of the human brain. Some redundant and nonexistent fibers, however, are inevitably generated in whole brain tractography. We propose to add one critical step while constructing the networks to remove these fibers using the linear fascicle evaluation (LiFE) method, and study the differences between the networks with and without LiFE optimization. For a cohort of nine healthy adults and for 9 out of the 35 subjects from Human Connectome Project, the T 1-weighted images and dMRI data are analyzed. Each brain is parcellated into 90 regions-of-interest, whilst a probabilistic tractography algorithm is applied to generate the original connectome. The elimination of redundant and nonexistent fibers from the original connectome by LiFE creates the optimized connectome, and the random selection of the same number of fibers as the optimized connectome creates the non-optimized connectome. The combination of parcellations and these connectomes leads to the optimized and non-optimized networks, respectively. The optimized networks are constructed with six weighting schemes, and the correlations of different weighting methods are analyzed. The fiber length distributions of the non-optimized and optimized connectomes are compared. The optimized and non-optimized networks are compared with regard to edges, nodes and networks, within a sparsity range of 0.75-0.95. It has been found that relatively more short fibers exist in the optimized connectome. About 24.0% edges of the optimized network are significantly different from those in the non-optimized network at a sparsity of 0.75. About 13.2% of edges are classified as false positives or the possible missing edges. The strength and betweenness centrality of some nodes are significantly different for the non-optimized and optimized networks, but not the node efficiency. The normalized clustering coefficient, the normalized characteristic path length and the small-worldness are higher in the optimized network weighted by the fiber number than in the non-optimized network. These observed differences suggest that LiFE optimization can be a crucial step for the construction of more reasonable and more accurate structural brain networks.

Spatiotemporal mapping of interictal epileptiform discharges in human absence epilepsy: A MEG study.

PURPOSE: Although absence epilepsy is considered to be a prototypic type of generalized epilepsy, it is still under debate whether generalized 3 Hz spike-and-wave discharges (SWDs) might have a cortical focal origin. Here it is investigated whether focal interictal epileptiform discharges (IEDs), which typically occur in the electro- (EEG) and magnetoencephalogram (MEG) in case of focal epilepsy, are present in the MEG of children with absence epilepsy. Next, the location of the sources of the IEDs is established, and it is investigated whether the location is concordant to the earlier established focal cortical regions involved in the generalized SWDs of these children. METHODS: Whole head MEG recordings of seven children with absence epilepsy were reviewed with respect to the presence of IEDs (spikes and sharp waves). These IEDs were grouped into distinct clusters, in which each contribution to a cluster yields a comparable magnetic field distribution. Source localization was then performed onto the average signal of each cluster using an equivalent current dipole model and a realistic head model of the cortical surface. RESULTS: IEDs were detected in 6 out of 7 patients. Source reconstruction indicated most often frontal, central or parietal origins of the IED in either the left and or right hemisphere. Spatiotemporal assessment of the IEDs indicated a stable location of the averages of these discharges, indicating a single underlying cortical source. DISCUSSION: The outcome of this pilot study shows that MEG is well suited for the detection of IEDs and suggests that their estimated sources coincide rather well with the cortical regions involved during the spikes of the SWDs. It is discussed whether the presence of IEDs, classically seen as a marker of focal epilepsies, indicate that absence epilepsy should be considered as a focal type of epilepsy, in which changes in the network are evolving rapidly.

The influence of construction methodology on structural brain network measures: A review.

Structural brain networks based on diffusion MRI and tractography show robust attributes such as small-worldness, hierarchical modularity, and rich-club organization. However, there are large discrepancies in the reports about specific network measures. It is hypothesized that these discrepancies result from the influence of construction methodology. We surveyed the methodological options and their influences on network measures. It is found that most network measures are sensitive to the scale of brain parcellation, MRI gradient schemes and orientation model, and the tractography algorithm, which is in accordance with the theoretical analysis of the small-world network model. Different network weighting schemes represent different attributes of brain networks, which makes these schemes incomparable between studies. Methodology choice depends on the specific study objectives and a clear understanding of the pros and cons of a particular methodology. Because there is no way to eliminate these influences, it seems more practical to quantify them, optimize the methodologies, and construct structural brain networks with multiple spatial resolutions, multiple edge densities, and multiple weighting schemes.

Network analysis of EEG related functional MRI changes due to medication withdrawal in focal epilepsy.

Anti-epileptic drugs (AEDs) have a global effect on the neurophysiology of the brain which is most likely reflected in functional brain activity recorded with EEG and fMRI. These effects may cause substantial inter-subject variability in studies where EEG correlated functional MRI (EEG-fMRI) is used to determine the epileptogenic zone in patients who are candidate for epilepsy surgery. In the present study the effects on resting state fMRI are quantified in conditions with AED administration and after withdrawal of AEDs. EEG-fMRI data were obtained from 10 patients in the condition that the patient was on the steady-state maintenance doses of AEDs as prescribed (condition A) and after withdrawal of AEDs (condition B), at the end of a clinically standard pre-surgical long term video-EEG monitoring session. Resting state networks (RSN) were extracted from fMRI. The epileptic component (ICE) was identified by selecting the RSN component with the largest overlap with the EEG-fMRI correlation pattern. Changes in RSN functional connectivity between conditions A and B were quantified. EEG-fMRI correlation analysis was successful in 30% and 100% of the cases in conditions A and B, respectively. Spatial patterns of ICEs are comparable in conditions A and B, except for one patient for whom it was not possible to identify the ICE in condition A. However, the resting state functional connectivity is significantly increased in the condition after withdrawal of AEDs (condition B), which makes resting state fMRI potentially a new tool to study AED effects. The difference in sensitivity of EEG-fMRI in conditions A and B, which is not related to the number of epileptic EEG events occurring during scanning, could be related to the increased functional connectivity in condition B.

Are Epilepsy-Related fMRI Components Dependent on the Presence of Interictal Epileptic Discharges in Scalp EEG?

Spatial independent component analysis (ICA) is increasingly being used to extract resting-state networks from fMRI data. Previous studies showed that ICA also reveals independent components (ICs) related to the seizure onset zone. However, it is currently unknown how these epileptic ICs depend on the presence of interictal epileptic discharges (IEDs) in the EEG. The goal of this study was to explore the relation between ICs obtained from fMRI epochs during the occurrence of IEDs in the EEG and those without IEDs. fMRI data sets with co-registered EEG were retrospectively selected of patients from whom the location of the epileptogenic zone was confirmed by outcome of surgery (n = 8). The fMRI data were split into two epochs: one with IEDs visible in scalp EEG and one without. Spatial ICA was applied to the fMRI data of each part separately. The maps of all resulting components were compared to the resection area and the EEG-fMRI correlation pattern by computing a spatial correlation coefficient to detect the epilepsy-related component. For all patients, except one, there was a remarkable resemblance between the epilepsy-related components selected during epochs with IEDs and those without IEDs. These findings suggest that epilepsy-related ICs are not dependent on the presence of IEDs in scalp EEG. Since these epileptic ICs showed partial overlap with resting-state networks of healthy volunteers (n = 10), our study supports the need for new ways to classify epileptic ICs.

Evaluating contextual processing in diffusion MRI: application to optic radiation reconstruction for epilepsy surgery.

Diffusion MRI and tractography allow for investigation of the architectural configuration of white matter in vivo, offering new avenues for applications like presurgical planning. Despite the promising outlook, there are many pitfalls that complicate its use for (clinical) application. Amongst these are inaccuracies in the geometry of the diffusion profiles on which tractography is based, and poor alignment with neighboring profiles. Recently developed contextual processing techniques, including enhancement and well-posed geometric sharpening, have shown to result in sharper and better aligned diffusion profiles. However, the research that has been conducted up to now is mainly of theoretical nature, and so far these techniques have only been evaluated by visual inspection of the diffusion profiles. In this work, the method is evaluated in a clinically relevant application: the reconstruction of the optic radiation for epilepsy surgery. For this evaluation we have developed a framework in which we incorporate a novel scoring procedure for individual pathways. We demonstrate that, using enhancement and sharpening, the extraction of an anatomically plausible reconstruction of the optic radiation from a large amount of probabilistic pathways is greatly improved in three healthy controls, where currently used methods fail to do so. Furthermore, challenging reconstructions of the optic radiation in three epilepsy surgery candidates with extensive brain lesions demonstrate that it is beneficial to integrate these methods in surgical planning.

High density wireless EEG prototype: Design and evaluation against reference equipment.

A high density wireless electroencephalographic (EEG) platform has been designed. It is able to record up to 64 EEG channels with electrode to tissue impedance (ETI) monitoring. The analog front-end is based on two kinds of low power ASICs implementing the active electrodes and the amplifier. A power efficient compression algorithm enables the use of continuous wireless transmission of data through Bluetooth for real-time monitoring with an overall power consumption of about 350 mW. EEG acquisitions on five subjects (one healthy subject and four patients suffering from epilepsy) have been recorded in parallel with a reference system commonly used in clinical practice and data of the wireless prototype and reference system have been processed with an automatic tool for seizure detection and localization. The false alarm rates (0.1-0.5 events per hour) are comparable between the two system and wireless prototype also detected the seizure correctly and allowed its localization.

EEG-fMRI correlation patterns in the presurgical evaluation of focal epilepsy: a comparison with electrocorticographic data and surgical outcome measures.

EEG-correlated functional MRI (EEG-fMRI) visualizes brain regions associated with interictal epileptiform discharges (IEDs). This technique images the epileptiform network, including multifocal, superficial and deeply situated cortical areas. To understand the role of EEG-fMRI in presurgical evaluation, its results should be validated relative to a gold standard. For that purpose, EEG-fMRI data were acquired for a heterogeneous group of surgical candidates (n=16) who were later implanted with subdural grids and strips (ECoG). The EEG-fMRI correlation patterns were systematically compared with brain areas involved in IEDs ECoG, using a semi-automatic analysis method, as well as to the seizure onset zone, resected area, and degree of seizure freedom. In each patient at least one of the EEG-fMRI areas was concordant with an interictally active ECoG area, always including the early onset area of IEDs in the ECoG data. This confirms that EEG-fMRI reflects a pattern of onset and propagation of epileptic activity. At group level, 76% of the BOLD regions that were covered with subdural grids, were concordant with interictally active ECoG electrodes. Due to limited spatial sampling, 51% of the BOLD regions were not covered with electrodes and could, therefore, not be validated. From an ECoG perspective it appeared that 29% of the interictally active ECoG regions were missed by EEG-fMRI and that 68% of the brain regions were correctly identified as inactive with EEG-fMRI. Furthermore, EEG-fMRI areas included the complete seizure onset zone in 83% and resected area in 93% of the data sets. No clear distinction was found between patients with a good or poor surgical outcome: in both patient groups, EEG-fMRI correlation patterns were found that were either focal or widespread. In conclusion, by comparison of EEG-fMRI with interictal invasive EEG over a relatively large patient population we were able to show that the EEG-fMRI correlation patterns are spatially accurate at the level of neurosurgical units (i.e. anatomical brain regions) and reflect the underlying network of IEDs. Therefore, we expect that EEG-fMRI can play an important role for the determination of the implantation strategy.

Novel artefact removal algorithms for co-registered EEG/fMRI based on selective averaging and subtraction.

Co-registered EEG and functional MRI (EEG/fMRI) is a potential clinical tool for planning invasive EEG in patients with epilepsy. In addition, the analysis of EEG/fMRI data provides a fundamental insight into the precise physiological meaning of both fMRI and EEG data. Routine application of EEG/fMRI for localization of epileptic sources is hampered by large artefacts in the EEG, caused by switching of scanner gradients and heartbeat effects. Residuals of the ballistocardiogram (BCG) artefacts are similarly shaped as epileptic spikes, and may therefore cause false identification of spikes. In this study, new ideas and methods are presented to remove gradient artefacts and to reduce BCG artefacts of different shapes that mutually overlap in time. Gradient artefacts can be removed efficiently by subtracting an average artefact template when the EEG sampling frequency and EEG low-pass filtering are sufficient in relation to MR gradient switching (Gonçalves et al., 2007). When this is not the case, the gradient artefacts repeat themselves at time intervals that depend on the remainder between the fMRI repetition time and the closest multiple of the EEG acquisition time. These repetitions are deterministic, but difficult to predict due to the limited precision by which these timings are known. Therefore, we propose to estimate gradient artefact repetitions using a clustering algorithm, combined with selective averaging. Clustering of the gradient artefacts yields cleaner EEG for data recorded during scanning of a 3T scanner when using a sampling frequency of 2048 Hz. It even gives clean EEG when the EEG is sampled with only 256 Hz. Current BCG artefacts-reduction algorithms based on average template subtraction have the intrinsic limitation that they fail to deal properly with artefacts that overlap in time. To eliminate this constraint, the precise timings of artefact overlaps were modelled and represented in a sparse matrix. Next, the artefacts were disentangled with a least squares procedure. The relevance of this approach is illustrated by determining the BCG artefacts in a data set consisting of 29 healthy subjects recorded in a 1.5 T scanner and 15 patients with epilepsy recorded in a 3 T scanner. Analysis of the relationship between artefact amplitude, duration and heartbeat interval shows that in 22% (1.5T data) to 30% (3T data) of the cases BCG artefacts show an overlap. The BCG artefacts of the EEG/fMRI data recorded on the 1.5T scanner show a small negative correlation between HBI and BCG amplitude. In conclusion, the proposed methodology provides a substantial improvement of the quality of the EEG signal without excessive computer power or additional hardware than standard EEG-compatible equipment.

A semi-automatic method to determine electrode positions and labels from gel artifacts in EEG/fMRI-studies.

The analysis of simultaneous EEG and fMRI data is generally based on the extraction of regressors of interest from the EEG, which are correlated to the fMRI data in a general linear model setting. In more advanced approaches, the spatial information of EEG is also exploited by assuming underlying dipole models. In this study, we present a semi automatic and efficient method to determine electrode positions from electrode gel artifacts, facilitating the integration of EEG and fMRI in future EEG/fMRI data models. In order to visualize all electrode artifacts simultaneously in a single view, a surface rendering of the structural MRI is made using a skin triangular mesh model as reference surface, which is expanded to a "pancake view". Then the electrodes are determined with a simple mouse click for each electrode. Using the geometry of the skin surface and its transformation to the pancake view, the 3D coordinates of the electrodes are reconstructed in the MRI coordinate frame. The electrode labels are attached to the electrode positions by fitting a template grid of the electrode cap in which the labels are known. The correspondence problem between template and sample electrodes is solved by minimizing a cost function over rotations, shifts and scalings of the template grid. The crucial step here is to use the solution of the so-called "Hungarian algorithm" as a cost function, which makes it possible to identify the electrode artifacts in arbitrary order. The template electrode grid has to be constructed only once for each cap configuration. In our implementation of this method, the whole procedure can be performed within 15 min including import of MRI, surface reconstruction and transformation, electrode identification and fitting to template. The method is robust in the sense that an electrode template created for one subject can be used without identification errors for another subject for whom the same EEG cap was used. Furthermore, the method appears to be robust against spurious or missing artifacts. We therefore consider the proposed method as a useful and reliable tool within the larger toolbox required for the analysis of co-registered EEG/fMRI data.

Space-time network connectivity and cortical activations preceding spike wave discharges in human absence epilepsy: a MEG study.

To describe the spatial and temporal profiles of connectivity networks and sources preceding generalized spike-and-wave discharges (SWDs) in human absence epilepsy. Nonlinear associations of MEG signals and cluster indices obtained within the framework of graph theory were determined, while source localization in the frequency domain was performed in the low frequency bands with dynamic imaging of coherent sources. The results were projected on a three-dimensional surface rendering of the brain using a semi-realistic head model and MRI images obtained for each of the five patients studied. An increase in clustering and a decrease in path length preceding SWD onset and a rhythmic pattern of increasing and decreasing connectivity were seen during SWDs. Beamforming showed a consistent appearance of a low frequency frontal cortical source prior to the first generalized spikes. This source was preceded by a low frequency occipital source. The changes in the connectivity networks with the onset of SWDs suggest a pathologically predisposed state towards synchronous seizure networks with increasing connectivity from interictal to preictal and ictal state, while the occipital and frontal low frequency early preictal sources demonstrate that SWDs are not suddenly arising but gradually build up in a dynamic network.

Comparison of analytical strategies for EEG-correlated fMRI data in patients with epilepsy.

The simultaneous recording of electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) can be used to localize interictal epileptiform discharges (IEDs). Previous studies have reported varying degrees of concordance of EEG-fMRI with electroclinical findings. The aim of the present study is to evaluate to what extent this variability is determined by the analytical strategy or by the properties of the EEG data. For that purpose, 42 IED sets obtained in 29 patients with epilepsy were reanalyzed using a finite impulse response approach, which estimates the hemodynamic response function (HRF) from the data and allows non-causal effects. Cardiac effects were treated as additional confounders in the model. This approach was compared to the classical approach assuming a fixed HRF for each voxel in the brain. The performance of each method was assessed by comparing the fMRI results to the EEG focus. The flexible model revealed more significantly activated voxels, which resulted in more activated brain regions concordant with the EEG focus (26 vs. 16). Correction for cardiac effects improved the results in 7 out of the 42 data sets. Furthermore, design theory for event-related experiments was applied in order to determine the influence of the number of IEDs and their temporal distribution on the success of an experiment. It appeared that this success is highly dependent upon the number of IEDs present during the recording and less on their temporal spacing. We conclude that the outcome of EEG-fMRI can be improved by using an optimized analytical strategy, but also depends on the number of IEDs occurring during the recording.

Correction for pulse height variability reduces physiological noise in functional MRI when studying spontaneous brain activity.

EEG correlated functional MRI (EEG-fMRI) allows the delineation of the areas corresponding to spontaneous brain activity, such as epileptiform spikes or alpha rhythm. A major problem of fMRI analysis in general is that spurious correlations may occur because fMRI signals are not only correlated with the phenomena of interest, but also with physiological processes, like cardiac and respiratory functions. The aim of this study was to reduce the number of falsely detected activated areas by taking the variation in physiological functioning into account in the general linear model (GLM). We used the photoplethysmogram (PPG), since this signal is based on a linear combination of oxy- and deoxyhemoglobin in the arterial blood, which is also the basis of fMRI. We derived a regressor from the variation in pulse height (VIPH) of PPG and added this regressor to the GLM. When this regressor was used as predictor it appeared that VIPH explained a large part of the variance of fMRI signals acquired from five epilepsy patients and thirteen healthy volunteers. As a confounder VIPH reduced the number of activated voxels by 30% for the healthy volunteers, when studying the generators of the alpha rhythm. Although for the patients the number of activated voxels either decreased or increased, the identification of the epileptogenic zone was substantially enhanced in one out of five patients, whereas for the other patients the effects were smaller. In conclusion, applying VIPH as a confounder diminishes physiological noise and allows a more reliable interpretation of fMRI results.