Europe PMC

This website requires cookies, and the limited processing of your personal data in order to function. By using the site you are agreeing to this as outlined in our privacy notice and cookie policy.

Abstract 


Persistently active neurons during mnemonic periods have been regarded as the mechanism underlying working memory maintenance. Alternatively, neuronal networks could instead store memories in fast synaptic changes, thus avoiding the biological cost of maintaining an active code through persistent neuronal firing. Such "activity-silent" codes have been proposed for specific conditions in which memories are maintained in a nonprioritized state, as for unattended but still relevant short-term memories. A hallmark of this "activity-silent" code is that these memories can be reactivated from silent, synaptic traces. Evidence for "activity-silent" working memory storage has come from human electroencephalography (EEG), in particular from the emergence of decodability (EEG reactivations) induced by visual impulses (termed pinging) during otherwise "silent" periods. Here, we reanalyze EEG data from such pinging studies. We find that the originally reported absence of memory decoding reflects weak statistical power, as decoding is possible based on more powered analyses or reanalysis using alpha power instead of raw voltage. This reveals that visual pinging EEG "reactivations" occur in the presence of an electrically active, not silent, code for unattended memories in these data. This crucial change in the evidence provided by this dataset prompts a reinterpretation of the mechanisms of EEG reactivations. We provide 2 possible explanations backed by computational models, and we discuss the relationship with TMS-induced EEG reactivations.

Free full text 


Logo of plosbiolLink to Publisher's site
PLoS Biol. 2021 Oct; 19(10): e3001436.
Published online 2021 Oct 21. https://doi.org/10.1371/journal.pbio.3001436
PMCID: PMC8641864
PMID: 34673775

Pinging the brain with visual impulses reveals electrically active, not activity-silent, working memories

Joao Barbosa, Conceptualization, Formal analysis, Funding acquisition, Investigation, Methodology, Software, Writing – original draft, Writing – review & editing, 1 , 2 Diego Lozano-Soldevilla, Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Software, 1 , 3 , 4 and Albert Compte, Conceptualization, Funding acquisition, Investigation, Methodology, Supervision, Writing – review & editing 1 ,*
Frank Tong, Academic Editor

Associated Data

Supplementary Materials
Data Availability Statement

Abstract

Persistently active neurons during mnemonic periods have been regarded as the mechanism underlying working memory maintenance. Alternatively, neuronal networks could instead store memories in fast synaptic changes, thus avoiding the biological cost of maintaining an active code through persistent neuronal firing. Such “activity-silent” codes have been proposed for specific conditions in which memories are maintained in a nonprioritized state, as for unattended but still relevant short-term memories. A hallmark of this “activity-silent” code is that these memories can be reactivated from silent, synaptic traces. Evidence for “activity-silent” working memory storage has come from human electroencephalography (EEG), in particular from the emergence of decodability (EEG reactivations) induced by visual impulses (termed pinging) during otherwise “silent” periods. Here, we reanalyze EEG data from such pinging studies. We find that the originally reported absence of memory decoding reflects weak statistical power, as decoding is possible based on more powered analyses or reanalysis using alpha power instead of raw voltage. This reveals that visual pinging EEG “reactivations” occur in the presence of an electrically active, not silent, code for unattended memories in these data. This crucial change in the evidence provided by this dataset prompts a reinterpretation of the mechanisms of EEG reactivations. We provide 2 possible explanations backed by computational models, and we discuss the relationship with TMS-induced EEG reactivations.

Introduction

A hallmark of the activity-silent working memory framework [1] is that memories stored silently in synaptic traces through short-term synaptic plasticity can be reactivated through nonspecific stimuli [15]. Evidence supporting activity-silent working memory has recently emerged from human electroencephalography (EEG) [6,7], in particular from EEG reactivations of unattended memories induced by visual impulses [7]—the so-called visual pinging. Despite their relevance for upcoming memory-guided behavior, currently unattended memories could not be robustly decoded from raw EEG voltage traces [6,7] (Fig 1A, red). In view of this, unattended memories resemble memories rendered behaviorally irrelevant by a contextual cue (discarded, Fig 1B, dashed lines), but they differ from attended memories with similar upcoming behavioral requirements, which are represented in sustained, active codes [68] (Fig 1A and 1B, solid lines). This observation has been key in interpreting EEG reactivations in pinging studies as evidence for activity-silent storage (see, e.g., recent reviews [814] or [24,15] for explicit simulations of this interpretation of the data). Intriguingly, pinging-induced increase in EEG decodability occurred exclusively for items that remained relevant for future, memory-guided behavior, suggesting that only unattended but still relevant items were kept in activity-silent traces. The mechanisms for such selective reactivation of activity-silent traces are unclear, as in existing computational models of activity-silent storage [1,2,4,5,16,17] short-term plasticity changes are induced by neuronal activity, regardless of its behavioral relevance. Here, we reanalyzed EEG recordings of these influential pinging studies [7,18]. We found that unattended memories, previously shown to be inaccessible from scalp EEG voltages despite remaining behaviorally relevant [7], could in fact be decoded from raw EEG voltage and robustly from alpha power signals. This reanalysis demonstrates that the original claim of a silent representation of unattended memories is not supported by the data, which instead show an active code and thus calls for a reinterpretation of pinging-induced increases in EEG decodability. Finally, we argue that the increase in stimulus decodability following an unspecific stimulus, seen in human [6,7,1820] and monkey electrophysiological experiments [21,22], can be explained by network models without short-term plasticity based on ongoing active, not silent, neural representations.

An external file that holds a picture, illustration, etc.
Object name is pbio.3001436.g001.jpg
Decoding from alpha power reveals an underlying, active working memory code.

(A) Strength of stimulus decoding from raw voltage traces as in [7]. As in the original study, unattended or (B) discarded memories cannot be decoded from raw voltage traces. (C, D) Same as (A) and (B), but decoding from alpha power (Methods), which reveals a sustained representation of the unattended stimulus. In (A) and (C), we analyze data from experiment 2 [7], while in (B) and (D), data from experiment 1 [7]. Light gray bars mark stimulus presentation periods. Notice that data immediately preceding pinging stimulus presentation are shown in this figure. Dashed lines mark the periods in which memories are irrelevant for upcoming behavior, following an instruction cue (dark green). All error bars are bootstrapped SEM, and color bars on the top mark the periods where bootstrapped 95% CI was above zero. Data from Wolff and colleagues (2017) [7].

Results

Attended and unattended working memories are robustly decoded from alpha power

We realized that an earlier study had reported that (attended) spatial memories were decoded more reliably from EEG total alpha power than from evoked activity [23]. Thus, we analyzed alpha-power information content for attended, unattended, and no longer relevant orientation memories in the publicly available dataset of the original publication by Wolff and colleagues [7]. We found that a sustained alpha power code tracks the orientation of the items that remain relevant for future behavior, whether attended or unattended (Fig 1C and 1D, solid lines, S1 Fig). This shows that working memory contents are maintained in an electrically active neural code, even for items outside of the current attentional focus [24]. However, while attended memories were decodable both in alpha power and voltage traces, unattended memories were only robustly detected in alpha power. While this could reflect qualitative differences in what these 2 neural signals represent [2527], we will explore here the parsimonious possibility that this stems from differential sensitivity of these 2 measures to the same underlying neural activity. Indeed, EEG voltage is known to lose decodability shortly following a reference baselining due to slow electrical drifts [28], while oscillations could be more robust to these baseline drifts. In this view, attended items would be represented by strong neural signals (represented both in voltage and in alpha power), while unattended items would be kept in analogous but weaker neural signals (picked only by alpha power). These neural dynamics could result from competitive interactions between prioritized and unprioritized memory items, in line with competing attractors in networks without activity-silent mechanisms [20,24,29].

Lack of statistical power suggests spurious evidence for silent representations of unattended memories

In line with the hypothesis of an active but weaker representation of the unattended items, recent studies [30,31] show that lack of decodability for unattended working memories can be overcome by increasing statistical power (e.g., sample size). We wondered if a similar strategy could improve decodability of unattended orientations from raw voltages in this dataset. We addressed this tentatively in the current datasets with 2 complementary approaches. First, we attempted to increase the statistical power by (1) smoothing the voltage traces with a 32-ms kernel prior to decoding (instead of smoothing instantaneous decoding accuracies as done in the original study [7]; Methods); (2) averaging decoding accuracies over an interval of 200 ms before the pinging impulse; and (3) pooling trials from all sessions and subjects (Methods). We found that discarded memories were still not decodable from raw voltage (p = 0.48, t = −0.04, one-sided t test), but the decoding of unattended memories almost reached the typical statistical threshold (p = 0.06, t = 1.55, one-sided t test). We also varied the number of sessions and trials included in the analyses to visualize how decoding depended on sample size (Fig 2). This further suggested that our analyses were underpowered and motivated our second approach, in which we reasoned that the low signal-to-noise ratio in this cohort could be due to specific sessions with overall low decodability. Sessions with low decodability could reflect technical issues during that particular session (e.g., EEG sensor placement) or specific subject characteristics, such as skull thickness or hair density. We thus divided our full dataset using cross-validation in high and low decoding sessions, based on the average decoding accuracy during the early delay (“split period” in Fig 3, Methods). We found that unattended memories could be robustly decoded during the whole delay (0.25–1.2 s, p = 0.002 randomization test, Methods) and in particular immediately before pinging (250 ms window, p = 0.039, randomization test, Methods) from high-decoding sessions, while discarded memories could not (both p>0.45, Fig 3). Finally, we show that sessions with high early-delay decoding (“split period” in Fig 3) are also those that have “reactivations” in voltage seen in the original publication (Fig 3, inset). Note that we used one-sided statistical tests (Figs (Figs22 and and3),3), since negative decoding strengths are not expected. Additionally, one-sided statistical tests represent a conservative approach when claiming lack of decoding. These analyses, together with previous studies showing robust decoding of unattended memory items [22,30,3236], suggest that also in this dataset, unattended items are not stored in activity-silent traces.

An external file that holds a picture, illustration, etc.
Object name is pbio.3001436.g002.jpg
Average voltage decoding during 200 ms prior to the impulse as a function of trial and session numbers [37].

Attended items could be robustly decoded from voltage, but not discarded items. Decoding of unattended items suggested a possible underlying signal (p = 0.06, t = 1.55, one-sided t test, upper-right corner in second panel). Subsampling of sessions and trials was done by randomly subsampling (n = 5,000) without repetition from the full dataset (Methods). In white, contour lines for different p-value levels (one-sided t test). Data from Wolff and colleagues (2017) [7].

An external file that holds a picture, illustration, etc.
Object name is pbio.3001436.g003.jpg

Sessions with high early-delay (split period, Methods) voltage decoding have a sustained code for unattended memories (left, red), but not for discarded memories (right). Replotting the reactivation period (inset), separately for high and low early-delay decoding sessions shows that “reactivations” only occur for sessions with a sustained code (left, red). Note that at time 0, the decoding strength is not actually zero (inset vs arrow). This is an artifact of baselining during the delay period (see also S3 Fig). Error bars are sem. Decoding strengths from high-decoding sessions were compared to the shuffle predictor (top black bars mark significant deviation, one-sided p<0.05, Methods). Time course and data are similar to Fig 1A and 1B. Data from Wolff and colleagues (2017) [7].

Two plausible explanations for the increase in decodability that do not require activity-silent mechanisms

If there is an active EEG code for both attended and unattended stimuli prior to the visual impulse, as our analyses suggest, then what is the interpretation of the observed increase in EEG decodability [7,18]? We reasoned that EEG reactivation events may emerge from either an increase in the signal about the stimulus (as assumed in the activity-silent interpretation) or through a reduction in the across-trial variability (S3C Fig). In the data, we found that pinging reduces across-trial variability of EEG voltage (Fig 4A), as expected for neural responses to sensory stimuli [38]. In addition, we found that trials with stronger EEG decodability showed lower across-trial variability than trials with weaker EEG decodability during pinging (Fig 4B), demonstrating a link between trial-by-trial EEG variability and pinging-induced increase in EEG decodability. We argue that a reduction of variability with an otherwise intact active memory representation (Fig 3) is a parsimonious interpretation of the visual pinging effect. Alternatively, there is another interpretation of pinging-induced increases in EEG decodability consistent with our findings. Recent modeling work has shown how the enhancement of active representations [20,21] is expected when pinging recurrent neural networks with no need for activity-silent mechanisms [29,39]. An existing representation maintained in an attractor supported by recurrent and competitive interactions enhances its tuning when it is stimulated unspecifically (attractor-boost model, S2 Fig). Also, this mechanism would be consistent with these data, as it shows reduced variability (Fig 4C, gray), concomitant with boosted attractor tuning (S2 Fig). While both of these possible interpretations do not exclude an interplay of active representations with activity-silent mechanisms [5,17,40,41], they offer a parsimonious view that renders activity-silent working memory an inadequate framework to understand increases in decodability induced by nonspecific stimuli [7,18,20,21]. To further support this, we sought to evaluate variability predictions from a computational model where reactivations occur because of factual memory reactivation from silent, synaptic traces. We tested an available biophysical network model for (continuous) activity-silent working memory [5], which is an extension of the canonical but discrete model of activity-silent working memory [1] (Methods). In these simulations, a nonspecific input induced reactivations in some trials, causing an increase of across-trial variability (Fig 4C, black). This is because reactivations in such attractor networks are an all-or-none phenomenon, and great variability is expected when triggering them from weak, decaying activity-silent traces in noisy spiking networks. In sum, pinging reveals an underlying active memory, perhaps by reducing noise (Figs 4A, 4B, and S3C) in the presence of an active code (Figs (Figs11 and and3)3) or by enhancing tuning in an active representation (S2 Fig), but not by reactivating stimulus signals from silent traces.

An external file that holds a picture, illustration, etc.
Object name is pbio.3001436.g004.jpg
Impulse-induced, across-trial variability change in human EEG and computational models.

(A) Percentage of variability change, relative to 0.2 s before the impulse, computed when the impulse was a visual stimulus in Wolff and colleagues (2017) [7] in gray and Wolff and colleagues (2015) [18] in black. (B) Difference in variances computed across trials with low vs high stimulus decoding (computed at the time of maximal decodability, black arrow) in Wolff and colleagues (2015) [18]. Trials with strong memory decoding showed significantly lower across-trial variance than trials with weak or absent memory decoding. We did not find a significant correlation in Wolff and colleagues (2017) [7] possibly because of a weaker pinging stimulus, which may have contributed to weaker increase in EEG decodability, not visible without baselining the data during the delay period (see S3 Fig). (C) Simulations of the activity-silent working memory model with short-term plasticity (dark) predict an increase of across-trial variability (Fano factor) following reactivations induced by a nonspecific drive. Simulations of the bump-attractor model without short-term plasticity (light) predict a decrease in the Fano factor following a nonspecific drive. See also S3 Fig. (D) Same as (A) but when the impulse was a single-pulse TMS (data from Rose and colleagues, 2016 [6]). Solid bars mark where the change in variability was significant (two-sided t test, p < 0.005), and error bars are bootstrapped 95% CI of the mean. Six out of 54 sessions had outlier TMS artifact variance (i.e., extremely high variance at the time of the impulse) and were removed from this analysis. EEG, electroencephalography; TMS, transcranial magnetic stimulation.

Differential mechanisms for pinging- and TMS-induced reactivations

Increases of EEG decodability at the presentation of nonspecific impulses have been shown to occur not only for visual impulses [7] but also for external perturbation with single-pulse transcranial magnetic stimulation (TMS) [6]. Despite apparent similarities, TMS perturbations impact working memory performance [5,6,4244], while pinging does not [7,45] (see also chapter 4 of [46]). This suggested that these approaches could be interacting with fundamentally different neural mechanisms [7,47]. Indeed, we found through the reanalysis of the data of [6] that single-pulse TMS increases across-trial EEG variability (Fig 4D), in contrast to the reduction observed upon visual pinging (Fig 4A). Such increase in across-trial variability is in accordance with the activity-silent working memory model presented before (Fig 4C, black), thus potentially supporting the interpretation of TMS EEG reactivations as signals recovered from activity-silent traces [5,6]. However, a note of caution is in order: The difficulty in precisely locating the TMS coil in different trials may contribute to increased EEG variability by virtue of the long-lasting effect of the TMS pulse on neural excitability [48]. This could mask EEG signals reflecting TMS-induced reactivations.

Discussion

Through the reanalysis of existing datasets, we provide here evidence that working memory “reactivations” by visual pinging, considered prime evidence for activity-silent working memory [24,815], occur in the presence of active (not silent) ongoing memory representations in the delay period. In addition, we verify that representations for unattended items are notably weaker than for attended items, consistent with biased competition between active memories [29]. Based on this substrate for working memory, visual pinging may increase EEG decodability through (1) ping-induced reduction in across-trial variability; and/or (2) ping-induced boosting of attractor tuning. We further compare visual pinging with TMS perturbations, and we find qualitative differences suggesting different underlying mechanisms. Based on the difference in behavioral impact of these 2 perturbation protocols (visual pinging does not affect working memory behavior, but TMS does), we speculate that visual pinging may increase EEG decodability via reduced across-trial variability or by transient boosting of active attractors, while TMS-induced reactivations would be supported by activity-silent mechanisms. Note that temporarily boosting an active attractor should not have a strong impact on behavior beyond the boosting period (unless additional long-lasting cellular mechanisms are engaged), while true reactivations from activity-silent stores should have a long-lasting impact, as the silent trace is refreshed.

We show here that relevant memories are stored in active codes and are thus decodable, while irrelevant memories could not be decoded, therefore potentially discarded. This interpretation sheds light on the intriguing ineffectiveness of pinging on discarded memories, while being effective in “reactivating” unattended but still relevant ones [7]. The reactivation effectiveness of pinging appears now to depend on whether memories are maintained in active neuronal representations, decodable from EEG. In line with this, a recent study shows that pinging reveals diffusing dynamics [45], a hallmark of active memories [45,49], instead of decaying dynamics, as expected in the activity-silent framework [1,45]. Our data reconcile the influential study by Wolff and colleagues [7] with recent works showing irrelevant or unattended memories actively encoded in scalp EEG [3234], in the activity of cortical association areas using large-sample fMRI analyses [30] or intracranial recordings in monkeys [22,50], and in neural activity in visual areas of rodents [36]. While there is extensive evidence for long-lasting cellular and synaptic mechanisms in cortical neurons (“silent” mechanisms, e.g., [51]) that must coexist [52] with, or even support [17,40,5355], active representations such as those reported in this study, there is more scant evidence that working memory can be volitionally stored without spiking activity [1,56].

Explicit evidence for activity-silent processes is difficult to obtain but is particularly confounded in the presence of active representations. A possible approach is to seek evidence for activity-silent traces from previously memorized but already irrelevant items [5,57], for which chance-decoding is expected in principle. There, too, it is hard to discard low-powered designs, so positive evidence must be sought. Recently, it has been reported that, between consecutive trials—when the previous memory should be discarded, similarly to uncued memories in [7]—neurons fire more synchronously after having been engaged in active working memory storage [5], suggesting that discarded memories can leave involuntary silent traces. Importantly, selecting sessions based on good overall decodability (as in Fig 3) did not reveal decodability of these across-trial discarded memories (Fig 1 in [5]), supporting a true activity-silent substrate. We argue that the effect of induced reactivations must be validated during similar conditions in which memories are demonstrably discarded, presumably leaving an activity-silent trace. Current evidence shows across-trial behavioral effects of TMS, as is expected by reactivating previous, discarded memories [5,44], but not of visual pinging during intertrial periods (chapter 4 of [46]). Also, the evidence in Fig 3 suggests that visual pinging cannot reactivate putative activity-silent traces following a discarded item.

Finally, our results suggest that voltage and alpha power encode similar working memory content. Previous work, however, shows that alpha and raw voltage play different roles in working memory [2527], in particular that alpha power tracks spatial attention instead of the actual memory content. In principle, with the experimental design of Wolff and colleagues [7], subjects could recode orientation as a spatial location, which would be tracked by alpha power and by sustained voltages [25]. Future work with other experimental designs including independent variation of orientation and attentional location [25] could clarify this point further. Regardless of whether EEG is tracking orientation or a location recodification, we argue that in these data, both signals are carrying analogous contents, with voltage being noisier.

In sum, our results add to previous literature showing robust decoding of unattended working memories from electrophysiological signals [22,30,3236,50]. Our analyses reinforce the idea that interpreting null decoding as evidence for storage in silent traces is not straightforward, because null results might result from weak signals in insufficiently powered analyses [30,31].

Methods

EEG experiments

We analyzed 2 available datasets of visual pinging [7,18]. For decoding and EEG variability analyses, we focused on both experiments from [7]. In experiment 1 (n = 30), subjects were cued for which item was going to be probed (cued item, here also called attended, or uncued item, here called discarded memory). In experiment 2 (n = 19), subjects had to alternate their attention between 2 items (their early/late, here attended/unattended memory item). Experiment 1 consisted of 1 session, while experiment 2 consisted of 2 sessions (separated by approximately 1 to 2 weeks) on the same set of subjects. For variability analyses, we also analyzed the experiment of [18] (n = 24). In this experiment, the subjects had to memorize 1 item, thus always within the focus of attention. Importantly, the item decodability from raw voltage never dropped to chance. Additionally, we also analyze the voltage variance of the experiment 2 (n = 6) of a TMS study [6]. We refer the reader to the original studies for extra details [6,7,18]. All these datasets were made available in a fully anonymized format and had been approved by the corresponding institutional review boards, as indicated in the original publications.

Data preprocessing

The data available online was epoched and baselined relative to the beginning of each epoch (S3 Fig). To revert this baselining, we computed trial-by-trial voltage difference between consecutive epochs. We then added this voltage difference to the beginning of each baselined epoch, effectively reverting all baselining effects. Additionally, for the variability analyses (see below), we remove any signal drift caused by, e.g., moving electrodes using the python function scipy.signal.detrend on each subject’s variability. Finally, for the decoding analysis (see below), we also computed the alpha power. For this purpose, the data were Hilbert-transformed (using the FieldTrip function “ft_freqanalysis.m”) to extract frequencies in the alpha-band (8 to 12 Hz), and total power was calculated as the squared complex magnitude of the signal.

Decoding analyses

We used freely available code to perform these analyses, so we will only briefly describe the methodology here. For a detailed description of the decoding methods, please refer to the original study [7]. As in the original study, we decoded from all the 17 posterior channels (P7, P5, P3, P1, PO7, PO3, and O1 versus P8, P5, P6, P4, P2, PO8, PO4, and O2). Briefly, we collected the Mahalanobis distance between all possible pairwise combinations of the orientations and thus form a representational dissimilarity function. Finally, the decoding strength was calculated as the vector strength of this function. We decoded from raw voltage or alpha power with the exact same code. The decoding strengths (Fig 1) were smoothed over time with a gaussian kernel (SD = 10 ms).

Across-trial variability analyses

We computed variability as the variance (var) across trials of the raw voltage traces (trials × sensors × time). Before averaging variances across sensors, we detrended them using the function scipy.signal.detrend to account for any drift in the signal. Finally, we computed the percentage of variability change (Δvar) relative to the baseline period of 2 s before the pinging stimulus (b): Δvar = (var − b) / b * 100. This referencing to the baseline ensures that changes in variability can be attributed solely to the pinging and not to other factors that are common to both pre and after pinging, such as varying stimulus orientations. We computed across-trial variability for each session separately, and then averaged variances across sessions. This is important for the TMS experiment, in which TMS location was held fixed. This way, our variance analysis is not capturing aspects of the design that vary from session to session (e.g., spatial attention, TMS target location).

Fano factor

To compute the variability drop in the simulated spiking activity, we used the Fano factor [38], which is defined as the variance of spike counts in a given window (100 ms) divided by their mean. We then computed ΔFF, as the difference relative to the baseline period of 2 s before pinging stimulus.

Phenomenological simulations of EEG trials

To study how single-trial baseline correction impacts pinging-induced increases in EEG decodability, we applied our decoders to synthetic EEG data generated by a model where spurious EEG reactivations are caused by a reduction in noise variability. First, we simulated 2 hypothetical delay maintenance EEG time series representing 2 independent experimental conditions (i.e., grating oriented 0 0 versus 45 0; n = 200 each condition) using the following Gamma function (f) as a single-trial waveform:

f(t|a,b)=ea1[b(a1)]a1ta1et/b

For each trial and condition, parameters a and b were drawn from a gaussian distribution (condition 1, μa = 2, μb = 130 ms; condition 2, μa = 3, μb = 80 ms; all conditions with σa = 0.2, σb = 0.5 ms). Each waveform was then scaled by 0.5 and 0.25, respectively, and the time onset was set to time point 0.1 s.

Finally, for each single-trial and condition, gaussian noise was added (μ = 0, σ = 0.75). During the impulse period (1.3 s to 1.5 s), for each trial and experimental condition, we reduced 30% the variability of the noise (σ = 0.525).

To use the same decoding method as the original publication, we generated a multichannel time series. We created 2 sets of dipoles located around left (position [1.5–8.6 1.5] cm, orientation [–111]) and right (position [1.5–8.6 1.5] cm, orientation [111]) primary visual cortex, and the simulated time courses were projected to the scalp via a forward model [58]. S3C Fig plots decoding of these signals upon different conditions of EEG voltage baselining. The Matlab code for these simulations is available on https://github.com/comptelab/reactivations.

Activity-silent network model

We used a previously proposed computational model [5] to simulate memory reactivations. The model consists of a network of interconnected 2,048 excitatory and 512 inhibitory leaky integrate-and-fire neurons [59]. This network was organized according to a ring structure: Excitatory and inhibitory neurons were spatially distributed on a ring so that nearby neurons encoded nearby spatial locations [60]. Excitatory connections were all to all and spatially tuned, so that nearby neurons with similar preferred directions had stronger than average connections, while distant neurons had weaker connections. All inhibitory connections were all to all and untuned. Network parameters were taken from [5]. Simulation of “activity-silent” mechanisms was done by simulating 2 presynaptic variables x and u, as described in [5]. Reactivations were accomplished stimulating all excitatory neurons with a nonspecific external stimulus [5].

Attractor-boost network model

For the attractor-boost model, we used a bump attractor model similar to the activity-silent model described above, but without short-term plasticity. As in this other model, attractor boosting was achieved with a nonspecific external stimulus to all excitatory neurons. The Brian [61] code for this model is available on https://github.com/comptelab/reactivations.

Improving statistical power

Smoothing

To improve signal-to-noise ratio for our decoding analyses, we smoothed the voltage traces using a gaussian kernel with σ = 32 ms. This was in contrast with the original study that used instantaneous, nonsmoothed voltages for decoding and smoothed the resulting decoding accuracies with σ = 16 ms.

Pooling all trials, across different sessions

We also pooled all trials across sessions and subjects. Because all subjects and sessions consisted of a similar number of trials, we are not biasing our analyses toward a specific subject. Note that we only pooled trials across sessions for the analyses in Fig 2. We averaged across sessions for the other analyses.

Cross-validated median split

To simulate an increase of signal-to-noise ratio, we removed sessions with low decodability. Importantly, to avoid circularity in our analysis, we cross-validated this selection in the following way. For each session, we randomly split the single-trial decoding accuracies in two-halves. With the first half, we sorted the sessions by their average decoding accuracy during early delay [0.25 to 0.45] s, and we selected high and low decoding sessions (median split). We then computed the average decoding accuracy in the second half for low- and high-decoding sessions defined in the first half. We repeated this procedure 2,000 times for different random half-splits of the data (split-folds). Finally, we established the chance-level decoding accuracy (shuffle predictor) for high-decoding sessions by averaging the decoding strengths of 1,000 shuffled decoders (i.e. permuting orientation labels) analyzed with the same procedure (200 split-folds). For high-decoding sessions, we considered the average decoding strength significant if at least 95% of the split-folds were higher than the shuffle predictor.

Supporting information

S1 Fig

A sustained alpha-power code tracks behaviorally relevant orientation memories.

Data from Wolff and colleagues (2017) [7].

(TIF)

S2 Fig

In an attractor-boost network model, a nonspecific stimulus increases tuning without reactivation.

Two example stimulations of a bump-attractor with (weak drive) and without (no drive) a nonspecific drive at the end of the trial. Importantly, we did not include short-term plasticity in either simulation; thus, reactivations are not possible. Right, tuning during the last 0.5 s is higher for the trials in which a nonspecific drive was delivered (green), compared to when no drive was delivered (black). Similar models have been put forward in previous publications [29] (see also [39]).

(TIF)

S3 Fig

Effect of EEG baselining on stimulus decoding analyses.

(A) Original data from Wolff and colleagues (2017) [7] was baselined twice. (B) For our analyses, we de-baselined the second baselining of each trial so we could get continuity in EEG voltage traces through the whole trial (early baseline). Importantly, the exact time of the baseline affects the strength of “reactivations.” Note that data with early baselining do not show any visible increase in EEG decoding. (C) Top: diagram outlining the computer simulation that generated the EEG synthetic data, as if 2 current dipoles were placed within the visual cortex (yellow) to recreate artificial trials. Two different sets of trials (n = 200, blue and red), corresponding to 2 different stimuli (e.g., 0°, 45°) were generated (Methods). We simulated 2 event-related potentials (electrode “Oz”; white) followed by a drop in across-trial variability. Bottom: illustration of how a similar signal-to-noise ratio increase (shown in “decoding” for the drop in variance case) can result from a drop in variance or actual reactivation (signal increase). (D) Through simulations, we show that the baselining procedure introduces spurious reactivations in data without any true reactivation signal. By applying the decoding methods of Wolff and colleagues [7], we observe that (as in the data) spurious reactivations are barely visible with a distant baseline (early baseline) but are amplified for more proximal baselines (late baseline). These analyses illustrate that baselining is problematic and should be avoided, as it has been previously pointed out [28]. Darker lines mark the impulse period in the data and drop in variance in the simulations. Data from Wolff and colleagues (2017) [7].

(TIF)

Acknowledgments

We thank Nathan Rose for sharing part of Rose and colleagues (2016) [6] dataset privately and the authors of Wolff and colleagues (2015, 2017) [7,18] for sharing the full dataset publicly. We thank Michael Wolff and Mark Stokes for constructive discussions.

Funding Statement

This work was funded by the Spanish Ministry of Science, Innovation and Universities and European Regional Development Fund (https://www.ciencia.gob.es/, Refs: BFU2015-65315-R and RTI2018-094190-B-I00) to JB, DLS, AC; by the Institute Carlos III, Spain (https://eng.isciii.es/eng.isciii.es/Paginas/Inicio.html, grant PIE 16/00014) to AC, DLS and grant AC20/00071 to AC; by the Cellex Foundation to DLS, AC; by the Generalitat de Catalunya (AGAUR, https://agaur.gencat.cat/en/inici/index.html, 2014SGR1265, 2017SGR01565) to JB, DLS, AC; and by the CERCA Programme/Generalitat de Catalunya (https://cerca.cat/en/) to AC. JB was supported by the Spanish Ministry of Economy and Competitiveness (FPI program) and by the Bial Foundation (https://www.bial.com/com/bial-foundation, Ref: 356/18). This work was developed at the building Centro Esther Koplowitz, Barcelona. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Data Availability

All the code used for the analyses can be found at https://github.com/comptelab/reactivations. The data and data decoder scripts are available from the original publications (Wolff et al., 2015 [18]; Wolff et al., 2017 [7]; Rose et al. 2016 [6]).

References

1. Mongillo G, Barak O, Tsodyks M. Synaptic theory of working memory. Science. 2008;319:1543–6. 10.1126/science.1150769 [Abstract] [CrossRef] [Google Scholar]
2. Masse NY, Yang GR, Song HF, Wang X-J, Freedman DJ. Circuit mechanisms for the maintenance and manipulation of information in working memory. Nat Neurosci. 2019;22:1159–67. 10.1038/s41593-019-0414-3 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
3. Pals M, Stewart TC, Akyurek EG, Borst JP. A functional spiking-neuron model of activity-silent working memory in humans based on calcium-mediated short-term synaptic plasticity. PLoS Comput Biol. 2020;16:e1007936. 10.1371/journal.pcbi.1007936 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
4. Manohar SG, Zokaei N, Fallon SJ, Vogels TP, Husain M. Neural mechanisms of attending to items in working memory. Neurosci Biobehav Rev. 2019;101:1–12. 10.1016/j.neubiorev.2019.03.017 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
5. Barbosa J, Stein H, Martinez RL, Galan-Gadea A, Li S, Dalmau J, et al.. Interplay between persistent activity and activity-silent dynamics in the prefrontal cortex underlies serial biases in working memory. Nat Neurosci. 2020;23:1016–24. 10.1038/s41593-020-0644-4 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
6. Rose NS, LaRocque JJ, Riggall AC, Gosseries O, Starrett MJ, Meyering EE, et al.. Reactivation of latent working memories with transcranial magnetic stimulation. Science. 2016;354:1136–9. 10.1126/science.aah7011 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
7. Wolff MJ, Jochim J, Akyürek EG, Stokes MG. Dynamic hidden states underlying working-memory-guided behavior. Nat Neurosci. 2017;20:864–71. 10.1038/nn.4546 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
8. Masse NY, Rosen MC, Freedman DJ. Reevaluating the Role of Persistent Neural Activity in Short-Term Memory. Trends Cogn Sci (Regul Ed). 2020;24:242–58. 10.1016/j.tics.2019.12.014 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
9. Beukers AO, Buschman TJ, Cohen JD, Norman KA. Is activity silent working memory simply episodic memory? Trends Cogn Sci (Regul Ed). 2021;25:284–93. 10.1016/j.tics.2021.01.003 [Abstract] [CrossRef] [Google Scholar]
10. Oberauer K. Working Memory and Attention—A Conceptual Analysis and Review. J Cogn. 2019;2:36. 10.5334/joc.58 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
11. Scimeca JM, Kiyonaga A, D’Esposito M. Reaffirming the sensory recruitment account of working memory. Trends Cogn Sci (Regul Ed). 2018;22:190–2. 10.1016/j.tics.2017.12.007 [Abstract] [CrossRef] [Google Scholar]
12. Sreenivasan KK, D’Esposito M. The what, where and how of delay activity. Nat Rev Neurosci. 2019;20:466–81. 10.1038/s41583-019-0176-7 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
13. Mallett R, Lewis-Peacock JA. Behavioral decoding of working memory items inside and outside the focus of attention. Ann N Y Acad Sci. 2018;1424:256–67. 10.1111/nyas.13647 [Abstract] [CrossRef] [Google Scholar]
14. Rademaker RL, Serences JT. Pinging the brain to reveal hidden memories. Nat Neurosci. 2017;20:767–9. 10.1038/nn.4560 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
15. Orhan AE, Ma WJ. A diverse range of factors affect the nature of neural representations underlying short-term memory. Nat Neurosci. 2019;22:275–83. 10.1038/s41593-018-0314-y [Abstract] [CrossRef] [Google Scholar]
16. Bliss DP, D’Esposito M. Synaptic augmentation in a cortical circuit model reproduces serial dependence in visual working memory. PLoS ONE. 2017;12:e0188927. 10.1371/journal.pone.0188927 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
17. Carter E, Wang X-J. Cannabinoid-mediated disinhibition and working memory: dynamical interplay of multiple feedback mechanisms in a continuous attractor model of prefrontal cortex. Cereb Cortex. 2007;17(Suppl 1):i16–26. 10.1093/cercor/bhm103 [Abstract] [CrossRef] [Google Scholar]
18. Wolff MJ, Ding J, Myers NE, Stokes MG. Revealing hidden states in visual working memory using electroencephalography. Front Syst Neurosci. 2015;9:123. 10.3389/fnsys.2015.00123 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
19. LaRocque JJ, Lewis-Peacock JA, Drysdale AT, Oberauer K, Postle BR. Decoding attended information in short-term memory: an EEG study. J Cogn Neurosci. 2013;25:127–42. 10.1162/jocn_a_00305 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
20. Sprague TC, Ester EF, Serences JT. Restoring latent visual working memory representations in human cortex. Neuron. 2016;91:694–707. 10.1016/j.neuron.2016.07.006 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
21. Stokes MG, Kusunoki M, Sigala N, Nili H, Gaffan D, Duncan J. Dynamic coding for cognitive control in prefrontal cortex. Neuron. 2013;78:364–75. 10.1016/j.neuron.2013.01.039 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
22. Watanabe K, Funahashi S. Neural mechanisms of dual-task interference and cognitive capacity limitation in the prefrontal cortex. Nat Neurosci. 2014;17:601–11. 10.1038/nn.3667 [Abstract] [CrossRef] [Google Scholar]
23. Foster JJ, Sutterer DW, Serences JT, Vogel EK, Awh E. The topography of alpha-band activity tracks the content of spatial working memory. J Neurophysiol. 2016;115:168–77. 10.1152/jn.00860.2015 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
24. Barbosa J. Working memories are maintained in a stable code. J Neurosci. 2017;37:8309–11. 10.1523/JNEUROSCI.1547-17.2017 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
25. Bae G-Y, Luck SJ. Dissociable Decoding of Spatial Attention and Working Memory from EEG Oscillations and Sustained Potentials. J Neurosci. 2018;38:409–22. 10.1523/JNEUROSCI.2860-17.2017 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
26. Hakim N, Adam KCS, Gunseli E, Awh E, Vogel EK. Dissecting the Neural Focus of Attention Reveals Distinct Processes for Spatial Attention and Object-Based Storage in Visual Working Memory. Psychol Sci. 2019;30:526–40. 10.1177/0956797619830384 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
27. Wolff MJ, Akyurek EG, Stokes MG. What is the functional role of delay-related alpha oscillations during working memory? PsyArXiv. 2021. 10.31234/osf.io/z6y5b [CrossRef] [Google Scholar]
28. van Driel J, Olivers CNL, Fahrenfort JJ. High-pass filtering artifacts in multivariate classification of neural time series data. J Neurosci Methods. 2021;352:109080. 10.1016/j.jneumeth.2021.109080 [Abstract] [CrossRef] [Google Scholar]
29. Schneegans S, Bays PM. Restoration of fMRI Decodability Does Not Imply Latent Working Memory States. J Cogn Neurosci. 2017;29:1977–94. 10.1162/jocn_a_01180 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
30. Christophel TB, Iamshchinina P, Yan C, Allefeld C, Haynes J-D. Cortical specialization for attended versus unattended working memory. Nat Neurosci. 2018;21:494–6. 10.1038/s41593-018-0094-4 [Abstract] [CrossRef] [Google Scholar]
31. Iamshchinina P, Christophel TB, Gayet S, Rademaker RL. Essential considerations for exploring visual working memory storage in the human brain. Vis Cogn. 2021:1–12. 10.1080/13506285.2020.1848949 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
32. Wan Q, Cai Y, Samaha J, Postle BR. Tracking stimulus representation across a 2-back visual working memory task. R Soc Open Sci. 2020;7:190228. 10.1098/rsos.190228 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
33. van Loon AM, Olmos-Solis K, Fahrenfort JJ, Olivers CN. Current and future goals are represented in opposite patterns in object-selective cortex. Elife. 2018;7. 10.7554/eLife.38677 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
34. Yu Q, Teng C, Postle BR. Different states of priority recruit different neural representations in visual working memory. PLoS Biol. 2020;18:e3000769. 10.1371/journal.pbio.3000769 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
35. Spaak E, Watanabe K, Funahashi S, Stokes MG. Stable and dynamic coding for working memory in primate prefrontal cortex. J Neurosci 2017;37:6503–16. 10.1523/JNEUROSCI.3364-16.2017 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
36. Libby A, Buschman TJ. Rotational dynamics reduce interference between sensory and memory representations. Nat Neurosci. 2021;24:715–26. 10.1038/s41593-021-00821-9 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
37. Hajonides JE, Nobre AC, van Ede F, Stokes MG. Decoding visual colour from scalp electroencephalography measurements. Neuroimage. 2021;237:118030. 10.1016/j.neuroimage.2021.118030 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
38. Churchland MM, Yu BM, Cunningham JP, Sugrue LP, Cohen MR, Corrado GS, et al.. Stimulus onset quenches neural variability: a widespread cortical phenomenon. Nat Neurosci. 2010;13:369–78. 10.1038/nn.2501 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
39. Edin F, Klingberg T, Johansson P, McNab F, Tegnér J, Compte A. Mechanism for top-down control of working memory capacity. Proc Natl Acad Sci U S A. 2009;106:6802–7. 10.1073/pnas.0901894106 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
40. Seeholzer A, Deger M, Gerstner W. Stability of working memory in continuous attractor networks under the control of short-term plasticity. PLoS Comput Biol. 2019;15:e1006928. 10.1371/journal.pcbi.1006928 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
41. Kilpatrick ZP. Synaptic mechanisms of interference in working memory. Sci Rep. 2018;8:7879. 10.1038/s41598-018-25958-9 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
42. Fulvio JM, Postle BR. Cognitive control, not time, determines the status of items in working memory. J Cogn. 2020;3:8. 10.5334/joc.98 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
43. Zokaei N, Manohar S, Husain M, Feredoes E. Causal evidence for a privileged working memory state in early visual cortex. J Neurosci. 2014;34:158–62. 10.1523/JNEUROSCI.2899-13.2014 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
44. de Azevedo Neto RM, Bartels A. Disrupting short-term memory maintenance in premotor cortex serial dependence in visuomotor integration. J Neurosci. 2021; JN-RM-0380–21. 10.1523/JNEUROSCI.0380-21.2021 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
45. Wolff MJ, Jochim J, Akyürek EG, Buschman TJ, Stokes MG. Drifting codes within a stable coding scheme for working memory. PLoS Biol. 2020;18:e3000625. 10.1371/journal.pbio.3000625 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
46. Barbosa J. Neural network mechanisms of working memory interference [doctoral dissertation]. University of Barcelona; 2019. Available from: http://hdl.handle.net/2445/166717 [Google Scholar]
47. Stokes MG, Muhle-Karbe PS, Myers NE. Theoretical distinction between functional states in working memory and their corresponding neural states. Vis Cogn. 2020;28:420–32. 10.1080/13506285.2020.1825141 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
48. Moliadze V, Zhao Y, Eysel U, Funke K. Effect of transcranial magnetic stimulation on single-unit activity in the cat primary visual cortex. J Physiol (Lond). 2003;553:665–79. 10.1113/jphysiol.2003.050153 [Abstract] [CrossRef] [Google Scholar]
49. Wimmer K, Nykamp DQ, Constantinidis C, Compte A. Bump attractor dynamics in prefrontal cortex explains behavioral precision in spatial working memory. Nat Neurosci. 2014;17:431–9. 10.1038/nn.3645 [Abstract] [CrossRef] [Google Scholar]
50. Panichello MF, Buschman TJ. Shared mechanisms underlie the control of working memory and attention. Nature. 2021;592:601–5. 10.1038/s41586-021-03390-w [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
51. Zucker RS, Regehr WG. Short-term synaptic plasticity. Annu Rev Physiol. 2002;64:355–405. 10.1146/annurev.physiol.64.092501.114547 [Abstract] [CrossRef] [Google Scholar]
52. Zylberberg J, Strowbridge BW. Mechanisms of persistent activity in cortical circuits: possible neural substrates for working memory. Annu Rev Neurosci. 2017;40:603–27. 10.1146/annurev-neuro-070815-014006 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
53. Itskov V, Hansel D, Tsodyks M. Short-Term Facilitation may Stabilize Parametric Working Memory Trace. Front Comput Neurosci. 2011;5:40. 10.3389/fncom.2011.00040 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
54. Tegnér J, Compte A, Wang X-J. The dynamical stability of reverberatory neural circuits. Biol Cybern 2002;87:471–81. 10.1007/s00422-002-0363-9 [Abstract] [CrossRef] [Google Scholar]
55. Stein H, Barbosa J, Rosa-Justicia M, Prades L, Morató A, Galan-Gadea A, et al.. Reduced serial dependence suggests deficits in synaptic potentiation in anti-NMDAR encephalitis and schizophrenia. Nat Commun. 2020;11:4250. 10.1038/s41467-020-18033-3 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
56. Stokes MG. “Activity-silent” working memory in prefrontal cortex: a dynamic coding framework. Trends Cogn Sci (Regul Ed). 2015;19:394–405. 10.1016/j.tics.2015.05.004 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
57. Bae G-Y, Luck SJ. Reactivation of previous experiences in a working memory task. Psychol Sci. 2019;30:587–95. 10.1177/0956797619830398 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
58. Oostenveld R, Praamstra P. The five percent electrode system for high-resolution EEG and ERP measurements. Clin Neurophysiol. 2001;112:713–9. 10.1016/s1388-2457(00)00527-7 [Abstract] [CrossRef] [Google Scholar]
59. Tuckell HC. Introduction to Theoretical Neurobiology: Volume 2, Nonlinear and Stochastic Theories. Cambridge University Press; 1988. [Google Scholar]
60. Compte A, Brunel N, Goldman-Rakic PS, Wang XJ. Synaptic mechanisms and network dynamics underlying spatial working memory in a cortical network model. Cereb Cortex. 2000;10:910–23. 10.1093/cercor/10.9.910 [Abstract] [CrossRef] [Google Scholar]
61. Goodman D, Brette R. Brian: a simulator for spiking neural networks in python. Front Neuroinform. 2008;2:5. 10.3389/neuro.11.005.2008 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
2021 Oct; 19(10): e3001436.
Published online 2021 Oct 21. 10.1371/journal.pbio.3001436.r001

Decision Letter 0

Gabriel Gasque, Senior Editor

2 Apr 2021

Dear Albert,

Thank you for submitting your manuscript entitled "Unattended short-term memories are maintained in active neural representations" for consideration as a Short Report by PLOS Biology.

Your manuscript has now been evaluated by the PLOS Biology editorial staff, as well as by an academic editor with relevant expertise, and I am writing to let you know that we would like to send your submission out for external peer review.

However, before we can send your manuscript to reviewers, we need you to complete your submission by providing the metadata that is required for full assessment. To this end, please login to Editorial Manager where you will find the paper in the 'Submissions Needing Revisions' folder on your homepage. Please click 'Revise Submission' from the Action Links and complete all additional questions in the submission questionnaire.

Please re-submit your manuscript within two working days, i.e. by Apr 06 2021 11:59PM.

Login to Editorial Manager here: https://www.editorialmanager.com/pbiology

During resubmission, you will be invited to opt-in to posting your pre-review manuscript as a bioRxiv preprint. Visit http://journals.plos.org/plosbiology/s/preprints for full details. If you consent to posting your current manuscript as a preprint, please upload a single Preprint PDF when you re-submit.

Once your full submission is complete, your paper will undergo a series of checks in preparation for peer review. Once your manuscript has passed all checks it will be sent out for review.

Given the disruptions resulting from the ongoing COVID-19 pandemic, please expect delays in the editorial process. We apologise in advance for any inconvenience caused and will do our best to minimize impact as far as possible.

Feel free to email us at gro.solp@ygoloibsolp if you have any queries relating to your submission.

Kind regards,

Gabriel Gasque, Ph.D.,

Senior Editor

PLOS Biology

    2021 Oct; 19(10): e3001436.
    Published online 2021 Oct 21. 10.1371/journal.pbio.3001436.r002

    Decision Letter 1

    Gabriel Gasque, Senior Editor

    6 May 2021

    Dear Albert,

    Thank you very much for submitting your manuscript "Unattended short-term memories are maintained in active neural representations" for consideration as a Short Report at PLOS Biology. Your manuscript has been evaluated by the PLOS Biology editors, by an Academic Editor with relevant expertise, and by three independent reviewers. You will note that reviewers 1 and 3 have revealed their identities.

    In light of the reviews (below), we will not be able to accept the current version of the manuscript, but we would welcome re-submission of a revised version that takes into account the reviewers' comments. We cannot make any decision about publication until we have seen the revised manuscript and your response to the reviewers' comments. Your revised manuscript is also likely to be sent for further evaluation by the reviewers.

    We expect to receive your revised manuscript within 3 months.

    Please email us (gro.solp@ygoloibsolp) if you have any questions or concerns, or would like to request an extension. At this stage, your manuscript remains formally under active consideration at our journal; please notify us by email if you do not intend to submit a revision so that we may end consideration of the manuscript at PLOS Biology.

    **IMPORTANT - SUBMITTING YOUR REVISION**

    Your revisions should address the specific points made by each reviewer. As you will see, all reviewers agree the work is timely and important, and most of their concerns can be addressed with textual changes --we think-- regarding issues of interpretation and providing appropriate context and discussion of other related work. In addition, reviewers 2 and 3 also suggest some additional analyses that could bolster your conclusions. Having discussed these comments with the academic editor, we think it is important that you include those new analyses in the current manuscript.

    Please submit the following files along with your revised manuscript:

    1. A 'Response to Reviewers' file - this should detail your responses to the editorial requests, present a point-by-point response to all of the reviewers' comments, and indicate the changes made to the manuscript.

    *NOTE: In your point by point response to the reviewers, please provide the full context of each review. Do not selectively quote paragraphs or sentences to reply to. The entire set of reviewer comments should be present in full and each specific point should be responded to individually, point by point.

    You should also cite any additional relevant literature that has been published since the original submission and mention any additional citations in your response.

    2. In addition to a clean copy of the manuscript, please also upload a 'track-changes' version of your manuscript that specifies the edits made. This should be uploaded as a "Related" file type.

    *Re-submission Checklist*

    When you are ready to resubmit your revised manuscript, please refer to this re-submission checklist: https://plos.io/Biology_Checklist

    To submit a revised version of your manuscript, please go to https://www.editorialmanager.com/pbiology/ and log in as an Author. Click the link labelled 'Submissions Needing Revision' where you will find your submission record.

    Please make sure to read the following important policies and guidelines while preparing your revision:

    *Published Peer Review*

    Please note while forming your response, if your article is accepted, you may have the opportunity to make the peer review history publicly available. The record will include editor decision letters (with reviews) and your responses to reviewer comments. If eligible, we will contact you to opt in or out. Please see here for more details:

    https://blogs.plos.org/plos/2019/05/plos-journals-now-open-for-published-peer-review/

    *PLOS Data Policy*

    Please note that as a condition of publication PLOS' data policy (http://journals.plos.org/plosbiology/s/data-availability) requires that you make available all data used to draw the conclusions arrived at in your manuscript. If you have not already done so, you must include any data used in your manuscript either in appropriate repositories, within the body of the manuscript, or as supporting information (N.B. this includes any numerical values that were used to generate graphs, histograms etc.). For an example see here: http://www.plosbiology.org/article/info%3Adoi%2F10.1371%2Fjournal.pbio.1001908#s5

    *Blot and Gel Data Policy*

    We require the original, uncropped and minimally adjusted images supporting all blot and gel results reported in an article's figures or Supporting Information files. We will require these files before a manuscript can be accepted so please prepare them now, if you have not already uploaded them. Please carefully read our guidelines for how to prepare and upload this data: https://journals.plos.org/plosbiology/s/figures#loc-blot-and-gel-reporting-requirements

    *Protocols deposition*

    To enhance the reproducibility of your results, we recommend that if applicable you deposit your laboratory protocols in protocols.io, where a protocol can be assigned its own identifier (DOI) such that it can be cited independently in the future. Additionally, PLOS ONE offers an option for publishing peer-reviewed Lab Protocol articles, which describe protocols hosted on protocols.io. Read more information on sharing protocols at https://plos.org/protocols?utm_medium=editorial-email&utm_source=authorletters&utm_campaign=protocols

    Thank you again for your submission to our journal. We hope that our editorial process has been constructive thus far, and we welcome your feedback at any time. Please don't hesitate to contact us if you have any questions or comments.

    Sincerely,

    Gabriel Gasque

    Senior Editor

    PLOS Biology

    gro.solp@euqsagg

    *****************************************************

    REVIEWS:

    Reviewer #1, Brad Postle: "Unattended short-term memories are maintained in active neural representations," Barbosa, Lozano-Soldevilla, and Compte. This manuscript present the reanalysis of data from two publications (from other groups) about which, as the authors write, "The[ir] interpretation of reactivations from absent decoding as evidence for 'activity-silent' storage has had a strong impact in the field." The major findings are threefold. First, reanalysis of the data from Wolfe et al. (2017) replicated the original finding of failing to decode evidence for an active representation of the unprioritized memory item (UMI) when the EEG data were analyzed in their voltage format, but SUCCEEDED in decoding evidence for an active representation of the UMI from the alpha-band component of the signal, after these same data had been spectrally transformed. Second, they show with simulations, then empirically, that the effect of the visual "ping" is to decrease trial-to-trial variability in the signal. Third, they show with simulations, then with empirical analysis of the data from Rose et al., (2016), that the effect of a pulse of TMS is the opposite of that of a ping: TMS increases trial-to-trial variability. The significance of (2) and (3) is that (2) may explain why the effect of the ping on time-voltage data is to briefly rescue decoding of the UMI, whereas (3) would seem to indicate that the activity-silent interpretation of TMS-evoked rescue of decoding of the UMI remains viable (although, as this review will highlight, this interpretation is equivocal in the manuscript as currently written). This manuscript will be of considerable interest to many who study working memory from a variety of perspectives (computation, behavior, intra- and extracranial electrophysiology, fMRI, neurostimulation…), because of its important implications at both empirical and theoretical levels. (As the senior author of Rose et al. (2016) I can't not forgo anonymity and still write a comprehensive review.)

    This work is important and will get a lot of attention, but I worry that some of the framing fails to emphasize some of the most important elements and pushes a narrative that's both too strong and somewhat out of date. Let's start with the title - it's stating proposition that has already been described on several occasions, including in one paper that they cite, but not in this context (Christophel et al. 2018), and several that they don't (van Loon et al. 2018, Wan et al. 2020, Yu et al. 2020, Libby and Buschman 2021). This latter group all document a phenomenon that the decoding methods used here can't distinguish, but that the authors might consider investigating with a different method, which is that the active representation of the UMI is transformed from its format when it is a prioritized (P)MI. (Two of these papers are from my group, and I want to be clear that I have no expectations that the authors need to cite them.) Furthermore there's an implication here and in many parts of the manuscript that the UMI in the Wolfe et al. studies is NOT also stored in an activity-silent format, but of course the authors can't know that. (Finally, with regard to the title, the authors of both (Wolff et al. 2017) and (Rose et al. 2016) refer to their tasks as "working memory" tasks, and the authors might consider using this label, which is much more widely used in the field.) This false dichotomy (active vs. activity-silent) is particularly dissonant in the final paragraph, in which the authors seem to ignore their own (very elegant) work that showed compelling positive evidence for activity-silent representations in the PFC of the monkey. Also important to consider are different results from the monkey PFC in which a "cognitive ping" reveals otherwise undetectable (and so 'probably' activity-silent) representations of stimuli (Stokes et al. 2013). Finally, as the authors (seemingly grudgingly) acknowledge, the most straightforward conclusion from their reanalysis of the data from Rose et al. (2016) is the UMI in that study most likely was held via activity-silent mechanism. Finally, it need not be the case that the idea of activity-silent representation "stands in contradiction with computational models of 'activity-silent' storage, where short-term plasticity changes are induced by neuronal activity" if one allows for existence of cognitive control. Indeed, Jacqueline Fulvio in my group has show behaviorally that the 'behavioral reactivation' TMS effect is subject to control (Fulvio and Postle 2020). (Again, I'm not asking for citations, but it is the case that I work in this area …) Masse et al. (Masse et al. 2019) and Manohar et al. (Manohar et al. 2019) have demonstrated that the two formats can, in principle, co-exist, and if the authors embraced this idea they could preempt readers like me getting distracted by an issue that's not the main point of the results.

    I realize that I spilled a lot of ink to make this point, but it is really the only big-ish concern that I have with this otherwise excellent paper. I'll make more specific comments in the order in which they appear.

    2nd paragraph: Rose el al. (2016) also reported an inability to decode the UMI from spectrally transformed data (Figure S5, which also show that the PMI is decodable from alpha and the UMI reactivation from beta).

    Results and Discussion: The first time the key result is mentioned "We found that a sustained alpha power code tracks the items that remain relevant for future behavior" it'd be helpful to specify that it's tracking orientation; during my first read-through I was uncertain whether it was orientation or location-on-the-screen that was being decoded.

    "Furthermore, it challenges the current view on the role of alpha power during working memory maintenance [18], which would suppress immediately irrelevant memories [19]." This is unclear, because 18 argues against a suppression/inhibition for alpha?

    This clause needs more unpacking "possibly reflecting the prioritization of strong competition between actively held memories in attractor networks [20,21]" because someone unfamiliar with the details of 20, 21 won't necessarily understand what is meant here.

    Is the work of Bae and Luck (recovery of decoding of previous trial's stimulus during the next trial) also relevant here, perhaps as part of a more detailed consideration of the implications of (Barbosa et al. 2020)?

    Data preprocessing: it's unclear what is meant by the words "to revert this baselining." Additionally, a sentence or two explaining in more detail how this was done, and why it was important to do, would be helpful.

    "funcion" is not English.

    Figure 1: this is picky, but "decoding from alpha power (Methods), which reveals a strong sustained code of the unattended stimulus" strikes me as imprecise. What's being revealed is a stimulus representation, and from that one infers that there is a code that supports this representation, right?

    Figure 2: "A signal-to-noise ratio increase can reflect a drop in variance …" Isn't it meant that an SNR increase can result from a drop in variance?

    Signed, Brad Postle

    Barbosa, J., H. Stein, R. L. Martinez, A. Galan-Gadea, S. Li, J. Dalmau, K. C. S. Adam, J. Valls-Solé, C. Constantinidis and A. Compte (2020). "Interplay between persistent activity and activity-silent dynamics in the prefrontal cortex underlies serial biases in working memory." Nature Neuroscience 23: 1016-1024.

    Christophel, T. B., P. Iamshchinina, C. Yan, C. Allefeld and J.-D. Haynes (2018). "Cortical specialization for attended versus unattended working memory." Nature Neuroscience 21: 494-496.

    Fulvio, J. M. and B. R. Postle (2020). "Cognitive control, not time, determines the status of items in working memory." Journal of Cognition 3: 1-8.

    Libby, A. and T. J. Buschman (2021). "Rotational dynamics reduce interference between sensory and memory representations." Nature Neuroscience.

    Manohar, S. G., N. Zokaei, S. J. Fallon, T. P. Vogels and M. Husain (2019). "Neural mechanisms of attending to items in working memory." Neuroscience and Biobehavioral Reviews 101: 1-12.

    Masse, N. Y., G. R. Yang, H. F. Song, X.-J. Wang and D. J. Freedman (2019). "Circuit mechanisms for the maintenance and manipulation of information in working memory." Nature Neuroscience 22: 1159-1167.

    Rose, N., J. J. Larocque, A. C. Riggall, O. Gosseries, M. J. Starrett, E. Meyering and B. R. Postle (2016). "Reactivation of latent working memories with transcranial magnetic stimulation." Science 354: 1136-1139.

    Stokes, M. G., M. Kusunoki, N. Sigala, H. Nili, D. Gaffan and J. Duncan (2013). "Dynamic coding for cognitive control in prefrontal cortex." Neuron 78(2): 364-375.

    van Loon, A. M., K. Olmos-Solis, J. J. Fahrenfort and C. N. L. Olivers (2018). "Current and future goals are represented in opposite patterns in object-selective cortex." eLife 7: e38677.

    Wan, Q., Y. Cai, J. Samaha and B. R. Postle (2020). "Tracking stimulus representation across a 2-back visual working memory task." Royal Society Open Science 7: 190228

    Wolff, M. J., J. Jochim, E. G. Akyürek and M. G. Stokes (2017). "Dynamic hidden states underlying working-memory-guided behavior." Nature Neuroscience.

    Yu, Q., C. Teng and B. R. Postle (2020). "Different states of priority recruit different neural codes in visual working memory." PLoS Biology 18: e3000769.

    Reviewer #2: Summary: The authors investigate whether they can decode information that was previously thought to be maintained in silent working memory. The authors ask this question by using alpha power at posterior electrodes to decode working memory representations. In at least one experiment, they are able to decode working memory information. Furthermore, the authors then suggest that visual "pinging" may change the signal-to-noise ratio of EEG activity, thus increasing decoding accuracy of already active neural representations (by increasing SNR), rather than re-activating latent traces. Many have pointed out already that absence of evidence is weak evidence for absence of decodable memory activity, and this study is an important example of this point. Given the strong influence of the "activity silent" models of working memory storage, the present findings provide a critical alternative explanation for one of the more prominent studies arguing in favor of activity-silent modes of storage in working memory. The present authors also report simulations that support their hypothesis regarding the SNR effects of visual "pinging" (i.e., presentation of an irrelevant visual stimuli), showing that pinging may reduce across trial variability and thus increase SNR. This adds strength to the authors' speculations that the putative "silent" representations may have simply been masked by noise rather than truly silent.

    However, the authors' account doesn't provide a direct explanation for why it was the *relevant* and not the irrelevant memory representation that was decodable in the raw voltages after visual pinging. That is, while it is clear that they have successfully decoded using alpha power, why would visual pinging only resurrect the relevant item's representation if the authors are correct about the effects of pinging on across trial variability? Those variability effects should influence decoding of both relevant and irrelevant items, shouldn't they?

    The authors assert that "…alpha power tightly tracks working memory contents, *regardless* of their immediate behavioral relevance." (p. 3, results and discussion), but I thought figure 1 was ambiguous on this point. 1C shows a trend towards higher decoding strength for the relevant item, especially at the end of the delay period. Figure 1D also seems to show a divergence of the relevant and irrelevant decoding strengths at 1sec, but I wasn't sure why the irrelevant item line was dashed instead of solid (perhaps this is indicated somewhere in the manuscript, but I couldn't find it. I'd recommend making this more clear in the figure caption).

    This then leaves the mystery of why only the relevant item appears to be tracked by EEG voltage following the ping. The authors do state, "However, while attended memories are decodable both in alpha power and voltage traces, unattended memories are only detected in alpha power, *possibly reflecting the prioritization of strong competition between actively held memories in attractor networks.*" Does this sentence imply that the traces reflected in EEG voltage were indeed "silent" prior to the pinging? I thought this was unclear from the manuscript. The later arguments about pinging reducing across trial variability seem to leave open the possibility that both attended and unattended items were represented in EEG voltage, but with SNR too low to confirm it with the original analysis pipeline. I think this points merits careful clarification in the manuscript.

    The authors also point out that TMS reactivation may be qualitatively different, because simulation suggest that TMS *increases* across trial variability, and that this falls in line with the predictions of "a biophysical network model of memory reactivation from silent, synaptic traces". I thought this was an interesting point, but it seems to sparsely described to really have a strong impact in the paper. Moreover, while being "potentially consistent" with that biophysical model is interesting, how strong is this evidence? I would be surprised if models that denied the role of "activity-silent" representations in working memory could not also be "consistent" with the finding that across trial variability is increased following TMS. But if the authors have a compelling argument that I should be surprised here, I think it should be clearly spelled out in the paper. If there is no strong argument, then I'd recommend tempering this particular conclusion. My view is that the really clear evidence in this study is the positive findings with alpha power in the pinging studies, and the simulations of how changes in SNR could explain the original findings in the pinging study. The other speculations regarding the potential causes of the TMS findings are not yet as convincing.

    Minor points:

    1. The authors use the same analysis pipeline as the originally published papers, with the inclusion of alpha power instead of raw EEG amplitude. This is a powerful approach. However, the original analysis pipeline only included 17 posterior electrodes. It is possible that information about working memory could also be represented in other electrodes on the scalp. Therefore, decoding accuracy could actually be higher if the authors deviated from the original methods and included all electrodes in their models. It may be useful to include an additional analysis that addresses whether a model that includes all electrodes increases decoding accuracy of working memory representations.

    2. In the current manuscript, the authors clearly show that alpha power tracks working memory representations in the time periods where the prior work suggested that there were no neurally active representations. This is the critical empirical pattern. That said, recent work suggests that alpha power and EEG voltage may play qualitatively different roles in working memory tasks:

    Bae, G. Y., & Luck, S. J. (2018). Dissociable decoding of spatial attention and working memory from EEG oscillations and sustained potentials. Journal of Neuroscience, 38(2), 409-422.

    Hakim, N., Adam, K. C., Gunseli, E., Awh, E., & Vogel, E. K. (2019). Dissecting the neural focus of attention reveals distinct processes for spatial attention and object-based storage in visual working memory. Psychological Science, 30(4), 526-540.

    So, it may be worthy of some discussion whether the active neural signal the present authors have identified might reflect a different aspect of working memory maintenance than do the patterns of activity in EEG voltage. This could be an interesting compromise between the original framing of the prior reactivation studies and the present framing.

    3. The published data was baselined, and the authors reversed this baselining for their analyses. The authors should provide further explanation for why they reversed this baseline and should additionally discuss whether they were able to decode working memory representations with the published, baselined data.

    4. The authors introduce fano-factors in Figure 3, but do not include a discussion of fano-factor anywhere else in the paper. Given the prominence of "[big up triangle, open]fano-factor" in their figure, the authors should include a brief description of fano-factors in the Method's section of their manuscript. This description will help bridge the literature that typically uses fano-factors to the literature that typically investigates alpha power activity in human EEG.

    Conclusions:

    The authors make two main conclusions (1) working memory representations that were previously thought to be maintained in an activity silent state have been shown to be tracked via an active neural trace in the alpha band, and (2) visual "pinging" changes the ratio of signal-to-noise in EEG signals, thereby allowing active signals that may have been masked by noise to be detected. This work provides an important new interpretation of a highly influential set of studies, and I believe it would have a strong positive impact in the literature that will generate vigorous follow-up work. Their conclusions, however, would be strengthened by addressing the above comments.

    Reviewer #3, Thomas Christophel and Vivien Chopurian: The authors report a short reanalysis of EEG data from recent studies investigating reactivation (via 'pinging') of supposedly 'activity-silent' mnemonic traces during working memory. This reanalysis shows that alpha band activity carries an active trace of memorized content, in contrast to this prior work. Additional reanalyses concern the nature of the 'pinging' effects reported in this and other prior work. The manuscript argues for differential effects of 'pinging' using either sensory mask-like stimuli and TMS on the reduction and increase of across-trial EEG variability, respectively.

    The finding that data in a core study seemingly supporting 'silent' working memory contains evidence for active neural representations alone is critical to the field and beyond. This finding is robust and consistent with previous work and without any doubt deserves the attention given here. There are opportunities to elaborate more on this finding and its relationship to other work, but it is convincing as is. The additional reanalyses and their interpretations highlight some explanations for 'pinging' effects, but several concerns let me doubt the conclusions drawn here.

    We will outline our suggestions, concerns, and recommendations in the following sections. Beforehand, we would like to emphasize that our expertise primarily lies in related fMRI work.

    Major:

    * Our concerns solely relate to the second part of the manuscript which argues that "visual pinging reveals an underlying active code by quenching EEG noise". As the authors outline this is indeed a plausible explanation, but it is questionable whether the data analyses reported are sufficient to support the hypothesis in a substantive way.

    As the authors outline, a reduction in trial-by-trial variability is the expected outcome of introducing a stimulus into a system that is constant across trials. Simply put, stimulation (like a 'ping') can be seen to replace endogenous by with exogenous activity. If this exogenous activity is constant across trials, variability across trials is reduced. It is equally plausible that in a decoding analysis, the response related to a memorized item is more easily decoded when noise is reduced. The 'noise' in this decoding analysis and the 'trial-by-trial variability' in the analysis provided, however, are not the same. On the contrary, the trial-by-trial variability is composed of both the noise and the signal used for the decoding analysis. This becomes obvious as one considers the fact that subjects memorize different orientations on different trials, making the signal itself vary across trials. How different components (ping signal, mnemonic signal, and residual noise) are aggregated to form the overall EEG signal is unknown, but here we see little indication that trial-by-trial variability can be used as an exclusive indicator of the residual noise component. For these reasons, we do not see how the trial-by-trial variability analyses provide evidence in favor or against the 'quenching noise' hypothesis.

    * This problem becomes even clearer when considering the trial-by-trial variability in Rose et al: In Rose, not only does the content vary across trials, but the ping is variable across trials, too (in the EEG Exp 2). Considering that this means that on different trials the TMS coil is pointed at different parts of the brain, one might expect to see it in the absence of any mnemonic or cognitive activity (e.g. in a lifeless neuronal substrate). Notably, this also provides a relatively straight forward explanation for the 'reactivation' in Rose et al. (in their Exp 2, at least): The TMS itself already carries information about the memorized content, which makes decoding the content after the pulse a trivial finding.

    * What It is essential here,, is that a critical finding (active representations in the Wolff data) is not hampered in impact, by jumping to conclusions in the second part of the manuscript. While the 'quenching noise' interpretation is a plausible one, alternative interpretations need to be considered. One simple interpretation is that 'pings' increase overall alertness or modality-specific attention which in turn increases the signal of any memorized content. This could be seen as a mechanism to counter interference by distracting stimulation (see Bettencourt and Xu as well as Rademaker et al., in particular, the difference between expected and unexpected distraction in Bettencourt and Xu). This highlights that a robust explanation of 'pinging' effects must consider both enhancing ('pinging') and disrupting ('distracting') effects of stimulus presentations during working memory delays which are likely to jointly affect neural decoding in different ways across studies.

    Minor:

    * The main finding of the paper is a rather simple but important one: If you run the same analysis as Wolff et al. (2017) with the same data but using Alpha power rather than raw data, then you can decode unattended items, which were not decodable before. One potential avenue to make this manuscript even more insightful would be to provide further analyses that help explaining this finding.

    As the authors mention, raw EEG power already carries information about a memorized item in Wolff et al. (2015) when combining 'pinging' and 'no-pinging' trials prior to the 'ping'. The same analyses split half for 'pinging' trials results in an occasional null result ('long trials only'). This highlights the trivial but apparently forgotten idea that statistical power might be the essential determinant of whether one is likely to find a positive result in a given study. The change from raw signal to Alpha power might simply constitute an increase in per trial effect size (like other forms of filtering and smoothing). In other words, the question is whether the difference between raw and Alpha -based analyses represents a qualitative difference in the underlying signals (between UMI and AMI) or simply a quantitative difference in the strength of representation.

    We suggest quantifying this effect by running raw-EEG and Alpha decoding for both studies and 'simulating' studies with different sample sizes (by randomly removing subjects from a given study) or even different number of trials (e.g. split half). We also suggest running analyses combining decoding accuracy across time-points (e.g. averaging across the delay) on a given trial and/or combining trials prior to the decoding analyses (similar to run-wise beta decoding in fMRI) to increase statistical power for raw analyses. We want to stress that these are suggestions to better explain an already relevant findings rather than analyses necessary to support the finding. Not all these suggestions might be feasible for the datasets available.

    * One final recommendation is to investigate the differential involvement of different electrode positions in the decoding of AMIs and UMIs. We want to be very clear, however, that is more driven by the reviewers' curiosity than anything else.

    * "[…] moving electrodes using the python function signal.detrend on each subject variability." We assume this refers to the function in the scipy package and the linear detrending option thereof, please specify. There also might be a typo here ("each subject's").

    * "Mahalanobis distance between all possible pairwise combinations of the orientations and thus form a tuning curve." Please specify how these pairwise distances are turned into a tuning curve. An array of pairwise distances is typically referred to as a representational dissimilarity matrix.

    * "We realized that an earlier study had reported that attended spatial memories were decoded more reliably from EEG alpha power than from raw voltages [16]." This might be related to our primary expertise in fMRI data, but when reading Foster et al. (2016), our understanding is that they differentiate between evoked and total power not raw voltages and Alpha power. See David et al. (2006): "In short, evoked responses can be characterized as the power of the average; while induced responses are the average power that cannot be explained by the power of the average." 'Total' power would then be both components combined which (in our understanding) is what is used, here.

    Naturally, there is a relationship between evoked power and evoked responses in preprocessed (rather than raw) EEG data, so the results by Foster et al. still can serve as a motivation for the reanalysis in the current manuscript. The authors should however clarify what prior work found.

    * "Despite their relevance for upcoming memory-guided behavior, currently unattended memories cannot be robustly decoded from raw EEG voltage traces [6,7] (Figure 1a, red)." This might be somewhat misleading as prior work only showed that they were not able to decode from raw voltages which might simply be a false negative. The prior work does not show that UMIs "cannot" - in principle - be decoded.

    * "Furthermore, it challenges the current view on the role of alpha power during working memory maintenance [18], which would suppress immediately irrelevant memories [19]." Here, it might be worth noting that suppression might lead to inverted tuning which could be decodable equally well than non-inverted tuning. This means that being able to decode an item does not necessarily mean that it isn't suppressed or that the signal one decodes isn't suppressive in nature. (See the different works by Lorenc and Postle and van Loon). Notably, there is plenty of debate about these inverted tuned representations, so it is unclear to me whether to include the debate here.

    * "Despite apparent similarities, TMS reactivations impact working memory performance [5,6], while pinging does not [7]." Here it is worth noting that (by our count) in Rose et al. only one out of three TMS experiments show a behavioral effect of the TMS (one out of six possible effects in exp. 4 at p = 0.01, which is not reanalyzed in the current manuscript). Overall effects of interfering stimulation with distractor or TMS are rather scarce, so we caution to interpret this apparent difference more carefully.

    * "Decoding from raw voltage or alpha power was done with the exact same code, but preprocessing the data differently.": The second sentence seems to miss a "with" or we suggest changing the voice from passive to active

    * "Each 'trial' trace was simulated as a slope α, different for each stimulus (α1 = 1 and α1 = 1/2) on top of noise sampled from a normal distribution ξ, with mean 0 and standard deviation 1." 'α1 = ½' should probably be 'α2 = ½'. It would be helpful to explain the selection of these model parameters (why 1 and 1/2)? Further it might be helpful to spell out the simulation a little (e.g. What are the two items, what is the data being simulated, where does the data start, where is 0 ….), it took us more time than necessary to get what the authors are doing here and why. Arguably, however, the whole simulation can be seen as demonstrating the trivial fact that SNR means Signal to Noise ratio, but we'll leave it to the authors judgement whether explaining this is necessary.

    We wish you all the best for this fascinating project.

    Thomas Christophel and Vivien Chopurian

      2021 Oct; 19(10): e3001436.
      Published online 2021 Oct 21. 10.1371/journal.pbio.3001436.r003

      Author response to Decision Letter 1

      3 Aug 2021

      Attachment

      Submitted filename:

        2021 Oct; 19(10): e3001436.
        Published online 2021 Oct 21. 10.1371/journal.pbio.3001436.r004

        Decision Letter 2

        Gabriel Gasque, Senior Editor

        30 Aug 2021

        Dear Albert,

        Thank you for submitting your revised Short Report entitled "Pinging reveals active, not silent, working memories" for publication in PLOS Biology. I have now obtained advice from the original reviewers and have discussed their comments with the Academic Editor. 

        Based on the reviews, we will probably accept this manuscript for publication, provided you satisfactorily address the remaining points raised by reviewer 3. Please also make sure to address the data and other policy-related requests listed below my signature.

        As you address these items, please take this last chance to review your reference list to ensure that it is complete and correct. If you have cited papers that have been retracted, please include the rationale for doing so in the manuscript text, or remove these references and replace them with relevant current references. Any changes to the reference list should be mentioned in the cover letter that accompanies your revised manuscript.

        We expect to receive your revised manuscript within two weeks.

        To submit your revision, please go to https://www.editorialmanager.com/pbiology/ and log in as an Author. Click the link labelled 'Submissions Needing Revision' to find your submission record. Your revised submission must include the following:

        -  a cover letter that should detail your responses to any editorial requests, if applicable, and whether changes have been made to the reference list

        -  a Response to Reviewers file that provides a detailed response to the reviewers' comments (if applicable)

        -  a track-changes file indicating any changes that you have made to the manuscript. 

        NOTE: If Supporting Information files are included with your article, note that these are not copyedited and will be published as they are submitted. Please ensure that these files are legible and of high quality (at least 300 dpi) in an easily accessible file format. For this reason, please be aware that any references listed in an SI file will not be indexed. For more information, see our Supporting Information guidelines:

        https://journals.plos.org/plosbiology/s/supporting-information  

        *Published Peer Review History*

        Please note that you may have the opportunity to make the peer review history publicly available. The record will include editor decision letters (with reviews) and your responses to reviewer comments. If eligible, we will contact you to opt in or out. Please see here for more details:

        https://blogs.plos.org/plos/2019/05/plos-journals-now-open-for-published-peer-review/

        *Early Version*

        Please note that an uncorrected proof of your manuscript will be published online ahead of the final version, unless you opted out when submitting your manuscript. If, for any reason, you do not want an earlier version of your manuscript published online, uncheck the box. Should you, your institution's press office or the journal office choose to press release your paper, you will automatically be opted out of early publication. We ask that you notify us as soon as possible if you or your institution is planning to press release the article.

        *Protocols deposition*

        To enhance the reproducibility of your results, we recommend that if applicable you deposit your laboratory protocols in protocols.io, where a protocol can be assigned its own identifier (DOI) such that it can be cited independently in the future. Additionally, PLOS ONE offers an option for publishing peer-reviewed Lab Protocol articles, which describe protocols hosted on protocols.io. Read more information on sharing protocols at https://plos.org/protocols?utm_medium=editorial-email&utm_source=authorletters&utm_campaign=protocols

        Please do not hesitate to contact me should you have any questions.

        Sincerely,

        Gabriel Gasque, Ph.D.,

        Senior Editor,

        gro.solp@euqsagg,

        PLOS Biology

        ------------------------------------------------------------------------

        TITLE:

        We would like to decompress your title a little, to make it more accesible and appealing to our broad readership. We recommend:

        Pinging the brain with visual impulses reveals an electrically active code for unattended working memories.

        However, we would be happy to discuss alternatives, if you think our suggestion is inaccurate or misrepresents your findings.

        ------------------------------------------------------------------------

        BLURB:

        Please provide a blurb, which will be included in our weekly and monthly Electronic Table of Contents, sent out to readers of PLOS Biology, and may be used to promote your article in social media. The blurb should be about 30-40 words long and is subject to editorial changes. It should, without exaggeration, entice people to read your manuscript. It should not be redundant with the title and should not contain acronyms or abbreviations.

        ------------------------------------------------------------------------

        DATA POLICY:

        Please include in your Data Statement in the submission system an explicit statement on where the original data can be found (“the original publications.”), while keeping the following “All the code used for the analyses can be found at https://github.com/comptelab/reactivations”

        Please also ensure that each figure legend in your manuscript includes information on where the underlying data can be found.

        ------------------------------------------------------------------------

        DATA NOT SHOWN:

        Please note that per journal policy, we do not allow the mention of "data not shown", "personal communication", "manuscript in preparation" or other references to data that is not publicly available or contained within this manuscript. Please either remove mention of these data or provide figures presenting the results and the data underlying the figure(s).

        ------------------------------------------------------------------------

        Reviewer remarks:

        Reviewer #1, Brad Postle: That authors have thoroughly and satisfactorily addressed all the points raised in my initial review.

        Reviewer #2: I think the authors have done a strong job with the revision, and I believe the paper is ready for publication.

        Reviewer #3, Thomas Christophel & Vivien Chopurian: The authors provide an extensive response to our comments. The revised manuscript includes several additional analyses and simulation which critically add to the prior work. The updated main finding shows that unattended memory items can be decoded from both raw and frequency transformed EEG data. Moreover, 'reactivations' by 'pinging' appear to depend on the presence of an active raw EEG trace of the unattended item. Finally, as an alternative explanation for pinging effects, they introduce an attractor-boost model which appear to equally explain their results. The authors have broadened their interpretation of their initial findings in the light of these additional findings and discuss possible underlying mechanisms.

        These changes greatly improve the quality of an already noteworthy study. Some open points remain. These concerns predominantly relate to the interpretation of Rose et al., which play a lesser role in the revised manuscript. We will outline our suggestions, concerns, and recommendations in the following sections.

        Minor:

        1. Data from Rose et al. (2016) is reanalyzed for the current study and plays a part in the conclusions drawn by the authors (e.g. in the discussion: 'We further compare visual pinging with TMS perturbations, and we find qualitative differences suggesting different underlying mechanisms'). The authors however did not analyze the TMS-EEG data in the same way as the Pinging data and conclude different mechanisms based on one analysis looking at difference in variation measures between the two sets of experiments.

        The authors now acknowledge that variation in the TMS site might result in an increase in variability but appear to contend that reactivations (meaning an interaction of the TMS with the respective conditions) are the more probable cause of this increase, as the TMS itself should average out due to the block-wise nature of their analysis. When looking at the data directly, however, we see that the by far largest increase in variability (so large that it evades the y-axis on the respective plot) occurs at 0ms relative to the pulse. In other words, the largest change in the EEG measurements is the pulse itself which - according to the authors - should have been averaged out. Hence, we would attribute the differences in variability found between pinging and TMS to variations in the TMS itself not a difference in the underlying mechanism of the reactivation which might be sensibly explained by boosted active representations or even by noise-quenching (which would be overshadowed TMS-related noise). We suggest a more careful interpretation of these differences in variability and possible underlying mechanisms.

        2. For the across-trial variability analyses, please clarify the block-wise nature of this analysis in the methods section.

        3. Please extend the y-axis for the variability analyses such that the reader can see the extend of the peak for the data from Rose et al. (2016).

        4. "We then tested the influence of sample size on decoding estimates in each condition, by varying the number of sessions and trials included in the analyses. This showed a striking difference between unattended and discarded memories: while increasing the sample size in the unattended condition resulted in a monotonic increase of t-value, it did not for discarded memories (see diagonal in Fig. 2). This result suggests that increasing the number of sessions would lead to decodability of unattended memories, but not of discarded memories." This monotonic increase might be less striking than one might think. In a bootstrapping procedure like this, the closer a randomly drawn sample of trials gets to the full set the closer it will approximate the results with all trials included. For unattended items this results appears to be p < 0.1 whereas for discarded it is p ~ 0.48, which appears to be the main difference here. Notably, p < 0.1 would constitute a significant result if tested using a one-sided test (p < 0.05) and we see little reason to perform a two-sided test (because below chance classification is neither hypothesized nor plausible unless confounds affect the cross-validation, see Görgen et al., 2018, Neuroimage). Please provide a detailed explanation of this analysis in the methods section. Furthermore, it is unclear to me what time-points are used for these plots (i.e., what p < 0.1 refers to).

        5. Throughout the manuscript the work by Panichello and Buschman (2021) is cited as evidence for active representations of unattended items found in intracranial recordings. In our view, this study cannot be easily interpreted as such evidence as (1) these items can be discarded after the cue and (2) these signals were only observed until 500 ms after the cue when they were in decline (i.e., there might have been too little time for them to reach baseline), At least in our reading, Panichello and Buschman might not agree with such a characterization of their results.

        Best of Luck

        Thomas Christophel & Vivien Chopurian

          2021 Oct; 19(10): e3001436.
          Published online 2021 Oct 21. 10.1371/journal.pbio.3001436.r005

          Author response to Decision Letter 2

          1 Oct 2021

          Attachment

          Submitted filename:

            2021 Oct; 19(10): e3001436.
            Published online 2021 Oct 21. 10.1371/journal.pbio.3001436.r006

            Decision Letter 3

            Gabriel Gasque, Senior Editor

            4 Oct 2021

            Dear Albert,

            On behalf of my colleagues and the Academic Editor, Frank Tong, I am pleased to say that we can in principle offer to publish your Short Report "Pinging the brain with visual impulses reveals electrically active, not activity-silent working memories" in PLOS Biology, provided you address any remaining formatting and reporting issues. These will be detailed in an email that will follow this letter and that you will usually receive within 2-3 business days, during which time no action is required from you. Please note that we will not be able to formally accept your manuscript and schedule it for publication until you have made the required changes.

            Please take a minute to log into Editorial Manager at http://www.editorialmanager.com/pbiology/, click the "Update My Information" link at the top of the page, and update your user information to ensure an efficient production process.

            PRESS

            We frequently collaborate with press offices. If your institution or institutions have a press office, please notify them about your upcoming paper at this point, to enable them to help maximise its impact. If the press office is planning to promote your findings, we would be grateful if they could coordinate with gro.solp@sserpygoloib. If you have not yet opted out of the early version process, we ask that you notify us immediately of any press plans so that we may do so on your behalf.

            We also ask that you take this opportunity to read our Embargo Policy regarding the discussion, promotion and media coverage of work that is yet to be published by PLOS. As your manuscript is not yet published, it is bound by the conditions of our Embargo Policy. Please be aware that this policy is in place both to ensure that any press coverage of your article is fully substantiated and to provide a direct link between such coverage and the published work. For full details of our Embargo Policy, please visit http://www.plos.org/about/media-inquiries/embargo-policy/.

            Thank you again for choosing PLOS Biology for publication and supporting Open Access publishing. We look forward to publishing your study. 

            Sincerely, 

            Gabriel Gasque, Ph.D. 

            Senior Editor 

            PLOS Biology

            gro.solp@euqsagg


              Articles from PLOS Biology are provided here courtesy of PLOS

              Citations & impact 


              Impact metrics

              Jump to Citations

              Citations of article over time

              Alternative metrics

              Altmetric item for https://www.altmetric.com/details/115565600
              Altmetric
              Discover the attention surrounding your research
              https://www.altmetric.com/details/115565600

              Article citations


              Go to all (14) article citations

              Data 


              Data behind the article

              This data has been text mined from the article, or deposited into data resources.

              Funding 


              Funders who supported this work.

              Departament d&apos;Innovació, Universitats i Empresa, Generalitat de Catalunya (2)

              Fundación Cellex

                Fundação Bial (1)

                Instituto de Salud Carlos III (1)

                Ministerio de Ciencia, Innovación y Universidades (2)

                instituto de salud carlos iii (1)

                ministerio de economía y competitividad (1)