
==== Front
Commun Biol
Commun Biol
Communications Biology
2399-3642
Nature Publishing Group UK London

6780
10.1038/s42003-024-06780-8
Article
Dynamic layer-specific processing in the prefrontal cortex during working memory
http://orcid.org/0009-0003-8682-3920
Degutis Jonas Karolis j.karolis.degutis@maxplanckschools.de

123
http://orcid.org/0000-0003-4532-3084
Chaimow Denis 4
http://orcid.org/0000-0003-2310-5086
Haenelt Daniel 4
Assem Moataz 5
http://orcid.org/0000-0002-9695-2764
Duncan John 5
Haynes John-Dylan 123678
http://orcid.org/0000-0001-5239-1881
Weiskopf Nikolaus 14910
Lorenz Romy romy.lorenz@tuebingen.mpg.de

411
1 https://ror.org/01hcx6992 grid.7468.d 0000 0001 2248 7639 Department of Psychology, Humboldt-Universität zu Berlin, Berlin, Germany
2 grid.484013.a 0000 0004 6879 971X Bernstein Center for Computational Neuroscience Berlin and Berlin Center for Advanced Neuroimaging, Charité Universitätsmedizin Berlin, corporate member of the Freie Universität Berlin, Humboldt-Universität zu Berlin, and Berlin Institute of Health, Berlin, Germany
3 grid.4372.2 0000 0001 2105 1091 Max Planck School of Cognition, Leipzig, Germany
4 https://ror.org/0387jng26 grid.419524.f 0000 0001 0041 5028 Department of Neurophysics, Max Planck Institute for Human Cognitive and Brain Sciences, Leipzig, Germany
5 grid.5335.0 0000000121885934 MRC Cognition and Brain Sciences Unit, University of Cambridge, Cambridge, Cambridgeshire UK
6 https://ror.org/01hcx6992 grid.7468.d 0000 0001 2248 7639 Research Training Group “Extrospection” and Berlin School of Mind and Brain, Humboldt-Universität zu Berlin, Berlin, Germany
7 https://ror.org/03v4gjf40 grid.6734.6 0000 0001 2292 8254 Research Cluster of Excellence “Science of Intelligence”, Technische Universität Berlin, Berlin, Germany
8 https://ror.org/042aqky30 grid.4488.0 0000 0001 2111 7257 Collaborative Research Center “Volition and Cognitive Control”, Technische Universität Dresden, Dresden, Germany
9 https://ror.org/03s7gtk40 grid.9647.c 0000 0004 7669 9786 Felix Bloch Institute for Solid State Physics, Faculty of Physics and Earth Sciences, Leipzig University, Leipzig, Germany
10 grid.83440.3b 0000000121901201 Wellcome Centre for Human Neuroimaging, Institute of Neurology, University College London, London, UK
11 https://ror.org/026nmvv73 grid.419501.8 0000 0001 2183 0052 Max Planck Institute for Biological Cybernetics, Tübingen, Germany
14 9 2024
14 9 2024
2024
7 114022 11 2023
26 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
The dorsolateral prefrontal cortex (dlPFC) is reliably engaged in working memory (WM) and comprises different cytoarchitectonic layers, yet their functional role in human WM is unclear. Here, participants completed a delayed-match-to-sample task while undergoing functional magnetic resonance imaging (fMRI) at ultra-high resolution. We examine layer-specific activity to manipulations in WM load and motor response. Superficial layers exhibit a preferential response to WM load during the delay and retrieval periods of a WM task, indicating a lamina-specific activation of the frontoparietal network. Multivariate patterns encoding WM load in the superficial layer dynamically change across the three periods of the task. Last, superficial and deep layers are non-differentially involved in the motor response, challenging earlier findings of a preferential deep layer activation. Taken together, our results provide new insights into the functional laminar circuitry of the dlPFC during WM and support a dynamic account of dlPFC coding.

Layer-specific imaging of the human dorsolateral prefrontal cortex reveals distinct laminar responses to working memory load and dynamic coding of working memory trial phases.

Subject terms

Cognitive control
Working memory
https://doi.org/10.13039/501100007316 Klaus Tschira Stiftung (Klaus Tschira Foundation) GSO/KT18 Degutis Jonas Karolis issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

The prefrontal cortex (PFC) is critical for a diverse range of higher-level cognitive processes including working memory (WM)1. Early work on WM focused on the neural instantiation of stimulus-specific WM activity in the PFC of non-human primates2,3, with later studies finding content-specific signals across multiple areas of the human cortex4–7. In a parallel line of work, instead of elucidating WM contents, studies examined neural activity related to WM load8–11 (the number of remembered items) and WM manipulation (the reordering of stored information)12. These studies found a positive relationship between PFC activation and task demand for WM as well as other cognitive processes, suggesting that the PFC plays a critical role in the multiple-demand system, a frontoparietal network that commonly responds to a particular diverse range of different cognitively challenging tasks13,14. More recent studies have shown that multivariate activity patterns of the PFC dynamically change depending on WM demands15, with these changes potentially being driven by altered coupling between the PFC and other key brain regions that process WM load16,17.

Despite the PFC’s importance to WM, the role of its laminar circuitry in relation to human cognition remains unclear. Evidence from non-human primates indicates that the PFC has cytoarchitectonic layers that play distinct roles in WM18,19. Superficial layers have been hypothesized to underlie the maintenance of WM20, while deep layers are considered the output layers that send signals away from the PFC towards motor and premotor areas21,22. Recent advances in high spatial resolution functional MRI (fMRI) at ultra-high field strength have enabled the study of lamina-specific responses in human participants in a non-invasive manner with studies predominantly investigating sensory regions23–28, and only recently the left dorsolateral PFC (dlPFC)29. In a double-dissociation, Finn et al. 29. found that superficial layers of the dlPFC were preferentially activated by the manipulation of verbal WM compared to its maintenance, while deep layers responded to a motor response manipulation.

Here, we expand on these human laminar fMRI findings by addressing several open questions about the superficial and deep layers of frontal cortex, while focusing on the aforementioned left dlPFC. First, we manipulate WM load to see whether the superficial layer result from Finn et al. 29. generalizes to other demanding high-level cognitive tasks. Second, we seek to replicate the preferential activation of deep layers to a motor response. Third, drawing from prior research on the impact of WM load on the strength and directionality of coupling between the dlPFC and other brain areas16,17, we explore whether WM load induces layer-specific neural activity pattern alterations by investigating the multivariate code underlying WM load in superficial and deep layers. Fourth, we investigate the stability of this code across the entire duration of the trial, aiming to discern whether well-established dynamic coding properties of the PFC30–34 localize to specific layers.

While dynamic coding more commonly refers to changes of the coding format of WM content4,30,31,34, a broader definition of the concept has been used in the past to investigate dynamic changes in multivariate activity underlying different trial types31 or to characterize different WM control processes at different trial stages15. Here, we follow this broader notion of dynamic coding and employ multivariate decoding of WM load to investigate WM control processes across three periods of the WM task (encoding, delay, retrieval) with the aim to understand whether the multivariate patterns underlying these control processes change in a layer-specific manner.

To preview, we found that WM load preferentially activated superficial layers compared to deep layers during the delay period, and we observed higher load decoding accuracy in superficial layers during the retrieval period. We did not find a significant difference between deep and superficial layers to the motor response manipulation. Finally, we observed that multivariate patterns underlying load at different periods of the WM task changed across time in the superficial layers of the left dlPFC.

Results

Nine participants were scanned at 7 T field strength using a GE-BOLD sequence and used their right hand to complete a WM delayed match-to-sample task, where in two of the four runs we manipulated WM load. Participants had to remember either four or one item and had to always respond during the retrieval period (Fig. 1a). In the other half of the runs the motor response was manipulated: participants had to always remember four items and during the retrieval period they were asked to respond or abstain from responding to a presented probe (Fig. 1a). Participants accurately remembered the stimulus in both the low load (mean = 95.5% correct, standard deviation (SD) = 6.34%) and high load condition (mean = 74.3% correct, SD = 10.6%; Supplementary Fig. 1, left panel) with the high load trials having significantly lower accuracy (p < 0.001, t(8) = 9.44, CI95 = [16.8, 27.7], d = 3.15, two-tailed paired t-test). During the motor manipulation, the participants followed instructions by responding (mean = 99.3% correct, SD = 2.08%) and abstaining to respond (mean = 97.9% correct, SD = 3.13%; Supplementary Fig. 1, right panel) with no significant difference between motor conditions (p = 0.35, t(8) = 1.00, CI95 = [-1.81, 4.59], d = 0.33, two-tailed paired t-test).Fig. 1 Trial design and ROIs.

a Trial design. In load runs (left), participants performed a delayed match-to-sample task, where they had to remember a presented high or low load stimulus (four vs. one item) and had to indicate whether the probe presented at the end of the trial was a part of the stimulus array or not. In motor runs (right), participants were always presented with four items. If an ‘X’ was presented during the probe, they had to abstain from answering. High load stimulus depicts four outdoor scenes. Low load stimulus depicts one outdoor scene and three masks. b Regions of interest. The left dlPFC was defined as four parcels from the frontoparietal network: 8 C, IFJp, IFSa, p9-46v76,83. The left control regions were defined as four parcels from the cingulo-opercular network anatomically adjacent to the left dlPFC: FEF, 6r, 46, FOP576,83 (for full description of ROIs see Methods). c Three equidistant gray matter layers (only superficial and deep layers were used for analyses) defined in the space of the functional images were projected onto an anatomical T1 and functional T2* image of an example participant.

To allow a systematic comparison with the findings reported by Finn et al.29, all our analyses focused on the left dlPFC (Fig. 1b). To assess the specificity of our results, we additionally looked at a set of anatomically adjacent frontal regions in the left hemisphere from the cingulo-opercular network as a control (Fig. 1b). We also report results from the right dlPFC in the Supplementary Material due to its involvement in WM processing and as it was covered in our fMRI acquisition slab. As a control for the right dlPFC, we also report results from right-lateralized frontal regions from the cingulo-opercular network. Please note though, that we a priori decided to specifically investigate the left dlPFC and as such optimized our signal in the left hemisphere (see Methods). We found no significant difference in the number of voxels per layer within a given ROI (Supplementary Fig. 4), allowing comparability of our results across different layers and regions.

All results in the main text are presented for trial periods that are adjusted for a hemodynamic delay of 6 s while “true” trial timings (without taking into account the hemodynamic delay) are additionally indicated in Figs. 2 and 3.Fig. 2 Layer-specific univariate BOLD responses.

a Univariate trial time courses averaged across all voxels from the superficial (top) and deep (bottom) layers in the left dlPFC for load trials (high and low conditions; left plots) and motor trials (response and abstain conditions; right plots). The vertical gray transparent rectangles indicate the two trial periods of interest (adjusted for the hemodynamic delay of 6 s): the delay (11.3–15.1 s) and the retrieval period (20.7–24.5 s). The shaded area indicates standard errors across participants. The horizontal colored transparent rectangles at the bottom of the figure indicate “true” trial periods (without taking hemodynamic delay into account) in the following order: stimulus presentation (light blue), delay period (blue), probe presentation and retrieval period (gray). b The layer-specific load effect (top) and motor effect (bottom) in the superficial and deep layers of the left dlPFC for the delay and retrieval time periods. The load effect was calculated per subject by subtracting the low from the high load time course and averaging time points associated with the trial periods of interest. The motor effect was calculated per subject by subtracting the abstain from the motor response time course and averaging time points in the retrieval period. The stars indicate p < 0.01 (two-tailed paired t-test). Error bars indicate ±SEM. c same as (a) but for left control regions. d Same as (b) but for left control regions.

Fig. 3 Layer-specific decoding accuracy.

a Top: decoding accuracy (linear SVM) of high vs. low load conditions across time from the superficial (gray) and deep layers (red) of the left dlPFC. Bottom: decoding accuracy of response vs. abstain motor conditions across time from the superficial (gray) and deep (red) layers of the dlPFC. Diamonds indicate above-chance decoding significance of superficial and deep layers (in their corresponding color) and comparison between layers (black). Filled diamonds depict significance of p < 0.05, (cluster-permutation test). Hollow diamonds depict p < 0.1, (cluster-permutation test). The shaded rectangles indicate the trial periods in the following order: stimulus presentation (light blue), delay period (blue), probe presentation and retrieval period (gray). “True” trial periods (without accounting for hemodynamic delay) are displayed at the bottom and trial periods adjusted for a 6 s hemodynamic delay are displayed on top. b same as (a) but for left control regions. c Temporal cross-decoding of load from the superficial and deep layers of the dlPFC. Top: superficial layers. Bottom: deep layers. A classifier trained on each time point was tested on itself and all other time points. The matrix diagonal corresponds to the classification accuracy of load in (a). The black outline denotes all above-chance clusters within a layer (nonparametric cluster-permutation test against null, p < 0.05). The red outline signifies all dynamic clusters (conjunction of two nonparametric cluster-permutation tests, p < 0.05; see Methods: Temporal cross-decoding). d Selected columns from (c) visualized in 2D as decoding time courses: accuracy across the trial when trained on the encoding (top plot; 10 s), delay (middle plot; 16 s), and retrieval (bottom plot; 24 s) trial periods and tested on all time points in turn. Columns at 10 s, 16 s, and 24 s in (c) correspond to the top, middle, and bottom time courses, respectively. Transparent gray vertical bars indicate the time points that the classifier was trained on. Shaded area depicts ± SEM. To increase signal-to-noise, all decoding points were averaged with their preceding one in a moving average manner.

Superficial layers preferentially activate to higher WM load during the delay period

Our first goal was to examine the influence of WM load on the superficial and deep layer activity of the left dlPFC (Fig. 1b, c). We predicted that superficial layers would be preferentially influenced by WM load. Trial time courses were estimated for each load condition in each layer in the dlPFC and frontal control regions (Fig. 2a, c). To examine the load effect, we subtracted the low load time course from the high load for a given layer in a given region and examined two trial periods of interest: middle of the delay period (11.3–15.1 s), which was after the encoding peak and before the response probe, and retrieval (20.7–24.5 s; indicated as two gray shaded rectangles in Fig. 2a, c), which was after the presentation of the probe.

We found a higher load effect across participants in superficial as compared to deep layers during the delay period within the left dlPFC (p = 0.0072, t(8) = 3.58, CI95 = [0.040, 0.188], d = 1.19, two-tailed paired t-test), but no significant difference during the response (p = 0.072, t(8) = 2.07, CI95 = [-0.024, 0.46], d = 0.69) (Fig. 2b). In left control regions we found no difference between layers for the load effect in either the delay (p = 0.55, t(8) = 0.630, CI95 = [-0.14, 0.24], d = 0.21) or the retrieval period (p = 0.83, t(8) = -0.223, CI95 = [-0.20, 0.16], d = -0.07) (Fig. 2d). Preferential superficial layer activation to high WM load during the delay and retrieval periods was significant in the right dlPFC (Supplementary Fig. 2b) and right control regions (Supplementary Fig. 2d).

No differential layer activation to motor response in retrieval period

Additionally, we examined the motor effect during the retrieval period (20.7–24.5 s). We expected to replicate higher activation of the deep layers29. Similarly to the load effect, we subtracted the abstain from the response time course for each layer in each ROI (Fig. 2a, c). However, we found no significant difference between superficial and deep layers either in the dlPFC (p = 0.30, t(8) = -1.11, CI95 = [-0.268, 0.094], d = -0.37) (Fig. 2b) or in the control regions (p = 0.18, t(8) = -1.452, CI95 = [-0.20, 0.045], d = -0.48) (Fig. 2d). We also found no significant motor effect difference between layers in the right dlPFC or right control regions (Supplementary Fig. 2b, d).

Superficial layers preferentially code for WM load in retrieval period

Univariate analyses demonstrate changes in average signal between layers; however, individual voxel responses within each layer are not necessarily homogeneous. To be more sensitive to such patterns we used multivariate decoding analyses (for examples in GE-BOLD see refs. 25,35,36) using trial normalization (see Methods) that mitigated the effect of differences in mean univariate response between the conditions which might drive the decoding results. This analysis allowed us to investigate multivariate activation patterns of different layers to load and motor conditions as well as potential differences in response patterns over the course of the trial. Therefore, unlike in the univariate analysis where we compared responses averaged over delay and retrieval periods of the task, here we ran the classification in an exploratory manner across the length of the trial, encompassing the encoding, delay, and retrieval periods.

Previous studies have used multivariate decoding to investigate control processes of WM memory using load15,31,37–40. Here we followed these studies and trained a binary classifier to differentiate between high and low load trials independently for each layer within each ROI and statistically compared decoding results against permutation-derived null distributions as well as between layers. For superficial layers of the left dlPFC, we found two above-chance temporal decoding clusters during the encoding and early delay period (8-12 s, p = 0.0035, dcluster X̄ = 0.86, one-tailed permutation test) as well as at retrieval (22-24 s, p = 0.0027, dcluster X̄ = 1.34, one-tailed permutation test) (Fig. 3a, filled gray diamonds) while two above-chance time points during the late delay failed to reach significance on the temporal cluster-level (16-18 s, p = 0.077, dcluster X̄ = 0.67, one-tailed permutation test) (Fig. 3a, hollow gray diamonds). For deep layers of the dlPFC, two above-chance time points during the encoding and early delay period trended towards temporal cluster significance (8-10 s, p = 0.058, dcluster X̄ = 0.75, one-tailed permutation test; Fig. 3a, hollow red diamonds). For between layer comparison, we found two time points during the encoding and early delay period (10-12 s) where the superficial layer showed stronger WM load decoding than the deep layer, with this time period trending towards significance at a temporal cluster level (p = 0.059, dcluster X̄ = 0.85, two-tailed permutation test; Fig. 3a, black hollow diamonds). Interestingly, we found one cluster that showed a significantly higher decoding of WM load in the superficial compared to the deep layer during the retrieval period (22-24 s, p = 0.023, dcluster X̄ = 1.08, two-tailed permutation test; Fig. 3a, black filled diamonds) even though the bottom-up sensory input at retrieval was completely identical between the two conditions and consisted of a single probed image (Fig. 1a, response probe). These results were not driven by differences in accuracy (Supplementary Fig. 3a) and reaction time (Supplementary Fig. 3b) between low and high load trials. Moreover, we found a significant interaction between the three task periods of encoding (8–12 s), delay (14–20 s), and retrieval (22–24 s) and superficial and deep layers of the left dlPFC (p = 0.034, η2 = 0.35, non-parametric repeated-measures ANOVA).

For the left control regions we found above-chance decoding clusters in the superficial (8–10 s, p = 0.046, dcluster X̄ = 0.77, one-tailed permutation test) and deep layers (10-12 s, p = 0.0043, dcluster X̄ = 0.95, one-tailed permutation test) during the encoding and early delay period (Fig. 3b), but no differences between layers nor an interaction between the three trial periods and layers (p = 0.70, η2 = 0.04, non-parametric repeated-measures ANOVA). We also decoded WM load from the encoding and early delay period in the superficial layers and deep layers of the right dlPFC and found a difference between layers with higher superficial layer decoding, but no interaction between the three trials periods (p = 0.97, η2 = 0.005, non-parametric repeated-measures ANOVA) (Supplementary Fig. 2e). We also decoded load from the right control regions during the encoding and early delay period but found no interaction between the three trial periods and layers (p = 0.083, η2 = 0.27, non-parametric repeated-measures ANOVA) (Supplementary Fig. 2f).

Both superficial and deep layers non-preferentially code for motor response in retrieval period

To detect patterns associated with the motor effect, we trained a binary classifier to distinguish between response and abstain trials. Since the participants only found out about the type of trial during the presentation of the probe, as expected, there was no significant decoding prior, either in the left dlPFC or left control regions. We found a single above-chance cluster each in the superficial (p = 0.0054, dcluster X̄ = 1.27, one-tailed permutation test) and deep layers (p = 0.015, dcluster X̄ = 1.18, one-tailed permutation test) during the retrieval period (24-26 s; Fig. 3a, gray and red diamonds). However, there was no significant difference between the two layers. These results agree with our univariate findings showing no difference between superficial and deep layers during the motor retrieval period. The right dlPFC and right control region showed similar motor decoding results. Motor trial decoding was significant in both layers of the right dlPFC and right control regions during the retrieval period (Supplementary Fig. 2f).

Dynamic coding of WM load in superficial layers

We sought to determine whether the multivariate pattern sustaining load information is stable across time. To do so, we used temporal cross-decoding, as previously done to investigate the stability of WM content in humans and non-human primates30,31,41,42. We trained on a given time point and tested on all time points in turn. This analysis provided us with two distinct types of elements within the cross-decoding matrix: the diagonal, which indicates the amount of information coding for load within a given time point, and the off-diagonal, which provides the generalization between two time points. If the trained classifier generalizes to other time points (i.e., high decoding accuracies of the off-diagonal elements of the decoding matrix), the analysis reveals that the multivariate pattern sustaining load information (i.e., “coding”) is stable across time. However, when the diagonal of the matrix, representing training and testing on the same time points, shows higher decoding accuracies than the off-diagonal elements, then we can conclude from this lack of cross-generalization of the decoder, that the multivariate activity is dynamic and varies across time. We examined the generalization of multivariate patterns in response to load in the superficial and deep layers of the dlPFC.

Within the superficial layers of the dlPFC, we identified two above-chance decoding clusters (Fig. 3c, black dashed line). The first decoding cluster spanned both the encoding and delay period (p = 0.001, dcluster X̄ = 1.04, one-tailed permutation test), while the second cluster was isolated to the retrieval period (23–25 s; p = 0.003, dcluster X̄ = 1.15, one-tailed permutation test). Looking more closely into the first decoding cluster, we observed that the classifier trained on the encoding period (8–10 s) generalized to the middle of the delay period (16 s), and the classifier trained on the delay period (14–18 s) generalized to the late delay period (up to 22 s). Importantly though, the classifier trained on the encoding period did not generalize to the late delay period and vice versa, indicating a differential multivariate pattern sustaining load information in those two periods of the trial. To visualize these results, we plotted the columns from Fig. 3c as time courses (Fig. 3d, top panel), which depict that the classifier trained on 10 s had high decoding during the encoding and early delay, but the decoding accuracy waned in the second half of the trial. Similarly, the classifier trained on 16 s had the highest decoding during the late delay period (Fig. 3d, middle panel). With respect to the second decoding cluster identified at retrieval, we observe that it only generalizes to its surrounding time points, indicating that the multivariate pattern reflecting load during the retrieval period was distinct from encoding and delay (bottom time course in Fig. 3d).

Within the deep layers of the dlPFC, a single above-chance cluster (Fig. 3c, black dashed line, p = 0.002, dcluster X̄ = 1.10, one-tailed permutation test) was found during the encoding period (8-10 s) which generalized to the delay (18 s). Plotting the column of Fig. 3c as a time course (Fig. 3d, top panel), we see that the deep layer results follow a similar trend as the superficial layer results: the classifier trained on 10 s had high decoding during the encoding and early delay, but the decoding accuracy waned in the second half of the trial. Similarly, the classifier trained on 16 s had the highest decoding during the late delay period (Fig. 3d, middle panel).

To empirically check the lack of generalization between the above-chance clusters, we identified dynamic elements within the cross-decoding matrices. Dynamic elements were calculated as a conjunction of above-chance elements surviving two cluster permutation tests (see Methods for details). Three clusters were found in the superficial layers of the left dlPFC: two off-diagonal elements at the intersection of the encoding and early delay (train) and retrieval (test) periods, two elements at the intersection of the middle delay (train) and encoding and early delay (test) periods, and four elements at the intersection of retrieval (train) and encoding and early delay (test) periods (Fig. 3c).

To summarize, our cross-decoding results showed three differential multivariate response patterns sustaining load in the superficial layers during the encoding, delay, and retrieval periods, indicative of dynamic coding between trial periods.

Discussion

We scanned participants at ultra-high field strength performing a delayed match-to-sample task where we manipulated WM load and motor response while measuring responses from the superficial and deep layers of the dlPFC. In the left dlPFC, we found a higher superficial compared to deep layer activation in response to increased working memory (WM) load during the delay period; however, we did not find a differential laminar response to a motor manipulation. To examine the multivariate response across the length of a trial, we decoded load and motor trial types from the layers of the left dlPFC. We found that superficial layers preferentially code for difference in WM load during the retrieval period. With respect to the motor trials, we found that both superficial and deep layers code for the motor response in the retrieval period but failed to find any significant laminar difference. Finally, we also found that the multivariate code underlying WM load in the superficial layers dynamically changed between the encoding, delay, and retrieval periods.

Our results replicate early non-human primate studies showing persistent spiking in the lateral PFC in all periods of a WM task2 by showing sustained univariate activity in both superficial and deep layers during high load trials. Our work also corroborates findings showing higher activation of the lateral PFC to higher task demands, as previously seen in univariate human fMRI work on WM load8,9 and the manipulation of WM information12. Additionally, our results expand on previous WM literature that has used multivariate decoding to investigate neural mechanisms of WM5,6,41–44.

We also provide insights into the mesoscale functional architecture of the dlPFC in response to WM load. Greater neural processing associated with WM load in superficial layers compared to deep throughout the delay and retrieval periods might indicate activity related to more general task demands; these demands increase as a result of higher WM load, possibly by both recruiting new neural populations (e.g., in retrieval) and by increasing univariate activity (e.g., in delay). These WM control processes might correspond to the updating of content during encoding, the monitoring or contents throughout the delay period45, or the compression or reorganization of WM contents as a result of increased load demand46. Interestingly, there was a preferential univariate response of the superficial compared to deep layers to high load during the delay period, but no laminar difference in decoding. Superficial layers might utilize the same neural population underlying one or several control processes to a higher degree. This would result in higher univariate activity, but lower decodability since univariate contributions were minimized through trial normalization. Taking into account results from Finn et al.29., preferential response in the superficial layers to both WM manipulation and higher WM load might point towards a lamina-specific (i.e., superficial layer) activation of the frontoparietal multiple demand system13 for heightened task demands more generally. However, it is essential to underscore the necessity for future investigations to systematically examine the layer-specific activation of the frontoparietal network across a diverse spectrum of cognitive tasks with varying demands; as can be accomplished employing novel real-time fMRI approaches utilizing neuroadaptive Bayesian optimization47–49.

Previous studies have established that recurrent connections within superficial layers of the PFC support sensory input-specific gamma-band activity, while deep layer alpha and beta oscillations modulate this superficial layer activity18,33. In turn, the interaction between deep and superficial layers might act as a gating mechanism for sensory information to enter the PFC and WM storage50,51. In our study, both superficial layers and deep layers had above-chance temporal cross-decoding clusters during the encoding of WM information. Deep layer coding for WM load during the encoding period might correspond to the opening of the gate, which allows for the encoding of new WM information within the superficial layers. Based on our results, we can only speculate whether deep layers drive superficial layer activity. It should be noted that deep layer coding was significant in the right dlPFC and only trended towards significance in our temporal cluster analysis of the left dlPFC. This might indicate a lateralization in the WM encoding control process to the right hemisphere52; however, further studies are needed to corroborate these results and unravel any causal laminar relationships in both hemispheres.

For the retrieval period, we found significantly higher superficial compared to deep layer decoding of WM load which might indicate increased demands during the probe’s comparison to the stored memory (1 vs. 4 items) despite the presentation of identical response probes in low and high WM load conditions. The comparison process seems to happen within the superficial layers, which might multiplex both the novel sensory probe stimulus and the remembered memoranda; to achieve this at high load, additional neural populations within the superficial layers might be recruited, thus driving the multivariate decoding. Superficial layer decoding relates to previous findings showing ramping gamma-band activity at the end of memory delays53, as a consequence of its disinhibition resulting from the drop of alpha/beta power to allow a successful read-out of information from WM54. Intriguingly, load decoding peaked prior to motor decoding, which might denote the sequence of computations in the dlPFC: a comparison of the probe to the WM set before the initialization of a motor response.

Our results show multivariate laminar activity patterns related to load that are specific to the dlPFC. For the left dlPFC, we find that superficial layers show a significant decoding during the encoding and retrieval period and trend towards significance during the delay period. We find a higher superficial compared to deep layer decoding in the response period as well as trending towards significance in the encoding period. When comparing this to the left control regions, we only find above-chance decoding during the encoding period, but no effects for other trial periods. More importantly, no layer-specific results were identified in the left control region. We also find a significant interaction effect between the three trial periods and layers in the left dlPFC. The results from the right dlPFC show a preferential involvement of superficial layers during the encoding period. Despite the similar result seen in the right control regions, the laminar difference emerges before the onset of the encoding period and seems to be driven by the single negative time point in deep layers, which likely negates its significance.

Our findings of WM load-induced layer-specific changes in neural activity patterns in the left dlPFC can be contextualized by previous research in non-human primates and human fMRI studies. More specifically, it has been shown that increased WM load changes the strength and directionality of neuronal coupling between the prefrontal cortex, the frontal eye fields, and lateral intraparietal cortex16, and influences both between-network and within-network coupling among frontoparietal, ventral attention, and default mode networks17, which are reflected in a change in multivariate pattern. This altered network coupling in the left dlPFC may result in subtle variations in neural activity patterns when WM demands change, which has been observed before15 and which we localize to the superficial layers of the dlPFC. This observation suggests a network-coding perspective of WM and aligns with the established intra-prefrontal and cortico-cortical connections of the superficial layers in the prefrontal cortex18,33. In addition, changes in the multivariate patterns of the superficial layers in response to WM load may also be attributed to activation spread beyond initially activated regions, potentially recruiting new neural populations, although such signal spread might also manifest in univariate activation differences. Equally, the multivariate pattern in the dlPFC could suggest content-independent pointers38 that project to parietal and sensory cortices where the maintained memoranda are thought to be stored7,55,56.

The canonical model of the WM suggests that that information coded within the PFC is stable and firing is persistent across the delay period3,57. However, evidence has emerged that WM is both rhythmical and non-persistent in activation and the format of the coded information dynamically varies across time33,58. Studies using temporal cross-decoding have observed a dynamic code between the encoding period and delay activity in the PFC of non-human primates and humans30,31,34, and across the visual hierarchy in humans4,42. We find three non-generalizing sub-clusters within the span of the trial; each sub-cluster corresponding to either the encoding, delay, or retrieval periods. Previous studies argued that dynamic coding indicates a point where transient sensory stimuli are transformed into a more stable representation30; however, since we were decoding load and not the contents of WM, the decoding sub-clusters might underlie separate WM control processes that occur on the WM content within the same layer.

In contrast to superficial layers of the PFC, the deep layers are considered to be the output - or top-down - layers of the cortical area and have neurons modulated by D2 dopamine receptors21,22. Finn et al.29. found that a motor modulation only activated the deep layers of the prefrontal cortex, while superficial layers stayed at baseline. We do not find a preferential activation of either superficial or deep layers of the left dlPFC to the motor manipulation either in univariate activity or in multivariate decoding. The discrepancy between Finn et al.29. and our results may be related to a difference in paradigm, a difference in scanning sequences, or a smaller sample size.

Finn et al.29. contrasted trials in which participants either had to respond after the presentation of a probe or had to not respond when a probe was not presented. This resulted in a change of two factors between conditions: the presence or absence of the probe stimulus and the motor response. Our paradigm always included the presentation of the probe and participants were asked to either respond or abstain from responding based on a presented cue. If the deep layer modulates the information that enters the superficial layer of the dlPFC, then that computation could provide an additional deep layer response. Thus, higher activation of the deep layers in Finn et al.29. might not be a result of the output, but rather of the opening of the gate by the deep layers to let in the probe for subsequent comparison in the superficial layers. Hence, this probe input signal would have been absent when contrasting between conditions in our paradigm, since participants could not predict whether a given trial required a response or not, they would have maintained the memoranda and automatically compared it to the presented probe. Based on our results we might conclude that the superficial and deep layers are both non-differentially involved in the output of a motor signal. Further research is needed to distinguish between the input and output hypotheses by possibly extending the time between the presentation of the probe and motor initiation.

An additional difference between the two studies was the pulse sequence used to measure signals from the layers. We used conventional GE-BOLD, while Finn et al.29. used both BOLD and cerebral blood volume (CBV), which resulted in an SS-SI-vascular space occupancy (VASO) contrast59,60. The BOLD signal has poorer spatial specificity, but higher sensitivity, as a larger part of the signal comes from larger veins present in the pial surface61 resulting in a broader point-spread function of the signal in both columns62 and layers63 compared to more spatially precise scanning sequences. VASO has a smaller contrast-to-noise ratio, but localizes activity to both superficial and deep layers, with less bias from large draining veins59. The main concern in GE-BOLD are interdependencies in the signal between different layers of the cortex, as a result of the blood draining from deep layers towards superficial layers64,65. The superficial layers might thus contain signals from the deep layers as well. One strategy of mitigating the effects of draining veins is to compare differences of tasks, as in the case of our study (high load compared to low load and response compared to abstain trials) and as done in previous GE-BOLD studies; a subtraction of conditions would theoretically remove a linear draining effect23,24,26,28,66. Despite employing this method, both superficial and deep layers activated to the motor manipulation. Thus, the linear detrend might have not completely removed the draining vein signal, thus future work can employ more sophisticated models to account for the superficial layer blurring63,67; however, preferential deep layer activation has been previously observed in GE-BOLD studies that have not employed these models23,24,26,68,69.

In summary, we show that superficial layers of the dlPFC are not only crucial for WM maintenance, but in fact are involved in various WM subprocesses by dynamically adapting to current task demands. In addition, we highlight that deep layers seem to have a more complex role in WM than previously understood. Our findings offer new insights into the neural mechanisms of WM and try to bridge the gap between the human and non-human primate WM literature.

Methods

Participants

We scanned nine participants (ages 23-36 years, three female) that were recruited from the Max Planck Institute for Human Cognitive and Brain Sciences (MPI CBS) subject database. The sample size was comparable to previous laminar fMRI studies25,66,70. All participants gave written informed consent and received monetary compensation for their participation. The study was approved by the ethics committee of the University of Leipzig (441/20-ek). All ethical regulations relevant to human research participants were followed.

Stimuli and procedure

Following the design of Finn et. al.29, participants completed four runs of a delayed match-to-sample task (Fig. 1a). In two runs we manipulated the task difficulty by changing the WM load (from now on referred to as the load runs). In the other two runs we manipulated the motor response during the retrieval period (from now on referred to as the motor runs). Each trial began with the presentation of the sample for 3.5 s which consisted of a four-square array where each square was equidistant from the center. A fixation cross was presented centrally. In the load runs either one (low load) or four (high load) items were presented within the squares as the to-be-remembered sample. The items consisted of faces (subcategory: male and female) or scenes (subcategory: indoor and outdoor). The face and scene stimuli were retrieved from the Face-Place database71. Faces had neutral emotional expressions. The scenes included indoor and outdoor pictures of houses. All stimuli were converted to grayscale. For each trial, the category and subcategory of the items within the array were the same (e.g., only female faces). In the case of a low load condition, the remaining three squares consisted of Fourier-scrambled items from the same category, in order to preserve the low-level visual features of the overall sample array. Visual masks consisting of Fourier-scrambled images were presented for 200 ms at the location of the four squares after sample presentation (not shown in Fig. 1a). All scrambled images were from the same category as the sample. Motor runs always had four items (i.e., high load) as the sample.

A WM delay period of 13 s followed during which the participants had to remember the sample. At the end of the delay period a probe image (same size as one of the sample squares) appeared centrally with a question mark placed above for 1.2 s. The probe was always from the same category and subcategory as the sample. Participants had 4.2 s to respond using their index and middle finger of their right hand to indicate whether the presented probe was part of the remembered sample. During motor runs participants had to maintain high-load stimuli for each trial (i.e., high-load WM maintenance, no load manipulation) while we manipulated whether participants had to make a motor response (response trial) or had to abstain from their motor response (abstain trials) during the retrieval period. Identical to the load runs, during the response trials in the motor runs participants had to indicate using their index and middle fingers whether the presented probe was part of the remembered sample. During an abstain trial the probe was still presented; however, an ‘X’ was placed above instead of a question mark cueing the participants to abstain from responding. If a subject responded during this period on an abstain trial, this was marked as a miss-trial. An inter-trial-interval of 10 s followed the retrieval period.

Stimuli (across both sample and probe) did not repeat within a run. The types of trials within a run were counterbalanced. Categories were counterbalanced for high and low load for each load run. The position of the single item was counterbalanced across the four positions in the low load trials for each load run. The sequence of trials was pseudorandomized for each run; a maximum of three trials of the same type could be subsequently presented. Each trial was 30.7 seconds long. A run had 16 trials and lasted 8 min and 24 s.

Prior to entering the scanner, participants trained by completing eight load trials. To proceed to the scanner participants had to respond correctly to at least 70% of the trials, otherwise, the practice load run was repeated.

fMRI protocol

MRI data were acquired using a Siemens MAGNETOM Terra 7 T MRI system (Siemens, Erlangen, Germany) using an 8Tx/32Rx head coil (Nova Medical Inc., Wilmington, MA, USA). Anatomical images were acquired using an MP2RAGE sequence (Marques et al., 2010) (TR = 5000 ms, TE = 2.27 ms, 0.75 mm isotropic voxels) at two inversion times (TI of 900 ms, 2700 ms with a flip angle of 3°, 5°, respectively) that were combined to yield a T1-weighted image.

High-resolution functional data were acquired using a T2*-weighted 2D gradient-echo EPI sequence (TR = 2000 ms, TE = 23 ms, 0.8 mm isotropic voxels, flip angle = 80°, 35 slices). Prior to running the high-resolution scans, a whole-brain functional localizer scan using a T2*-weighted 2D gradient-echo EPI sequence (TR = 2000 ms, TE = 18 ms, 2 mm isotropic voxels, flip angle = 75°) was acquired and analyzed to subsequently position the high-resolution partial-brain slab. During this localizer run participants performed an additional load run (see Stimuli and procedure). The resulting peak activation defined by the contrast between delay activity and baseline informed the positioning of the high-resolution slab around the expected anatomical location of the dlPFC.

The procedure of the scanning was the following: participants were first scanned using the low-resolution localizer after which the anatomical image was acquired. After semi-manual shimming, participants performed two load runs followed by two motor runs (except for S01 who performed load and motor runs interleaved).

Preprocessing and co-registration

To obtain accurate reconstructions of gray and white matter surfaces from structural MP2RAGE data, a Freesurfer-based pipeline (Chaimow et al. in preparation,72) was applied. First, a T1-weighted image that is spatially homogeneous yet free of extracerebral noise was generated by bias-correcting the INV2 MP2RAGE image before multiplying with the UNI image (MPRAGEize)73. Next, a high-quality brainmask was computed by passing this T1-weighted image to CAT1274 for segmentation and merging the resulting gray and white matter components. Finally, the T1-weighted image was processed using Freesurfer’s recon-all pipeline (version 7.1) using the high-resolution flag and substituting Freesurfer’s auto-generated brain mask with the CAT12-derived brain mask.

For subsequent ROI selection (see Definition of ROI), the segmented structural surfaces of each participant were brought into fsLR164k CIFTI space using ciftify_recon_all, which is part of the Python-based ciftify package75 with expert settings to register to the high-resolution 0.5 mm FSL MNI152_T1 template. Ciftify adapts the post-Freesurfer portion of the Human Connectome Project’s (HCP) minimal preprocessing pipeline to non-HCP acquired data76 and includes surface based alignment (MSMSulc)77 to HCPs fs_LR space.

Registration between the Freesurfer processed T1-weighted image and the mean image of the third functional GE-BOLD run was calculated using ANTs78. The structural image was initially registered using a linear affine transform and subsequently a symmetric normalization (SyN) algorithm generated warp field for distortion correction79.

Functional scans were motion corrected using AFNI’s 3dVolreg using Fourier interpolation to the volume which had the minimum outlier fraction of voxels as calculated by 3dToutcount across all runs (as per standard AFNI preprocessing)80.

We calculated three equidistant depth bins across the anatomical cortical ribbon in functional space using LAYNII81. For smoother layer estimation, prior to layer generation, the anatomical ribbon was upsampled by a factor of five. The upsampled ribbon was then used to generate three equidistant layers. Finally, the layers were downsampled back to the functional images’ 0.8 mm-isotropic resolution. Even though we were interested in only two layers (superficial and deep), the middle layer was estimated as a ‘buffer’ between the deepest and most superficial depths, thus minimizing partial volume effects by excluding voxels where there might have been a large overlap between superficial and deep layers. We applied the affine transformation matrix and warp field generated by ANTs to register the layers to functional space.

Definition of ROI

We generated two regions of interest using the Human Connectome Project Multi-Modal Parcellation version 1.0 (HCP MMP 1.0) atlas76. This surface-based cortical parcellation has been generated using high-resolution, multi-modal data from 210 participants. Surface parcels were transformed from fsLR space to single subject functional space following a combined transformation, estimated by structural-functional registration and MSMSulc77 surface-based alignment to atlas space (fsLR) in each individual subject. This subject-specific surface-based approach has been shown to substantially improve cortical area localization compared to traditional volumetric approaches and is also superior to other surface-based approaches82. We focused all primary analyses on the left dlPFC, following Finn el al.29. We defined the left dlPFC by selecting four left-lateralized frontal parcels (8 C, IFJp, IFSa, p9-46v) for which we had sufficient coverage and that correspond to the frontoparietal network83. To assess whether our findings were specific to the dlPFC, we selected an anatomically and functionally similar control set of four left-lateralized frontal parcels (FEF, 6r, 46, FOP5) from the cingulo-opercular network (COP)83. As an exploratory analysis, we also looked at the dlPFC and the control region in the right frontal cortex, which were defined using the same corresponding parcels in the right hemisphere.

We focused our main analyses on the left dlPFC as this followed Finn et al.29. allowing us to have a direct comparison of our results to theirs. Since we decided a priori to focus on the left hemisphere, we also aimed to minimize the distortion of the signal, specifically in that area. The results from the right dlPFC are presented for transparency, as they were part of the acquisition slab and are involved in WM processing.

Univariate analysis

For each layer, trial time courses were extracted by finite-impulse response modeling using AFNI’s 3dDeconvolve with the ‘TENTzero’ basis function model. A total of 17 regressors were fit to 32 s after the onset of the stimulus. The first and last were zero basis functions, while the middle fifteen free basis functions. We used ‘TENTzero’ to remove any possible pre-stimulus differences between trial types. In addition, we fit polynomial nuisance regressors up to the fifth order to detrend the signal. One time course was estimated for each trial-type in each type of run (high and low load from the load runs; response and abstain from the motor runs). We used the constant predictor to calculate the percent signal change for each voxel in each trial-type.

Multivariate decoding

We performed two decoding analyses: high vs. low load (load decoder) using the load runs, and response vs. abstain (motor decoder) using the motor runs. Decoding analyses were run separately on the superficial and deep layers of the left dlPFC and on the left frontal control regions and corresponding ROIs in the right hemisphere.

We ran the classification across the entire duration of a trial. Time courses of each voxel were high-pass filtered and then temporally z-scored for each run independently. Since the onset of the trial was not TR-locked, we chose the nearest rounded-down TR closest to the onset of the trial and the following 16 TRs (32 s) as a single trial (resulting in a total of 17 TRs). Trial normalization and feature z-scoring was applied to training and testing data independently, thus reducing the classifier’s ability to decode based on overall univariate activity differences between classes. Additionally, to increase the signal-to-noise ratio, we averaged the data from two adjacent TRs for each decoding time point. We used a linear support vector machine classifier (LSVMs84; libsvm: www.csie.ntu.edu.tw/~cjlin/libsvm/) with a cost parameter of c = 1.

Within a cross-validation fold, we trained on the same number of trials in order to account for the trial transition matrix, which corrected for possible carry-over effects from the previous trial. For example, in the motor decoding, we trained on the same number of retrieval trials where the previous trial was a retrieval or an abstain trial, thus creating four types of trials: retrieval-abstain, retrieval-retrieval, abstain-retrieval, and abstain-abstain. Despite the four types of trials, we ran a binary classification where the current trial determined the label name. The trials chosen for a given fold were randomly selected from all available trials across both runs. The number of cross-validation folds was determined by taking the smallest number of trial types and multiplying them by four. Despite some trials being trained on more than once across folds, each fold had a unique combination of the trials trained and tested on.

Temporal cross-decoding

To investigate the temporal stability of the multivariate decoding across the delay period, we ran a temporal cross-decoding classification analysis where we decoded the load condition from the superficial and deep layers of the dlPFC. This analysis followed the same classification procedure as the temporal decoding analysis, except each classifier trained on a given time point was tested on data from all time points in turn. This analysis allowed us to see whether the multivariate code at one time point generalizes to other time points and thus to determine whether the underlying code is dynamic or stable. There was no informational overlap between the training and testing data, since the same trials were never included in the two subsets, both when training and testing on the same and different time points.

Statistics and Reproducibility

Statistical inference for the multivariate analyses was done on empirically generated permuted null distributions. For each decoding analysis, time point, and layer a null distribution was generated by running the same procedure of classification on training data with permuted labels. The analysis was run 250 times for each participant generating a classification accuracy on each iteration. A population null distribution (for each analysis, time point, and layer) was generated by drawing a single sample from each subject’s null distribution (resulting in nine samples) and taking the average resulting in a mean t-value. This was done 10,000 times. Significance was then determined using a cluster-permutation test85 where the summed t-value of the empirical data clusters was compared to the summed t-value of clusters within the generated null distribution. One-tailed tests were run to test for above-chance decoding significance, while two-tailed tests were used for the comparison between layers.

Above-chance decoding of the temporal cross-decoding matrix was tested using a cluster-permutation test85. We calculated the summed t-value of all above-chance elements within a cluster. Within a matrix, there could have been more than one cluster. A null distribution was generated by a sign permutation test by randomly sampling one of the decoding accuracies within a given element of the matrix. This was run 10,000 times. The summed t-value of the largest above-chance cluster within the permuted data was taken as a statistic for the null distribution. The empirical clusters were then compared to the distribution; all clusters where p < 0.05 were deemed to be above-chance.

In addition to calculating an above-chance cluster, we also identified dynamic elements within the cross-decoding matrix. Following previous research30,42, dynamic elements were defined as off-diagonal elements that had lower decoding accuracies than their two corresponding diagonal elements (e.g., aij < aii ∧ aij < ajj). Two cluster-permutation tests were run following the above-chance cluster sign permutation method; in this case all aij elements were initially subtracted by aii in one test and ajj in the second test. Only elements that were within clusters from both tests were deemed dynamic. Also, the diagonal entries corresponding to a dynamic element had to both be within an above-chance cluster.

Supplementary information

Peer review file

Supplementary material

Supplementary information

The online version contains supplementary material available at 10.1038/s42003-024-06780-8.

Acknowledgements

J.K.D. was supported by the Max Planck Society and BMBF. R.L. was funded by the Wellcome Trust (209139/Z/17/Z). M.A. and J.D. were funded by the UK Medical Research Council (MC_UU_00030/7). J.D.H. was supported by the Deutsche Forschungsgemeinschaft (DFG, Exzellenzcluster Science of Intelligence). N.W. was supported by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) – project no. 347592254 (WE 5046/4-2); the Federal Ministry of Education and Research (BMBF) under support code 01ED2210; the European Union’s Horizon 2020 research and innovation programme under the grant agreement No 681094. This project also received funding from the Klaus Tschira Stiftung gGmbH foundation (support code GSO/KT18). Stimulus images courtesy of Michael J. Tarr, Center for the Neural Basis of Cognition and Department of Psychology, Carnegie Mellon University, http://www.tarrlab.org/. We thank Rob Mok, Karla Matic, and Laurentius Huber for valuable discussions. The authors acknowledge support by the Open Access Publication Fund of Humboldt-Universität zu Berlin.

Author contributions

J.K.D., R.L, and D.C. conceived the study; J.K.D. analyzed the data; J.K.D., D.C., D.H., M.A., J.D, and R.L. developed methods; J.K.D. visualized the data; R.L., N.W., and J-D.H, supervised the work; J.K.D. and R.L. drafted the manuscript. All authors reviewed the last version of the manuscript.

Peer review

Peer review information

Communications biology thanks the anonymous reviewers for their contribution to the peer review of this work. Primary Handling Editors: Michel Thiebaut de Schotten and Benjamin Bessieres. A peer review file is available.

Funding

Open Access funding enabled and organized by Projekt DEAL.

Data availability

Preprocessed data for all participants and corresponding analysis code are available on the OSF platform (https://osf.io/dtkne/).

Competing interests

The authors declare the following competing interests: The Max Planck Institute for Human Cognitive and Brain Sciences and Wellcome Centre for Human Neuroimaging have institutional research agreements with Siemens Healthcare. N.W. holds a patent on acquisition of MRI data during spoiler gradients (US 10,401,453 B2). N.W. was a speaker at an event organized by Siemens Healthcare and was reimbursed for the travel expenses.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Curtis CE D’Esposito M Persistent activity in the prefrontal cortex during working memory Trends Cogn. Sci. 2003 7 415 423 10.1016/S1364-6613(03)00197-9 12963473
Curtis, C. E. & D’Esposito, M. Persistent activity in the prefrontal cortex during working memory. Trends Cogn. Sci. 7, 415–423 (2003).12963473 10.1016/S1364-6613(03)00197-9
2. Funahashi S Bruce CJ Goldman-Rakic PS Mnemonic coding of visual space in the monkey’s dorsolateral prefrontal cortex J. Neurophysiol. 1989 61 331 349 10.1152/jn.1989.61.2.331 2918358
Funahashi, S., Bruce, C. J. & Goldman-Rakic, P. S. Mnemonic coding of visual space in the monkey’s dorsolateral prefrontal cortex. J. Neurophysiol. 61, 331–349 (1989).2918358 10.1152/jn.1989.61.2.331
3. Fuster JM Alexander GE Neuron activity related to short-term memory Science 1971 173 652 654 10.1126/science.173.3997.652 4998337
Fuster, J. M. & Alexander, G. E. Neuron activity related to short-term memory. Science 173, 652–654 (1971).4998337 10.1126/science.173.3997.652
4. Sreenivasan KK Vytlacil J D’Esposito M Distributed and dynamic storage of working memory stimulus information in extrastriate Cortex. J. Cogn. Neurosci. 2014 26 1141 1153 10.1162/jocn_a_00556 24392897
Sreenivasan, K. K., Vytlacil, J. & D’Esposito, M. Distributed and dynamic storage of working memory stimulus information in extrastriate. Cortex. J. Cogn. Neurosci. 26, 1141–1153 (2014).24392897 10.1162/jocn_a_00556
5. Harrison SA Tong F Decoding reveals the contents of visual working memory in early visual areas Nature 2009 458 632 635 10.1038/nature07832 19225460
Harrison, S. A. & Tong, F. Decoding reveals the contents of visual working memory in early visual areas. Nature 458, 632–635 (2009).19225460 10.1038/nature07832
6. Christophel TB Hebart MN Haynes J-D Decoding the contents of visual short-term memory from human visual and parietal cortex J. Neurosci. 2012 32 12983 12989 10.1523/JNEUROSCI.0184-12.2012 22993415
Christophel, T. B., Hebart, M. N. & Haynes, J.-D. Decoding the contents of visual short-term memory from human visual and parietal cortex. J. Neurosci. 32, 12983–12989 (2012).22993415 10.1523/JNEUROSCI.0184-12.2012
7. Christophel TB Klink PC Spitzer B Roelfsema PR Haynes J-D The distributed nature of working memory Trends Cogn. Sci. 2017 21 111 124 10.1016/j.tics.2016.12.007 28063661
Christophel, T. B., Klink, P. C., Spitzer, B., Roelfsema, P. R. & Haynes, J.-D. The distributed nature of working memory. Trends Cogn. Sci. 21, 111–124 (2017).28063661 10.1016/j.tics.2016.12.007
8. Druzgal TJ D’Esposito M Dissecting contributions of prefrontal cortex and fusiform face area to face working memory J. Cogn. Neurosci. 2003 15 771 784 10.1162/089892903322370708 14511531
Druzgal, T. J. & D’Esposito, M. Dissecting contributions of prefrontal cortex and fusiform face area to face working memory. J. Cogn. Neurosci. 15, 771–784 (2003).14511531 10.1162/089892903322370708
9. Emrich SM Riggall AC LaRocque JJ Postle BR Distributed patterns of activity in sensory cortex reflect the precision of multiple items maintained in visual short-term memory J. Neurosci. 2013 33 6516 6523 10.1523/JNEUROSCI.5732-12.2013 23575849
Emrich, S. M., Riggall, A. C., LaRocque, J. J. & Postle, B. R. Distributed patterns of activity in sensory cortex reflect the precision of multiple items maintained in visual short-term memory. J. Neurosci. 33, 6516–6523 (2013).23575849 10.1523/JNEUROSCI.5732-12.2013
10. Todd JJ Marois R Capacity limit of visual short-term memory in human posterior parietal cortex Nature 2004 428 751 754 10.1038/nature02466 15085133
Todd, J. J. & Marois, R. Capacity limit of visual short-term memory in human posterior parietal cortex. Nature 428, 751–754 (2004).15085133 10.1038/nature02466
11. Ranganath C DeGutis J D’Esposito M Category-specific modulation of inferior temporal activity during working memory encoding and maintenance Cogn. Brain Res. 2004 20 37 45 10.1016/j.cogbrainres.2003.11.017
Ranganath, C., DeGutis, J. & D’Esposito, M. Category-specific modulation of inferior temporal activity during working memory encoding and maintenance. Cogn. Brain Res. 20, 37–45 (2004).10.1016/j.cogbrainres.2003.11.017
12. D’Esposito M Postle BR Ballard D Lease J Maintenance versus manipulation of information held in working memory: an event-related fMRI study Brain Cogn 1999 41 66 86 10.1006/brcg.1999.1096 10536086
D’Esposito, M., Postle, B. R., Ballard, D. & Lease, J. Maintenance versus manipulation of information held in working memory: an event-related fMRI study. Brain Cogn 41, 66–86 (1999).10536086 10.1006/brcg.1999.1096
13. Duncan J The Structure of Cognition: Attentional Episodes in Mind and Brain Neuron 2013 80 35 50 10.1016/j.neuron.2013.09.015 24094101
Duncan, J. The Structure of Cognition: Attentional Episodes in Mind and Brain. Neuron 80, 35–50 (2013).24094101 10.1016/j.neuron.2013.09.015
14. Duncan J The multiple-demand (MD) system of the primate brain: mental programs for intelligent behaviour Trends Cogn. Sci. 2010 14 172 179 10.1016/j.tics.2010.01.004 20171926
Duncan, J. The multiple-demand (MD) system of the primate brain: mental programs for intelligent behaviour. Trends Cogn. Sci. 14, 172–179 (2010).20171926 10.1016/j.tics.2010.01.004
15. Soreq E Leech R Hampshire A Dynamic network coding of working-memory domains and working-memory processes Nat. Commun. 2019 10 936 10.1038/s41467-019-08840-8 30804436
Soreq, E., Leech, R. & Hampshire, A. Dynamic network coding of working-memory domains and working-memory processes. Nat. Commun. 10, 936 (2019).30804436 10.1038/s41467-019-08840-8
16. Pinotsis DA Buschman TJ Miller EK Working Memory Load Modulates Neuronal Coupling Cereb. Cortex 2019 29 1670 1681 10.1093/cercor/bhy065 29608671
Pinotsis, D. A., Buschman, T. J. & Miller, E. K. Working Memory Load Modulates Neuronal Coupling. Cereb. Cortex 29, 1670–1681 (2019).29608671 10.1093/cercor/bhy065
17. Eryilmaz H Working memory load-dependent changes in cortical network connectivity estimated by machine learning NeuroImage 2020 217 116895 10.1016/j.neuroimage.2020.116895 32360929
Eryilmaz, H. et al. Working memory load-dependent changes in cortical network connectivity estimated by machine learning. NeuroImage 217, 116895 (2020).32360929 10.1016/j.neuroimage.2020.116895
18. Bastos AM Loonis R Kornblith S Lundqvist M Miller EK Laminar recordings in frontal cortex suggest distinct layers for maintenance and control of working memory Proc. Natl. Acad. Sci. USA 2018 115 1117 1122 10.1073/pnas.1710323115 29339471
Bastos, A. M., Loonis, R., Kornblith, S., Lundqvist, M. & Miller, E. K. Laminar recordings in frontal cortex suggest distinct layers for maintenance and control of working memory. Proc. Natl. Acad. Sci. USA 115, 1117–1122 (2018).29339471 10.1073/pnas.1710323115
19. Bastos AM Lundqvist M Waite AS Kopell N Miller EK Layer and rhythm specificity for predictive routing Proc. Natl. Acad. Sci. 2020 117 31459 31469 10.1073/pnas.2014868117 33229572
Bastos, A. M., Lundqvist, M., Waite, A. S., Kopell, N. & Miller, E. K. Layer and rhythm specificity for predictive routing. Proc. Natl. Acad. Sci. 117, 31459–31469 (2020).33229572 10.1073/pnas.2014868117
20. Goldman-Rakic PS Cellular basis of working memory Neuron 1995 14 477 485 10.1016/0896-6273(95)90304-6 7695894
Goldman-Rakic, P. S. Cellular basis of working memory. Neuron 14, 477–485 (1995).7695894 10.1016/0896-6273(95)90304-6
21. Arnsten AFT Wang MJ Paspalas CD Neuromodulation of thought: flexibilities and vulnerabilities in prefrontal cortical network synapses Neuron 2012 76 223 239 10.1016/j.neuron.2012.08.038 23040817
Arnsten, A. F. T., Wang, M. J. & Paspalas, C. D. Neuromodulation of thought: flexibilities and vulnerabilities in prefrontal cortical network synapses. Neuron 76, 223–239 (2012).23040817 10.1016/j.neuron.2012.08.038
22. Wang M Vijayraghavan S Goldman-Rakic PS Selective D2 receptor actions on the functional circuitry of working memory Science 2004 303 853 856 10.1126/science.1091162 14764884
Wang, M., Vijayraghavan, S. & Goldman-Rakic, P. S. Selective D2 receptor actions on the functional circuitry of working memory. Science 303, 853–856 (2004).14764884 10.1126/science.1091162
23. Lawrence SJD Laminar organization of working memory signals in human visual cortex Curr. Biol. 2018 28 3435 3440.e4 10.1016/j.cub.2018.08.043 30344121
Lawrence, S. J. D. et al. Laminar organization of working memory signals in human visual cortex. Curr. Biol. 28, 3435–3440.e4 (2018).30344121 10.1016/j.cub.2018.08.043
24. Lawrence SJD Norris DG de Lange FP Dissociable laminar profiles of concurrent bottom-up and top-down modulation in the human visual cortex eLife 2019 8 e44422 10.7554/eLife.44422 31063127
Lawrence, S. J. D., Norris, D. G. & de Lange, F. P. Dissociable laminar profiles of concurrent bottom-up and top-down modulation in the human visual cortex. eLife 8, e44422 (2019).31063127 10.7554/eLife.44422
25. Muckli L Contextual Feedback to Superficial Layers of V1 Curr. Biol. 2015 25 2690 2695 10.1016/j.cub.2015.08.057 26441356
Muckli, L. et al. Contextual Feedback to Superficial Layers of V1. Curr. Biol. 25, 2690–2695 (2015).26441356 10.1016/j.cub.2015.08.057
26. Kok P Bains LJ van Mourik T Norris DG de Lange FP Selective Activation of the Deep Layers of the Human Primary Visual Cortex by Top-Down Feedback Curr. Biol. 2016 26 371 376 10.1016/j.cub.2015.12.038 26832438
Kok, P., Bains, L. J., van Mourik, T., Norris, D. G. & de Lange, F. P. Selective Activation of the Deep Layers of the Human Primary Visual Cortex by Top-Down Feedback. Curr. Biol. 26, 371–376 (2016).26832438 10.1016/j.cub.2015.12.038
27. De Martino F Frequency preference and attention effects across cortical depths in the human primary auditory cortex Proc. Natl. Acad. Sci. 2015 112 16036 16041 10.1073/pnas.1507552112 26668397
De Martino, F. et al. Frequency preference and attention effects across cortical depths in the human primary auditory cortex. Proc. Natl. Acad. Sci. 112, 16036–16041 (2015).26668397 10.1073/pnas.1507552112
28. Aitken F Prior expectations evoke stimulus-specific activity in the deep layers of the primary visual cortex PLOS Biol 2020 18 e3001023 10.1371/journal.pbio.3001023 33284791
Aitken, F. et al. Prior expectations evoke stimulus-specific activity in the deep layers of the primary visual cortex. PLOS Biol 18, e3001023 (2020).33284791 10.1371/journal.pbio.3001023
29. Finn ES Huber L Jangraw DC Molfese PJ Bandettini PA Layer-dependent activity in human prefrontal cortex during working memory Nat. Neurosci. 2019 22 1687 1695 10.1038/s41593-019-0487-z 31551596
Finn, E. S., Huber, L., Jangraw, D. C., Molfese, P. J. & Bandettini, P. A. Layer-dependent activity in human prefrontal cortex during working memory. Nat. Neurosci. 22, 1687–1695 (2019).31551596 10.1038/s41593-019-0487-z
30. Spaak E Watanabe K Funahashi S Stokes MG Stable and dynamic coding for working memory in primate prefrontal cortex J. Neurosci 2017 37 6503 6516 10.1523/JNEUROSCI.3364-16.2017 28559375
Spaak, E., Watanabe, K., Funahashi, S. & Stokes, M. G. Stable and dynamic coding for working memory in primate prefrontal cortex. J. Neurosci 37, 6503–6516 (2017).28559375 10.1523/JNEUROSCI.3364-16.2017
31. Stokes MG Dynamic Coding for Cognitive Control in Prefrontal Cortex Neuron 2013 78 364 375 10.1016/j.neuron.2013.01.039 23562541
Stokes, M. G. et al. Dynamic Coding for Cognitive Control in Prefrontal Cortex. Neuron 78, 364–375 (2013).23562541 10.1016/j.neuron.2013.01.039
32. Cavanagh SE Towers JP Wallis JD Hunt LT Kennerley SW Reconciling persistent and dynamic hypotheses of working memory coding in prefrontal cortex Nat. Commun. 2018 9 3498 10.1038/s41467-018-05873-3 30158519
Cavanagh, S. E., Towers, J. P., Wallis, J. D., Hunt, L. T. & Kennerley, S. W. Reconciling persistent and dynamic hypotheses of working memory coding in prefrontal cortex. Nat. Commun. 9, 3498 (2018).30158519 10.1038/s41467-018-05873-3
33. Miller EK Lundqvist M Bastos AM Working Memory 2.0 Neuron 2018 100 463 475 10.1016/j.neuron.2018.09.023 30359609
Miller, E. K., Lundqvist, M. & Bastos, A. M. Working Memory 2.0. Neuron 100, 463–475 (2018).30359609 10.1016/j.neuron.2018.09.023
34. Meyers EM Freedman DJ Kreiman G Miller EK Poggio T Dynamic Population Coding of Category Information in Inferior Temporal and Prefrontal Cortex J. Neurophysiol. 2008 100 1407 1419 10.1152/jn.90248.2008 18562555
Meyers, E. M., Freedman, D. J., Kreiman, G., Miller, E. K. & Poggio, T. Dynamic Population Coding of Category Information in Inferior Temporal and Prefrontal Cortex. J. Neurophysiol. 100, 1407–1419 (2008).18562555 10.1152/jn.90248.2008
35. Iamshchinina, P. et al. Benchmarking GE-BOLD, SE-BOLD, and SS-SI-VASO sequences for depth-dependent separation of feedforward and feedback signals in high-field MRI (Neuroscience) (2021).
36. Iamshchinina P Perceived and mentally rotated contents are differentially represented in cortical depth of V1 Commun. Biol. 2021 4 1069 10.1038/s42003-021-02582-4 34521987
Iamshchinina, P. et al. Perceived and mentally rotated contents are differentially represented in cortical depth of V1. Commun. Biol. 4, 1069 (2021).34521987 10.1038/s42003-021-02582-4
37. Adam KCS Vogel EK Awh E Multivariate analysis reveals a generalizable human electrophysiological signature of working memory load Psychophysiology 2020 57 e13691 10.1111/psyp.13691 33040349
Adam, K. C. S., Vogel, E. K. & Awh, E. Multivariate analysis reveals a generalizable human electrophysiological signature of working memory load. Psychophysiology 57, e13691 (2020).33040349 10.1111/psyp.13691
38. Thyer, W., et al Storage in Visual Working Memory Recruits a Content-Independent Pointer System. Psychol. Sci., 095679762210909. (2022).
39. Majerus S Cross-Modal Decoding of Neural Patterns Associated with Working Memory: Evidence for Attention-Based Accounts of Working Memory Cereb. Cortex 2016 26 166 179 10.1093/cercor/bhu189 25146374
Majerus, S. et al. Cross-Modal Decoding of Neural Patterns Associated with Working Memory: Evidence for Attention-Based Accounts of Working Memory. Cereb. Cortex 26, 166–179 (2016).25146374 10.1093/cercor/bhu189
40. Majerus S Péters F Bouffier M Cowan N Phillips C The Dorsal Attention Network Reflects Both Encoding Load and Top–down Control during Working Memory J. Cogn. Neurosci. 2018 30 144 159 10.1162/jocn_a_01195 28984526
Majerus, S., Péters, F., Bouffier, M., Cowan, N. & Phillips, C. The Dorsal Attention Network Reflects Both Encoding Load and Top–down Control during Working Memory. J. Cogn. Neurosci. 30, 144–159 (2018).28984526 10.1162/jocn_a_01195
41. 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 32119658
Wolff, M. J., Jochim, J., Akyürek, E. G., Buschman, T. J. & Stokes, M. G. Drifting codes within a stable coding scheme for working memory. PLOS Biol 18, e3000625 (2020).32119658 10.1371/journal.pbio.3000625
42. Li H-H Curtis CE Neural population dynamics of human working memory Curr. Biol. 2023 33 3775 3784.e4 10.1016/j.cub.2023.07.067 37595590
Li, H.-H. & Curtis, C. E. Neural population dynamics of human working memory. Curr. Biol. 33, 3775–3784.e4 (2023).37595590 10.1016/j.cub.2023.07.067
43. Christophel TB Iamshchinina P Yan C Allefeld C Haynes J-D Cortical specialization for attended versus unattended working memory Nat. Neurosci. 2018 21 494 496 10.1038/s41593-018-0094-4 29507410
Christophel, T. B., Iamshchinina, P., Yan, C., Allefeld, C. & Haynes, J.-D. Cortical specialization for attended versus unattended working memory. Nat. Neurosci. 21, 494–496 (2018).29507410 10.1038/s41593-018-0094-4
44. Rademaker RL Chunharas C Serences JT Coexisting representations of sensory and mnemonic information in human visual cortex Nat. Neurosci. 2019 22 1336 1344 10.1038/s41593-019-0428-x 31263205
Rademaker, R. L., Chunharas, C. & Serences, J. T. Coexisting representations of sensory and mnemonic information in human visual cortex. Nat. Neurosci. 22, 1336–1344 (2019).31263205 10.1038/s41593-019-0428-x
45. D’Esposito M Postle BR Rypma B Prefrontal cortical contributions to working memory: evidence from event-related fMRI studies Exp. Brain Res. 2000 133 3 11 10.1007/s002210000395 10933205
D’Esposito, M., Postle, B. R. & Rypma, B. Prefrontal cortical contributions to working memory: evidence from event-related fMRI studies. Exp. Brain Res. 133, 3–11 (2000).10933205 10.1007/s002210000395
46. Rypma B Berger JS D’Esposito M The Influence of Working-Memory Demand and Subject Performance on Prefrontal Cortical Activity J. Cogn. Neurosci 2002 14 721 731 10.1162/08989290260138627 12167257
Rypma, B., Berger, J. S. & D’Esposito, M. The Influence of Working-Memory Demand and Subject Performance on Prefrontal Cortical Activity. J. Cogn. Neurosci 14, 721–731 (2002).12167257 10.1162/08989290260138627
47. Lorenz R The Automatic Neuroscientist: A framework for optimizing experimental design with closed-loop real-time fMRI NeuroImage 2016 129 320 334 10.1016/j.neuroimage.2016.01.032 26804778
Lorenz, R. et al. The Automatic Neuroscientist: A framework for optimizing experimental design with closed-loop real-time fMRI. NeuroImage 129, 320–334 (2016).26804778 10.1016/j.neuroimage.2016.01.032
48. Lorenz R Dissociating frontoparietal brain networks with neuroadaptive Bayesian optimization Nat. Commun. 2018 9 1227 10.1038/s41467-018-03657-3 29581425
Lorenz, R. et al. Dissociating frontoparietal brain networks with neuroadaptive Bayesian optimization. Nat. Commun. 9, 1227 (2018).29581425 10.1038/s41467-018-03657-3
49. Lorenz R A Bayesian optimization approach for rapidly mapping residual network function in stroke Brain 2021 144 2120 2134 10.1093/brain/awab109 33725125
Lorenz, R. et al. A Bayesian optimization approach for rapidly mapping residual network function in stroke. Brain 144, 2120–2134 (2021).33725125 10.1093/brain/awab109
50. Miller EK Cohen JD An Integrative Theory of Prefrontal Cortex Function Annu. Rev. Neurosci. 2001 24 167 202 10.1146/annurev.neuro.24.1.167 11283309
Miller, E. K. & Cohen, J. D. An Integrative Theory of Prefrontal Cortex Function. Annu. Rev. Neurosci. 24, 167–202 (2001).11283309 10.1146/annurev.neuro.24.1.167
51. O’Reilly RC Frank MJ Making working memory work: a computational model of learning in the prefrontal cortex and basal ganglia Neural Comput 2006 18 283 328 10.1162/089976606775093909 16378516
O’Reilly, R. C. & Frank, M. J. Making working memory work: a computational model of learning in the prefrontal cortex and basal ganglia. Neural Comput 18, 283–328 (2006).16378516 10.1162/089976606775093909
52. Friedman NP Robbins TW The role of prefrontal cortex in cognitive control and executive function Neuropsychopharmacology 2022 47 72 89 10.1038/s41386-021-01132-0 34408280
Friedman, N. P. & Robbins, T. W. The role of prefrontal cortex in cognitive control and executive function. Neuropsychopharmacology 47, 72–89 (2022).34408280 10.1038/s41386-021-01132-0
53. Hussar C Pasternak T Trial-to-trial variability of the prefrontal neurons reveals the nature of their engagement in a motion discrimination task Proc. Natl. Acad. Sci. USA 2010 107 21842 21847 10.1073/pnas.1009956107 21098286
Hussar, C. & Pasternak, T. Trial-to-trial variability of the prefrontal neurons reveals the nature of their engagement in a motion discrimination task. Proc. Natl. Acad. Sci. USA 107, 21842–21847 (2010).21098286 10.1073/pnas.1009956107
54. Buschman TJ Miller EK Working memory is complex and dynamic, like your thoughts J. Cogn. Neurosci 2022 35 17 23 10.1162/jocn_a_01940 36322832
Buschman, T. J. & Miller, E. K. Working memory is complex and dynamic, like your thoughts. J. Cogn. Neurosci 35, 17–23 (2022).36322832 10.1162/jocn_a_01940
55. Curtis CE Sprague TC Persistent activity during working memory from front to back Front. Neural Circuits 2021 15 696060 10.3389/fncir.2021.696060 34366794
Curtis, C. E. & Sprague, T. C. Persistent activity during working memory from front to back. Front. Neural Circuits 15, 696060 (2021).34366794 10.3389/fncir.2021.696060
56. Sreenivasan KK D’Esposito M The what, where and how of delay activity Nat. Rev. Neurosci. 2019 20 466 481 10.1038/s41583-019-0176-7 31086326
Sreenivasan, K. K. & D’Esposito, M. The what, where and how of delay activity. Nat. Rev. Neurosci. 20, 466–481 (2019).31086326 10.1038/s41583-019-0176-7
57. Constantinidis C Persistent spiking activity underlies working memory J. Neurosci 2018 38 7020 7028 10.1523/JNEUROSCI.2486-17.2018 30089641
Constantinidis, C. et al. Persistent spiking activity underlies working memory. J. Neurosci 38, 7020–7028 (2018).30089641 10.1523/JNEUROSCI.2486-17.2018
58. Stokes MG ‘Activity-silent’ working memory in prefrontal cortex: a dynamic coding framework Trends Cogn. Sci. 2015 19 394 405 10.1016/j.tics.2015.05.004 26051384
Stokes, M. G. ‘Activity-silent’ working memory in prefrontal cortex: a dynamic coding framework. Trends Cogn. Sci. 19, 394–405 (2015).26051384 10.1016/j.tics.2015.05.004
59. Huber L Cortical lamina-dependent blood volume changes in human brain at 7 T NeuroImage 2015 107 23 33 10.1016/j.neuroimage.2014.11.046 25479018
Huber, L. et al. Cortical lamina-dependent blood volume changes in human brain at 7 T. NeuroImage 107, 23–33 (2015).25479018 10.1016/j.neuroimage.2014.11.046
60. Huber L Techniques for blood volume fMRI with VASO: From low-resolution mapping towards sub-millimeter layer-dependent applications NeuroImage 2018 164 131 143 10.1016/j.neuroimage.2016.11.039 27867088
Huber, L. et al. Techniques for blood volume fMRI with VASO: From low-resolution mapping towards sub-millimeter layer-dependent applications. NeuroImage 164, 131–143 (2018).27867088 10.1016/j.neuroimage.2016.11.039
61. Polimeni JR Fischl B Greve DN Wald LL Laminar analysis of 7T BOLD using an imposed spatial activation pattern in human V1 NeuroImage 2010 52 1334 1346 10.1016/j.neuroimage.2010.05.005 20460157
Polimeni, J. R., Fischl, B., Greve, D. N. & Wald, L. L. Laminar analysis of 7T BOLD using an imposed spatial activation pattern in human V1. NeuroImage 52, 1334–1346 (2010).20460157 10.1016/j.neuroimage.2010.05.005
62. Chaimow D Yacoub E Uğurbil K Shmuel A Spatial specificity of the functional MRI blood oxygenation response relative to neuronal activity NeuroImage 2018 164 32 47 10.1016/j.neuroimage.2017.08.077 28882632
Chaimow, D., Yacoub, E., Uğurbil, K. & Shmuel, A. Spatial specificity of the functional MRI blood oxygenation response relative to neuronal activity. NeuroImage 164, 32–47 (2018).28882632 10.1016/j.neuroimage.2017.08.077
63. Markuerkiaga I Marques JP Gallagher TE Norris DG Estimation of laminar BOLD activation profiles using deconvolution with a physiological point spread function J. Neurosci. Methods 2021 353 109095 10.1016/j.jneumeth.2021.109095 33549635
Markuerkiaga, I., Marques, J. P., Gallagher, T. E. & Norris, D. G. Estimation of laminar BOLD activation profiles using deconvolution with a physiological point spread function. J. Neurosci. Methods 353, 109095 (2021).33549635 10.1016/j.jneumeth.2021.109095
64. Heinzle J Koopmans PJ den Ouden HEM Raman S Stephan KE A hemodynamic model for layered BOLD signals NeuroImage 2016 125 556 570 10.1016/j.neuroimage.2015.10.025 26484827
Heinzle, J., Koopmans, P. J., den Ouden, H. E. M., Raman, S. & Stephan, K. E. A hemodynamic model for layered BOLD signals. NeuroImage 125, 556–570 (2016).26484827 10.1016/j.neuroimage.2015.10.025
65. Koopmans PJ Barth M Norris DG Layer-specific BOLD activation in human V1 Hum. Brain Mapp. 2010 31 1297 1304 10.1002/hbm.20936 20082333
Koopmans, P. J., Barth, M. & Norris, D. G. Layer-specific BOLD activation in human V1. Hum. Brain Mapp. 31, 1297–1304 (2010).20082333 10.1002/hbm.20936
66. de Hollander G van der Zwaag W Qian C Zhang P Knapen T Ultra-high field fMRI reveals origins of feedforward and feedback activity within laminae of human ocular dominance columns NeuroImage 2021 228 117683 10.1016/j.neuroimage.2020.117683 33385565
de Hollander, G., van der Zwaag, W., Qian, C., Zhang, P. & Knapen, T. Ultra-high field fMRI reveals origins of feedforward and feedback activity within laminae of human ocular dominance columns. NeuroImage 228, 117683 (2021).33385565 10.1016/j.neuroimage.2020.117683
67. Havlicek M Uludağ K A dynamical model of the laminar BOLD response NeuroImage 2020 204 116209 10.1016/j.neuroimage.2019.116209 31546051
Havlicek, M. & Uludağ, K. A dynamical model of the laminar BOLD response. NeuroImage 204, 116209 (2020).31546051 10.1016/j.neuroimage.2019.116209
68. Bergmann J., Morgan A. T., & Muckli L. Two distinct feedback codes in V1 for ‘real’ and ‘imaginary’ internal experiences (Neuroscience). (2019).
69. Sharoh D Laminar specific fMRI reveals directed interactions in distributed networks during language processing Proc. Natl. Acad. Sci. USA 2019 116 21185 21190 10.1073/pnas.1907858116 31570628
Sharoh, D. et al. Laminar specific fMRI reveals directed interactions in distributed networks during language processing. Proc. Natl. Acad. Sci. USA 116, 21185–21190 (2019).31570628 10.1073/pnas.1907858116
70. Huang, P. et al. Correcting for Superficial Bias in 7T Gradient Echo fMRI. Front. Neurosci. 15 (2021).
71. Righi G Peissig JJ Tarr MJ Recognizing disguised faces Vis. Cogn. 2012 20 143 169 10.1080/13506285.2012.654624
Righi, G., Peissig, J. J. & Tarr, M. J. Recognizing disguised faces. Vis. Cogn. 20, 143–169 (2012).10.1080/13506285.2012.654624
72. Fischl B FreeSurfer NeuroImage 2012 62 774 781 10.1016/j.neuroimage.2012.01.021 22248573
Fischl, B. FreeSurfer. NeuroImage 62, 774–781 (2012).22248573 10.1016/j.neuroimage.2012.01.021
73. Sriranga Kashyap (2021). srikash/3dMPRAGEise: ondu. Version 1.0 (Zenodo). 10.5281/ZENODO.4626825.
74. Gaser, C., Dahnke, R., Thompson, P. M., Kurth, F., Luders, E., and Alzheimer’s Disease Neuroimaging Initiative. CAT – A Computational Anatomy Toolbox for the Analysis of Structural MRI Data (Neuroscience). 10.1101/2022.06.11.495736. (2022)
75. Dickie EW Ciftify: A framework for surface-based analysis of legacy MR acquisitions NeuroImage 2019 197 818 826 10.1016/j.neuroimage.2019.04.078 31091476
Dickie, E. W. et al. Ciftify: A framework for surface-based analysis of legacy MR acquisitions. NeuroImage 197, 818–826 (2019).31091476 10.1016/j.neuroimage.2019.04.078
76. Glasser MF A multi-modal parcellation of human cerebral cortex Nature 2016 536 171 178 10.1038/nature18933 27437579
Glasser, M. F. et al. A multi-modal parcellation of human cerebral cortex. Nature 536, 171–178 (2016).27437579 10.1038/nature18933
77. Robinson EC MSM: A new flexible framework for Multimodal Surface Matching NeuroImage 2014 100 414 426 10.1016/j.neuroimage.2014.05.069 24939340
Robinson, E. C. et al. MSM: A new flexible framework for Multimodal Surface Matching. NeuroImage 100, 414–426 (2014).24939340 10.1016/j.neuroimage.2014.05.069
78. Avants BB A reproducible evaluation of ANTs similarity metric performance in brain image registration NeuroImage 2011 54 2033 2044 10.1016/j.neuroimage.2010.09.025 20851191
Avants, B. B. et al. A reproducible evaluation of ANTs similarity metric performance in brain image registration. NeuroImage 54, 2033–2044 (2011).20851191 10.1016/j.neuroimage.2010.09.025
79. Haenelt, D., Chaimow, D., Nasr, S., Weiskopf, N., & Trampel, R. Decoding of columnar-level organization across cortical depth using BOLD- and CBV-fMRI at 7 T (Neuroscience) (2023).
80. Cox RW AFNI: Software for Analysis and Visualization of Functional Magnetic Resonance Neuroimages Comput. Biomed. Res. 1996 29 162 173 10.1006/cbmr.1996.0014 8812068
Cox, R. W. AFNI: Software for Analysis and Visualization of Functional Magnetic Resonance Neuroimages. Comput. Biomed. Res. 29, 162–173 (1996).8812068 10.1006/cbmr.1996.0014
81. Huber, L. et al. LayNii: A software suite for layer-fMRI. NeuroImage 237, 118091 (2021).
82. Coalson, T. S., Van Essen, D. C., and Glasser, M. F. The impact of traditional neuroimaging methods on the spatial localization of cortical areas. Proc. Natl. Acad. Sci. USA 115. 10.1073/pnas.1801582115. (2018)
83. Ji JL Mapping the human brain’s cortical-subcortical functional network organization NeuroImage 2019 185 35 57 10.1016/j.neuroimage.2018.10.006 30291974
Ji, J. L. et al. Mapping the human brain’s cortical-subcortical functional network organization. NeuroImage 185, 35–57 (2019).30291974 10.1016/j.neuroimage.2018.10.006
84. Chang C-C Lin C-J LIBSVM: A library for support vector machines. ACM Trans Intell. Syst. Technol. 2011 2 1 27 10.1145/1961189.1961199
Chang, C.-C. & Lin, C.-J. LIBSVM: A library for support vector machines. ACM Trans. Intell. Syst. Technol. 2, 1–27 (2011).10.1145/1961189.1961199
85. Maris E Oostenveld R Nonparametric statistical testing of EEG- and MEG-data J. Neurosci. Methods 2007 164 177 190 10.1016/j.jneumeth.2007.03.024 17517438
Maris, E. & Oostenveld, R. Nonparametric statistical testing of EEG- and MEG-data. J. Neurosci. Methods 164, 177–190 (2007).17517438 10.1016/j.jneumeth.2007.03.024
