Complementary structural and functional abnormalities to localise epileptogenic tissue

Abstract

Background: When investigating suitability for epilepsy surgery, people with drug-refractory focal epilepsy may have intracranial EEG (iEEG) electrodes implanted to localise seizure onset. Diffusion-weighted magnetic resonance imaging (dMRI) may be acquired to identify key white matter tracts for surgical avoidance. Here, we investigate whether structural connectivity abnormalities, inferred from dMRI, may be used in conjunction with functional iEEG abnormalities to aid localisation of the epileptogenic zone (EZ), improving surgical outcomes in epilepsy.

Methods: We retrospectively investigated data from 43 patients with epilepsy who had surgery following iEEG. Twenty-five patients (58%) were free from disabling seizures (ILAE 1 or 2) at one year. Interictal iEEG functional, and dMRI structural connectivity abnormalities were quantified by comparison to a normative map and healthy controls. We explored whether the resection of maximal abnormalities related to improved surgical outcomes, in both modalities individually and concurrently. Additionally, we suggest how connectivity abnormalities may inform the placement of iEEG electrodes pre-surgically using a patient case study.

Findings: Seizure freedom was 15 times more likely in patients with resection of maximal connectivity and iEEG abnormalities (p=0.008). Both modalities separately distinguished patient surgical outcome groups and when used simultaneously, a decision tree correctly separated 36 of 43 (84%) patients.

Interpretation: Our results suggest that both connectivity and iEEG abnormalities may localise epileptogenic tissue, and that these two modalities may provide complementary information in pre-surgical evaluations.

Funding: This research was funded by UKRI, CDT in Cloud Computing for Big Data, NIH, MRC, Wellcome Trust and Epilepsy Research UK.

Keywords: Epilepsy; diffusion-weighted MRI; intracranial EEG; machine learning; multi-modal analysis; surgical prediction.

Figure 1:. Abnormality calculation pipeline.

Connectivity pipeline is shown in panels A-D and iEEG abnormality pipeline is shown in panels E-H. A) Connectivity matrices were generated for each patient using average FA between each pair of regions for the Lausanne-60 brain atlas. B) Connectivity matrices were harmonized across the two sites using ComBat. Known biological effects, age and sex, were regressed out. C) Each connection in each patient was z-scored against healthy controls to get connection abnormalities. D) Connection abnormalities involving each region were averaged (mean) to obtain region-level connectivity abnormalities. E) For each patient, 70s of interictal iEEG recording were analysed. F) The relative band power was calculated for five frequency bands for each electrode contact. G) The relative band power was computed for each region by averaging contacts assigned to that region. H) iEEG abnormalities were calculated for a region by z-scoring the relative band power in each frequency band to a normative map and taking the maximum abnormality.

Figure 2:. Connectivity and iEEG abnormalities in two example patients.

The brain plots show connectivity (left) and iEEG abnormalities (right) for the regions implanted with iEEG electrodes in two example patients (Panel A: good outcome; Panel B: poor outcome). Regions are coloured depending on whether they were resected (red) or spared (blue) in surgery. Increasing point size relates to increasing abnormality. The same information is shown on plots with the x-axis indicating whether a region was resected or spared and the y-axis indicating abnormality size.

Figure 3:. Patients with resection of maximal abnormalities were more likely to be seizure-free.

A) The iEEG and connectivity abnormalities from two example patients. Each point is either a resected (red) or spared (blue) region. Abnormalities were not correlated in all patients. For example, Patient 1 had some high-connectivity and low-iEEG abnormalities (single red arrow), and some low-connectivity and high-iEEG abnormalities (double red arrows). Within each patient, a SVM separated the regions into resected and spared zones based on the size of abnormalities. If the top right of the plot was in the resected (red) zone, then that patient had maximal abnormalities resected. The SVM successfully separated the abnormalities into two zones in 28 out of 43 patients. B) For these patients, surgical outcome was related to resection of maximal abnormalities.

Figure 4:. Connectivity and iEEG abnormality distribution in resected versus spared tissue explains post-surgical seizure freedom.

Both A) connectivity D_RS and B) iEEG D_RS were used to separate patients based on surgical outcome. The top plots in each panel show regional abnormalities, indicated with circular points, in example patients. The bottom plots in each panel show patient D_RS values, indicated with diamond points. C) A decision tree was fit to both D_RS values simultaneously to classify patient outcome, achieving an accuracy of 84%.

Figure 5:. Incorporating dMRI into pre-surgical evaluations - a case study.

A) The connectivity abnormalities are shown for an example patient (iEEG D_RS = 0.76, connectivity D_RS = 0.11). Regions with a black triangle indicate that the region was sampled by iEEG implantation and a red outline indicates that region was resected in surgery. Despite a relatively widespread implantation, the regions with the largest connectivity abnormalities were not implanted, but were resected. This patient was seizure-free following surgery. B) Incorporating connectivity abnormalities from diffusion-weighted imaging into the pre-surgical evaluation may allow for more targeted iEEG implantation or avoiding the need for implantation altogether.