- 1Epilepsy Center, Medical Center–University of Freiburg, Faculty of Medicine, University of Freiburg, Freiburg im Breisgau, Germany
- 2Spemann Graduate School of Biology and Medicine (SGBM), University of Freiburg, Freiburg im Breisgau, Germany
- 3Faculty of Biology, University of Freiburg, Freiburg im Breisgau, Germany
- 4Department of Neurology, Paracelsus Medical University Salzburg, Salzburg, Austria
Human High-Frequency-Oscillations (HFO) in the ripple band are oscillatory brain activity in the frequency range between 80 and 250 Hz. HFOs may comprise different subgroups that either play a role in physiologic or pathologic brain functions. An exact differentiation between physiologic and pathologic HFOs would help elucidate their relevance for cognitive and epileptogenic brain mechanisms, but the criteria for differentiating between physiologic and pathologic HFOs remain controversial. In particular, the separation of pathologic HFOs from physiologic HFOs could improve the identification of epileptogenic brain regions during the pre-surgical evaluation of epilepsy patients. In this study, we performed intracranial electroencephalography recordings from the hippocampus of epilepsy patients before, during, and after the patients completed a spatial navigation task. We isolated hippocampal ripples from the recordings and categorized the ripples into the putative pathologic group iesRipples, when they coincided with interictal spikes, and the putative physiologic group isolRipples, when they did not coincide with interictal spikes. We found that the occurrence of isolRipples significantly decreased during the task as compared to periods before and after the task. The rate of iesRipples was not modulated by the task. In patients who completed the spatial navigation task on two consecutive days, we furthermore examined the occurrence of ripples in the intervening night. We found that the rate of ripples that coincided with sleep spindles and were therefore putatively physiologic correlated with the performance improvement on the spatial navigation task, whereas the rate of all ripples did not show this relationship. Together, our results suggest that the differentiation of HFOs into putative physiologic and pathologic subgroups may help identify their role for spatial memory and memory consolidation processes. Conversely, excluding putative physiologic HFOs from putative pathologic HFOs may improve the HFO-based identification of epileptogenic brain regions in future studies.
Introduction
High Frequency Oscillations (HFO) are an electrographic manifestation of hyper-synchronized neurons and are subdivided into Ripples and Fast-Ripples according to the frequency range of the oscillations (80–250 and 250–500 Hz, respectively) (1, 2). In the field of epilepsy, Ripples and Fast-Ripples were initially considered improved biomarkers of epileptogenic networks (3–8). However, recent research has drawn a more complex picture (9–11) and has highlighted the importance of being able to differentiate between physiologic and pathologic HFOs.
Interictal Epileptic Spikes (IES) are another type of epileptic biomarker. IES are common in patients with epilepsy, have a waveform of a fast transient, are commonly generated in epileptic cortex and reflect a hyper-excitability of neural networks (12). IES are very sensitive but not specific to epileptogenic areas (13, 14). Ripples are known to coincide temporally and spatially with IES to some extent. These IES coincident Ripples appear to have different amplitude and waveform characteristics when compared to Ripples associated with physiologic events such as sleep spindles (15). IES coincident Ripples may be more sensitive to the seizure-onset-zone than Fast-Ripples and also more specific to it than Ripples occurring in isolation from IES (16). IES coincident Ripples may better predict post-surgical outcomes than Ripples not coinciding with IES (17) and Ripples coinciding with IES showed the highest correspondence with the resected volume in seizure-free patients as compared to other HFO subgroups (18). Moreover, a combined marker composed of IES and HFO occurrence rates appeared to be useful for estimating the epileptogenic zone (11). Together, the coincidence with IES may constitute a good criterion for separating pathologic Ripples from physiologic Ripples.
Sleep spindles are a third type of electrographic pattern which is observed in the human electroencephalogram (EEG) recorded with scalp or implanted electrodes (19). Sleep-spindle events have a distinct oscillatory waveform with durations between 0.5 and 3.0 s and frequencies between 11 and 16 Hz (20–22). Sleep spindles are generated and controlled by thalamic networks, with several hypotheses linking them to gating functions of sensory information flow. However, so far, the complete and definitive functional meaning of sleep spindles remains to be explored (23). Amongst others, sleep spindles may be relevant for memory consolidation during sleep, particularly when coupled to hippocampal Ripples (24–32). The occurrence rate of sleep spindle-coupled Ripples during sleep may thus reflect the intensity of memory consolidation.
Based on this prior knowledge, we hypothesized in the current study that hippocampal ripples could be differentiated into a putatively pathologic subgroup (iesRipples) and a putatively physiologic subgroup (isolRipples) based on their temporal and spatial coincidence with IES: iesRipples should coincide with an IES temporally (i.e., occurring within the duration of an IES) and spatially (i.e., when recorded on the same channel as an IES), whereas isolRipples should occur in isolation. Moreover, we hypothesized that ripples temporally co-occurring with ipsilateral hippocampal sleep spindles (spindleRipples) could serve as a marker for memory-consolidation processes. In our analyses, we therefore assessed whether the occurrence rates of iesRipples and of isolRipples were altered during the spatial navigation task. The spatial navigation task required the patients to form associative memories between objects and their corresponding locations and thus imposed an increased cognitive demand on the patients. Since physiological ripples are associated with cognitive functioning, we hypothesized that the rate of isolRipples should be altered during the task, whereas the activity of iesRipples should be unaffected. Additionally, in patients who performed the spatial navigation task on two consecutive days, we assessed the correlation between the occurrence rate of spindleRipples in the intervening night and the performance improvement of the patients between both days. Based on the proposed role of sleep spindle-coupled ripples in memory consolidation, we hypothesized that a higher rate of spindleRipples should be associated with a greater performance improvement.
Methods
Patient Selection
Participating patients (Table 1) suffered from pharmaco-resistant focal epilepsy and underwent intracranial EEG (iEEG) recordings from the hippocampus to identify their seizure onset zone at the Freiburg Epilepsy Center, Germany. The clinical decision for the implantation of iEEG electrodes was made individually for each patient in cases when the epileptogenic zone remained unclear using non-invasive methods. All patients gave their informed consent to participate in a study aiming at the identification of electrophysiological correlates of cognitive processing, including a spatial navigation task. A total of 19 patients performing the spatial navigation task were included in the current study, 9 of which completed the spatial navigation task on two subsequent days. All patients gave their written informed consent and the study was approved by the Ethics Committee of the Freiburg University Medical Center.
Spatial Navigation Paradigm
The spatial navigation paradigm was adapted from previous studies (33–35). In this paradigm, the patients performed an object-location memory task navigating freely in a circular virtual environment. The environment comprised a grassy plain bounded by a cylindrical cliff. Two mountains, a sun, and several clouds provided patients with distal orientation cues. Patients completed the task on a laptop using the arrow keys for moving forward and turning left and right and the spacebar to indicate their response. Patients were asked to complete up to 160 trials but were instructed to pause or quit the task whenever they wanted. At the very beginning, patients collected eight everyday objects (randomly drawn from a total number of 12 potential objects) from different locations in the arena. Objects appeared one after the other. Afterward, patients completed variable numbers of trials, depending on compliance. Each trial consisted of four different phases (Figure 1). First, one of the eight objects was presented for 2 s (cue presentation). Afterwards, patients were asked to navigate to the associated object location within the virtual environment (retrieval). After patients had indicated their response via a button press at the assumed object location, they received feedback depending on response accuracy. Response accuracy was measured as the distance between the remembered location and the correct location (drop error). Last, the object was presented in the correct location, and patients had to collect the object from there to further improve their associative memory between the object and its location. After each trial, a fixation crosshair was shown for a variable duration of 3–5 s. Triggers were detected using a phototransistor attached to the screen marking onsets and offsets of the cue presentation phase. The intra-paradigm period of the iEEG was then delimited by the first and last phototransistor triggers. We calculated the patients' performance in each trial as the ratio between the drop error and the largest possible drop error (maximum random drop error):
The largest possible drop error in a given trial was determined by creating one million random locations within the virtual environment and then selecting the location with the largest Euclidean distance to the correct object location. The performance for the entire paradigm was calculated as the median performance across all trials.
Figure 1. Associative object-location memory task during virtual spatial navigation. At the beginning of the experiment, patients collected eight different objects from different locations within the virtual environment. Afterward, patients completed variable numbers of retrieval trials, during which they were first presented with one of the eight objects serving as cue (cue presentation). Patients then navigated to the remembered location of that object (retrieval) and made a response. Following this response, patients received feedback via an emoticon (feedback) and had to collect the object from its correct location (re-encoding).
Identification of Hippocampal and White Matter Channels
Preimplantation and post-implantation MRI scans were available for all patients. Electrode localization was performed using FSL (https://fsl.fmrib.ox.ac.uk/fsl/fslwiki/FSL) and PyLocator (http://pylocator.thorstenkranz.de/). The post-implantation MR image was coregistered with the preimplantation MR image. Next, the preimplantation MR image was skull-stripped and normalized to MNI space, and the same normalization matrix was applied to the post-implantation MR image. Normalized post-implantation images were visually inspected using PyLocator, and channel locations were manually identified. For the analyses in this study, we only considered electrode channels located in the right or left hippocampus.
Detection of Ripples, isolRipples, iesRipples, and IES
Intracranial EEG data were recorded using “Profusion EEG Software” (Compumedics Limited, Abbotsford Victoria, Australia). Original signals were low-pass filtered at 800 Hz and sampled at 2 kHz using Cz as a hardware reference. All signal analyses were performed using bipolar montages.
The detection of Ripples and IES on the hippocampal iEEG signals was accomplished by using automatic detectors (36), which were based on multivariate classifications of iEEG epochs using kernelized support-vector-machines. The spectral analysis of these detections was performed using the Morlet wavelet transform, the maximum and minimum frequencies clustering the events' power corresponded to those frequency-bins with a power contribution within that one of the spectral peak (i.e., the frequency with the maximum power contribution) +/– the standard deviation of the power across all frequency bins.
The automatic detectors were run on time-selected segments of the hippocampal iEEG signals. These segments corresponded to the pre-, intra- and post-paradigm phases. The intra-paradigm phase was delimited by the first and last cue-presentation triggers. The duration of the pre- and post-paradigm phases was the same as for the intra-paradigm phase. The pre-paradigm phase ended 30 min before the start of the intra-paradigm phase and the post-paradigm phase started 30 min after the end of the intra-paradigm phase.
The automatic detectors provided discrete events of the classes Ripples and IES. Each event comprised a start time and an end time. We used custom scripts in MATLAB 2018b that determined the temporal and spatial coincidence of the Ripples and IES in order to identify ripples belonging to the class iesRipples. Ripples which were not coincident with IES formed the class isolRipples. For each class of ripples, we then calculated their occurrence rate per minute in each of the hippocampal channels. If a patient had more than one hippocampal electrode, we averaged the occurrence rates across the different hippocampal channels (Figure 2).
Figure 2. Determining the hippocampal activity for Patient 16 and each event class. Firstly, the automatic detections of Ripples, iesRipples, isolRipples, and IES were used to characterize each channel (y-axes) with their occurrence rate per minute (heat-maps). The specific brain region Hippocampus (red channel labels) was then characterized by the average occurrence rate across all hippocampal channels. The average hippocampal activity for each event class and the pre-, intra- and post-paradigm segments is shown in parentheses on the x-axis.
Analysis of the Ripples, isolRipples, iesRipples, and IES Activity
We tested for the effect of the spatial navigation paradigm on the occurrence rate of Ripples, isolRipples, iesRipples, and IES using the data from all 19 patients. In patients who performed the task on two consecutive days, we only used the data from the first day for this analysis in order to avoid that these patients had a stronger effect on the statistical results. A mixed ANOVA was conducted to test for main effects and interactions between the factors “time period” (pre-, intra- and post-paradigm phase) and “ripple class” (isolRipples vs. iesRipples) on the occurrence rate of hippocampal ripples. In a post-hoc analysis, we firstly used a two-tailed, non-parametric Wilcoxon signed rank test to analyze whether the occurrence rates differed (i) between the pre-paradigm and the intra-paradigm phase, (ii) between the intra-paradigm phase and the post-paradigm phase, and (iii) between the control phases pre-paradigm and post-paradigm; at this stage no correction for multiple comparisons was applied since the families (i.e., pre vs. intra, intra vs. post, pre vs. post) had no repeated analyses (i.e., α = 0.05) (37). Finally, as part of the same post-hoc analysis, we performed either a subsequent left-tailed or a right-tailed Wilcoxon signed rank test to analyze if the activity-difference between phases corresponded to an increase or a decrease; at this stage a Bonferroni correction was applied to the significance threshold (i.e., α = 0.025) since two null hypotheses were considered (two-tailed and one-tailed Wilcoxon signed rank tests).
Detection of spindleRipples
The spindleRipple events corresponded to those Ripples which were temporally and spatially coincident with hippocampal sleep spindles. The detection of hippocampal sleep spindles was achieved using an automatic detector (38) based on multivariate classifications of iEEG epochs using kernelized support-vector-machines. We used a custom MATLAB 2018b script to determine coinciding ripples and sleep spindles, which then composed the class spindleRipples. Only ripples that were completely within the start and end time of a sleep spindle were considered to be temporally coincident. The spatial coincidence of a ripple and a sleep spindle was present if they were both hippocampal and ipsilateral. All ripples complying with these rules of temporal and spatial coincidence comprised the spindleRipples class.
We estimated the occurrence rate of spindleRipples in patients who performed the spatial navigation task on two consecutive days. For all these patients, we used the data between 10:00 pm of the first day and 6:00 am of the second day to estimate the occurrence rate of spindleRipples. This time period was selected with the aim of maximizing the inclusion of non-rapid eye movement N2 sleep stages, since the occurrence and power of sleep spindles is highest during this sleep stage (21, 23, 39, 40). If a patient had more than one hippocampal electrode channel, we averaged the channel-specific spindleRipple rates across the different channels to obtain one overall occurrence rate, which quantified the number of spindleRipples per minute.
Analysis of spindleRipples and Their Correlation With Spatial Navigation Performance
This analysis was only performed with the data from the patients who performed the spatial navigation task on two consecutive days. To test whether the occurrence rates of spindleRipples could be a marker of memory consolidation of the associative object-location memories that the patients formed during the spatial navigation task, we calculated the correlation between spindleRipple rates and the difference in performance obtained on days 1 and 2 (performance Δ). A non-paired, non-parametric, left-tailed Mann–Whitney U test was applied to each patient using all their trials from days 1 and 2 to test if the performance Δ was significant. The correlation between the hippocampal spindleRipple rates and the performance Δ was measured using Spearman's rank correlation coefficient. As a control, we quantified this relationship while controlling for the potential effect of the number of trials on day 1 using a partial correlation.
Results
Navigation Paradigm
A total of 19 patients completed the spatial navigation task (day 1). A subgroup of nine patients completed the spatial navigation task also on the subsequent day (day 2). The performance from each patient was measured as the median performance across all trials. For day 1, the maximum, minimum and mean performance corresponded to 86% (patient 9), 48% (patient 7) and 67 ± 11% (mean ± SD), respectively (Figure 3A).
Figure 3. Performance in the spatial navigation paradigm. (A) Median performance across trials obtained by the 19 patients on their first (or only) day of conducting the spatial-navigation task. (B) The line plots depict the median performance across trials for each of the nine patients that conducted the spatial navigation paradigm on 2 consecutive days (day 1 and 2). The bar plots depict the average across patients.
In the subgroup of the 9 patients conducting the paradigm on two consecutive days, the maximum, minimum and mean performance on day 1 corresponded to 86% (patient 9), 48% (patient 7), and 72 ± 12%, respectively. On day 2, the maximum, minimum, and mean performance corresponded to 93% (patient 9), 53% (patient 7), and 79 ± 13%, respectively (Figure 3B). All patients showed a significant performance improvement, when comparing day 1 with day 2 (left-tailed Mann–Whitney U, all p < 0.025).
Detection of isolRipples and iesRipples
The automatic detection on the hippocampal channels generated an average of 1,762 Ripple events and 1,042 IES events per patient across the pre-, intra- and post-paradigm phases. Of these Ripples, 74% occurred isolated from IES and were thus labeled isolRipples; the remaining 26% of the Ripples occurred spatially and temporally coinciding with an IES and were thus labeled iesRipples. A depiction of example waveforms and the corresponding spectrograms of isolRipples and iesRipples is provided in Figure 4.
Figure 4. Detected hippocampal ripples and spikes. (A) Selected events from each event type: the- first row shows the waveform from the raw EEG, the second row shows the waveform from the bandpassed EEG (80–250 Hz), the third row shows the spectrogram from each event-type obtained using a Morlet wavelet transform. (B) Average event-types from patient 19. (C) Grand average of the three event-types from all patients.
A challenge when detecting HFOs is the correct handling of fast transients which, when filtered, can produce artifacts resembling authentic HFOs (41, 42). A depiction of the waveform and spectrogram of the IES detections is provided in Figure 4 in order to show the handling of fast transients by the automatic detectors and to allow the comparison with the spectrogram from the detected Ripples. The isolRipples showed a distinct increase in power within a narrow band which resembled a blob in the spectrogram, the tuples (Hz) consisting of the spectral-peak, blob-lower-frequency and blob-higher-frequency for the single-event, averaged-events and patients-average were (117, 82, 144), (109, 82, 144), and (109, 88, 144), respectively. This blob-like power increase was also shown by iesRipples when considering only patient 19. The same tuples for iesRipples were (82, 82, 102), (102, 82, 144), and (109, 82, 165), respectively.
The spectrogram of the IES depicted in the third column of Figure 4 showed typical spectral characteristics of a fast-transient, which much like a single pulse are represented in all frequency bins when analyzed by a time-frequency transformation. The key characteristic that differentiates a filtering artifact from a real HFO event is then this blob-like, narrow-band power increase which is shown by the spectrograms of both isolRipples and iesRipples.
The events forming the averages in Figure 4 were centered using their maximum peak. Hence a peak was formed at the center of the raw averages while surrounding samples were canceled because of varying pre- and post-event waveforms across events.
Modulation of Ripple Activity by Increased Cognitive Demands During the Spatial Navigation Task
The mixed ANOVA test showed a non-significant main effect of time period (pre-, intra-, post-paradigm) (F = 0.679, p = 0.510), a significant main effect of the ripple class (isolRipples, iesRipples) (F = 8.948, p < 0.005), and a significant interaction between both factors (F = 4.069, p < 0.025). To further understand this interaction, we performed post-hoc analyses (Wilcoxon signed rank tests).
We found no significant differences in the occurrence rates of allRipples for the pair-wise comparisons (pre-paradigm vs. intra-paradigm, p < 0.717; intra-paradigm vs. post-paradigm, p < 0.243; pre-paradigm vs. post-paradigm, p < 0.314). Hence, the increased cognitive load exerted by the spatial-navigation paradigm did not modulate the activity from the allRipples event class.
For isolRipples, the comparison pre-paradigm vs. intra-paradigm showed a significant difference (p < 0.043), the subsequent right-tailed test showed a significance of p < 0.022, hence the isolRipples presented an activity decrease when transitioning from the pre- to the intra-paradigm phase. Similarly, the comparison of isolRipples from the intra-paradigm vs. post-paradigm phase showed a significant difference (p < 0.007), the subsequent left-tailed test showed a significance of p < 0.004. Hence, the isolRipples presented an activity increase when transitioning from the intra- to the post-paradigm phase. When comparing the occurrence rates of isolRipples during the control phases pre-paradigm and post-paradigm, no significant difference was found (p > 0.314). In summary, the increased cognitive load exerted by the spatial navigation paradigm did modulate the activity from the isolRipples in a way that their activity was reduced during the paradigm (Figure 5), additionally, the isolRipples activity returned to the pre-paradigm control values after the phase of increased cognitive load.
Figure 5. Modulation of the event-classes by the increased cognitive load. The bar plots show pairs of activity differences between phases (Activity Δ), i.e., difference between the intra- and pre-paradigm phase (orange) and difference between the intra- and post-paradigm phase (blue), for each patient and using each of the event-classes. The bar plots with a negative value indicate a decrease in activity when entering or leaving the period of increased cognitive load (i.e., intra-paradigm phase). The bar plots with a positive value indicate an increase in activity during the paradigm when compared to the pre- or post-paradigm phases. The hippocampal isolRipples were the only event-class to show a significant modulation by the navigation task, presenting an activity decrease when comparing the intra-paradigm occurrence rates with those from both the preceding (pre-paradigm, p < 0.043) and succeeding (post-paradigm, p < 0.007) control phases.
We found no significant differences in the occurrence rates of iesRipples for the pair-wise comparisons (pre-paradigm vs. intra-paradigm, p > 0.277; intra-paradigm vs. post-paradigm, p > 0.295; pre-paradigm vs. post-paradigm, p > 0.260). Thus, the increased cognitive load exerted by the spatial navigation paradigm did not modulate the activity from the iesRipples event class.
For IES, we found no significant differences for the pair-wise comparisons pre-paradigm vs. intra-paradigm (p > 0.778) and intra-paradigm vs. post-paradigm phase (p > 0.260); however, the IES activity from the control phase pre-paradigm was significantly different from the post-paradigm activity (p < 0.022), the subsequent left-tailed test showed a significance of p < 0.012. These results suggest that the IES activity was not modulated by the increased cognitive load, however the IES activity from the post-paradigm control phase was higher than the IES activity during the pre-paradigm control phase.
The results described in this section are also shown in Table 2 and summarized in Figure 5.
Improvement in Performance and the Correlation With spindleRipples
An example of a detected sleep spindle is presented in Figure 6. The performances obtained by the patients when conducting the spatial navigation task on days 1 and 2 are presented in Table 3.
As previously mentioned, all the patients' performance for the navigation task improved on day 2 when comparing it with day 1 (left-tailed Mann–Whitney U all p < 0.025). During the night between the two sessions of the spatial navigation task, hippocampal Ripples were detected and their occurrence rate (incidences per minute) was obtained. We found that the activity from the allRipples event class showed no correlation with the patients' performance improvement when repeating the spatial navigation task (rho = 0.13, p = 0.74; controlling for the number trials on day 1: rho = 0.05, p = 0.90). In contrast, the occurrence rate of the detected spindleRipples showed a significant positive correlation with performance improvement (rho = 0.73, p = 0.03; controlling for the number of trials on day 1, rho = 0.77, p < 0.03) (Figure 7).
Figure 7. Occurrence Rates of the event-class all hippocampal ripples and spindleRipples and their correlation with the patients' difference in performance when comparing day 1 and 2 (i.e., performance Δ).
As a control, we also computed correlations between the occurrence rate of ripples and the performance values on both days in order to present evidence that the association was presumably due to the learning process and not simply related to the patients' general performance level. These correlations were not significant (for performance on day 1 and allRipples rate, rho = 0.17, p = 0.68; performance on day 1 and spindleRipples rate, rho = −0.14, p = 0.74; performance on day 2 and allRipples rate, rho = 0.32, p = 0.41; performance on day 2 and spindleRipples rate, rho = 0.05, p = 0.91).
Discussion
In this study, we examined the activity of putatively physiologic and pathologic ripples in the human hippocampus during a spatial navigation task. We differentiated between the two groups of ripples by analyzing their coincidence with interictal spikes. We found that the rate of the putatively physiologic ripples, isolRipples, decreased during the task as compared to pre- and post-task phases, whereas the putatively pathologic ripples, iesRipples, did not show this effect. In addition, the ripples associated with hippocampal sleep spindles, spindleRipples, showed a positive correlation with the performance improvement of patients who completed the spatial navigation task on two consecutive days.
In summary, the increased cognitive demands from the spatial navigation task exerted differential effects on iesRipples and isolRipples; these results [in addition to previously published evidence showing iesRipples as the most and isolRipples as the least accurate estimator of the epileptogenic-zone (18)] hence provide further evidence to support the putatively physiologic and pathologic nature initially attributed to isolRipples and iesRipples. Furthermore, our study suggests that ripples associated with sleep spindles may constitute a marker of memory-consolidation processes.
Our analysis was focused on the hippocampus, as it plays a major role for both spatial navigation (24–30) and declarative memory formation (43–50) and is a candidate region for generating epileptic activity (51, 52) and is therefore often assessed for the decision making and planning for surgical epilepsy therapy.
Our first goal was to analyze the activity of putative physiologic and putative pathologic Ripples during periods when the cognitive load was higher than normal and then compare this activity with controls coming from periods with a lower demand on cognitive functioning. During these periods with a lower cognitive demand, all patients stayed in their hospital bed in the fairly quiet and stable environment of their room. It can thus be assumed that cognitive demands during these periods were indeed lower than during the spatial navigation task. Furthermore, all hippocampal activities reported corresponded to the average across the duration of the pre-, intra-, and post-paradigm phases; short increases of cognitive-load during the control phases would have been averaged out. In contrast, during the intra-paradigm phase it is conceivable that the reported average was derived from a period of a constant and increased cognitive-load.
The periods selected as controls before and after the paradigm showed no difference in the activity from either allRipples, isolRipples, or iesRipples, which in the case of the allRipples and iesRipples, as previously mentioned, evidences a lack of modulation exerted by the period of increased cognitive-load. In the case of isolRipples however, this equivalence between controls also shows that the isolRipples activity was modulated by the increased cognitive-load and then returned to the levels previous to the conducted navigation task.
Differentiation of Pathologic and Physiologic Ripples
Numerous publications have reported, for both the human and non-human brain, on the occurrence of Ripples during wakefulness and while conducting cognitive tasks (43–50); despite this, to our knowledge, only the studies from (47) and those from the Brázdil1 group (48–50), have analyzed the effects of cognitive processes on the activity of putatively physiologic and putatively pathologic Ripples during the awake state and in humans. The study from (50) presents the largest patient cohort with 36 patients and will thus be considered for further discussion. This study explored if the effect of cognitive load in the form of different tasks (visual oddball, Go/NoGo, Ultimatum Game, Mismatch Negativity) on putative-pathologic ripples (i.e., ripples from epileptic hippocampi, hereinafter referenced as pathoBrazRipples) differed from the effect on putative-physiologic ripples (i.e., from non-epileptic hippocampi, hereinafter referred to as physioBrazRipples). Both pathoBrazRipples and physioBrazRipples were reported to show a significant activity reduction when transitioning from the pre- to the intra-paradigm phase, however this reduction was more significant for the physioBrazRipples than for the pathoBrazRipples when averaging the activity across the analyzed hippocampal channels.
Our results agree with (50) in that both isolRipples and physioBrazRipples showed a decrease in their activity during the intra-paradigm phase when compared to either the resting states pre-paradigm or post-paradigm. Our results disagree with (50) in that in contrast to the decrease in activity of pathoBrazRipples during the task, our iesRipples did not show any modulation exerted by the increased cognitive load. The differences between our results and those from (50) are likely due to the fact that neither approach is exhaustive, i.e., both approaches are likely to increase the proportion of physiologic to pathologic Ripples but it is still possible that this formed sub-groups of Ripples are not exclusively physiologic or pathologic.
Another recent study, (47), examined ripple-occurrence rates across two cognitive tasks and a resting state during wakefulness. This study detected ripples firstly in the time domain by thresholding the power in the ripple band, and secondly by only accepting those detections with spectral power bursts narrowed down to the ripple range. This procedure used for the removal of potential artifactual ripples will produce detected events with spectral characteristics resembling those from our isolRipples class (Figure 4). Interestingly, in agreement with our findings for isolRipples, the results from (47) showed their ripple occurrence-rate to be higher during the resting state than during the cognitive tasks.
Modulation of IES Activity
Interactions between epileptic activity and cognition have been discussed for many years (53). We compare our results with those obtained by other studies analyzing IES activity changes during cognitive tasks. The studies from (54–56) give evidence that cognitive tasks and movements can change the properties of epileptogenic networks and thus the occurrence of IES, these studies however provide disagreeing conclusions on the activity patterns followed by the reported IES. The work from (54) reported a reduction of the spike rate during successful encoding while conducting a visual recognition memory task in amygdala, hippocampus, and temporal cortex. In agreement with the latter study, (55) showed a decrease of IES activity during movement in two patients with a focal cortical dysplasia in the pre- and/or post-central gyrus. The more recent study by (56) presented an increase of temporal lobe interictal spikes in the hippocampus during a spatial memory task and both in hippocampus and lateral temporal lobe during an episodic memory task.
In contrast with the mentioned conflicting studies, our results did not show any significant modulation, whether increasing or decreasing, of the IES activity during the period of increased cognitive load. We did find, however, a difference between the control phases pre-paradigm and post-paradigm. The difference in activity from our control periods then calls for a further exploration of the importance of selecting a control period, which can then allow the comparison of results between studies.
Correlation of spindleRipples and Memory Consolidation
The fact that all patients improved their performance when repeating the spatial-navigation paradigm provides evidence that the used paradigm did in fact exert a cognitive load which lead to the learning of newly acquired information. Interestingly, we found a strong correlation between spindleRipple activity and performance improvement.
To our knowledge only one other study has analyzed Ripple-rates and their correlation with cognitive-performance. This study from (24) analyzed the rates of ripples in the hippocampus and rhinal cortex during a short nap of 1 h, a set of images was presented pre- and post-nap and then again at the control stage, where patients had to distinguish known from novel images. Their results showed firstly, that the ripple events were circumscribed to the frequency range between 80 and 120 Hz. Secondly, that the ripple rate in the hippocampus was on average 1.90/min. Thirdly, that only rhinal, but not hippocampal ripples were correlated with the number of correctly recognized items. Our results differ from (24) in that our average rate of hippocampal ripples is higher (Table 3, allRipples: 13.35/min, spindleRipples: 1.14/min), which can be explained by the different detection methods used (amplitude thresholding vs. multivariate analysis). Our results slightly differ with (24) in that the frequency range of the detected ripples was circumscribed to a broader frequency range spanning between 88 and 144 Hz for the isolRipples, and 82 to 165 Hz for the iesRipples. An important agreement between our results and those from (24) is that hippocampal ripples, when undifferentiated (i.e., allRipples), do not present a correlation with performance improvement (measured by the difference in performances obtained pre-sleep and post-sleep).
The strong correlation shown by the spindleRipples with the performance improvement provides further evidence for their involvement in memory consolidation processes, moreover, these findings may contribute to the separation of physiological and non-physiological high frequency oscillations in the human hippocampus.
Limitations
This study presents a grand average of the Ripple activity during cognitive load and does not look into more local phenomena which could arise at specific time points, e.g., the isolRipples' activity dynamics at specific time intervals after cue-presentation.
The selection of spindleRipples was based on mere co-occurrence, however previous research has shown that Ripples strongly cluster around the troughs of the sleep spindles (15, 25, 26). A selection of spindleRipples while considering their clustering around the spindle trough could provide a more depurated sub-set of Ripples promoting the memory-consolidation mechanism.
Conclusions
In conclusion, the proposed method for the differentiation of physiological and pathological Ripples could help to understand the neural processes that allow the brain to execute cognitive functions such as spatial navigation and may also help to identify specific forms of ripples as biomarkers of epileptogenicity and ictogenicity. We also presented evidence supporting the role of sleep spindle-coincident ripples in memory consolidation processes, which may contribute to better understand the neural interactions allowing the storage of newly acquired information in the brain.
Data Availability Statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author/s.
Ethics Statement
The studies involving human participants were reviewed and approved by Ethics Committee of the Freiburg University Medical Center. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
DL-P and LK conducted the navigation tasks with the patients, analyzed the data, and co-wrote the manuscript. AB conducted parts of the navigation tasks with the patients and provided assistance with the EEG measurements. MD provided assistance with the EEG measurements and co-wrote the manuscript. AT helped with the data analysis. AS-B co-wrote the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the Medical Infomatics Platform for Stereoencephalography (SEEGMIP), grant agreement No.: 650003 Human Brain Project and the BrainLinks-BrainTools Cluster of Excellence funded by the German Research Foundation (DFG; EXC 1086). LK was supported by the Federal Ministry of Education and Research (BMBF; grant number 01GQ1705A), National Science Foundation (NSF) grant BCS-1724243, NIH/NINDS grant U01 NS1113198-01, and by a travel grant from the Boehringer Ingelheim Fonds (Mainz, Germany). AT was partially financed by the Austrian Science Fund (FWF), project KLI657-B31. The article processing charge was funded by the Baden-Württemberg Ministry of Science, Research and Art and the University of Freiburg in the funding programme Open Access Publishing.
Conflict of Interest
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.
References
1. Buzsáki G, Silva FL da. High frequency oscillations in the intact brain. Progr Neurobiol. (2012) 98:241–9. doi: 10.1016/j.pneurobio.2012.02.004
2. Bragin A, Engel J, Wilson C L, Fried I, Buzsáki G. High-frequency oscillations in human brain. Hippocampus. (1999) 9:137–42. doi: 10.1002/(SICI)1098-1063(1999)9:2<137::AID-HIPO5>3.0.CO;2-0
3. Jacobs J, LeVan P, Châtillon C-É, Olivier A, Dubeau F, Gotman J. High frequency oscillations in intracranial EEGs mark epileptogenicity rather than lesion type. Brain. (2009) 132:1022–37. doi: 10.1093/brain/awn351
4. Jacobs J, Zijlmans M, Zelmann R, Chatillon C-E, Hall J, Olivier A, et al. High-frequency electroencephalographic oscillations correlate with outcome of epilepsy surgery. Ann Neurol. (2010) 67:209–20. doi: 10.1002/ana.21847
5. Crépon B, Navarro V, Hasboun D, Clemenceau S, Martinerie J, Baulac M, et al. Mapping interictal oscillations >200 Hz recorded with intracranial macroelectrodes in human epilepsy. Brain. (2010) 133:33–45. doi: 10.1093/brain/awp277
6. Malinowska U, Boatman-Reich D. Cross-frequency coupling during auditory perception in human cortex. (2016). In: 38th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC). Orlando, FL (2016). p. 5521–4. doi: 10.1109/EMBC.2016.7591977
7. Worrell G A, Parish L, Cranstoun SD, Jonas R, Baltuch G, Litt B. High-frequency oscillations and seizure generation in neocortical epilepsy. Brain. (2004) 127:1496–506. doi: 10.1093/brain/awh149
8. Zijlmans M, Jacobs J, Kahn YU, Zelmann R, Dubeau F, Gotman J. Ictal and interictal high frequency oscillations in patients with focal epilepsy. Clin Neurophysiol. (2011) 122:664–71. doi: 10.1016/j.clinph.2010.09.021
9. Jacobs J, Wu JY, Perucca P, Zelmann R, Mader M, Dubeau F, et al. Removing high-frequency oscillations: a prospective multicenter study on seizure outcome. Neurology. (2018) 91:e1040–52. doi: 10.1212/WNL.0000000000006158
10. Gliske SV, Irwin ZT, Chestek C, Hegeman GL, Brinkmann B, Sagher O, et al. Variability in the location of high frequency oscillations during prolonged intracranial EEG recordings. Nat Commun. (2018) 9:2155. doi: 10.1038/s41467-018-04549-2
11. Roehri N, Pizzo F, Lagarde S, Lambert I, Nica A, McGonigal A, et al. High-frequency oscillations are not better biomarkers of epileptogenic tissues than spikes. Ann Neurol. (2018) 83:84–97. doi: 10.1002/ana.25124
12. Staley KJ, Dudek FE. Interictal spikes and epileptogenesis. Epilepsy Curr. (2006) 6:199–202. doi: 10.1111/j.1535-7511.2006.00145.x
13. Ward AA. The epileptic spike. Epilepsia. (1959) 1:600–6. doi: 10.1111/j.1528-1157.1959.tb04293.x
14. Hufnagel A, Dümpelmann M, Zentner J, Schijns O, Elger CE. Clinical relevance of quantified intracranial interictal spike activity in presurgical evaluation of epilepsy. Epilepsia. (2000) 41:467–78. doi: 10.1111/j.1528-1157.2000.tb00191.x
15. Bruder JC, Dümpelmann M, Lachner Piza D, Mader M, Schulze-Bonhage A, Julia J-L. Physiological ripples associated with sleep spindles differ in waveform morphology from epileptic ripples. Int J Neural Syst. (2016) 27:1750011. doi: 10.1142/S0129065717500113
16. Wang S, Wang IZ, Bulacio JC, Mosher JC, Gonzalez-Martinez J, Alexopoulos AV, et al. Ripple classification helps to localize the seizure-onset zone in neocortical epilepsy. Epilepsia. (2013) 54:370–6. doi: 10.1111/j.1528-1167.2012.03721.x
17. Wang S, So NK, Jin B, Wang IZ, Bulacio JC, Enatsu R, et al. Interictal ripples nested in epileptiform discharge help to identify the epileptogenic zone in neocortical epilepsy. Clin Neurophysiol. (2017) 128:945–51. doi: 10.1016/j.clinph.2017.03.033
18. Lachner-Piza D, Jacobs J, Schulze-Bonhage A, Stieglitz T, Dümpelmann M. Estimation of the epileptogenic-zone with HFO sub-groups exhibiting various levels of epileptogenicity*. In: 41st Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC). Berlin (2019). p. 2543–6. doi: 10.1109/EMBC.2019.8856733
19. Andrillon T, Nir Y, Staba RJ, Ferrarelli F, Cirelli C, Tononi G, et al. Sleep spindles in humans: insights from intracranial EEG, unit recordings. J Neurosci. (2011) 31:17821–34. doi: 10.1523/JNEUROSCI.2604-11.2011
20. Iber C, Ancoli-Israel S, Chesson AL Jr, Quan SF. The AASM Manual for the Scoring of Sleep and Associated Events: Rules, Terminology, and Technical Specifications, 1st edn. Westchester, IL: American Academy of Sleep Medicine. (2007).
21. Rechtschaffen A, Kales A, University of California LABIS. NINDB Neurological Information Network (U.S.). In: Allan Rechtschaffen and Anthony Kales, editors. A Manual of Standardized Terminology, Techniques and Scoring System for Sleep Stages of Human Subjects. Bethesda, MD: U. S. National Institute of Neurological Diseases and Blindness, Neurological Information Network. (1968). p. 6.
22. Silber MH, Ancoli-Israel S, Bonnet MH, Chokroverty S, Grigg-Damberger MM, Hirshkowitz M, et al. The visual scoring of sleep in adults. J Clin Sleep Med. (2007) 3:121–31. doi: 10.5664/jcsm.26814
23. De Gennaro L, Ferrara M. Sleep spindles: an overview. Sleep Med Rev. (2003) 7:423–40. doi: 10.1053/smrv.2002.0252
24. Axmacher N, Elger CE, Fell J. Ripples in the medial temporal lobe are relevant for human memory consolidation. Brain. (2008) 131:1806–17. doi: 10.1093/brain/awn103
25. Clemens Z, Mölle M, Eross L, Jakus R, Rásonyi G, Halász P, et al. Fine-tuned coupling between human parahippocampal ripples and sleep spindles: ripple-spindle events in human sleep. Eur J Neurosci. (2011) 33:511–20. doi: 10.1111/j.1460-9568.2010.07505.x
26. Staresina BP, Bergmann TO, Bonnefond M, Meij R van der, Jensen O, Deuker L, et al. Hierarchical nesting of slow oscillations, spindles and ripples in the human hippocampus during sleep. Nat Neurosci. (2015) 18:1679–86. doi: 10.1038/nn.4119
27. Rasch B, Born J. About sleep's role in memory. Physiol Rev. (2013) 93:681–766. doi: 10.1152/physrev.00032.2012
28. Diekelmann S, Born J. The memory function of sleep. Nat Rev Neurosci. (2010) 11:114–26. doi: 10.1038/nrn2762
29. Mölle M, Bergmann T O, Marshall L, Born J. Fast and slow spindles during the sleep slow oscillation: disparate coalescence and engagement in memory processing. Sleep. (2011) 34:1411–21. doi: 10.5665/SLEEP.1290
30. Buzsáki G. Hippocampal sharp wave-ripple: a cognitive biomarker for episodic memory and planning. Hippocampus. (2015) 25:1073–188. doi: 10.1002/hipo.22488
31. Siapas AG, Wilson MA. Coordinated interactions between hippocampal ripples and cortical spindles during slow-wave sleep. Neuron. (1998) 21:1123–8. doi: 10.1016/S0896-6273(00)80629-7
32. Buzsáki G. The hippocampo-neocortical dialogue. Cerebral Cortex. (1996) 6:81–92. doi: 10.1093/cercor/6.2.81
33. Kunz L, Schröder TN, Lee H, Montag C, Lachmann B, Sariyska R, et al. Reduced grid-cell–like representations in adults at genetic risk for Alzheimer's disease. Science. (2015) 350:430–3. doi: 10.1126/science.aac8128
34. Chen D, Kunz L, Wang W, Zhang H, Wang W-X, Schulze-Bonhage A, et al. Hexadirectional modulation of theta power in human entorhinal cortex during spatial navigation. Curr Biol. (2018) 28:3310–5.e4. doi: 10.1016/j.cub.2018.08.029
35. Kunz L, Wang L, Lachner-Piza D, Zhang H, Brandt A, Dümpelmann M, et al. Hippocampal theta phases organize the reactivation of large-scale electrophysiological representations during goal-directed navigation. Sci Adv. (2019) 5:eaav8192. doi: 10.1126/sciadv.aav8192
36. Lachner-Piza D, Jacobs J, Bruder JC, Schulze-Bonhage A, Stieglitz T, Duempelmann M. Automatic detection of high-frequency-oscillations and their sub-groups co-occurring with interictal-epileptic-spikes. J Neural Eng. (2019) 17:016030. doi: 10.1088/1741-2552/ab4560
37. Läuter J. HochbergY, Ajit C. Tamhane: multiple comparison procedures. John Wiley & Sons, New York – Chichester – Brisbane – Toronto – Singapore 1987, XXII, 450 S., £ 38.95. Biometr J. (1989) 31:122. doi: 10.1002/bimj.4710310115
38. Lachner-Piza D, Epitashvili N, Schulze-Bonhage A, Stieglitz T, Jacobs J, Dümpelmann M. A single channel sleep-spindle detector based on multivariate classification of EEG epochs: MUSSDET. J Neurosci Methods. (2018) 297:31–43. doi: 10.1016/j.jneumeth.2017.12.023
39. Himanen S-L, Virkkala J, Huhtala H, Hasan J. Spindle frequencies in sleep EEG show U-shape within first four NREM sleep episodes. J Sleep Res. (2002) 11:35–42. doi: 10.1046/j.1365-2869.2002.00273.x
40. Dijk D-J, Hayes B, Czeisler CA. Dynamics of electroencephalographic sleep spindles and slow wave activity in men: effect of sleep deprivation. Brain Res. (1993) 626:190–9. doi: 10.1016/0006-8993(93)90579-C
41. Bénar C G, Chauvire L, Bartolomei F, Wendling F. Pitfalls of high-pass filtering for detecting epileptic oscillations: a technical note on false? Ripples. Clin Neurophysiol. (2010) 121:301–10. doi: 10.1016/j.clinph.2009.10.019
43. Joo HR, Frank LM. The hippocampal sharp wave–ripple in memory retrieval for immediate use and consolidation. Nat Rev Neurosci. (2018) 19:744–57. doi: 10.1038/s41583-018-0077-1
44. Nokia MS, Mikkonen JE, Penttonen M, Wikgren J. Disrupting neural activity related to awake-state sharp wave-ripple complexes prevents hippocampal learning. Front Behav Neurosci. (2012) 6:84. doi: 10.3389/fnbeh.2012.00084
45. Papale AE, Zielinski MC, Frank LM, Jadhav SP, Redish AD. Interplay between hippocampal sharp-wave-ripple events and vicarious trial and error behaviors in decision making. Neuron. (2016) 92:975–82. doi: 10.1016/j.neuron.2016.10.028
46. Kucewicz MT, Cimbalnik J, Matsumoto JY, Brinkmann BH, Bower MR, Vasoli V, et al. High frequency oscillations are associated with cognitive processing in human recognition memory. Brain. (2014) 137:2231–44. doi: 10.1093/brain/awu149
47. Chen YY, Yoshor D, Sheth SA, Foster BL. Stability of ripple events during task engagement in human hippocampus. bioRxiv. (2020) 2020.10.17.342881. doi: 10.1101/2020.10.17.342881
48. Brázdil M, Cimbálník J, Roman R, Shaw DJ, Stead MM, Daniel P, et al. Impact of cognitive stimulation on ripples within human epileptic and non-epileptic hippocampus. BMC Neurosci. (2015) 16:47. doi: 10.1186/s12868-015-0184-0
49. Pail M, Cimbálník J, Roman R, Daniel P, Shaw DJ, Chrastina J, et al. High frequency oscillations in epileptic and non-epileptic human hippocampus during a cognitive task. Sci Rep. (2020) 10:18147. doi: 10.1038/s41598-020-74306-3
50. Cimbalnik J, Pail M, Klimes P, Travnicek V, Roman R, Vajcner A, et al. Cognitive processing impacts high frequency intracranial EEG activity of human hippocampus in patients with pharmacoresistant focal epilepsy. Front Neurol. (2020) 11:578571. doi: 10.3389/fneur.2020.578571
51. Engel J. Mesial temporal lobe epilepsy: what have we learned? Neuroscientist. (2001) 7:340–52. doi: 10.1177/107385840100700410
52. Elger CE, Helmstaedter C, Kurthen M. Chronic epilepsy and cognition. Lancet Neurol. (2004) 3:663–72. doi: 10.1016/S1474-4422(04)00906-8
53. Aarts JH, Binnie CD, Smit AM, Wilkins AJ. Selective cognitive impairment during focal and generalized epileptiform EEG activity. Brain. (1984) 107:293–308. doi: 10.1093/brain/107.1.293
54. Matsumoto JY, Stead M, Kucewicz MT, Matsumoto AJ, Peters PA, Brinkmann BH, et al. Network oscillations modulate interictal epileptiform spike rate during human memory. Brain. (2013) 136:2444–56. doi: 10.1093/brain/awt159
55. Yanagisawa T, Hirata M, Kishima H, Goto T, Saitoh Y, Oshino S, et al. Movement induces suppression of interictal spikes in sensorimotor neocortical epilepsy. Epilepsy Res. (2009) 87:12–7. doi: 10.1016/j.eplepsyres.2009.07.002
Keywords: high-frequency oscillations, ripples, interictal epileptiform spikes, sleep spindles, hippocampus, cognition, memory consolidation, spatial memory
Citation: Lachner-Piza D, Kunz L, Brandt A, Dümpelmann M, Thomschewski A and Schulze-Bonhage A (2021) Effects of Spatial Memory Processing on Hippocampal Ripples. Front. Neurol. 12:620670. doi: 10.3389/fneur.2021.620670
Received: 26 October 2020; Accepted: 04 February 2021;
Published: 05 March 2021.
Edited by:
Maeike Zijlmans, University Medical Center Utrecht, NetherlandsReviewed by:
Jan Cimbalnik, International Clinical Research Center (FNUSA-ICRC), CzechiaEleonora Tamilia, Boston Children's Hospital and Harvard Medical School, United States
Copyright © 2021 Lachner-Piza, Kunz, Brandt, Dümpelmann, Thomschewski and Schulze-Bonhage. 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) and the copyright owner(s) 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: Daniel Lachner-Piza, ZGFuaWVsLmxhY2huZXImI3gwMDA0MDt1bmlrbGluaWstZnJlaWJ1cmcuZGU=