Head model and electrical source imaging: a study of 38 epileptic patients

Dysfunctional gaze processing in bipolar disorder

Gaze conveys emotional information, and humans present sensitivity to its direction from the earliest days of life. Bipolar disorder is a disease characterized by fluctuating states of emotional and cognitive dysregulation. To explore the role of attentional control on face processing in bipolar patients (BP) we used gaze direction as an emotion modulation parameter in a two-back Working Memory (WM) task while high-density EEG data were acquired. Since gaze direction influences emotional attributions to faces with neutral expressions as well, we presented neutral faces with direct and averted gaze. Nineteen euthymic BP and a sample of age- and gender-matched controls were examined.

In BP we observed diminished P200 and augmented P300 evoked responses, differentially modulated by non-repeated or repeated faces, as well as by gaze direction. BP showed a reduced P200 amplitude, significantly stronger for faces with direct gaze than averted gaze. Source localization of P200 indicated decreased activity in sensory-motor regions and frontal areas suggestive of abnormal affective processing of neutral faces.

The present study provides neurophysiological evidence for abnormal gaze processing in BP and suggests dysfunctional processing of direct eye contact as a prominent characteristic of bipolar disorder.

•

This ERP study identified abnormalities in gaze processing in bipolar patients.

•

We observed functional anomalies in the P200 and P300 evoked responses.

•

BP showed a strong suppression of the P200 for faces with direct gaze.

•

Source localization indicated decreased activity in sensory-motor regions.

Localizing value of electrical source imaging: Frontal lobe, malformations of cortical development and negative MRI related epilepsies are the best candidates

OBJECTIVE

We aimed to prospectively assess the anatomical concordance of electric source localizations of interictal discharges with the epileptogenic zone (EZ) estimated by stereo-electroencephalography (SEEG) according to different subgroups: the type of epilepsy, the presence of a structural MRI lesion, the aetiology and the depth of the EZ.

METHODS

In a prospective multicentric observational study, we enrolled 85 consecutive patients undergoing pre-surgical SEEG investigation for focal drug-resistant epilepsy. Electric source imaging (ESI) was performed before SEEG. Source localizations were obtained from dipolar and distributed source methods. Anatomical concordance between ESI and EZ was defined according to 36 predefined sublobar regions. ESI was interpreted blinded to- and subsequently compared with SEEG estimated EZ.

RESULTS

74 patients were finally analyzed. 38 patients had temporal and 36 extra-temporal lobe epilepsy. MRI was positive in 52. 41 patients had malformation of cortical development (MCD), 33 had another or an unknown aetiology. EZ was medial in 27, lateral in 13, and medio-lateral in 34. In the overall cohort, ESI completely or partly localized the EZ in 85%: full concordance in 13 cases and partial concordance in 50 cases. The rate of ESI full concordance with EZ was significantly higher in (i) frontal lobe epilepsy (46%; p = 0.05), (ii) cases of negative MRI (36%; p = 0.01) and (iii) MCD (27%; p = 0.03). The rate of ESI full concordance with EZ was not statistically different according to the depth of the EZ.

SIGNIFICANCE

We prospectively demonstrated that ESI more accurately estimated the EZ in subgroups of patients who are often the most difficult cases in epilepsy surgery: frontal lobe epilepsy, negative MRI and the presence of MCD.

Automated long-term EEG analysis to localize the epileptogenic zone

Objective

We investigated the performance of automatic spike detection and subsequent electroencephalogram (EEG) source imaging to localize the epileptogenic zone (EZ) from long-term EEG recorded during video-EEG monitoring.

Methods

In 32 patients, spikes were automatically detected in the EEG and clustered according to their morphology. The two spike clusters with most single events in each patient were averaged and localized in the brain at the half-rising time and peak of the spike using EEG source imaging. On the basis of the distance from the sources to the resection and the known patient outcome after surgery, the performance of the automated EEG analysis to localize the EZ was quantified.

Results

In 28 out of the 32 patients, the automatically detected spike clusters corresponded with the reported interictal findings. The median distance to the resection in patients with Engel class I outcome was 6.5 and 15 mm for spike cluster 1 and 27 and 26 mm for cluster 2, at the peak and the half-rising time of the spike, respectively. Spike occurrence (cluster 1 vs. cluster 2) and spike timing (peak vs. half-rising) significantly influenced the distance to the resection (p < 0.05). For patients with Engel class II, III, and IV outcomes, the median distance increased to 36 and 36 mm for cluster 1. Localizing spike cluster 1 at the peak resulted in a sensitivity of 70% and specificity of 100%, positive prediction value (PPV) of 100%, and negative predictive value (NPV) of 53%. Including the results of spike cluster 2 led to an increased sensitivity of 79% NPV of 55% and diagnostic OR of 11.4, while the specificity dropped to 75% and the PPV to 90%.

Significance

We showed that automated analysis of long-term EEG recordings results in a high sensitivity and specificity to localize the epileptogenic focus.

SEEG dipole source localization based on an empirical Bayesian approach taking into account forward model uncertainties

Electromagnetic brain source localization consists in the inversion of a forward model based on a limited number of potential measurements. A wide range of methods has been developed to regularize this severely ill-posed problem and to reduce the solution space, imposing spatial smoothness, anatomical constraint or sparsity of the activated source map. This last criteria, based on physiological assumptions stating that in some particular events (e.g., epileptic spikes, evoked potential) few focal area of the brain are simultaneously actives, has gained more and more interest. Bayesian approaches have the ability to provide sparse solutions under adequate parametrization, and bring a convenient framework for the introduction of priors in the form of probabilistic density functions. However the quality of the forward model is rarely questioned while this parameter has undoubtedly a great influence on the solution. Its construction suffers from numerous approximation and uncertainties, even when using realistic numerical models. In addition, it often encodes a coarse sampling of the continuous solution space due to the computational burden its inversion implies. In this work we propose an empirical Bayesian approach to take into account the uncertainties of the forward model by allowing constrained variations around a prior physical model, in the particular context of SEEG measurements. We demonstrate on simulations that the method enhance the accuracy of the source time-course estimation as well as the sparsity of the resulting source map. Results on real signals prove the applicability of the method in real contexts.

Directed Functional Brain Connectivity Based on EEG Source Imaging: Methodology and Application to Temporal Lobe Epilepsy

OBJECTIVE

The importance of functional brain connectivity to study physiological and pathological brain activity has been widely recognized. Here, we aimed to 1) review a methodological pipeline to investigate directed functional connectivity between brain regions using source signals derived from high-density EEG; 2) elaborate on some methodological challenges; and 3) apply this pipeline to temporal lobe epilepsy (TLE) patients and healthy controls to investigate directed functional connectivity differences in the theta and beta frequency bands during EEG epochs without visible pathological activity.

METHODS

The methodological pipeline includes: EEG acquisition and preprocessing, electrical-source imaging (ESI) using individual head models and distributed inverse solutions, parcellation of the gray matter in regions of interest, fixation of the dipole orientation for each region, computation of the spectral power in the source space, and directed functional connectivity estimation using Granger-causal modeling. We specifically analyzed how the signal-to-noise ratio (SNR) changes using different approaches for the dipole orientation fixation. We applied this pipeline to 20 left TLE patients, 20 right TLE patients, and 20 healthy controls.

RESULTS

Projecting each dipole to the predominant dipole orientation leads to a threefold SNR increase as compared to the norm of the dipoles. By comparing connectivity in TLE versus controls, we found significant frequency-specific outflow differences in physiologically plausible regions.

CONCLUSION

The results suggest that directed functional connectivity derived from ESI can help better understand frequency-specific resting-state network alterations underlying focal epilepsy.

SIGNIFICANCE

EEG-based directed functional connectivity could contribute to the search of new biomarkers of this disorder.

In-vivo measurements of human brain tissue conductivity using focal electrical current injection through intracerebral multicontact electrodes

In-vivo measurements of human brain tissue conductivity at body temperature were conducted using focal electrical currents injected through intracerebral multicontact electrodes. A total of 1,421 measurements in 15 epileptic patients (age: 28 ± 10) using a radiofrequency generator (50 kHz current injection) were analyzed. Each contact pair was classified as being from healthy (gray matter, n = 696; white matter, n = 530) or pathological (epileptogenic zone, n = 195) tissue using neuroimaging analysis of the local tissue environment and intracerebral EEG recordings. Brain tissue conductivities were obtained using numerical simulations based on conductivity estimates that accounted for the current flow in the local brain volume around the contact pairs (a cube with a side length of 13 mm). Conductivity values were 0.26 S/m for gray matter and 0.17 S/m for white matter. Healthy gray and white matter had statistically different median impedances (P < 0.0001). White matter conductivity was found to be homogeneous as normality tests did not find evidence of multiple subgroups. Gray matter had lower conductivity in healthy tissue than in the epileptogenic zone (0.26 vs. 0.29 S/m; P = 0.012), even when the epileptogenic zone was not visible in the magnetic resonance image (MRI) (P = 0.005). The present in-vivo conductivity values could serve to create more accurate volume conduction models and could help to refine the identification of relevant intracerebral contacts, especially when located within the epileptogenic zone of an MRI-invisible lesion. Hum Brain Mapp 38:974-986, 2017. © 2016 Wiley Periodicals, Inc.

The Role of Skull Modeling in EEG Source Imaging for Patients with Refractory Temporal Lobe Epilepsy

We investigated the influence of different skull modeling approaches on EEG source imaging (ESI), using data of six patients with refractory temporal lobe epilepsy who later underwent successful epilepsy surgery. Four realistic head models with different skull compartments, based on finite difference methods, were constructed for each patient: (i) Three models had skulls with compact and spongy bone compartments as well as air-filled cavities, segmented from either computed tomography (CT), magnetic resonance imaging (MRI) or a CT-template and (ii) one model included a MRI-based skull with a single compact bone compartment. In all patients we performed ESI of single and averaged spikes marked in the clinical 27-channel EEG by the epileptologist. To analyze at which time point the dipole estimations were closer to the resected zone, ESI was performed at two time instants: the half-rising phase and peak of the spike. The estimated sources for each model were validated against the resected area, as indicated by the postoperative MRI. Our results showed that single spike analysis was highly influenced by the signal-to-noise ratio (SNR), yielding estimations with smaller distances to the resected volume at the peak of the spike. Although averaging reduced the SNR effects, it did not always result in dipole estimations lying closer to the resection. The proposed skull modeling approaches did not lead to significant differences in the localization of the irritative zone from clinical EEG data with low spatial sampling density. Furthermore, we showed that a simple skull model (MRI-based) resulted in similar accuracy in dipole estimation compared to more complex head models (based on CT- or CT-template). Therefore, all the considered head models can be used in the presurgical evaluation of patients with temporal lobe epilepsy to localize the irritative zone from low-density clinical EEG recordings.

Brain imaging in the assessment for epilepsy surgery

Brain imaging has a crucial role in the presurgical assessment of patients with epilepsy. Structural imaging reveals most cerebral lesions underlying focal epilepsy. Advances in MRI acquisitions including diffusion-weighted imaging, post-acquisition image processing techniques, and quantification of imaging data are increasing the accuracy of lesion detection. Functional MRI can be used to identify areas of the cortex that are essential for language, motor function, and memory, and tractography can reveal white matter tracts that are vital for these functions, thus reducing the risk of epilepsy surgery causing new morbidities. PET, SPECT, simultaneous EEG and functional MRI, and electrical and magnetic source imaging can be used to infer the localisation of epileptic foci and assist in the design of intracranial EEG recording strategies. Progress in semi-automated methods to register imaging data into a common space is enabling the creation of multimodal three-dimensional patient-specific datasets. These techniques show promise for the demonstration of the complex relations between normal and abnormal structural and functional data and could be used to direct precise intracranial navigation and surgery for individual patients.

Electrical source imaging of interictal spikes using multiple sparse volumetric priors for presurgical epileptogenic focus localization

Electrical source imaging of interictal spikes observed in EEG recordings of patients with refractory epilepsy provides useful information to localize the epileptogenic focus during the presurgical evaluation. However, the selection of the time points or time epochs of the spikes in order to estimate the origin of the activity remains a challenge. In this study, we consider a Bayesian EEG source imaging technique for distributed sources, i.e. the multiple volumetric sparse priors (MSVP) approach. The approach allows to estimate the time courses of the intensity of the sources corresponding with a specific time epoch of the spike. Based on presurgical averaged interictal spikes in six patients who were successfully treated with surgery, we estimated the time courses of the source intensities for three different time epochs: (i) an epoch starting 50 ms before the spike peak and ending at 50% of the spike peak during the rising phase of the spike, (ii) an epoch starting 50 ms before the spike peak and ending at the spike peak and (iii) an epoch containing the full spike time period starting 50 ms before the spike peak and ending 230 ms after the spike peak. To identify the primary source of the spike activity, the source with the maximum energy from 50 ms before the spike peak till 50% of the spike peak was subsequently selected for each of the time windows. For comparison, the activity at the spike peaks and at 50% of the peaks was localized using the LORETA inversion technique and an ECD approach. Both patient-specific spherical forward models and patient-specific 5-layered finite difference models were considered to evaluate the influence of the forward model. Based on the resected zones in each of the patients, extracted from post-operative MR images, we compared the distances to the resection border of the estimated activity. Using the spherical models, the distances to the resection border for the MSVP approach and each of the different time epochs were in the same range as the LORETA and ECD techniques. We found distances smaller than 23 mm, with robust results for all the patients. For the finite difference models, we found that the distances to the resection border for the MSVP inversions of the full spike time epochs were generally smaller compared to the MSVP inversions of the time epochs before the spike peak. The results also suggest that the inversions using the finite difference models resulted in slightly smaller distances to the resection border compared to the spherical models. The results we obtained are promising because the MSVP approach allows to study the network of the estimated source-intensities and allows to characterize the spatial extent of the underlying sources.

•

A Bayesian ESI technique is evaluated to localize interictal spike activity.

•

Averaged spikes in six patients were used that were seizure free after surgery.

•

We compared the technique with the LORETA an ECD technique.

•

We evaluated both spherical and 5-layered finite difference forward models.

•

Our approach is potentially useful to delineate the irritative zone.

Early recurrence and ongoing parietal driving during elementary visual processing

Visual stimuli quickly activate a broad network of brain areas that often show reciprocal structural connections between them. Activity at short latencies (<100 ms) is thought to represent a feed-forward activation of widespread cortical areas, but fast activation combined with reciprocal connectivity between areas in principle allows for two-way, recurrent interactions to occur at short latencies after stimulus onset. Here we combined EEG source-imaging and Granger-causal modeling with high temporal resolution to investigate whether recurrent and top-down interactions between visual and attentional brain areas can be identified and distinguished at short latencies in humans. We investigated the directed interactions between widespread occipital, parietal and frontal areas that we localized within participants using fMRI. The connectivity results showed two-way interactions between area MT and V1 already at short latencies. In addition, the results suggested a large role for lateral parietal cortex in coordinating visual activity that may be understood as an ongoing top-down allocation of attentional resources. Our results support the notion that indirect pathways allow early, evoked driving from MT to V1 to highlight spatial locations of motion transients, while influence from parietal areas is continuously exerted around stimulus onset, presumably reflecting task-related attentional processes.

Intrinsic coupling modes in source-reconstructed electroencephalography

Intrinsic coupling of neuronal assemblies constitutes a key feature of ongoing brain activity, yielding the rich spatiotemporal patterns observed in neuroimaging data and putatively supporting cognitive processes. Intrinsic coupling has been investigated in electrophysiological recordings using two types of functional connectivity measures: amplitude and phase coupling. These two coupling modes differ in their likely causes and functions, and have been proposed to provide complementary insights into intrinsic neuronal interactions. Here, we investigate the relationship between amplitude and phase coupling in source-reconstructed electroencephalography (EEG). Volume conduction is a key obstacle for connectivity analysis in EEG-we therefore also test the envelope correlation of orthogonalized signals and the phase lag index. Functional connectivity between six seed source regions (bilateral visual, sensorimotor, and auditory cortices) and all other cortical voxels was computed. For all four measures, coupling between homologous sensory areas in both hemispheres was significantly higher than with other voxels at the same physical distance. The frequency of significant coupling differed between sensory areas: 10 Hz for visual, 30 Hz for auditory, and 40 Hz for sensorimotor cortices. By contrasting envelope correlations and phase locking values, we observed two distinct clusters of voxels showing a different relationship between amplitude and phase coupling. Large clusters contiguous to the seed regions showed an identity (1:1) relationship between amplitude and phase coupling, whereas a cluster located around the contralateral homologous regions showed higher phase than amplitude coupling. These results show a relationship between intrinsic coupling modes that is distinct from the effect of volume conduction.

Brain Source Imaging in Preclinical Rat Models of Focal Epilepsy using High-Resolution EEG Recordings

Electroencephalogram (EEG) has been traditionally used to determine which brain regions are the most likely candidates for resection in patients with focal epilepsy. This methodology relies on the assumption that seizures originate from the same regions of the brain from which interictal epileptiform discharges (IEDs) emerge. Preclinical models are very useful to find correlates between IED locations and the actual regions underlying seizure initiation in focal epilepsy. Rats have been commonly used in preclinical studies of epilepsy; hence, there exist a large variety of models for focal epilepsy in this particular species. However, it is challenging to record multichannel EEG and to perform brain source imaging in such a small animal. To overcome this issue, we combine a patented-technology to obtain 32-channel EEG recordings from rodents and an MRI probabilistic atlas for brain anatomical structures in Wistar rats to perform brain source imaging. In this video, we introduce the procedures to acquire multichannel EEG from Wistar rats with focal cortical dysplasia, and describe the steps both to define the volume conductor model from the MRI atlas and to uniquely determine the IEDs. Finally, we validate the whole methodology by obtaining brain source images of IEDs and compare them with those obtained at different time frames during the seizure onset.

Dynamic directed interictal connectivity in left and right temporal lobe epilepsy

OBJECTIVE

There is increasing evidence that epileptic activity involves widespread brain networks rather than single sources and that these networks contribute to interictal brain dysfunction. We investigated the fast-varying behavior of epileptic networks during interictal spikes in right and left temporal lobe epilepsy (RTLE and LTLE) at a whole-brain scale using directed connectivity.

METHODS

In 16 patients, 8 with LTLE and 8 with RTLE, we estimated the electrical source activity in 82 cortical regions of interest (ROIs) using high-density electroencephalography (EEG), individual head models, and a distributed linear inverse solution. A multivariate, time-varying, and frequency-resolved Granger-causal modeling (weighted Partial Directed Coherence) was applied to the source signal of all ROIs. A nonparametric statistical test assessed differences between spike and baseline epochs. Connectivity results between RTLE and LTLE were compared between RTLE and LTLE and with neuropsychological impairments.

RESULTS

Ipsilateral anterior temporal structures were identified as key drivers for both groups, concordant with the epileptogenic zone estimated invasively. We observed an increase in outflow from the key driver already before the spike. There were also important temporal and extratemporal ipsilateral drivers in both conditions, and contralateral only in RTLE. A different network pattern between LTLE and RTLE was found: in RTLE there was a much more prominent ipsilateral to contralateral pattern than in LTLE. Half of the RTLE patients but none of the LTLE patients had neuropsychological deficits consistent with contralateral temporal lobe dysfunction, suggesting a relationship between connectivity changes and cognitive deficits.

SIGNIFICANCE

The different patterns of time-varying connectivity in LTLE and RTLE suggest that they are not symmetrical entities, in line with our neuropsychological results. The highest outflow region was concordant with invasive validation of the epileptogenic zone. This enhanced characterization of dynamic connectivity patterns could better explain cognitive deficits and help the management of epilepsy surgery candidates.

Mapping epileptic activity: sources or networks for the clinicians?

Epileptic seizures of focal origin are classically considered to arise from a focal epileptogenic zone and then spread to other brain regions. This is a key concept for semiological electro-clinical correlations, localization of relevant structural lesions, and selection of patients for epilepsy surgery. Recent development in neuro-imaging and electro-physiology and combinations, thereof, have been validated as contributory tools for focus localization. In parallel, these techniques have revealed that widespread networks of brain regions, rather than a single epileptogenic region, are implicated in focal epileptic activity. Sophisticated multimodal imaging and analysis strategies of brain connectivity patterns have been developed to characterize the spatio-temporal relationships within these networks by combining the strength of both techniques to optimize spatial and temporal resolution with whole-brain coverage and directional connectivity. In this paper, we review the potential clinical contribution of these functional mapping techniques as well as invasive electrophysiology in human beings and animal models for characterizing network connectivity.