Andrew P. Bagshaw

Combining EEG and fMRI: A multimodal tool for epilepsy research

Patients with epilepsy often present in their electroencephalogram (EEG) short electrical potentials (spikes or spike‐wave bursts) that are not accompanied by clinical manifestations but are of important diagnostic significance. They result from a population of abnormally hyperactive and hypersynchronous neurons. It is not easy to determine the location of the cerebral generators and the other brain regions that may be involved as a result of this abnormal activity. The possibility to combine EEG recording with functional MRI (fMRI) scanning opens the opportunity to uncover the regions of the brain showing changes in the fMRI signal in response to epileptic spikes seen in the EEG. These regions are presumably involved in the abnormal neuronal activity at the origin of epileptic discharges. This paper reviews the methodology involved in performing such studies, particularly the challenge of recording a good quality EEG inside the MR scanner while scanning is taking place, and the methods required for the statistical analysis of the combined EEG and fMRI time series. We review the results obtained in patients with different types of epileptic disorders and discuss the difficult theoretical problems raised by the interpretation of an increase (activation) and decrease (deactivation) in blood oxygen level dependent (BOLD) signal, both frequently seen in response to spikes. J. Magn. Reson. Imaging 2006.

EEG-fMRI of focal epileptic spikes: Analysis with multiple haemodynamic functions and comparison with gadolinium-enhanced MR angiograms

Combined EEG‐fMRI has recently been used to explore the BOLD responses to interictal epileptiform discharges. This study examines whether misspecification of the form of the haemodynamic response function (HRF) results in significant fMRI responses being missed in the statistical analysis. EEG‐fMRI data from 31 patients with focal epilepsy were analysed with four HRFs peaking from 3 to 9 sec after each interictal event, in addition to a standard HRF that peaked after 5.4 sec. In four patients, fMRI responses were correlated with gadolinium‐enhanced MR angiograms and with EEG data from intracranial electrodes. In an attempt to understand the absence of BOLD responses in a significant group of patients, the degree of signal loss occurring as a result of magnetic field inhomogeneities was compared with the detected fMRI responses in ten patients with temporal lobe spikes. Using multiple HRFs resulted in an increased percentage of data sets with significant fMRI activations, from 45% when using the standard HRF alone, to 62.5%. The standard HRF was good at detecting positive BOLD responses, but less appropriate for negative BOLD responses, the majority of which were more accurately modelled by an HRF that peaked later than the standard. Co‐registration of statistical maps with gadolinium‐enhanced MRIs suggested that the detected fMRI responses were not in general related to large veins. Signal loss in the temporal lobes seemed to be an important factor in 7 of 12 patients who did not show fMRI activations with any of the HRFs. Hum. Brain Mapp. 22:179–192, 2004.

Negative BOLD responses to epileptic spikes

Simultaneous electroencephalogram/functional magnetic resonance imaging (EEG‐fMRI) during interictal epileptiform discharges can result in positive (activation) and negative (deactivation) changes in the blood oxygenation level‐dependent (BOLD) signal. Activation probably reflects increased neuronal activity and energy demand, but deactivation is more difficult to explain. Our objective was to evaluate the occurrence and significance of deactivations related to epileptiform discharges in epilepsy. We reviewed all EEG‐fMRI studies from our database, identified those with robust responses (P = 0.01, with ≥5 contiguous voxels with a |t| > 3.1, including ≥1 voxel at |t| > 5.0), and divided them into three groups: activation (A = 8), deactivation (D = 9), and both responses (AD = 43). We correlated responses with discharge type and location and evaluated their spatial relationship with regions involved in the “default” brain state (Raichle et al. [ 2001 ]: Proc Natl Acad Sci 98:676–682]. Deactivations were seen in 52/60 studies (AD+D): 26 related to focal discharges, 12 bilateral, and 14 generalized. Deactivations were usually distant from anatomical areas related to the discharges and more frequently related to polyspike‐ and spike‐and‐slow waves than to spikes. The “default” pattern occurred in 10/43 AD studies, often associated with bursts of generalized discharges. In conclusion, deactivations are frequent, mostly with concomitant activation, for focal and generalized discharges. Discharges followed by a slow wave are more likely to result in deactivation, suggesting neuronal inhibition as the underlying phenomenon. Involvement of the “default” areas, related to bursts of generalized discharges, provides evidence of a subclinical effect of the discharges, temporarily suspending normal brain function in the resting state. Hum Brain Mapp, 2005.

Effect of sleep stage on interictal high-frequency oscillations recorded from depth macroelectrodes in patients with focal epilepsy.

Purpose:  To investigate the effect of sleep stage on the properties of high‐frequency oscillations (HFOs) recorded from depth macroelectrodes in patients with focal epilepsy.

Electrical impedance tomography of human brain function using reconstruction algorithms based on the finite element method

Electrical impedance tomography (EIT) is a recently developed technique which enables the internal conductivity of an object to be imaged using rings of external electrodes. In a recent study, EIT during cortical evoked responses showed encouraging changes in the raw impedance measurements, but reconstructed images were noisy. A simplified reconstruction algorithm was used which modelled the head as a homogeneous sphere. In the current study, the development and validation of an improved reconstruction algorithm are described in which realistic geometry and conductivity distributions have been incorporated using the finite element method. Data from computer simulations and spherical or head-shaped saline-filled tank phantoms, in which the skull was represented by a concentric shell of plaster of Paris or a real human skull, have been reconstructed into images. There were significant improvements in image quality as a result of the incorporation of accurate geometry and extracerebral layers in the reconstruction algorithm. Image quality, assessed by blinded subjective expert observers, also improved significantly when data from the previous evoked response study were reanalysed with the new algorithm. In preliminary images collected during epileptic seizures, the new algorithm generated EIT conductivity changes which were consistent with the electrographic ictal activity. Incorporation of realistic geometry and conductivity into the reconstruction algorithm significantly improves the quality of EIT images and lends encouragement to the belief that EIT may provide a low-cost, portable functional neuroimaging system in the foreseeable future.

Temporal and Extratemporal BOLD Responses to Temporal Lobe Interictal Spikes

Summary:  Purpose: Simultaneous EEG and functional MRI (fMRI) allows measuring metabolic changes related to interictal spikes. Our objective was to investigate blood oxygenation level–dependent (BOLD) responses to temporal lobe (TL) spikes by using EEG‐fMRI recording.

BOLD changes occur prior to epileptic spikes seen on scalp EEG

This study examined BOLD changes prior to interictal discharges in the EEG of patients with epilepsy. From a database of 143 EEG-fMRI studies, we selected the 16 data sets that showed both strong fMRI activation in the original analysis and only a single spike type in the EEG. Scans were then analyzed using seven model HRFs, peaking 3 or 1 s before the event, or 1, 3, 5, 7, or 9 s after it. An HRF was calculated using a deconvolution method for all activations seen in each analysis. The results showed that seven data sets had HRFs that peaked 1 s after the event or earlier, indicating a BOLD change starting prior to the spike seen on the scalp EEG. This is surprising since the BOLD change is expected to result from the spike. For most of the data sets with early peaking HRFs, the maximum activation in all of the statistical maps was when the model HRF peaked 1 s after the event, suggesting that the early activation was at least as important as any later activation. We suggest that this early activity is the result of neuronal changes occurring several seconds prior to a surface EEG event, but that these changes are not visible on the scalp. This is the first report of a BOLD response occurring several seconds prior to an interictal event seen on the scalp and could have important implications for our understanding of the generation of epileptic discharges.

Hemodynamic and metabolic responses to activation, deactivation and epileptic discharges.

To investigate the coupling between the hemodynamic and metabolic changes following functional brain activation as well as interictal epileptiform discharges (IEDs), blood oxygenation level dependent (BOLD), perfusion and oxygen consumption responses to a unilateral distal motor task and interictal epileptiform discharges (IEDs) were examined via continuous EEG-fMRI. Seven epilepsy patients performed a periodic (1 Hz) right-hand pinch grip using approximately 8% of their maximum voluntary contraction, a paradigm previously shown to produce contralateral MI neuronal excitation and ipsilateral MI neuronal inhibition. A multi-slice interleaved pulsed arterial spin labeling and T(2)*-weighted gradient echo sequence was employed to quantify cerebral blood flow (CBF) and BOLD changes. EEG was recorded throughout the imaging session and reviewed to identify the IEDs. During the motor task, BOLD, CBF and cerebral metabolic rate of oxygen consumption (CMR(O(2))) signals increased in the contra- and decreased in the ipsilateral primary motor cortex. The relative changes in CMR(O(2)) and CBF were linearly related, with a slope of 0.46 +/- 0.05. The ratio of contra- to ipsilateral CBF changes was smaller in the present group of epilepsy patients than in the healthy subjects examined previously. IEDs produced both increases and decreases in BOLD and CBF signals. In the two case studies for which the estimation criteria were met, the coupling ratio between IED-induced CMR(O(2)) and CBF changes was estimated at 0.48 +/- 0.17. These findings provide evidence for a preserved coupling between hemodynamic and metabolic changes in response to both functional activation and, for the two case studies available, in response to interictal epileptiform activity.

Functional Connectivity of the Posteromedial Cortex

As different areas within the PMC have different connectivity patterns with various cortical and subcortical regions, we hypothesized that distinct functional modules may be present within the PMC. Because the PMC appears to be the most active region during resting state, it has been postulated to play a fundamental role in the control of baseline brain functioning within the default mode network (DMN). Therefore one goal of this study was to explore which components of the PMC are specifically involved in the DMN. In a sample of seventeen healthy volunteers, we performed an unsupervised voxelwise ROI-based clustering based on resting state functional connectivity. Our results showed four clusters with different network connectivity. Each cluster showed positive and negative correlations with cortical regions involved in the DMN. Progressive shifts in PMC functional connectivity emerged from anterior to posterior and from dorsal to ventral ROIs. Ventral posterior portions of PMC were found to be part of a network implicated in the visuo-spatial guidance of movements, whereas dorsal anterior portions of PMC were interlinked with areas involved in attentional control. Ventral retrosplenial PMC selectively correlated with a network showing considerable overlap with the DMN, indicating that it makes essential contributions in self-referential processing, including autobiographical memory processing. Finally, ventral posterior PMC was shown to be functionally connected with a visual network. The paper represents the first attempt to provide a systematic, unsupervised, voxelwise clustering of the human posteromedial cortex (PMC), using resting-state functional connectivity data. Moreover, a ROI-based parcellation was used to confirm the results.

Intrinsic epileptogenicity in polymicrogyric cortex suggested by EEG-fMRI BOLD responses

Polymicrogyria (PMG) is a widespread cortical malformation frequently associated with seizures and EEG spikes. Its epileptogenicity is poorly understood. Nine patients with simultaneous EEG and fMRI were studied to assess the blood oxygenation level-dependent response to spikes. Sixteen of 18 studies showed responses, with maximum activation involving the lesion in 61.5%, but often limited to a small fraction of that lesion, suggesting intrinsic epileptogenicity in small areas of the PMG cortex.