- 1School of Medicine, Chongqing University, Chongqing, China
- 2College of Mathematics and Statistics, Chongqing University, Chongqing, China
- 3Department of Psychological and Cognitive Sciences, Tsinghua University, Beijing, China
- 4Department of Physics, South China University of Technology, Guangzhou, China
Sequential working memory (SWM), referring to the temporary storage and manipulation of information in order, plays a fundamental role in brain cognitive functions. The serial position effect refers to the phenomena that recall accuracy of an item is associated to the order of the item being presented. The neural mechanism underpinning the serial position effect remains unclear. The synaptic mechanism of working memory proposes that information is stored as hidden states in the form of facilitated neuronal synapse connections. Here, we build a continuous attractor neural network with synaptic short-term plasticity (STP) to explore the neural mechanism of the serial position effect. Using a delay recall task, our model reproduces the the experimental finding that as the maintenance period extends, the serial position effect transitions from the primacy to the recency effect. Using both numerical simulation and theoretical analysis, we show that the transition moment is determined by the parameters of STP and the interval between presented stimulus items. Our results highlight the pivotal role of STP in processing the order information in SWM.
1 Introduction
Sequential working memory (SWM), a function responsible for temporarily storing and manipulating information in a specific order (Stephan and MJB, 1989; Jensen and Lisman, 2005; Logan, 2021), plays a fundamental role in brain cognitive functions, such as reasoning, comprehension and learning (Alan, 2003; Curtis and Lee, 2010; Potagas et al., 2011; Calmels et al., 2012; Ru et al., 2022). SWM supports human mental processes by providing an interface between perception, long-term memory and actions (Tsetsos et al., 2012). The memory recall paradigm is widely utilized to investigate the storage of multiple items in SWM (Endel and CFI, 2000; Pantelis et al., 2008; Henson, 2013), where subjects are required to retrieve previously presented information. A large volume of psychological experiments has demonstrated that the retrieval performances of human subjects are associated with the order of presented items, displaying a serial position effect, namely, subjects exhibit better performances for memory items appearing at the beginning or at end of a sequence, called the primacy or recency effect, respectively (Simon, 1962; Postman and Phillips, 1965; Glanzer and Cunitz, 1966). This serial position effect is observed in various types of working memory systems, including visual (Kiani et al., 2008), auditory (Hurlstone et al., 2014; Borderie et al., 2024), and spatial working memories (Groeger et al., 2008). The serial position effect is a well-established phenomenon in memory research, yet its underlying neural mechanism, contextual variation, and functional implication remain largely unresolved.
A number of psychophysical experiments have indicated that the serial position effect of SWM is affected by factors related to attention, context and interference. For instance, the attentional gradient, i.e., a gradual decrease in attention level as different items are presented during the encoding period, affects the primacy effect. Contents can also affect participants' retrieval performances, with the recency effect observed when the recall cue is the item order, while the primacy effect observed when the recall cue is the relative size of a specific feature of the item (Cowan et al., 2002; Li et al., 2021). Interference between items affects the recency effect (Gorgoraptis et al., 2011). Additionally, some variations in the experimental paradigm can affect the serial position effect, such as, increasing the inter-stimulus interval during the encoding can weaken the primacy effect but not the recency effect (Glanzer and Cunitz, 1966); prolonging the maintenance period in a visual sequence working memory task can shift participants' performances from the recency to the primacy effect (Knoedler et al., 1999), and this recency-primacy shift were observed in experiments including auditory, verbal, and text materials (Knoedler et al., 1999; Storm and Bjork, 2016). Finally, distractions at different time points during the maintenance period can lead to fluctuations in participants' retrieval accuracy (Lui et al., 2023).
Up to now, the neural mechanism underlying the serial position effect in SWM remains largely unclear. The perspective of “limited resource” proposed that the differential allocation of memory resources across multiple items governs their relative recall precision, thereby leading to the primacy and recency effects as observed in SWM tasks (Gorgoraptis et al., 2011; Ma et al., 2014; Lee et al., 2020; Wang et al., 2021). Another study proposed an attractor network model with firing rate adaptation, which explains the power law of recall capacity (Naim et al., 2020), as well as the primacy and recency effects in human free recall (Boboeva et al., 2021). Nevertheless, these studies did not explain how the detailed dynamics of a neural circuit account for the serial position effect. Recent studies have suggested that working memory is mediated by rapid transitions in “activity-silent” neural states (Wolff et al., 2017; Barbosa et al., 2020), and the strength of hidden-state representation predicts the accuracy of working memory-guided behavior, including recall precision, i.e., the primacy and recency effects (Stokes, 2015; Wolff et al., 2015, 2017; Katkov et al., 2017; Naim et al., 2020). The synaptic mechanism of working memory posits that information is encoded in the facilitated synaptic connections between neurons, rather than in the persistent responses of neurons (Mongillo et al., 2008; Mi et al., 2017). These works mainly studied the neural mechanism for storing and manipulating a single memory item. However, the neural circuit dynamics underlying the serial position effect during the storage of multiple (>2) memory items has not been investigated, which is the focus of the present study.
In this work, we adopt the view that information resides in hidden states of a neural circuit (Stokes, 2015) and is expressed by facilitated synapses between neurons (Mongillo et al., 2008). Specifically, we develop a model of continuous attractor neural network (CANN) with short-term synaptic plasticity (STP). Utilizing a delay task paradigm, we investigate the serial position effect in SWM. In the delay task, the whole period is divided into stimulus encoding, maintenance, and retrieval/response phases. During the encoding phase, participants sequentially receive and encode multiple items into their working memory. After a maintenance period, they are prompted to recall the items, with each item's recall performance indicating the precision of the corresponding memory representation. Typically, participants are required to retain a specific attribute of each item in the sequence, such as visual orientation, direction, or spatial location. Our model shows that with the prolongation of the maintenance period, the serial position effect gradually shifts from a significant primacy effect to a significant recency effect, with the recency effect diminishing in significance over time. This agrees well with the experimental finding. We further analyze that the transition moment of the serial position effect is predominantly determined by the STP dynamics and the inter-item interval of presenting stimuli. Our study highlights the important role of STP plays in processing the order information in SWM.
2 The model
To elucidate the neural mechanism underlying the temporal dynamics of the serial position effect in SWM, we adopted a continuous attractor neural network (CANN) with short-term synaptic plasticity effect (STP). CANNs are a canonical model for neural information storage and representation (Wu et al., 2013, 2016) (Figure 1A), which has been successfully applied to describe the encoding of continuous features in neural systems, such as orientation (Ben-Yishai et al., 1995), moving direction (Georgopoulos et al., 1986), head direction (Taube et al., 1990), and spatial location of objects (Bottomley, 1987). Additionally, CANNs has been extensively used to model the neural mechanism of working memory (Mi et al., 2017; Li et al., 2021). STP is a ubiquitous phenomenon in neural systems, referring to the property that synaptic efficacy between neurons dynamically changes over time in a way that reflects the firing history of the pre-synaptic neuron (Figure 1B). Based on the property of STP, Mongillo et al. proposed a synaptic mechanism of working memory, stating that a neural circuit need not to maintain energy-intensive firings during the entire period of the task for memorizing stimuli, rather the neural circuit can utilize facilitated synaptic connections to retain information (Mongillo et al., 2008). The alteration in synaptic strength induced by STP is a relatively slow process that temporarily modifies the network's connectivity pattern, leading to the network's computation relying on the history of external inputs. Combining CANNs with STP, computational models have elucidated the maintenance and manipulation of working memory (Mi et al., 2017), and the phase precession phenomenon in hippocampus (Chu et al., 2022).
Figure 1. A continuous attractor neural network (CANN) model with short-term plasticity (STP). (A) The schematic diagram of the CANN. Excitatory neurons are arranged in a ring based on their preferred visual orientations θ (θ∈[−π2, π2)). The connection strength between two excitatory neurons at θ and θ′ is denoted as J(θ, θ′), which depends only on |θ−θ′| (the varying shades of gray lines represent the connection strength) and is translation-invariant in the feature space. All excitatory neurons in the network are connected to a global inhibitory neuronal pool (the gray node). The network generates a Gaussian-shaped bump (red curve) to represent the external stimulus (Iext(θ, t)) using population coding strategy. (B) The schematic diagram of STP. u(θ, t), x(θ, t) represent the release probability and the fraction of available neurotransmitters of neurons θ at time t, respectively. The firing rate of neurons at θ(rE(θ, t)) instantaneously increases (the red solid line) upon receiving external signals, leading to an increase in release probability of neurotransmitters (u(θ, t)) and a decrease in the fraction of available resources (x(θ, t)). Following the removal of the external stimulus, u(θ, t) decays to 0 within the time of τf, and x(θ, t) returns to 1 within τd. Due to the dominance of the synaptic short-term facilitation effect, u(θ, t) remains at a high level for an extended duration.
2.1 Continuous attractor neural networks
Consider a one-dimensional continuous stimulus θ, such as the visual orientation, is encoded by an ensemble of neurons. All excitatory neurons in the CANN (the red circles in Figure 1A) are aligned in a ring according to their preference under the periodic boundary condition, i.e., θ∈[−π/2, +π/2), and they are all reciprocally connected to a global inhibitory neuronal pool (the black circle in Figure 1A). Denotes hE(θ, t) as the synaptic inputs at time t of excitatory neurons at θ. The dynamics of hE(θ, t) is determined by a decay term, the recurrent current from other neurons, the inhibitory input from the global inhibitory neuronal pool and the external input, which is written as,
where τ denotes the time constant of neurons, ρ the neuronal density, I0+σ0η0(θ, t) the background input, with η0(θ, t) the Gaussian white noise of zero mean and unit variance and σ0 the corresponding noise strength. Iext(θ, t) refers to the external input, such as the visual stimulus during the encoding period and the cue during the recalling period. JEI is the synaptic strength from the global inhibitory neuronal pool to excitatory neurons. rE(θ, t) is the firing rate of neurons with preference at θ, and its relationship with the synaptic current is given by rE(θ, t) = αln [1+exp(hE(θ, t)/α)], which is a smoothed threshold-linear function.
J(θ, θ′) is the synaptic connection between neurons at θ and θ′, as set in Equation 2,
where J(θ, θ′) is a function of the difference in neuronal preferences (i.e., (θ−θ′)), which is translation-invariant in the feature space. Due to this characteristic topological structure, a CANN can hold a continuous family of stationary states, metaphorically understood as a valley of local minima in the network's energy landscape. J and J0 determine the synaptic connection strength between neurons, while B controls the synaptic interaction range. Neurons with similar preferences have stronger synaptic connections, while those with significantly different preferences have weaker connections.
The synaptic input to the global inhibitory neuronal pool is denoted as hI, with rI the corresponding firing rate, and rI(t) = αln [1+exp(hI(t)/α)]. The dynamics of the global inhibitory neuronal pool is written as,
where τ denotes the time constant of the inhibitory neuronal pool, JIE the connection strength from excitatory neurons in the ring to the inhibitory neuronal pool. The global inhibitory neural pool plays a crucial role in maintaining a balanced state between excitation and inhibition in the network, thereby preventing excessive neuronal firing. Additionally, it fosters competition among different groups of excitatory neurons, ensuring that only one memory item is represented at a given moment.
2.2 Short-term synaptic plasticity
Two types of STP, known as short-term facilitation (STF) and short-term depression (STD), have been observed in various cortical areas. STF is caused by the influx of calcium into the synaptic terminal of the pre-synaptic neuron following spike generation, which increases the release probability of neurotransmitters. On the other hand, STD is caused by the depletion of neurotransmitters at the synaptic terminal of the pre-synaptic neuron after spike generation.
In the model proposed by Mongillo et al. (2008), the STF effect is modeled by u(θ, t) (u∈[0, 1]), which indicates the release probability of neurotransmitters from pre-synaptic neurons at θ, and STD is modeled by x(θ, t)(x∈[0, 1]), indicating the fraction of available neurotransmitters in pre-synaptic neurons. The dynamics of STF and STD are given in Equation 4,
where τf and τd denote the time constants of STF and STD, respectively, and U0 the increment of u caused by spiking of the pre-synaptic neuron. When a neuron at θ receives external inputs, its firing rate (rE(θ, t)) increases. The increase in firing rate results in an increase in the release probability of neurotransmitter u(θ, t) (with the increment determined by U0), leading to the STF effect, while the proportion of available neurotransmitter x(θ, t) decreases, leading to the STD effect. Subsequently, u(θ, t) decays to its baseline of 0 with a time constant τf, and x(θ, t) returns to its baseline of 1 with a time constant τd, as illustrated in Figure 1B. The product of u(θ, t) and x(θ, t) represents the instantaneous synaptic efficacy at time t, i.e., Ju(θ, t)x(θ, t), which reflects the strength of memory representation in the network (Mi et al., 2017; Li et al., 2021).
To elucidate the neural mechanism of SWM, we selected parameters consistent with the synaptic connectivity between neurons in the prefrontal cortex (PFC), which is the primary crotical region involved in working memory (Wang et al., 2006). We adopted the STF dominant parameters as in the model proposed by Mongillo et al. (2008), i.e., τd≪τf and a smaller U0. This implies that after neuron firing, the synaptic connection efficacy is maintained at a high value for an extended period to sustain memory information.
In our numerical simulation, we model the CANN by considering N neurons uniformly distributed in the range of [−π/2, π/2) in the feature space. The integration in Equation (1) is computed by,
The parameters are given in Supplementary Tables 1, 2.
3 Results
3.1 The serial position effect in SWM
Based on the above model of CANN with STP, we investigated the serial position effect in SWM. We first studied the case of two memory items and later generalized the study to the case of multiple items. We utilized the same paradigm as in the psychophysical experiments for SWM (Li et al., 2021), and investigated the recall accuracy of items based on their visual orientations, as illustrated in Figure 2A.
Figure 2. The serial position effect in SWM task for two items. (A) Schematic diagram of the psychophysical experiment paradigm. In the encoding phase, two distinct visual stimuli are sequentially presented. Each stimulus is presented for a duration denoted by Tencode, separated by an time interval of Tgap between them. Following a delay period of Tmaintain, a recall cue lasting for Trecall is presented. Participates are then tasked with recalling the orientation value of either the first or second visual stimulus based on the given cue. (B, C) The temporal dynamics of neural activity pattern in an individual trial simulation example. (B) Two visual stimuli, denoted as θi (i=1,2), are presented sequentially during the encoding phase. In response to each stimulus, the network generates successive bump-shaped neural activity patterns centered at θi, respectively. During the delay period, the neural activity gradually decays to a silent state. (C) The temporal course of synaptic efficacies of neural groups encoding two visual stimuli throughout the trail. When a stimulus is presented, the synaptic efficacy of the corresponding neural group rapidly decreases. After the stimulus is removed, the synaptic efficacy gradually recovers to a maximum value, denoted as , within a certain period of time. Subsequently, it remains at a high level for an extended duration. (D) The recalling performance at varying Tmaintain. Tc denotes the critical moment of recall performance shift from primacy to recency effect, and Tmaintain∈{0.1s, 0.2s, 0.3s, 0.4s, 0.5s, 1s, 2s, 4s, 6s, 8s, 10s}. (Top) The normalized target probability of the ith presented item, denoted as Pi for i = 1, 2, at varying Tmaintain. (Bottom) The normalized target probability difference between the 1st and 2nd item (P2−P1) at varying Tmaintain. (E) The relative synaptic efficacy (ΔJux(t)±SEM) of the neuronal groups encoding the two visual stimuli during maintaining period. denotes the critical moment at which ΔJux(t)≡0. The parameters settings see Supplementary Tables 1, 2. (n.s.: p > 0.05, *: 0.01 <p <0.05, **: 0.001 <p <0.01, ***: p <0.001).
In each trail, two stimuli with different orientations (referred to as θ1, θ2) are sequentially presented in the encoding period. Following a delay period, a visual recall cue lasting for Trecall is presented, and the network retrieves either the first or second visual item based on the recalling cue. Let Tencode denote the duration of presenting each item, Tgap the time interval between two items, and Tmaintain the duration of the delay period. The orientation values of two stimuli are set as: θ1 is randomly selected from the range [−π2, π2), and θ2 = θ1+Δθ, where Δθ is the difference between two stimuli, randomly selected from the data set [±17°, ±24°, ±38°, ±52°, ±66°, ±80°] (Li et al., 2021). The visual stimulus in the encoding period and the cues in the recalling period are denoted as Iext(θ, t), which are written as,
where θtype(type = encode, recall) represents the visual stimulus during different periods. The parameters atype(t) and Btype regulate the strength and accuracy of external signals, respectively. A larger atype(t) and Btype result in more precise encoding of orientation information from the stimulus. As the recalling signals are unrelated to the task, the parameters are set as follows: aencode≫arecall, Bencode>Brecall (Li et al., 2021).
When two visual stimuli are presented sequentially, the neural network generates successive bump-shaped neural activity patterns. The peaks of these bumps correspond to θ1 and θ2, respectively, as illustrated in Figure 2B. Due to the strong interactions among neurons with similar preferences and weak interactions among those with significantly different preferences in the network, we define a neuronal group Gi as the ensemble of neurons whose preferred values satisfy |θ−θi| ≤ Δ (i = 1, 2), and this group of neurons primarily encodes the ith item. The corresponding neural activity and synaptic strength are calculated to be and , respectively, with mi representing the number of neurons in Gi. After removing stimuli, the firing rates of both neuronal groups gradually decrease to zero (Stokes, 2015; Wolff et al., 2015). However, their synaptic strengths remain at high levels due to STF, which maintain the stimulus information (see Figure 2C). During the recalling period, the network generates a weak bump-shaped activity pattern in response to the recalling cue, and the retrieved orientation (denoted as , i = 1, 2) is decoded using the population vector method, with details given in the Supplementary material 1.1.
We investigated how the maintenance duration Tmaintain affects the recall performance. We set Tmaintain as a variable ranging from 0 to 10 s and selected 11 values within this range. For each chosen value of Tmaintain, we evolved the network for 50 times (corresponding to 50 different participants in a psychophysical experiment), each consisting of 300 trials. We utilized the normalized target probability method (Bays et al., 2009; Schneegans and Bays, 2016) to calculate the recall performance of each item (i.e.,, i = 1, 2). More details see Supplementary material 1.2. We found that (as shown in Figure 2D):
• When Tmaintain is smaller than a critical value denoted as Tc, i.e., Tmaintain<Tc, the recall performance exhibits the primacy effect, indicating that participates memorize the first item more accurately. Moreover, the primacy effect becomes more pronounced as the value of Tmaintain decreases.
• When Tmaintain>Tc, the recall performance shifts to the recency effect, indicating that participates memorize the second item more accurately. The significance of the recency effect gradually decreases as Tmaintain increases.
We further utilized the methods of Circular Variance (CV) and Circular Kurtosis (CK; Berens, 2009) to calculate the accuracy of recall performance. The statistical results are consistent with those shown in Figure 2D (more details see Supplementary material 1.2 and Supplementary Figure 1). In conclusion, with the increase of the maintenance period, the serial position effect in SWM dynamically shifts from the primacy effect to the recency effect.
To further reveal the neural mechanism underlying the dynamical change of the serial position effects, we calculated the relative synaptic efficacy of two neuronal groups encoding two stimuli (θ1 and θ2) over time, denoted as ΔJux(t) = Jux2(t)−Jux1(t) hereafter. The synaptic mechanism of WM posits that information is maintained in the facilitated synaptic interactions between neurons, with the synaptic efficacy determining the accuracy of the memorized item (Teyler and Discenna, 1984; Henry and Misha, 1996; Mongillo et al., 2008). For example, when two items (θ1 and θ2) are presented sequentially in a trail (Figure 2C), the transient synchronous firing of a neuronal group (G1 or G2, respectively) leads to rapid decrease in synaptic efficacy, due to the depletion of available neurotransmitters in neurons. After the visual stimulus disappears, the synaptic efficacy of the neural group recovers to the maximum value and maintains at a high level for an extended period. We calculated the synaptic efficacy between two neuron groups ΔJux(t) during the maintaining period (Figure 2E) and found that:
• Since the second item is presented later, its synaptic efficacy Jux2(t) is smaller than that of the first one Jux1(t) before it recovers to the maximum value. Therefore, ΔJux(t) <0 when , where is the moment when ΔJux(t)≡0.
• As the time t approaches , both ΔJux(t) and its variance approach zero. When , ΔJux(t) first increases and then gradually diminishes over time.
The critical moment at which ΔJux(t) = 0 coincides with the value of Tc at which the recall performance transfers from the primacy effect to the recency effect, as shown in Figures 2D, E. Specifically,
• When the maintenance period Tmaintain is smaller than a critical value, i.e., & Tc, the recall performance exhibits the primacy effect; and the greater the value of ΔJux(Tmaintain) is, the more significant the primacy effect becomes.
• As Tmaintain approaches the critical value and ΔJux(Tmaintain) approaches 0, the recall performance switches to the recency effect and is no longer significant (Figure 2D).
• When Tmaintain is much larger than the critical value and meanwhile ΔJux(Tmaintain)≫0, the recall performance displays the recency effect. The greater the value of ΔJux(Tmaintain) is, the more pronounced the recency effect becomes.
It is worth noting that the temporal shift of the serial position effect is independent of the orientation difference between two visual stimuli, which is consistent with the results in Figures 2D, E. More details see Supplementary material and Supplementary Figure 3.
We further investigated the detailed dependence of the transition from the primacy to the recency effect on the model and experimental parameters, including the time interval between two stimuli (Tgap) and the parameters of STP (i.e., τf and τd; Figure 3). For each set of parameters {τf, τd, Tgap}, we calculated the memory accuracy of participants when the recall cue was given at different times (i.e., Tmaintain). For each given Tmaintain, we simulated the network 50 times, each consisting of 300 trails. We then calculated the transition moment from the primacy to the recency effect (i.e., Tc) and the critical moment () when ΔJux(t)≡0. Tc is calculated using different statistical methods, such as CK, CV and P.
Figure 3. The dependence of Tc, on different variables τd, τf and Tgap. The calculations for Tc employ various statistical methods, including normalized target probability (indicated by orange asterisks), Circular Variance (green crosses), and Circular Kurtosis (red hollow circles). The theoretical analyzes represented by the blue solid line, while numerical simulations denoted by purple triangles. (A) Tc and increase with increasing τf; (B) Tc and increase with increasing τd; (C) Tc and decrease with increasing Tgap. in (A–C). For each given parameter set { τd, τf, Tgap }, we computed the memory accuracy for different maintenance times Tmaintain when the recall cue was presented. For every selected Tmaintain, the network was simulated 50 times, with each simulation comprising 300 trials. More details see Supplementary material 4.
We found that both Tc and increase with τf and τd, respectively, as shown in Figures 3A, B, and they both decrease with Tgap (Figure 3C). Furthermore, Tc is approximately equal to for each given parameter set, suggesting that the relative synaptic efficacy between two neural groups (G1 and G2) determines the recall performance. In conclusion, the shift of the serial position effect is determined by STP (i.e., τf, τd) and the inter-stimulus interval (Tgap). Notably, as depicted in Figure 3, the transition from the primacy to the recency effect coincides with the time constant of STD (τd), precisely aligning with its time order.
3.2 Theoretical analysis
In the above, we have utilized a simplified mean-field model to elucidate the neural mechanism of SWM, there are still many variables and parameters involved, including the time constants of STF and STD (τf, τd), the time constant of a single neuron (τ), the connection strength between neurons (i.e., J,J0,JEI,JIE, etc.), the duration of loading each stimulus (Tencode), and the time interval between adjacent stimuli (Tgap). If we continue using numerical simulations, it will be very time-consuming to explore how the recall performance depends on all these variables. We therefore conducted theoretical analysis to elucidate how the critical moment Tc depends on various variables that lead to the shift of the serial position effect. The advantage of theoretical analysis lies in its power of prediction, which can be validated with experiments. Specifically, we focused on examining the dependence of Tc on the variables τf, τd, and Tgap.
To carry our theoretical analysis, we simplified the model of CANN with STP (i.e., Equations 1, 3, 5) into a model composed of multiple neuronal groups. In this simplified model, each ith neuronal group (Gi, i = 1, ⋯ , M with M the number of items in the SWM task) encodes the ith item, and the interaction and overlap between different neuronal groups are ignored. This is because that the orientation difference between two stimuli does not impact much the recall performance, as shown in Supplementary Figure 3. Note that, the value of the visual orientation difference Δθ between memory items reflects the extent of interaction between neural groups encoding them. Using the experimental paradigm with two items as an example, we showed that the recall performances (Supplementary Figure 3) for Δθ being a random number in the range of [−π, +π] are consistent with that when Δθ takes a large value (Figures 2, 3, Δθ∈±17o, ±24o, ±38o, ±52o, ±66o, ±80o), indicating that the ignoring the interactions between neural groups is a proper approximation for theoretical analysis. Indeed, representation of too many items will introduce unignorable overlaps between neural groups. However, consider the limited capacity of working memory (~4 items, Zhang and Luck, 2008), ignoring the interactions between neuron groups is feasible in the theoretical analysis. The STP effect is considered in each neuronal group. Moreover, Figures 2, 3 have shown that the firing rates of different neuronal groups during the maintaining period approach zero (i.e., ri(t) → 0Hz), and this allows us to disregard the change in firing rates over time. On the other hand, the relative memory accuracy of participates for multiple items is mainly determined by the relative synaptic efficacy of neuronal groups during the maintaining period. Therefore, we only need to consider the dynamics of STP of different neuronal groups (Gi) during the maintaining period, which is written as:
where ui(t) and xi(t) denote the STF and STD effects of the ith neuronal groups at time t, respectively, and t is aligned to the presence of the sencond item.
We first examined the case of two items (i.e., M = 2). In accordance with Equation (6), the synaptic efficacy of neuronal group Gi is calculated as , where u0 and x0 represent the values of STF and STD of the ith neuronal group when the stimulus is removed. Thus, the relative synaptic efficacy among neuronal groups during the maintaining period is expressed as:
where . According to Equation (7), the critical moment for ΔJux(t)≡0 is derived as,
For details, see Supplementary material 3. Since it takes an amount of time for the cue in the recalling period to trigger the activity of the corresponding neuronal group, a constant bias (denoted as tb, tb∈(0, Trecal]) is considered into Equation (8). Meanwhile, since τf≫τd, and has the same time scale as τf, we further simplify Equation (8) to be:
Based on the theoretical predictions in Equation (9), we found that the theoretical results for (represented by the blue colored line in Figure 3) are consistent with both the critical moment (depicted by the purple triangle in Figure 3) and Tc (shown by the blue line in Figure 3). It is important to highlight that our theoretical analysis effectively captures the dynamic patterns of storage for multiple items qualitatively during the maintenance period. Therefore, the shift of the serial position effect is positively correlated with the time constants of STD and STF, which implies that a larger time constant of STD enables the working memory system to maintain the primacy effect for a longer duration. This shift of the serial position effect is inversely correlated with the sum of presentation durations of items and the inter-stimulus time intervals (referred to as t*), which indicates that the larger the value of t*, the more difficulty it is for the working memory system to keep the primacy effect.
To further validate the theoretical results, we conducted three different tasks to investigate the temporal dynamics of relative accuracy for multiple memory information (i.e., and Tc). We explored the serial position effect on various parameters of the network and the design of the experiment, as illustrated in Figure 4, where is computed based on Equation (9).
In the first task, we calculated Tc using the same experimental parameters as depicted in Figure 2, but altered the values of τf and τd of the network, as shown in Figure 2A. We selected 21 different values of τf uniformly from the range of [2, 8] s and 21 different values of τd from the range of [0.1, 0.4] s. For each parameter set { τf, τd }, we simulated the network and assessed recall performances of participants at different moments (i.e., Tmaintain), and obtained the critical moment Tc when the recall performance shifted from the primacy to the recency effect (left panel). The selection of Tmaintain followed the same procedure as depicted Figure 3.
In the second task, we fixed τd and varied the variable Tgap in the psychophysical experiment, as well as τf, as shown in Figure 4B. We selected 21 different values of τf uniformly from the range of [2, 8] s and 21 different values of Tgap from the range of [0, 2] s. The calculation method for Tc is the same as that in Figure 4A.
Figure 4. The dependence of temporal dynamics of memory accuracy on different variables and parameters, particularly τf, τd, and Tgap. (left panel) the numerical simulation results of Tc, (right panel) the theoretical predictions (). (A) Both and Tc increases with increasing τf, τd. (B) Both and Tc increases with increasing τf, and decreases with increasing Tgap. (C) Both and Tc increases with increasing τd, and decreases with increasing Tgap.
In the third task, we fixed τf and varied the variable Tgap in the psychophysical experiment, along with τd, as shown in Figure 4C. We uniformly selected 21 different values of τd from the range of [0.1, 0.4] s and 21 different values of Tgap from the range of [0, 2] s. The calculation method for Tc is the same as in Figure 4A.
In summary, we found that: (1) the theoretical analysis of the critical moment () is qualitatively consistent with the results obtained by numerical simulations (Tc). (2) Tc increases with τf and τd, as shown in Figure 4A, (3) Tc decreases with Tgap, as shown in Figures 4B, C.
3.3 Model prediction: the serial position effect in SWM for multiple items
In this section, we demonstrate that the results obtained in SWM with two items can be extended to cases involving multiple items. We continued using the experimental paradigm depicted in Figure 2A for simulations, and studied the case of loading three items successively into the working memory system during the encoding period.
In the WM task, we considered that three items with different orientations are presented sequentially in each trail. Following a maintaining period lasting Tmaintain, a recalling signal with duration of Trecall is presented, which triggers the network to retrieve the task-related feature of the corresponding item. Denote Tencode the loading duration of each item, Tgap the time interval between two adjacent items, and θi, (i = 1, 2, 3) the orientation of each item. The value of θi(i = 1, 2, 3) is determined as follows. θ1 is randomly selected from the range of (−π/2, π/2], and θi is determined by θi = θ1+Δθ(i = 2, 3), where Δθ denotes the orientation difference between the ith (i = 2, 3) and the first items. The value of Δθ is chosen according to the relevant psychophysical experiments (Huang et al., 2023), specifically, Δθ∈[±12°, ±24°, ±36°, ±48°, ±60°, ±72°, ±84°], with Δθ being randomly selected in each trail. For each Tmaintain, we conducted 50 runs, each consisting of 500 trials, and used three methods of normalized target probability, Circular Variance, and Circular Kurtosis to measure the recall performance of each trail. More details are given in Supplementary material 1.2 and Supplementary Figure 2.
In each experimental trial, three items are sequentially loaded to the network, and the network generates three bump-shaped population firing patterns to represent the corresponding items (Figure 5A), respectively. The synaptic efficacy of each neuronal group (Juxi(t), i = 1, 2, 3) decreases rapidly due to the neurotransmitter depletion. After all items are removed, Juxi(t) (i = 1, 2, 3) recovers to its maximum value (Figure 5B) and then retains at a relatively high level to preserve the item information. We calculated the memory performance of participants at different Tmaintain, as shown in Figure 5C.
Figure 5. The temporal pattern of memory accuracy in a SWM task for three items. (A, B) The temporal dynamics of neural activity pattern in an individual trial simulation example. (A) Three items are presented sequentially loaded to the network, and the network generates three bump-shaped population firing pattern to represented the corresponding items, respectively. After the removal of items, the neural activity decays to a slient state. (B) The synaptic efficacy of each neuronal group (Juxi(t) for i = 1, 2, 3) rapidly decays when the corresponding items is presented, and then recovers to a maximum value, denoted as , within a certain time and then remains at a high level for an extended duration. (C) The recall performance at varying Tmaintain. and present the critical moment of recall performance shift from the primacy to recency effects. (Top) The normalized target probability of the ith presented item, denoted as Pi, i = 1, 2, 3, at varying Tmaintain. (Middle) The normalized target probability difference (denoted as P1−P2) and relative synaptic efficacy (Jux1−Jux2, the red line) between the 1st and 2nd items. (Bottom) The normalized target probability difference (denoted as P2−P3) and relative synaptic efficacy (Jux2−Jux3, the red line) between the 2nd and 3rd items. (n.s.: p >0.05, *: 0.01<p <0.05, **: 0.001<p <0.01, ***: p <0.001).
We extended the theoretical analysis in Section 4 to the SWM task involving three items. By comparing the relative synaptic efficacy of neuronal groups encoding different items, we can deduce the network's recalling performance. Due to the neglect of connections and overlaps between different neuronal groups, and the assumption that all neuronal groups receive the same intensity and duration of external signals, the relationship of synaptic efficiency between any two different neuronal groups (for example, Juxi(t), Juxj(t), and i<j, for the ith and j loaded items) during the maintaining period is solved to be , with . Thus, the relative synaptic efficacy between the ith and jth neuronal groups (ΔJuxij(t)) are given by,
The critical moments (Equation 10) can be theoretically resolved based on ΔJux12(t)≡0 and ΔJux23(t)≡0, which are:
where denote the response times of neuronal groups to recalling signals, which are approximately in the time order of τ. According to the above theoretical analysis, if , then Jux1(t)>Jux2(t); otherwise, Jux1(t) <Jux2(t). Similarly, if , then Jux2(t)>Jux3(t); otherwise, Jux2(t) <Jux3(t).
The above theoretical analysis is in perfect alignment with the numerical simulation results, as shown in Equation 11 and Figure 5C.
Firstly, when , the network exhibits the primacy effect (see Figure 5C middle panel). This is because that after the removal of multiple items, Jux1 firstly recovers to the maximum valve, while the synaptic efficacy of other neuronal groups still remains low values due to neurotransmitter depletion (Figure 5B). The larger the relative synaptic efficacy ΔJux12(t), the more significant the primacy effect. Meanwhile, both the value of ΔJux12(t) and the significance of the primacy effect decrease with Tmaintain over time. Thus, the recalling performance indeed shifts around the critical moment (Figure 5C).
Secondly, when , the network exhibits the recency effect (see Figure 5C bottom panel). This is because the synaptic efficacy of the third item (Jux3(t)) recovers to its maximum valve, which is larger than that of other neuronal groups. Furthermore, the synaptic efficacy of neuronal groups gradually decrease over time, and the recency effect becomes insignificant (Figure 5C, bottom panel).
Thirdly, the critical moments ( and ) are primarily influenced by the time constants of STF and STD (τf and τd), as well as the time interval between adjacent items (Tgap). Specifically, the values of and increase with τf and τd, and decrease with Tgap.
4 Discussion
In this study, we built a CANN incorporating STP to investigate the neural mechanism underlying the shift of the serial position effect in SWM. We found that with the elongation of the delay period, the participates' recall performance undergoes a shift from the pronounced primacy effect to the significant recency effect. Additionally, the prominence of the recency effect gradually wanes with the extension of the delay period. Furthermore, we show that the transition moment of the serial position effect is predominantly determined by STP, the time interval between adjacent stimuli, and the duration of stimulus presentation. We carried out theoretical analysis to confirm the simulation results and made predictions to be validated by future experiments. Overall, our study indicates that STP gives us insights into understanding how the ordinal information is processed in working memory.
Utilizing a CANN to depict the encoding, maintenance, and retrieval processes in SWM is biologically plausible. CANNs have been extensively applied to elucidate the neural mechanisms of information processing in working memory. In these tasks, the stimuli are typically represented by continuous variables' including visual orientation (Ben-Yishai et al., 1995), spatial location (Bottomley, 1987), and motion direction (Georgopoulos et al., 1986) in visual working memory, as well as auditory frequency in auditory working memory (Borderie et al., 2024). Additionally, CANNs have been widely used to describe the neural representation and storage of continuous variables, such as head orientation (Stringer et al., 2002; Wang and Kang, 2022), visual orientation (Li et al., 2021), motion direction (Fung et al., 2010; Mi et al., 2014), and spatial location (Samsonovich and McNaughton, 1997; Yoon et al., 2013), and they align with experimental data. A core feature of CANNs is that the synaptic connections between neurons are solely dependent on their difference in preferred continuous variables (Wu et al., 2013, 2016), implying that the synaptic connections in their feature space exhibit spatial translation invariance. Moreover, the network employs population coding (forming bump-shaped neural activity patterns) to represent external stimuli, where the peak value of the bump corresponds to the continuous value of the external stimulus. This structural characteristics of CANNs and its mode of representing external stimuli have been empirically validated.
In this study, we adopted the view of the synaptic mechanism of working memory which considers that information is stored in the facilitated synaptic connections (Mongillo et al., 2008). The prefrontal cortex (PFC) is a crucial cortical region for the execution of working memory. A significant body of empirical evidence demonstrates that the connections between neurons in the PFC exhibit STP and are dominated by STF (Wang et al., 2006; Masse et al., 2020; Bocincova et al., 2022). Consequently, after the removal of external stimuli, information can be stored in the temporally enhanced synaptic connections between neurons, rather than in the sustained firings of neurons (Rainer and Miller, 2002; Shafi et al., 2007), i.e., the memory is residing in the activity-silent hidden state.
The synaptic mechanism of working memory postulates that the memory information is primarily stored in the facilitated synaptic interactions between neurons, and the synaptic efficacy determines the accuracy of the memorized item (Stokes, 2015; Wolff et al., 2015). However, other studies suggest that STD also plays a critical role in the storage and manipulation processes of memory information. Regarding the storage of memory information, the limited capacity of working memory is directly proportional to the time constant of STD (Mi et al., 2017). In term of memory manipulation, for instance, external dynamic perturbations that are irrelevant to the task but weakly related to the attributes of the remembered items can alter the relative synaptic connection strengths between neuronal groups encoding different items through STD, thereby changing the relative accuracy of participates' recollection of different items in visual working memory in real time, a transition from the recency effect to the primacy effect (Knoedler et al., 1999; Li et al., 2021). Moreover, the time constant of STD determines the effective time window for dynamically manipulating working memory. Compared to previous works (particularly Li et al., 2021), the contributions of our study include: (1) We theoretically calculated the critical moment (Tc) at which the relative memory accuracy of two or more items in a SWM task undergoes a transition. We showed that Tc depends on the time constants of short-term facilitation (τf) and short-term depression (τd), and the sum of the duration of presenting the items and the inter-item intervals (t*). (2) Previous studies focused on manipulating two memory items. In this study, we investigated the relative memory accuracy changes of multiple items and derived the critical moments for when these changes occur. Furthermore, our study indicates that the efficacy of synaptic connections (Jux) encoding memory items in neuron groups not only determines the accuracy of item retrieval but is also related to the ordinal information of items, suggesting the importance of STP in processing order information in SWM.
In our model, the efficacy of synaptic connections within neuronal groups, denoted as Jux, determines the accuracy of memory storage. On one hand, u represents the short-term synaptic facilitation effect. According to the synaptic computational theory of working memory, information is retained in the facilitated synaptic connections between neurons within each neuronal group. Therefore, the accuracy of information storage gradually diminishes as u decreases. On the other hand, x represents the short-term synaptic depression effect, which provides the network with a slow negative feedback effect. This negative feedback effect can induce the mobility of neuronal activity bumps within the network (Fung et al., 2012), consequently leading to the phenomenon where neural activity bumps drift in the feature space of the network during the information maintenance period [i.e., the memory drift phenomena reported by Funahashi et al. (1989)]. Such drifting of activity bumps weaken the accuracy of memory representation (Seeholzer et al., 2019). Therefore, during the information maintenance period, a decrease in u within the neuronal groups leads to memory accuracy decay; simultaneously, if the STD effect in the network is sufficiently strong, this effect can also induce the drifting of weak neuronal activities during the maintenance period, further reducing memory accuracy decay.
The design of the psychophysical experimental paradigm also affects the shift of the serial position effect. Our theoretical analysis and simulation results indicate that an increase in the duration of stimulus and the inter-stimulus interval can affect the relative synaptic efficacy between neuronal groups encoding different items. Consequently, this reduces the prominence of the primacy effect and its significance, and concurrently induces the shift of the serial position effect. These theoretical insights into the effects of the psychophysical experimental design on the serial position phenomena pave the way for further investigation into the control of SWM with multiple items.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.
Author contributions
JZ: Writing – original draft, Conceptualization, Formal analysis, Resources, Software, Validation, Visualization. LG: Writing – original draft, Data curation, Formal analysis, Methodology. XH: Writing – original draft. CM: Writing – original draft. YM: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing - original draft, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by the National Natural Science Foundation of China (No: T2122016, YM), National Science and Technology Innovation 2030 Major Program (Nos. 2021ZD0203700/2021ZD0203705, YM), Guangdong Basic and Applied Basic Research Foundation (Grant No. 2021A1515010500, XH), the National Natural Science Foundation of China (No: 12271064, CM), the Chongqing Talent Support program (Grant No. cstc2022ycjh-bgzxm0169), Natural Science Foundation of Chongqing (Grant No. cstc2021jcyj-msxmX1051), the Fundamental Research Funds for the Central Universities (Grant Nos. 2022CDJJCLK002,2020CDJQY-2001), Key Laboratory of Nonlinear Analysis and Applications, Ministry of Education, and Chongqing Key Laboratory of Analytic Mathematics and Applications.
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/fncom.2024.1430244/full#supplementary-material
References
Alan, B. (2003). Working memory: looking back and looking forward. Nat. Rev. Neurosci. 4:829–839. doi: 10.1038/nrn1201
Barbosa, J., Stein, H., Martinez, R. L., Galan-Gadea, A., Li, S., Dalmau, J., et al. (2020). Interplay between persistent activity and activity-silent dynamics in the prefrontal cortex underlies serial biases in working memory. Nat. Neurosci. 23:1016–1024. doi: 10.1038/s41593-020-0644-4
Bays, P. M., Catalao, R. F. G., and Husain, M. (2009). The precision of visual working memory is set by allocation of a shared resource. J. Vis. 9, 7.1–11. doi: 10.1167/9.10.7
Ben-Yishai, R., Bar-Or, R. L., and Sompolinsky, H. (1995). Theory of orientation tuning in visual cortex. Proc. Natl. Acad. Sci. U. S. A. 92, 3844–3848.
Berens, P. (2009). CircStat: a matlab toolbox for circular statistics. J. Stat. Softw. 31, 1–21. doi: 10.18637/jss.v031.i10
Boboeva, V., Pezzotta, A., and Clopath, C. (2021). Free recall scaling laws and short-term memory effects in a latching attractor network. Proc. Natl. Acad. Sci. U. S. A. 118:e2026092118. doi: 10.1073/pnas.2026092118
Bocincova, A., Buschman, T. J., Stokes, M. G., and Manohar, S. G. (2022). Neural signature of flexible coding in prefrontal cortex. Proc. Natl. Acad. Sci. U. S. A. 119:e2200400119. doi: 10.1073/pnas.2200400119
Borderie, A., Caclin, A., Lachaux, J. P., Perrone-Bertollotti, M., Hoyer, R. S., Kahane, P., et al. (2024). Cross-frequency coupling in cortico-hippocampal networks supports the maintenance of sequential auditory information in short-term memory. PLoS Biol. 22, 1–24. doi: 10.1371/journal.pbio.3002512
Bottomley, P. (1987). Spatial localization in nmr spectroscopy in vivo. Ann. N. Y. Acad. Sci. 508, 333–348.
Calmels, C., Foutren, M., and Stam, C. (2012). Beta functional connectivity modulation during the maintenance of motion information in working memory: importance of the familiarity of the visual context. Neuroscience 212, 49–58. doi: 10.1016/j.neuroscience.2012.03.045
Chu, T., Ji, Z., Zuo, J., Zhang, W., Huang, T., Mi, Y., et al. (2022). "Oscillatory tracking of continuous attractor neural networks account for phase precession and procession of hippocampal place cells," in Advances in Neural Information Processing Systems, Vol. 35, eds. S. Koyejo, S. Mohamed, A. Agarwal, D. Belgrave, K. Cho, and A. Oh (Red Hook, NY: Curran Associates, Inc.), 33159–33172.
Cowan, N., Saults, J., Elliott, E. M., and Moreno, M. V. (2002). Deconfounding serial recall. J. Mem. Lang. 46, 153–177. doi: 10.1006/jmla.2001.2805
Curtis, C. E., and Lee, D. (2010). Beyond working memory: the role of persistent activity in decision making. Trends Cogn. Sci. 14, 216–222. doi: 10.1016/j.tics.2010.03.006
Endel, T., and CFI, M. (2000). The Oxford Handbook of Memory. New York, NY: Oxford University Press.
Funahashi, S., Bruce, C. J., and Goldman-Rakic, P. S. (1989). Mnemonic coding of visual space in the monkey's dorsolateral prefrontal cortex. J. Neurophysiol. 61, 331–349.
Fung, C. A., Wong, K. M., Wang, H., and Wu, S. (2012). Dynamical synapses enhance neural information processing: gracefulness, accuracy, and mobility. Neural Comput. 24, 1147–1185. doi: 10.48550/arXiv.1104.0305
Fung, C. C. A., Wong, K. Y. M., and Wu, S. (2010). A moving bump in a continuous manifold: a comprehensive study of the tracking dynamics of continuous attractor neural networks. Neural Comput. 22, 752–792. doi: 10.1162/neco.2009.07-08-824
Georgopoulos, A. P., Schwartz, A. B., and Kettner, R. E. (1986). Neuronal population coding of movement direction. Science 233, 1416–1419.
Glanzer, M., and Cunitz, A. R. (1966). Two storage mechanisms in free recall. J. Verb. Learn. Verb. Behav. 5, 351–360.
Gorgoraptis, N., Catalao, R. F. G., Bays, P. M., and Husain, M. (2011). Dynamic updating of working memory resources for visual objects. J. Neurosci. 31, 8502–8511. doi: 10.1523/JNEUROSCI.0208-11.2011
Groeger, J. A., Banks, A. P., and Simpson, P. J. (2008). Serial memory for sound-specified locations: effects of spatial uncertainty and motor suppression. Quart. J. Exp. Psychol. 61, 248–262. doi: 10.1080/17470210601138746
Henry, M., and Misha, T. (1996). Redistribution of synaptic efficacy between neocortical pyramidal neurons. Nature 382, 807–810.
Henson, R. N. (2013). Short-term memory for serial order: the start-end model. Cogn. Psychol. 36, 73–137.
Huang, Q., Luo, M., Mi, Y., and Luo, H. (2023). “leader–follower” dynamic perturbation manipulates multi-item working memory in humans. eNeuro 10:2023. doi: 10.1523/ENEURO.0472-22.2023
Hurlstone, M. J., Hitch, G. J., and Baddeley, A. D. (2014). Memory for serial order across domains: an overview of the literature and directions for future research. Psychol. Bullet. 140, 339–373. doi: 10.1037/a0034221
Jensen, O., and Lisman, J. E. (2005). Hippocampal sequence-encoding driven by a cortical multi-item working memory buffer. Trends Neurosci. 28, 67–72.
Katkov, M., Romani, S., and Tsodyks, M. (2017). Memory retrieval from first principles. Neuron 94, 1027–1032. doi: 10.1016/j.neuron.2017.03.048
Kiani, R., Hanks, T. D., and Shadlen, M. N. (2008). Bounded integration in parietal cortex underlies decisions even when viewing duration is dictated by the environment. J. Neurosci. 28, 3017–3302. doi: 10.1523/JNEUROSCI.4761-07.2008
Knoedler, A. J., Hellwig, K. A., and Neath, I. (1999). The shift from recency to primacy with increasing delay. J. Exp. Psychol. 25, 474–484.
Lee, H., Choi, W., Park, Y., and Paik, S.-B. (2020). Distinct role of flexible and stable encodings in sequential working memory. Neural Netw. 121, 419–429. doi: 10.1016/j.neunet.2019.09.034
Li, J., Huang, Q., Han, Q., Mi, Y., and Luo, H. (2021). Temporally coherent perturbation of neural dynamics during retention alters human multi-item working memory. Progr. Neurobiol. 201:102023. doi: 10.1016/j.pneurobio.2021.102023
Logan, G. D. (2021). Serial order in perception, memory, and action. Psychol. Rev. 128, 1–44. doi: 10.1037/rev0000253
Lui, T. K.-Y., Obleser, J., and Wöstmann, M. (2023). Slow neural oscillations explain temporal fluctuations in distractibility. Progr. Neurobiol. 226:102458. doi: 10.1016/j.pneurobio.2023.102458
Ma, W. J., Husain, M., and Bays, P. M. (2014). Changing concepts of working memory. Nat. Neurosci. 17, 347–356. doi: 10.1038/nn.3655
Masse, N. Y., Rosen, M. C., and Freedman, D. J. (2020). Reevaluating the role of persistent neural activity in short-term memory. Trends Cogn. Sci. 24, 242–258. doi: 10.1016/j.tics.2019.12.014
Mi, Y., Fung, C. C. A., Wong, K. Y. M., and Wu, S. (2014). "Spike frequency adaptation implements anticipative tracking in continuous attractor neural networks," in Z. Ghahramani, M. Welling, C. Cortes, N. Lawrence, and K. Weinberger Advances in Neural Information Processing Systems (Red Hook, NY: Curran Associates, Inc.), 27.
Mi, Y., Katkov, M., and Tsodyks, M. (2017). Synaptic correlates of working memory capacity. Neuron 93, 323–330. doi: 10.1016/j.neuron.2016.12.004
Mongillo, G., Barak, O., and Tsodyks, M. (2008). Synaptic theory of working memory. Science 319, 1543–1546. doi: 10.1126/science.1150769
Naim, M., Katkov, M., Romani, S., and Tsodyks, M. (2020). Fundamental law of memory recall. Phys. Rev. Lett. 124:e018101. doi: 10.1103/PhysRevLett.124.018101
Pantelis, P. C., van Vugt, M. K., Sekuler, R., Wilson, H. R., and Kahana, M. J. (2008). Why are some people's names easier to learn than others? the effects of face similarity on memory for face-name associations. Mem. Cogn. 36, 1182–1195. doi: 10.3758/MC.36.6.1182
Postman, L., and Phillips, L. W. (1965). Short-term temporal changes in free recall. Quart. J. Exp. Psychol. 17, 132–138.
Potagas, C., Kasselimis, D., and Evdokimidis, I. (2011). Short-term and working memory impairments in aphasia. Neuropsychologia 49, 2874–2878. doi: 10.1016/j.neuropsychologia.2011.06.013
Rainer, G., and Miller, E. K. (2002). Timecourse of object-related neural activity in the primate prefrontal cortex during a short-term memory task. Eur. J. Neurosci. 15, 1244–1254. doi: 10.1046/j.1460-9568.2002.01958.x
Ru, X., He, K., Lyu, B., Li, D., Xu, W., Gu, W., et al. (2022). Intracranial brain-computer interface spelling using localized visual motion response. NeuroImage 258:119363. doi: 10.1016/j.neuroimage.2022.119363
Samsonovich, A., and McNaughton, B. L. (1997). Path integration and cognitive mapping in a continuous attractor neural network model. J. Neurosci. 17, 5900–5920.
Schneegans, S., and Bays, P. M. (2016). No fixed item limit in visuospatial working memory. Cortex 83, 181–193. doi: 10.1016/j.cortex.2016.07.021
Seeholzer, A., Deger, M., and Gerstner, W. (2019). Stability of working memory in continuous attractor networks under the control of short-term plasticity. PLoS Comput. Biol. 15:e1006928. doi: 10.1371/journal.pcbi.1006928
Shafi, M., Zhou, Y., Quintana, J., Chow, C., Fuster, J., and Bodner, M. (2007). Variability in neuronal activity in primate cortex during working memory tasks. Neuroscience 146, 1082–1108. doi: 10.1016/j.neuroscience.2006.12.072
Stokes, M. G. (2015). "Activity-silent" working memory in prefrontal cortex: a dynamic coding framework. Trends Cogn. Sci. 19, 394–405. doi: 10.1016/j.tics.2015.05.004
Storm, B. C., and Bjork, R. A. (2016). Do learners predict a shift from recency to primacy with delay? Mem. Cogn. 44, 1204–1214. doi: 10.3758/s13421-016-0632-9
Stringer, S.M., Trappenberg, T. P. E. R., and Araujo, I. (2002). Self-organizing continuous attractor networks and path integration: one-dimensional models of head direction cells. Network 13, 217–242. doi: 10.1080/net.13.2.217.242
Taube, J., Muller, R., and JB Ranck, J. (1990). Head-direction cells recorded from the postsubiculum in freely moving rats. I. Description and quantitative analysis. J. Neurosci. 10, 420–435.
Teyler, T. J., and Discenna, P. (1984). Long-term potentiation as a candidate mnemonic device. Brain Res. Rev. 7, 15–28.
Tsetsos, K., Gao, J., McClelland, J. L., and Usher, M. (2012). Using time-varying evidence to test models of decision dynamics: bounded diffusion vs. the leaky competing accumulator model. Front. Neurosci. 6:79. doi: 10.3389/fnins.2012.00079
Wang, C., Deng, H., Dong, Y., Zhang, X., and Wang, D. H. (2021). The capacity and cognitive processing of vibrotactile working memory for frequency. Curr. Psychol. 2021, 1–11. doi: 10.1007/s12144-021-02212-6
Wang, R., and Kang, L. (2022). Multiple bumps can enhance robustness to noise in continuous attractor networks. PLoS Comput. Biol. 18, 1–38. doi: 10.1371/journal.pcbi.1010547
Wang, Y., Markram, H., Goodman, P. H., Berger, T. K., Ma, J., and Goldman-Rakic, P. S. (2006). Heterogeneity in the pyramidal network of the medial prefrontal cortex. Nat. Neurosci. 9, 534–542. doi: 10.1038/nn1670
Wolff, M. J., Ding, J., Myers, N. E., and Stokes, M. G. (2015). Revealing hidden states in visual working memory using electroencephalography. Front. Syst. Neurosci. 9:123. doi: 10.3389/fnsys.2015.00123
Wolff, M. J., Jochim, J., Akyürek, E. G., and Stokes, M. G. (2017). Dynamic hidden states underlying working-memory-guided behavior. Nat. Neurosci. 20, 864–871. doi: 10.1038/nn.4546
Wu, S., Hamaguchi, K., and Amari, S. I. (2013). Dynamics and computation of continuous attractors. Neural Comput. 20, 994–1025. doi: 10.1162/neco.2008.10-06-378
Wu, S., Wong, K., Fung, C., Mi, Y., and Zhang, W. (2016). Continuous attractor neural networks: candidate of a canonical model for neural information representation [version 1; peer review: 2 approved]. F1000Research 5:1. doi: 10.12688/f1000research.7387.1
Yoon, K., Buice, M. A., Barry, C., Hayman, R., Burgess, N., and Fiete, I. R. (2013). Specific evidence of low-dimensional continuous attractor dynamics in grid cells. Nat. Neurosci. 16, 1077–1084. doi: 10.1038/nn.3450
Keywords: sequential working memory, serial position effect, short-term plasticity effect, continuous attractor neural networks, the primacy and recency effect
Citation: Zhou J, Gong L, Huang X, Mu C and Mi Y (2024) The synaptic correlates of serial position effects in sequential working memory. Front. Comput. Neurosci. 18:1430244. doi: 10.3389/fncom.2024.1430244
Received: 09 May 2024; Accepted: 26 June 2024;
Published: 15 July 2024.
Edited by:
Yu-Guo Yu, Fudan University, ChinaReviewed by:
Daqing Guo, University of Electronic Science and Technology of China, ChinaDa-Hui Wang, Beijing Normal University, China
Copyright © 2024 Zhou, Gong, Huang, Mu and Mi. 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: Yuanyuan Mi, bWl5dWFueXVhbiYjeDAwMDQwO3RzaW5naHVhLmVkdS5jbg==; Chunlai Mu, Y2xtdTIwMDUmI3gwMDA0MDsxNjMuY29t
†These authors share first authorship