Skip to main content

Main menu

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

User menu

  • Log out
  • Log in
  • My Cart

Search

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

Advanced Search

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

Stimulus-Specific Visual Working Memory Representations in Human Cerebellar Lobule VIIb/VIIIa

James A. Brissenden, Sean M. Tobyne, Mark A. Halko and David C. Somers
Journal of Neuroscience 3 February 2021, 41 (5) 1033-1045; https://doi.org/10.1523/JNEUROSCI.1253-20.2020
James A. Brissenden
1Department of Psychology, University of Michigan, Ann Arbor, Michigan 48109
2Department of Psychological and Brain Sciences, Boston University, Boston, Massachusetts 02215
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for James A. Brissenden
Sean M. Tobyne
3Graduate Program for Neuroscience, Boston University, Boston, Massachusetts 02215
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Mark A. Halko
4Psychotic Disorders Division, McLean Hospital, Belmont, MA, 02478; Harvard Medical School, Boston, MA, 02115
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Mark A. Halko
David C. Somers
2Department of Psychological and Brain Sciences, Boston University, Boston, Massachusetts 02215
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for David C. Somers
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF
Loading

Abstract

fMRI research has revealed that cerebellar lobule VIIb/VIIIa exhibits load-dependent activity that increases with the number of items held in visual working memory (VWM). However, it remains unclear whether these cerebellar responses reflect processes specific to VWM or more general visual attentional mechanisms. To investigate this question, we examined whether cerebellar activity during the delay period of a VWM task is selective for stimuli held in working memory. A sample of male and female human subjects performed a VWM continuous report task in which they were retroactively cued to remember the direction of motion of moving dot stimuli. Cerebellar lobule VIIb/VIIIa delay-period activation accurately decoded the direction of the remembered stimulus, as did frontal and parietal regions of the dorsal attention network. Arguing against a motor explanation, no other cerebellar area exhibited stimulus specificity, including the oculomotor vermis, a key area associated with eye movement control. Finer-scale analysis revealed that the medial portion of lobule VIIb and to a lesser degree the lateral most portion of lobules VIIb and VIIIa, which exhibit robust resting state connectivity with frontal and parietal regions of the dorsal attention network, encoded the identity of the remembered stimulus, while intermediate portions of lobule VIIb/VIIIa did not. These findings of stimulus-specific coding of VWM within lobule VIIb/VIIIa indicate for the first time that the distributed network responsible for the encoding and maintenance of mnemonic representations extends to the cerebellum.

SIGNIFICANCE STATEMENT There is considerable debate concerning where in the brain the contents of visual working memory (VWM) are stored. To date, this literature has primarily focused on the role of regions located within cerebral cortex. There is growing evidence for cerebellar involvement in higher-order cognitive functions including working memory. While the cerebellum has been previously shown to be recruited by VWM paradigms, it is unclear whether any portion of cerebellum actively encodes and maintains mnemonic representations. The present study demonstrates that cerebellar lobule VIIb/VIIIa activity patterns are selective for remembered stimuli and that this selectivity persists in the absence of perceptual input. These findings provide novel evidence for the participation of cerebellar structures in the persistent storage of visual information.

  • cerebellum
  • cortico-cerebellar networks
  • delay activity
  • fMRI
  • stimulus selectivity
  • visual working memory

Introduction

Visual working memory (VWM) can be broadly defined as the temporary maintenance of visual information independent of the constant influx of perceptual input. Investigating the neural mechanisms of VWM is of considerable interest to the neuroscience community as performance on VWM tasks has been shown to be highly correlated with broader measures of cognitive ability such as fluid intelligence, reading comprehension, and scholastic aptitude (Engle et al., 1999; Engle, 2002; Fukuda et al., 2010). The question of where working memory contents are stored in the brain is the subject of ongoing debate. Based on electrophysiological recordings in non-human primates and neuroimaging in humans, it has long been asserted that parietal and prefrontal cortices support working memory maintenance (Fuster and Alexander, 1971; Funahashi et al., 1989; Chafee and Goldman-Rakic, 1998; Courtney et al., 1998; Mendoza-Halliday et al., 2014). On the other hand, the sensory recruitment hypothesis posits that working memory storage is mediated by the same areas involved in the initial sensory processing of stimuli and that parietal and prefrontal areas instead serve as a source of top-down biasing signals (Pasternak and Greenlee, 2005; D'Esposito and Postle, 2015). Recently, it has been suggested that working memory contents are distributed across number of cerebral cortical areas including both sensory and prefrontal regions (Ester et al., 2015; Lee and Baker, 2016; Serences, 2016; Christophel et al., 2017). Despite the progress that has been made in characterizing the neural substrates of working memory encoding and storage, the involvement of subcortical structures has received little consideration. A number of studies have reported that portions of the cerebellum are actively involved in attention and working memory task performance (Allen et al., 1997; Desmond et al., 1997; Chen and Desmond, 2005a; Kirschen et al., 2005; Hayter et al., 2007; Baier et al., 2010; Stoodley et al., 2012; Striemer et al., 2015; Brissenden et al., 2016, 2018; Guell et al., 2018a; King et al., 2019). However, as these studies either did not include a VWM maintenance component or employed a block design in which sample, delay, and probe responses could not be distinguished, it is unclear whether the cerebellar contribution to these tasks reflects working memory storage processes or rather a more generalized attentional control mechanism. The cerebellum has been proposed to play a predictive role in working memory (Sheu et al., 2019), mirroring cerebellar models of motor control (Wolpert et al., 1998). In order to support a predictive function in WM, the cerebellum would need to represent the item-specific contents of working memory. In the current study, we sought to directly investigate whether specific portions of the cerebellum contribute to the encoding and maintenance of mnemonic representations.

Early research aimed at determining where working memory contents are stored in the brain was primarily informed by either lesion-induced deficits in working memory task performance or the observation of sustained, elevated activity over extended delays (Funahashi et al., 1989; Courtney et al., 1998; Postle et al., 2000; Wager and Smith, 2003; Sreenivasan and D'Esposito, 2019). However, working memory deficits following lesioning or elevated delay period activity could simply indicate that an area influences storage in other areas and does not imply that an area contains a representation of working memory contents (D'Esposito and Postle, 2015; Lee and Baker, 2016; Serences, 2016; Christophel et al., 2017). More recently, researchers have employed multivariate methods to determine whether an area exhibits activity patterns that are selective for the specific stimulus stored in working memory (Serences et al., 2009; Harrison and Tong, 2009; Christophel et al., 2012, 2015, 2018; Riggall and Postle, 2012; Lee et al., 2013; Pratte and Tong, 2014; Ester et al., 2015). As a result, content specificity is now argued to be the critical marker of information storage in VWM (Christophel et al., 2017). Here, we used an event-related delayed recall paradigm and a multivariate encoding model to investigate whether this distributed network extends to the cerebellum.

Materials and Methods

Participants

A total of 16 healthy adult volunteers (nine female) participated in this study. The Institutional Board of Boston University approved the study. All subjects gave written informed consent and were compensated to participate in this study. Subjects were recruited from Boston University and the Greater Boston area. Our sample included six laboratory members (including one author), each of which had extensive experience as an experimental subject in neuroimaging studies. The remaining participants were not lab members but had prior experience performing psychophysical tasks in an fMRI scanner. All subjects possessed normal or corrected-to-normal vision. All subjects completed a ∼30-min behavioral session before being scanned to familiarize them with the task. One subject's behavioral accuracy was found to not significantly differ from chance performance and was excluded from further analysis, leaving us with 15 subjects (nine female).

Magnetic resonance image acquisition

Data were acquired from a 3 Tesla Siemens Prisma scanner located at the Cognitive Neuroimaging Center at Boston University in Boston, Massachusetts using a 64-channel head coil. High-resolution T1-weighted multiecho MPRAGE (sagittal; TR = 2780 ms; TE = 1.32, 3.19, 5.11, 7.03 ms; FA = 7°; 0.8-mm isotropic voxels; 224 slices; FOV = 256 × 256 × 180 mm; in-plane GRAPPA acceleration 2) and T2-weighted (sagittal; TR = 3200 ms; TE = 564 ms; 0.8-mm isotropic voxels; 224 slices; FOV = 256 × 256 × 180 mm; in-plane GRAPPA acceleration 2) structural images were acquired. Functional data were acquired using a multiband gradient-echo echoplanar pulse sequence (Feinberg et al., 2010; Moeller et al., 2010; Xu et al., 2013) with the following acquisition parameters: TR = 2000 ms; TE = 35 ms; FA = 80°; 2.2-mm isotropic voxels; FOV = 207 × 207 × 152 mm; multiband factor (SMS) = 3. A total of 69 slices were acquired with 0% skip, fully covering the cerebral cortex and cerebellum. Image reconstruction was performed using the split-slice GRAPPA method, which has been shown to dramatically reduce the incidence of false positive activation because of signal leakage (Cauley et al., 2014; Todd et al., 2016). Spin echo field maps were also acquired with opposite phase encoding directions (anterior-to-posterior; posterior-to-anterior) and matching parameters to the gradient-echo EPI fMRI timeseries for subsequent EPI distortion correction.

Visual stimuli and experimental paradigm

Stimuli were generated and presented using Python with the PsychoPy software package (Peirce, 2007, 2008) and were projected with a PROPixx digital light processing LED projector (VPixx Technologies Inc.) onto a screen within the scanner bore. Each trial consisted of a trial start cue period (1 s), sample period (1 s), retro-cue and delay period (10 s), and probe period (5 s; Fig. 1A). The beginning of each trial was denoted by the white fixation cross briefly changing to red for 500 ms. A 500-ms fixation interval separated the trial start cue and the sample period. During the sample period, subjects were presented with two circular random dot kinematograms (one in each hemifield) each spanning 8° in diameter. Dots within each patch were 0.05° in size and possessed a density of ∼2 dots per square degree. The dot motion patches were centered 6° to the left or right from fixation along the horizontal meridian. Dots moved with 100% coherence at 1.5°/s and had a limited lifetime of 10 frames (167 ms). The direction of coherent motion on a trial within each patch was drawn from a uniform distribution over 0−315° in 45° increments. A small ±1−10° angular jitter was added randomly to the motion direction on each trial to minimize verbal coding of stimuli; 250 ms following the offset of the sample, an arrow appeared pointing either to the left or right. This “retro-cue” indicated which of the two motion directions subjects should maintain over the subsequent 9.75-s delay period. The use of a retro-cue was a key element of the experimental design as it permitted us to distinguish between initial stimulus encoding and maintenance in VWM by comparing delay period activity between the remembered and non-remembered stimuli. On 20% of trials, an “X” was presented indicating that both items should be dropped. Following the delay period, subjects were presented with a probe display (5 s). The probe display consisted of a circle (8° diameter) centered at fixation and a line extending from the center of the display to a location on the circumference of the circle. Subjects were instructed to rotate the line so that it aligned with the remembered direction of motion. Subjects pressed one key to incrementally rotate the line 10° clockwise, another key to rotate the line 10° counterclockwise, and a third key to flip the line 180°. The initial direction of the line segment was selected randomly among the 36 possible offsets that are an integer multiple of 10° of the to-be-remembered motion direction. The randomization of the probe stimulus ensured that subjects were unable to prospectively plan a specific motor action (i.e., direction of rotation) during the delay period. On drop trials, subjects were instructed to make a random number of button-press responses. Trials were separated by a 7-s intertrial interval. Each fMRI run comprised 20 trials (16 store trials + 4 drop trials; 482 s total). Subjects completed eight runs during the session. In-scanner eye position tracking was planned, but because of hardware issues was not possible during this experiment. Stimulus direction and location (left or right visual hemifield) of the cued stimulus were fully crossed within a run. Additionally, the motion directions of cued and non-cued stimuli were fully crossed across all eight runs.

fMRI data preprocessing and regions of interest (ROIs)

Functional task data first underwent the Human Connectome Project (HCP)'s “minimal” preprocessing pipeline, which comprises gradient nonlinearity distortion correction, motion correction, EPI image distortion correction, co-registration with the subject's T1-weighted image, and nonlinear transformation to MNI152 space (Glasser et al., 2013). The transforms involved in each step of this pipeline were concatenated into a single nonlinear transformation and performed as a single resampling step to reduce interpolation related blurring (Glasser et al., 2013). Voxel timeseries within each run were then standardized and detrended with a third order polynomial. Multivariate analyses were performed on the average response of individual voxels across three TRs beginning 6, 8, and 10 s (TRs 4, 5, and 6) after the start of each trial, consistent with previous working memory decoding studies (Harrison and Tong, 2009; Ester et al., 2015). Additionally, we generated stimulus reconstructions across time by applying our decoding model to a sliding window of the average response across three consecutive TRs (0–4, 2–6 s, etc.)

Resting-state data underwent additional preprocessing using custom scripts in MATLAB. The following preprocessing steps were performed: linear interpolation across high-motion time-points (>0.5 mm FD; Power et al., 2012; Carp, 2013), application of a fourth-order Butterworth temporal bandpass filter to extract frequencies between 0.009 and 0.08 Hz, mean “grayordinate” signal regression (MGSR; Burgess et al., 2016) and censoring of high-motion time-points (Power et al., 2012).

ROIs were defined from publicly available atlases. Cortical ROIs corresponding with a combined IPS0-IPS3 ROI [left (L): 670 vertices; right (R): 606 vertices] and a combined V1-V3 ROI (L: 1008 vertices; R: 960 vertices) were defined in each hemisphere from a probabilistic retinotopy atlas (Wang et al., 2015). A putative VWM region in dorsolateral frontal cortex was defined from the 6a label from the Glasser multimodal parcellation in both hemispheres (L: 360 vertices; R: 339 vertices; Glasser et al., 2016); area 6a corresponds with the intersection of the superior frontal sulcus and the superior branch of the precentral sulcus, and exhibits robust working memory activation in the HCP VWM (2-back vs 0-back) task contrast (Glasser et al., 2016). This ROI will be referred to as superior frontal junction (SFJ) throughout the remainder of the manuscript. The cerebellar spatially unbiased infratentorial template (SUIT) atlas was used to define an ROI for each lobule within the cerebellum in each hemisphere (L: 133–1866 voxels; R: 140–1765 voxels; Diedrichsen et al., 2009). To exclude any potential signal spill-over (resulting from the nonlinear spatial transformation procedure) from ventral visual cortex into cerebellum (Buckner et al., 2011; van Es et al., 2019), we applied a dilated (6 mm) cortical ribbon mask to dorsal cerebellar lobule ROIs (lobules I–IV, V, and VI and Crus I). As our prior work indicates that an area spanning both lobule VIIb and lobule VIIIa is recruited by attention and working memory tasks along with fronto-parietal areas (Brissenden et al., 2016, 2018), we created a combined lobule VIIb/VIIIa ROI (L: 1660 voxels; R: 1731 voxels). An oculomotor vermis ROI (L: 99 voxels; R: 55 voxels) was defined using a publicly available cerebellar retinotopy atlas produced by an analysis of the HCP retinotopy dataset (van Es et al., 2019). We further created ROIs subdividing VIIb/VIIIa using resting-state functional connectivity. To do so, a resting-state correlation with each cortical ROI (V1-3, IPS0-3, SFJ) was computed for each voxel within lobule VIIb/VIIIa. As each cerebral cortical hemisphere is anatomically connected with the contralateral hemisphere of cerebellar cortex (Strick, 1985; Schmahmann and Pandya, 1997), left hemisphere cortical seeds were correlated with right hemisphere lobule VIIb/VIIIa voxels and vice versa. We then created ROIs for each subject corresponding with the top 25% of correlations within lobule VIIb/VIIIa for each cortical seed. A percentile threshold was used to account for potential SNR differences between subjects and to ensure that ROIs were of similar size across subjects. This analysis revealed a differential pattern of connectivity within lobule VIIb/VIIIa for each cortical seed. To quantify the location of each cortical area's connectivity within lobule VIIb/VIIIa across subjects, we summed 25th percentile ROIs across subjects to create a probabilistic subject count map for each cortical ROI. The top 25% of resting-state correlations for IPS0-3 and SFJ were found to be concentrated in the medial and lateral portions of lobule VIIb/VIIIa for the majority of subjects. In contrast, the top 25% of V1-3 correlations was more evenly distributed throughout lobule VIIb/VIIIa across subjects. We then subdivided lobule VIIb/VIIIa on the basis of this differential connectivity. Medial (L: 376 voxels; R: 322 voxels) and lateral (L: 105 voxels; R: 128 voxels) VIIb/VIIIa ROIs were defined by taking the maximum across the IPS0-3 and SFJ subject count maps and then thresholding at six or more (i.e., >33.3% of) subjects. A central lobule VIIb/VIIIa ROI (L: 161 voxels; R: 140 voxels) intermediate to the medial and lateral ROIs was defined by thresholding the V1-3 subject count map at six or more subjects and then excluding any voxel that already belonged to either the medial or lateral ROI.

Statistical analysis

Behavioral analysis

To assess each participant's behavioral accuracy, we computed the mean absolute value of the angular difference between reported and actual motion directions. We further fit several different mixture models to each participant's response error distribution using the MemToolBox package (Suchow et al., 2013). We fit the following models: standard mixture model (Zhang and Luck, 2008), swap model (Bayset al., 2009), variable precision model (van den Berg et al., 2012; Fougnie et al., 2012), and a swap + variable precision model (van den Berg et al., 2014). The standard mixture model fits participants' response errors with a mixture of two distributions, a uniform distribution and a von Mises distribution, reflecting (1) trials in which items were not stored in working memory and (2) the precision of items successfully represented in working memory [parameters: guess rate (g) and precision (SD)]. The swap model is an extension of the standard mixture model that includes an additional parameter representing the proportion of trials for which participants incorrectly reported the identity of the non-cued item [parameters: g, SD, swap rate (b)]. The variable precision model, on the other hand, assumes that the precision of a participant's reports is itself drawn from a higher-order variability distribution [parameters: g, mean precision (SDµ), precision variability (SDσ)]. The swap + variable precision combines swap and variable precision models and includes parameters associated with both non-target reports and precision variability (parameters: g, b, SDµ, and SDσ). Each model additionally included a bias parameter to allow the central tendency of the fitted distribution to differ from zero. To assess the relative goodness of fit between models, we computed the Akaike Information Criterion (AIC; Akaike, 1974) for each model: Embedded Image where L is the model log-likelihood and k is the number of free parameters. As AIC assumes infinite data and will tend to overfit for smaller datasets, we computed a corrected version of AIC (AICc; Hurvich and Tsai, 1989): Embedded Image where n represents the number of working memory trials used to fit the model. The variable precision model was determined to best fit our data (mean AICc difference ± SEM: 20.17 ± 5.29 vs standard mixture model, 22.25 ± 5.27 vs swap model, 2.43 ± 0.47 vs swap + variable precision model). Consequently, we report mean parameters for the best fitting variable precision model across subjects associated with guess rate (g), mean precision (SDµ), and precision variability (SDσ).

Generative model

We employed a generative model-based approach developed by van Bergen et al. (2015). This model assumes voxel responses can be modeled as the linear weighted sum of 8 motion direction selective neuronal populations or channels. Each channel was represented as a half-wave rectified (co)sinusoid raised to the fifth power: Embedded Image where θk indicates the preferred direction of the kth population or channel. Channel basis functions were maximally tuned at one of eight equally spaced angles (0°, 45°, 90°, 135°, 180°, 225°, 270°, or 315°).

Model fitting and assessment was performed using an eightfold leave-one-run-out cross-validation scheme, which iteratively partitioned the data into a training set (B1) and a test set (B2). The training set was expressed as an n × m matrix B1, where n is the number of trials in the training set and m is the number of voxels. To account for retinotopic preferences, the training and test sets were restricted to trials in which the stimulus appeared in the ROI's preferred hemifield (contralateral for cortical ROIs and ipsilateral for cerebellar ROIs; Brissenden et al., 2018; van Es et al., 2019). An n × k hypothetical channel output matrix C1 (where k is the number of channels) represented the idealized channel responses for each trial given the remembered motion direction. A k × m channel weight matrix W related the observed voxel responses to the hypothetical channel responses via a general linear model of the form: Embedded Image

The weight matrix was estimated via ordinary least squares as follows: Embedded Image

Then for each trial in the test set, we used Bayes rule to obtain a posterior probability distribution over motion direction indicating which stimulus value was most probable given the observed pattern of BOLD responses. Embedded Image where s denotes a particular stimulus value and b2 indicates a single trial/row of test set BOLD matrix (B2). The prior, p(s), was set to 1 for all angles as each stimulus value was presented equally often and we did not want to bias decoding toward particular motion directions (van Bergen et al., 2015). To avoid numerical underflow, we computed log-probability rather than probability. The conditional log-probability of a voxel activation pattern given a specific stimulus value was defined as: Embedded Image where c2 denotes a vector of hypothetical channel responses corresponding with the stimulus value s and Σ denotes the m × m voxel covariance matrix estimated from the training set. As computing the above conditional probability requires computing the inverse of the voxel covariance matrix, we performed a procedure described in Naselaris et al. (2009) if Σ was singular. First, we performed principal component analysis on the predicted responses (Embedded Image) for the training set. p components were chosen with the constraint that each component explained at least 5% of the total variance. Using the resulting m × p projection matrix, we then projected the predicted (Embedded Image) and actual (B1) responses from the training set onto the first p principal components. These dimensionality-reduced responses were then normalized to unit length: Embedded Image Embedded Image where P represents the computed projection matrix, Embedded Image and Embedded Image respectively denote a single row from the predicted (Embedded Image) and actual (B1) response matrix, and Embedded Image and Embedded Image, respectively, denote a single row from the dimensionality-reduced n × p predicted (Embedded Image) and actual (Embedded Image response matrix. The voxel covariance matrix was then computed as: Embedded Image where Σ* is a p × p matrix. During the testing phase, the same procedure (projection onto the first p components using the projection matrix obtained from the training set followed by normalization) was performed on the predicted and actual responses for the test set. Thus, the conditional log-probability of the dimensionality-reduced pattern of voxel responses given a specific stimulus value was re-expressed as: Embedded Image

To assess whether an area contains stimulus-specific information across trials, we then circularly shifted trial log-posterior distributions so that the center value corresponded with the presented stimulus on that trial. Probability densities, as normalized ratios, must be averaged using the geometric mean rather than the arithmetic mean (Fleming and Wallace, 1986; Nelson, 2017). This is equivalent to computing the arithmetic mean of log-probability and then exponentiating the result. Consequently, we averaged the shifted log-probability distributions across trials to produce a subject-level log-posterior probability distribution. We then averaged subject log-posterior distributions to form a group-level log-posterior probability distribution. As the training and testing procedure was performed separately for each hemisphere, a group-level log-posterior probability distribution was first generated for each hemisphere and then these distributions were averaged across hemispheres. The bilateral group-level log-posterior probability distribution was then exponentiated and divided by a normalization constant to compute the geometric mean and yield a probability distribution. To assess the degree to which ROI activity patterns contained information about the presented stimulus, we performed non-parametric permutation tests comparing the probability of the presented stimulus [corresponding with p(0°|b) for circularly shifted distributions] to an empirical null distribution (Golland and Fischl, 2003). To estimate a null distribution, we permuted the motion direction labels within each run before model fitting and testing. This procedure was repeated 1000 times for each subject, yielding a group-level probability distribution for each permutation. We then compared the empirical null distribution to the probability obtained with the true class labeling (α < 0.05, one-sided). A p value was computed as [(# of permutation probabilities ≥ actual probability) + 1]/(N permutations + 1) (Phipson and Smyth, 2010). Correction for multiple comparisons across ROIs was performed using the Holm–Bonferroni procedure (Holm, 1979). Inspired by cluster correction approaches for group-level fMRI analyses (Heller et al., 2006; Woo et al., 2014; Eklund et al., 2016), time-series decoding comparisons were corrected using a nonparametric cluster correction procedure to account for the dependency across sliding TR windows. A cluster-defining threshold was set at p < 0.05. Adjacent timepoints with below-threshold p values (<0.05) were considered a cluster. Permutation probabilities were also converted to p values by comparing each permutation to the remaining permutations and the actual probability, i.e., count number of times permuted probabilities or the actual probability is greater than or equal to the current permutation and then divide by the number of permutations plus 1. We then created a maximum null cluster size distribution by thresholding and recording the maximum cluster size for each permutation across TR windows (Nichols and Holmes, 2002). A p value for the observed cluster size was computed as [(# of permutation cluster sizes ≥ actual cluster size) + 1]/(N permutations + 1) (Phipson and Smyth, 2010). It should be noted that permutation tests impose a minimum p value that is determined by the number of permutations performed. Here, the minimum possible p value before correction for multiple comparisons was 1/1001 = 0.000999. The minimum p value was greater following multiple comparison correction. We additionally examined whether the posterior probability of the presented stimulus given an area's activity patterns was correlated between ROIs across trials. To do so, we computed pairwise correlations between trial-wise estimates of the posterior probability of the presented stimulus for lobule VIIb/VIIIa, IPS0-3, SFJ, and V1-3. This analysis can be likened to a β-series correlation analysis (Rissman et al., 2004) that examines the correspondence between areas in terms of representational content rather than task-evoked BOLD response. Permutation tests were performed by repeatedly shuffling trial-wise posterior probabilities within each subject before computing the correlation between areas. We further investigated the relationship between area decoding correlations for each subject and parameters derived from the variable precision mixture model fitted to each subject's response error distribution. We again performed permutation tests by independently shuffling subject-wise decoding correlations and model parameters before computing the correlation between the two.

Results

To investigate VWM stimulus specificity in the cerebellum, we had participants perform a delayed motion direction recall task while measuring BOLD responses using fMRI (Fig. 1A). On each trial, participants were simultaneously presented with two patches of dot motion (100% coherence), presented to the right and left of fixation. After stimulus offset, participants were retroactively cued to maintain the direction of motion of one of the patches over a blank delay interval (9.75 s) or drop both items. Following the delay period, participants then adjusted a probe stimulus to match the remembered direction of motion. Behavioral performance was assessed by computing a distribution of recall errors over trials for each subject (Fig. 1B). Trial errors were found to be clustered around the cued direction of motion, indicating that participants were able to successfully maintain the direction of motion of the cued item. The average mean absolute value of the angular difference between reported and actual motion directions was 16.75 ± 1.09° SEM. We additionally fit a computational model to participants' response error distribution that enabled us to estimate for each subject the proportion of trials in which they guessed, the mean precision with which they remembered stimuli, and the degree to which precision was variable across trials (see methods; van den Berg et al., 2012; Fougnie, Suchow and Alvarez, 2012). The average guess rate (g) across subjects was 9.29 ± 2.64% SEM. The average precision (SDµ) across subjects was 8.04 ± 1.41° and the average standard deviation of precision (SDσ) was 16.67 ± 2.28°.

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

Behavioral task paradigm and behavioral performance. A, Participants were presented with two patches of moving dots. 250 ms following the offset of the stimulus presentation period participants were postcued to remember the direction of motion of the moving dots. Following a 10-s delay, participants were given 5 s to adjust an oriented line segment to match the remembered direction of motion with key presses. The initial orientation of the line segment was randomized with respect to the cued direction of motion on each trial. B, Kernel density estimates of recall error for all participants. The colored lines denote density estimates of recall error for each individual participant and the black line reflects the density estimate of recall error across all participants.

To ensure that participants were sensitive to the actual stimulus direction on each trial and not simply storing a categorical representation (i.e., 0°, 45°, 90°, 135°, 180°, 225°, 270°, or 315°), we investigated the relationship between the directional jitter applied on each trial (±1–10°) and response error. To do so, we examined recall errors for five categories of trials based on jitter direction and size [zero jitter, small clockwise jitter, small counterclockwise jitter (<5°), large clockwise jitter, and large counterclockwise jitter (>5°)]. We did not find an overall effect of jitter size/direction on signed error, F(4,1912) = 1.96, p = 0.097. Pairwise comparisons between each jitter size/direction and zero jitter further found no significant differences (all p > 0.37), indicating that participants did not exhibit any directional bias in their responses relative to non-jittered trials. We additionally found no differences in absolute error between jitter conditions (all p > 0.43), which showed that recall precision did not depend on the direction or size of jitter applied and therefore that subjects were encoding the actual direction of motion rather than a categorical representation of the non-jittered direction.

Stimulus encoding model

To quantify feature selectivity during working memory, we employed a generative model-based decoding approach. This analysis uses information about voxel motion direction preferences to produce a posterior probability distribution over motion direction for each trial given the observed pattern of voxel responses during the delay period (see methods). To produce a group-level reconstruction of the probability of a particular stimulus given the observed responses, we circularly shifted trial posterior probability distributions so 0° corresponded with the presented motion direction and then averaged the shifted distributions across trials and participants. If an ROI contains feature-selective information about the remembered direction of motion then this averaged posterior distribution should peak at or near 0°. We first present our decoding results for the last three TRs of the delay period (TRs beginning 6–10 s after trial onset), consistent with previous working memory decoding studies (Harrison and Tong, 2009; Ester et al., 2015). Centered and averaged posterior distributions for V1-3, IPS0-3, and SFJ are shown in Figure 2. Each cortical ROI was found to contain a robust representation of the direction of motion maintained in working memory (IPS0-3: p = 0.004 corrected; SFJ: p = 0.004 corrected; V1-3: p = 0.004 corrected; minimum possible p value following correction was 0.004; Fig. 2A–C). The representation of the non-remembered stimulus in these areas also appears to persist into the delay period but in a weaker form (IPS0-3: p = 0.004 corrected, SFJ: p = 0.044 corrected; V1-3: p = 0.027 corrected; Fig. 2D–F). The difference between remembered and non-remembered stimuli for the average across the last three TRs of the delay period is significant for SFJ (p = 0.027 corrected). Although numerically greater the difference for IPS0-3 and V1-3 does not reach significance (p = 0.131 uncorrected and p = 0.129 uncorrected). However, a difference does emerge at the following TR window (average across TRs 5, 6, and 7; IPS: p = 0.079 corrected; V1-3: p = 0.002 corrected). These results replicate prior reports of working memory content specificity in visual, parietal and frontal cortex (Christophel et al., 2012, 2017; Riggall and Postle, 2012; Lee et al., 2013, 2013; Sprague et al., 2014; Ester et al., 2015).

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

Group-level cortical encoding model results. Shifted and averaged posterior probability distributions over motion direction indicating which stimulus is most probable given the observed pattern of vertex responses during the delay period for (A) IPS0-3, (B) SFJ, and (C) V1-3. Posterior probability over motion direction for non-remembered stimuli for (D) IPS0-3, (E) SFJ, and (F) V1-3. Shaded ribbon reflects a bootstrap estimate of SEM.

We next examined whether cerebellar lobule VIIb/VIIIa delay-period responses were also selective for a stimulus stored in working memory. lobule VIIb/VIIIa was preselected for analysis as our prior work reported VWM task activation localized to this portion of cerebellum (Brissenden et al., 2016, 2018). A group-level reconstruction of the probability of a stimulus given the pattern of voxel responses in lobule VIIb/VIIIa was found to peak at the remembered direction of motion (p = 0.004 corrected; Fig. 3A). We also found a significant, albeit weak, representation of the non-remembered stimulus in lobule VIIb/VIIIa (p = 0.039 uncorrected; Fig. 3B). To test whether the strength of representation was stronger for remembered or non-remembered stimuli we compared the observed difference in probability for remembered and non-remembered items to a null probability difference distribution. The probability of the remembered stimulus given the observed pattern of lobule VIIb/VIIIa voxel responses was significantly greater than the probability of the non-remembered stimulus (p = 0.0025 corrected). This difference indicates that the representation of the remembered stimulus in VIIb/VIIIa during the delay period is working memory specific and cannot be explained by residual stimulus-driven responses. We also examined the hemispheric specificity of working memory representations in lobule VIIb/VIIIa by performing the same analysis for stimuli presented in the non-preferred hemifield (contralateral hemifield for cerebellum). We found a significant representation of the remembered item in the contralateral hemifield in lobule VIIb/VIIIa (p = 0.004). However, this representation did not significantly differ from the non-cued item (p = 0.788). Thus, in contrast to items presented in the preferred hemifield (ipsilateral) we are unable to claim that the VIIb/VIIIa representation of items in the non-preferred hemifield is specific to VWM.

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

Lobule VIIb/VIIIa encoding model results. A, Posterior probability over motion direction for remembered stimuli. B, Posterior probability over motion direction for non-remembered stimuli. Shaded ribbon reflects bootstrap SEM. C, Time-resolved encoding model results. Each row of the heat map represents an averaged and shifted probability distribution generated by applying an encoding model to an average of three consecutive TRs. tr123 – TRs 1, 2, 3; tr234 – TRs 2, 3, 4; tr345 – TRs 3, 4, 5; tr456 – TRs 4, 5, 6; tr567 – TRs 5, 6, 7; tr678 – TRs 6, 7, 8.

We next investigated the relationship between decoding in lobule VIIb/VIIIa and cerebral cortical regions. To do so, we extracted the probability of the presented stimulus from shifted posterior distributions for each trial [i.e., p(0°|b) for circularly shifted distributions] and then correlated these values between areas for each subject. lobule VIIb/VIIIa decoding was positively correlated with decoding in each cortical ROI and permutation tests indicated that correlations averaged across subjects were greater than what would be expected by chance (VIIb/VIIIa–IPS0-3: r = 0.171 ± 0.052, p = 0.006 corrected; VIIb/VIIIa–SFJ: r = 0.142 ± 0.050, p = 0.006 corrected; VIIb/VIIIa–V1-3: r = 0.148 ± 0.038, p = 0.006 corrected; 0.006 was the minimum possible corrected p value for 6 comparisons). Furthermore, the observed trial-wise correlations between lobule VIIb/VIIIa and cortical ROIs were comparable to those between cortical ROIs (IPS0-3–SFJ: r = 0.281 ± 0.037, p = 0.006 corrected; IPS0-3–V1-3: r = 0.161 ± 0.039, p = 0.006 corrected; SFJ–V1-3: r = 0.101 ± 0.025, p = 0.006 corrected). Taken together, these results indicate that fluctuations in information content within lobule VIIb/VIIIa and canonical VWM areas are linked on trial-wise basis.

We also found that the correlation in decoding between certain areas predicted behavioral recall precision. The correlation between IPS0-3 and lobule VIIb/VIIIa and the correlation between IPS0-3 and SFJ for each participant were found to be significantly correlated with the mean precision parameter from the variable precision model fitted to each participant's error distribution [cor(IPS0-3–VIIb/VIIIa, SDµ): r = −0.59, p = 0.042 corrected; cor(IPS0-3–SFJ, SDµ): r = −0.61, p = 0.04995 corrected]. No significant correlations were identified between the other model parameters (g and SDσ) and the correlation in decoding between areas. These results indicate that subjects for which trial-by-trial fluctuations in information content were more tightly linked between IPS and both cerebellum and frontal cortex remembered stimuli with greater precision.

Time-resolved lobule VIIb/VIIIa stimulus probability reconstructions for both remembered and non-remembered items were produced by applying the encoding model to a sliding window consisting of three consecutive TRs (Fig. 3C). A non-parametric cluster-based correction procedure, a method commonly used to correct for multiple correlated comparisons in group-level fMRI analyses (Heller et al., 2006; Woo et al., 2014; Eklund et al., 2016), revealed a significant effect from the second to fifth TR window for remembered items (cluster p = 0.003). On the other hand, the non-remembered item was not significantly represented across time (p = 0.224). Furthermore, the difference in probability between remembered and non-remembered items was significantly different from the second to fifth TR window (p = 0.002). Thus, cerebellar lobule VIIb/VIIIa delay period responses appear to be content-specific and this specificity is sustained throughout the delay period.

Working memory stimulus specificity was further shown to be restricted to lobule VIIb/VIIIa within the cerebellum. We were unable to reconstruct the remembered stimulus from any other cerebellar lobule (all p > 0.05 uncorrected), with the exception of Crus II (p = 0.027 uncorrected), which has been previously implicated in working memory (Guell et al., 2018a), exhibits robust functional connectivity with the cognitive control network (Buckner et al., 2011; Guell et al., 2018b), and is adjacent to lobule VIIb/VIIIa (Fig. 4). However, Crus II only weakly represented the remembered stimulus and this effect did not survive correction for multiple comparisons (p = 0.216 corrected). We additionally examined whether an oculomotor vermis ROI, an area known to be highly responsive to saccades (Noda and Fujikado, 1987; Ohtsuka and Noda, 1995; Takagi et al., 1998; Thier et al., 2002; Soetedjo and Fuchs, 2006), could decode the remembered stimulus to determine whether saccadic eye movements could explain the observed stimulus selective activity. We were unable to recover a significant representation of the presented stimulus from oculomotor vermis delay period activity (p = 0.098 uncorrected; Fig. 5). Collectively, these results indicate that the encoding of stimulus-specific mnemonic representations is limited to lobule VIIb/VIIIa within the cerebellum. Furthermore, the lack of evidence for feature-selective representations elsewhere in the cerebellum argues against a purely motor explanation for lobule VIIb/VIIIa content specificity, as motor processes would be expected to drive areas of the cerebellum previously implicated in motor functions such as the anterior lobe (lobules I–V) and the oculomotor vermis.

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

Cerebellar lobule encoding model results. Group-level posterior probability distributions for cerebellar lobules outside of Lobule VIIb/VIIIa. Inset shows flatmap representation of cerebellum (Diedrichsen and Zotow, 2015). Colors and dotted lines demarcate cerebellar lobules. Shaded ribbon reflects bootstrap SEM.

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

Oculomotor vermis encoding model results. Averaged and shifted group-level posterior probability distribution for oculomotor vermis. Inset depicts oculomotor vermis ROI overlaid on cerebellar flatmap. Shaded ribbon reflects bootstrap SEM.

We next sought to determine whether a particular subdivision of lobule VIIb/VIIIa was driving the observed working memory stimulus specificity. To do so, we performed a seed-based resting-state functional connectivity analysis using V1-3, IPS0-3, and SFJ ROIs as seeds. As each cerebral cortical hemisphere is preferentially connected with the contralateral hemisphere of cerebellum, these analyses were performed separately for each hemisphere (right hemisphere cerebral cortical ROI to left hemisphere cerebellum and vice versa). For each subject, we then identified the top 25% of correlations within lobule VIIb/VIIIa for each cortical ROI. We then computed the overlap across subjects by summing individual percentile ROIs across subjects for each resting-state analysis (Fig. 6). Cortical ROIs exhibited differential connectivity across lobule VIIb/VIIIa. IPS0-3 was preferentially connected with the medial portion of lobule VIIb across subjects, and to a lesser extent the lateral portion of lobule VIIb/VIIIa (Fig. 6A), while largely sparing the central portion. SFJ connectivity, on the other hand, was equally concentrated in medial VIIb and lateral VIIb/VIIIa, again sparing the central portion, across subjects (Fig. 6B). V1-3 connectivity was more evenly distributed across all of lobule VIIb/VIIIa across subjects (Fig. 6C). Using these subject count maps, we then subdivided lobule VIIb/VIIIa to create ROIs localized to the medial, central, and lateral portions of lobule VIIb/VIIIa bilaterally (see Materials and Methods; Fig. 7A). We then applied the encoding model to each of these ROIs' delay-period activity to determine which portion of lobule VIIb/VIIIa exhibited feature-selection activation. Medial lobule VIIb/VIIIa was shown to encode feature-selective information about the remembered stimulus (p = 0.004 corrected; Fig. 7B), but little information about the non-remembered stimulus (p = 0.09 uncorrected; Fig. 7C). Additionally, medial lobule VIIb/VIIIa exhibited a significant difference in decoding between remembered and non-remembered stimuli (p = 0.02 corrected). Lateral lobule VIIb/VIIIa also exhibited significant decoding of the remembered stimulus (p = 0.004 corrected; Fig. 7F), but was also shown to significantly decode the non-remembered stimulus (p = 0.036 corrected; Fig. 7G). The difference in decoded probability between remembered and non-remembered stimuli was not significant (p = 0.169 uncorrected), indicating that the stimulus representation in lateral lobule VIIb/VIIIa may be in part driven by a perceptual trace rather than a working memory trace. Nevertheless, combining the medial and lateral VIIb/VIIIa ROIs into a single ROI yielded significant decoding of the remembered stimulus (p = 0.004 corrected; Fig. 7H), but not the non-remembered stimulus (p = 0.122 uncorrected; Fig. 7I), and an even more robust difference between remembered and non-remembered stimuli (p = 0.0002 corrected). In contrast, the central portion in each hemisphere of lobule VIIb/VIIIa did not exhibit evidence for a representation of either the remembered or non-remembered stimulus (remembered: p = 0.882 uncorrected; non-remembered: p = 0.111 uncorrected; difference: p = 0.953 uncorrected; Fig. 7D,E). In sum, portions of lobule VIIb/VIIIa that are preferentially connected with fronto-parietal cortex exhibit content-specific activation during the delay-period, while portions of lobule VIIb/VIIIa that are preferentially connected with visual cortex do not.

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

Cerebellar Lobule VIIb/VIIIa functional connectivity with cortical ROIs. The top 25% of correlations across Lobule VIIb/VIIIa with either (A) IPS0-3, (B) SFJ, or (C) V1-3 were identified for each participant. Values within Lobule VIIb/VIIIa reflect the number of subjects for which that vertex is included in the top 25% of correlations.

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

Cerebellar Lobule VIIb/VIIIa subdivision encoding model results. A, Medial, central, and lateral VIIb/VIIIa ROIs defined using probabilistic functional connectivity subject count maps depicted in Figure 6. B, C, Medial Lobule VIIb/VIIIa ROI group-level posterior probability distribution for (B) remembered and (C) non-remembered items. D, E, Central Lobule VIIb/VIIIa posterior distribution for (D) remembered and (E) non-remembered items. F, G, Lateral Lobule VIIb/VIIIa posterior distribution for (F) remembered and (G) non-remembered items. H, I, Combined medial and lateral Lobule VIIb/VIIIa posterior distribution for (H) remembered and (I) non-remembered items. Shaded ribbon reflects bootstrap SEM.

Discussion

The findings presented here demonstrate for the first time that cerebellar lobule VIIb/VIIIa encodes VWM representations. Prior research had shown that lobule VIIb/VIIIa is recruited by VWM and attention tasks (Chen and Desmond, 2005b; Kirschen et al., 2005; Brissenden et al., 2016, 2018), yet the extent to which lobule VIIb/VIIIa is involved in working memory storage processes or rather more generalized attentional processes was previously unclear. Here, we used a retro-cue paradigm to examine working memory storage independent of perception-related, attention-related, and retrieval-related processes. An encoding model of motion direction revealed that lobule VIIb/VIIIa robustly represented the remembered stimulus. Moreover, the representation of the remembered stimulus was significantly stronger than the representation of the non-remembered stimulus, indicating that this representation can be attributed to a working memory trace rather than a perceptual one. Delay-period lobule VIIb/VIIIa stimulus-selectivity was further shown to be robustly correlated with selectivity in cortical regions previously implicated in VWM. Notably, the correlation in decoding between IPS and lobule VIIb/VIIIa was strongly related to behavioral recall precision. Among cerebellar structures, stimulus representation during the delay period was unique to lobule VIIb/VIIIa. We were unable to recover the remembered direction of motion from other portions of the cerebellum. Further analysis demonstrated that specific subdivisions of lobule VIIb/VIIIa exhibit stimulus-specific responses. Medial VIIb and lateral VIIb/VIIIa were shown to be preferentially connected with frontoparietal cortex, while the intermediate area was found to be preferentially connected with early visual cortex. Subsequent analysis demonstrated that both medial and lateral portions of lobule VIIb/VIIIa encoded remembered items, while the intervening central portion did not. It is important to note that these same medial and lateral portions have been shown to contain topographic representations of the visual field (Brissenden et al., 2018; van Es et al., 2019). Taken together, these findings indicate that stimuli stored in VWM are represented by the cerebellar component of the dorsal attention network (Brissenden et al., 2016).

The current results provide further evidence for cerebellar contributions to cognitive function. It has become increasingly clear that cerebellar function is not limited to the motor domain and that substantial functional heterogeneity exists within cerebellar cortex (Stoodley and Schmahmann, 2009; Buckner et al., 2011; Stoodley et al., 2012; Brissenden et al., 2016, 2018; Guell et al., 2018a,b; van Es et al., 2019; King et al., 2019). Along with recent work demonstrating the existence of visual field representations in lobule VIIb/VIIIa (Brissenden et al., 2018; van Es et al., 2019), this study demonstrates that the cerebellum possesses considerable representational specificity for items within the focus of attention or working memory. It remains to be seen whether areas of the cerebellum associated with other cognitive domains exhibit similar specificity in terms of their representational content.

A number of prior studies have provided evidence for the recruitment of the cerebellum by both visual and verbal working memory (Desmond et al., 1997; Chen and Desmond, 2005b; Kirschen et al., 2005; Hayter et al., 2007; Stoodley and Schmahmann, 2009; Stoodley et al., 2012; Brissenden et al., 2016; Guell et al., 2018a; King et al., 2019). These studies revealed activation across a broad swath of cerebellar cortex including lobule VI, Crus I, Crus II, lobule VIIb, and lobule VIIIa in response to a variety of working memory paradigms. This body of research includes our prior work in which we showed that the portion of lobule VIIb/VIIIa that is functionally connected with cortical dorsal attention network areas responds in a load-dependent manner to a VWM change detection task (Brissenden et al., 2016). However, as the methods used in these prior studies did not lend themselves to revealing content-specific activity, it was previously unclear whether the cerebellum could be specifically implicated in the storage of VWM representations. In the present study, we employed an event-related retro-cue paradigm, which enabled the dissociation of activity related to remembered and non-remembered items during the maintenance period. Consequently, we were able to show that specific cerebellar areas encoded the identity of a stimulus during the delay-period and that this selectivity was substantially more robust for remembered items. These findings extend the prior literature on cerebellar contributions to working memory to provide novel evidence that the cerebellum is actively involved in the persistent storage of information in working memory.

A substantial body of work has attempted to identify the locus or loci of a working memory storage buffer. Our results extend recent proposals that working memory contents are distributed across a number of areas (Lee and Baker, 2016; Serences, 2016; Christophel et al., 2017) by showing that working memory representations are additionally encoded in the cerebellum. Central to the distributed working memory network hypothesis is the notion that storage across different areas reflects different levels of representational abstraction (Christophel et al., 2017). Recent work has proposed a predictive role for cerebellum in verbal working memory (Sheu et al., 2019). Motor control studies suggest that cerebellum serves as a “forward-model” that predicts actions (Wolpert et al., 1998). Given the notable homogeneity of cerebellar circuit architecture, similar computations may be conducted by cognitive portions of the cerebellum (Ramnani, 2006). Follow-up work will need to more finely probe the nature of representations in cerebellum to elucidate its unique role in VWM storage processes.

As the cerebellum has been extensively implicated in oculomotor processes (Ron and Robinson, 1973; Fujikado and Noda, 1987; Voogd and Barmack, 2006; Baier et al., 2009), a candidate alternative explanation for our findings is that VIIb/VIIIa activity patterns reflect eye movement processes rather than a VWM representation. This interpretation is unlikely for several reasons. First, we could not decode remembered motion direction from oculomotor vermis, an area that has been extensively implicated in the control of eye movements (Noda and Fujikado, 1987; Ohtsuka and Noda, 1995; Takagi et al., 1998; Thier et al., 2002; Soetedjo and Fuchs, 2006). Second, participants were experienced observers who extensively practiced maintaining fixation before scanning. Although hardware problems prevented us from tracking eye position in the present study, in our prior work with lateralized VWM targets, subjects exhibited no difference in either eye position or movement between “remember-left” and “remember-right” conditions (Brissenden et al., 2018). Third, our findings cannot be explained by the two most obvious saccade strategies. One strategy would be to make eye movements during stimulus presentation along the direction of motion. However, as participants did not know which dot motion patch would be cued, there is no reason that eye movements would be selective for the direction of motion of the subsequently cued item and not the non-cued item. Another strategy would be to make a saccade to the cued hemifield following the presentation of the cue. Yet, as motion direction and hemifield were dissociated from one another, a saccade toward the cued hemifield would not be predictive of stimulus identity. A less obvious strategy would be to refrain from making eye movements during the presentation of coherent motion and then make a saccade parallel to the direction of motion of the cued item following the presentation of the cue and the offset of the motion stimulus. lobule VIIb/VIIIa stimulus specificity could potentially be explained by this scenario. Concurrent eye-tracking during fMRI scanning would allow future studies to explicitly rule out this possibility.

A recent study in mice reported that the persistent representation of information in frontal cortex depends on the cerebellum (Gao et al., 2018). Head-fixed mice were presented with a sample stimulus that cued one of two actions (left or right lick) to be made following a delay period. Neurons in both frontal cortex and the cerebellar fastigial nucleus were shown to exhibit selectivity for the cued action during the delay period. Critically, optogenetic silencing of fastigial neurons abolished selectivity in frontal cortex and resulted in incorrect behavioral choices. Genetic manipulation of different cell types in the cortico-cerebellar circuit has the potential to further advance our understanding of delay activity in the cerebellum (Low et al., 2018). In the current study, we explicitly dissociated motor planning from VWM storage by randomizing the start position of the probe stimulus relative to the remembered stimulus. As a consequence, subjects were unable to prospectively plan their responses. Thus, our results suggest that the cerebellum is involved in the sustained representation of information in working memory even when stored items are not related with specific actions. Stimulation protocols such as TMS could be used to investigate whether inactivation of specific cerebellar areas in humans similarly interferes with cortical representations of remembered stimuli, as well as behavioral performance (Halko et al., 2014; Esterman et al., 2017).

Cortico-centric models of cognition are pervasive in cognitive neuroscience. Our findings highlight the shortcomings of a narrow focus on cerebral cortex in characterizing the neural mechanisms of working memory storage. The more expansive characterization of the neural substrates of working memory storage suggested by the current results could provide new insights into a wide range of goal-directed behaviors that are known to rely on working memory.

Footnotes

  • This work was supported by National Institutes of Health Grants R21-EY027703 (to D.C.S.), R01-MH116170 and R01-MH111868 (to M.A.H.), and F31-NS103306 (to S.M.T.), as well as the National Science Foundation Grant DGE-1247312 (to J.A.B.). We thank Dr. Stephanie McMains and Dr. Ross Mair for technical advice.

  • The authors declare no competing financial interests.

  • Correspondence should be addressed to David C. Somers at somers{at}bu.edu or James A. Brissenden at brissend{at}umich.edu

SfN exclusive license.

References

  1. ↵
    1. Akaike H
    (1974) A new look at the statistical model identification. IEEE Trans Automat Contr 19:716–723. doi:10.1109/TAC.1974.1100705
    OpenUrlCrossRef
  2. ↵
    1. Allen G,
    2. Buxton RB,
    3. Wong EC,
    4. Courchesne E
    (1997) Attentional activation of the cerebellum independent of motor involvement. Science 275:1940–1943. doi:10.1126/science.275.5308.1940 pmid:9072973
    OpenUrlAbstract/FREE Full Text
  3. ↵
    1. Baier B,
    2. Stoeter P,
    3. Dieterich M
    (2009) Anatomical correlates of ocular motor deficits in cerebellar lesions. Brain J Brain 132:2114–2124. doi:10.1093/brain/awp165 pmid:19567703
    OpenUrlCrossRefPubMed
  4. ↵
    1. Baier B,
    2. Dieterich M,
    3. Stoeter P,
    4. Birklein F,
    5. Müller NG
    (2010) Anatomical correlate of impaired covert visual attentional processes in patients with cerebellar lesions. J Neurosci 30:3770–3776. doi:10.1523/JNEUROSCI.0487-09.2010 pmid:20220011
    OpenUrlAbstract/FREE Full Text
  5. ↵
    1. Bays PM,
    2. Catalao RFG,
    3. Husain M
    (2009) The precision of visual working memory is set by allocation of a shared resource. J Vision 9:1–11.
    OpenUrlAbstract
  6. ↵
    1. Brissenden JA,
    2. Levin EJ,
    3. Osher DE,
    4. Halko MA,
    5. Somers DC
    (2016) Functional evidence for a cerebellar node of the dorsal attention network. J Neurosci 36:6083–6096. doi:10.1523/JNEUROSCI.0344-16.2016 pmid:27251628
    OpenUrlAbstract/FREE Full Text
  7. ↵
    1. Brissenden JA,
    2. Tobyne SM,
    3. Osher DE,
    4. Levin EJ,
    5. Halko MA,
    6. Somers DC
    (2018) Topographic cortico-cerebellar networks revealed by visual attention and working memory. Curr Biol 28:3364–3372.e5. doi:10.1016/j.cub.2018.08.059 pmid:30344119
    OpenUrlCrossRefPubMed
  8. ↵
    1. Buckner RL,
    2. Krienen FM,
    3. Castellanos A,
    4. Diaz JC,
    5. Yeo BTT
    (2011) The organization of the human cerebellum estimated by intrinsic functional connectivity. J Neurophysiol 106:2322–2345. doi:10.1152/jn.00339.2011 pmid:21795627
    OpenUrlCrossRefPubMed
  9. ↵
    1. Burgess GC,
    2. Kandala S,
    3. Nolan D,
    4. Laumann TO,
    5. Power JD,
    6. Adeyemo B,
    7. Harms MP,
    8. Petersen SE,
    9. Barch DM
    (2016) Evaluation of denoising strategies to address motion-correlated artifacts in resting-state functional magnetic resonance imaging data from the human connectome project. Brain Connect 6:669–680. doi:10.1089/brain.2016.0435 pmid:27571276
    OpenUrlCrossRefPubMed
  10. ↵
    1. Carp J
    (2013) Optimizing the order of operations for movement scrubbing: comment on Power et al. Neuroimage 76:436–438. doi:10.1016/j.neuroimage.2011.12.061 pmid:22227884
    OpenUrlCrossRefPubMed
  11. ↵
    1. Cauley SF,
    2. Polimeni JR,
    3. Bhat H,
    4. Wald LL,
    5. Setsompop K
    (2014) Interslice leakage artifact reduction technique for simultaneous multislice acquisitions. Magnet Reson Med 72:93–102.
    OpenUrl
  12. ↵
    1. Chafee MV,
    2. Goldman-Rakic PS
    (1998) Matching patterns of activity in primate prefrontal area 8a and parietal area 7ip neurons during a spatial working memory task. J Neurophysiol 79:2919–2940. doi:10.1152/jn.1998.79.6.2919 pmid:9636098
    OpenUrlCrossRefPubMed
  13. ↵
    1. Chen SHA,
    2. Desmond JE
    (2005a) Temporal dynamics of cerebro-cerebellar network recruitment during a cognitive task. Neuropsychologia 43:1227–1237. doi:10.1016/j.neuropsychologia.2004.12.015 pmid:15949507
    OpenUrlCrossRefPubMed
  14. ↵
    1. Chen SHA,
    2. Desmond JE
    (2005b) Cerebrocerebellar networks during articulatory rehearsal and verbal working memory tasks. Neuroimage 24:332–338. doi:10.1016/j.neuroimage.2004.08.032 pmid:15627576
    OpenUrlCrossRefPubMed
  15. ↵
    1. Christophel TB,
    2. Hebart MN,
    3. Haynes JD
    (2012) Decoding the contents of visual short-term memory from human visual and parietal cortex. J Neurosci 32:12983–12989. doi:10.1523/JNEUROSCI.0184-12.2012 pmid:22993415
    OpenUrlAbstract/FREE Full Text
  16. ↵
    1. Christophel TB,
    2. Cichy RM,
    3. Hebart MN,
    4. Haynes JD
    (2015) Parietal and early visual cortices encode working memory content across mental transformations. Neuroimage 106:198–206. doi:10.1016/j.neuroimage.2014.11.018 pmid:25463456
    OpenUrlCrossRefPubMed
  17. ↵
    1. Christophel TB,
    2. Klink PC,
    3. Spitzer B,
    4. Roelfsema PR,
    5. Haynes JD
    (2017) The distributed nature of working memory. Trends Cogn Sci 21:111–124. doi:10.1016/j.tics.2016.12.007 pmid:28063661
    OpenUrlCrossRefPubMed
  18. ↵
    1. Christophel TB,
    2. Iamshchinina P,
    3. Yan C,
    4. Allefeld C,
    5. Haynes JD
    (2018) Cortical specialization for attended versus unattended working memory. Nat Neurosci 21:494–496. doi:10.1038/s41593-018-0094-4 pmid:29507410
    OpenUrlCrossRefPubMed
  19. ↵
    1. Courtney SM,
    2. Petit L,
    3. Maisog JM,
    4. Ungerleider LG,
    5. Haxby JV
    (1998) An area specialized for spatial working memory in human frontal cortex. Science 279:1347–1351. doi:10.1126/science.279.5355.1347 pmid:9478894
    OpenUrlAbstract/FREE Full Text
  20. ↵
    1. Desmond JE,
    2. Gabrieli JDE,
    3. Wagner AD,
    4. Ginier BL,
    5. Glover GH
    (1997) Lobular patterns of cerebellar activation in verbal working-memory and finger-tapping tasks as revealed by functional MRI. J Neurosci 17:9675–9685. doi:10.1523/JNEUROSCI.17-24-09675.1997 pmid:9391022
    OpenUrlAbstract/FREE Full Text
  21. ↵
    1. D'Esposito M,
    2. Postle BR
    (2015) The cognitive neuroscience of working memory. Annu Rev Psychol 66:1–28.
    OpenUrlCrossRefPubMed
  22. ↵
    1. Diedrichsen J,
    2. Zotow E
    (2015) Surface-based display of volume-averaged cerebellar imaging data. PLoS One 10:e0133402. doi:10.1371/journal.pone.0133402 pmid:26230510
    OpenUrlCrossRefPubMed
  23. ↵
    1. Diedrichsen J,
    2. Balsters J,
    3. Flavell J,
    4. Cussans E,
    5. Ramnani N
    (2009) A probabilistic MR atlas of the human cerebellum. Neuroimage 47:S122. doi:10.1016/S1053-8119(09)71166-8
    OpenUrlCrossRef
  24. ↵
    1. Eklund A,
    2. Nichols TE,
    3. Knutsson H
    (2016) Cluster failure: why fMRI inferences for spatial extent have inflated false-positive rates. Proc Natl Acad Sci USA 113:7900–7905. doi:10.1073/pnas.1602413113 pmid:27357684
    OpenUrlAbstract/FREE Full Text
  25. ↵
    1. Engle RW
    (2002) Working memory capacity as executive attention. Curr Dir Psychol Sci 11:19–23. doi:10.1111/1467-8721.00160
    OpenUrlCrossRef
  26. ↵
    1. Engle RW,
    2. Tuholski SW,
    3. Laughlin JE,
    4. Conway ARA
    (1999) Working memory, short-term memory, and general fluid intelligence: a latent-variable approach. J Exp Psychol Gen 128:309–331. doi:10.1037//0096-3445.128.3.309 pmid:10513398
    OpenUrlCrossRefPubMed
  27. ↵
    1. Ester EF,
    2. Sprague TC,
    3. Serences JT
    (2015) Parietal and frontal cortex encode stimulus-specific mnemonic representations during visual working memory. Neuron 87:893–905. doi:10.1016/j.neuron.2015.07.013 pmid:26257053
    OpenUrlCrossRefPubMed
  28. ↵
    1. Esterman M,
    2. Thai M,
    3. Okabe H,
    4. DeGutis J,
    5. Saad E,
    6. Laganiere SE,
    7. Halko MA
    (2017) Network-targeted cerebellar transcranial magnetic stimulation improves attentional control. Neuroimage 156:190–198. doi:10.1016/j.neuroimage.2017.05.011 pmid:28495634
    OpenUrlCrossRefPubMed
  29. ↵
    1. Feinberg DA,
    2. Moeller S,
    3. Smith SM,
    4. Auerbach E,
    5. Ramanna S,
    6. Gunther M,
    7. Glasser MF,
    8. Miller KL,
    9. Ugurbil K,
    10. Yacoub E
    (2010) Multiplexed echo planar imaging for sub-second whole brain FMRI and fast diffusion imaging. PLoS One 5:e15710. doi:10.1371/journal.pone.0015710 pmid:21187930
    OpenUrlCrossRefPubMed
  30. ↵
    1. Fleming PJ,
    2. Wallace JJ
    (1986) How not to lie with statistics: the correct way to summarize benchmark results. Commun Acm 29:218–221. doi:10.1145/5666.5673
    OpenUrlCrossRef
  31. ↵
    1. Fougnie D,
    2. Suchow JW,
    3. Alvarez GA
    (2012) Variability in the quality of visual working memory. Nat Commun 3:1229. doi:10.1038/ncomms2237 pmid:23187629
    OpenUrlCrossRefPubMed
  32. ↵
    1. Fujikado T,
    2. Noda H
    (1987) Saccadic eye movements evoked by microstimulation of lobule VII of the cerebellar vermis of macaque monkeys. J Physiol 394:573–594. doi:10.1113/jphysiol.1987.sp016885 pmid:3443975
    OpenUrlCrossRefPubMed
  33. ↵
    1. Fukuda K,
    2. Vogel E,
    3. Mayr U,
    4. Awh E
    (2010) Quantity, not quality: the relationship between fluid intelligence and working memory capacity. Psychon Bull Rev 17:673–679. doi:10.3758/17.5.673 pmid:21037165
    OpenUrlCrossRefPubMed
  34. ↵
    1. Funahashi S,
    2. Bruce CJ,
    3. Goldman-Rakic PS
    (1989) Mnemonic coding of visual space in the monkey's dorsolateral prefrontal cortex. J Neurophysiol 61:331–349. doi:10.1152/jn.1989.61.2.331 pmid:2918358
    OpenUrlCrossRefPubMed
  35. ↵
    1. Fuster JM,
    2. Alexander GE
    (1971) Neuron activity related to short-term memory. Science 173:652–654. doi:10.1126/science.173.3997.652 pmid:4998337
    OpenUrlAbstract/FREE Full Text
  36. ↵
    1. Gao Z,
    2. Davis C,
    3. Thomas AM,
    4. Economo MN,
    5. Abrego AM,
    6. Svoboda K,
    7. Zeeuw CID,
    8. Li N
    (2018) A cortico-cerebellar loop for motor planning. Nature 563:113–116. doi:10.1038/s41586-018-0633-x pmid:30333626
    OpenUrlCrossRefPubMed
  37. ↵
    1. Glasser MF,
    2. Sotiropoulos SN,
    3. Wilson JA,
    4. Coalson TS,
    5. Fischl B,
    6. Andersson JL,
    7. Xu J,
    8. Jbabdi S,
    9. Webster M,
    10. Polimeni JR,
    11. Essen DCV,
    12. Jenkinson M
    ; WU-Minn HCP Consortium (2013) The minimal preprocessing pipelines for the human connectome project. Neuroimage 80:105–124. doi:10.1016/j.neuroimage.2013.04.127 pmid:23668970
    OpenUrlCrossRefPubMed
  38. ↵
    1. Glasser MF,
    2. Coalson TS,
    3. Robinson EC,
    4. Hacker CD,
    5. Harwell J,
    6. Yacoub E,
    7. Ugurbil K,
    8. Andersson J,
    9. Beckmann CF,
    10. Jenkinson M,
    11. Smith SM,
    12. Essen DCV
    (2016) A multi-modal parcellation of human cerebral cortex. Nature 536:171–178. doi:10.1038/nature18933 pmid:27437579
    OpenUrlCrossRefPubMed
  39. ↵
    1. Golland P,
    2. Fischl B
    (2003) Permutation tests for classification: towards statistical significance in image-based studies. Inf Process Med Imaging Proc 18:330–341.
    OpenUrl
  40. ↵
    1. Guell X,
    2. Gabrieli JDE,
    3. Schmahmann JD
    (2018a) Triple representation of language, working memory, social and emotion processing in the cerebellum: convergent evidence from task and seed-based resting-state fMRI analyses in a single large cohort. bioRxiv 254110.
  41. ↵
    1. Guell X,
    2. Schmahmann JD,
    3. Gabrieli JD,
    4. Ghosh SS
    (2018b) Functional gradients of the cerebellum. Elife 7:e36652. doi:10.7554/eLife.36652
    OpenUrlCrossRefPubMed
  42. ↵
    1. Halko MA,
    2. Farzan F,
    3. Eldaief MC,
    4. Schmahmann JD,
    5. Pascual-Leone A
    (2014) Intermittent theta-burst stimulation of the lateral cerebellum increases functional connectivity of the default network. J Neurosci 34:12049–12056. doi:10.1523/JNEUROSCI.1776-14.2014 pmid:25186750
    OpenUrlAbstract/FREE Full Text
  43. ↵
    1. Harrison SA,
    2. Tong F
    (2009) Decoding reveals the contents of visual working memory in early visual areas. Nature 458:632–635. doi:10.1038/nature07832 pmid:19225460
    OpenUrlCrossRefPubMed
  44. ↵
    1. Hayter AL,
    2. Langdon DW,
    3. Ramnani N
    (2007) Cerebellar contributions to working memory. Neuroimage 36:943–954. doi:10.1016/j.neuroimage.2007.03.011 pmid:17468013
    OpenUrlCrossRefPubMed
  45. ↵
    1. Heller R,
    2. Stanley D,
    3. Yekutieli D,
    4. Rubin N,
    5. Benjamini Y
    (2006) Cluster-based analysis of FMRI data. Neuroimage 33:599–608. doi:10.1016/j.neuroimage.2006.04.233 pmid:16952467
    OpenUrlCrossRefPubMed
  46. ↵
    1. Holm S
    (1979) A simple sequentially rejective multiple test procedure. Scand J Stat 6:65–70.
    OpenUrl
  47. ↵
    1. Hurvich CM,
    2. Tsai C-L
    (1989) Regression and time series model selection in small samples. Biometrika 76:297–307. doi:10.1093/biomet/76.2.297
    OpenUrlCrossRef
  48. ↵
    1. King M,
    2. Hernandez-Castillo CR,
    3. Poldrack RA,
    4. Ivry RB,
    5. Diedrichsen J
    (2019) Functional boundaries in the human cerebellum revealed by a multi-domain task battery. Nat Neurosci 22:1371–1378. doi:10.1038/s41593-019-0436-x pmid:31285616
    OpenUrlCrossRefPubMed
  49. ↵
    1. Kirschen MP,
    2. Chen SHA,
    3. Schraedley-Desmond P,
    4. Desmond JE
    (2005) Load- and practice-dependent increases in cerebro-cerebellar activation in verbal working memory: an fMRI study. Neuroimage 24:462–472. doi:10.1016/j.neuroimage.2004.08.036 pmid:15627588
    OpenUrlCrossRefPubMed
  50. ↵
    1. Lee SH,
    2. Baker CI
    (2016) Multi-voxel decoding and the topography of maintained information during visual working memory. Front Syst Neurosci 10:2.
    OpenUrlCrossRefPubMed
  51. ↵
    1. Lee SH,
    2. Kravitz DJ,
    3. Baker CI
    (2013) Goal-dependent dissociation of visual and prefrontal cortices during working memory. Nat Neurosci 16:997–999. doi:10.1038/nn.3452 pmid:23817547
    OpenUrlCrossRefPubMed
  52. ↵
    1. Low AYT,
    2. Thanawalla AR,
    3. Yip AKK,
    4. Kim J,
    5. Wong KLL,
    6. Tantra M,
    7. Augustine GJ,
    8. Chen AI
    (2018) Precision of discrete and rhythmic forelimb movements requires a distinct neuronal subpopulation in the interposed anterior nucleus. Cell Rep 22:2322–2333. doi:10.1016/j.celrep.2018.02.017 pmid:29490269
    OpenUrlCrossRefPubMed
  53. ↵
    1. Mendoza-Halliday D,
    2. Torres S,
    3. Martinez-Trujillo JC
    (2014) Sharp emergence of feature-selective sustained activity along the dorsal visual pathway. Nat Neurosci 17:1255–1262. doi:10.1038/nn.3785 pmid:25108910
    OpenUrlCrossRefPubMed
  54. ↵
    1. Moeller S,
    2. Yacoub E,
    3. Olman CA,
    4. Auerbach E,
    5. Strupp J,
    6. Harel N,
    7. Uğurbil K
    (2010) Multiband multislice GE-EPI at 7 tesla, with 16-fold acceleration using partial parallel imaging with application to high spatial and temporal whole-brain fMRI. Magn Reson Med 63:1144–1153. doi:10.1002/mrm.22361 pmid:20432285
    OpenUrlCrossRefPubMed
  55. ↵
    1. Naselaris T,
    2. Prenger RJ,
    3. Kay KN,
    4. Oliver M,
    5. Gallant JL
    (2009) Bayesian reconstruction of natural images from human brain activity. Neuron 63:902–915. doi:10.1016/j.neuron.2009.09.006 pmid:19778517
    OpenUrlCrossRefPubMed
  56. ↵
    1. Nelson K
    (2017) Assessing probabilistic inference by comparing the generalized mean of the model and source probabilities. Entropy 19:286. doi:10.3390/e19060286
    OpenUrlCrossRef
  57. ↵
    1. Nichols TE,
    2. Holmes AP
    (2002) Nonparametric permutation tests for functional neuroimaging: a primer with examples. Hum Brain Mapp 15:1–25. doi:10.1002/hbm.1058 pmid:11747097
    OpenUrlCrossRefPubMed
  58. ↵
    1. Noda H,
    2. Fujikado T
    (1987) Topography of the oculomotor area of the cerebellar vermis in macaques as determined by microstimulation. J Neurophysiol 58:359–378. doi:10.1152/jn.1987.58.2.359 pmid:3655873
    OpenUrlCrossRefPubMed
  59. ↵
    1. Ohtsuka K,
    2. Noda H
    (1995) Discharge properties of Purkinje cells in the oculomotor vermis during visually guided saccades in the macaque monkey. J Neurophysiol 74:1828–1840. doi:10.1152/jn.1995.74.5.1828 pmid:8592177
    OpenUrlCrossRefPubMed
  60. ↵
    1. Pasternak T,
    2. Greenlee MW
    (2005) Working memory in primate sensory systems. Nat Rev Neurosci 6:97–107. doi:10.1038/nrn1603 pmid:15654324
    OpenUrlCrossRefPubMed
  61. ↵
    1. Peirce JW
    (2007) PsychoPy—psychophysics software in Python. J Neurosci Methods 162:8–13. doi:10.1016/j.jneumeth.2006.11.017 pmid:17254636
    OpenUrlCrossRefPubMed
  62. ↵
    1. Peirce JW
    (2008) Generating stimuli for neuroscience using PsychoPy. Front Neuroinform 2:10. doi:10.3389/neuro.11.010.2008 pmid:19198666
    OpenUrlCrossRefPubMed
  63. ↵
    1. Phipson B,
    2. Smyth GK
    (2010) Permutation p-values should never be zero: calculating exact p-values when permutations are randomly drawn. Stat Appl Genet Mol Biol 9:Article39. doi:10.2202/1544-6115.1585 pmid:21044043
    OpenUrlCrossRefPubMed
  64. ↵
    1. Postle BR,
    2. Zarahn E,
    3. D'Esposito M
    (2000) Using event-related fMRI to assess delay-period activity during performance of spatial and nonspatial working memory tasks. Brain Res Brain Res Protoc 5:57–66. doi:10.1016/S1385-299X(99)00053-7 pmid:10719266
    OpenUrlCrossRefPubMed
  65. ↵
    1. Power JD,
    2. Barnes KA,
    3. Snyder AZ,
    4. Schlaggar BL,
    5. Petersen SE
    (2012) Spurious but systematic correlations in functional connectivity MRI networks arise from subject motion. Neuroimage 59:2142–2154. doi:10.1016/j.neuroimage.2011.10.018 pmid:22019881
    OpenUrlCrossRefPubMed
  66. ↵
    1. Pratte MS,
    2. Tong F
    (2014) Spatial specificity of working memory representations in early visual cortex. J Vision 14:1–12.
    OpenUrlAbstract/FREE Full Text
  67. ↵
    1. Ramnani N
    (2006) The primate cortico-cerebellar system: anatomy and function. Nat Rev Neurosci 7:511–522. doi:10.1038/nrn1953 pmid:16791141
    OpenUrlCrossRefPubMed
  68. ↵
    1. Riggall AC,
    2. Postle BR
    (2012) The relationship between working memory storage and elevated activity as measured with functional magnetic resonance imaging. J Neurosci 32:12990–12998. doi:10.1523/JNEUROSCI.1892-12.2012 pmid:22993416
    OpenUrlAbstract/FREE Full Text
  69. ↵
    1. Rissman J,
    2. Gazzaley A,
    3. D'Esposito M
    (2004) Measuring functional connectivity during distinct stages of a cognitive task. Neuroimage 23:752–763. doi:10.1016/j.neuroimage.2004.06.035 pmid:15488425
    OpenUrlCrossRefPubMed
  70. ↵
    1. Ron S,
    2. Robinson DA
    (1973) Eye movements evoked by cerebellar stimulation in the alert monkey. J Neurophysiol 36:1004–1022. doi:10.1152/jn.1973.36.6.1004 pmid:4202613
    OpenUrlCrossRefPubMed
  71. ↵
    1. Schmahmann JD,
    2. Pandya DN
    (1997) The cerebrocerebellar system. Int Rev Neurobiol 41:31–60. doi:10.1016/s0074-7742(08)60346-3 pmid:9378595
    OpenUrlCrossRefPubMed
  72. ↵
    1. Serences JT
    (2016) Neural mechanisms of information storage in visual short-term memory. Vision Res 128:53–67. doi:10.1016/j.visres.2016.09.010 pmid:27668990
    OpenUrlCrossRefPubMed
  73. ↵
    1. Serences JT,
    2. Ester EF,
    3. Vogel EK,
    4. Awh E
    (2009) Stimulus-specific delay activity in human primary visual cortex. Psychol Sci 20:207–214. doi:10.1111/j.1467-9280.2009.02276.x pmid:19170936
    OpenUrlCrossRefPubMed
  74. ↵
    1. Sheu YS,
    2. Liang Y,
    3. Desmond JE
    (2019) Disruption of cerebellar prediction in verbal working memory. Front Hum Neurosci 13:61. doi:10.3389/fnhum.2019.00061 pmid:30846935
    OpenUrlCrossRefPubMed
  75. ↵
    1. Soetedjo R,
    2. Fuchs AF
    (2006) Complex spike activity of Purkinje cells in the oculomotor vermis during behavioral adaptation of monkey saccades. J Neurosci 26:7741–7755. doi:10.1523/JNEUROSCI.4658-05.2006 pmid:16855102
    OpenUrlAbstract/FREE Full Text
  76. ↵
    1. Sprague TC,
    2. Ester EF,
    3. Serences JT
    (2014) Reconstructions of information in visual spatial working memory degrade with memory load. Curr Biol 24:2174–2180. doi:10.1016/j.cub.2014.07.066 pmid:25201683
    OpenUrlCrossRefPubMed
  77. ↵
    1. Sreenivasan KK,
    2. D'Esposito M
    (2019) The what, where and how of delay activity. Nat Rev Neurosci 20:466–481. doi:10.1038/s41583-019-0176-7 pmid:31086326
    OpenUrlCrossRefPubMed
  78. ↵
    1. Stoodley CJ,
    2. Schmahmann JD
    (2009) Functional topography in the human cerebellum: a meta-analysis of neuroimaging studies. Neuroimage 44:489–501. doi:10.1016/j.neuroimage.2008.08.039 pmid:18835452
    OpenUrlCrossRefPubMed
  79. ↵
    1. Stoodley CJ,
    2. Valera EM,
    3. Schmahmann JD
    (2012) Functional topography of the cerebellum for motor and cognitive tasks: an fMRI study. Neuroimage 59:1560–1570. doi:10.1016/j.neuroimage.2011.08.065 pmid:21907811
    OpenUrlCrossRefPubMed
  80. ↵
    1. Strick PL
    (1985) How do the basal ganglia and cerebellum gain access to the cortical motor areas? Behav Brain Res 18:107–123. doi:10.1016/0166-4328(85)90067-1 pmid:3913443
    OpenUrlCrossRefPubMed
  81. ↵
    1. Striemer CL,
    2. Cantelmi D,
    3. Cusimano MD,
    4. Danckert JA,
    5. Schweizer TA
    (2015) Deficits in reflexive covert attention following cerebellar injury. Front Hum Neurosci 9:428. doi:10.3389/fnhum.2015.00428 pmid:26300756
    OpenUrlCrossRefPubMed
  82. ↵
    1. Suchow JW,
    2. Brady TF,
    3. Fougnie D,
    4. Alvarez GA
    (2013) Modeling visual working memory with the MemToolbox. J Vision 13:1–8.
    OpenUrlAbstract/FREE Full Text
  83. ↵
    1. Takagi M,
    2. Zee DS,
    3. Tamargo RJ
    (1998) Effects of lesions of the oculomotor vermis on eye movements in primate: saccades. J Neurophysiol 80:1911–1931. doi:10.1152/jn.1998.80.4.1911 pmid:9772249
    OpenUrlCrossRefPubMed
  84. ↵
    1. Thier P,
    2. Dicke PW,
    3. Haas R,
    4. Thielert C-D,
    5. Catz N
    (2002) The role of the oculomotor vermis in the control of saccadic eye movements. Ann NY Acad Sci 978:50–62. doi:10.1111/j.1749-6632.2002.tb07555.x
    OpenUrlCrossRefPubMed
  85. ↵
    1. Todd N,
    2. Moeller S,
    3. Auerbach EJ,
    4. Yacoub E,
    5. Flandin G,
    6. Weiskopf N
    (2016) Evaluation of 2D multiband EPI imaging for high-resolution, whole-brain, task-based fMRI studies at 3T: Sensitivity and slice leakage artifacts. Neuroimage 124:32–42.
    OpenUrl
  86. ↵
    1. van Bergen RS,
    2. Ma WJ,
    3. Pratte MS,
    4. Jehee JFM
    (2015) Sensory uncertainty decoded from visual cortex predicts behavior. Nat Neurosci 18:1728–1730. doi:10.1038/nn.4150 pmid:26502262
    OpenUrlCrossRefPubMed
  87. ↵
    1. van den Berg R,
    2. Shin H,
    3. Chou W-C,
    4. George R,
    5. Ma WJ
    (2012) Variability in encoding precision accounts for visual short-term memory limitations. Proc Natl Acad Sci USA 109:8780–8785. doi:10.1073/pnas.1117465109 pmid:22582168
    OpenUrlAbstract/FREE Full Text
  88. ↵
    1. van den Berg R,
    2. Awh E,
    3. Ma WJ
    (2014) Factorial comparison of working memory models. Psychol Rev 121:124–149. doi:10.1037/a0035234 pmid:24490791
    OpenUrlCrossRefPubMed
  89. ↵
    1. van Es DM,
    2. van der Zwaag W,
    3. Knapen T
    (2019) Topographic maps of visual space in the human cerebellum. Curr Biol 29:1689–1694.e3. doi:10.1016/j.cub.2019.04.012 pmid:31080082
    OpenUrlCrossRefPubMed
  90. ↵
    1. Voogd J,
    2. Barmack NH
    (2006) Oculomotor cerebellum. Prog Brain Res 151:231–268. doi:10.1016/S0079-6123(05)51008-2 pmid:16221591
    OpenUrlCrossRefPubMed
  91. ↵
    1. Wager TD,
    2. Smith EE
    (2003) Neuroimaging studies of working memory. Cogn Affect Behav Neurosci 3:255–274. doi:10.3758/cabn.3.4.255 pmid:15040547
    OpenUrlCrossRefPubMed
  92. ↵
    1. Wang L,
    2. Mruczek REB,
    3. Arcaro MJ,
    4. Kastner S
    (2015) Probabilistic maps of visual topography in human cortex. Cereb Cortex 25:3911–3931.
    OpenUrlCrossRefPubMed
  93. ↵
    1. Wolpert DM,
    2. Miall RC,
    3. Kawato M
    (1998) Internal models in the cerebellum. Trends Cogn Sci 2:338–347. doi:10.1016/s1364-6613(98)01221-2 pmid:21227230
    OpenUrlCrossRefPubMed
  94. ↵
    1. Woo CW,
    2. Krishnan A,
    3. Wager TD
    (2014) Cluster-extent based thresholding in fMRI analyses: pitfalls and recommendations. Neuroimage 91:412–419. doi:10.1016/j.neuroimage.2013.12.058
    OpenUrlCrossRefPubMed
  95. ↵
    1. Xu J,
    2. Moeller S,
    3. Auerbach EJ,
    4. Strupp J,
    5. Smith SM,
    6. Feinberg DA,
    7. Yacoub E,
    8. Uğurbil K
    (2013) Evaluation of slice accelerations using multiband echo planar imaging at 3 T. Neuroimage 83:991–1001. doi:10.1016/j.neuroimage.2013.07.055 pmid:23899722
    OpenUrlCrossRefPubMed
  96. ↵
    1. Zhang W,
    2. Luck SJ
    (2008) Discrete fixed-resolution representations in visual working memory. Nature 453:233–235. doi:10.1038/nature06860 pmid:18385672
    OpenUrlCrossRefPubMed
Back to top

In this issue

The Journal of Neuroscience: 41 (5)
Journal of Neuroscience
Vol. 41, Issue 5
3 Feb 2021
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Index by author
  • Ed Board (PDF)
Email

Thank you for sharing this Journal of Neuroscience article.

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

Enter multiple addresses on separate lines or separate them with commas.
Stimulus-Specific Visual Working Memory Representations in Human Cerebellar Lobule VIIb/VIIIa
(Your Name) has forwarded a page to you from Journal of Neuroscience
(Your Name) thought you would be interested in this article in Journal of Neuroscience.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Print
View Full Page PDF
Citation Tools
Stimulus-Specific Visual Working Memory Representations in Human Cerebellar Lobule VIIb/VIIIa
James A. Brissenden, Sean M. Tobyne, Mark A. Halko, David C. Somers
Journal of Neuroscience 3 February 2021, 41 (5) 1033-1045; DOI: 10.1523/JNEUROSCI.1253-20.2020

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Respond to this article
Request Permissions
Share
Stimulus-Specific Visual Working Memory Representations in Human Cerebellar Lobule VIIb/VIIIa
James A. Brissenden, Sean M. Tobyne, Mark A. Halko, David C. Somers
Journal of Neuroscience 3 February 2021, 41 (5) 1033-1045; DOI: 10.1523/JNEUROSCI.1253-20.2020
Twitter logo Facebook logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

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

Keywords

  • cerebellum
  • cortico-cerebellar networks
  • delay activity
  • fMRI
  • stimulus selectivity
  • visual working memory

Responses to this article

Respond to this article

Jump to comment:

No eLetters have been published for this article.

Related Articles

Cited By...

More in this TOC Section

Research Articles

  • Regional Excitatory–Inhibitory Balance Relates to Self-Reference Effect on Recollection via the Precuneus/Posterior Cingulate Cortex–Medial Prefrontal Cortex Connectivity
  • Attention Alters Population Spatial Frequency Tuning
  • Sequence Diversity and Encoded Enzymatic Differences of Monocistronic L1 ORF2 mRNA Variants in the Aged Normal and Alzheimer's Disease Brain
Show more Research Articles

Behavioral/Cognitive

  • Attention Alters Population Spatial Frequency Tuning
  • Complex Impact of Stimulus Envelope on Motor Synchronization to Sound
  • The Molecular Substrates of Second-Order Conditioned Fear in the Basolateral Amygdala Complex
Show more Behavioral/Cognitive
  • Home
  • Alerts
  • Follow SFN on BlueSky
  • Visit Society for Neuroscience on Facebook
  • Follow Society for Neuroscience on Twitter
  • Follow Society for Neuroscience on LinkedIn
  • Visit Society for Neuroscience on Youtube
  • Follow our RSS feeds

Content

  • Early Release
  • Current Issue
  • Issue Archive
  • Collections

Information

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

About

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

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

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