Skip to main content

Face-induced gamma oscillations and event-related potentials in patients with epilepsy: an intracranial EEG study

Abstract

Background

To examine the pathological effect of a mesial temporal seizure onset zone (SOZ) on local and inter-regional response to faces in the amygdala and other structures of the temporal lobe.

Methods

Intracranial EEG data was obtained from the amygdala, hippocampus, fusiform gyrus and parahippocampal gyrus of nine patients with drug-refractory epilepsy during visual stimulation with faces and mosaics. We analyzed event-related potentials (ERP), gamma frequency power, phase-amplitude coupling and phase-slope-index and compared the results between patients with versus without a mesial temporal SOZ.

Results

In the amygdala and fusiform gyrus, faces triggered higher ERP amplitudes compared to mosaics in both patient groups and higher gamma power in patients without a mesial temporal SOZ. In the hippocampus, famous faces triggered higher gamma power for both groups combined but did not affect ERPs in either group. The differentiated ERP response to famous faces in the parahippocampal gyrus was more pronounced in patients without a mesial temporal SOZ. Phase-amplitude coupling and phase-slope-index results yielded bidirectional modulation between amygdala and fusiform gyrus, and predominately unidirectional modulation between parahippocampal gyrus and hippocampus.

Conclusions

A mesial temporal SOZ was associated with an impaired response to faces in the amygdala, fusiform gyrus and parahippocampal gyrus in our patients. Compared to this, the response to faces in the hippocampus was impaired in patients with, as well as without, a mesial temporal SOZ. Our results support existing evidence for face processing deficits in patients with a mesial temporal SOZ and suggest the pathological effect of a mesial temporal SOZ on the amygdala to play a pivotal role in this matter in particular.

Peer Review reports

Background

A mesial temporal seizure onset zone (SOZ) in temporal lobe epilepsy (TLE) has consistently been associated with deficits in social cognition [1,2,3,4]. One of the most basic levels of social cognition is constituted by the ability to recognize a face as a face, e.g. as a prerequisite to recognize and perceive face emotions [5]. Typically, the evaluation of this ability involves measuring face-selective responses in a distributed network of brain areas thought to be involved in face processing summarized as the “face network” [6]. Inside the face network, the brain area most strongly associated with face emotion recognition is the amygdala [7,8,9]. While the amygdala has a long history of being linked with face emotion recognition, its response to neutral faces compared to non-face stimuli has been gaining attention in the more recent years as well. Some studies showed the amygdala’s face-selectivity to be on a similar scale as, or sometimes even exceeding, the fusiform face area [10,11,12]. This face-selective response in the amygdala was suggested to be a key component in the amygdala’s contribution to face processing [13].

In connectivity analyses, the amygdala and hippocampus form a distinct cluster associated with emotional processes, memory formation and face discrimination [14, 15]. The hippocampus not only responds to faces [11, 12, 16] but also shows a differentiated response to familiar versus non-familiar faces [17, 18] and emotional faces [19]. From a pathological standpoint, the epileptogenic hippocampus represents one of the most common causes for focal epilepsy [20]. Inside the face network, evidence exists for the close connectivity between the amygdala and hippocampus to the fusiform and parahippocampal gyrus respectively [21,22,23,24]. The fusiform gyrus with its functionally localized fusiform face area [25] was suggested to receive direct influence from the amygdala in fearful face processing [26], and a face-specific connectivity between them was shown using neutral face stimuli [11]. The role of the parahippocampal gyrus was proposed to be that of ‘contextual processing’, meaning it provides the hippocampus with the information as to “where” and “when” an input was associated with [27]. Multiple studies have shown activation of both the hippocampus and parahippocampal gyrus during associative tasks involving faces (for review, see Ref. [28]).

While intracranially implanted electrodes ensure a high level of spatial and temporal precision, their use remains a feature in a fraction of all EEG studies in humans due to the invasive implantation procedure. Two established methods to analyze intracranial EEG consist of averaging the post-stimulus epochs to compute event-related potentials (ERP) and conducting a time–frequency analysis [18, 29,30,31]. In time–frequency analyses inside the face network, gamma frequencies have been studied most extensively out of all frequency bands [32, 33]. It has been proposed that synchronization in the gamma frequencies affect communication between neuronal groups and thus represent a fundamental process in neuronal connections subserving higher cognitive functions [34]. A number of studies compared ERPs and gamma power elicited by faces in selective sites of the visual network [35,36,37,38]. The two methods show different response characteristics to faces, and it has been speculated that gamma power is associated with elaborate processing of faces, while ERPs may represent synchronization within regions of the face processing network [39]. So far it is unclear if and how a SOZ in the examined area affects ERPs and gamma frequencies differently. Most intracranial EEG studies aim to examine healthy brain functionality and involve patients with a SOZ outside of their regions of interest only [40,41,42]. To our knowledge, no study examining face processing in the temporal lobe has compared epileptogenic versus non-epileptogenic functionality using intracranial EEG.

It was therefore the goal of our study to evaluate the effect of a SOZ on the EEG response to faces inside the temporal lobe. To accomplish this, intracranial ERPs and gamma power were measured in response to faces versus mosaics, and between different categories of faces, in the amygdala, hippocampus, fusiform and parahippocampal gyrus in patients with drug-refractory epilepsy. Different categories of neutral faces, including famous public figures and wooden masks as an abstract representation of faces, were chosen to test for the robustness of the response to faces. A linear mixed-effects model was implemented to compare the results between patients with a mesial temporal SOZ, Group M, and patients without a mesial temporal SOZ, Group O. As fMRI studies examining face processing in patients with a mesial temporal SOZ demonstrated decreased activity in the amygdala, hippocampus and parahippocampal gyrus on the side of seizure onset [43,44,45] and decreased functional connectivity in face-processing regions of the temporal lobe centered around the amygdala [46], we hypothesized to see a reduced ERP and gamma frequency response to faces in the examined regions of Group M compared to Group O.

In a second step, we analyzed the functional connectivity between the amygdala, hippocampus, fusiform gyrus, and parahippocampal gyrus. It has been proposed that gamma frequencies contribute to communication between neuronal groups through cross-frequency-coupling with lower frequency bands [34]. Out of all cross-frequency-couplings, phase-amplitude coupling (PAC) has gained interest in particular [15, 47,48,49]. This is due to the assumption that lower frequency phase modulates the time window in which neuronal spiking activity increases, which is then reflected by gamma frequency amplitude accordingly [50, 51]. Recent research suggests that when PAC is calculated between inter-regional intracranial electrodes, gamma frequency amplitude as an output of the higher order, “driver” area will be related to the low frequency phase of the lower order, “receiver” area [52]. Following this line of thinking, we hypothesized to find significant PAC between gamma frequency amplitude in the amygdala and low frequency (3–20 Hz) phase in the fusiform gyrus, and between gamma frequency amplitude in the hippocampus and low frequency phase in the parahippocampal gyrus. For further analysis of the modulation direction, we calculated the phase-slope-index (PSI), which also considers the time lag that is bound to occur when a signal modulates another in a remote area [53]. As our PSI calculations represent cumulative directionality over a span of 750 ms, we hypothesized to see bidirectional connectivity between the amygdala and fusiform gyrus as reported before using fMRI and Granger causality [54], as well as between the hippocampus and parahippocampal gyrus.

Materials and methods

Patients

Nine patients with drug-refractory epilepsy participated in this study (aged 18 to 54 years, 5 female and 4 male), all of whom were stereotactically implanted with intracranial electrodes for pre-surgical diagnostics independent of this study (see Fig. 1 for examples of electrode localizations). Using intracranial EEG data, patients were divided into two groups according to their SOZ (Table 1). Group O (patients 1 to 4) included patients with an extratemporal SOZ, and Group M (patients 5 to 9) with a mesial temporal SOZ. For bilaterally implanted patients presenting a unilateral SOZ, the unaffected contralateral hemisphere was attributed to Group O (patients 5 and 6). All patients gave their written informed consent to participate in this study. The study protocol was approved by the institutional review board.

Fig. 1
figure 1

Examples of intracranial electrode localizations. Coronal MRI images (T1-weighted) of two patients. Red dots: Electrode contacts inside target structure. Orange dots: Electrode contacts outside target structure. Only data derived from electrode contacts inside target structure were analyzed. A Amygdala electrode contacts in patient 5, right hemisphere (Group M). B Hippocampus electrode contacts in patient 6, right hemisphere (Group O)

Table 1 Patients and clinical data

Experiment procedure

We chose facial portraits in the categories Caucasian faces, dark-skinned faces, faces of famous public figures, veiled faces, and wooden masks of African or oceanic decent, each containing 50 pictures, respectively. Caucasian faces and faces of famous public figures each contained 25 female and 25 male models, while dark-skinned and veiled faces comprised of 50 female models. Caucasian, dark-skinned and veiled faces did not contain any person familiar to the participating study patients. We added wooden masks as a category to evaluate the response to abstract representations of faces. Veiled faces were included to test for differences between whole faces and faces where the mouth and nose regions are hidden. Famous public figures were chosen out of an online database ranking the most famous public figures in each age group [55]. In all pictures, facial expressions range from neutral to smiling, with none of the pictures showing an obvious sign of negative emotion. Additionally, a non-face category of 50 colored mosaics was compiled, bringing the total sum of pictures to 300. All pictures were rated by 11 participants in a separate rating procedure according to valence and arousal using the self-assessment manikin [56], to ensure that the pictures are emotionally neutral. Valence was rated on a scale from 1 for unpleasant to 5 for pleasant and resulted in a mean value of 2.64 with a standard deviation of 0.80. Arousal was rated on a scale from 1 for calm to 5 for excited and resulted in a mean value of 2.83 with a standard deviation of 0.71.

During the visual stimulation experiment, patients were seated in a darkened room, facing a 66 × 44 cm monitor with a face-screen distance of 1 m. The interstimulus interval varied randomly between 2800 and 3000 ms, after which a black screen with a central white cross was shown for 300 ms as a visual focus. Each picture was scaled to 600 × 800 pixels and shown for 1500 ms. Patients were asked to passively view the images.

EEG recording

Electrode localizations were determined using co-registered (Compumedics® Neuroscan™, CURRY 7) pre-implantation 3D-MRI-datasets (3-Tesla, T1-weighted) and post-implantation CT scans. Subsequently each contact was attributed to a respective anatomical structure according to an atlas of anatomy [57]. Only electrode contacts implanted into one of our four regions of interest were selected for data analysis: Amygdala, hippocampus, fusiform gyrus and parahippocampal gyrus.

EEG recordings were sampled at 1024 or 2048 Hz. All EEG data and statistical analyses were carried out using Matlab R2011a. Data sampled at 2048 Hz was down-sampled to 1024 Hz.

Event-related potentials

ERPs were computed by averaging 1500 ms post-stimulus epochs and subtracting the mean of a pre-stimulus baseline of 1000 ms. Surface electrode Fz was used as the reference. A bandpass-filter between 0.5 and 20 Hz was applied. To remove outliers, epochs with a standard deviation greater than 2.5 times mean standard deviation across all epochs were rejected. This removal of outliers was performed for each electrode contact in each patient. Overall, the mean number of epochs that went into the statistical analyses per patient per electrode localization per category was 82 with a mean standard deviation of 13.1. Expected peak locations were defined to be around 110, 240 and 360 ms after stimulus onset based on visual inspection of the computed ERPs and information from published studies with a study design similar to ours: We mainly followed the ERP categorization of Barbeau et al. [18], who conducted intracranial ERPs during face recognition in multiple areas inside the temporal lobe including three out of four of our own regions of interest. This study identified an early ERP component peaking at 110 ms after visual stimulation in the fusiform gyrus, followed by stages of widespread parallel processing in different areas of the visual network at 240 and 360 ms post stimulus. In a different study, intracranial ERPs in the amygdala elicited a first negative peak at around 250 ms when viewing faces [36]. We hypothesize this peak to be comparable to the N240 described by Barbeau et al. [18]. For each of the above-mentioned peak locations, mean ERP amplitude was calculated in a window of 80 ms, which were marked in grey in Fig. 2. All graphs were plotted to show the ERP pointing upwards for easier viewing, i.e., multiplied by -1. All amplitudes mentioned below were adjusted accordingly.

Fig. 2
figure 2

Event-related potentials. Grey areas mark the windows in which ERP amplitude was measured. Windows containing a significant categorical difference are marked with an asterisk (*). A ERPs for faces versus mosaics. In AM, N240 was higher for faces than mosaics in both groups, while N360 was higher only in Group O. Both groups showed higher N240 for faces in FU and higher N360 for faces in FU and PH. B ERPs for all face categories. In AM of Group O, N360 was higher for famous faces than wooden masks and veiled faces. In AM of Group M, N360 was higher for famous and dark-skinned faces than Caucasian faces, while N110 was lower for wooden masks than dark-skinned or veiled faces. In HI of Group M, N360 was lower for wooden masks than dark-skinned faces. In FU, wooden masks triggered lower amplitudes than all other face categories for N240 in both groups and N360 in Group O. In PH, N240 was higher for famous faces than dark-skinned faces and wooden masks in Group O. N360 was higher for famous faces than all other face categories in Group O and higher than wooden masks, Caucasian and veiled faces in Group M

Gamma frequency analysis

Gamma frequency analysis was computed using a bipolar reference. The same method of artefact rejection using standard deviation as in the ERP calculation was applied here as well. We subtracted the ERP from each epoch to remove any effect that is phase-locked to the ERP. A time–frequency representation was computed using fast Fourier transform and wavelet convolution with complex Morlet wavelets in 40 logarithmically increasing frequencies between 1 and 200 Hz. Baseline normalization was performed via conversion into decibel using a pre-stimulus period of 1100 ms. We then extracted the mean power between 45 and 150 Hz during 1–750 ms (early-onset gamma) and 250 – 1000 ms (late-onset gamma) and applied the linear mixed-effects model. In order to find significant clusters of gamma power, we first subtracted the time–frequency representation for mosaics from that of faces (Fig. 3 A3–D3, A4–D4). We then compared the face minus mosaic time–frequency representation between Group O and Group M by performing a two-tailed two-sample t-test for each pixel of the time–frequency representation to create a map of test statistic values (Fig. 3A5–D5). Then a permutation test was performed by randomly dividing all patients into two groups and computing the same aforementioned two-sample t-test, which was repeated for 3000 times. In each permuted map of t-statistics, as well as in the real map of t-statistics, all voxels beneath a threshold corresponding to a p-value of 0.05 were removed. The sum of t-statistics in each remaining cluster was calculated and the maximum and minimum sum of each permuted map registered. We determined the 95th percentile from the distribution of maximum values and the 5th percentile from the distribution of minimum values. In the real map of t-statistics, all clusters with a sum of t-statistics above the 5th and beneath the 95th percentile were removed (Fig. 3A6–D6) [58, 59].

Fig. 3
figure 3

Gamma frequency analysis. A1D1, A2D2 Power map in response to faces and mosaics in Group O and Group M. Mean gamma power was measured during 0–750 ms in AM, FU and PH and 250–1000 ms in HI after stimulus onset and applied to the linear mixed-effects model. Gamma power was higher for faces than mosaics in AM (A1) and FU (C1) of Group O. A3D3, A4D4 Power map showing the response to faces minus mosaics in Group O and Group M. A5D5 Unthresholded t-statistics comparing the power map for faces minus mosaics in Group O and Group M. Positive t-statistics indicate higher power in Group O and negative t-statistics indicate higher power in Group M. A6D6 Remaining clusters of t-statistics after voxel-based and cluster-based thresholding at p ≤ 0.05. Large clusters in the gamma frequency range remained in AM and FU, indicating higher gamma power in Group O

Linear Mixed-Effects Models and Statistical Analyses

All statistical analyses were carried out using R. Using the lme4 package [60], we built two linear mixed-effects models with ERP amplitude defined as the dependent variable and differing fixed effects. The fixed effects for Model 1 were group (O and M), picture category (face and mosaic), hemisphere (left and right), peaks (N110, N240 and N360) and the interaction between group and picture category. For each of the four electrode localizations, all trials by all patients conducted in the respective electrode localization were pooled together and entered into Model 1. The fixed effects for Model 2 were group, all face categories (Caucasian, dark-skinned, famous, veiled and masks), hemisphere, peaks and the interaction between group and face categories. As we were only interested in the face categories in Model 2, we removed all trials containing the mosaics category from the pooled trials and entered them into this model. Both models included patient as a random effect to account for non-independence of trials and were fit by restricted maximum likelihood. Visual inspection of residual plots was used to ensure that assumptions on linearity, homogeneity of variance and normality of residuals were met. P-values were calculated using the lmerTest package [61], which applies Satterthwaite’s method to generate p-values. If the effect of the interaction between group and category was significant, we followed up with post hoc pairwise comparisons using the emmeans package [62] with built-in Tukey method to adjust for multiple comparisons. For each electrode localization and each peak, we also conducted pairwise comparisons for picture or face category using the emmeans package. To evaluate gamma power, we used both Model 1 and 2 without peaks as a fixed effect. All p-values excluding those already adjusted in the post hoc analyses were corrected for multiple comparisons using the Holm-Bonferroni method. The corrected p-values were considered significant at p < 0.05. Confidence intervals reported under Results were not corrected for multiple comparisons.

Phase-amplitude coupling and phase-slope-index

PAC between low frequency phase (3–20 Hz) and gamma frequency amplitude (45–200 Hz) was computed for a chosen range of electrode pairs between (1) the amygdala and fusiform gyrus, (2) the hippocampus and parahippocampal gyrus, and (3) the amygdala and hippocampus. Eight out of nine patients were implanted in both amygdala and fusiform gyrus, six patients in both hippocampus and parahippocampal gyrus, and seven patients in both amygdala and hippocampus. For each of these patients, we visually inspected the time frequency representation of every electrode inside the aforementioned structures, to choose one electrode with the least amount of artifacts for each structure. Inter-regional PAC was then calculated for each patient’s electrode pair by multiplying the gamma frequency amplitude time series of one structure with the low frequency (3–20 Hz) phase time series of the other, and vice versa, over 750 ms after stimulus presentation [50]. PAC results were converted into Z-scores and considered significant at z ≥ 2.

Following this, we computed PSI for all electrode pairs showing a significant amount of PAC. PSI was calculated between low frequency phase and gamma power envelope using cross frequency coherence and a segment length of 750 ms [63]. Only the PSI results within the phase and amplitude range revealed as significant according to each electrode pair’s respective PAC results were considered. All PSI results were converted into Z-scores and considered significant at z ≥ 2 or ≤ − 2.

Results

For an overview of the statistical results for ERP and gamma frequency analysis, see Tables 2 and 3.

Table 2 Linear mixed-effects model results for ERP and gamma frequency analysis
Table 3 Linear mixed effects model post-hoc analysis for ERP peaks

Responsive versus non-responsive electrode contacts

Each of the 51 electrode contacts in Group O and 53 electrode contacts in Group M were examined individually. In Group O, 40 (78%) electrode contacts had an ERP in response to faces, 27 (53%) had enhanced gamma power in response to faces, and 26 (51%) had both ERPs and gamma power. In Group M, ERP was seen in 44 (83%), gamma power in 23 (43%), and both ERPs and gamma power in 22 (42%) electrode contacts in response to faces. Combining both groups together, two (4%) out of 50 electrode contacts with face-induced gamma power did not show a face-induced ERP, while 36 (43%) out of 84 electrode contacts with face-induced ERP did not show face-induced gamma power.

Response to faces versus mosaics

In the amygdala, faces triggered higher ERP amplitudes than mosaics, which was visible in N240 in both groups and N360 in Group O only. For gamma amplitude, the linear mixed-effects model did not show an effect of picture category for both groups combined, but the interaction between group and picture category was significant. To study this interaction effect, pairwise comparisons were performed, which revealed faces to trigger higher gamma power compared to mosaics in Group O (β = 0.70, 95% CI [0.44, 0.96], p = 0.017; Fig. 3A1) but not in Group M (β = − 0.11, 95% CI [− 0.45, 0.23], p > 0.05; Fig. 3A2). The gamma cluster analysis confirmed this result by displaying large clusters in the gamma frequency range with positive test statistic values indicating higher gamma power in Group O than Group M (Fig. 3A6). The earliest face-induced gamma cluster was shown to be around 200 ms after stimulus onset. These results combined demonstrate that, while both groups show a differentiated response to faces in the amygdala, Group O’s response is stronger on a statistically significant level. From the interaction effect between group and picture category, which was significant in the linear mixed-effects model for gamma power but not ERPs, we can also conclude that this difference between Group O and Group M was more pronounced in the gamma analysis.

The fusiform gyrus exhibited the most prominent difference between faces and mosaics out of all examined regions. ERP amplitude was higher in response to faces in N240 and N360 in both groups. While the linear mixed-effects model yielded an interaction effect between group and picture category, pairwise comparisons showed faces to trigger higher ERPs than mosaics in Group O (β = 14.4, 95% CI [11.6, 17.3], p < 0.001) as well as Group M (β = 6.25, 95% CI [4.42, 8.08], p < 0.001). Comparing the estimated marginal means for faces and mosaics in each group thus revealed that the interaction effect stems from Group O showing an even higher difference between faces and mosaics than Group M. The linear mixed-effects model for gamma power analysis also displayed a difference between faces and mosaics for both groups combined. An interaction effect between group and picture category could not be shown after correcting for multiple comparisons. As the p-value was significant before applying the Holm-Bonferroni method, pairwise comparisons were still performed. These resulted in higher gamma power in response to faces in Group O (β = 0.90, 95% CI [0.30, 1.49], p < 0.001; Fig. 3C1) but not in Group M (β = − 0.05, 95% CI [− 0.59, 0.49], p > 0.05; Fig. 3C2). In accordance to this, our gamma cluster analysis demonstrated higher gamma power in Group O than Group M, which lasted until almost 1000 ms after stimulus onset (Fig. 3C6).

In the hippocampus, neither ERP amplitude, nor gamma power could be shown to give rise to a distinct response to faces compared to mosaics, with no differences between Group O and Group M either. In the parahippocampal gyrus, N360 amplitude was higher for faces in both groups, while gamma power analysis did not reveal a difference between faces and mosaics in either group. Gamma cluster analysis also did not display any significant clusters comparing faces and mosaics in neither the hippocampus (Fig. 3B6), nor the parahippocampal gyrus (Fig. 3D6). But visual inspection of the gamma power plots presented an important difference between the hippocampus and parahippocampal gyrus. While gamma power in the hippocampus was almost non-existent for both picture categories and groups (Fig. 3B1, B2), the parahippocampal gyrus displayed a noticeable gamma response to both faces and mosaics in both groups (Fig. 3D1, D2). The fact that the response to mosaics was just as pronounced as the response to faces in the parahippocampal gyrus, led to the lack of a measurable difference between the two stimuli there.

Face categories

Conclusively, the linear mixed-effects models revealed ERP amplitudes to differ between face categories in each examined region. When we studied the category differences across all regions, two face categories stood out. Firstly, famous faces were recorded to trigger higher ERPs compared to other categories most often, i.e., N240 in the fusiform gyrus of both groups, N360 in the amygdala and parahippocampal gyrus of both groups, N240 in the parahippocampal gyrus of Group O, and N360 in the fusiform gyrus of Group O. While the hippocampus was the only region where this distinct ERP response to famous faces was not measured, it was the only region that did show a difference between face categories in the gamma analysis. Pairwise comparisons revealed gamma power in the hippocampus to be highest for famous faces and significantly higher compared to dark-skinned faces (β = 0.52, 95% CI [0.06, 0.98], p = 0.017), wooden masks (β = 0.49, 95% CI [0.05, 0.92], p = 0.022) and veiled faces (β = 0.43, 95% CI [0.01, 0.85], p = 0.049). There was no interaction effect between group and face categories, meaning that this distinct reaction to famous faces can be attributed to both groups. Secondly, wooden masks triggered lower ERPs more often than other categories, i.e., N240 in the fusiform gyrus of both groups, N360 in the parahippocampal gyrus of both groups, N360 in the fusiform gyrus of Group O, N240 in the parahippocampal gyrus of Group O, N360 in the hippocampus of Group M, and N110 in the amygdala of Group M. Gamma power, on the other hand, was not shown to be lower for wooden masks compared to other categories.

PAC and PSI

Between the amygdala and the fusiform gyrus, three out of eight electrode pairs showed PAC between low frequency phase in the fusiform gyrus and gamma amplitude in the amygdala, while showing no PAC in the reverse constellation (i.e. low frequency phase in the amygdala and gamma amplitude in the fusiform gyrus; see Fig. 4A1–4 for an example). One out of eight electrode pairs showed PAC in both of the aforementioned constellations.

Fig. 4
figure 4

Phase-amplitude coupling and phase-slope index. A1A4 PAC between AM and FU in patient 6, right hemisphere (Group O). Enhanced PAC (Z-score ≥ 2) was observed between FU theta/alpha phase and AM low gamma frequency amplitude when viewing faces. No enhanced PAC was observed when viewing mosaics. B Overlap between PSI and the area of enhanced PAC. Red area indicates modulation directionality from FU theta phase towards AM low gamma amplitude. Blue area indicates modulation directionality from AM low gamma amplitude towards FU alpha phase. C1–4 PAC between HI and PH in patient 9, right hemisphere (Group M). Enhanced PAC (Z-score ≥ 2) was calculated between PH alpha/beta phase and HI high gamma amplitude when viewing faces. No enhanced PAC was observed when viewing mosaics. D Overlap between PSI and the area of enhanced PAC. Red area indicates modulation directionality from PH alpha/beta phase towards HI high gamma amplitude. Blue areas indicate modulation directionality from HI high gamma amplitude towards PH alpha phase

Between the hippocampus and the parahippocampal gyrus, three out of six electrode pairs showed PAC between low frequency phase in the parahippocampal gyrus and gamma amplitude in the hippocampus, with no coupling in the reverse constellation (i.e., low frequency phase in the hippocampus and gamma amplitude in the parahippocampal gyrus) (see Fig. 4C1–4 for an example).

Between the amygdala and hippocampus, one out of seven electrode pairs showed PAC between low frequency phase in the amygdala and gamma amplitude in the hippocampus, with no coupling in the reverse constellation.

The remaining electrode pairs showed no enhanced PAC.

Overall, eight out of nine significant PAC computations resulted in considerably smaller to non-existent couplings for mosaics compared to faces.

PSI between low frequency phase in the fusiform gyrus and gamma amplitude in the amygdala was calculated for four electrode pairs, with phase and amplitude range specified according to each pair’s respective PAC. All four electrode pairs showed both positive and negative PSI values, indicating bidirectional modulation between low frequency phase in the fusiform gyrus and gamma amplitude in the amygdala (see Fig. 4B for an example).

Between low frequency phase in the parahippocampal gyrus and gamma amplitude in the hippocampus, PSI was calculated for three electrode pairs. While all three pairs revealed negative PSI values, only one of them additionally showed a positive PSI value as well (see Fig. 4D for an example). Negative PSI values indicate modulation directionality from gamma amplitude in the hippocampus towards low frequency phase in the parahippocampal gyrus, while positive PSI values indicate the opposite direction.

PSI between low frequency phase in the amygdala and gamma amplitude in the hippocampus was calculated for the single electrode pair showing enhanced PAC and revealed both positive and negative values, indicating bidirectional modulation.

Discussion

ERP and gamma frequency analysis

We set out to compare the response to faces in specific regions of the temporal lobe between patients with and without a mesial temporal SOZ. Conclusively, we saw predominantly preserved, non-epileptogenic functionality in the amygdala, fusiform and parahippocampal gyrus in patients without a mesial temporal SOZ (Group O), while patients with a mesial temporal SOZ (Group M) displayed impaired functionality in all of the examined regions, which we will illustrate in the following.

Previous studies examining non-epileptogenic amygdala functionality have shown faces to trigger a negative intracranial ERP peak comparable to N240 and a subsequent lower peak comparable to N360 [8, 36], as well as intracranial gamma power activity starting from 200 ms after stimulus onset in the amygdala [36]. Compared to this, our amygdala results can be summarized as follows: While Group O presented both higher N240 and N360 for faces, as well as the first face-induced gamma cluster around 200 ms, Group M was only able to replicate a higher N240 for faces compared to mosaics. Sato et al. [36] speculated this response to faces starting from 200 ms in the amygdala to be related to conscious perceptual processing of faces, referencing the “visual awareness negativity”, a prominent ERP peak observed around 200–300 ms in scalp EEG for consciously detected visual stimuli only [64, 65]. As Group O was able to replicate non-epileptogenic amygdala functionality while Group M was not, we interpret from these results that conscious perception of faces may be impaired in patients with a mesial temporal SOZ. Multiple behavioral studies examining humans and primates associate impaired amygdala functionality with impaired perception of faces: It was reported that monkeys with bilateral amygdala lesions had no viewing preference for faces or face-like objects and also exhibited changed face-viewing behavior, i.e. no advantage for the eyes or mouth region [66]. This lack of fixation on the eye region was also seen in a human patient with bilateral amygdala damage [67]. In a patient with unilateral amygdala damage, gaze shift towards the eye and mouth regions was impaired in brief (150 ms) stimulus presentations, while it was normal during longer (5000 ms) presentations [68]. These results indicate that the amygdala is involved in the rapid detection of facial features, which in turn could explain why patients with a mesial temporal SOZ were reported to score lower than healthy controls in tests involving the recognition of face emotions and familiar faces [4, 69, 70]. Studies on unfamiliar face recognition on the other hand revealed conflicting results [71, 72], demonstrating the necessity for conducting further electrophysiological studies.

The hippocampal ERPs of our patients, while presenting visible N240 and N360 peaks, did not show a distinction between faces and mosaics. Although Halgren et al. [73] in their intracranial EEG study recorded a face-specific N240 in a portion of their hippocampal electrode contacts and a non-specific N360, they and also Barbeau et al. [18] described a later positive peak, P480, to be the dominant ERP component in the hippocampus. Contrary to this, visual inspection of our own ERPs implies that P480 was not a prominent factor in the hippocampi of our patients. When individual face categories were compared, none of the ERP peaks showed a distinction to famous faces either, contrary to results from studies examining non-epileptogenic hippocampal functionality [18, 30].

In the fusiform gyrus, ERPs altogether closely replicated the sequence of N110-P160-N240 defined by Barbeau et al. [18], which is also comparable to the sequence of N130-P180-N240 in the study by Halgren et al. [73]. Both studies not only observed the N110-P160-N240 in the fusiform gyrus in response to faces but attributed the role of the principal generator of the sequence to the fusiform gyrus, exerting a causal influence on other regions of the temporal lobe. In particular, they postulated the N110 to represent a rapid feed-forward signal triggering further processing: Halgren et al. [73] theorized the N110-P160-N240 to represent a bottom-up sequence with the fusiform gyrus receiving a “primitive sketch” from the lower-level visual cortices during N110 and projecting face-specific encoding to higher cortical areas during P160 and N240. P160 possibly corresponds to N170, a well-studied ERP component that multiple source localization studies reported to originate from face-selective regions in the fusiform gyrus and occipitotemporal cortex (for review, see [74]). While N170 or P160 was repeatedly demonstrated to be larger in response to faces compared to non-face stimuli [75, 76], P160 in the fusiform gyrus of our patients revealed lower amplitudes for faces than mosaics in both groups, for which the reason is unclear. In our gamma analysis, both groups exhibited early face-induced gamma power (Fig. 3C1, C2), but the gamma cluster analysis indicated higher gamma power for Group O that lasted until > 500 ms after stimulus onset (Fig. 3C6), similar to results from studies exploring healthy fusiform gyrus functionality [77, 78].

ERPs in the parahippocampal gyrus, while not replicating the clear triphasic N240-P300-N360 observed by Barbeau et al. [18], did show higher N360 in response to faces than mosaics in both groups. Interestingly, the parahippocampal gyrus in Group O showed the most distinguished ERP response to famous faces out of all regions examined. Studies using fMRI and MEG have shown the parahippocampal cortex to be associated with extracting familiar contextual associations from visual stimuli [79], with an increase in phase synchrony between the parahippocampal cortex and other regions of the contextual associations network starting at 150–250 ms after stimulus onset [80]. Combined with these insights, our results suggest that the parahippocampal N240 and N360 are associated with the presence of contextual information of a face. Furthermore, our connectivity analyses detected PAC between low frequency phase in the parahippocampal gyrus and gamma frequency amplitude in the hippocampus, while a directed modulation of the hippocampus by the parahippocampal gyrus could not be shown. This reflects a feedforward projection from the parahippocampal gyrus to the hippocampus during face processing and also suggests that the source of the impaired response to famous faces seem to be rooted in the hippocampus, as opposed to the parahippocampal gyrus.

Thus, we can conclude that a mesial temporal SOZ was associated with an impaired ERP and gamma power response to faces in the amygdala, fusiform gyrus and parahippocampal gyrus in our patients. While this part of our hypothesis was confirmed, we did not expect both our patient groups to display an impaired response to faces in the hippocampus. In a similar fashion, Halgren et al. [73] noted that a large number of their hippocampal contacts did not show distinct ERPs and the ones who did displayed extremely variable waveforms, while amygdala contacts generally showed clear ERPs even in patients with no hippocampal ERPs. They attributed the lack of hippocampal ERPs to structural damage arising from hippocampal sclerosis. When we break down the patients in Group O, one third of the patients had a SOZ ipsilateral extratemporal, and another third in the contralateral mesial temporal lobe in the form of hippocampal sclerosis. Evidence for structural damage in the hippocampus exists for both of these etiologies: In patients with extratemporal epilepsy with normal hippocampal MRI, magnetic resonance spectroscopic imaging could demonstrate the presence of structural damage in the hippocampus [81], while previously T2 relaxometry had revealed abnormalities [82]. Evidence also exists for structural anomalies in the contralateral hippocampus at the presence of unilateral mesial-TLE [83, 84]. A possible explanation for the involvement of the contralateral hippocampus is seizure propagation through the dorsal hippocampal commissure [85, 86], although the functional relevance of this structure is disputed [87, 88]. These findings could explain why hippocampal EEG responses were similarly impaired in both groups of our study.

A limitation of our study is the fact that two patients are incorporated in both patient groups, as the hemisphere containing the SOZ was attributed to Group M and the unaffected hemisphere to Group O. Our linear mixed-effects model accounts for this by assigning a different random intercept for each patient in addition to the fixed effects of group and hemisphere. This way the non-independence of observations is taken into consideration. Nevertheless, the limited number of patients prevented us from being able to make group comparisons for PAC and PSI results. Further research involving a higher number of patients is necessary for better assessment of the pathological effect of a SOZ on the connectivity of specific temporal lobe regions during face processing.

PAC and PSI

Due to the restricted number of electrode pairs that we examined for our PAC and PSI calculations, we refrained from making comparative statements between Group O and M. We expected to find significant PAC between gamma frequency amplitude in the amygdala and low frequency phase in the fusiform gyrus, and between gamma frequency amplitude in the hippocampus and low frequency phase in the parahippocampal gyrus. This turned out to be the case in the majority of the examined electrode pairs, indicating feedback modulation from the amygdala to the fusiform gyrus, and from the hippocampus to the parahippocampal gyrus. Keeping the interpretational differences between local and inter-regional PAC in mind, our PAC results are still in accordance with the theory of “predictive coding” [89], which says that brain areas simultaneously communicate through lower and higher frequencies representing feedback predictions and feedforward prediction errors. Our results indicate that the amygdala and the hippocampus represent “driver” areas that send feedback predictions to the fusiform gyrus and the parahippocampal gyrus, which represent “receiver” areas and in turn send feedforward prediction errors.

The PSI results in the hippocampus and parahippocampal gyrus predominately indicate unidirectional modulation of the parahippocampal gyrus by the hippocampus. In the amygdala and the fusiform gyrus, however, PSI results indicate bidirectional modulation in all examined electrode pairs. Previously, bidirectional connectivity between the amygdala and fusiform gyrus was shown using fMRI and Granger analysis [54] as well as dynamic causal modeling [23] during visual processing. Considering our PSI results represent cumulative results over 750 ms, a strictly unidirectional modulation seems less plausible than bidirectional. Therefore, this discrepancy possibly represents further evidence for impaired hippocampal functionality in our patients.

Conclusions

In summary, a mesial temporal SOZ was associated with an impaired ERP and gamma power response to faces in the amygdala, fusiform gyrus and parahippocampal gyrus in our patients. Compared to this, the response to faces in the hippocampus was impaired in patients with, as well as without, a mesial temporal SOZ. Our connectivity analyses showed bidirectional modulation between the amygdala and fusiform gyrus when viewing faces, and predominately unidirectional modulation between the hippocampus and the parahippocampal gyrus, which further suggests overall higher levels of functional impairment in the hippocampi compared to the amygdalae in our patients. Our results support existing fMRI and behavioral studies reporting deficits in face processing in patients with a mesial temporal SOZ, with the added benefit of the spatial and temporal accuracy of intracranial EEG. Our results further suggest the pathological effect of a mesial temporal SOZ on the amygdala to play a pivotal role in this matter in particular.

Availability of data and materials

The raw datasets analyzed for the current study are not publicly available because we did not obtain the consent of participants to provide them to third parties for publication, but they are available from the corresponding author on reasonable request.

Abbreviations

AM:

Amygdala

ERP:

Event-related potential

FU:

Fusiform gyrus

HI:

Hippocampus

PAC:

Phase-amplitude coupling

PH:

Parahippocampal gyrus

PSI:

Phase-slope-index

SOZ:

Seizure onset zone

TLE:

Temporal lobe epilepsy

References

  1. Bala A, Okruszek Ł, Piejka A, Głębicka A, Szewczyk E, Bosak K, et al. Social perception in mesial temporal lobe epilepsy: interpreting social information from moving shapes and biological motion. J Neuropsychiatry Clin Neurosci. 2018;30(3):228–35. https://doi.org/10.1176/appi.neuropsych.17080153.

    Article  PubMed  Google Scholar 

  2. Ives-Deliperi VL, Jokeit H. Impaired social cognition in epilepsy: a review of what we have learnt from neuroimaging studies. Front Neurol. 2019;10:940. https://doi.org/10.3389/fneur.2019.00940.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Schacher M, Winkler R, Grunwald T, Kraemer G, Kurthen M, Reed V, Jokeit H. Mesial temporal lobe epilepsy impairs advanced social cognition. Epilepsia. 2006;47(12):2141–6. https://doi.org/10.1111/j.1528-1167.2006.00857.

    Article  PubMed  Google Scholar 

  4. Broicher SD, Kuchukhidze G, Grunwald T, Krämer G, Kurthen M, Jokeit H. “Tell me how do I feel”–emotion recognition and theory of mind in symptomatic mesial temporal lobe epilepsy. Neuropsychologia. 2012;50(1):118–28. https://doi.org/10.1016/j.neuropsychologia.2011.11.005.

    Article  PubMed  Google Scholar 

  5. Steiger BK, Jokeit H. Why epilepsy challenges social life. Seizure. 2017;44:194–8. https://doi.org/10.1016/j.seizure.2016.09.008.

    Article  PubMed  Google Scholar 

  6. Ishai A. Let’s face it: it’s a cortical network. Neuroimage. 2008. https://doi.org/10.1016/j.neuroimage.2007.10.040.

    Article  PubMed  Google Scholar 

  7. Morris JS, Frith CD, Perrett DI, Rowland D, Young AW, Calder AJ, Dolan RJ. A differential neural response in the human amygdala to fearful and happy facial expressions. Nature. 1996;383(6603):812–5.

    Article  CAS  Google Scholar 

  8. Krolak-Salmon P, Hénaff MA, Vighetto A, Bertrand O, Mauguière F. Early amygdala reaction to fear spreading in occipital, temporal, and frontal cortex: a depth electrode ERP study in human. Neuron. 2004;42(4):665–76.

    Article  CAS  Google Scholar 

  9. Méndez-Bértolo C, Moratti S, Toledano R, Lopez-Sosa F, Martínez-Alvarez R, Mah YH, et al. A fast pathway for fear in human amygdala. Nat Neurosci. 2016;19(8):1041–9. https://doi.org/10.1038/nn.4324.

    Article  CAS  PubMed  Google Scholar 

  10. Young LR, Yu W, Holloway M, Rodgers BN, Chapman SB, Krawczyk DC. Amygdala activation as a marker for selective attention toward neutral faces in a chronic traumatic brain injury population. Neuropsychologia. 2017;104:214–22. https://doi.org/10.1016/j.neuropsychologia.2017.08.026.

    Article  PubMed  Google Scholar 

  11. Mende-Siedlecki P, Verosky SC, Turk-Browne NB, Todorov A. Robust selectivity for faces in the human amygdala in the absence of expressions. J Cogn Neurosci. 2013;25(12):2086–106. https://doi.org/10.1162/jocn_a_00469.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Rossion B, Hanseeuw B, Dricot L. Defining face perception areas in the human brain: a large-scale factorial fMRI face localizer analysis. Brain Cogn. 2012;79(2):138–57. https://doi.org/10.1016/j.bandc.2012.01.001.

    Article  PubMed  Google Scholar 

  13. Todorov A. The role of the amygdala in face perception and evaluation. Motiv Emot. 2012;36(1):16–26. https://doi.org/10.1007/s11031-011-9238-5.

    Article  PubMed  Google Scholar 

  14. Amft M, Bzdok D, Laird AR, Fox PT, Schilbach L, Eickhoff SB. Definition and characterization of an extended social-affective default network. Brain Struct Funct. 2015;220(2):1031–49. https://doi.org/10.1007/s00429-013-0698-0.

    Article  PubMed  Google Scholar 

  15. Zheng J, Anderson KL, Leal SL, Shestyuk A, Gulsen G, Mnatsakanyan L, et al. Amygdala-hippocampal dynamics during salient information processing. Nat Commun. 2017;8:14413. https://doi.org/10.1038/ncomms14413.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Ishai A, Schmidt CF, Boesiger P. Face perception is mediated by a distributed cortical network. Brain Res Bull. 2005;67(1–2):87–93.

    Article  Google Scholar 

  17. Dietl T, Trautner P, Staedtgen M, Vannucci M, Mecklinger A, Grunwald T, et al. Processing of famous faces and medial temporal lobe event-related potentials: a depth electrode study. Neuroimage. 2005;25(2):401–7.

    Article  CAS  Google Scholar 

  18. Barbeau EJ, Taylor MJ, Regis J, Marquis P, Chauvel P, Liégeois-Chauvel C. Spatio temporal dynamics of face recognition. Cereb Cortex. 2008;18(5):997–1009.

    Article  Google Scholar 

  19. Wang S, Tudusciuc O, Mamelak AN, Ross IB, Adolphs R, Rutishauser U. Neurons in the human amygdala selective for perceived emotion. Proc Natl Acad Sci USA. 2014;111(30):E3110–9. https://doi.org/10.1073/pnas.1323342111.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Blumcke I, Spreafico R, Haaker G, Coras R, Kobow K, Bien CG, et al. Histopathological findings in brain tissue obtained during epilepsy surgery. N Engl J Med. 2017;377(17):1648–56. https://doi.org/10.1056/NEJMoa1703784.

    Article  PubMed  Google Scholar 

  21. Burwell RD. The parahippocampal region: corticocortical connectivity. Ann NY Acad Sci. 2000;911:25–42.

    Article  CAS  Google Scholar 

  22. Fairhall SL, Ishai A. Effective connectivity within the distributed cortical network for face perception. Cereb Cortex. 2007;17(10):2400–6.

    Article  Google Scholar 

  23. Herrington JD, Taylor JM, Grupe DW, Curby KM, Schultz RT. Bidirectional communication between amygdala and fusiform gyrus during facial recognition. Neuroimage. 2011;56(4):2348–55. https://doi.org/10.1016/j.neuroimage.2011.03.072.

    Article  PubMed  Google Scholar 

  24. Libby LA, Ekstrom AD, Ragland JD, Ranganath C. Differential connectivity of perirhinal and parahippocampal cortices within human hippocampal subregions revealed by high-resolution functional imaging. J Neurosci. 2012;32(19):6550–60. https://doi.org/10.1523/JNEUROSCI.3711-11.2012.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Kanwisher N, Yovel G. The fusiform face area: a cortical region specialized for the perception of faces. Philos Trans R Soc Lond B Biol Sci. 2006;361(1476):2109–28. https://doi.org/10.1098/rstb.2006.1934.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Vuilleumier P, Richardson MP, Armony JL, Driver J, Dolan RJ. Distant influences of amygdala lesion on visual cortical activation during emotional face processing. Nat Neurosci. 2004;7(11):1271–8. https://doi.org/10.1038/nn1341.

    Article  CAS  PubMed  Google Scholar 

  27. Aminoff EM, Kveraga K, Bar M. The role of the parahippocampal cortex in cognition. Trends Cogn Sci. 2013;17(8):379–90. https://doi.org/10.1016/j.tics.2013.06.009.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Eichenbaum H, Yonelinas AP, Ranganath C. The medial temporal lobe and recognition memory. Annu Rev Neurosci. 2007;30:123–52.

    Article  CAS  Google Scholar 

  29. Niessing J, Ebisch B, Schmidt KE, Niessing M, Singer W, Galuske RAW. Hemodynamic signals correlate tightly with synchronized gamma oscillations. Science. 2005. https://doi.org/10.1126/science.1110948.

    Article  PubMed  Google Scholar 

  30. Trautner P, Dietl T, Staedtgen M, Mecklinger A, Grunwald T, Elger CE, Kurthen M. Recognition of famous faces in the medial temporal lobe—an invasive ERP study. Neurology. 2004. https://doi.org/10.1212/01.WNL.0000140487.55973.D7.

    Article  PubMed  Google Scholar 

  31. Gray CM, König P, Engel AK, Singe W. Oscillatory responses in cat visual cortex exhibit inter-columnar synchronization which reflects global stimulus properties. Nature. 1989;338:334–7.

    Article  CAS  Google Scholar 

  32. Lachaux JP, George N, Tallon-Baudry C, Martinerie J, Hugueville L, Minotti L, et al. The many faces of the gamma band response to complex visual stimuli. Neuroimage. 2005;25(2):491–501. https://doi.org/10.1016/j.neuroimage.2004.11.052.

    Article  PubMed  Google Scholar 

  33. Tsuchiya N, Kawasaki H, Oya H, Howard MA, Adolphs R. Decoding face information in time, frequency and space from direct intracranial recordings of the human brain. PLoS ONE. 2008;3(12): e3892. https://doi.org/10.1371/journal.pone.0003892.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Fries P. Rhythms for cognition: communication through coherence. Neuron. 2015;88(1):220–35. https://doi.org/10.1016/j.neuron.2015.09.034.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Vidal JR, Ossandón T, Jerbi K, Dalal SS, Minotti L, Ryvlin P, et al. Category-specific visual responses: an intracranial study comparing gamma, beta, alpha, and ERP response selectivity. Front Hum Neurosci. 2010;4:195. https://doi.org/10.3389/fnhum.2010.00195.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Sato W, Kochiyama T, Uono S, Matsuda K, Usui K, Inoue Y, Toichi M. Temporal profile of amygdala γ oscillations in response to faces. J Cogn Neurosci. 2012;24(6):1420–33. https://doi.org/10.1162/jocn_a_00142.

    Article  PubMed  Google Scholar 

  37. Engell AD, McCarthy G. Selective attention modulates face-specific induced gamma oscillations recorded from ventral occipitotemporal cortex. J Neurosci. 2010;30(26):8780–6. https://doi.org/10.1523/JNEUROSCI.1575-10.2010.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Fisch L, Privman E, Ramot M, Harel M, Nir Y, Kipervasser S, et al. Neural “ignition”: enhanced activation linked to perceptual awareness in human ventral stream visual cortex. Neuron. 2009;64(4):562–74.

    Article  CAS  Google Scholar 

  39. Engell AD, McCarthy G. The relationship of gamma oscillations and face-specifiC ERPs recorded subdurally from occipitotemporal cortex. Cereb Cortex. 2011;21(5):1213–21. https://doi.org/10.1093/cercor/bhq206.

    Article  PubMed  Google Scholar 

  40. Meletti S, Cantalupo G, Benuzzi F, Mai R, Tassi L, Gasparini E, et al. Fear and happiness in the eyes: an intra-cerebral event-related potential study from the human amygdala. Neuropsychologia. 2012;50(1):44–54. https://doi.org/10.1016/j.neuropsychologia.2011.10.020.

    Article  PubMed  Google Scholar 

  41. Pourtois G, Spinelli L, Seeck M, Vuilleumier P. Temporal precedence of emotion over attention modulations in the lateral amygdala: intracranial ERP evidence from a patient with temporal lobe epilepsy. Cogn Affect Behav Neurosci. 2010;10(1):83–93. https://doi.org/10.3758/CABN.10.1.83.

    Article  PubMed  Google Scholar 

  42. Huijgen J, Dinkelacker V, Lachat F, Yahia-Cherif L, El Karoui I, Lemaréchal J-D, et al. Amygdala processing of social cues from faces: an intracrebral EEG study. Soc Cogn Affect Neurosci. 2015;10(11):1568–76. https://doi.org/10.1093/scan/nsv048.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Labudda K, Mertens M, Steinkroeger C, Bien CG, Woermann FG. Lesion side matters - an fMRI study on the association between neural correlates of watching dynamic fearful faces and their evaluation in patients with temporal lobe epilepsy. Epilepsy Behav. 2014;31:321–8. https://doi.org/10.1016/j.yebeh.2013.10.014.

    Article  PubMed  Google Scholar 

  44. Benuzzi F, Meletti S, Zamboni G, Calandra-Buonaura G, Serafini M, Lui F, et al. Impaired fear processing in right mesial temporal sclerosis: a fMRI study. Brain Res Bull. 2004;63(4):269–81. https://doi.org/10.1016/j.brainresbull.2004.03.005.

    Article  PubMed  Google Scholar 

  45. Schacher M, Haemmerle B, Woermann FG, Okujava M, Huber D, Grunwald T, et al. Amygdala fMRI lateralizes temporal lobe epilepsy. Neurology. 2006;66(1):81–7. https://doi.org/10.1212/01.wnl.0000191303.91188.00.

    Article  CAS  PubMed  Google Scholar 

  46. Steiger BK, Muller AM, Spirig E, Toller G, Jokeit H. Mesial temporal lobe epilepsy diminishes functional connectivity during emotion perception. Epilepsy Res. 2017;134:33–40. https://doi.org/10.1016/j.eplepsyres.2017.05.004.

    Article  PubMed  Google Scholar 

  47. Axmacher N, Henseler MM, Jensen O, Weinreich I, Elger CE, Fell J. Cross-frequency coupling supports multi-item working memory in the human hippocampus. Proc Natl Acad Sci USA. 2010;107(7):3228–33. https://doi.org/10.1073/pnas.0911531107.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Heusser AC, Poeppel D, Ezzyat Y, Davachi L. Episodic sequence memory is supported by a theta-gamma phase code. Nat Neurosci. 2016;19(10):1374–80. https://doi.org/10.1038/nn.4374.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Canolty RT, Knight RT. The functional role of cross-frequency coupling. Trends Cogn Sci. 2010;14(11):506–15. https://doi.org/10.1016/j.tics.2010.09.001.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Canolty RT, Edwards E, Dalal SS, Soltani M, Nagarajan SS, Kirsch HE, et al. High gamma power is phase-locked to theta oscillations in human neocortex. Science. 2006;313(5793):1626–8.

    Article  CAS  Google Scholar 

  51. Lisman JE, Jensen O. The theta-gamma neural code. Neuron. 2013;77(6):1002–16. https://doi.org/10.1016/j.neuron.2013.03.007.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Nandi B, Swiatek P, Kocsis B, Ding M. Inferring the direction of rhythmic neural transmission via inter-regional phase-amplitude coupling (ir-PAC). Sci Rep. 2019;9(1):6933. https://doi.org/10.1038/s41598-019-43272-w.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Nolte G, Ziehe A, Nikulin VV, Schlögl A, Krämer N, Brismar T, Müller KR. Robustly estimating the flow direction of information in complex physical systems. Phys Rev Lett. 2008;100(23): 234101.

    Article  Google Scholar 

  54. Frank DW, Costa VD, Averbeck BB, Sabatinelli D. Directional interconnectivity of the human amygdala, fusiform gyrus, and orbitofrontal cortex in emotional scene perception. J Neurophysiol. 2019;122:1530–7.

    Article  Google Scholar 

  55. Korioth D. Berühmte Personen nach Lebensalter. 2010. https://geboren.am/geburtstage. Accessed 30 Sep 2019.

  56. Bradley MM, Lang PJ. Measuring emotion: the self-assessment manikin and the semantic differential. J Behav Ther Exp Psychiatry. 1994;25(1):49–59.

    Article  CAS  Google Scholar 

  57. Tamraz J, Comair Y. Atlas of regional anatomy of the brain using MRI: with functional correlations. Berlin: Springer-Verlag; 2006.

    Google Scholar 

  58. Maris E, Oostenveld R. Nonparametric statistical testing of EEG- and MEG-data. J Neurosci Methods. 2007;164(1):177–90. https://doi.org/10.1016/j.jneumeth.2007.03.024.

    Article  PubMed  Google Scholar 

  59. Cohen MX. Analyzing neural time series data: theory and practice. 1st ed. Cambridge: The MIT Press; 2014.

    Book  Google Scholar 

  60. Bates D, Maechler M, Bolker B, Walker S. Fitting linear mixed-effects models using lme4. J Stat Softw. 2015;67(1):1–48. https://doi.org/10.18637/jss.v067.i01.

    Article  Google Scholar 

  61. Kuznetsova A, Brockhoff PB, Christensen RHB. lmerTest package: tests in linear mixed effects models. J Stat Softw. 2017;82(13):1–26. https://doi.org/10.18637/jss.v082.i13.

    Article  Google Scholar 

  62. Lenth RV. emmeans: Estimated Marginal Means, aka Least-Squares Means. 2020. https://CRAN.R-project.org/package=emmeans. Accessed 11 May 2021.

  63. Jiang H, Bahramisharif A, Van Gerven MA, Jensen O. Measuring directionality between neuronal oscillations of different frequencies. Neuroimage. 2015;118:359–67. https://doi.org/10.1016/j.neuroimage.2015.05.044.

    Article  PubMed  Google Scholar 

  64. Koivisto M, Lähteenmäki M, Sørensen TA, Vangkilde S, Overgaard M, Revonsuo A. The earliest electrophysiological correlate of visual awareness? Brain Cogn. 2008;66(1):91–103. https://doi.org/10.1016/j.bandc.2007.05.010.

    Article  PubMed  Google Scholar 

  65. Wilenius-Emet M, Revonsuo A, Ojanen V. An electrophysiological correlate of human visual awareness. Neurosci Lett. 2004;354(1):38–41. https://doi.org/10.1016/j.neulet.2003.09.060.

    Article  CAS  PubMed  Google Scholar 

  66. Taubert J, Flessert M, Wardle SG, Basile BM, Murphy AP, Murray EA, Ungerleider LG. Amygdala lesions eliminate viewing preferences for faces in rhesus monkeys. Proc Natl Acad Sci USA. 2018;115(31):8043–8. https://doi.org/10.1073/pnas.1807245115.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  67. Adolphs R, Gosselin F, Buchanan TW, Tranel D, Schyns P, Damasio AR. A mechanism for impaired fear recognition after amygdala damage. Nature. 2005;433(7021):68–72.

    Article  CAS  Google Scholar 

  68. Gamer M, Schmitz AK, Tittgemeyer M, Schilbach L. The human amygdala drives reflexive orienting towards facial features. Curr Biol. 2013;23(20):R917–8. https://doi.org/10.1016/j.cub.2013.09.008.

    Article  CAS  PubMed  Google Scholar 

  69. Meletti S, Benuzzi F, Rubboli G, Cantalupo G, Stanzani Maserati M, Nichelli P, Tassinari CA. Impaired facial emotion recognition in early-onset right mesial temporal lobe epilepsy. Neurology. 2003;60(3):426–31. https://doi.org/10.1212/wnl.60.3.426.

    Article  CAS  PubMed  Google Scholar 

  70. Seidenberg M, Griffith R, Sabsevitz D, Moran M, Haltiner A, Bell B, et al. Recognition and identification of famous faces in patients with unilateral temporal lobe epilepsy. Neuropsychologia. 2002;40(4):446–56. https://doi.org/10.1016/s0028-3932(01)00096-3.

    Article  PubMed  Google Scholar 

  71. Volfart A, Jonas J, Maillard L, Busigny T, Rossion B, Brissart H. Typical visual unfamiliar face individuation in left and right mesial temporal epilepsy. Neuropsychologia. 2020;147: 107583. https://doi.org/10.1016/j.neuropsychologia.2020.107583.

    Article  PubMed  Google Scholar 

  72. Gomez-Ibañez A, Urrestarazu E, Viteri C. Recognition of facial emotions and identity in patients with mesial temporal lobe and idiopathic generalized epilepsy: an eye-tracking study. Seizure. 2014;23(10):892–8. https://doi.org/10.1016/j.seizure.2014.08.012.

    Article  PubMed  Google Scholar 

  73. Halgren E, Baudena P, Heit G, Clarke JM, Marinkovic K, Clarke M. Spatio-temporal stages in face and word processing. I. Depth-recorded potentials in the human occipital, temporal and parietal lobes. J Physiol Paris. 1994;88(1):1–50. https://doi.org/10.1016/0928-4257(94)90092-2.

    Article  CAS  PubMed  Google Scholar 

  74. Gao C, Conte S, Richards JE, Xie W, Hanayik T. The neural sources of N170: understanding timing of activation in face-selective areas. Psychophysiology. 2019;56(6): e13336. https://doi.org/10.1111/psyp.13336.

    Article  PubMed  PubMed Central  Google Scholar 

  75. Rossion B. Understanding face perception by means of human electrophysiology. Trends Cogn Sci. 2014;18(6):310–8. https://doi.org/10.1016/j.tics.2014.02.013.

    Article  PubMed  Google Scholar 

  76. Bentin S, Allison T, Puce A, Perez E, McCarthy G. Electrophysiological studies of face perception in humans. J Cogn Neurosci. 1996;8(6):551–65.

    Article  Google Scholar 

  77. Ghuman AS, Brunet NM, Li Y, Konecky RO, Pyles JA, Walls SA, et al. Dynamic encoding of face information in the human fusiform gyrus. Nat Commun. 2014;5:5672. https://doi.org/10.1038/ncomms6672.

    Article  CAS  PubMed  Google Scholar 

  78. Engell AD, McCarthy G. Face, eye, and body selective responses in fusiform gyrus and adjacent cortex: an intracranial EEG study. Front Hum Neurosci. 2014;8:642. https://doi.org/10.3389/fnhum.2014.00642.

    Article  PubMed  PubMed Central  Google Scholar 

  79. Bar M, Aminoff E. Cortical analysis of visual context. Neuron. 2003;38(2):347–58. https://doi.org/10.1016/s0896-6273(03)00167-3.

    Article  CAS  PubMed  Google Scholar 

  80. Kveraga K, Ghuman AS, Kassam KS, Aminoff EA, Hämäläinen MS, Chaumon M, Bar M. Early onset of neural synchronization in the contextual associations network. Proc Natl Acad Sci USA. 2011;108(8):3389–94. https://doi.org/10.1073/pnas.1013760108.

    Article  PubMed  PubMed Central  Google Scholar 

  81. Mueller SG, Laxer KD, Cashdollar N, Lopez RC, Weiner MW. Spectroscopic evidence of hippocampal abnormalities in neocortical epilepsy. Eur J Neurol. 2006;13(3):256–60.

    Article  CAS  Google Scholar 

  82. Scott RC, Cross JH, Gadian DG, Jackson GD, Neville BG, Connelly A. Abnormalities in hippocampi remote from the seizure focus: a T2 relaxometry study. Brain. 2003;126(Pt 9):1968–74.

    Article  Google Scholar 

  83. Bernasconi N, Bernasconi A, Caramanos Z, Antel SB, Andermann F, Arnold DL. Mesial temporal damage in temporal lobe epilepsy: a volumetric MRI study of the hippocampus, amygdala and parahippocampal region. Brain. 2003;126(Pt 2):462–9.

    Article  CAS  Google Scholar 

  84. Yu O, Mauss Y, Namer IJ, Chambron J. Existence of contralateral abnormalities revealed by texture analysis in unilateral intractable hippocampal epilepsy. Magn Reson Imaging. 2001;19(10):1305–10.

    Article  CAS  Google Scholar 

  85. Gloor P, Salanova V, Olivier A, Quesney LF. The human dorsal hippocampal commissure. An anatomically identifiable and functional pathway. Brain. 1993;116(Pt 5):1249–73.

    Article  Google Scholar 

  86. Rosenzweig I, Beniczky S, Brunnhuber F, Alarcon G, Valentin A. The dorsal hippocampal commissure: when functionality matters. J Neuropsychiatry Clin Neurosci. 2011;23(3):E45–8. https://doi.org/10.1176/appi.neuropsych.23.3.E45.

    Article  PubMed  Google Scholar 

  87. Lacruz ME, García Seoane JJ, Valentin A, Selway R, Alarcón G. Frontal and temporal functional connections of the living human brain. Eur J Neurosci. 2007;26(5):1357–70.

    Article  CAS  Google Scholar 

  88. Wilson CL, Isokawa M, Babb TL, Crandall PH, Levesque MF, Engel J, JR,. Functional connections in the human temporal lobe. II. Evidence for a loss of functional linkage between contralateral limbic structures. Exp Brain Res. 1991;85(1):174–87.

    Article  CAS  Google Scholar 

  89. Bastos AM, Usrey WM, Adams RA, Mangun GR, Fries P, Friston KJ. Canonical microcircuits for predictive coding. Neuron. 2012;76(4):695–711. https://doi.org/10.1016/j.neuron.2012.10.038.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We are grateful to all of the study participants for their patience and cooperation.

Funding

Open Access funding enabled and organized by Projekt DEAL. This study received no specific grant from any funding agency in the public, commercial, or not-for-profit sectors.

Author information

Authors and Affiliations

Authors

Contributions

JWK conceptualized the study, analyzed and interpreted the data, drafted and revised the manuscript. KEB, CB, WH and EKE played a major role in the acquisition of data and revised the manuscript for intellectual content. MLVQ analyzed and interpreted the data. MS and SRGS played a major role in the acquisition of data. JMH and ML supervised the study, analyzed and interpreted the data and revised the manuscript for intellectual content; both authors contributed equally to this work. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Ji-Won Kim.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the institutional review board of the Protestant Hospital Alsterdorf, Hamburg, Germany. The study was performed in accordance with the ethical standards as laid down in the 1964 Declaration of Helsinki and its later amendments or comparable ethical standards. Informed consent was obtained from all individual participants included in the study.

Consent for publication

Not applicable.

Competing interests

SRGS reports personal fees from UCB Pharma, GW Pharmaceuticals, Eisai Pharma and Desitin Pharma outside the submitted work. There are no competing interests to report for JWK, KEB, CB, WH, MS, MLVQ, EKE, JMH and ML.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Kim, JW., Brückner, K.E., Badenius, C. et al. Face-induced gamma oscillations and event-related potentials in patients with epilepsy: an intracranial EEG study. BMC Neurosci 23, 36 (2022). https://doi.org/10.1186/s12868-022-00715-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12868-022-00715-x

Keywords