Research Paper Volume 15, Issue 24 pp 14574—14590

Age-related alterations in the oscillatory dynamics serving verbal working memory processing

Seth D. Springer1,2, , Hannah J. Okelberry1,2, , Madelyn P. Willett1,2, , Hallie J. Johnson1,2, , Chloe E. Meehan1,3, , Mikki Schantell1,2, , Christine M. Embury1, , Maggie P. Rempe1,2, , Tony W. Wilson1,2,4, ,

  • 1 Institute for Human Neuroscience, Boys Town National Research Hospital, Boys Town, NE 68010, USA
  • 2 College of Medicine, University of Nebraska Medical Center, Omaha, NE 68198, USA
  • 3 Department of Neurology, Washington University School of Medicine, St. Louis, MO 63110, USA
  • 4 Department of Pharmacology and Neuroscience, Creighton University, Omaha, NE 68131, USA

Received: July 21, 2023       Accepted: November 15, 2023       Published: December 27, 2023      

https://doi.org/10.18632/aging.205403
How to Cite

Copyright: © 2023 Springer et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Abstract

Working memory (WM) is a foundational cognitive function involving the temporary storage of information. Unfortunately, WM is also one of the most sensitive cognitive functions to the detrimental effects of aging. Expanding the field’s understanding of age-related WM changes is critical to advancing the development of strategies to mitigate age-related WM declines. In the current study, we investigated the neural mechanisms serving WM function in seventy-eight healthy aging adults (range: 20.2–65.2 years) using magnetoencephalography (MEG) and a Sternberg WM task with letter stimuli. Neural activity during the different phases of the WM task (i.e., encoding, maintenance, and retrieval) were imaged using a time-frequency resolved beamformer and whole-brain statistics were performed. We found stronger increases in theta activity and stronger decreases in alpha and beta activity (i.e., more negative relative to baseline) as a function of healthy aging. Specifically, age-related increases in theta activity were detected during the encoding period in the primary visual and left prefrontal cortices. Additionally, alpha and beta oscillations were stronger (i.e., more negative) during both encoding and maintenance in the left prefrontal cortex in older individuals. Finally, alpha and beta oscillations during the retrieval phase were stronger (i.e., more negative) in older participants within the prefrontal, parietal, and temporal cortices. Together, these results indicate that healthy aging strongly modulates the neural oscillatory dynamics serving WM function.

Introduction

Working memory (WM) is a fundamental executive process which involves the storage, maintenance, and manipulation of relevant information. WM processes are generally grouped into three phases: information encoding, maintenance, and retrieval. The maintenance phase, which involves the rehearsal of the encoded information and the active inhibition of any new sensory stimuli, is the defining component which differentiates WM from other types of memory. Multiple studies have shown that WM processes are supported by a combination of neural spiking and oscillatory activity in the frontal and posterior cortices [18]. Historically, the neurophysiological underpinnings of WM function have been primarily investigated using functional MRI (fMRI). Such studies have revealed a left dominant, bilateral WM network with converging nodes in the dorsolateral prefrontal cortex (dlPFC), inferior frontal gyri, parieto-temporal cortices, and occipital regions [911]. More recently, electrophysiological studies have helped illuminate the specific oscillatory dynamics underlying WM function. Specifically, theta band activity during WM performance has been associated with the encoding of visual stimuli in primary visual cortices and higher-order cognitive control in the prefrontal cortex (PFC) [8, 1214]. Similarly, alpha and beta band activity has been associated with higher-order visual processing in the posterior cortices and cognitive control in the PFC [1417].

Previous literature has demonstrated that WM processes are vulnerable to age-related decline [1820], with these changes being particularly evident at higher cognitive loads [2123]. Across the adult lifespan, it has been shown that individuals tend to recruit the same brain regions when performing WM tasks [16, 24]. However, differential activation patterns become apparent when the cognitive load of the WM tasks is varied. For example, at lower cognitive loads, older individuals maintain WM performance comparable to younger individuals by recruiting more neural resources (i.e., stronger activations and/or more sites of activity) [2426]. Nevertheless, at higher cognitive loads, this neural compensation in older individuals faulters, resulting in behavioral decrements accompanied by weakened neural activity [21, 24, 2729]. This pattern of age-related overactivation, plateau, and subsequent underactivation is known as the compensation-related utilization of neural circuits hypothesis (CRUNCH) [30]. The CRUNCH phenomenon has been most commonly observed in PFC regions [21, 24, 27, 3133], which are thought to be responsible for exerting executive control over WM processes [10, 11]. Further, these executive prefrontal areas have shown age-related differentiation in the lateralization of neural activity. Specifically, younger individuals consistently show stronger left hemispheric neural activity during WM performance, compared to a more bilateral pattern of activity in older adults [16, 24, 29, 31, 34, 35].

In addition to the plethora of past fMRI and positron-emission tomography (PET) work investigating age-related WM changes, more recent electrophysiologic work has been conducted to quantify the neural dynamics of these WM changes with age. One such study by Proskovec and colleagues [16] found that older adults had stronger alpha/beta oscillations (i.e., decreases in power relative to baseline), in right frontal regions during WM encoding and maintenance, and in the right superior temporal gyrus later in the maintenance phase. Conversely, it was shown that younger adults had stronger alpha/beta oscillations in the right parieto-temporal regions early in the encoding phase [16]. Age-related changes in alpha activity during WM processing was later shown by Tran and associates [36], who demonstrated that older adults had lower pre-trial alpha phase consistency, and that this lower phase consistency predicted lower WM task accuracy in older participants. Finally, event-related potential (ERP) analyses have revealed that, overall, older adults have reduced amplitudes in the fronto-central positivity (i.e., P200) and parietal positivity (i.e., P300) [37]. However, when high- and low-performing older adults were separately analyzed, high-performing adults had stronger P200 and P300 amplitudes, more similar to younger adults [38]. This interesting pattern of ERP results are in agreement with CRUNCH, with the high-performing older adults seemingly showing compensatory hyperactivation and the low-performing adults having reached a “resource ceiling” and decompensated.

Despite previous work characterizing age-related changes in WM, no previous studies have examined spectrally-specific age-related changes in neural oscillatory activity throughout all phases of WM (i.e., encoding, maintenance, and retrieval). Thus, in the current study we used the spatiotemporal precision of MEG to examine the effects of healthy aging on the spectrally-resolved neural oscillatory dynamics underlying all three phases of WM processing. We hypothesized that older adults would require stronger engagement of key left hemispheric frontal and parieto-occipital WM hubs. Additionally, we expected that prefrontal activity lateralization (i.e., stronger left hemispheric activity) during WM performance would diminish as a function of age, with older individuals tending to utilize a more bilaterally distributed WM network.

Methods

Participants

Seventy-eight adults with a mean age of 45.10 (SD = 12.76) years were selected for inclusion in this study. The age range for males was 20.2 to 65.2 years and that for females was 21.4 to 62.2 years. These participants were chosen from a larger-scale study of accelerated aging in persons with HIV [39, 40], with only the HIV-negative participants included in this investigation of healthy aging. Of the 78 adults, 95% were right-handed, 79% were male, 10% were African-American, 6% were Asian, 77% were Caucasian, 4% were more than one race, and the remaining 3% preferred not to answer. This distribution corresponds closely to the racial demographics of the surrounding region. Notably, there was no effect of age on years of education (r74 = −.073, p = .531). Exclusionary criteria included any medical illness affecting CNS function (e.g., HIV/AIDS, Lupus, etc.), any neurological or psychiatric disorder, cognitive impairment, history of head trauma, current substance abuse, and the MEG laboratory’s standard exclusionary criteria (e.g., ferromagnetic implants). The Institutional Review Board reviewed and approved this investigation. Each participant provided written informed consent following a detailed description of the study.

Experimental paradigm

Participants were shown a centrally-presented fixation cross embedded in a 3 × 2 grid for 1.3 s (Figure 1). An array of six consonants then appeared at fixed locations within the grid for 2.0 s (i.e., encoding phase). Following the encoding phase, the letters disappeared and the empty grid remained on the screen for 3.0 s (i.e., maintenance phase). Finally, during the retrieval phase, a single probe consonant then appeared in the grid for 0.9 s, and the participant was instructed to respond via button press with their right index or middle finger as to whether the probe was in or out of the previous array of letters. A total of 128 trials were completed, equally split and pseudorandomized between in- and out-of-set trials, for a total run time of about 15 min. This task has been utilized in several previous studies from our laboratory [16, 17, 4144].

Modified Sternberg working memory paradigm. Participants were required to encode and maintain six letter stimuli. During the retrieval stage, participants were then required to indicate whether the probe letter was or was not present in the initial six letter set.

Figure 1. Modified Sternberg working memory paradigm. Participants were required to encode and maintain six letter stimuli. During the retrieval stage, participants were then required to indicate whether the probe letter was or was not present in the initial six letter set.

MEG data acquisition

All recordings were conducted in a magnetically-shielded room with active shielding engaged. Neuromagnetic responses were sampled continuously at 1 kHz, with an acquisition bandwidth of 0.1-330 Hz, using a MEGIN Vectorview MEG system with 306 magnetic sensors (MEGIN, Helsinki, Finland). During data acquisition, participants were monitored via real-time audio-visual feeds from inside the shielded room. Subject-wise MEG data were corrected for head motion and subjected to external noise reduction using signal space separation with a temporal extension [45].

Structural MRI processing and MEG coregistration

Preceding MEG measurement, four head position indicator (HPI) coils were attached to the participant’s head and localized, together with three fiducial points and at least 100 scalp surface points, with a 3D digitizer (Fastrak 3SF0002, Polhemus Navigator Sciences, Colchester, VT, USA). Once in the MEG, electrical currents with unique frequencies (e.g., 322 Hz) were fed into each of the HPI coils. These HPI coil currents induced measurable magnetic fields, allowing the position of the HPI coils to be activity tracked relative to the MEG sensors throughout the recording. Since HPI coil locations are known in head coordinates, all MEG measurements could be transformed into a common coordinate system. With this coordinate system, MEG data were coregistered with each individual’s high-resolution structural T1-weighted MRI data prior to source reconstruction using Brain Electrical Source Analysis MRI (BESA MRI, Version 2.0, BESA GmbH, Gräfelfing, Germany). Individual structural MRI data for each participant was acquired using a Siemens Prisma 3T scanner (Siemens Medical Solutions) with a 64-channel head coil. An MP-RAGE sequence was utilized with the following parameters: TR: 2300 ms; TE = 2.98 ms; flip angle = 9°; FOV = 256 mm; slice thickness = 1.00 mm; voxel size = 1 × 1 × 1 mm. Structural MRI data were transformed into standardized space (i.e., Talairach space) and aligned parallel to the anterior and posterior commissures. Following source analysis, each participant’s MEG functional images were also transformed into standardized space and spatially resampled.

MEG preprocessing, time-frequency transformation, and sensor-level statistics

Blink and cardiac artifacts were removed from the raw data using an adaptive artifact correction method in which brain activity is selectively separated from artifactual activities [46]. This adaptive artifact correction is accounted for during subsequent source analysis. The continuous magnetic time series was divided into 7200 ms epochs, with the baseline period being defined as the 400 ms prior to the encoding phase (i.e., −400 to 0 ms). Subsequently, epochs containing artifacts were removed based on a fixed threshold method, supplemented with visual inspection. Briefly, the amplitude and gradient distributions across all trials were determined per participant, and those trials containing the highest amplitude and/or gradient values relative to this distribution were rejected based on participant-specific thresholds. This approach was employed to minimize the impact of individual differences in sensor proximity to the brain and overall head size, which strongly affect MEG signal amplitude. Importantly, there were no age-related changes in either the amplitude (r74 = −.182, p = .115) or gradient (r74 = −.135, p = .245) thresholds used for artifact rejection. Artifact-free epochs were then transformed into the time-frequency domain using complex demodulation [4749], with a resolution of 1 Hz and 50 ms between 2 and 100 Hz. Following time-frequency transformation, spectral power estimates per sensor were averaged across trials to generate plots of mean spectral density per sensor. These sensor-level data were then normalized to the baseline power within each frequency bin, which was calculated as the mean power for that 1 Hz bin during the −400 to 0 ms time period.

The significant time-frequency windows used for source imaging were determined by statistical analysis of the sensor-level spectrograms across the entire array of gradiometers. Briefly, each pixel per spectrograms was initially evaluated using a mass univariate approach based on the general linear model, followed by cluster-based permutation testing to address the problem of multiple comparisons [50, 51]. Specifically, a two-stage procedure was utilized to minimize false positive results while maintaining sensitivity. The first stage consisted of performing paired-sample t-tests against baseline on each pixel per spectrogram and thresholding the output spectrograms of t-values at p < .05 to define time-frequency bins containing potentially significant oscillatory deviations from baseline. Bins that survived thresholding (at p < .05) were clustered with temporally and/or spectrally neighboring bins that also survived, and cluster values were derived by summing all t-values within each cluster. In stage two, nonparametric permutation testing was used to derive a distribution of cluster-values and the significance level of the cluster(s) from stage one was tested directly using this permuted distribution, which was the result of 10,000 permutations. Based on this cluster-based permutation analysis, only the time-frequency windows that contained significant oscillatory deviations from baseline at the p < .001, corrected, threshold across all participants were subjected to source imaging (i.e., beamforming).

MEG source imaging and statistics

Cortical networks were imaged through a time-frequency-resolved extension of the linearly constrained minimum variance (LCMV) beamformer [5254]. The images were derived from the cross spectral densities of all combinations of MEG gradiometers averaged over the time-frequency range of interest, and the solution of the forward problem for each location on a grid specified by input voxel space. In principle, the beamformer operator generates a spatial filter for each grid point that passes signals without attenuation from a given neural region, while suppressing activity in all other brain areas. The filter properties arise from the forward solution (i.e., lead field matrix) for each location on a volumetric grid specified by input voxel space, and from the MEG cross spectral density matrix. Basically, for each voxel, a set of beamformer weights is determined, which amounts to each MEG sensor being allocated a sensitivity weighting for activity in the particular voxel. Following convention, the source power in these images was normalized per participant using a pre-stimulus period (i.e., baseline) of equal duration and bandwidth [55]. Such images are typically referred to as pseudo-t maps, with units (pseudo-t) that reflect noise-normalized power differences (i.e., active vs. passive) per voxel. MEG pre-processing and imaging used BESA (version 7.0) software.

After imaging, average whole-brain maps were computed across all participants for the selected time-frequency windows. These 3D maps of brain activity were used to assess the neuroanatomical basis of the significant oscillatory responses identified through the sensor-level analysis. Finally, these source images were subjected to correlation analyses with age to investigate how the oscillatory processes underlying WM encoding, maintenance, and retrieval change as a function of age. To control the multiple comparisons inherent to whole-brain voxel-wise statistics, cluster-based permutation testing was again utilized [50, 51, 56]. This whole-brain cluster-based permutation testing was virtually identical to the approach used for sensor-level statistics (see above), with the exception that the input data was voxel-level neural responses and we used a different statistical threshold following permutation testing (i.e., stage one threshold of p < .001; threshold of p < .05 following permutation testing for multiple comparisons correction).

Statistical analyses

Statistical analyses were performed using custom R [57] and MATLAB [58] scripts. Correlation analyses were used to test differences as a function of healthy aging in education, accuracy, reaction time, total correct trials, amplitude cutoffs, and gradient cutoffs. Permutation testing of the sensor-level spectrograms and whole-brain images was performed using BESA Statistics (Version 2.1, BESA GmbH, Gräfelfing, Germany). Finally, to reduce the impact of outliers on statistical analyses, participants with values 2.5 SDs above or below the group mean were excluded for each analysis.

Data availability

The data that support the findings of this study are available from the corresponding author, upon reasonable request. Specifically, the anonymized data will be available for non-commercial research purposes and responses will occur at earliest convenience with a goal of delivering the data within 1 month.

Results

Behavioral results

Of the 78 participants that performed the verbal WM task, two were removed due to poor task performance (i.e., <65% correct). While there was no effect of age on accuracy (r74 = .012, p = .919), there was an effect of age on reaction time (r74 = .453, p < .001) such that older participants were slower than younger participants. Across all participants, the mean accuracy was 84.33% (SD = 7.04%) and the mean reaction time was 895.4 ms (SD = 193.06 ms).

Spectrally-specific neural oscillations

Sensor-level time-frequency analysis across all participants revealed significant clusters (p < .001, corrected) of theta, alpha, and beta band oscillatory activity during the encoding, maintenance, and retrieval phases (Figure 2). During the encoding phase, theta band (3-6 Hz) activity sharply increased immediately following stimulus presentation (i.e., 0 ms), continued through the encoding phase, and slowly dissipated over the first 800 ms of the maintenance phase (i.e., 3-6 Hz, 0-2800 ms; p < .001, corrected). Additionally, decreases in alpha band (9-15 Hz) activity began around 200 ms after the onset of the encoding grid and were sustained throughout the encoding phase, terminating near the beginning of the maintenance phase (i.e., 9-15 Hz, 200-1800 ms; p < .001, corrected). During the maintenance phase, there was a significant cluster of increased alpha/beta activity in a slightly higher frequency band (12-17 Hz; 2500-2900 ms; p < .001, corrected). During the retrieval phase, an increase in theta activity followed retrieval probe presentation (3-6 Hz; 5000-5350 ms; p < .001, corrected). Finally, overlapping decreases in alpha (8-12 Hz; 5350-5750 ms) and beta (14-19 Hz; 5250-5650 ms) band activity were observed during the retrieval phase (both ps < .001, corrected).

Grand-averaged time-frequency spectrogram from a sensor near parieto-occipital areas (MEG1922) with time (ms) shown on the x-axis and frequency (Hz) denoted on the y-axis. A color scale bar is shown to the right of the spectrogram which shows the percent power change relative to the baseline period (−400 to 0 ms).

Figure 2. Grand-averaged time-frequency spectrogram from a sensor near parieto-occipital areas (MEG1922) with time (ms) shown on the x-axis and frequency (Hz) denoted on the y-axis. A color scale bar is shown to the right of the spectrogram which shows the percent power change relative to the baseline period (−400 to 0 ms).

Beamformer analysis

In order to determine the neural regions involved in verbal WM, the aforementioned sensor-level time-frequency bins of interest for theta, alpha, and beta band oscillatory activity were imaged using a time-frequency-resolved beamformer. Considering that beamforming windows must be the same length as the baseline period (i.e., 400 ms), the oscillatory deviations from baseline were imaged in 400 ms increments or less. For the temporally extended oscillatory responses (i.e., those over 400 ms), participant-level image averaging was performed for each neural response as follows: theta encoding (3-6 Hz; 0-2000 ms), theta maintenance (3-6 Hz; 2000-2800 ms), and alpha encoding (9-15 Hz; 200-1800 ms). These whole-brain average maps per participant were then averaged across all individuals to allow for visualization of the neural oscillations during the different WM phases (i.e., encoding, maintenance, and retrieval). These images revealed differential patterns of neural activity in the theta, alpha, and beta bands throughout all phases of WM performance (Figure 3). During the encoding phase, a strong increase in theta activity relative to the baseline was observed in the bilateral primary visual cortices, while decreased alpha activity (i.e., more negative relative to baseline) was observed in the lateral occipital cortices bilaterally, and left frontal cortex. During the maintenance phase, increased theta activity was observed in the right dorsolateral prefrontal cortex (dlPFC). Further, the decreased alpha/beta oscillations (i.e., more negative) observed in the frontal cortices during the encoding phase were sustained through the maintenance phase, along with decreased alpha activity in the lateral occipital cortices bilaterally. During the retrieval phase, theta band activity was shown to strongly increase in the right posterior cingulate gyrus and the bilateral primary visual cortices. Decreased alpha activity (i.e., more negative) during the retrieval phase were observed in the lateral occipital cortices, left supramarginal gyrus (SMG), and left angular gyrus. Similarly, decreased beta oscillations (i.e., more negative) during the retrieval phase were found in the lateral occipital cortices, bilateral primary motor cortices, and bilateral superior parietal cortices.

Grand-averaged beamformer images (pseudo-t) across all participants for theta (left) and alpha/beta (right) bands. During the encoding phase (top), there were increases in theta band amplitude relative to baseline in the primary visual cortex and decreases in alpha band activity (i.e., more negative relative to baseline) in the lateral occipital cortices, parietal cortices, and left hemisphere frontal cortex (including language areas). During the maintenance phase (middle), theta increases were limited to the right dlPFC and decreased alpha/beta activity in the left hemisphere frontal cortex (relative to baseline) was largely sustained. Additionally, a decrease in alpha/band activity relative to the baseline was observed in the lateral occipital cortices during the maintenance phase. During the retrieval phase (bottom), theta band activity increased in the right posterior cingulate gyrus and bilateral primary visual cortices, while alpha activity decreased (i.e., became more negative) in the lateral occipital cortices, left SMG, and left angular gyrus. Finally, decreased beta activity (i.e., more negative relative to baseline) in the bilateral primary motor cortices, lateral occipital cortices, and right superior parietal lobe was observed.

Figure 3. Grand-averaged beamformer images (pseudo-t) across all participants for theta (left) and alpha/beta (right) bands. During the encoding phase (top), there were increases in theta band amplitude relative to baseline in the primary visual cortex and decreases in alpha band activity (i.e., more negative relative to baseline) in the lateral occipital cortices, parietal cortices, and left hemisphere frontal cortex (including language areas). During the maintenance phase (middle), theta increases were limited to the right dlPFC and decreased alpha/beta activity in the left hemisphere frontal cortex (relative to baseline) was largely sustained. Additionally, a decrease in alpha/band activity relative to the baseline was observed in the lateral occipital cortices during the maintenance phase. During the retrieval phase (bottom), theta band activity increased in the right posterior cingulate gyrus and bilateral primary visual cortices, while alpha activity decreased (i.e., became more negative) in the lateral occipital cortices, left SMG, and left angular gyrus. Finally, decreased beta activity (i.e., more negative relative to baseline) in the bilateral primary motor cortices, lateral occipital cortices, and right superior parietal lobe was observed.

To statistically examine the effects of healthy aging on the neural dynamics serving WM function, we subjected the frequency-specific whole-brain encoding, maintenance, and retrieval participant-level average maps to correlation analysis. These whole-brain correlation analyses were then followed by cluster-based permutation testing to correct for multiple comparisons. These analyses revealed that theta oscillations during encoding became significantly stronger as a function of age in the left primary visual cortex (p = .037, corrected) and left dlPFC (p = .045, corrected; Figure 4), and that there were no age-related theta effects during the maintenance or retrieval phases. In contrast, the decreases in alpha and beta became stronger (i.e., more negative) with increasing age throughout all phases of WM processing (i.e., encoding, maintenance, and retrieval; Figure 5). During the encoding phase, decreases in alpha became stronger (i.e., more negative) as a function of age in the left middle frontal gyrus (p = .023, corrected) and the left anterior cingulate gyrus (p = .033, corrected). This age-related change in activity in the left middle frontal gyrus was sustained through the maintenance phase (p = .029, corrected), along with stronger decreases in alpha/band responses (i.e., more negative) in the left dlPFC (p = .031, corrected). During the retrieval phase, decreases in alpha oscillations also became stronger (i.e., more negative) as a function of age in the left insula (p = .009, corrected), left anterior cingulate (p = .020, corrected), and left middle temporal gyrus (p = .031, corrected). Finally, the decreases in beta oscillations (14-19 Hz) during the retrieval phase were found to get stronger (i.e., more negative) with increasing age in the left primary motor cortex (p < .001, corrected), right primary somatosensory cortex (p < .001, corrected), bilateral superior parietal cortices (ps < .001, corrected), left dlPFC (p = .007, corrected), and left anterior superior temporal gyrus (p = .020, corrected). Alpha and beta retrieval phase scatterplots that were not included in Figure 5 can be found in the Supplementary Material (Supplementary Figure 1). A complete summary of all significant effects can be found in Table 1.

Effect of age on theta oscillations during the encoding phase. Whole-brain linear regression analysis revealed that theta activity increased in the left primary visual cortex and left dlPFC with older age. Linear regression plots of peak voxel pseudo-t values are shown as a function of age. Lines of best-fit and 95% CI (shaded area) are overlaid.

Figure 4. Effect of age on theta oscillations during the encoding phase. Whole-brain linear regression analysis revealed that theta activity increased in the left primary visual cortex and left dlPFC with older age. Linear regression plots of peak voxel pseudo-t values are shown as a function of age. Lines of best-fit and 95% CI (shaded area) are overlaid.

Effect of age on alpha and beta oscillations during the encoding, maintenance, and retrieval phases. Whole-brain linear regression analyses testing for age-related changes in alpha (encoding and retrieval), alpha/beta (maintenance), and beta (retrieval) are shown during the encoding (top), maintenance (middle), and retrieval bottom) phases. During encoding (top), decreases in alpha activity became significantly stronger (i.e., more negative relative to baseline), as a function of age, in prefrontal cortical regions and anterior cingulate. Additionally, age-related alterations in alpha and beta activity during the maintenance period (middle) were found in the left prefrontal cortices, while age-related effects during retrieval (bottom) were observed across a diverse set of brain regions that have been implicated in working memory processing (i.e., frontal, temporal, and parietal cortices). Beyond statistical maps, linear regression plots of peak voxel pseudo-t values are shown as a function of age for encoding and maintenance peaks. Lines of best-fit and 95% CI (shaded area) are overlaid.

Figure 5. Effect of age on alpha and beta oscillations during the encoding, maintenance, and retrieval phases. Whole-brain linear regression analyses testing for age-related changes in alpha (encoding and retrieval), alpha/beta (maintenance), and beta (retrieval) are shown during the encoding (top), maintenance (middle), and retrieval bottom) phases. During encoding (top), decreases in alpha activity became significantly stronger (i.e., more negative relative to baseline), as a function of age, in prefrontal cortical regions and anterior cingulate. Additionally, age-related alterations in alpha and beta activity during the maintenance period (middle) were found in the left prefrontal cortices, while age-related effects during retrieval (bottom) were observed across a diverse set of brain regions that have been implicated in working memory processing (i.e., frontal, temporal, and parietal cortices). Beyond statistical maps, linear regression plots of peak voxel pseudo-t values are shown as a function of age for encoding and maintenance peaks. Lines of best-fit and 95% CI (shaded area) are overlaid.

Table 1. Coordinates of the peak response in each significant age correlation cluster.

Region/Effect of interestFrequencyXYZdfStatistic
Encodingr
Left primary visual cortexθ−2−851070.42
Left dlPFCθ−34243870.40
Left middle frontal gyrusα−42−55474−.41
Left anterior cingulate gyrusα−10−54274−.40
Maintenancer
Left middle frontal gyrusα/β−42−15472−.42
Left dlPFCα/β−42284272−.42
Retrievalr
Left insulaα−30−131072−.41
Left anterior cingulate gyrusα−6−54272−.40
Left middle temporal gyrusα−58−37-372−.39
Left primary motor cortexβ−42−135071−.45
Right primary somatosensory cortexβ42−214271−.46
Left superior parietal cortexβ−14−455471−.45
Right superior parietal cortexβ30−535871−.45
Left dlPFCβ−50322271−.41
Left anterior superior temporal gyrusβ−5024-1171−.41
All test statistics are significant at p < .001 level and survive cluster-based permutation testing. All coordinates are in Talairach space. Abbreviations: dlPFC: dorsolateral prefrontal cortex; df: degrees of freedom.

Discussion

In the present study, we investigated the effect of healthy aging on the neural oscillatory dynamics serving verbal WM processing using MEG and advanced source reconstruction. Similar to previous studies utilizing a modified Sternberg WM task in healthy aging adults, we found that reaction time increased with age and that there were no significant aging effects on task accuracy [16, 59]. Across all participants, our results also showed that WM processing is supported by regionally distinct neural oscillations in the theta, alpha, and beta bands throughout all phases (i.e., encoding, maintenance, and retrieval) and that such activity is broadly affected by healthy aging. Briefly, theta oscillations were shown to become stronger in the primary visual cortices and in the frontal cortex with increasing age. Furthermore, throughout all WM phases and in several frontal, parietal, and temporal regions, alpha, beta, and alpha/beta activity became stronger (i.e., more negative) as a function of age. Critically, these results represent the first oscillatory analysis of verbal WM function in a healthy aging sample in which neural activity was examined throughout all phases of processing, including the retrieval stage. Below, we discuss the implications of these novel findings on our understanding of how healthy aging affects verbal WM processing.

Across all frequencies, our findings demonstrate widespread and robust cortical activation patterns during verbal working memory task performance, as well as age-related changes in these neural activations.

These task-related activations and age-related changes are most prevalent in the occipital, parietal, and frontal cortices; regions which have been shown to be critical to verbal working memory performance [9, 10]. According to the theories of working memory function originally proposed by Baddeley and colleagues, the neural subsystems serving working memory processing include the visuospatial sketchpad, phonological loop, episodic buffer, and the central executive [60, 61]. The visuospatial sketchpad is responsible for processing the visual and spatial components of a stimulus, and is served by the occipital and parieto-temporal cortices [9, 10]. The phonological loop deals with the processing and rehearsal of auditory and language information. Further, the phonological loop is a predominately left hemispheric network comprised of the phonological store in the parietal lobe [9, 62], which briefly stores fleeting memory traces, and the articulatory process in the inferior frontal gyrus and superior temporal gyrus (e.g., Broca’s area) [6365], responsible for the rehearsal and manipulation of information in the phonological store (e.g., through subvocal repetition). The episodic buffer is a more distributed subsystem which links working memory and long-term memory, facilitating the integration of information from separate modalities (i.e., visual, spatial, auditory, language, etc.) and binding information into coherent episodes [61]. Finally, the central executive component exerts top-down control over the other working memory subsystems, via the allocation of attentional resources and direction of information flow to the other subsystems, and is thought be housed in the prefrontal cortex [10, 11]. Considering the letter-based stimuli utilized in our modified Sternberg working memory task, we expected robust activation of each of the aforementioned subsystems.

One of our most interesting findings was the broad strengthening of alpha and beta oscillations (i.e., decreases from baseline) with increasing age. The local suppression of alpha and beta activity, relative to baseline levels, has been shown to represent active cortical engagement (i.e., active visual processing in higher-order visual areas) [16, 17, 6668]. This is further supported by multimodal studies which have demonstrated co-localization of decreases in alpha and beta activity with increased BOLD fMRI signal during cognitive tasks [69, 70]. During the encoding phase, older individuals showed significantly stronger decreases in alpha activity (i.e., more negative relative to baseline) in the left middle frontal gyrus (MFG) and the left anterior cingulate. This significant age-related shift in left MFG activity extended into the maintenance period, with the addition of age-related effects in alpha/beta activity in the left dlPFC. These prefrontal cortical regions are thought to act as central executive nodes during WM task performance, allowing for the allocation of WM subsystem resources and the direction of attention [10, 11]. These findings of age-related increases in frontal regions are in widespread agreement with the aging literature, which has demonstrated that increased prefrontal neural recruitment can offset declines in cognitive performance as a function of healthy aging (i.e., CRUNCH) [21, 24, 27, 3033]. Of note, we found no significant age-related changes during the encoding and maintenance phases in alpha/beta activity in either visual processing or language areas, which have historically been linked to the visuospatial sketchpad and phonological loop, respectively [9, 10, 62]. In sum, the present work is in agreement with our hypotheses that older adults would require stronger engagement of left hemispheric frontal WM hubs and reinforces previous aging WM literature in demonstrating how crucial increased neural activity in prefrontal cortices is for older individuals during WM processing.

Beyond the encoding and maintenance phases, we found significant alpha and beta oscillatory changes as a function of age during the retrieval phase. Though previous studies have extensively characterized neural activity during WM encoding and maintenance, studies analyzing the neural dynamics during the retrieval phase are sparse. During the retrieval phase, older participants showed stronger decreases in alpha activity (i.e., more negative) in the left posterior middle temporal gyrus (MTG), left insula, and left anterior cingulate. Additionally, beta activity during the retrieval phase decreased (i.e., more negative) with increasing age in the left primary motor cortex, bilateral superior parietal cortices, left dlPFC, left anterior superior temporal gyrus, and other regions. Thus, unlike the age-related alpha/beta effects observed during the encoding and maintenance phases, which were largely confined to prefrontal cortices, age-related changes in alpha/beta activity during the retrieval phase were found to extend to other WM subsystems. Specifically, in the retrieval phase, older individuals tended to more strongly engage neural regions linked to higher-order visual processing (i.e., visuospatial sketchpad; left posterior MTG), language processing (i.e., left STG), phonological storage (i.e., superior parietal lobes), central executive processes (i.e., left dlPFC), and motor execution (i.e., primary motor cortex). These retrieval phase findings are in line with our hypothesis that older adults would more strongly recruit left hemispheric frontal and parieto-occipital WM nodes, relative to younger participants. As noted above, only a few studies have examined the retrieval stage, Guran and associates [71] found that older (i.e., 53-80 years) compared to younger (i.e., 18-30 years) individuals had stronger alpha/beta decreases relative to baseline in sensors overlying frontal, central, and parietal regions [71], which is consistent with our findings. Finally, our finding of stronger decreases in beta (i.e., more negative) in the primary motor cortex in older participants is consistent with work from the motor control literature [72, 73]. Thus, this effect of aging on motor beta activity is likely task independent.

We also found robust age-related effects in the theta range. Theta oscillations in the primary visual cortex have been widely implicated in basic visual processing, while frontal theta activity has been shown to be important for the higher-order organization and maintenance of stimulus information [8, 1214, 66, 74]. Though theta activity has been shown to be critical for WM function, age-related changes in these theta band dynamics have only been weakly characterized. Age-related theta alterations in the current study were found to be confined to the encoding phase, with activity in the primary visual cortex and left dlPFC getting stronger as a function of increasing age. One comparable study found that older participants (i.e., 62–71 years) were more susceptible to WM disruption and that frontal midline theta was stronger in older individuals (i.e., 19–29 years) [75]. Our findings of age-related increases in theta activity, particularly those in the dlPFC, are consistent with our alpha and beta band findings and the neural compensation hypothesis (i.e., CRUNCH). Additionally, theta band activity in the PFC of healthy young adults during WM performance has been shown to increase stepwise with increasing task load [12, 14, 66], further supporting the idea that the increased theta activity with increasing age may represent these older individuals needing to exert more neural resources to complete this WM task (i.e., less efficient processing). Altogether, in conjunction with our alpha and beta findings, we have demonstrated that the widely cited age-related increases in neural recruitment are not frequency specific, with age-related effects spanning multiple frequency bands (i.e., theta, alpha, and beta) and working memory subsystems.

Interestingly, our results did not show the hypothesized and commonly reported age-related change in the lateralization of prefrontal activation. Specifically, previous work has demonstrated that younger individuals more strongly utilize left hemispheric prefrontal regions during working memory performance, compared to older individuals who exhibit more bilateral prefrontal activity [16, 24, 29, 31, 34, 35]. This discrepancy may be due to study design and analysis differences. Firstly, our sample of 78 health adults is nearly double the size of most previous healthy aging WM studies and our whole-brain data were modeled using age as a continuous variable, compared to dichotomizing age groups as has commonly been done [16, 24, 31, 35]. Secondly, the current study utilized cluster-based permutation testing of the source images to strongly control for false positives, compared to the much less stringent multiple comparison correction approaches that have previously been used (e.g., uncorrected or cluster-size corrections) [16, 31, 35]. Finally, the oldest participants in the current study were younger (i.e., 65 years) than many aging WM studies (e.g., 75–82 years) [24, 31, 34, 35]. Thus, our relatively younger participants may not have reached a level of cognitive dysfunction that necessitates bilateral prefrontal recruitment, but future studies are needed to confirm these findings.

Before closing, it is important to note the limitations of this study. First, as mentioned above, the oldest participants in our sample were slightly younger (i.e., 65 years) than some previous WM studies of aging (e.g., 75–82 years). For future studies, increasing the age of the oldest participants would allow for the hypothesized transition from age-related compensation (i.e., increased neural activity and normal accuracy) to age-related decompensation (i.e., decreased neural activity and decreased accuracy) to be more thoroughly studied. Finally, the modified Sternberg WM paradigm that was utilized for the current study required that participants remember six letter stimuli per trial and previous work has shown that varying the cognitive load of WM tasks leads to an interesting pattern of effects in older individuals (i.e., greater neural activity at low working memory loads with decompensation and decreased neural activity at higher working memory loads) [21, 24, 2729]. Thus, future aging WM research would benefit from having multiple memory load conditions of varying difficulty.

Taken together, these data support the previous literature showing greater recruitment of the prefrontal cortex in older individuals during WM performance. Importantly, our findings extend this literature by demonstrating that age-related hyperactivation of prefrontal regions occurs in all phases of WM (i.e., encoding, maintenance, and retrieval) and is multispectral involving theta, alpha, and beta oscillatory activity. Further, though the age-related changes in alpha and beta band activity were restricted to prefrontal executive control regions during the encoding and maintenance phases, these age-related changes emerged in a more widespread network that included language (i.e., phonological loop) and motor regions during the retrieval phase. Developing a better understanding of these age-related alterations in WM function may have important implications for developing targeted interventions aimed at improving cognitive function and promoting healthy aging across the lifespan.

Supplementary Materials

Supplementary Figure 1

Author Contributions

Seth D. Springer: Data curation, Formal analysis, Writing – Original draft, Writing – Review and Editing; Hannah J. Okelberry: Data curation, Project administration, Participant accrual; Madelyn P. Willett: Project administration, Participant accrual; Hallie J. Johnson: Project administration, Participant accrual; Chloe E. Meehan: Data curation; Mikki Schantell: Data curation; Christine M. Embury: Formal analysis; Maggie P. Rempe: Formal analysis; Tony W. Wilson: Study conceptualization, Writing – Review and Editing, Supervision, Project administration, Funding acquisition.

Conflicts of Interest

The authors declare no conflicts of interest related to this study.

Ethical Statement and Consent

The authors assert that all procedures contributing to this work comply with the ethical standards of the relevant national and institutional committees on human experimentation and with the Helsinki Declaration of 1975, as revised in 2008. The Institutional Review Board reviewed and approved this investigation (IRB protocol number 20-29xp). Each participant provided written informed consent following detailed description of the study.

Funding

This study was supported by the National Institutes of Health through grants R01-MH116782 (TWW), R01-MH118013 (TWW), R01-DA047828 (TWW), R01-DA056223 (TWW), P20-GM144641 (TWW), F30-AG076259 (SDS), and F31-DA056296 (MS). The funders had no role in the study design, collection, analysis, or interpretation of data, nor did they influence writing the report or the decision to submit this work for publication.

References

  • 1. Miller EK, Lundqvist M, Bastos AM. Working Memory 2.0. Neuron. 2018; 100:463–75. https://doi.org/10.1016/j.neuron.2018.09.023 [PubMed]
  • 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–49. https://doi.org/10.1152/jn.1989.61.2.331 [PubMed]
  • 3. Fuster JM, Alexander GE. Neuron activity related to short-term memory. Science. 1971; 173:652–4. https://doi.org/10.1126/science.173.3997.652 [PubMed]
  • 4. Salazar RF, Dotson NM, Bressler SL, Gray CM. Content-specific fronto-parietal synchronization during visual working memory. Science. 2012; 338:1097–100. https://doi.org/10.1126/science.1224000 [PubMed]
  • 5. Lundqvist M, Rose J, Herman P, Brincat SL, Buschman TJ, Miller EK. Gamma and Beta Bursts Underlie Working Memory. Neuron. 2016; 90:152–64. https://doi.org/10.1016/j.neuron.2016.02.028 [PubMed]
  • 6. Honkanen R, Rouhinen S, Wang SH, Palva JM, Palva S. Gamma Oscillations Underlie the Maintenance of Feature-Specific Information and the Contents of Visual Working Memory. Cereb Cortex. 2015; 25:3788–801. https://doi.org/10.1093/cercor/bhu263 [PubMed]
  • 7. Roux F, Wibral M, Mohr HM, Singer W, Uhlhaas PJ. Gamma-band activity in human prefrontal cortex codes for the number of relevant items maintained in working memory. J Neurosci. 2012; 32:12411–20. https://doi.org/10.1523/JNEUROSCI.0421-12.2012 [PubMed]
  • 8. van Ede F, Jensen O, Maris E. Supramodal Theta, Gamma, and Sustained Fields Predict Modality-specific Modulations of Alpha and Beta Oscillations during Visual and Tactile Working Memory. J Cogn Neurosci. 2017; 29:1455–72. https://doi.org/10.1162/jocn_a_01129 [PubMed]
  • 9. Rottschy C, Langner R, Dogan I, Reetz K, Laird AR, Schulz JB, Fox PT, Eickhoff SB. Modelling neural correlates of working memory: a coordinate-based meta-analysis. Neuroimage. 2012; 60:830–46. https://doi.org/10.1016/j.neuroimage.2011.11.050 [PubMed]
  • 10. Baddeley A. Working memory: looking back and looking forward. Nat Rev Neurosci. 2003; 4:829–39. https://doi.org/10.1038/nrn1201 [PubMed]
  • 11. Riley MR, Constantinidis C. Role of Prefrontal Persistent Activity in Working Memory. Front Syst Neurosci. 2016; 9:181. https://doi.org/10.3389/fnsys.2015.00181 [PubMed]
  • 12. Jensen O, Tesche CD. Frontal theta activity in humans increases with memory load in a working memory task. Eur J Neurosci. 2002; 15:1395–9. https://doi.org/10.1046/j.1460-9568.2002.01975.x [PubMed]
  • 13. Brookes MJ, Wood JR, Stevenson CM, Zumer JM, White TP, Liddle PF, Morris PG. Changes in brain network activity during working memory tasks: a magnetoencephalography study. Neuroimage. 2011; 55:1804–15. https://doi.org/10.1016/j.neuroimage.2010.10.074 [PubMed]
  • 14. Proskovec AL, Heinrichs-Graham E, Wilson TW. Load modulates the alpha and beta oscillatory dynamics serving verbal working memory. Neuroimage. 2019; 184:256–65. https://doi.org/10.1016/j.neuroimage.2018.09.022 [PubMed]
  • 15. 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 U S A. 2018; 115:1117–22. https://doi.org/10.1073/pnas.1710323115 [PubMed]
  • 16. Proskovec AL, Heinrichs-Graham E, Wilson TW. Aging modulates the oscillatory dynamics underlying successful working memory encoding and maintenance. Hum Brain Mapp. 2016; 37:2348–61. https://doi.org/10.1002/hbm.23178 [PubMed]
  • 17. Heinrichs-Graham E, Wilson TW. Spatiotemporal oscillatory dynamics during the encoding and maintenance phases of a visual working memory task. Cortex. 2015; 69:121–30. https://doi.org/10.1016/j.cortex.2015.04.022 [PubMed]
  • 18. Park DC, Lautenschlager G, Hedden T, Davidson NS, Smith AD, Smith PK. Models of visuospatial and verbal memory across the adult life span. Psychol Aging. 2002; 17:299–320. [PubMed]
  • 19. Bopp KL, Verhaeghen P. Aging and verbal memory span: a meta-analysis. J Gerontol B Psychol Sci Soc Sci. 2005; 60:P223–33. https://doi.org/10.1093/geronb/60.5.p223 [PubMed]
  • 20. Salthouse TA, Babcock RL, Shaw RJ. Effects of adult age on structural and operational capacities in working memory. Psychol Aging. 1991; 6:118–27. https://doi.org/10.1037/0882-7974.6.1.118 [PubMed]
  • 21. Mattay VS, Fera F, Tessitore A, Hariri AR, Berman KF, Das S, Meyer-Lindenberg A, Goldberg TE, Callicott JH, Weinberger DR. Neurophysiological correlates of age-related changes in working memory capacity. Neurosci Lett. 2006; 392:32–7. https://doi.org/10.1016/j.neulet.2005.09.025 [PubMed]
  • 22. Heinzel S, Lorenz RC, Duong QL, Rapp MA, Deserno L. Prefrontal-parietal effective connectivity during working memory in older adults. Neurobiol Aging. 2017; 57:18–27. https://doi.org/10.1016/j.neurobiolaging.2017.05.005 [PubMed]
  • 23. Billig AR, Feng NC, Behforuzi H, McFeeley BM, Nicastri CM, Daffner KR. Capacity-limited resources are used for managing sensory degradation and cognitive demands: Implications for age-related cognitive decline and dementia. Cortex. 2020; 133:277–94. https://doi.org/10.1016/j.cortex.2020.09.005 [PubMed]
  • 24. Cappell KA, Gmeindl L, Reuter-Lorenz PA. Age differences in prefontal recruitment during verbal working memory maintenance depend on memory load. Cortex. 2010; 46:462–73. https://doi.org/10.1016/j.cortex.2009.11.009 [PubMed]
  • 25. Reuter-Lorenz PA, Park DC. Human neuroscience and the aging mind: a new look at old problems. J Gerontol B Psychol Sci Soc Sci. 2010; 65:405–15. https://doi.org/10.1093/geronb/gbq035 [PubMed]
  • 26. Schneider-Garces NJ, Gordon BA, Brumback-Peltz CR, Shin E, Lee Y, Sutton BP, Maclin EL, Gratton G, Fabiani M. Span, CRUNCH, and beyond: working memory capacity and the aging brain. J Cogn Neurosci. 2010; 22:655–69. https://doi.org/10.1162/jocn.2009.21230 [PubMed]
  • 27. Smith EE, Geva A, Jonides J, Miller A, Reuter-Lorenz P, Koeppe RA. The neural basis of task-switching in working memory: effects of performance and aging. Proc Natl Acad Sci U S A. 2001; 98:2095–100. https://doi.org/10.1073/pnas.98.4.2095 [PubMed]
  • 28. Rypma B, Prabhakaran V, Desmond JE, Gabrieli JD. Age differences in prefrontal cortical activity in working memory. Psychol Aging. 2001; 16:371–84. https://doi.org/10.1037/0882-7974.16.3.371 [PubMed]
  • 29. Rypma B, Eldreth DA, Rebbechi D. Age-related differences in activation-performance relations in delayed-response tasks: a multiple component analysis. Cortex. 2007; 43:65–76. https://doi.org/10.1016/s0010-9452(08)70446-5 [PubMed]
  • 30. Reuter-Lorenz PA, Cappell KA. Neurocognitive aging and the compensation hypothesis. Curr Dir Psychol Sci. 2008; 17:177–82. https://doi.org/10.1111/j.1467-8721.2008.00570.x
  • 31. Emery L, Heaven TJ, Paxton JL, Braver TS. Age-related changes in neural activity during performance matched working memory manipulation. Neuroimage. 2008; 42:1577–86. https://doi.org/10.1016/j.neuroimage.2008.06.021 [PubMed]
  • 32. Morcom AM, Henson RNA. Increased Prefrontal Activity with Aging Reflects Nonspecific Neural Responses Rather than Compensation. J Neurosci. 2018; 38:7303–13. https://doi.org/10.1523/JNEUROSCI.1701-17.2018 [PubMed]
  • 33. Fakhri M, Sikaroodi H, Maleki F, Ali Oghabian M, Ghanaati H. Age-related frontal hyperactivation observed across different working memory tasks: an fMRI study. Behav Neurol. 2012; 25:351–61. https://doi.org/10.3233/BEN-2012-120280 [PubMed]
  • 34. Cabeza R, Daselaar SM, Dolcos F, Prince SE, Budde M, Nyberg L. Task-independent and task-specific age effects on brain activity during working memory, visual attention and episodic retrieval. Cereb Cortex. 2004; 14:364–75. https://doi.org/10.1093/cercor/bhg133 [PubMed]
  • 35. Reuter-Lorenz PA, Jonides J, Smith EE, Hartley A, Miller A, Marshuetz C, Koeppe RA. Age differences in the frontal lateralization of verbal and spatial working memory revealed by PET. J Cogn Neurosci. 2000; 12:174–87. https://doi.org/10.1162/089892900561814 [PubMed]
  • 36. Tran TT, Hoffner NC, LaHue SC, Tseng L, Voytek B. Alpha phase dynamics predict age-related visual working memory decline. Neuroimage. 2016; 143:196–203. https://doi.org/10.1016/j.neuroimage.2016.08.052 [PubMed]
  • 37. Lubitz AF, Niedeggen M, Feser M. Aging and working memory performance: Electrophysiological correlates of high and low performing elderly. Neuropsychologia. 2017; 106:42–51. https://doi.org/10.1016/j.neuropsychologia.2017.09.002 [PubMed]
  • 38. Morrison C, Kamal F, Taler V. The influence of working memory performance on event-related potentials in young and older adults. Cogn Neurosci. 2019; 10:117–28. https://doi.org/10.1080/17588928.2019.1570104 [PubMed]
  • 39. Spooner RK, Taylor BK, Moshfegh CM, Ahmad IM, Dyball KN, Emanuel K, Schlichte SL, Schantell M, May PE, O'Neill J, Kubat M, Bares SH, Swindells S, et al. Neuroinflammatory profiles regulated by the redox environment predicted cognitive dysfunction in people living with HIV: A cross-sectional study. EBioMedicine. 2021; 70:103487. https://doi.org/10.1016/j.ebiom.2021.103487 [PubMed]
  • 40. Spooner RK, Taylor BK, Ahmad IM, Dyball K, Emanuel K, O'Neill J, Kubat M, Swindells S, Fox HS, Bares SH, Stauch KL, Zimmerman MC, Wilson TW. Mitochondrial redox environments predict sensorimotor brain-behavior dynamics in adults with HIV. Brain Behav Immun. 2023; 107:265–75. https://doi.org/10.1016/j.bbi.2022.10.004 [PubMed]
  • 41. Wiesman AI, Heinrichs-Graham E, McDermott TJ, Santamaria PM, Gendelman HE, Wilson TW. Quiet connections: Reduced fronto-temporal connectivity in nondemented Parkinson's Disease during working memory encoding. Hum Brain Mapp. 2016; 37:3224–35. https://doi.org/10.1002/hbm.23237 [PubMed]
  • 42. Wilson TW, Proskovec AL, Heinrichs-Graham E, O'Neill J, Robertson KR, Fox HS, Swindells S. Aberrant Neuronal Dynamics during Working Memory Operations in the Aging HIV-Infected Brain. Sci Rep. 2017; 7:41568. https://doi.org/10.1038/srep41568 [PubMed]
  • 43. Koshy SM, Wiesman AI, Spooner RK, Embury C, Rezich MT, Heinrichs-Graham E, Wilson TW. Multielectrode Transcranial Electrical Stimulation of the Left and Right Prefrontal Cortices Differentially Impacts Verbal Working Memory Neural Circuitry. Cereb Cortex. 2020; 30:2389–400. https://doi.org/10.1093/cercor/bhz246 [PubMed]
  • 44. Embury CM, Wiesman AI, Proskovec AL, Heinrichs-Graham E, McDermott TJ, Lord GH, Brau KL, Drincic AT, Desouza CV, Wilson TW. Altered Brain Dynamics in Patients With Type 1 Diabetes During Working Memory Processing. Diabetes. 2018; 67:1140–8. https://doi.org/10.2337/db17-1382 [PubMed]
  • 45. Taulu S, Simola J. Spatiotemporal signal space separation method for rejecting nearby interference in MEG measurements. Phys Med Biol. 2006; 51:1759–68. https://doi.org/10.1088/0031-9155/51/7/008 [PubMed]
  • 46. Ille N, Berg P, Scherg M. Artifact correction of the ongoing EEG using spatial filters based on artifact and brain signal topographies. J Clin Neurophysiol. 2002; 19:113–24. https://doi.org/10.1097/00004691-200203000-00002 [PubMed]
  • 47. Kovach CK, Gander PE. The demodulated band transform. J Neurosci Methods. 2016; 261:135–54. https://doi.org/10.1016/j.jneumeth.2015.12.004 [PubMed]
  • 48. Hoechstetter K, Bornfleth H, Weckesser D, Ille N, Berg P, Scherg M. BESA source coherence: a new method to study cortical oscillatory coupling. Brain Topogr. 2004; 16:233–8. https://doi.org/10.1023/b:brat.0000032857.55223.5d [PubMed]
  • 49. Papp N, Ktonas P. Critical evaluation of complex demodulation techniques for the quantification of bioelectrical activity. Biomed Sci Instrum. 1977; 13:135–45. [PubMed]
  • 50. Maris E, Oostenveld R. Nonparametric statistical testing of EEG- and MEG-data. J Neurosci Methods. 2007; 164:177–90. https://doi.org/10.1016/j.jneumeth.2007.03.024 [PubMed]
  • 51. Ernst MD. Permutation Methods: A Basis for Exact Inference. Statist Sci. 2004; 19: 676–85. https://doi.org/10.1214/088342304000000396
  • 52. Dalal SS, Sekihara K, Nagarajan SS. Modified beamformers for coherent source region suppression. IEEE Trans Biomed Eng. 2006; 53:1357–63. https://doi.org/10.1109/TBME.2006.873752 [PubMed]
  • 53. Gross J, Kujala J, Hamalainen M, Timmermann L, Schnitzler A, Salmelin R. Dynamic imaging of coherent sources: Studying neural interactions in the human brain. Proc Natl Acad Sci U S A. 2001; 98:694–9. https://doi.org/10.1073/pnas.98.2.694 [PubMed]
  • 54. Van Veen BD, van Drongelen W, Yuchtman M, Suzuki A. Localization of brain electrical activity via linearly constrained minimum variance spatial filtering. IEEE Trans Biomed Eng. 1997; 44:867–80. https://doi.org/10.1109/10.623056 [PubMed]
  • 55. Hillebrand A, Singh KD, Holliday IE, Furlong PL, Barnes GR. A new approach to neuroimaging with magnetoencephalography. Hum Brain Mapp. 2005; 25:199–211. https://doi.org/10.1002/hbm.20102 [PubMed]
  • 56. Nichols TE, Holmes AP. Nonparametric permutation tests for functional neuroimaging: a primer with examples. Hum Brain Mapp. 2002; 15:1–25. https://doi.org/10.1002/hbm.1058 [PubMed]
  • 57. R Core Team. R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing; 2020. https://www.R-project.org/.
  • 58. MATLAB version: 9.5.0 (R2018b), Natick. Massachusetts: The MathWorks Inc.; 2018.
  • 59. Sghirripa S, Graetz L, Merkin A, Rogasch NC, Semmler JG, Goldsworthy MR. Load-dependent modulation of alpha oscillations during working memory encoding and retention in young and older adults. Psychophysiology. 2021; 58:e13719. https://doi.org/10.1111/psyp.13719 [PubMed]
  • 60. Baddeley AD, Hitch G. Working Memory. In: Bower GH, ed. Psychology of Learning and Motivation. Academic Press. 1974; 47–89.
  • 61. Baddeley A. The episodic buffer: a new component of working memory? Trends Cogn Sci. 2000; 4:417–23. https://doi.org/10.1016/s1364-6613(00)01538-2 [PubMed]
  • 62. Vallar G, Di Betta AM, Silveri MC. The phonological short-term store-rehearsal system: patterns of impairment and neural correlates. Neuropsychologia. 1997; 35:795–812. https://doi.org/10.1016/s0028-3932(96)00127-3 [PubMed]
  • 63. Friederici AD, Rüschemeyer SA, Hahne A, Fiebach CJ. The role of left inferior frontal and superior temporal cortex in sentence comprehension: localizing syntactic and semantic processes. Cereb Cortex. 2003; 13:170–7. https://doi.org/10.1093/cercor/13.2.170 [PubMed]
  • 64. Awh E, Jonides J, Smith EE, Schumacher EH, Koeppe RA, Katz S. Dissociation of Storage and Rehearsal in Verbal Working Memory: Evidence From Positron Emission Tomography. Psychological Science. 1996; 7:25–31. https://doi.org/10.1111/j.1467-9280.1996.tb00662.x
  • 65. Paulesu E, Frith CD, Frackowiak RS. The neural correlates of the verbal component of working memory. Nature. 1993; 362:342–5. https://doi.org/10.1038/362342a0 [PubMed]
  • 66. Pavlov YG, Kotchoubey B. Oscillatory brain activity and maintenance of verbal and visual working memory: A systematic review. Psychophysiology. 2022; 59:e13735. https://doi.org/10.1111/psyp.13735 [PubMed]
  • 67. Neuper C, Pfurtscheller G. Event-related dynamics of cortical rhythms: frequency-specific features and functional correlates. Int J Psychophysiol. 2001; 43:41–58. https://doi.org/10.1016/s0167-8760(01)00178-7 [PubMed]
  • 68. Hanslmayr S, Staudigl T, Fellner MC. Oscillatory power decreases and long-term memory: the information via desynchronization hypothesis. Front Hum Neurosci. 2012; 6:74. https://doi.org/10.3389/fnhum.2012.00074 [PubMed]
  • 69. Scheeringa R, Fries P, Petersson KM, Oostenveld R, Grothe I, Norris DG, Hagoort P, Bastiaansen MC. Neuronal dynamics underlying high- and low-frequency EEG oscillations contribute independently to the human BOLD signal. Neuron. 2011; 69:572–83. https://doi.org/10.1016/j.neuron.2010.11.044 [PubMed]
  • 70. Hermes D, Nguyen M, Winawer J. Neuronal synchrony and the relation between the blood-oxygen-level dependent response and the local field potential. PLoS Biol. 2017; 15:e2001461. https://doi.org/10.1371/journal.pbio.2001461 [PubMed]
  • 71. Guran CA, Herweg NA, Bunzeck N. Age-Related Decreases in the Retrieval Practice Effect Directly Relate to Changes in Alpha-Beta Oscillations. J Neurosci. 2019; 39:4344–52. https://doi.org/10.1523/JNEUROSCI.2791-18.2019 [PubMed]
  • 72. Heinrichs-Graham E, McDermott TJ, Mills MS, Wiesman AI, Wang YP, Stephen JM, Calhoun VD, Wilson TW. The lifespan trajectory of neural oscillatory activity in the motor system. Dev Cogn Neurosci. 2018; 30:159–68. https://doi.org/10.1016/j.dcn.2018.02.013 [PubMed]
  • 73. Heinrichs-Graham E, Wilson TW. Is an absolute level of cortical beta suppression required for proper movement? Magnetoencephalographic evidence from healthy aging. Neuroimage. 2016; 134:514–21. https://doi.org/10.1016/j.neuroimage.2016.04.032 [PubMed]
  • 74. Onton J, Delorme A, Makeig S. Frontal midline EEG dynamics during working memory. Neuroimage. 2005; 27:341–56. https://doi.org/10.1016/j.neuroimage.2005.04.014 [PubMed]
  • 75. Arnau S, Wascher E, Küper K. Age-related differences in reallocating cognitive resources when dealing with interruptions. Neuroimage. 2019; 191:292–302. https://doi.org/10.1016/j.neuroimage.2019.02.048 [PubMed]