Abstract
Introduction:
Spatial working memory is a kind of short-term memory that allows temporarily storing and manipulating spatial information. Evidence suggests that spatial working memory is processed through three distinctive phases: Encoding, maintenance, and retrieval. Though the medial prefrontal cortex (mPFC) and mediodorsal thalamus (MD) are involved in memory retrieval, how the functional interactions and information transfer between mPFC and MD remains largely unclear.
Methods:
We recorded local field potentials (LFPs) from mPFC and MD while mice performed a spatial working memory task in T-maze. The temporal dynamics of functional interactions and bidirectional information flow between mPFC and MD was quantitatively assessed by using directed transfer function.
Results:
Our results showed a significantly elevated information flow from mPFC to MD, varied in time and frequency (theta in particular), accompanying successful memory retrieval.
Discussion:
Elevated theta information flow, a feature that was absent on error trials, indicates an important role of the directional information transfer from mPFC to MD for memory retrieval.
Introduction
Working memory is defined as a brain system that provides temporary storage and manipulation of the information necessary for complex cognitive tasks (, ; ). Spatial working memory, especially, refers to the temporary online maintenance and/or manipulation of spatial information to be used toward a specific goal (; ).
Cognitive neuroscientists have employed empirical and theoretical approaches to study the components of spatial working memory in attempts to localize and characterize the neural implementation. The process of spatial working memory is recognized typically divided into three phases: Encoding, maintenance, and retrieval (; ; ; ). Specially, encoding involves the loading of information into working memory, while maintenance encompasses the active rehearsal of that information for a brief period of time. Finally, retrieval refers to the recall and application of the information to achieve a cognitive goal, which is also bound up with the updating of memory information in working memory and evaluation or monitoring of this information.
Extensive evidence has demonstrated the medial prefrontal cortex (mPFC) and mediodorsal thalamus (MD) is involved in working memory retrieval (; ; ; ; ). The reciprocal MD-mPFC activity is recognized to be crucial for successful execution of spatial working memory task (, ; ). Though anatomical evidence has shown the direct projections between MD and ventral-medial PFC (; ), the question of how the information is parsed and dynamically transferred between the two areas to support working memory remains largely unclear.
Recently, brain networks have attracted increasing attention due to the potential to better characterize brain dynamics in cognition (; ; ; ; ). The methods with brain network analysis can be useful tools to estimate the directional functional interactions among brain regions and explore the mechanism for integration in the brain (; ; ). So far, brain network analysis method has been used to explore the mechanism underlying working memory in human and revealed the altered connectivity in functional networks during working memory performance persists following conclusion of that performance (; ; ). Previous studies have also shown the enhanced directional information transfer from hippocampus to prefrontal cortex is required for successful execution of a delayed non-match-to-place (DNMTP) task in rat (). Therefore, the network analysis approaches could offer crucial insight into the nature of functional interactions and information transfer in the mPFC-MD network correlates of spatial working memory performance.
In the present study, we recorded the LFPs from mPFC and MD of mice as they performed a spatial working memory task. Then the bidirectional information flow between the mPFC and MD were quantitatively evaluated during the task epochs. This study is expected to reveal the directional information transfer between mPFC and MD underlying spatial working memory retrieval.
Materials and methods
All experimental procedures were conducted in accordance with the Guide for Care and Use of Laboratory Animals and approved by the Tianjin Medical University Animal Care and Use Committee (license number: TMUaMEC 2021060).
Subjects
Male adult C57BL/6 mice (aged 10–12 weeks) were used for the current study and provided by the Experimental Animal Center of Tianjin Medical University (Tianjin, China). Mice were group-housed (4–5 per cage) under a 12-h light-dark cycle. Except when food-restricted for the purpose of behavioral training and testing, all mice were given ad libitum access to food and water.
Delayed non-match-to-place task
Animals were trained on a DNMTP task in T-maze, a typical spatial working memory task for rodents, as described previously (; ). The schema of a single trial of the DNMTP task can be found in Figure 1A. A total of 3 days before the training, mice were gradually food restricted to maintain 85% of their body weight. Mice were then given 2 days of habituation to the maze, which consisted of 10 min free exploration and foraging with all doors open. On the subsequent 2 days mice underwent behavioral shaping consisting of 10 min of running to baited goal arms in alternating directions and back to the start box. Mice then underwent training on the T-maze task until criterion performance (consisting of 70% correct trials on 2 out of 3 consecutive days) was achieved (). In this experiment, intra-trial delay was 10 s and inter-trial delay was 20 s. The performances of the mice in the training sessions were measured by the accuracy of responses. The behavioral results showed that the mice (n = 8) learned the task gradually, improving from 58.8 ± 5.4% to 82.5 ± 3.9% correct choice over 6 days of training (Figure 1B).
FIGURE 1
Electrodes array construction and implantation
Multi-electrode array for extracellular recording was constructed with polyimide insulated, tungsten wire (diameter: 13.97 μm, impedance <1 MΩ). All electrodes were referenced to a common electrode and common ground. For simultaneous recording of mPFC and MD, chamber was designed for accurate targeting to these two brain regions. For implantation, mice were deeply anesthetized with analgesics (Sodium pentobarbital, 10 mg/mL, 10 mg/kg i.p.) before surgery. Dura was carefully removed for the insertion of electrodes. According to the mouse brain atlas in stereotaxic coordinates, the stereotrodes were implanted in the mPFC (1.75 mm anterior, 0.4 mm lateral, and 1.8 mm ventral) and MD (1.2 mm posterior, 0.35 mm lateral, and 3.2 mm ventral). Finally, the electrodes were secured with dental cement. The recording sites of electrodes were further verified with histology after experiments.
Neurophysiological recording
After electrodes implantation, mice were recovered for 1 week before neurophysiological recording. Recordings were made from mPFC and MD while the mice again performed the spatial working memory task (Figure 1C). Extracellular signals were recorded with the using the Cerebus Acquisition System (Blackrock Microsystems Inc., UT, USA), and the performance of the mice was obtained simultaneously by the overhead camera. For LFP recording, the signals were filtered with a pass-band of 0.5–120 Hz, further amplified and digitized at 1 kHz. Then polynomial fitting was used to remove the baseline drifts from LFPs and a 50-Hz notch filter was used to remove the powerline artifact. The behavioral data was also analyzed and the mouse’s initial position was defined as the point a/d, and the turn began was denoted as the point b/e, and the reward position was denoted as the point c/fduring sample/choice runs (marked in Figure 1A). For comparison, the neuronal activity in resting stage was also recorded and analyzed. In the present study, the resting stage was defined as the pre-task period in a rest box (
Time-frequency spectral analysis
Spectral analysis was used to assess the dominant frequencies in the LFPs during the task. For averaged spectral density estimation during start-box, choice-point and reward occupancy, the multi-taper method was used (
Functional connectivity in network and quantification of information flow
Directed transfer function (DTF) method, a multi-channel parametric method of analysis based on an autoregressive model, has been applied to a number of neurobiological systems for analysis of causal connectivity and flow pattern of activity in the frequency domain (
According to the multivariate autoregressive (MVAR) model, multi-channel LFPs can be represented as a data vector X of N source signals:
The MVAR model can then be expressed as (
Where p is the model order, calculated by the Bayesian information criterion. An is the coefficient matrix of the MVAR model. E(t) is the vector of multivariate zero mean uncorrelated white noise at time t.
f denotes a specific frequency, H(f) is the transfer function matrix, and A(0) = −I, I is an identity matrix.
γij(f) represents the ratio between the effect from channel j to channel iand the combined effect from all other nodes to node j, k is the number of nodes, and H is the transfer matrix of the system.
DTFij refers to the connectivity strength from channel i to channel j. N refers to the number of channels. G refers to the set of channels in the mPFC-MD causal network. In order to measure the information transfer in the network, the bidirectional information flow (IF) between mPFC and MD was calculated, respectively:
IFmPFC→MD refers to the IF from the mPFC to the MD, and IFMD→mPFC indicates the IF from the MD to the mPFC. NmPFC and NMD refer to the number of channels in the mPFC and MD, respectively.
Histology
After the experiments, the mice were deeply anesthetized and transcardially perfused with PBS (phosphate buffered solution) followed by 4% PFA (paraformaldehyde). Brain sections were mounted on slides to visualize and photograph lesions. Recording sites of the mPFC and MD were verified histologically (Figure 1D) and shown overlaid on a representative drawing [taken from the atlas of
Data analysis and statistics
We simultaneously recorded LFPs from the mPFC and MD of eight mice as they performed the DNMTP task. In total, we describe 412 trials (including 295 correct trials and 117 incorrect trials) in the present article. The distribution of trials per subject is shown in Table 1. Data in the text and figures are expressed as the mean ± SEM. Statistical differences were evaluated by using Mann–Whitney test, two-way ANOVA, one-way ANOVA, and Bonferroni’s test for post-hoc analyses. P-values are marked statistically significant as follows: *P < 0.05, **P < 0.01, and ***P < 0.001. n.s., not significant.
TABLE 1
| Mouse | Trials | |
| Correct | Incorrect | |
| 1 | 42 | 18 |
| 2 | 48 | 18 |
| 3 | 22 | 12 |
| 4 | 50 | 13 |
| 5 | 34 | 13 |
| 6 | 41 | 18 |
| 7 | 29 | 12 |
| 8 | 29 | 13 |
| Sum | 295 | 117 |
Distribution of trials per subject.
Results
Increased theta oscillatory activity in medial prefrontal cortex and mediodorsal thalamus during memory retrieval
We simultaneously recorded LFPs from the mPFC and MD of eight mice using custom-made electrode arrays as they performed the DNMTP task that requires thalamic-prefrontal interactions. The examples of LFP traces recorded from the mPFC and MD were exhibited in Figure 1. Spectral analysis was used to assess the dominant frequencies in the LFPs during the task epochs. Figures 1E,F showed the time-frequency power spectrum across correct trials for a single mouse. A strong and transient increase in theta (4–12 Hz) activity in both mPFC and MD was observed preceded an animal’s correct choice during the choice phase. During the sample runs, however, the increased theta power was only found in mPFC but not MD. The time-frequency power spectrum reveals the prominent activity in the theta frequency band during working memory task as well as the mPFC and MD are involved in memory retrieval.
Increased information flow in theta frequency band from medial prefrontal cortex to mediodorsal thalamus for correct choice
Information flows between mPFC and MD were assessed to reveal the nature of information transfer in the mPFC-MD network underlying spatial working memory. Firstly, the information flow across frequency bands from mPFC to MD for individual subject was reliably evaluated during the choice phase (Figure 2A and Supplementary Figure 1). To statistically quantify this result, we made a comparison of the average information flow between the correct and incorrect trials and the information flow in theta frequency band showed a the dominant role on correct trials (Figure 2A and Supplementary Tables 1,2). Moreover, the information flow in theta shows a remarkably increase on correct trials; however, no consistent difference among different frequency bands on incorrect trials across subjects (Figure 2A and Supplementary Figure 1).
FIGURE 2

Theta-band information flow from mPFC to MD is prominent during the choice phase. (A) The average information flow across different frequencies during the choice phase. (Left) Information flow as a function of frequency across subjects (n = 8 mice). The red and blue curves represent the information flow on correct and incorrect trials, respectively. (Right) Comparison of information flow across different frequencies on correct and incorrect trials (correct: two-way AONVA, F = 55.02, P < 0.001; incorrect: two-way ANOVA, P > 0.05). (B) Same as (A), but for the sample phase (correct: two-way AONVA, P > 0.05; incorrect: two-way ANOVA, P > 0.05). (C) Running speed and duration on correct and incorrect trials showed no significant difference (Mann–Whitney test, speed: P > 0.05; duration: P > 0.05). (D) Running speed and duration between the choice and sample phases showed no significant difference (Mann–Whitney test, speed: P > 0.05; duration: P > 0.05). ***P < 0.001, ns, not significant.
To further investigate whether the enhanced theta-band information flow we observed is specific to the choice phase when memory retrieval is presumed to occur, we also analyzed the information flow during the sample phase. The results showed that there was no significant increase in the information flow between the mPFC and MD (Figure 2B). These findings suggested that the theta-band information transfer from mPFC to MD may serve the retrieval of working memory information.
In order to test whether the movements differ in the different trial types, we compared the running speed and the duration to choice between the correct and incorrect trials. The results showed no significant difference in the running speed and duration between the two the correct and incorrect trials (Figure 2C). The results indicated that the theta-band information flow were not influenced by animals’ locomotion we further compared the running speed between the choice run and sample sun. The results showed no significant difference in the running speed between the choice and sample phase (Figure 2D). These results suggest a role for theta-specific information transfer from mPFC to MD for successful memory retrieval.
We next subsampled the number of correct trials used for this analysis to match the number of incorrect trials for each animal and repeated this analysis to rule out potential differences to the number of trials in each condition. The increased theta-band information flow was also found in the re-sampled dataset. Moreover, the sub-sampled dataset of trials did not differ in running speed or duration (Supplementary Figure 4). The results showed that the difference in the numbers of correct and incorrect trials did not affect the analysis results.
Dynamic changes in theta-specific information transfer from medial prefrontal cortex to mediodorsal thalamus during memory retrieval
To gain more insights into how mPFC-MD information flow varies during memory retrieval, we next analyzed the temporal dynamics of information flow in theta band from mPFC to MD throughout choice phase. The information flow was found to experience the pattern of increase, peak, and decline on correct trials. Notably, the information flow was increased before the mice turned left or right during the choice phase on correct trials (Supplementary Figure 2A and Figure 3A). By contrast, no significant change in information flow was found on incorrect choices. We then compared the information flow in theta band between the trials with correct and incorrect choices. Averaged information flows on correct trials were statistically higher than those on incorrect trials. Besides, the measures of information flow during resting state were significantly lower than those in other conditions (Supplementary Figure 2B and Figure 3A, correct: 0.0155 ± 0.0009, incorrect: 0.0086 ± 0.0007, and rest: 0.0051 ± 0.0005, P < 0.001). These findings indicate that the temporal dynamics of theta-specific information transfer from mPFC to MD is required for successful memory retrieval. Moreover, the information flow in theta shows a remarkably increase during the choice phase on correct trials; however, no consistent difference among different frequency bands on incorrect trials across subjects (Figure 3B). For a better link the information flow we observed to a memory retrieval process, we further analyzed the information flow during the sample runs. The results showed that there was no significant increase in the information flow between the mPFC and MD (Figures 3C,D). These findings suggested that the theta-specific information transfer from mPFC to MD may serve the retrieval of working memory information.
FIGURE 3

Enhanced theta-band information flow from mPFC to MD correlates with the task epochs. (A) The average theta-band information flow from mPFC to MD during the choice phase across subjects (n = 8 mice). (left) Changes in information flow over position [from start box (location “d”) to arrival at the reward port (location “f”)]. (Right) Comparison of theta-band information flow on correct, incorrect trials and rest stage (two-way ANOVA, F = 312.9, P < 0.001). (B) Theta-band information flow shows a remarkably increase on correct trials [ΔIFmPFC→MD is defined as peak of IF (the maximum of IF during the choice period: d→e) minus the initial value of IF (at location “d”)]. Correct: two-way AONVA, F = 43.31, P < 0.001; incorrect: two-way ANOVA, P > 0.05. (C) same as (A), but for sample phase. The information flow showed no difference in the different trial types and the values on correct and incorrect trials were apparently higher than that at rest stage (two-way ANOVA, F = 274.8, P < 0.001). (D) Same as (B), but for sample phase [ΔIFmPFC→MD is defined as peak of IF (the maximum of IF during the sample period: a→b) minus the initial value of IF (at location “a”)]. Correct: two-way AONVA, F = 12.13, P < 0.001; incorrect: two-way ANOVA, F = 44.29, P < 0.001. *P < 0.05, **P < 0.01, and ***P < 0.001. ns, not significant.
Directional medial prefrontal cortex-mediodorsal thalamus information flow is required during memory retrieval
In order to test the directionality of the coordinated information processing between the mPFC and MD, we further compared the bidirectional information flow (in theta frequency band) between mPFC and MD during the task epochs. The results showed that the information flow from mPFC to MD is differentially modulated during the choice phase. Specially, on correct trials, the information flow from mPFC to MD showed a markedly increase prior to the choice-points. By contrast, there was no distinct change in the reverse direction (Figure 4A and Supplementary Figure 3A). Statistical results showed the information flow from mPFC to MD was significantly higher than the reverse (Figure 4A and Supplementary Figure 3B, mPFC to MD: 0.0157 ± 0.0006, MD to mPFC: 0.011 ± 0.0001, Mann–Whitney test, P < 0.001). However, during the sample phase, there was no statistically significant change in the information flow between the mPFC and MD (Figure 4B). These findings demonstrate that the directional information transfer from mPFC to MD may serve the retrieval of working memory information.
FIGURE 4

Directional information flow from mPFC to MD is distinct during the choice phase. (A) The bidirectional information flow between mPFC and MD across subjects (n = 8 mice) during the choice phase on correctly performed trials. (Left) The red and green curve represents the information flow from mPFC to MD and from MD to mPFC, respectively. (Right) The information flow from mPFC to MD was significantly higher than the reverse (Mann–Whitney test: P < 0.001). (B) Same as (A), but for sample (Mann–Whitney test: P > 0.05). ***P < 0.001, ns, not significant.
Discussion
In this study, we analyzed LFPs from mPFC and MD of mice performing a spatial working memory task and examined the temporal dynamics of functional interactions and directional information transfer between mPFC and MD during memory retrieval processing. Our results showed a significant increased information flow from mPFC to MD, varied in time and frequency, accompanying correct choice. Elevated theta information flow, a feature that was absent on error trials, indicates an important role of the directional information transfer from mPFC to MD for memory retrieval.
It is well known that successful spatial working memory performance depends on organized communication across expansive brain networks. The PFC, as a principal neural substrate of spatial working memory (
Emerging evidence indicates the thalamus has an extensive anatomical connectivity with the PFC. In particular, the MD (mediodorsal nucleus of the thalamus) shares dense, direct and reciprocal excitatory projections with various PFC subregions (
Notably, hippocampus is another substrate which has robust long-range anatomical connectivity and functional interactions with the PFC (
Strikingly, in recent years, neuroscientists have made significant advances toward an understanding of the dynamic communication of the hippocampal-thalamo-PFC network that supports spatial working memory. The three different phases (encoding, maintenance, and retrieval) in working memory exhibit distinct functional interactions among brain networks. Specially, the PFC is generally involved in goal-directed actions, path planning, and strategy switching. The vHPC inputs to the mPFC support spatial encoding (
Another point worth discussing in this study is the theta-specific pattern for successful memory retrieval. Previous studies have revealed that theta oscillations play an essential role in working memory. In human, theta activity can be engaged by the working memory task (
Recently, cross-frequency coupling (CFC) has been proposed as an index neural interactions of information gating and communication within and between broad networks in cognition (
Statements
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The animal study was reviewed and approved by the Tianjin Medical University Animal Care and Use Committee (license number: TMUaMEC 2021060).
Author contributions
XT, XZ, and WB designed the experiment. JW, SZ, and TL carried out the experiments. JW analyzed the data. JW and WB wrote the manuscript. All authors read and approved the final manuscript.
Funding
This work was supported by the National Natural Science Foundation of China (61971311, 61976157, and 61501329).
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.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnins.2022.1055986/full#supplementary-material
Supplementary Figure 1Information flow from mPFC to MD varies with frequency during the choice phase. Each row represents the data for individual subject. Information flow in theta frequency band was the highest among all frequency bands on correct trials (one-way ANOVA, Mouse 1: F = 104.3, P < 0.001; Mouse 2: F = 242.4, P < 0.001; Mouse 3: F = 85.60, P < 0.001; Mouse 4: F = 52.19, P < 0.001; Mouse 5: F = 72.81, P < 0.001; Mouse 6: F = 31.55, P < 0.001; Mouse 7: F = 304.9, P < 0.001; and Mouse 8: F = 97.53, P < 0.001). No consistent difference was found among different frequency bands on incorrect trials (one-way ANOVA, Mouse 1: F = 60.19, P < 0.001; Mouse 2: F = 38.60, P < 0.001; Mouse 3: F = 2.515, P > 0.05; Mouse 4: F = 9.953, P < 0.001; Mouse 5: F = 9.089, P < 0.001; Mouse 6: F = 32.25, P < 0.001; Mouse 7: F = 76.41, P < 0.001; and Mouse 8: F = 366.1, P < 0.001).
Supplementary Figure 2Theta-band information flow from mPFC to MD varies during the choice phase. (A) Information flow from mPFC to MD during the choice phase. The curves show the results from start box (location “d”) to arrival at the reward port (location “f”) and each panel represents individual subject. (B) Comparison of the information flow in the different conditions (One-way ANOVA, Mouse 1: F = 43.32, P < 0.001; Mouse 2: F = 96.72, P < 0.001; Mouse 3: F = 275.1, P < 0.001; Mouse 4: F = 149.9, P < 0.001; Mouse 5: F = 121.4, P < 0.001; Mouse 6: F = 271.6, P < 0.001; Mouse 7: F = 87.67, P < 0.001; and Mouse 8: F = 463.4, P < 0.001).
Supplementary Figure 3Directional theta-band information flow from mPFC to MD is prominent during the choice phase. (A) Bidirectional information flow between mPFC and MD on correct trials. (B) Comparison of the information flow in the two directions (t-test, Mouse 1: P < 0.01; Mouse 2: P < 0.01; Mouse 3: P < 0.05; Mouse 4: P < 0.001; Mouse 5: P < 0.01; Mouse 6: P < 0.01; Mouse 7: P < 0.01; and Mouse 8: P < 0.001).
Supplementary Figure 4Information flow from mPFC to MD for the subsampled dataset. (A) Information flow across different frequencies during the choice phase. (Left) Information flow as a function of frequency across subjects (n = 8 mice, 295 correct trials and 117 incorrect trials). (Right) Comparison of the information flow across different frequencies on correct and incorrect trials (correct: two-way AONVA, F = 55.02, P < 0.001; incorrect: two-way ANOVA, P > 0.05). (B) Same as (A), but for subsampled dataset (n = 8 mice, 117 correct trials and 117 incorrect trials; correct: two-way AONVA, F = 69.87, P < 0.001; right, incorrect: two-way ANOVA, P > 0.05). (C) Comparison of running speed and duration during the choice phase (Mann–Whitney test: Speed P > 0.05; duration P > 0.05). (D) Same as (C), but for subsampled dataset (Mann–Whitney test: Speed P > 0.05; duration: P > 0.05).
References
1
AbubakerM.Al QasemW.KvasnakE. (2021). Working memory and cross-frequency coupling of neuronal oscillations.Front. Psychol.12:756661. 10.3389/fpsyg.2021.756661
2
AlekseichukI.TuriZ.Amador de LaraG.AntalA.PaulusW. (2016). Spatial working memory in humans depends on theta and high gamma synchronization in the prefrontal cortex.Curr. Biol.261513–1521. 10.1016/j.cub.2016.04.035
3
AlexinskyT. (2001). Differential effect of thalamic and cortical lesions on memory systems in the rat.Behav. Brain Res.122:16. 10.1016/s0166-4328(01)00182-6
4
BaddeleyA. (1992). Working memory.Science255556–559. 10.1126/science.1736359
5
BaddeleyA. (2010). Working memory.Curr. Biol.20R136–R140. 10.1016/j.cub.2009.12.014
6
BahnerF.DemanueleC.SchweigerJ.GerchenM. F.ZamoscikV.UeltzhofferK.et al (2015). Hippocampal-dorsolateral prefrontal coupling as a species-conserved cognitive mechanism: A human translational imaging study.Neuropsychopharmacology401674–1681. 10.1038/npp.2015.13
7
BeatyR. E.BenedekM.SilviaP. J.SchacterD. L. (2016). Creative cognition and brain network dynamics.Trends Cogn. Sci.2087–95. 10.1016/j.tics.2015.10.004
8
BolkanS. S.StujenskeJ. M.ParnaudeauS.SpellmanT. J.RauffenbartC.AbbasA. I.et al (2017). Thalamic projections sustain prefrontal activity during working memory maintenance.Nat. Neurosci.20987–996. 10.1038/nn.4568
9
BraunU.HarneitA.PergolaG.MenaraT.SchaferA.BetzelR. F.et al (2021). Brain network dynamics during working memory are modulated by dopamine and diminished in schizophrenia.Nat. Commun.12:3478. 10.1038/s41467-021-23694-9
10
BresslerS. L.KumarA.SingerI. (2021). Brain synchronization and multivariate autoregressive (MVAR) modeling in cognitive neurodynamics.Front. Syst. Neurosci.15:638269. 10.3389/fnsys.2021.638269
11
ChenJ.TamA.KebetsV.OrbanC.OoiL. Q. R.AsplundC. L.et al (2022). Shared and unique brain network features predict cognitive, personality, and mental health scores in the ABCD study.Nat. Commun.13:2217. 10.1038/s41467-022-29766-8
12
ColeM. W.YarkoniT.RepovsG.AnticevicA.BraverT. S. (2012). Global connectivity of prefrontal cortex predicts cognitive control and intelligence.J. Neurosci.328988–8999. 10.1523/JNEUROSCI.0536-12.2012
13
D’EspositoM.PostleB. R. (2015). The cognitive neuroscience of working memory.Annu. Rev. Psychol.66115–142. 10.1146/annurev-psych-010814-015031
14
FarrokhiA.TafakoriS.DaliriM. R. (2022). Dynamic theta-modulated high frequency oscillations in rat medial prefrontal cortex during spatial working memory task.Physiol. Behav.254:113912. 10.1016/j.physbeh.2022.113912
15
FellJ.AxmacherN. (2011). The role of phase synchronization in memory processes.Nat. Rev. Neurosci.12105–118. 10.1038/nrn2979
16
FogelH.Levy-LamdanO.ZifmanN.HillerT.EfratiS.SuzinG.et al (2021). Brain network integrity changes in subjective cognitive decline: A possible physiological biomarker of dementia.Front. Neurol.12:699014. 10.3389/fneur.2021.699014
17
FunahashiS. (2017). Working memory in the prefrontal cortex.Brain Sci.7:49. 10.3390/brainsci7050049
18
GilsonM.Zamora-LopezG.PallaresV.AdhikariM. H.SendenM.CampoA. T.et al (2020). Model-based whole-brain effective connectivity to study distributed cognition in health and disease.Netw. Neurosci.4338–373. 10.1162/netn_a_00117
19
GordonE. M.BreedenA. L.BeanS. E.VaidyaC. J. (2014). Working memory-related changes in functional connectivity persist beyond task disengagement.Hum. Brain Mapp.351004–1017. 10.1002/hbm.22230
20
GriffinA. L. (2021). The nucleus reuniens orchestrates prefrontal-hippocampal synchrony during spatial working memory.Neurosci. Biobehav. Rev.128415–420. 10.1016/j.neubiorev.2021.05.033
21
HallockH. L.WangA.GriffinA. L. (2016). Ventral midline thalamus is critical for hippocampal-prefrontal synchrony and spatial working memory.J. Neurosci.368372–8389. 10.1523/JNEUROSCI.0991-16.2016
22
HelfrichR. F.KnightR. T. (2016). Oscillatory dynamics of prefrontal cognitive control.Trends Cogn. Sci.20916–930. 10.1016/j.tics.2016.09.007
23
HerwegN. A.SolomonE. A.KahanaM. J. (2020). Theta oscillations in human memory.Trends Cogn. Sci.24208–227. 10.1016/j.tics.2019.12.006
24
HooverW. B.VertesR. P. (2007). Anatomical analysis of afferent projections to the medial prefrontal cortex in the rat.Brain Struct. Funct.212149–179. 10.1007/s00429-007-0150-4
25
HorstN. K.LaubachM. (2009). The role of rat dorsomedial prefrontal cortex in spatial working memory.Neuroscience164444–456. 10.1016/j.neuroscience.2009.08.004
26
JadhavS. P.RothschildG.RoumisD. K.FrankL. M. (2016). Coordinated excitation and inhibition of prefrontal ensembles during awake hippocampal sharp-wave ripple events.Neuron90113–127. 10.1016/j.neuron.2016.02.010
27
JarvisM. R.MitraP. P. (2001). Sampling properties of the spectrum and coherency of sequences of action potentials.Neural Comput.13717–749. 10.1162/089976601300014312
28
JayachandranM.LinleyS. B.SchlechtM.MahlerS. V.VertesR. P.AllenT. A. (2019). Prefrontal pathways provide top-down control of memory for sequences of events.Cell Rep.28640–654.e6. 10.1016/j.celrep.2019.06.053
29
KaminskiM. J.BlinowskaK. J. (1991). A new method of the description of the information flow in the brain structures.Biol. Cybern.65203–210. 10.1007/BF00198091
30
KupferschmidtD. A.GordonJ. A. (2018). The dynamics of disordered dialogue: Prefrontal, hippocampal and thalamic miscommunication underlying working memory deficits in schizophrenia.Brain Neurosci. Adv.21–15. 10.1177/2398212818771821
31
MillerE. K.CohenJ. D. (2001). An integrative theory of prefrontal cortex function.Annu. Rev. Neurosci.24167–202. 10.1146/annurev.neuro.24.1.167
32
MillerR. L. A.FrancoeurM. J.GibsonB. M.MairR. G. (2017). Mediodorsal thalamic neurons mirror the activity of medial prefrontal neurons responding to movement and reinforcement during a dynamic DNMTP task.eNeuro41–18. 10.1523/ENEURO.0196-17.2017
33
MitchellA. S. (2015). The mediodorsal thalamus as a higher order thalamic relay nucleus important for learning and decision-making.Neurosci. Biobehav. Rev.5476–88. 10.1016/j.neubiorev.2015.03.001
34
MoraschiM.MascaliD.TommasinS.GiliT.HassanI. E.FratiniM.et al (2020). Brain network modularity during a sustained working-memory task.Front. Physiol.11:422. 10.3389/fphys.2020.00422
35
O’NeillP. K.GordonJ. A.SigurdssonT. (2013). Theta oscillations in the medial prefrontal cortex are modulated by spatial working memory and synchronize with the hippocampus through its ventral subregion.J. Neurosci.3314211–14224. 10.1523/JNEUROSCI.2378-13.2013
36
OsadaT.OgawaA.SudaA.NakajimaK.TanakaM.OkaS.et al (2021). Parallel cognitive processing streams in human prefrontal cortex: Parsing areal-level brain network for response inhibition.Cell Rep.36:109732. 10.1016/j.celrep.2021.109732
37
ParkH. J.FristonK. (2013). Structural and functional brain networks: From connections to cognition.Science342:1238411. 10.1126/science.1238411
38
ParnaudeauS.BolkanS. S.KellendonkC. (2018). The mediodorsal Thalamus: An essential partner of the prefrontal cortex for cognition.Biol. Psychiatry83648–656. 10.1016/j.biopsych.2017.11.008
39
ParnaudeauS.O’NeillP. K.BolkanS. S.WardR. D.AbbasA. I.RothB. L.et al (2013). Inhibition of mediodorsal thalamus disrupts thalamofrontal connectivity and cognition.Neuron771151–1162. 10.1016/j.neuron.2013.01.038
40
PaxinosG.CharlesW. (2005). The rat brain in stereotaxic coordinates.London: Elsevier Academic Press.
41
ProskovecA. L.Heinrichs-GrahamE.WilsonT. W. (2016). Aging modulates the oscillatory dynamics underlying successful working memory encoding and maintenance.Hum. Brain Mapp.372348–2361. 10.1002/hbm.23178
42
RatcliffeO.ShapiroK.StaresinaB. P. (2022). Fronto-medial theta coordinates posterior maintenance of working memory content.Curr. Biol.322121–2129.e3. 10.1016/j.cub.2022.03.045
43
RiddleJ.ScimecaJ. M.CellierD.DhananiS.D’EspositoM. (2020). Causal evidence for a role of theta and alpha oscillations in the control of working memory.Curr. Biol.301748–1754.e4. 10.1016/j.cub.2020.02.065
44
RouxF.UhlhaasP. J. (2014). Working memory and neural oscillations: Alpha-gamma versus theta-gamma codes for distinct WM information?Trends Cogn. Sci.1816–25. 10.1016/j.tics.2013.10.010
45
SchmittL. I.WimmerR. D.NakajimaM.HappM.MofakhamS.HalassaM. M. (2017). Thalamic amplification of cortical connectivity sustains attentional control.Nature545219–223. 10.1038/nature22073
46
SederbergP. B.KahanaM. J.HowardM. W.DonnerE. J.MadsenJ. R. (2003). Theta and gamma oscillations during encoding predict subsequent recall.J. Neurosci.2310809–10814. 10.1523/JNEUROSCI.23-34-10809.2003
47
SethA. K. (2010). A MATLAB toolbox for Granger causal connectivity analysis.J. Neurosci. Methods186262–273. 10.1016/j.jneumeth.2009.11.020
48
SiegelM.DonnerT. H.EngelA. K. (2012). Spectral fingerprints of large-scale neuronal interactions.Nat. Rev. Neurosci.13121–134. 10.1038/nrn3137
49
SigurdssonT.StarkK. L.KarayiorgouM.GogosJ. A.GordonJ. A. (2010). Impaired hippocampal-prefrontal synchrony in a genetic mouse model of schizophrenia.Nature464763–767. 10.1038/nature08855
50
SpellmanT.RigottiM.AhmariS. E.FusiS.GogosJ. A.GordonJ. A. (2015). Hippocampal-prefrontal input supports spatial encoding in working memory.Nature522309–314. 10.1038/nature14445
51
StoutJ. J.GriffinA. L. (2020). Representations of on-going behavior and future actions during a spatial working memory task by a high firing-rate population of medial prefrontal cortex neurons.Front. Behav. Neurosci.14:151. 10.3389/fnbeh.2020.00151
52
TamuraM.SpellmanT. J.RosenA. M.GogosJ. A.GordonJ. A. (2017). Hippocampal-prefrontal theta-gamma coupling during performance of a spatial working memory task.Nat. Commun.8:2182. 10.1038/s41467-017-02108-9
53
van AsselenM.KesselsR. P.NeggersS. F.KappelleL. J.FrijnsC. J.PostmaA. (2006). Brain areas involved in spatial working memory.Neuropsychologia441185–1194. 10.1016/j.neuropsychologia.2005.10.005
54
VanceA.WintherJ. (2021). Spatial working memory performance in children and adolescents with major depressive disorder and dysthymic disorder.J. Affect. Disord.278470–476. 10.1016/j.jad.2020.09.093
55
VertesR. P. (2006). Interactions among the medial prefrontal cortex, Hippocampus and midline thalamus in emotional and cognitive processing in the rat.Neuroscience1421–20. 10.1016/j.neuroscience.2006.06.027
56
VertesR. P.HooverW. B.Do ValleA. C.ShermanA.RodriguezJ. J. (2006). Efferent projections of reuniens and rhomboid nuclei of the thalamus in the rat.J. Comp. Neurol.499768–796. 10.1002/cne.21135
57
VogelP.HahnJ.DuvarciS.SigurdssonT. (2022). Prefrontal pyramidal neurons are critical for all phases of working memory.Cell Rep.39:110659. 10.1016/j.celrep.2022.110659
58
WangX.WangR.LiF.LinQ.ZhaoX.HuZ. (2020). Large-scale granger causal brain network based on resting-state fMRI data.Neuroscience425169–180. 10.1016/j.neuroscience.2019.11.006
59
XiaM.LiuT.BaiW.ZhengX.TianX. (2019). Information transmission in HPC-PFC network for spatial working memory in rat.Behav. Brain Res.356170–178. 10.1016/j.bbr.2018.08.024
60
XieY.LiY.DuanH.XuX.ZhangW.FangP. (2021). Theta oscillations and source connectivity during complex audiovisual object encoding in working memory.Front. Hum. Neurosci.15:614950. 10.3389/fnhum.2021.614950
61
ZhangL.GlascherJ. (2020). A brain network supporting social influences in human decision-making.Sci. Adv.6:eabb4159. 10.1126/sciadv.abb4159
Summary
Keywords
spatial working memory, medial prefrontal cortex (mPFC), mediodorsal thalamus (MD), local field potentials (LFPs), information flow
Citation
Wang J, Zhang S, Liu T, Zheng X, Tian X and Bai W (2022) Directional prefrontal-thalamic information flow is selectively required during spatial working memory retrieval. Front. Neurosci. 16:1055986. doi: 10.3389/fnins.2022.1055986
Received
28 September 2022
Accepted
08 November 2022
Published
23 November 2022
Volume
16 - 2022
Edited by
Anirban Dutta, University of Lincoln, United Kingdom
Reviewed by
Scott S. Bolkan, Princeton University, United States; Yang Yang, The Pennsylvania State University, United States
Updates

Check for updates
Copyright
© 2022 Wang, Zhang, Liu, Zheng, Tian and Bai.
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: Wenwen Bai, baiwenwen@tmu.edu.cn
This article was submitted to Neural Technology, a section of the journal Frontiers in Neuroscience
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.