Skip to main content

Main menu

  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
    • Podcast
  • ALERTS
  • FOR AUTHORS
    • Information for Authors
    • Fees
    • Journal Clubs
    • eLetters
    • Submit
  • EDITORIAL BOARD
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
  • SUBSCRIBE

User menu

  • Log in
  • My Cart

Search

  • Advanced search
Journal of Neuroscience
  • Log in
  • My Cart
Journal of Neuroscience

Advanced Search

Submit a Manuscript
  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
    • Podcast
  • ALERTS
  • FOR AUTHORS
    • Information for Authors
    • Fees
    • Journal Clubs
    • eLetters
    • Submit
  • EDITORIAL BOARD
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
  • SUBSCRIBE
PreviousNext
Research Articles, Behavioral/Cognitive

Anterior–Posterior Hippocampal Dynamics Support Working Memory Processing

Jin Li (李瑾), Dan Cao (曹丹), Vasileios Dimakopoulos, Weiyang Shi (时维 阳), Shan Yu (余山), Lingzhong Fan (樊令 仲), Lennart Stieglitz, Lukas Imbach, Johannes Sarnthein and Tianzi Jiang (蒋田 仔)
Journal of Neuroscience 19 January 2022, 42 (3) 443-453; DOI: https://doi.org/10.1523/JNEUROSCI.1287-21.2021
Jin Li (李瑾)
1Brainnetome Center, Institute of Automation, Chinese Academy of Sciences, Beijing 100190, People's Republic of China
2National Laboratory of Pattern Recognition, Institute of Automation, Chinese Academy of Sciences, Beijing, 100190, People's Republic of China
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Jin Li (李瑾)
Dan Cao (曹丹)
1Brainnetome Center, Institute of Automation, Chinese Academy of Sciences, Beijing 100190, People's Republic of China
2National Laboratory of Pattern Recognition, Institute of Automation, Chinese Academy of Sciences, Beijing, 100190, People's Republic of China
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Vasileios Dimakopoulos
3Department of Neurosurgery, University Hospital Zurich, University of Zurich, 8091 Zurich, Switzerland
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Weiyang Shi (时维 阳)
1Brainnetome Center, Institute of Automation, Chinese Academy of Sciences, Beijing 100190, People's Republic of China
2National Laboratory of Pattern Recognition, Institute of Automation, Chinese Academy of Sciences, Beijing, 100190, People's Republic of China
4School of Artificial Intelligence, University of Chinese Academy of Sciences, Beijing 100049, People's Republic of China
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Shan Yu (余山)
1Brainnetome Center, Institute of Automation, Chinese Academy of Sciences, Beijing 100190, People's Republic of China
2National Laboratory of Pattern Recognition, Institute of Automation, Chinese Academy of Sciences, Beijing, 100190, People's Republic of China
4School of Artificial Intelligence, University of Chinese Academy of Sciences, Beijing 100049, People's Republic of China
5Center for Excellence in Brain Science and Intelligence Technology, Chinese Academy of Sciences, Beijing 100190, People's Republic of China
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Lingzhong Fan (樊令 仲)
1Brainnetome Center, Institute of Automation, Chinese Academy of Sciences, Beijing 100190, People's Republic of China
2National Laboratory of Pattern Recognition, Institute of Automation, Chinese Academy of Sciences, Beijing, 100190, People's Republic of China
4School of Artificial Intelligence, University of Chinese Academy of Sciences, Beijing 100049, People's Republic of China
5Center for Excellence in Brain Science and Intelligence Technology, Chinese Academy of Sciences, Beijing 100190, People's Republic of China
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Lennart Stieglitz
3Department of Neurosurgery, University Hospital Zurich, University of Zurich, 8091 Zurich, Switzerland
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Lukas Imbach
6Swiss Epilepsy Center, 8008 Zurich, Switzerland
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Johannes Sarnthein
3Department of Neurosurgery, University Hospital Zurich, University of Zurich, 8091 Zurich, Switzerland
7Zurich Neuroscience Center, ETH Zurich, 8057 Zurich, Switzerland
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Johannes Sarnthein
Tianzi Jiang (蒋田 仔)
1Brainnetome Center, Institute of Automation, Chinese Academy of Sciences, Beijing 100190, People's Republic of China
2National Laboratory of Pattern Recognition, Institute of Automation, Chinese Academy of Sciences, Beijing, 100190, People's Republic of China
4School of Artificial Intelligence, University of Chinese Academy of Sciences, Beijing 100049, People's Republic of China
5Center for Excellence in Brain Science and Intelligence Technology, Chinese Academy of Sciences, Beijing 100190, People's Republic of China
8The Queensland Brain Institute, University of Queensland, Brisbane, Queensland 4072, Australia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Tianzi Jiang (蒋田 仔)
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF
Loading

Abstract

The hippocampus is a locus of working memory (WM) with anterior and posterior subregions that differ in their transcriptional and external connectivity patterns. However, the involvement and functional connections between these subregions in WM processing are poorly understood. To address these issues, we recorded intracranial EEG from the anterior and the posterior hippocampi in humans (seven females and seven males) who maintained a set of letters in their WM. We found that WM maintenance was accompanied by elevated low-frequency activity in both the anterior and posterior hippocampus and by increased theta/alpha band (3–12 Hz) phase synchronization between anterior and posterior subregions. Cross-frequency and Granger prediction analyses consistently showed that the correct WM trials were associated with theta/alpha band-coordinated unidirectional influence from the posterior to the anterior hippocampus. In contrast, WM errors were associated with bidirectional interactions between the anterior and posterior hippocampus. These findings imply that theta/alpha band synchrony within the hippocampus may support successful WM via a posterior to anterior influence. A combination of intracranial recording and a fine-grained atlas may be of value in understanding the neural mechanisms of WM processing.

SIGNIFICANCE STATEMENT Working memory (WM) is crucial to everyday functioning. The hippocampus has been proposed to be a subcortical node involved in WM processes. Previous studies have suggested that the anterior and posterior hippocampi differ in their external connectivity patterns and gene expression. However, it remains unknown whether and how human hippocampal subregions are recruited and coordinated during WM tasks. Here, by recording intracranial electroencephalography simultaneously from both hippocampal subregions, we found enhanced power in both areas and increased phase synchronization between them. Furthermore, correct WM trials were associated with a unidirectional influence from the posterior to the anterior hippocampus, whereas error trials were correlated with bidirectional interactions. These findings indicate a long-axis specialization in the human hippocampus during WM processing.

  • hippocampus
  • intracranial EEG
  • longitudinal axis
  • working memory

Introduction

Working memory (WM), the “sketchpad of conscious thought” (Miller et al., 2018), is critical for organizing goal-directed behaviors (D'Esposito, 2007). It describes our capacity to temporarily hold information and adaptively translate it later into a behavioral response. The hippocampus (Hipp) has been proposed to be a subcortical node underlying WM. Direct electrophysiological recordings in the human hippocampus consistently found persistent single-neuron activity (Kamiński et al., 2017; Kornblith et al., 2017; Boran et al., 2019), power change (Brzezicka et al., 2019), and increased phase–amplitude coupling (PAC; Axmacher et al., 2010) during WM maintenance. These studies modeled the hippocampus as a single, homogeneous entity; however, the hippocampus is a complex, heterogeneous structure consisting of distinct, interacting subregions. Rodent studies suggested that the ventral (homologous to the anterior in humans) and dorsal (homologous to the posterior in humans) hippocampi differ in their connectivity and gene expression patterns (Fanselow and Dong, 2010; Poppenk et al., 2013). Such connectivity and transcriptomic heterogeneity have recently been validated in humans (Vos de Wael et al., 2018; Vogel et al., 2020) and provide a rationale for the functional differentiation along the long axis (Poppenk et al., 2013; Strange et al., 2014). Consistent with this hypothesis, Lin et al. (2017, 2019) reported that during episodic memory encoding, the positive subsequent memory effect in the slow-theta range is greater in the posterior hippocampus (Lin et al., 2017); later, they found that gamma activity was greater during retrieval than during encoding only in the posterior hippocampus (Lin et al., 2019). Given that mounting evidence shows functional organization along the longitudinal axis, treating the hippocampus as a functionally uniform structure in WM processing may disregard potentially useful information about its subregions.

Rodent experiments suggested a critical role of the dorsal but not the ventral hippocampus in WM processing, with a significant WM deficit caused by dorsal but not ventral hippocampal lesions (Pothuizen et al., 2004). A recent meta-analysis of functional magnetic resonance imaging (MRI) studies involving humans reported activation of both the anterior Hipp (aHipp) and the posterior Hipp (pHipp) during WM processing (Grady, 2020). The limited temporal and spatial resolution of this noninvasive MRI technique in human studies has made it difficult to capture the functional specialization of the hippocampus at a fine-grained scale. Thus, direct brain recordings precisely located in hippocampal subregions during WM tasks are required. In addition, the anatomic connectivity suggests a flow of information between the aHipp and the pHipp (Schmidt et al., 2013). Rodent studies showed increased connectivity between the dorsal and ventral hippocampi as they engage in WM processing (Schmidt et al., 2013). However, it remains unknown whether and how the two hippocampal subregions communicate during WM in humans. Human intracranial electroencephalography (iEEG) recordings, with both millimeter spatial precision and millisecond temporal resolution, even in deep brain regions (Johnson et al., 2020), offer the possibility of delineating activity at the hippocampal subregion level and of studying the information flow between the two hippocampal subregions in vivo (Lin et al., 2017).

Here, we recorded iEEG simultaneously from the aHipp and the pHipp in 14 presurgical epilepsy patients while they performed a modified Sternberg task. We addressed the following questions. (1) Does the hippocampus play a single role in WM processing or is there spatial heterogeneity along the long axis in its involvement in working memory? (2) Are the aHiPP and the pHipp independent or do they interact in WM? (3) Do these activities or interactions support WM maintenance?

Materials and Methods

Subjects.

We used data from 14 adult human subjects (mean ± SD age, 32 ± 12 years; age range, 18–56 years; seven males) in this study. All subjects were implanted with depth electrodes [diameter, 1.3 mm; 8 contacts of 1.6 mm length; 5 mm spacing between contact centers; Ad-Tech (www.adtechmedical.com)] in the medial temporal lobe for evaluation of the surgical treatment of epilepsy. Electrode placement was exclusively guided by clinical needs, and patient selection was based on MRI-confirmed depth electrode placement in the hippocampus contralateral to or outside of the epileptogenic region. Before testing, all the subjects provided written informed consent for the study, which had been approved by the relevant institutional ethics review board (Kantonale Ethikkommission Zürich, Approval PB 2016–02055). There were no seizures recorded during any of the epochs, and any epochs with interictal epileptiform activity were excluded from analysis.

Experimental design and statistical analyses.

WM was tested using a modified Sternberg task in which the encoding of memory contents, maintenance, and recall were temporally separated (Fig. 1A). Each trial started with a fixation period (1 s) followed by the stimulus for 2 s. The memory content consisted of a set of letters at the center of the screen. After the disappearance of the stimulus, the maintenance interval was followed by a fixation square (3 s). After that, the subjects responded with a button press (“IN” or “OUT”) as fast as possible when a probe was presented. They were instructed to recall whether the probe was part of the stimulus. After the response, the subjects were encouraged to relax, and the next trial was initiated. Each subject performed several sessions with 50 trials per session (226 ± 113 trials per subject), and the order of the stimuli was fully random across sessions.

We then statistically compared the z-scored power between the aHipp and the pHipp using a cluster-based permutation test, which computes statistics at the cluster level and corrects for multiple comparisons (Maris and Oostenveld, 2007). The differences between the aHipp and the pHipp were quantified by t values. The t values corresponding to uncorrected p values of ≤0.05 were clustered. The differences were considered significant when the maximum of the cluster-level summed t values in the true data exceeded the threshold of a Monte Carlo distribution, which was created by randomly shuffling the labels of the outcomes 1000 times, and p < 0.05 was considered as significant. Such a comparison was also applied in the phase–amplitude coupling analysis to assess the task-induced effects from the maintenance period over the baseline (p < 0.05). In addition, we examined the statistical significance of the phase locking values (PLVs) of the aHipp–pHipp electrode pairs by permuting the trials 200 times to create a null distribution and computed the 95th percentile threshold for each pair. The PLVs above the thresholds were considered to be significant (p < 0.05). A similar permutation test was also applied on the Granger causality (GC) analysis. We used mixed-effect models (MEMs) to compare the electrophysiological indexes for two directions of information transfer between the aHipp and the pHipp. All data analyses were conducted using MATLAB (MathWorks).

Depth electrode localization.

The electrodes were localized using postimplantation computed tomography (CT) scans and postimplantation structural T1-weighted MRI scans. For each subject, the CT scan was coregistered to the postimplantation scan, as implemented in FieldTrip (Oostenveld et al., 2011). The electrode contacts were visually marked on the coregistered CT-MR images. To determine the electrode locations, we used a high-resolution anatomic template (1 mm) developed in our laboratory (Fan et al., 2016), which parcellates the brain into 246 subregions based on their anatomic connectivity patterns. This template has been widely used in past studies (Stolk et al., 2018; Assaneo et al., 2019; Lee et al., 2021). In this atlas, the hippocampus is divided into anterior and posterior subregions along the long axis (http://atlas.brainnetome.org/).

The final dataset contained 94 contacts in the hippocampus including 32 contacts in the aHipp and 62 contacts in the pHipp across all subjects. There were 3.2 ± 1.8 contacts per subject in the aHipp (range, 1–6; 10 of 14 subjects) and 6.7 ± 2.2 contacts per subject in the pHipp (range, 1–8). Thus, the data from 10 subjects with electrodes in both the aHipp and the pHipp were included in the interregional connectivity analysis. Figure 1B presents the electrodes transformed into standard MNI space in the aHipp and pHipp across all subjects.

Data acquisition and preprocessing.

Intracranial data were acquired using a Neuralynx ATLAS recording system, sampled at 4 kHz, and analog filtered >0.5 Hz against a common intracranial reference. After data acquisition, neural recordings were downsampled to 1 kHz and bandpass filtered from 1 to 200 Hz using the zero-phase delay FIR (finite impulse response) filter with Hamming window. Line noise harmonics were removed using a discrete Fourier transform. The filtered data were manually inspected to mark any contacts or epochs containing epileptiform activity or artifacts for exclusion. The data were then rereferenced to the average of the signal over all the contacts. The contacts localized in the hippocampus were then selected for subsequent analysis. Channels were excluded when they had continuous artifacts, usually because the cable was broken. In our analysis, only one contact was excluded for this reason. In addition, for all the participants, 1.4% of the trials (n = 43/3050) were excluded because of large singular artifacts, usually because the cable moved.

The continuous data were segmented into 8 s trials with a 1 s fixation period as the baseline, 2 s encoding period, 3 s maintenance period, and 2 s retrieval period. We focused on the maintenance period. The trials with residual artifacts were rejected via visual inspection. The subjects' responses defined correct and incorrect trials. Preprocessing routines were performed using the FieldTrip (Oostenveld et al., 2011), EEGLAB (Delorme and Makeig, 2004) toolboxes and customized scripts for MATLAB.

Frequency decomposition.

Time–frequency representations of power were computed for correct and incorrect trials. For each trial for each electrode, we convolved the signal within the aHipp and the pHipp with complex-valued Morlet wavelets (six cycles) to obtain power information (Solomon et al., 2019). The frequency range was 2–128 Hz at 23 logarithmically spaced frequencies with a time resolution of 1 ms. The task-induced power was analyzed per trial using the same statistical bootstrapping procedure as used in previous studies (Johnson et al., 2018). Briefly, for each channel and frequency, the baseline power (500 ms pretrial) was pooled into a single time series; then r data points (r = number of trials of the subject) were randomly selected and averaged. This procedure was repeated 1000 times to create a null distribution. The raw power for each time point was then z-scored by comparing it to the entire null distribution. The task-induced power was considered significant when the z-scored power exceeded the value of 1.96 (p < 0.05). Then, the z-scored power for the aHipp and the pHipp was compared at each time–frequency point with a cluster-based permutation test (Maris and Oostenveld, 2007) for the correct and the incorrect trials, separately.

Interregional phase synchrony.

The strength of the interregional neural synchrony was quantified by the PLV, which calculates the average phase φ difference between two electrodes (a,b) across trials for a given time point t and frequency f, as follows: PLVa,b(t,f)=1Ntrials|∑n=1n=Ntrialsexp(i[φn,a(t,f)−φn,b(t,f)])|.(1)

The PLV indexes the degree of consistency of the phase relationship for each pair of electrodes independent of the absolute phases and amplitudes among repeated trials and ranges from 0 to 1, with values approaching 1 if the two signals show a constant phase relationship. We calculated the PLV for each electrode pair of the aHipp–pHipp from 1 to 30 Hz across the correct and the incorrect trials, separately. The trialwise mean was subtracted from each trial to minimize contamination from simultaneous voltage changes of the phase estimates before calculating the PLV (Brincat and Miller, 2015).

To test the statistical significance of the PLVs of all electrode pairs within the maintenance trials, we permuted the trials 200 times to create a null distribution and computed the 95th percentile threshold for each pair. The thresholds were then averaged across all electrode pairs and only the time–frequency PLVs above the corresponding thresholds were kept as significant PLVs. The values below the thresholds were set to zero. To further specify the frequency-specific PLVs, the spectral PLVs were obtained across time. We selected and extracted the peak PLVs (theta/alpha, 3–12 Hz) to compare the PLVs with the baseline using a paired t test to confirm that this peak was associated with the maintenance process.

Granger causality analysis.

To examine the directionality of the aHipp–pHipp synchronization within the theta/alpha band, we analyzed the effective aHipp–pHipp connectivity using spectral GC, which quantifies the prediction error of the signal in the frequency domain by introducing another time series. Based on observations from previous studies (Kota et al., 2020), we extracted the GC values within 3–12 Hz that showed task-related elevated PLV. The aHipp–pHipp electrode pairs with task-induced elevated connectivity (theta/alpha PLVmaintenance/PLVbaseline > 1) for the correct and incorrect trials were selected separately. Following this criterion, 49 of 84 electrode pairs for correct trials and 40 pairs for incorrect trials were selected separately for later analysis. For each electrode pair of the aHipp–pHipp, the trialwise mean was subtracted from each trial before being fit to an autoregressive model and computing the spectral GC. We then applied the Multivariate Granger Causality MATLAB Toolbox (Barnett and Seth, 2014) based on the Akaike information criterion to define the model order for each pair. The GC index of both directions (from aHipp to pHipp and the reverse direction) was computed during the maintenance period for the correct and incorrect trials, separately. The significance testing for the spectral GC was conducted by creating a null distribution by randomly swapping the signal between the electrodes 200 times. The electrode pairs with a GC above the threshold, defined as the 95th percentile of the null distribution, were selected for the subsequent comparisons.

The theta/alpha GC for each pair was extracted in both directions for the correct trials and the incorrect trials, separately. Then, two MEMs were used to compare the GC value from both directions, for the correct and incorrect trials. In each MEM, we treated the direction as a fixed factor, the electrode pair number and subjects as random factors, and the extracted theta/alpha GC values for each electrode pair as the dependent variables.

Given that the GC analysis is sensitive to the signal-to-noise ratio across frequency bands, we also assessed the theta/alpha band directionality using the phase slope index (PSI; Basti et al., 2018) from both directions for the correct and incorrect trials, separately. Again, the electrode pairs with elevated PLVs for the correct and incorrect trials were selected separately. For each trial, the data segments during maintenance were zero padded and multiplied with a Hann taper from 3 to 12 Hz with 1 Hz steps, from which we computed the theta/alpha PSI (Johnson et al., 2018). To correct for any spurious results, null distributions were generated by randomly shuffling the trials 1000 times. Raw PSI outputs were then z-scored on the null distributions to assess the significance of the directional effects. Significant PSI was thresholded at |z| > 1.96, in which the aHipp leads were defined as z >1.96 and the pHipp leads as z < −1.96.

Phase–amplitude coupling.

The cross-frequency coupling between the lower frequency (1–20 Hz) phase from contact a and the higher-frequency (20–100 Hz) amplitude from contact b was examined individually by mean vector length (Canolty et al., 2006). Phase–amplitude coupling was calculated according to the following equation: PAC=|∑t=1nAt,aeiφt,bn|,(2) where n is the number of data points, t is a data point, At,a is the amplitude, and φt,b is the phase angle.

Then, the raw PAC was assessed against surrogate distributions via bootstrapping as was done in previous studies (Canolty et al., 2006). Specifically, we shuffled data by inserting a randomly generated time lag between the time series of the low-frequency phase and the high-frequency amplitude and recomputed the PAC. This procedure was repeated 200 times to create surrogate distributions of PAC values. Then, the raw PAC was z-scored based on the null distributions to correct for any spurious results. The z-scored PACs exceeding a threshold value of 1.96 were considered as significant (Wang et al., 2021). To assess the task-induced effects, the z-scored PAC during the maintenance was compared with the baseline values by the cluster-based permutation test used in the previous analysis. A one-tailed test was applied to find the task-increased coupling relative to the baseline (p < 0.05).

We extracted the theta/alpha (3–12 Hz)–gamma (30–100 Hz) z-scored PAC from the following two directions: the theta/alpha phase of the aHipp modulating the gamma amplitude of the pHipp as well as the theta/alpha phase of the pHipp modulating the gamma amplitude of the aHipp for the correct and incorrect trials, separately. The outputs were analyzed using an MEM, with direction as a fixed factor, subject and electrode pair as random factors, and z-scored PAC for each electrode pair as a dependent variable.

Data and code availability.

The dataset has been analyzed and described previously (Boran et al., 2019, 2020, 2021) and is freely available for download at https://doi.gin.g-node.org/10.12751/g-node.d76994/. The task is freely available for download at http://www.neurobs.com/ex_files/expt_view?id=266. Links to updates and further datasets can be found at https://hfozuri.ch.

Results

Task and behavior

Fourteen presurgical epilepsy patients (7 males; Table 1) performed a modified Sternberg WM task (58 total sessions from 14 subjects). In this task, items were presented simultaneously rather than sequentially, thus separating the encoding stage from the maintenance stage. In each trial, the subject was asked to memorize a set of 4, 6, or 8 letters presented for 2 s (encoding). After a delay (maintenance) period of 3 s, the subjects were instructed to judge whether a probe letter was identical to one of the letters held in memory (retrieval; Fig. 1A).

View this table:
  • View inline
  • View popup
Table 1.

Subject characteristics

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1.

Working memory task, recording sites in the hippocampus, and time–frequency power. A, An example trial of the task. Each trial consisted of a set of consonants (encoding, 2 s), followed by a delay (maintenance, 3 s). After the delay, a probe letter was shown, and the subjects indicated whether the probe was or was not shown during the encoding period (retrieval, 2 s). B, Electrode location across subjects in Montreal Neurologic Institute (MNI152) space. Recording subregions are indicated by different colors (magenta, aHipp; cyan, pHipp). C, An example MRI and template for a single subject. D, Average task-induced power, grouped in the anterior hippocampus (left column) and posterior hippocampus (right column), across all subjects during the maintenance phase for the correct trials. Warmer colors denote task-increased power above the baseline, and colder colors denote task-decreased power. Both the aHipp and the pHipp showed sustained task-increased power in the low-frequency range (2–20 Hz; z > 1.96, p < 0.05) during the maintenance phase. E, Spectral z-scored power within the anterior (magenta) and posterior hippocampus (cyan) across all subjects (±SEM shown as shading around the mean trace) for the correct trials. During maintenance, no significant difference in the spectral power (p > 0.05, cluster-based permutation test) was found between the two subregions across the frequency band.

The average accuracy was 90.8 ± 1.2% (range, 78−94.9%). The mean response time was faster for correct than incorrect trials (1.32 ± 0.04 vs 1.81 ± 0.15 s; paired t test: t(13) = −4.0489, p = 0.0005). Hence, the subjects were able to perform the WM task.

Low-frequency activity is sustained during maintenance

Local field potentials were recorded from the depth electrodes implanted in the aHipp (32 contacts) and the pHipp (62 contacts; Fig. 1B). The average distance between the anterior and posterior contacts was 19.59 ± 0.80 mm (range, 4.69–36.67 mm). An example MRI and template were also presented for a single subject (Fig. 1C). We calculated the time–frequency power for the correct trials/for each contact in each of the subregions during the maintenance period, and the power outputs were z-scored against pretrial baseline distributions to assess the significance of the task-induced power effects per trial. Both the aHipp and the pHipp showed sustained activity in the low-frequency range (2–20 Hz; z >1.96, p < 0.05; Fig. 1D). Then we averaged the spectral power across time and compared it between the aHipp and the pHipp at each frequency. We found no aHipp–pHipp difference in the spectral power (p > 0.05, clustered-based permutation test) across the frequency bands for the correct trials (Fig. 1E). Together, both the aHipp and the pHipp engaged in the maintenance of the WM information, as reflected by the low-frequency oscillations.

Theta/alpha synchronization of aHipp–pHipp supports the maintenance of WM information

Task-induced increases in spectral power were found within the low-frequency bands for both the aHipp and the pHipp. This raised a further question about whether the two subregions of the hippocampus interact during WM maintenance. Given the role of low-frequency oscillations in supporting interregional communication, we calculated the low-frequency phase synchrony with PLVs between the aHipp and the pHipp for the correct trials, thus measuring the consistency of the phase relationship for each aHipp–pHipp electrode pair. The PLVs up to 30 Hz were computed across the time–frequency domain to reveal the dynamic fluctuations throughout the functional connectivity during the maintenance period. This analysis was done using the data from 10 participants who had electrodes in both the anterior and the posterior hippocampi.

As shown in Figure 2A, evident phase synchronization (p < 0.05, permutation test) between aHipp and pHipp was found throughout the entire maintenance period. Further, we obtained the spectral PLVs across the time domain and confirmed that the PLVs of aHipp–pHipp peaked in the theta/alpha band (3–12 Hz; Fig. 2B, light orange shadowed). To confirm whether this peak in the PLVs was maintenance related, we extracted the PLVs within the theta/alpha band for the maintenance period and the baseline fixation, and compared them using a paired t test. For all the aHipp–pHipp electrode pairs, the PLVs for the maintenance phase were significantly stronger than those for the baseline fixation (paired t test: t(83) = 3.97, p = 0.0002), confirming that the peak of the PLVs within the theta/alpha band was strongly associated with the maintenance process (Fig. 2C). Therefore, the subsequent analyses focused on the theta/alpha band.

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2.

Frequency-specific aHipp–pHipp phase synchrony during maintenance. A, Phase synchrony (PLV) between the anterior hippocampus and the posterior hippocampus was identified across all subjects, with greater theta/alpha synchrony during the maintenance for the correct trials. The PLVs ranged from 0 to 1, with warmer colors indicating greater PLVs. The PLV maps show the PLVs that survived the threshold at p < 0.05 (cluster-based permutation test). B, Spectral PLVs within 1–30 Hz between the aHipp and the pHipp across subjects for the correct trials (±SEM shown as shaded area around the mean trace) with peaks in the theta/alpha band (3–12 Hz, shading in light orange). The shaded area indicates the theta/alpha range used for the subsequent analyses. C, Extracted theta/alpha PLVs of the aHipp–pHipp across all subjects for the baseline and maintenance in the correct trials. The theta/alpha PLVs were elevated by the maintenance relative to the baseline (paired t test, t(83) = 3.97, p = 0.0002). D, The GC index for the correct trials was extracted within the theta/alpha band across subjects for two directions (magenta: aHipp to pHipp; cyan: pHipp to aHipp). Stronger GC was found for pHipp to aHipp (cyan) than for aHipp to pHipp (magenta; mixed-effect model, *p < 0.05).

pHipp drives aHipp in the theta/alpha band supporting WM processing

To further examine the directionality of the aHipp–pHipp synchronization within the theta/alpha band, we analyzed the effective aHipp–pHipp connectivity using GC in the frequency domain. GC estimates the degree to which the signal from one region can be predicted by incorporating information from another. Next, we compared the GC values in the pHipp-to-aHipp direction with those from the reverse direction. We found that the theta/alpha-driven GC values were greater in the pHipp-to-aHipp direction than those in the reverse direction (mixed-effect model: t = 2.10, p = 0.039; Fig. 2D). These results provided evidence that directional aHipp–pHipp interactions driven by the pHipp in the theta/alpha band supported successful WM processing.

In addition to the GC analysis, we assessed the theta/alpha band directionality between the aHipp and the pHipp by calculating the PSI within the theta/alpha band (3–12 Hz). We found a unidirectional pHipp-driven connectivity (z = −2.78, p = 0.0054) for the correct trials and bidirectional aHipp–pHipp interactions (z = 0.89, p = 0.37) for incorrect trials in the theta/alpha band. This result confirmed the result from the GC analysis, and together they indicated that pHipp-driven information flow supports successful WM processing.

Cross-frequency information transfer between the aHipp and the pHipp

PAC serves as an important mechanism for coordinating interregional information transfer in that high-frequency amplitudes are modulated by phases with a low frequency (Canolty et al., 2006). It has been suggested that PAC, especially theta/alpha-gamma PAC within the hippocampus, reflects a neural correlate for WM maintenance in animals and humans (Tort et al., 2008; Axmacher et al., 2010; Roux and Uhlhaas, 2014). We computed the cross-regional PAC in both phase–amplitude combinations between the aHipp and the pHipp for the correct trials. The interregional influence was quantified as the coupling of low-frequency phases (2–20 Hz; Martínez-Cancino et al., 2019) from one region with high-frequency amplitudes (20–100 Hz) to another. The raw PAC was z-scored against surrogate distributions via bootstrapping. The results showed that the theta/alpha phase of the aHipp entrained the gamma amplitude of the pHipp (∼30–100 Hz; peak, 78 Hz; Fig. 3A, top) and the theta/alpha phase of the pHipp modulated the gamma amplitude of the aHipp (∼30–100 Hz; peak, 80 Hz; Fig. 3A, bottom). To further compare the coupling strength between the two combinations, the cross-regional PAC of the theta/alpha–gamma bands (Fig. 3A, orange box) was extracted for both combinations. We found that the coupling strength between the pHipp phase and the aHipp amplitude was greater than that between the aHipp phase and the pHipp amplitude (mixed-effect model, t = 2.84, p = 0.006), as shown in Figure 3B.

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3.

Cross-frequency coupling between the aHipp and the pHipp. A, Average cross-frequency phase–amplitude coupling (z score) between the anterior hippocampus and posterior hippocampus across all subjects for correct trials. To assess the task-induced effects, the z-scored PAC during the maintenance was compared with the baseline values by the cluster-based permutation test. A one-tailed test was applied to find the task-increased coupling relative to the baseline (p < 0.05). For the theta–alpha phase (3–12 Hz)/gamma amplitude (30–100 Hz) coupling shown in the orange boxes, interregional modulations for both directions were significantly higher during maintenance relative to those in the baseline period (all p values < 0.05). B, The theta/alpha–gamma (orange box) PACs were calculated across subjects for two directions (magenta, aHipp–pHipp means aHipp phase modulating pHipp amplitudes; cyan, pHipp–aHipp means the reverse direction) for the correct trials. The modulation strength from the pHipp to the aHipp was greater than that from the aHipp to the pHipp (mixed-effect model, **p < 0.01). C, Average z-scored intrahippocampal PAC for the slow-theta (3–5 Hz)–gamma, fast-theta (5–9 Hz)–gamma, and alpha (9–12 Hz)–gamma PACs between the aHipp and the pHipp in both phase–amplitude combinations (low-frequency phase from the aHipp and alpha amplitude from the pHipp, and vice versa). For each frequency band, the coupling strength between the pHipp phase and the aHipp amplitude (pHipp–aHipp, cyan) was consistently greater than that between the aHipp phase and the pHipp amplitude (aHipp–pHipp, magenta). D, Computation of theta/alpha–gamma interactions. Phase–amplitude distributions were constructed for both directions (top: magenta, gamma amplitudes of the pHipp distributed across the theta/alpha phase bins of the aHipp, abbreviated as aHipp–pHipp; cyan, gamma amplitudes of the aHipp distributed across the theta/alpha phase bins of the pHipp, abbreviated as pHipp–aHipp). The difference distributions were obtained by subtracting the aHipp–pHipp from the pHipp–aHipp (bottom, orange). The gamma activity from both the aHipp (magenta) and the pHipp (cyan) was modulated by the theta/alpha oscillations around the trough (Rayleigh test, average pHipp gamma amplitude occurred at 57.3°, ppHipp = 4.31 × 10−8; average aHipp gamma amplitude occurred at 57.2°, paHipp = 4.31 × 10−8). The difference distributions were not uniform (Rayleigh test, p = 4.31 × 10−8).

Based on previously reported aHipp/pHipp differences in preferred frequencies of the hippocampal cortex PAC in episodic memory processing (Wang et al., 2021), we further divided the theta/alpha band (3–12 Hz) into three frequency bands in our PAC analysis. We separately extracted the phase information for the slow-theta (3–5 Hz), fast-theta (5–9 Hz), and alpha (9–12 Hz) frequency bands in one hippocampal subregion and the amplitude information for the gamma band in the other hippocampal subregion. Then, the raw PAC outputs were z-scored against surrogate distributions. The coupling strength between the pHipp phase and the aHipp amplitude was consistently greater than that between the aHipp phase and the pHipp amplitude for these frequency bands (Fig. 3C). These results further demonstrated that the pHipp slow oscillations consistently directed the aHipp gamma during working memory processing.

Low-frequency oscillations are thought to organize cell assemblies to maintain WM information, and recent studies indicate that neuronal spiking in the hippocampus is locked to distinct phases of slow oscillations in a behaviorally relevant manner (Kamiński et al., 2020). As a surrogate for neuronal spiking, we examined when gamma activity from the modulated signals occurred relative to the theta/alpha phases of the modulating signal. The distributions of the gamma amplitudes across 18 phase bins were calculated for both directions. As shown in Figure 3D, top, the gamma activity from both the aHipp (magenta) and the pHipp (cyan) was modulated by the theta/alpha oscillations around the trough (Rayleigh test: average pHipp gamma amplitude occurred at 57.3°, ppHipp = 4.3082 × 10−8; average aHipp gamma amplitude occurred at 57.2°, paHipp = 4.3085 × 10−8). Moreover, the distributions of the gamma amplitudes in the pHipp across the theta/alpha phase bins in the aHipp were subtracted from those of the gamma amplitudes of the aHipp across the theta/alpha phase bins of the pHipp, as shown in Figure 3D, bottom. The differences were not uniform (Rayleigh test, p = 4.308 × 10−8). In summary, WM maintenance increased cross-frequency coupling between the aHipp and the pHipp at frequency- and phase-specific oscillations, and the gamma amplitudes during the theta/alpha phase oscillations optimized the information transfer from the posterior to the anterior hippocampus.

Directionality between aHipp and pHipp is unbiased for the incorrect trials

To investigate whether the findings above were specific for successful WM maintenance, we also analyzed the directionality of the information transfer between the aHipp and the pHipp for the incorrect trials. Within the theta/alpha band, significant GC above the threshold was found within the incorrect trials and no difference in GC was found for either direction between the aHipp and the pHipp (mixed-effect model: t = 1.668, p = 0.12), as shown in Figure 4A. To exclude the possibility that the lack of effect of direction for the incorrect trials are because of different electrode pairs used in statistical analysis, we repeated our analysis by using all 84 electrode pairs to compare the GC values from both directions for the correct and incorrect trials, separately. The GC values from pHipp to aHipp were significantly higher than those from the opposite direction (mixed-effect model: t = 2.06, p = 0.043) for the correct trials, while there was no significant difference between the two directions (mixed-effect model: t = 0.46, p = 0.65) for the incorrect trials. These results confirmed our findings. Therefore, the directionality of information transfer driven by the pHipp to the aHipp supported WM maintenance (Fig. 4B, left), which was beneficial for correct WM performance. However, this bias in the directionality of information transfer was absent in the hippocampus in incorrect trials (Fig. 4B, right).

Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4.

Directionality of the aHipp–pHipp for the incorrect trials. A, The GC index for the incorrect trials was extracted within the theta/alpha band across subjects for two directions (magenta, aHipp to pHipp; cyan, pHipp to aHipp). No difference in GC was found between the aHipp and the pHipp for either direction (mixed-effect model, p > 0.05). B, Schematics of the maintenance of WM processing within the hippocampus for the correct trials (left) and the incorrect trials (right). Directional information flow driven by the pHipp to the aHipp within the theta/alpha band subserved the correct performance in the WM trial, but this bias with respect to the directionality of information transfer was absent in incorrect WM trials.

As we did for the analyses conducted for the correct trials, we also computed the time–frequency power during the maintenance period for the incorrect trials within the aHipp and pHipp separately. As shown in Figure 5A, task-induced power for the incorrect trials on both the aHipp and the pHipp was not sustained during the maintenance in the low-frequency range (2–20 Hz). The spectral power was not significant (z < 1.96) during the maintenance period and no difference (p > 0.05, clustered-based permutation test) was found between the two subregions across the frequency band (Fig. 5B). In addition, the PLV across the time–frequency domain was also calculated between the aHipp and the pHipp for the incorrect trials. However, only a few scattered PLV points survived significance testing during the maintenance of incorrect trials (Fig. 5C). As a result, we did not conduct further analyses on the PLVs for the incorrect trials as we did for the correct trials.

Figure 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 5.

Power and phase synchrony for the incorrect trials. A, Average task-induced power across all subjects during maintenance for the incorrect trials, grouped in the anterior (left column) and posterior hippocampus (right column). Warmer colors denote task-increased power from the baseline and colder colors denote task-decreased power. Task-induced effects on both the aHipp and the pHipp were not sustained during maintenance in the low-frequency range (2–20 Hz; z < 1.96) for the incorrect trials. B, Spectral z-scored power within the anterior (magenta) and posterior hippocampus (cyan) across all subjects (±SEM shown as shading around the mean trace) for the incorrect trials. The spectral power was not significant (z < 1.96) during maintenance, and no significant difference (p > 0.05, clustered-based permutation test) was found between the two subregions across the frequency band. C, Significant PLVs above the thresholds (p < 0.05. permutation test) across the time–frequency domain were present between the aHipp and the pHipp for the incorrect trials. Only a few scattered points survived testing for incorrect trials. Therefore, the PLV of incorrect trials was not further analyzed.

Volume conduction and recording montage

To rule out effects of volume conduction, we repeated our analysis of connectivity and Granger causality after rereferencing our recordings to a bipolar scheme (for a similar approach, see Kota et al., 2020). First, we found sustained phase synchronization (p < 0.05, clustered-based permutation test) between aHipp and pHipp in the low-frequency band throughout the entire maintenance period (Fig. 6A). Then, we found higher GC values from the pHipp to aHipp than those in the opposite direction, with a trend toward statistical significance (mixed-effect model, p = 0.092) for the correct trials (Fig. 6B) and no significant difference between the two directions (mixed-effect model, p > 0.05) for the incorrect trials. These findings are consistent with what we observed by using a common reference. In addition, the following findings speak against volume conduction: (1) there was a strong frequency dependence of the intrahippocampal PLVs (Fig. 2A) and those to scalp EEGs (Boran et al., 2019), while the transfer of signal through tissue by volume conduction is independent of frequency in the range of the iEEG studied here (Nunez and Cutillo, 1995; Miceli et al., 2017); (2) there was a strong task dependence of GC values (Figs. 2D, 4A); and (3) GC is highly anisotropic (Figs. 2D, 4A). In conclusion, while we cannot rule out a contribution of volume conduction in our data, it seems to be small compared with the effects in the PLVs and GC, which were induced by processing the cognitive task.

Figure 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 6.

PLV and GC analyses with a bipolar reference. A, Average phase synchronization (z-scored PLV) across electrode pairs for the correct trials, with a bipolar reference. Elevated PLV in the low-frequency band was sustained throughout the entire maintenance period. B, Granger causality between hippocampal subregions from both directions by using a bipolar reference for the correct trials. #p < 0.1. GC values from the pHipp to aHipp are higher than those in the opposite direction with a trend toward statistical significance (mixed-effect model, p = 0.092).

Recordings in the epileptic seizure onset zone

Our analysis included recordings from electrodes in the seizure onset zone (SOZ). To assess the possible effects of the SOZ on our findings, we compared the hippocampal electrodes within and outside the SOZ. For the power analysis, the task-induced power was not significantly different between the non-SOZ electrodes and the SOZ electrodes (p = 0.52). For the connectivity analysis, we separated the electrode pairs into SOZ (at least one SOZ electrode) and non-SOZ (non-SOZ to non-SOZ electrode) pairs. There was no significant difference in the task-related connectivity measures between the SOZ and non-SOZ electrode pairs (PLV analysis, p = 0.27; PAC analysis: aHipp–pHipp, p = 0.62; pHipp-aHipp, p = 0.66; Granger causality analysis: aHipp to pHipp, p = 0.14; pHipp to aHipp, p = 0.40). Our analysis implicated residual function in the epileptogenic hippocampus. This is consistent with our previous finding of single neurons in the SOZ that responded to this task (Boran et al., 2019) and more general findings suggesting normal physiological responses to cognitive stimuli in the epileptic sites (Liu and Parvizi, 2019).

Discussion

Together, our results indicated that both subregions of the hippocampus were critical for WM processing; they interacted at a specific frequency band and with a directionality pattern to support the maintenance of the WM information. Specifically, WM maintenance depends on coordinated neural oscillations between the aHipp and the pHipp in the theta/alpha band, with a pHipp-driven information transfer in this band, and a higher PAC between the theta/alpha phase of the pHipp entraining the gamma amplitude of the aHipp. The directionality bias in information transfer within the hippocampus occurred exclusively during the correct WM trials.

The task-induced power revealed sustained low-frequency oscillations (including theta and alpha bands) within both the aHipp and the pHipp during the delay interval in which information was to be remembered. The theta and alpha bands are well known frequency bands that play important roles during WM maintenance. The theta rhythm is thought to underlie the organization of WM items (Roux and Uhlhaas, 2014), and theta power increases in the hippocampus reflect the storage component of WM (Brzezicka et al., 2019). Alpha oscillations were suggested to protect working memory maintenance against anticipated distracters (Bonnefond and Jensen, 2012) and increased during WM maintenance (Boran et al., 2019). Our results echo these findings and suggest that increases in theta/alpha power could support hippocampal WM processing by organizing neural activity (theta) and selectively filtering task-relevant information to avoid distracting information (alpha).

We further observed oscillatory phase synchrony between the aHipp and the pHipp predominantly in the theta/alpha range (3–12 Hz) during the retention interval. Only a few studies have studied oscillatory communication along the long axis of the hippocampus by simultaneously recording both areas in awake unrestrained rodents. Specifically, rodent studies reported increased connectivity between the dorsoventral hippocampus during a WM task (Schmidt et al., 2013). Lesion studies showed a disruption in WM accompanied by reduced theta phase coherence in the dorsoventral hippocampus, and pharmacological treatment resulted in both restored WM task performance and reversed abnormal hippocampal dorsoventral connectivity (Cardoso-Cruz et al., 2014). Our results are thus consistent with the animal literature suggesting that the contribution of synchrony between hippocampal subregions along the long axis to WM processing and extended these findings to humans.

Next, we found that the theta/alpha phase of the pHipp entrained the gamma activity of the aHipp and that theta/alpha-driven information flow occurs from the pHipp to the aHipp during the retention interval with subsequent successful WM responses. Intriguingly, this directional pHipp-to-aHipp information flow was only observed in correct trials and was absent in error trials. These results together suggest that the posterior to anterior influence in the theta/alpha band contributes to the successful maintenance of WM. These interregional information transfer dynamics are consistent with the anatomic circuitry showing projections from the dorsal to the more ventral part of the hippocampus (Amaral and Witter, 1989; Ishizuka et al., 1990). The longitudinal axis may be a locus for interaction between various pieces of WM information with longitudinal interneurons orchestrating the information flow (Harland et al., 2018). Human fMRI studies also suggested hierarchical representations of events along the hippocampal long axis with the pHipp representing detail-level events and the aHipp representing overlapping higher-level, multievent representations (Collin et al., 2015; Schlichting et al., 2015). In our research, the participants needed to remember individual letters; therefore, we speculate that there was a greater requirement for detailed representation, so the pHipp led the information transfer to the aHipp.

Together, our findings suggest a hierarchical information flow from the pHipp to the aHipp in theta/alpha oscillations to successfully maintain WM information. By using human iEEG recording and precise hippocampal subregion segmentation, we demonstrated that working memory maintenance depends on the activity of both the aHipp and the pHipp, which is in line with human fMRI findings. In addition, we provide evidence for coordinated neural oscillations between them with a unidirectional influence from the pHipp to the aHipp via theta/alpha oscillations. This echoes findings from rodents in which lesions of the dorsal Hipp, but not the ventral Hipp, led to deficits in working memory (Pothuizen et al., 2004), as well as rodent studies suggesting preferential involvement of cognitive function in the dorsal Hipp and of emotion in the ventral Hipp (Fanselow and Dong, 2010). Together, our work establishes a bridge between cellular-level work in animals and the extensive literature using noninvasive imaging in humans on the functional specialization of subregions within the hippocampus. The observed directional functional connectivity between the dorsal and ventral subregions is well in line with proposed models of hippocampal long-axis organization within discrete domains (Strange et al., 2014).

There are several limitations of the current study. First, only a verbal WM task was used here. This has the advantage of excluding the possibility of explaining the pHipp-dominated information transfer during WM as a preferential involvement of the pHipp in spatial processing. However, the same strength of the current study limits the generalization of our results to WM processing apart from other information, such as space and time. Second, we did not further map the aHipp or pHipp into cytoarchitectural subfields. In addition to the organizational patterns along the long axis, the hippocampus is composed of cytoarchitecturally different subfields (CA1–4, dentate gyrus, subiculum; Amunts et al., 2005). However, the small number of electrode contacts in the current study (aHipp for each subject: mean, 3.2 ± 1.8; range, 0–6; pHipp for each subject: mean, 6.7 ± 2.2; range, 1–8) limited the investigation of subfield-specific oscillation-based communication. Future studies with a larger sample size are thus needed to consider the information flow between cytoarchitectural subfields.

In summary, our results linked the intrahippocampal theta/alpha band interaction to WM maintenance and showed that directional information flow occurred from the pHipp to the aHipp during the maintenance interval when the subsequent WM responses were correct. These results add to the growing body of literature about the functional organization along the longitudinal axis of the hippocampus. More broadly, this work allows a detailed mechanistic and more direct understanding of WM.

Footnotes

  • This work received support from the following sources: National Key Research and Development Program of China (2017YFA0105203); the Chinese Academy of Sciences Science and Technology Service Network Initiative (Grant KFJ-STS-ZDTP-078); National Natural Science Foundation of China (Grants 31300934 and 31620103905); Science Frontier Program of the Chinese Academy of Sciences (Grant QYZDJ-SSW-SMC019); the National Key R&D Program of China (Grant 2017YFA0105203); Open Research Fund of the State Key Laboratory of Cognitive Neuroscience and Learning (Grant CNLYB2004); Open Research Fund of the CAS Key Laboratory of Behavioral Science, the Institute of Psychology; and the Swiss National Science Foundation (Grant SNSF 176222 to J.S.). We thank Prof. Yuanye Ma from the Key Laboratory of External Drug Delivery System and Preparation Technology in the University of Yunnan Province, and the School of Chinese Materia Medica, Yunnan University of Chinese Medicine for suggestions. We also thank Dr. Rhoda E. Perozzi and Dr. Edmund F. Perozzi for English and content editing assistance.

  • The authors declare no competing financial interests.

  • Correspondence should be addressed to Tianzi Jiang at jiangtz{at}nlpr.ia.ac.cn or Johannes Sarnthein at johannes.sarnthein{at}usz.ch

SfN exclusive license.

References

  1. ↵
    1. Amaral DG,
    2. Witter MP
    (1989) The three-dimensional organization of the hippocampal formation: a review of anatomical data. Neuroscience 31:571–591. doi:10.1016/0306-4522(89)90424-7 pmid:2687721
    OpenUrlCrossRefPubMed
  2. ↵
    1. Amunts K,
    2. Kedo O,
    3. Kindler M,
    4. Pieperhoff P,
    5. Mohlberg H,
    6. Shah NJ,
    7. Habel U,
    8. Schneider F,
    9. Zilles K
    (2005) Cytoarchitectonic mapping of the human amygdala, hippocampal region and entorhinal cortex: intersubject variability and probability maps. Anat Embryol (Berl) 210:343–352. doi:10.1007/s00429-005-0025-5 pmid:16208455
    OpenUrlCrossRefPubMed
  3. ↵
    1. Assaneo MF,
    2. Ripollés P,
    3. Orpella J,
    4. Lin WM,
    5. de Diego-Balaguer R,
    6. Poeppel D
    (2019) Spontaneous synchronization to speech reveals neural mechanisms facilitating language learning. Nat Neurosci 22:627–632. doi:10.1038/s41593-019-0353-z pmid:30833700
    OpenUrlCrossRefPubMed
  4. ↵
    1. Axmacher N,
    2. Henseler MM,
    3. Jensen O,
    4. Weinreich I,
    5. Elger CE,
    6. Fell J
    (2010) Cross-frequency coupling supports multi-item working memory in the human hippocampus. Proc Natl Acad Sci U S A 107:3228–3233. doi:10.1073/pnas.0911531107 pmid:20133762
    OpenUrlAbstract/FREE Full Text
  5. ↵
    1. Barnett L,
    2. Seth AK
    (2014) The MVGC multivariate Granger causality toolbox: a new approach to Granger-causal inference. J Neurosci Methods 223:50–68. doi:10.1016/j.jneumeth.2013.10.018 pmid:24200508
    OpenUrlCrossRefPubMed
  6. ↵
    1. Basti A,
    2. Pizzella V,
    3. Chella F,
    4. Romani GL,
    5. Nolte G,
    6. Marzetti L
    (2018) Disclosing large-scale directed functional connections in MEG with the multivariate phase slope index. Neuroimage 175:161–175. doi:10.1016/j.neuroimage.2018.03.004 pmid:29524622
    OpenUrlCrossRefPubMed
  7. ↵
    1. Bonnefond M,
    2. Jensen O
    (2012) Alpha oscillations serve to protect working memory maintenance against anticipated distracters. Curr Biol 22:1969–1974. doi:10.1016/j.cub.2012.08.029 pmid:23041197
    OpenUrlCrossRefPubMed
  8. ↵
    1. Boran E,
    2. Fedele T,
    3. Klaver P,
    4. Hilfiker P,
    5. Stieglitz L,
    6. Grunwald T,
    7. Sarnthein J
    (2019) Persistent hippocampal neural firing and hippocampal-cortical coupling predict verbal working memory load. Sci Adv 5:eaav3687. doi:10.1126/sciadv.aav3687 pmid:30944858
    OpenUrlFREE Full Text
  9. ↵
    1. Boran E,
    2. Fedele T,
    3. Steiner A,
    4. Hilfiker P,
    5. Stieglitz L,
    6. Grunwald T,
    7. Sarnthein J
    (2020) Dataset of human medial temporal lobe neurons, scalp and intracranial EEG during a verbal working memory task. Sci Data 7:30. doi:10.1038/s41597-020-0364-3 pmid:31964868
    OpenUrlCrossRefPubMed
  10. ↵
    1. Boran E,
    2. Stieglitz L,
    3. Sarnthein J
    (2021) Epileptic high-frequency oscillations in intracranial EEG are not confounded by cognitive tasks. Front Hum Neurosci 15:613125.
    OpenUrl
  11. ↵
    1. Brincat SL,
    2. Miller EK
    (2015) Frequency-specific hippocampal-prefrontal interactions during associative learning. Nat Neurosci 18:576–581. doi:10.1038/nn.3954 pmid:25706471
    OpenUrlCrossRefPubMed
  12. ↵
    1. Brzezicka A,
    2. Kamiński J,
    3. Reed CM,
    4. Chung JM,
    5. Mamelak AN,
    6. Rutishauser U
    (2019) Working memory load-related theta power decreases in dorsolateral prefrontal cortex predict individual differences in performance. J Cogn Neurosci 31:1290–1307. doi:10.1162/jocn_a_01417 pmid:31037988
    OpenUrlCrossRefPubMed
  13. ↵
    1. Canolty RT,
    2. Edwards E,
    3. Dalal SS,
    4. Soltani M,
    5. Nagarajan SS,
    6. Kirsch HE,
    7. Berger MS,
    8. Barbaro NM,
    9. Knight RT
    (2006) High gamma power is phase-locked to theta oscillations in human neocortex. Science 313:1626–1628. doi:10.1126/science.1128115 pmid:16973878
    OpenUrlAbstract/FREE Full Text
  14. ↵
    1. Cardoso-Cruz H,
    2. Dourado M,
    3. Monteiro C,
    4. Matos MR,
    5. Galhardo V
    (2014) Activation of dopaminergic D2/D3 receptors modulates dorsoventral connectivity in the hippocampus and reverses the impairment of working memory after nerve injury. J Neurosci 34:5861–5873. doi:10.1523/JNEUROSCI.0021-14.2014 pmid:24760846
    OpenUrlAbstract/FREE Full Text
  15. ↵
    1. Collin SH,
    2. Milivojevic B,
    3. Doeller CF
    (2015) Memory hierarchies map onto the hippocampal long axis in humans. Nat Neurosci 18:1562–1564. doi:10.1038/nn.4138 pmid:26479587
    OpenUrlCrossRefPubMed
  16. ↵
    1. Delorme A,
    2. Makeig S
    (2004) EEGLAB: an open source toolbox for analysis of single-trial EEG dynamics including independent component analysis. J Neurosci Methods 134:9–21. doi:10.1016/j.jneumeth.2003.10.009 pmid:15102499
    OpenUrlCrossRefPubMed
  17. ↵
    1. D'Esposito M
    (2007) From cognitive to neural models of working memory. Philos Trans R Soc Lond B Biol Sci 362:761–772. doi:10.1098/rstb.2007.2086 pmid:17400538
    OpenUrlCrossRefPubMed
  18. ↵
    1. Fan L,
    2. Li H,
    3. Zhuo J,
    4. Zhang Y,
    5. Wang J,
    6. Chen L,
    7. Yang Z,
    8. Chu C,
    9. Xie S,
    10. Laird AR,
    11. Fox PT,
    12. Eickhoff SB,
    13. Yu C,
    14. Jiang T
    (2016) The human brainnetome atlas: a new brain atlas based on connectional architecture. Cereb Cortex 26:3508–3526. doi:10.1093/cercor/bhw157 pmid:27230218
    OpenUrlCrossRefPubMed
  19. ↵
    1. Fanselow MS,
    2. Dong HW
    (2010) Are the dorsal and ventral hippocampus functionally distinct structures? Neuron 65:7–19. doi:10.1016/j.neuron.2009.11.031
    OpenUrlCrossRefPubMed
  20. ↵
    1. Grady CL
    (2020) Meta-analytic and functional connectivity evidence from functional magnetic resonance imaging for an anterior to posterior gradient of function along the hippocampal axis. Hippocampus 30:456–471. doi:10.1002/hipo.23164 pmid:31589003
    OpenUrlCrossRefPubMed
  21. ↵
    1. Harland B,
    2. Contreras M,
    3. Fellous J-M
    (2018) A role for the longitudinal axis of the hippocampus in multiscale representations of large and complex spatial environments and mnemonic hierarchies. In: The hippocampus-plasticity and functions, p 67. London: IntechOpen.
  22. ↵
    1. Ishizuka N,
    2. Weber J,
    3. Amaral DG
    (1990) Organization of intrahippocampal projections originating from CA3 pyramidal cells in the rat. J Comp Neurol 295:580–623. doi:10.1002/cne.902950407 pmid:2358523
    OpenUrlCrossRefPubMed
  23. ↵
    1. Johnson EL,
    2. Adams JN,
    3. Solbakk AK,
    4. Endestad T,
    5. Larsson PG,
    6. Ivanovic J,
    7. Meling TR,
    8. Lin JJ,
    9. Knight RT
    (2018) Dynamic frontotemporal systems process space and time in working memory. PLoS Biol 16:e2004274. doi:10.1371/journal.pbio.2004274 pmid:29601574
    OpenUrlCrossRefPubMed
  24. ↵
    1. Johnson EL,
    2. Kam JWY,
    3. Tzovara A,
    4. Knight RT
    (2020) Insights into human cognition from intracranial EEG: a review of audition, memory, internal cognition, and causality. J Neural Eng 17:051001. doi:10.1088/1741-2552/abb7a5 pmid:32916678
    OpenUrlCrossRefPubMed
  25. ↵
    1. Kamiński J,
    2. Sullivan S,
    3. Chung JM,
    4. Ross IB,
    5. Mamelak AN,
    6. Rutishauser U
    (2017) Persistently active neurons in human medial frontal and medial temporal lobe support working memory. Nat Neurosci 20:590–601. doi:10.1038/nn.4509 pmid:28218914
    OpenUrlCrossRefPubMed
  26. ↵
    1. Kamiński J,
    2. Brzezicka A,
    3. Mamelak AN,
    4. Rutishauser U
    (2020) Combined phase-rate coding by persistently active neurons as a mechanism for maintaining multiple items in working memory in humans. Neuron 106:256–264.e3. doi:10.1016/j.neuron.2020.01.032 pmid:32084331
    OpenUrlCrossRefPubMed
  27. ↵
    1. Kornblith S,
    2. Quian Quiroga R,
    3. Koch C,
    4. Fried I,
    5. Mormann F
    (2017) Persistent single-neuron activity during working memory in the human medial temporal lobe. Curr Biol 27:1026–1032. doi:10.1016/j.cub.2017.02.013 pmid:28318972
    OpenUrlCrossRefPubMed
  28. ↵
    1. Kota S,
    2. Rugg MD,
    3. Lega BC
    (2020) Hippocampal theta oscillations support successful associative memory formation. J Neurosci 40:9507–9518. doi:10.1523/JNEUROSCI.0767-20.2020 pmid:33158958
    OpenUrlAbstract/FREE Full Text
  29. ↵
    1. Lee J-J,
    2. Kim HJ,
    3. Čeko M,
    4. Park B-Y,
    5. Lee SA,
    6. Park H,
    7. Roy M,
    8. Kim S-G,
    9. Wager TD,
    10. Woo C-W
    (2021) A neuroimaging biomarker for sustained experimental and clinical pain. Nat Med 27:174–182. doi:10.1038/s41591-020-1142-7 pmid:33398159
    OpenUrlCrossRefPubMed
  30. ↵
    1. Lin JJ,
    2. Rugg MD,
    3. Das S,
    4. Stein J,
    5. Rizzuto DS,
    6. Kahana MJ,
    7. Lega BC
    (2017) Theta band power increases in the posterior hippocampus predict successful episodic memory encoding in humans. Hippocampus 27:1040–1053. doi:10.1002/hipo.22751 pmid:28608960
    OpenUrlCrossRefPubMed
  31. ↵
    1. Lin JJ,
    2. Umbach G,
    3. Rugg MD,
    4. Lega B
    (2019) Gamma oscillations during episodic memory processing provide evidence for functional specialization in the longitudinal axis of the human hippocampus. Hippocampus 29:68–72. doi:10.1002/hipo.23016 pmid:30394594
    OpenUrlCrossRefPubMed
  32. ↵
    1. Liu S,
    2. Parvizi J
    (2019) Cognitive refractory state caused by spontaneous epileptic high-frequency oscillations in the human brain. Sci Transl Med 11:eaax7830.
    OpenUrlAbstract/FREE Full Text
  33. ↵
    1. Maris E,
    2. Oostenveld R
    (2007) Nonparametric statistical testing of EEG- and MEG-data. J Neurosci Methods 164:177–190. doi:10.1016/j.jneumeth.2007.03.024 pmid:17517438
    OpenUrlCrossRefPubMed
  34. ↵
    1. Martínez-Cancino R,
    2. Heng J,
    3. Delorme A,
    4. Kreutz-Delgado K,
    5. Sotero RC,
    6. Makeig S
    (2019) Measuring transient phase-amplitude coupling using local mutual information. Neuroimage 185:361–378. doi:10.1016/j.neuroimage.2018.10.034 pmid:30342235
    OpenUrlCrossRefPubMed
  35. ↵
    1. Miceli S,
    2. Ness TV,
    3. Einevoll GT,
    4. Schubert D
    (2017) Impedance spectrum in cortical tissue: implications for propagation of LFP signals on the microscopic level. eNeuro 4:ENEURO.0291-16.2016. doi:10.1523/ENEURO.0291-16.2016
    OpenUrlAbstract/FREE Full Text
  36. ↵
    1. Miller EK,
    2. Lundqvist M,
    3. Bastos AM
    (2018) Working memory 2.0. Neuron 100:463–475. doi:10.1016/j.neuron.2018.09.023 pmid:30359609
    OpenUrlCrossRefPubMed
  37. ↵
    1. Nunez PL,
    2. Cutillo BA
    (1995) Neocortical dynamics and human EEG rhythms. Oxford: Oxford UP.
  38. ↵
    1. Oostenveld R,
    2. Fries P,
    3. Maris E,
    4. Schoffelen JM
    (2011) FieldTrip: open source software for advanced analysis of MEG, EEG, and invasive electrophysiological data. Comput Intell Neurosci 2011:156869. doi:10.1155/2011/156869 pmid:21253357
    OpenUrlCrossRefPubMed
  39. ↵
    1. Poppenk J,
    2. Evensmoen HR,
    3. Moscovitch M,
    4. Nadel L
    (2013) Long-axis specialization of the human hippocampus. Trends Cogn Sci 17:230–240. doi:10.1016/j.tics.2013.03.005 pmid:23597720
    OpenUrlCrossRefPubMed
  40. ↵
    1. Pothuizen HH,
    2. Zhang WN,
    3. Jongen-Rêlo AL,
    4. Feldon J,
    5. Yee BK
    (2004) Dissociation of function between the dorsal and the ventral hippocampus in spatial learning abilities of the rat: a within-subject, within-task comparison of reference and working spatial memory. Eur J Neurosci 19:705–712. doi:10.1111/j.0953-816x.2004.03170.x pmid:14984421
    OpenUrlCrossRefPubMed
  41. ↵
    1. Roux F,
    2. Uhlhaas PJ
    (2014) Working memory and neural oscillations: α-γ versus θ-γ codes for distinct WM information? Trends Cogn Sci 18:16–25. doi:10.1016/j.tics.2013.10.010 pmid:24268290
    OpenUrlCrossRefPubMed
  42. ↵
    1. Schlichting ML,
    2. Mumford JA,
    3. Preston AR
    (2015) Learning-related representational changes reveal dissociable integration and separation signatures in the hippocampus and prefrontal cortex. Nat Commun 6:8151. doi:10.1038/ncomms9151 pmid:26303198
    OpenUrlCrossRefPubMed
  43. ↵
    1. Schmidt B,
    2. Hinman JR,
    3. Jacobson TK,
    4. Szkudlarek E,
    5. Argraves M,
    6. Escabí MA,
    7. Markus EJ
    (2013) Dissociation between dorsal and ventral hippocampal theta oscillations during decision-making. J Neurosci 33:6212–6224. doi:10.1523/JNEUROSCI.2915-12.2013 pmid:23554502
    OpenUrlAbstract/FREE Full Text
  44. ↵
    1. Solomon EA,
    2. Stein JM,
    3. Das S,
    4. Gorniak R,
    5. Sperling MR,
    6. Worrell G,
    7. Inman CS,
    8. Tan RJ,
    9. Jobst BC,
    10. Rizzuto DS,
    11. Kahana MJ
    (2019) Dynamic theta networks in the human medial temporal lobe support episodic memory. Curr Biol 29:1100–1111.e4. doi:10.1016/j.cub.2019.02.020 pmid:30905609
    OpenUrlCrossRefPubMed
  45. ↵
    1. Stolk A,
    2. Griffin S,
    3. van der Meij R,
    4. Dewar C,
    5. Saez I,
    6. Lin JJ,
    7. Piantoni G,
    8. Schoffelen JM,
    9. Knight RT,
    10. Oostenveld R
    (2018) Integrated analysis of anatomical and electrophysiological human intracranial data. Nat Protoc 13:1699–1723. doi:10.1038/s41596-018-0009-6 pmid:29988107
    OpenUrlCrossRefPubMed
  46. ↵
    1. Strange BA,
    2. Witter MP,
    3. Lein ES,
    4. Moser EI
    (2014) Functional organization of the hippocampal longitudinal axis. Nat Rev Neurosci 15:655–669. doi:10.1038/nrn3785 pmid:25234264
    OpenUrlCrossRefPubMed
  47. ↵
    1. Tort AB,
    2. Kramer MA,
    3. Thorn C,
    4. Gibson DJ,
    5. Kubota Y,
    6. Graybiel AM,
    7. Kopell NJ
    (2008) Dynamic cross-frequency couplings of local field potential oscillations in rat striatum and hippocampus during performance of a T-maze task. Proc Natl Acad Sci U S A 105:20517–20522. doi:10.1073/pnas.0810524105 pmid:19074268
    OpenUrlAbstract/FREE Full Text
  48. ↵
    1. Vogel JW,
    2. La Joie R,
    3. Grothe MJ,
    4. Diaz-Papkovich A,
    5. Doyle A,
    6. Vachon-Presseau E,
    7. Lepage C,
    8. Vos de Wael R,
    9. Thomas RA,
    10. Iturria-Medina Y,
    11. Bernhardt B,
    12. Rabinovici GD,
    13. Evans AC
    (2020) A molecular gradient along the longitudinal axis of the human hippocampus informs large-scale behavioral systems. Nat Commun 11:960. doi:10.1038/s41467-020-14518-3 pmid:32075960
    OpenUrlCrossRefPubMed
  49. ↵
    1. Vos de Wael R,
    2. Larivière S,
    3. Caldairou B,
    4. Hong SJ,
    5. Margulies DS,
    6. Jefferies E,
    7. Bernasconi A,
    8. Smallwood J,
    9. Bernasconi N,
    10. Bernhardt BC
    (2018) Anatomical and microstructural determinants of hippocampal subfield functional connectome embedding. Proc Natl Acad Sci U S A 115:10154–10159. doi:10.1073/pnas.1803667115 pmid:30249658
    OpenUrlAbstract/FREE Full Text
  50. ↵
    1. Wang DX,
    2. Schmitt K,
    3. Seger S,
    4. Davila CE,
    5. Lega BC
    (2021) Cross-regional phase amplitude coupling supports the encoding of episodic memories. Hippocampus 31:481–492. doi:10.1002/hipo.23309 pmid:33544408
    OpenUrlCrossRefPubMed
Back to top

In this issue

The Journal of Neuroscience: 42 (3)
Journal of Neuroscience
Vol. 42, Issue 3
19 Jan 2022
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Index by author
  • Ed Board (PDF)
Email

Thank you for sharing this Journal of Neuroscience article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Anterior–Posterior Hippocampal Dynamics Support Working Memory Processing
(Your Name) has forwarded a page to you from Journal of Neuroscience
(Your Name) thought you would be interested in this article in Journal of Neuroscience.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Print
View Full Page PDF
Citation Tools
Anterior–Posterior Hippocampal Dynamics Support Working Memory Processing
Jin Li (李瑾), Dan Cao (曹丹), Vasileios Dimakopoulos, Weiyang Shi (时维 阳), Shan Yu (余山), Lingzhong Fan (樊令 仲), Lennart Stieglitz, Lukas Imbach, Johannes Sarnthein, Tianzi Jiang (蒋田 仔)
Journal of Neuroscience 19 January 2022, 42 (3) 443-453; DOI: 10.1523/JNEUROSCI.1287-21.2021

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Respond to this article
Request Permissions
Share
Anterior–Posterior Hippocampal Dynamics Support Working Memory Processing
Jin Li (李瑾), Dan Cao (曹丹), Vasileios Dimakopoulos, Weiyang Shi (时维 阳), Shan Yu (余山), Lingzhong Fan (樊令 仲), Lennart Stieglitz, Lukas Imbach, Johannes Sarnthein, Tianzi Jiang (蒋田 仔)
Journal of Neuroscience 19 January 2022, 42 (3) 443-453; DOI: 10.1523/JNEUROSCI.1287-21.2021
del.icio.us logo Digg logo Reddit logo Twitter logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Introduction
    • Materials and Methods
    • Results
    • Discussion
    • Footnotes
    • References
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF

Keywords

  • hippocampus
  • intracranial EEG
  • longitudinal axis
  • working memory

Responses to this article

Respond to this article

Jump to comment:

No eLetters have been published for this article.

Related Articles

Cited By...

More in this TOC Section

Research Articles

  • Distinct features of interictal activity predict seizure localization and burden in a mouse model of childhood epilepsy
  • Neuromedin B-expressing neurons in the retrotrapezoid nucleus regulate respiratory homeostasis and promote stable breathing in adult mice
  • A visual pathway into central complex for high frequency motion-defined bars in Drosophila
Show more Research Articles

Behavioral/Cognitive

  • Posterior Parietal Cortex Plays a Causal Role in Abstract Memory-Based Visual Categorical Decisions
  • The Differential Weights of Motivational and Task Performance Measures on Medial and Lateral Frontal Neural Activity
  • Prior Movement of One Arm Facilitates Motor Adaptation in the Other
Show more Behavioral/Cognitive
  • Home
  • Alerts
  • Visit Society for Neuroscience on Facebook
  • Follow Society for Neuroscience on Twitter
  • Follow Society for Neuroscience on LinkedIn
  • Visit Society for Neuroscience on Youtube
  • Follow our RSS feeds

Content

  • Early Release
  • Current Issue
  • Issue Archive
  • Collections

Information

  • For Authors
  • For Advertisers
  • For the Media
  • For Subscribers

About

  • About the Journal
  • Editorial Board
  • Privacy Policy
  • Contact
(JNeurosci logo)
(SfN logo)

Copyright © 2023 by the Society for Neuroscience.
JNeurosci Online ISSN: 1529-2401

The ideas and opinions expressed in JNeurosci do not necessarily reflect those of SfN or the JNeurosci Editorial Board. Publication of an advertisement or other product mention in JNeurosci should not be construed as an endorsement of the manufacturer’s claims. SfN does not assume any responsibility for any injury and/or damage to persons or property arising from or related to any use of any material contained in JNeurosci.