ORIGINAL RESEARCH article
Altered fMRI connectivity dynamics in temporal lobe epilepsy might explain seizure semiology
- 1Department of Neurology and Brain Imaging Center, Goethe-University Frankfurt am Main, Frankfurt am Main, Germany
- 2Department of Neurology, University Hospital Schleswig-Holstein, Christian-Albrechts-University Kiel, Kiel, Germany
- 3National Hospital for Neurology and Neurosurgery, University College London, London, UK
- 4The Epilepsy Society, Chalfont St. Peter, UK
Temporal lobe epilepsy (TLE) can be conceptualized as a network disease. The network can be characterized by inter-regional functional connectivity, i.e., blood oxygen level-dependent (BOLD) signal correlations between any two regions. However, functional connectivity is not constant over time, thus computing correlation at a given time and then at some later time could give different results (non-stationarity). We hypothesized (1) that non-stationarities can be induced by epilepsy (e.g., interictal epileptic activity) increasing local signal variance and that (2) these transient events contribute to fluctuations in connectivity leading to pathological functioning, i.e., TLE semiology. We analyzed fMRI data from 27 patients with TLE and 22 healthy controls focusing on EEG-confirmed wake epochs only to protect against sleep-induced connectivity changes. Testing hypothesis (1), we identified brain regions where the BOLD signal variance was significantly greater in TLE than in controls: the temporal pole – including the hippocampus. Taking the latter as the seed region and testing hypothesis (2), we calculated the time-varying inter-regional correlation values (dynamic functional connectivity) to other brain regions and found greater connectivity variance in the TLE than the control group mainly in the precuneus, the supplementary and sensorimotor, and the frontal cortices. We conclude that the highest BOLD signal variance in the hippocampi is highly suggestive of a specific epilepsy-related effect. The altered connectivity dynamics in TLE patients might help to explain the hallmark semiological features of dyscognitive seizures including impaired consciousness (precuneus, frontal cortex), sensory disturbance, and motor automatisms (sensorimotor cortices, supplementary motor cortex). Accounting for the non-stationarity and state-dependence of functional connectivity are a prerequisite in the search for potential connectivity-derived biomarkers in TLE.
Epilepsy affects the brain both during seizures and interictally in the form of neurobehavioral problems. These are considered to be due to progressive structural and functional changes in the brain related to syndrome-specific “network” variations (1). In temporal lobe epilepsy (TLE), for example, interictal language and memory impairment are typical (2). In terms of seizure semiology, motor automatisms and consciousness impairment are characteristics of dyscognitive (complex partial) seizures typical for TLE. We hypothesized that, during seizure-free wakeful rest, brain activity is altered in syndrome-specific regions owing to so-called interictal epileptic activity – which may or may not be visible on scalp EEG (3) – and additionally any other sort of paroxysmal activity patterns differing from activity within a healthy brain, an “epileptic process.” Hence, as a first step, we tested the basic assumption that in TLE patients any such paroxysmal process should induce bursts of the blood oxygen level-dependent (BOLD) signal leading to increased regional variance compared to healthy controls. Among the primary candidate regions, we expected either one or both of those contributing to “interictal” cognitive impairment and those functionally related to temporal lobe seizure semiology, i.e., not only temporal (e.g., hippocampal) but also extratemporal (consciousness subserving and motor) regions.
Temporal lobe epilepsy can be considered as a network disease (4, 5) because widespread anatomic abnormality exists outside the primary epileptic zone affecting inter-regional connectivity linked to distributed cognitive impairments (6). Hence, a natural second step is to explore to which other brain regions – exhibiting “normal” signal fluctuations – the locations of increased variance are more tightly connected in the epileptic condition than in the healthy control subjects. By always taking healthy brain activity as the reference, we ensured that the observed typical effects were specific for the epileptic condition, in particular TLE.
However, it has been shown that functional connectivity is not static (7–11) but differs between different stages of wakefulness (12), and in addition fluctuates within each stage over time. “Dynamic functional connectivity” is here understood in the context of fluctuating, non-constant coupling between brain areas when the coupling is computed over short temporal windows. In contrast with the traditional practice of assuming a steady correlation between regions, the study of dynamic connectivity aims to study how connectivity evolves in time, how stable connections are and what the possible short-term connectivity motifs are between regions. Assuming a paroxysmal pathological process influencing brain activity in TLE, it is conceivable that accounting for the non-stationarity of inter-regional BOLD signal correlations both between and within wake and other sleep stages is mandatory for obtaining specific results with high sensitivity. Because of this, we selected epochs of data recorded during epochs of wakefulness and within these based all inferences on dynamic functional connectivity measures as obtained by means of sliding window analysis. Then, the variability of dynamical functional connectivity was computed (as the variance of the functional connectivity time series), which represents the intensity of the temporal fluctuations in the connectivity of the regions. A low value of variance of the dynamic connectivity time series corresponds to stable connections, in which coupling does not deviate much from the average value. In contrast, large values of the variance correspond to widely fluctuating connectivity, i.e., switching over time between high and low values of the coupling between regions.
Materials and Methods
All procedures were subject to the relevant local and national research ethics committees’ approval.
Subjects and Patients
Data from 22 healthy subjects (11 females) served as control in this study (age 35 ± 12 years), and data from 27 epilepsy patients [13 classified as left TLE (age 32 ± 9 years), 14 with right (age 33 ± 11 years); 16 females], who were selected from a larger dataset with the inclusion criteria of TLE based on electro-clinical information-informed expert classification (Table 1).
In this study, we assumed the presence of an “epileptic process” possibly reaching beyond epochs of interictal epileptic activity visible on scalp EEG. Accordingly, the EEG was not used for the detection of epileptiform EEG activity but to identify fMRI epochs of wakefulness in order to increase the sensitivity and specificity of our results (12).
EEG via a cap (modified Brain Cap MR, Easycap, Herrsching, Germany) was recorded continuously (sampling rate 5 kHz, low pass filter 1 kHz) during fMRI acquisition with an MRI-compatible EEG system (BrainAmp MR+ and Brain Vision Analyzer; Brain Products GmbH, Gilching, Germany) yielding two 20-min sessions consisting of 404 T2*-weighted single-shot gradient-echo echo-planar images (EPIs; echo time/repetition time, 30/3000 ms; flip angle, 90°; 43 2.5 mm interleaved slices; FOV, 24 cm × 24 cm; matrix 64 × 64) acquired continuously on a 3T Signa Excite HDX MRI scanner (General Electric, Milwaukee, WI, USA).
EEG Data Preprocessing
MRI and pulse artifact correction were performed based on the average artifact subtraction (AAS) method (13, 14) as implemented in Vision Analyzer 2 (BrainProducts, Germany) resulting in EEG with a sampling rate of 250 Hz. EEG was re-referenced to common average. Sleep stages were scored manually by an expert according to the AASM criteria (15). Epochs other than wakefulness were erased from the analysis (by excluding the corresponding BOLD time courses from both the variance and dynamic functional connectivity analyses). Epochs of wakefulness shorter than 2 min were not included in the analysis.
fMRI Data Preprocessing
Using Statistical Parametric Mapping (SPM 8, http://www.fil.ion.ucl.ac.uk/spm), EPI data were realigned, normalized (MNI space), and spatially smoothed (Gaussian kernel, 5 mm3 full width at half maximum). The time course of the average signal at the ventricles (CSF as given by FSL’s Ventricle Mask standard image) as well as motion-induced noise was regressed out. We did not remove the global brain signal to avoid the issue of induced anti-correlations (16). Also, since there are direct electrophysiological correlates of the global resting-state signal (17), we believe that its removal is arbitrary and results in a loss of information. fMRI data was band pass filtered in the range 0.01–0.1 Hz using a sixth order Butterworth filter as described previously (8). In all cases, both 20 min fMRI sessions were analyzed and then the results were averaged within every subject (to guarantee the assumption of statistical independence) prior to statistical testing.
Subjects displayed similar degrees of head displacement across groups. Controls: 0.06 ± 0.02 mm, Left TLE: 0.08 ± 0.04 mm, Right TLE: 0.08 ± 0.05 mm. No significant differences (p > 0.05) between groups were found. As an additional control, we reproduced the results after erasing volumes associated with high ( >0.3) head displacement amplitude (see Appendix).
Computation of temporal variability
The temporal variability (variance) of BOLD signals carries important information about the brain state and also about task performance. An optimal value is typically observed during rest, which is decreased locally during task performance and increased in other brain states such as sleep (11, 18, 19). The possibility of interictal spikes occurring in TLE patient leads naturally to the hypothesis of locally increased variability (i.e., increased variance) of BOLD activity time courses. To quantify precisely the changes in variability observed in the patient vs. the control population, we use the traditional definition of variance in the time domain:
where X is the BOLD signal, σ2 the variance, and <, > denotes temporal averaging. This analysis is performed on a voxel-wise basis, resulting in spatial maps in which every voxel contains the variance of its associated time series. The procedure is summarized in Figure 1A.
Figure 1. Procedure to study BOLD signal and functional connectivity temporal variability. (A) To compute the voxel-wise map of BOLD signal variability, the time series for every voxel in the EPI data is first extracted. Then, the variance of the signal is computed, resulting in the required map. (B) To compute the voxel-wise map of BOLD functional connectivity temporal variability with a seed, a region in the AAL template is first selected (in this work, the left and right hippocampi are used as seeds) and the average signal in the region is computed. Then, for every voxel in the EPI data, correlations over time are obtained using a sliding window (30 s) and the temporal dynamics of functional connectivity are computed. Note that in this example the dynamics are non-constant with moments of drastic loss of connectivity between regions. Finally, the temporal variance of the functional connectivity time series is computed, resulting in a spatial map (see Figure 4) encoding the variability in the interaction between the seed region and every voxel.
Computation of functional connectivity variability
The study of functional connectivity between brain regions usually neglects the possibility of changes of connectivity occurring over time [dynamical functional connectivity (8)]. Not only functional connectivity values computed over extended periods of time characterize the healthy, resting brain but also a normal switching of connectivity values, representing the engagement and disengagement of different brain networks. At a first level, the normal level of this fluctuation in connectivity can be computed in two steps, which are also summarized graphically in Figure 1B:
(1) A temporal time course of functional connectivity is computed by means of a sliding window analysis. Inside each window, the normal correlation coefficient between the signal at a seed region and the signal at every voxel in the brain is computed. Then, the window is displaced one time step and the analysis repeated, obtaining one connectivity estimate per time unit. In this work, the sliding window length was set to 30 s, and the seed region to the left and right hippocampus [as defined in the automated anatomical labeling atlas (20)].
(2) Once the temporal time courses of functional connectivity are obtained, they are collapsed into a single spatial map by computing their temporal variability or variance (using Eq. 1). Thus, this results in a map in which every voxel has a value encoding how widely over time its connectivity with a seed region of interest is fluctuating.
Both maps of temporal variability and of dynamical functional connectivity variability are compared between groups by means of mass univariate t-tests as implemented in the SPM8 software. All maps are reported at a level of p < 0.001 uncorrected with only clusters passing a threshold of p < 0.05 FWE corrected being shown.
On average, healthy controls slept more than left and right TLE patients (Figure 2). Controls were awake 71 ± 11% of the time, left TLE patients 85 ± 8%, and the right TLE patients 86 ± 9% of the time. This justifies the detection and deleting of BOLD data corresponding to sleep, considering the possibility of confounds due to comparing groups of subjects in different vigilance states (12).
Figure 2. Distribution of percent time spent in different sleep stages (wake, N1 and N2; no N3 sleep was observed) for healthy controls, left and right TLE patients (mean ± SEM).
The voxel-wise analysis of BOLD signal variance in TLE was greater than in controls in the anterior temporal lobe bilaterally (Table 2; Figure 3). Significantly reduced variance compared to controls was not found in any voxel. A significant difference between left TLE and right TLE was not observed.
Table 2. Regions, hemisphere, and statistical significance of areas with increased BOLD signal variance in TLE patients compared to healthy controls (local maxima at least 8 mm apart).
Figure 3. Spatial map (coronal, sagittal, axial slices) of significantly greater variance of the blood oxygen level-dependent signal in patients with temporal lobe epilepsy (pooled right and left) than in healthy controls. Color bar indicates p-value (thresholded at p < 0.001 for display, cluster survives family-wise error correction at p < 0.05). Left on figure is right in the brain (coronal and axial slices).
Dynamic functional connectivity from the hippocampi in TLE (left and right TLE pooled) was greater than in controls in the supplementary motor area, the pre- and postcentral gyri, the (pre-)cuneus, the calcarine/middle occipital gyrus, and the superior frontal gyrus (Tables 3 and 4; Figure 4). A significant difference between left TLE and right TLE was not observed (compare Tables 3 and 4 for seeding in the left and right hippocampus, respectively). Of the pooled analysis, only results for the left hippocampus as the seed region are displayed and differences between seeding in the left and the right hippocampus were not significant (compare Tables 3 and 4).
Table 3. Regions, hemisphere, and statistical significance of areas with increased variance of dynamic functional connectivity with the left hippocampus in patients compared to healthy controls.
Table 4. Regions, hemisphere, and statistical significance of areas with increased variance of dynamic functional connectivity with the right hippocampus in patients compared to healthy controls.
Figure 4. Spatial map (coronal, sagittal, axial slices) of significantly greater variance of hippocampal dynamic functional connectivity in patients with temporal lobe epilepsy (here: seed in left hippocampus) than in healthy controls. For differences between seeding in the left hippocampus and in the right, see Tables 3 and 4, of the pooled analysis, only results for the left hippocampus as the seed region are displayed (seed in right hippocampus yielded most similar results). Color bar indicates p-value (thresholded at p < 0.001 for display, cluster survives family-wise error correction at p < 0.05). Left on figure is right in the brain (coronal and axial slices).
Interictally, compared to healthy controls and during wakefulness, in the TLE group we found increased BOLD signal variance in the anterior temporal lobe regions overlapping with the left and right hippocampi and that these regions showed increased dynamic functional connectivity most prominently to sensory motor structures, the precuneus, and superior frontal cortices.
The identification of anterior temporal lobe structures in particular including the hippocampi is highly suggestive of a syndrome-specific effect because these regions have long been known as crucial in TLE (21). The anterior hippocampus is thought to be highly epileptogenic (22). In the context of interictal neuropsychological dysfunction, it is essential for memory processes (23), and with respect to seizure control, its surgical removal has been shown to be associated with a good outcome (24). In addition, our study design in which we try to keep constant as many parameters as possible between the control and the patient group (including the state of wakefulness) ensures that any significant group differences are specific to the condition TLE. In particular, we can identify significant activity alterations in sets of brain regions constituting intrinsic connectivity networks, which can be present both in the patient and the control group and which might evade other analysis techniques such as independent component analysis or conventional seed region-based functional connectivity analysis.
Interpretation of Increased Hippocampal BOLD Signal Variance
The first observation that in our group of TLE patients, we found the highest BOLD signal variance in the anterior temporal lobe including, of note, the hippocampus supports our hypothesis of the presence of paroxysmal interictal activity resulting in regionally increased BOLD signal variance. It is well-established that local field potentials (LFP) explain the largest portion of this variance (25) and that interictal epileptic spikes in the anterior temporal lobe are reflected in LFP changes (26). We hence propose interictal epileptic activity (“epileptic process”) as the most likely explanation for the high anterior temporal lobe variance common to our TLE group. Kobayashi and colleagues as well as ourselves previously identified hippocampal BOLD signal increases at the group level in TLE, when correlating by means of a general linear model (GLM) interictal scalp EEG activity with the BOLD signal (27, 28). For these GLM-based EEG-fMRI studies, the occurrence of a number of interictal spikes was essential, whereas in the present approach, we on purpose designed our analysis to be independent of scalp EEG interictal activity. We did this for two reasons: (1) the sensitivity of scalp EEG in TLE is compromised with respect to hippocampal interictal epileptic activity (3, 29), allowing the creation of a scalp EEG-based GLM with limited sensitivity only and (2) we sought to identify an fMRI marker for paroxysmal activity specific to the epileptic condition independently of the simultaneously recorded EEG. Such paroxysmal activity in addition to “true” interictal spikes (as can be demonstrated with intracranial EEG) might include other types of epileptic (neuronal) activity leading to bursts of energy consuming processes, e.g., owing to damaged cells in the hippocampus, alterations in the glial milieu, or pathology of blood flow parameters (30). The exclusive selection of epochs during which the patients and control subjects, respectively, were awake ensured increased sensitivity when comparing resting-state BOLD signal properties: it is well known that the functional architecture of the brain changes significantly depending on the level of wakefulness (12, 31) including the variance of the BOLD signal in different brain regions (11, 32), but also graph theoretical network measures (33, 34), and variance in the brain’s functional connectivity structure (8).
We note that a very similar result (increased variance of BOLD signals in bilateral hippocampi) was recently found for subjects under the influence of a psychedelic substance (psilocybin) (35). Furthermore, a related result is discussed in the context of REM sleep (36), and in the psilocybin study the intensity of the “dream-like” experience, which was rated subjectively by the subjects, correlated positively with the magnitude of the hippocampal BOLD signal fluctuations. It can be speculated that subjective alterations in the conscious awareness of the patients – more subtle than full-blown seizures, but which can escalate into the former and also into auras – are shared in TLE and in the other altered states of consciousness mentioned above.
Interpretation of Increased Hippocampal Dynamic Functional Connectivity
Increased variance of dynamic functional connectivity in TLE compared to healthy controls indicates region pairs, which exhibit stronger fluctuations of the correlation of their BOLD time series in the epileptic condition than in the healthy. While the measure does not assess an absolute difference in connectivity between the two cohorts, regions are identified, which link and unlink typically more in TLE than in the control group. Our study design motivates the hypothesis that the “epileptic process” is responsible for this coupling and uncoupling, possibly in the form of the interictal epileptiform activity in the hippocampus, which leads to increased functional connectivity variance. Further analyses should reveal whether such epileptic activity weakens or strengthens the functional connection between regions. It may even be that “erratic” hippocampal coupling (and uncoupling) with cortical regions disrupts the physiological interplay of the latter with yet further parts of the brain – including subcortical and that this eventually results in the observed dysfunctionality.
Van Paesschen and colleagues used single photon emission computed tomography (SPECT) to study patients with hippocampal sclerosis (HS) and observed both ictal hyper- and hypoperfusion in (lateralized) regions partly overlapping with those we describe here in our study, which also includes many cases of HS. Hyperperfusion was described in temporal lobe, middle frontal and central regions, and in the frontal lobes and the precuneus hypoperfusion was found (37). It hence appears that both a “gain” as well as a “loss of function” in different brain regions is linked with ictal dysfunctioning. For example, Chassagnon and colleagues studied patients with mesial TLE and found ictal–interictal hypoperfusion in the posterior cingulate and prefrontal regions, which might be interpreted as a loss of function in the sense of impaired consciousness (precuneus) and executive functioning (frontal regions) (38).
Looking at the particular brain regions in TLE with increased connectivity variance, we speculate that (1) the hippocampal activity interferes with language (39) and memory function (40), both interictally and ictally and that (2) the increased dynamic connectivity to the precuneus and frontal cortex is ictally associated with impaired consciousness (4, 41) and executive functioning (42). We also propose that (3) the increased dynamics in functional connectivity between the hippocampus and the sensorimotor cortices might pave the way for ictal sensory and motor dysfunction and – probably tightly related to the supplementary motor area – in particular motor automatisms (43). The superior frontal gyrus was also described relevant for introspection (44).
How are Interictal Connectivity Changes Linked with Ictal Semiology?
Our study does not offer any objective clues as to what the changes we describe in the TLE group interictally have to do with ictal changes in brain function. While it is well described that language and memory impairment are present interictally (45), sensorimotor dysfunction does not obviously occur interictally. One explanation could be that a qualitative difference between what is called interictal and ictal activity may not exist as such, but rather a quantitative one: Binnie a good decade ago pointed out that if only we tested carefully enough, transitory cognitive impairment could be related to “interictal” activity in many individuals (46). In our own EEG-fMRI study looking at BOLD signal changes related to formally interictal activity, we interpreted signal changes in regions of the so-called default mode network to explain reduced consciousness during dyscognitive seizures. However, we selected for the study cohort patients with very frequent interictal discharges on the EEG increasing the sensitivity of our discharge-correlated GLM analysis (27). With this in mind and taking Binnie’s idea forward, it is conceivable that with so-called interictal activity vastly the same set of brain regions (network) is recruited as is ictally with the difference that when behavioral alterations become obvious they are called seizures and hence define “ictal activity.” Of course, some additional features distinguish seizure from interictal activity going beyond “duration” alone but include spreading of epileptic activity. Still, such spreading of activity might occur along interictally pre-existing paths (47). Supporting our speculation further, in TLE, structural changes have been shown to progress over time and memory function was more closely related to structural hippocampal changes than the overt seizure frequency: the group of Bernasconi found neocortical thinning in TLE progressing over time in bilateral frontal (sensorimotor) and temporal (hippocampal, entorhinal, temporo-polar) regions – overlapping with the regions we report here (48), and Pacagnella and colleagues most recently presented data proposing that memory impairment is more influenced by hippocampal damage than by seizure frequency (40). In our limited cohort, we failed to identify a correlation between epilepsy duration and BOLD signal variance or dynamic functional connectivity (analysis not shown).
Unidirectional Differences and Lack of Significant Lateralization of Our Findings
We did not find any significantly lateralized or side-dependent results. Instead, significant differences between healthy subjects and controls were usually bilateral. This might be due to a lack of sensitivity of our analysis, which was not optimized for this purpose (balancing of handedness, type of left- and right-sided pathologies, and EEG abnormalities). We hence do not discuss our lateralized results. In general, however, on first sight, a lack of lateralized findings is surprising taking into account clinical practice and surgical success with unilateral resections and the efforts spent with non- and invasive video-EEG telemetry and imaging to identify in which temporal lobe (hippocampus) the epileptogenic zone is located (49). Although it is clinically not a contested issue that it is relevant indeed to operate on the correct side of the brain, we are not aware of any systematic review of epilepsy surgery cases in which – for whatever reason – retrospectively the wrong side was operated upon. It is likely that the vastly symmetrical results we present reflect secondary bilateral “network” effects of a lateralized primary cause.
It is equally interesting that we found exclusively variance increases – and not any decreases – of BOLD signal amplitude and hippocampal functional connectivity in patients with TLE compared to control subjects. However, we would like to note that this does not rule out the possibility of decreased absolute functional inter-regional connectivity in the epilepsy cohort, as this differs from the variance of functional connectivity, which we report in our study and the interpretation of which is discussed above. Still, a relationship between the two measures might exist because many of the regions in which we found altered dynamic functional connectivity were also reported by Haneef and colleagues. They compared TLE patients to controls and found that the classical measure of static hippocampal functional connectivity was greater to the bilateral temporal lobes, insula, fornix, frontal poles, angular gyrus, basal ganglia, thalamus, and cerebellum. They found reduced connectivity with the occipital pole, calcarine, lingual, precuneus, sensorimotor cortex, and parts of insula and frontal lobes as well as medial frontal areas (50). We report results of a “Standard Seed Correlation Analysis” of our data in the Appendix. Because we analyzed wakefulness epochs exclusively, comparability to other data remains limited.
We are aware of the many factors influencing resting-state – and any other – fMRI studies (51) but not aware of any study formally assessing the order of relevance of the numerous confounding effects. We do know, though, from our own data (12, 31) that sleep alters the neuronal resting-state brain architecture significantly; and motion is known to introduce BOLD signal changes of several magnitudes the size of those commonly induced neuronally (52). In comparison, e.g., effects of sex and age are less pronounced both in terms of extent and intensity, and only optimized analysis methods will reveal such differences (53, 54). Nevertheless, we tried to match both gender and age as much as possible between the examined cohorts. We regressed rigid body motion from the data and in an additional analysis accounting for motion-induced variance in a very conservative way showed that our results were robust w.r.t. motion-induced variance (see “Reproduction of Results Erasing High Amplitude Head Movements” in the Appendix; Figures 5 and 6).
Figure 5. Spatial map (coronal, sagittal, axial slices) of significantly greater variance of the blood oxygen level-dependent signal in patients with temporal lobe epilepsy (pooled right and left) than in healthy controls. An additional preprocessing step was performed by erasing volumes associated with large head displacements, as well as surrounding volumes. Color bar indicates p-value (thresholded at p < 0.001 for display, cluster survives family-wise error correction at p < 0.05). Left on figure is right in the brain (coronal and axial slices).
Figure 6. Spatial map (coronal, sagittal, axial slices) of significantly greater variance of hippocampal dynamic functional connectivity in patients with temporal lobe epilepsy (here: seed in left hippocampus) than in healthy controls. An additional preprocessing step was performed by erasing volumes associated with large head displacements, as well as six surrounding volumes. Color bar indicates p-value (thresholded at p < 0.001 for display, cluster survives family-wise error correction at p < 0.05). Left on figure is right in the brain (coronal and axial slices).
Regarding confounding fluctuations in wakefulness, to our knowledge, this is the first “wakeful rest” functional connectivity study controlling for and restricting the analysis to true, EEG-defined, awake epochs only and in addition accounting for non-stationarity of the functional connectivity even within the awake state. We advocate both as a desirable standard. To demonstrate independence of our results to the length of the sliding window used in the dynamic functional connectivity analysis, we performed an additional analysis with a different (shorter) window length (see “Robustness against Changes in Sliding Window Length” in the Appendix, Figure 7).
Figure 7. Spatial map (coronal, sagittal, axial slices) of significantly greater variance of hippocampal dynamic functional connectivity in patients with temporal lobe epilepsy (here: seed in left hippocampus) than in healthy controls. A sliding window length of 15 s was used for the computation of dynamic functional connectivity time series. Color bar indicates p-value (thresholded at p < 0.001 for display, cluster survives family-wise error correction at p < 0.05). Left on figure is right in the brain (coronal and axial slices).
We did not make a strong point about left TLE vs. right TLE comparisons, because we did not observe significant group differences. Evidently, absence of proof is not proof of absence but may be due to patient heterogeneity. Such is almost inevitable in patient studies like this one caused by a variety of factors, such as the range of dyscognitive seizure semiology (not any two patients have the same seizures), seizure frequency, and the potential occurrence of additional generalized tonic clonic seizures, type and location – or by current radiological standards even absence – of structural pathology visible on MRI, IQ, duration of epilepsy, and medication – just to name a few. Although epilepsy syndrome diagnosis relied on multidisciplinary experts reviewing extensive electro-clinical information, we cannot rule out classification errors. But differences within the TLE cohort between left and right TLE patients will not affect our main positive findings, i.e., the group differences between TLE and healthy controls. However, we note that any medication taken by the patients if systematically leading to, e.g., alterations of consciousness could have biased our results; but there is no evidence to date supporting this possibility for the specific patterns we observed.
In addition to these mere technical issues, like in any resting-state study, result interpretation conceptionally is limited given the absence of a task. We based the justification of the study design on the fact that interictal epileptic activity occurs spontaneously at rest as indicated by neuronal discharges measurable with EEG, although we assumed a more general “epileptic process” based on the well-established general clinical observation of interictal cognitive compromitation in TLE. In an attempt to tie our interpretation of our observations more closely to the results, we regressed IQ against BOLD signal variance and dynamic functional connectivity. We report this analysis only in “Correlations Between BOLD Signal Variance/Variance of Dynamic Connectivity Time Series with the Left Hippocampus and VIQ/PIQ” in the Appendix as it needs to be considered of anecdotal character owing to a lack of statistical significance. Causes of the latter include those discussed above.
We found evidence in support of our hypothesis that, interictally, brain activity is altered in syndrome-specific regions. Assuming that interictal processes like bursts of interictal epileptiform discharges will generate large changes in BOLD amplitude, we analyzed the variance of the BOLD signal and found this increased in anterior temporal regions, which suggest a TLE specific effect. Starting from the anterior temporal lobe, we found hippocampal dynamic connectivity increased in regions, which might explain the hallmark semiological features of complex partial seizures including impaired consciousness (precuneus, frontal cortex), sensory disturbance, and motor automatisms (sensorimotor cortices, supplementary motor area). Taking into account state of the art knowledge about the non-stationarity and state-dependence of functional connectivity, we sought to increase the sensitivity and specificity of our results. More generally, this work encourages the further development of connectivity-derived measures as potential functional imaging biomarkers in TLE.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
This work was funded by the Bundesministerium für Bildung und Forschung (grant 01 EV 0703), the Deutsche Forschungsgemeinschaft (grant LA 1452/3-1), and the LOEWE Neuronale Koordination Forschungsschwerpunkt Frankfurt (NeFF). We thank Verena Brodbeck, Astrid Morzelewski, and Sajitha Cannadathu for EEG processing including sleep scoring. We are indebted to Pam Thompson for neuropsychology assessments. We are grateful to the Wolfson Trust and the Epilepsy Society for supporting the Epilepsy Society MRI scanner. This work was undertaken at UCLH/UCL who received a proportion of funding from the Department of Health’s NIHR UCLH Biomedical Research Centre funding scheme. We thank all participants of this study for their contribution.
2. Sidhu MK, Stretton J, Winston GP, Bonelli S, Centeno M, Vollmar C, et al. A functional magnetic resonance imaging study mapping the episodic memory encoding network in temporal lobe epilepsy. Brain (2013) 136:1868–88. doi:10.1093/brain/awt099
8. Tagliazucchi E, Von Wegner F, Morzelewski A, Brodbeck V, Laufs H. Dynamic BOLD functional connectivity in humans and its electrophysiological correlates. Front Hum Neurosci (2012) 6:339. doi:10.3389/fnhum.2012.00339
10. Hutchison RM, Womelsdorf T, Gati JS, Everling S, Menon RS. Resting-state networks show dynamic functional connectivity in awake humans and anesthetized macaques. Hum Brain Mapp (2013) 34:2154–77. doi:10.1002/hbm.22058
11. Tagliazucchi E, Von Wegner F, Morzelewski A, Brodbeck V, Jahnke K, Laufs H. Breakdown of long-range temporal dependence in default mode and attention networks during deep sleep. Proc Natl Acad Sci U S A (2013) 110:15419–24. doi:10.1073/pnas.1312848110
12. Tagliazucchi E, Laufs H. Decoding wakefulness levels from typical fMRI resting-state data reveals reliable drifts between wakefulness and sleep. Neuron (2014) 82:695–708. doi:10.1016/j.neuron.2014.03.020
13. Allen PJ, Polizzi G, Krakow K, Fish DR, Lemieux L. Identification of EEG events in the MR scanner: the problem of pulse artifact and a method for its subtraction. Neuroimage (1998) 8:229–39. doi:10.1006/nimg.1998.0361
16. Murphy K, Birn RM, Handwerker DA, Jones TB, Bandettini PA. The impact of global signal regression on resting state correlations: are anti-correlated networks introduced? Neuroimage (2009) 44:893–905. doi:10.1016/j.neuroimage.2008.09.036
20. Tzourio-Mazoyer N, Landeau B, Papathanassiou D, Crivello F, Etard O, Delcroix N, et al. Automated anatomical labeling of activations in SPM using a macroscopic anatomical parcellation of the MNI MRI single-subject brain. Neuroimage (2002) 15:273–89. doi:10.1006/nimg.2001.0978
22. Bartolomei F, Chauvel P, Wendling F. Epileptogenicity of brain structures in human temporal lobe epilepsy: a quantified study from intracerebral EEG. Brain (2008) 131:1818–30. doi:10.1093/brain/awn111
23. Bonelli SB, Powell RH, Yogarajah M, Samson RS, Symms MR, Thompson PJ, et al. Imaging memory in temporal lobe epilepsy: predicting the effects of temporal lobe resection. Brain (2010) 133:1186–99. doi:10.1093/brain/awq006
26. Demont-Guignard S, Benquet P, Gerber U, Wendling F. Analysis of intracerebral EEG recordings of epileptic spikes: insights from a neural network model. IEEE Trans Biomed Eng (2009) 56:2782–95. doi:10.1109/TBME.2009.2028015
27. Laufs H, Hamandi K, Salek-Haddadi A, Kleinschmidt AK, Duncan JS, Lemieux L. Temporal lobe interictal epileptic discharges affect cerebral activity in “default mode” brain regions. Hum Brain Mapp (2007) 28:1023–32. doi:10.1002/hbm.20323
28. Kobayashi E, Grova C, Tyvaert L, Dubeau F, Gotman J. Structures involved at the time of temporal lobe spikes revealed by interindividual group analysis of EEG/fMRI data. Epilepsia (2009) 50:2549–56. doi:10.1111/j.1528-1167.2009.02180.x
29. Bettus G, Ranjeva JP, Wendling F, Benar CG, Confort-Gouny S, Regis J, et al. Interictal functional connectivity of human epileptic networks assessed by intracerebral EEG and BOLD signal fluctuations. PLoS One (2011) 6:e20071. doi:10.1371/journal.pone.0020071
31. Tagliazucchi E, Von Wegner F, Morzelewski A, Borisov S, Jahnke K, Laufs H. Automatic sleep staging using fMRI functional connectivity data. Neuroimage (2012) 63:63–71. doi:10.1016/j.neuroimage.2012.06.036
32. Fukunaga M, Horovitz SG, Van Gelderen P, De Zwart JA, Jansma JM, Ikonomidou VN, et al. Large-amplitude, spatially correlated fluctuations in BOLD fMRI signals during extended rest and early sleep stages. Magn Reson Imaging (2006) 24:979–92. doi:10.1016/j.mri.2006.04.018
33. Spoormaker VI, Schroter MS, Gleiser PM, Andrade KC, Dresler M, Wehrle R, et al. Development of a large-scale functional brain network during human non-rapid eye movement sleep. J Neurosci (2010) 30:11379–87. doi:10.1523/JNEUROSCI.2015-10.2010
34. Tagliazucchi E, Von Wegner F, Morzelewski A, Brodbeck V, Borisov S, Jahnke K, et al. Large-scale brain functional modularity is reflected in slow electroencephalographic rhythms across the human non-rapid eye movement sleep cycle. Neuroimage (2013) 70:327–39. doi:10.1016/j.neuroimage.2012.12.073
37. Van Paesschen W, Dupont P, Van Driel G, Van Billoen H, Maes A. SPECT perfusion changes during complex partial seizures in patients with hippocampal sclerosis. Brain (2003) 126:1103–11. doi:10.1093/brain/awg108
38. Chassagnon S, Namer IJ, Armspach JP, Nehlig A, Kahane P, Kehrli P, et al. SPM analysis of ictal-interictal SPECT in mesial temporal lobe epilepsy: relationships between ictal semiology and perfusion changes. Epilepsy Res (2009) 85:252–60. doi:10.1016/j.eplepsyres.2009.03.020
40. Pacagnella D, Lopes TM, Morita ME, Yasuda CL, Cappabianco FA, Bergo F, et al. Memory impairment is not necessarily related to seizure frequency in mesial temporal lobe epilepsy with hippocampal sclerosis. Epilepsia (2014) 55:1197–204. doi:10.1111/epi.12691
43. Baumgartner C, Flint R, Tuxhorn I, Van Ness PC, Kosalko J, Olbrich A, et al. Supplementary motor area seizures: propagation pathways as studied with invasive recordings. Neurology (1996) 46:508–14. doi:10.1212/WNL.46.2.508
47. Vulliemoz S, Lemieux L, Daunizeau J, Michel CM, Duncan JS. The combination of EEG source imaging and EEG-correlated functional MRI to map epileptic networks. Epilepsia (2010) 51:491–505. doi:10.1111/j.1528-1167.2009.02342.x
48. Bernhardt BC, Bernasconi N, Concha L, Bernasconi A. Cortical thickness analysis in temporal lobe epilepsy: reproducibility and relation to outcome. Neurology (2010) 74:1776–84. doi:10.1212/WNL.0b013e3181e0f80a
53. Mowinckel AM, Espeseth T, Westlye LT. Network-specific effects of age and in-scanner subject motion: a resting-state fMRI study of 238 healthy adults. Neuroimage (2012) 63:1364–73. doi:10.1016/j.neuroimage.2012.08.004
54. Smith DV, Utevsky AV, Bland AR, Clement N, Clithero JA, Harsch AE, et al. Characterizing individual differences in functional connectivity using dual-regression and seed-based approaches. Neuroimage (2014) 95:1–12. doi:10.1016/j.neuroimage.2014.03.042
55. Lemieux L, Salek-Haddadi A, Lund TE, Laufs H, Carmichael D. Modelling large motion events in fMRI studies of patients with epilepsy. Magn Reson Imaging (2007) 25:894–901. doi:10.1016/j.mri.2007.03.009
56. Power JD, Mitra A, Laumann TO, Snyder AZ, Schlaggar BL, Petersen SE. Methods to detect, characterize, and remove motion artifact in resting state fMRI. Neuroimage (2014) 84:320–41. doi:10.1016/j.neuroimage.2013.08.048
Reproduction of Results Erasing High Amplitude Head Movements
Following the work of Lemieux et al. and Power et al. (55, 56), we have scanned the time course of estimated head movement and erased those volumes associated with a displacement larger than 0.3 mm, as well as the three previous and following volumes. We then reproduced Figures 3 and 4 using this alternative preprocessing. Results are shown in Figure 5 (voxel-wise variance, controls vs. TLE) and Figure 6 (hippocampus-based dynamic functional connectivity variance, controls vs. TLE). It can be appreciated from both figures that the general composition of the patterns discussed in the main text is preserved.
Robustness Against Changes in Sliding Window Length
The computation of dynamic functional connectivity time series requires the specification of a sliding window length, which can be seen as a free parameter (even though its correct choice bears an obvious relationship with the time scale where the changes occur). We demonstrate that results are robust against a chance of this parameter by reproducing Figure 4 using a different, shorter window length (15 s). Results are shown in Figure 7. From this figure, it is evident that the main differences between controls and TLE patients still hold using this different window length.
Standard Seed Correlation Analysis
We performed a standard “static” linear correlation analysis seeded in the left hippocampus. Results of the controls vs. TLE comparison are presented in Figure A1. A decrease in seed connectivity was observed in the patient group relative to the control group in brain areas including visual cortex, precuneus, and sensorimotor cortex. Compared to the regions in which a change in dynamic connectivity variance was observed (Figure 4), these regions overlap but also include a less extended network than that observed in the dynamic analysis; in particular, frontal regions do not appear. To a large extent, for regions having strong baseline connectivity, increased variance of the dynamic connectivity time series implies decreased “static” correlation. The opposite, however, is not true. This points to an origin of the decreased “static” correlation that can be traced to a dynamic transient phenomenon, which we hypothesize corresponds to paroxysmal events.
Figure A1. Spatial map (coronal, sagittal, axial slices) of significantly decreased hippocampal signal linear correlation in patients with temporal lobe epilepsy (here: seed in left hippocampus) than in healthy controls. Color bar indicates p-value (thresholded at p < 0.001 for display, cluster survives family-wise error correction at p < 0.05). Left on figure is right in the brain (coronal and axial slices).
Correlations between BOLD Signal Variance/Variance of Dynamic Connectivity Time Series with the Left Hippocampus and VIQ/PIQ
In an attempt to link the results derived from resting-state data more tightly to our interpretation, we performed a correlation analysis of the behavioral measures PIQ/VIQ (obtained close in time to the fMRI experiment) and fMRI signal changes in the regions of increased BOLD signal variance and increased dynamic functional connectivity, respectively (Figure A2). We found positive correlations between the IQ measures and the BOLD signal variance, and negative correlations with the dynamic BOLD functional connectivity. With the caveat of uncorrected significance values – possibly a power or, alternatively, a conceptual problem – this argues for a relationship of our findings obtained at rest with these psychological measures, i.e., a connection between interictal cognitive impairment and the observed changes between the TLE and the control cohorts. That a lower IQ score is linked to higher functional connectivity could be interpreted as paroxysmal interference of the “epileptic process” with normally required and less pronounced variations in functional connectivity. Why at the same time an impaired IQ is associated with reduced variance in the BOLD signal time series in regions characterized by an overall increase in BOLD signal variance in TLE compared to controls – admittedly – is difficult to interpret. Although the functional role of anterior temporal lobe is still debated (57), it might be possible that increased “inertia” (i.e., reduced BOLD signal variance) of the anterior temporal lobe reflects deficits in semantic memory (58).
Figure A2. Correlations between BOLD signal variance/variance of dynamic connectivity time series with the left hippocampus and VIQ/PIQ scores for the TLE group (left + right). Statistical significance was determined at p < 0.05 with a cluster threshold of 10 voxels. Maps were masked with the regions where the metrics were significantly different between controls and TLE patients. Left: results for PIQ and VIQ correlation with variance. In all cases positive correlations were found. Right: results for PIQ and VIQ correlation with variance of connectivity time series with the left hippocampus. In all cases negative correlations were found.
Keywords: functional connectivity, non-stationarity, temporal lobe epilepsy, biomarker, EEG-fMRI, interictal epileptiform discharges, semiology, seizure
Citation: Laufs H, Rodionov R, Thornton R, Duncan JS, Lemieux L and Tagliazucchi E (2014) Altered fMRI connectivity dynamics in temporal lobe epilepsy might explain seizure semiology. Front. Neurol. 5:175. doi: 10.3389/fneur.2014.00175
Received: 08 July 2014; Accepted: 28 August 2014;
Published online: 11 September 2014.
Edited by:David Vaughan, The Florey Institute of Neuroscience and Mental Health, Australia
Reviewed by:Gaelle Eve Doucet, Thomas Jefferson University Hospital, USA
Zulfi Haneef, Baylor College of Medicine, USA
Copyright: © 2014 Laufs, Rodionov, Thornton, Duncan, Lemieux and Tagliazucchi. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Helmut Laufs, Theodor-Stern-Kai 7, 60590 Frankfurt am Main, Hessen, Germany e-mail: email@example.com