Cerebral cortex activation during the Sternberg verbal working memory task

Ilya S. Bakulin; Zabirova Ah; Kopnin Pn; Sinitsyn Do; Poydasheva Ag; Fedorov Mv; Gnedovskaya Ev; Suponeva Na; Piradov Ma

CEREBRAL CORTEX ACTIVATION DURING THE STERNBERG VERBAL WORKING MEMORY TASK

Bakulin IS1 Zabirova AH1, Kopnin PN1, Sinitsyn DO1, Poydasheva AG1, Fedorov MV2, Gnedovskaya EV1,2, Suponeva NA1, Piradov MA1

1 Research Center of Neurology, Moscow, Russia

2 Skolkovo Institute of Science and Technology, Moscow, Russia

Despite intensive study, the data regarding functional role of specific brain regions in the working memory processes still remain controversial. The study was aimed to determine the activation of cerebral cortex regions at different stages of the working memory task (information encoding, maintenance and retrieval). Functional magnetic resonance imaging (fMRI) with the modified Sternberg task was applied to 19 healthy volunteers. The objective of the task was to memorize and retain in memory the sequence of 7 letters with the subsequent comparison of one letter with the sequence. Activation was analyzed during three periods of the task compared to the rest period, as well as temporal dynamics of changes in BOLD signal intensity in three regions: left dorsolateral prefrontal, left posterior parietal and left occipital cortex. According to the results, significant activation of the regions in prefrontal and posterior parietal cortex was observed during all periods of the task р < 0.05), but there were changes in its localization and lateralization. The activation pattern during the maintenance period corresponded to the fronto-parietal control network components. According to the analysis of temporal dynamics of changes in BOLD signal intensity, the most prominent activation of the dorsolateral prefrontal cortex and parietal cortex was observed in the end of the encoding period, during the maintenance period and in the beginning of the retrieval period, which confirmed the role of those areas in the working memory processes. The maximum of occipital cortex activation was observed during encoding period. The study confirmed the functional role of the dorsolateral prefrontal cortex and posterior parietal cortex in the working memory mechanisms during all stages of the Sternberg task.

Keywords: working memory, functional magnetic resonance imaging, the Sternberg task, delayed-response task, prefrontal cortex, fronto-parietal control network

Author contribution: Bakulin IS, Fedorov MV, Gnedovskaya EV, Suponeva NA, Piradov MA — study planning and design; Bakulin IS, Zabirova AH — literature analysis; Bakulin IS, Poydasheva AG — data acquisition; Zabirova AH, Kopnin PN, Sinitsyn DO — data analysis; Bakulin IS, Zabirova AH, Sinitsyn DO, Fedorov MV Gnedovskaya EV, Suponeva NA, Piradov MA — data interpretation; Bakulin IS, Zabirova AH — manuscript draft preparation; Kopnin PN, Sinitsyn DO, Poydasheva AG, Fedorov MV Gnedovskaya EV, Suponeva NA, Piradov MA — manuscript draft editing; all the authors — final version of the article preparation.

Compliance with ethical standards: the study was approved by the Ethics Committee of Research Center of Neurology (protocol № 9-3/17 dated August 30, 2017). All volunteers submitted the informed consent to participation in the study.

[23 Correspondence should be addressed: Ilya S. Bakulin

Volokolamskoye shosse, 80, Moscow, 125367; [email protected]

Received: 22.01.2020 Accepted: 22.02.2020 Published online: 29.02.2020

DOI: 10.24075/brsmu.2020.013

АКТИВАЦИЯ КОРЫ ГОЛОВНОГО МОЗГА ПРИ ВЫПОЛНЕНИИ ЗАДАЧИ СТЕРНБЕРГА НА ВЕРБАЛЬНУЮ РАБОЧУЮ ПАМЯТЬ

И. С. Бакулин1 А. Х. Забирова1, П. Н. Копнин1, Д. О. Синицын1, А. Г. Пойдашева1, М. В. Федоров2, Е. В. Гнедовская1,2, Н. А. Супонева1, М. А. Пирадов1

1 Научный центр неврологии, Москва, Россия

2 Сколковский институт науки и технологий, Москва, Россия

Несмотря на интенсивное изучение, данные о функциональном значении отдельных областей головного мозга для процессов рабочей памяти (РП) остаются противоречивыми. Целью данного исследования было определить активацию участков коры на разных этапах выполнения задачи на РП (запечатление, удержание и извлечение информации). На 19 здоровых добровольцах была использована функциональная магнитная резонансная томография (фМРТ) с модифицированной задачей Стернберга. Суть задачи — в запоминании и удержании в памяти последовательности из семи букв с последующим сравнением с ней одной буквы. Проанализирована активация в трех периодах задачи по сравнению с периодом покоя, а также динамика изменения интенсивности BOLD-сигнала в зависимости от времени в трех областях: левой дорсолатеральной префронтальной коре (ДЛПФК), левой задней теменной коре и левой затылочной коре. По результатам анализа, на протяжении всех периодов задачи наблюдали статистически значимую активацию участков ПФК и задней теменной коры (р < 0,05), но с изменением их локализации и латерализации. При этом паттерн активации в периоде удержания соответствовал компонентам лобно-теменной сети контроля. По результатам анализа изменения интенсивности BOLD-сигнала наиболее выраженную активацию ДЛПФК и теменной коры наблюдали в конце периода запечатления, периоде удержания и начале периода извлечения, что подтверждает роль данных областей в процессах РП. Пик активации затылочной коры соответствовал периоду запечатления. Данное исследование подтверждает функциональное значение ДЛПФК и задней теменной коры для процессов РП на всех этапах выполнения задачи Стернберга.

Ключевые слова: рабочая память, функциональная магнитно-резонансная томография, задача Стернберга, задача с отсроченным воспроизведением, префронтальная кора, лобно-теменная сеть контроля

Вклад авторов: И. С. Бакулин, М. В. Федоров, Е. В. Гнедовская, Н. А. Супонева, М. А. Пирадов — планирование и дизайн исследования; И. С. Бакулин, А. Х. Забирова — анализ литературы; И. С. Бакулин, А. Г. Пойдашева — сбор данных; А. Х. Забирова, П. Н. Копнин, Д. О. Синицын — анализ данных; И. С. Бакулин, А. Х. Забирова, Д. О. Синицын, М. В. Федоров, Е. В. Гнедовская, Н. А. Супонева, М. А. Пирадов — интерпретация данных; И. С. Бакулин, А. Х. Забирова — подготовка черновика рукописи статьи; П. Н. Копнин, Д. О. Синицын, А. Г. Пойдашева, М. В. Федоров, Е. В. Гнедовская, Н. А. Супонева, М. А. Пирадов — внесение правок в черновик рукописи статьи; все авторы — подготовка финального варианта статьи.

Соблюдение этических стандартов: исследование одобрено этическим комитетом ФГБНУ «Научный центр неврологии» (протокол № 9-3/17 от 30 августа 2017 г.). Все участники исследования подписали добровольное информированное согласие на участие в исследовании.

123 Для корреспонденции: Илья Сергеевич Бакулин

Волоколамское шоссе, д. 80, г. Москва, 125367; [email protected]

Статья получена: 22.01.2020 Статья принята к печати: 22.02.2020 Опубликована онлайн: 29.02.2020 DOI: 10.24075/vrgmu.2020.013

Working memory is the ability to maintain and process a limited amount of information for a short time for its use in the ongoing cognitive activity. The working memory decline is associated with normal aging [1], it is also observed in patients with various neurological and psychiatric disorders [2], that confirms the importance of the working memory study.

Research data regarding working memory mechanisms remain controversial [3]. One of the brain areas, most investigated in context of working memory, is the dorsolateral prefrontal cortex. The persistent neural activity during the maintenance period of the working memory task was firstly observed in the prefrontal cortex of primates [4]. Later, the similar activity was detected in other regions. However, a number of studies confirmed the importance of the persistent neural activity in prefrontal cortex. For example, despite the simultaneous activation of prefrontal and temporal cortex during the task, inactivation of prefrontal cortex leaded to the inhibition of persistent neural activity in the temporal cortex of primates, while inactivation of the temporal cortex did not affect the prefrontal cortex [5]. The temporary inactivation of the prefrontal cortex led to a more pronounced working memory decline than inactivation of the posterior parietal cortex [6].

The role of the dorsolateral prefrontal cortex in the working memory in humans was confirmed by fMRI with different tasks (n-back, delayed match-to-sample task, Sternberg task, etc.) [7-9], as well as by transcranial magnetic stimulation (TMS) [10]. However, many cerebral cortex regions are activated during a working memory task: parietal lobe (areas around the intraparietal sulcus, superior parietal lobule), insula, anterior cingulate gyrus [7, 11, 12], visual cortex [13]. At the same time, the question of the link between activation and working memory mechanisms (information maintenance and manipulation) remains open, since activation can be caused by additional cognitive processes involved in the working memory task performance, for example, the perception of a stimulus [13]. The use of the delayed response tasks (the Sternberg task in particular) [14] allows to separate some working memory processes by analyzing different task stages: encoding, maintenance and retrieval.

The assessment of temporal dynamics of activation can be helpful in comparison of fMRI results with the studies using causal methods (for example, TMS) where the online effects at the different stages of the task execution were studied. Given the increasing interest in the use of non-invasive brain stimulation techniques for the improvement of cognitive functions [15, 16] the issue is also important for clinical neurology. fMRI results can help in determination of the TMS target which might increase the stimulation effectiveness [17]. The data of the cortical activity dynamics might also be used for combination of TMS protocols with cognitive load.

Despite the frequent use of the Sternberg task in neuroimaging studies of verbal working memory [11], the analysis of activation at different time was performed in a relatively small number of studies [9, 18-20]. Considering the limitations caused by the fMRI temporal resolution, the reproducibility of the results of such studies, which is necessary to confirm the possibility of using fMRI for the investigation of individual working memory mechanisms, still remains important.

The study was aimed to determine the activation of cerebral cortex regions at the different stages of the modified Sternberg task for verbal working memory, and to assess the dynamics of individual regions activation by the analysis of temporal dynamics of BOLD signal changes. The following cerebral cortex regions were chosen for assessment (according to their potential functional role in the Sternberg task performance,

including participation in the working memory processes): left dorsolateral prefrontal cortex [9, 10], left posterior parietal cortex [7, 11, 21] and left occipital cortex [13].

METHODS

Study participants

The study was performed in 2017-2019 at the Research Center of Neurology.

Inclusion criteria: age 18-55; dominant right hand (according to the Edinburgh Handedness Inventory [22]); normal (corrected to normal) vision. Exclusion criteria: contraindications for MRI [23]; neurological or psychiatric disorders; intake of drugs that affect the central nervous system at the time of the study; severe chronic somatic disorders.

In accordance with the inclusion criteria, 19 healthy volunteers (6 men) aged 22-31 were recruited, all of them right-handed. The study was performed as a part of the project "Effects of high-frequency repetitive transcranial magnetic stimulation over the left dorsolateral prefrontal cortex on working memory in healthy volunteers".

fMRI with the Sternberg task

fMRI with the modified Sternberg task for verbal working memory was used in the study [9, 14]. The participants received instructions about all stages of the investigation before performing the task in the MRI scanner. During the investigation, participants were wearing headphones to reduce the scanner noise. The participants could see the task screen located in the scanning room through a mirror, which was installed on the 8-channel head coil. Hand-held grips with buttons were given to each participant, which they pressed according to the instructions.

Task-fMRI used in our study had a block design and consisted of 24 task presentations (Fig. 1). Each presentation included the activation period (the Sternberg task) and the rest period. The activation period was divided into three blocks (periods): encoding, maintenance and retrieval. During the encoding period (1st-3rd seconds from the start of the task) the stimulus (7 latin consonants selected in random order using the algorithm created in the MATLAB R2017a (Mathworks, Natick, MA; USA)) was shown on the screen. Angular size of the letters was 1.5. During the maintenance period (4-12th second) the participant maintained the stimuli in memory, while the screen remained blank (without fixation point). During the retrieval period (13-15th second) a single latin letter appeared on the screen. The participant was instructed to press a button in the right hand-held grip if the letter presented during the retrieval period matched to the sequence, and a button in the left hand-held grip if the letter had not. Matching and non-matching stimuli were presented in random order, with 12 presentations of each stimulus type. The duration of activation period, as well as of the rest period without any stimulus (black screen), was 15 seconds. Participants were instructed to keep their eyes open during all periods. The experiment duration was 12 minutes.

Image acquisition and preprocessing

The study was performed using a 3 T scanner Magnetom Verio (Siemens; Germany). Anatomical data were acquired using the 3D-T1-gradient echo (T1-MPR) MRI sequence and consisted of 176 sagittal slices with the following parameters:

TR 1900 ms, TE 2.47 ms, slice thickness 1.0 mm, distance factor 50%; voxel size 1.0 ■ 1.0 ■ 1.0 mm3, field-of-view 250 mm. The functional data were acquired by investigation of the BOLD signal changes using T2*-gradient echo MRI sequence with the following parameters: TR 3000 ms, TE 30 ms, slice thickness 3 mm, distance factor 25%, voxel size 3.0 ■ 3.0 ■ 3.0 mm3, field-of-view 192 mm, 36 axial slices. To ensure dynamic magnetic equilibrium, first five images were excluded from the analysis. Preprocessing and statistical analysis were performed using the SPM12 software (Institute of Neurology, London; UK) based on the MATLAB R2017a (Mathworks, Natick, MA; USA).

Image preprocessing included following stages: 1) image realignment for correction of motion artifacts using the shifting and rotation transformations; 2) co-registration of anatomical and functional data; 3) segmentation of the obtained structural images and their spatial normalization to transform the individual images into the Montreal Neurological Institute (MNI) space; 4) spatial smoothing with a Gaussian kernel 8 mm FWMH to reduce the impact of anatomical variability.

The low-frequency drift was removed using the high-pass filter (time constant 128 s).

Statistical analysis

Statistical parametric mapping was performed using general linear models, estimated separately for the signal of each voxel. During the 1st level analysis, individual activation data were obtained for the encoding, maintenance and retrieval periods. The last 9 seconds of the rest period (22-30th seconds) were used as the control state for the activation level calculation. Head motion parameters were included as nuisance regressors. The data of 19 participants were included in the group analysis (separately for each period). Activation was assessed both at peak and cluster levels. The threshold for the inclusion of an individual voxel in the clusters was set at p < 0.0001. The clusters in the gray matter of cerebral hemispheres and cerebellum were selected for the analysis with a significance level of p < 0.05 (with the family-wise error, FWE correction for multiple comparisons ).

For regions showing significant activation during the maintenance period, the additional comparison with activation of the frontoparietal control network (FPCN) was performed, according to the resting state fMRI data [24]. The results were visualized using xjView [25] through overlay of the FPCN from the atlas [24] with the activation map of the maintenance period.

The temporal dynamics of the changes in BOLD signal intensity was analyzed in three regions of interest: left

dorsolateral prefrontal cortex, left posterior parietal cortex and left occipital cortex. To determine the specific regions that were active during the task performance , the difference between the signal intensity in the maintenance period and the rest period was used as contrast within the selected regions of interest. Voxels with activation above the threshold were selected in the regions of interest (p < 0.0001). Activation clusters containing selected voxel and restricted to spheres with the centers in a voxel with locally maximal activation were identified on group data. Radii were selected manually to obtain the region of intersection of the sphere and cluster with a sufficient number of voxels without spreading to adjacent anatomical regions (7 mm for dorsolateral prefrontal cortex, 8 mm for posterior parietal cortex, 10 mm for occipital cortex). At the individual level, BOLD signal was extracted from voxels in the determined cluster, as well as a signal adjusted for alignment effects and temporal non-sphericity in the framework of the autoregressive process AR(1). The extracted and simulated signals were averaged over the cluster voxels, and then over all participants. After that, the signal was averaged over the 24 task presentations and interpolated with periodic boundary conditions by the fast Fourier transformation method with a step of TR/10 (equal to 0.3 s). By the interpretation of the data, 5 s were subtracted from the onsets of BOLD events to account for the hemodynamic response delay.

RESULTS

The activation of the cerebral cortex regions at different stages of the modified Sternberg task is shown in Fig. 2.

During the encoding period significant activation was detected in the clusters within the right and left parietal lobes (supramarginal and postcentral gyri, superior parietal lobule), and also in the cerebellar cortex and occipital lobe (cuneus) bilaterally. Activation was observed in the right middle frontal and superior frontal gyri, and in the lower part of the left middle frontal gyrus spreading to the inferior frontal gyrus. A cluster of activation was also detected in the left rolandic operculum (Table 1, Fig. 2A).

During the maintenance period significant activation was found in the medial parts of the frontal lobes (bilaterally), prefrontal cortex (right middle frontal gyrus , left middle frontal and inferior frontal gyri), parietal lobe (supramarginal gyrus and superior parietal lobule bilaterally, right angular gyrus). In addition, activation was detected in the anterior parts of the insula (bilaterally) and in the left prefrontal cortex (middle frontal gyrus) (Table 2, Fig. 2B).

Encoding

24 presentations

Retrieval

Duration — 30 s

Fig. 1. fMRI with the modified Sternberg task for verbal memory (delayed item recognition task)

Fig. 2. Graphical depiction of significant activation regions on the cerebral cortex convexital surface during three periods of task compared with rest state (clusters containing more than 100 voxels are shown). А — encoding period; B — maintenance period; C — retrieval period. Left to right: top view, right view, left view. Color bars display the T-statistics level

Comparing the activation pattern in maintenance period with the FPCN, a partial overlap of the activation zones in the parietal and frontal lobes was noted, as well as in the medial parts of the frontal lobes and insular cortex (bilaterally) (Fig. 3).

During the retrieval period the regions of significant activation were detected in the right medial orbital and anterior orbital gyri, right middle temporal gyrus, in the cerebellar hemispheres bilaterally (spreading to the fusiform gyrus on the left). A cluster of activation was also detected in the right anterior insula and in the operculum. Within the parietal lobe clusters of activation were found in the right supramarginal gyrus, as well as in the left postcentral and supramarginal gyri. A cluster of activation was also identified in the right middle frontal gyrus. Activation was detected in the cerebellar vermis spreading to the lingual gyri bilaterally (Table 3, Fig. 2C).

A cluster of activation in the left middle frontal gyrus at the intersection with a sphere with a radius of 7 mm and center at the point of maximal activation [-36; 47; 5] (according to MNI space) containing 21 voxels was selected as the region of interest, in which temporal change of BOLD signal intensity was analyzed. According to the analysis results, an increase in signal intensity in the left middle frontal gyrus begins at the beginning of the task. The activation curve contains two peaks, 6 and 15 seconds from the beginning of the block. After the latter the signal intensity gradually decreases and reaches a plateau (Fig. 4).

The temporal change of BOLD signal intensity was analyzed for a cluster in the left inferior parietal lobule at the intersection with a sphere with a radius of 8 mm and the center at the point of maximum activation [-42; -49; 54] (according to MNI space)

Region Cluster level Peak level Peak coordinates, MNI [x; y; z]

Volume, voxels рFWE Т рFWE

Supramarginal gyrus, D 3587 <0.001 14.56 <0.001 [50; -32; 46]

Cerebellar hemisphere, D 8293 <0.001 12.55 0.001 [28; -58; -46]

Supramarginal gyrus, S 1502 <0.001 11.34 0.002 [-42 -34; 42]

Superior parietal lobule, S 11.28 0.002 [-28 -60; 46]

Postcentral gyrus, S 9.43 0.019 [-44 -36; 60]

Middle frontal gyrus, D 935 <0.001 10.49 0.005 [36 34; 26]

Cuneus, bilateral 569 <0.001 9.99 0.009 [2; -84; 34]

Cerebellar hemisphere, S 75 0.008 7.55 0.131 [-36; -46; -46]

Middle frontal gyrus, S Inferior frontal gyrus (pars triangularis), S 362 <0.001 8.58 0.038 [2; -22; 16]

Rolandic operculum, S 128 0.001 7.45 0.046 [34; -32; 6]

Middle and superior frontal gyri, D 48 0.037 6.15 0.516 [28; 8; 54]

Cerebellar hemisphere, S 152 <0.001 5.97 0.597 [-30; -62; -46]

Table 1. Regions of significant activation during the encoding period of the modified Sternberg task for verbal working memory

Note: MNI — Montreal Neurological Institute, FWE — family-wise error, T — T-statistics, D (dexter) — right, S (sinister) — left.

Table 2. Regions of significant activation during the maintenance period of the modified Sternberg task for verbal working memory

Region Cluster level Peak level Peak coordinates, MNI [x; y; z]

Volume, voxels рFWE Т рFWE

Superial frontal gyrus medial parts, bilateral 420 <0.001 14.12 <0.001 [0; 26; 44]

Supplementary motor area, bilateral 10.37 0.006 [2; 20; 52]

Cerebellar hemisphere, D 78 0.009 8.87 0.030 [24; -64; -50]

Middle frontal gyrus, S Inferior frontal gyrus (pars triangularis), S 207 <0.001 8.26 0.057 [-40; 44; 2]

Middle frontal gyrus, D 261 <0.001 8.1 0.067 [42; 30; 32]

Supramarginal gyrus, D 276 <0.001 7.86 0.086 [48 -42; 50]

Angular gyrus, D 7.1 0.192 [46 -50; 46]

Superior parietal lobule, D 7 0.213 [40 -46; 42]

Superior parietal lobule, S 317 <0.001 6.96 0.223 [-36; -48; 46]

Supramarginal gyrus, S 6.19 0.47 [-46; -38; 48]

Anterior insula, S 126 0.001 6.93 0.229 [-32; 20; -4]

Anterior insula D 96 0.004 6.89 0.238 [34; 22; -2]

Middle frontal gyrus, S 149 <0.001 6.86 0.245 [-44; 24; 30]

Note: MNI — Montreal Neurological Institute, FWE — family-wise error, T — T-statistics, D (dexter) — right, S (sinister) — left.

containing 23 voxels. The temporal pattern of the BOLD signal intensity change in this region is similar to the curve obtained for the dorsolateral prefrontal cortex. It contains two peaks, 6 and 18 seconds from the beginning of the block. Between the peaks, a decrease in activation was observed with a local minimum at 13 seconds from the beginning of the block. Similar to the left dorsolateral prefrontal cortex, an increase in signal intensity begins with the beginning of the block, and the minimum is detected 24-29 seconds from the beginning of the module (Fig. 4).

In the left middle occipital gyrus a cluster of activation was selected at the intersection with a sphere with a radius of 10 mm and the center at [-24; -82; 1] (MNI space), containing 33 voxels. In this region of interest an increase in the BOLD signal intensity is observed from the beginning of the task. The maximum is observed 5-6 seconds after the beginning of the block, then the signal level decreases. The minimum is registered 11 seconds from the beginning followed by a slight increase in signal intensity with a subsequent decrease in activation (Fig. 4).

DISCUSSION

In this study, the activation in regions of cerebral cortex at different stages of the modified Sternberg verbal working

memory task was assessed in 19 healthy volunteers using fMRI. Temporal pattern of changes in the left dorsolateral prefrontal, left posterior parietal and left occipital cortex activation was analyzed. During the encoding and retrieval periods activation in prefrontal cortex was observed (bilaterally and in the right hemisphere, respectively), as well as bilateral activation in parietal cortex, cerebellar hemispheres and occipital cortex. During the maintenance period the activation in FPCN components was found: symmetrical regions of the prefrontal and posterior parietal cortex, insular cortex, and medial parts of the superior frontal gyrus. The maximal activation of the left prefrontal and parietal cortex was observed during the encoding period and at the end of the maintenance period, between them it remained relatively constant, while for the occipital cortex, the peak of activation matched to the encoding period. Thus, according to the results of fMRI analysis, activation of prefrontal and posterior parietal cortex (components of FPCN) was detected during maintenance of verbal information in working memory.

Regions of activation in our study were generally consistent with the data of other researchers who used a similar modification of the Sternberg task for verbal working memory assessment [9]. Activation of the upper and lower parts of the left prefrontal cortex , parietal cortex bilaterally and the occipital cortex was observed during the encoding period. However, in our study, we detected additional activation in the right

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

M

Activation during the maintenance period

FPCN

Fig. 3. Comparison of activation map during the maintenance period with FPCN. FPCN Is the fronto-parietal control network (according to [24])

Table 3. Regions of significant activation during the retrieval period of the modified Sternberg task for verbal working memory

Region Cluster level Peak level Peak coordinates, MNI [x; y; z]

Volume, voxels рFWE Т рFWE

Anterior and medial orbital gyrus, D 1074 <0.001 11.36 0.002 [26; 44; -18]

Middle temporal gyrus, D 309 <0.001 9.24 0.018 [50; -40; -8]

Cerebellar hemisphere, D 206 <0.001 9 0.023 [24; -52; -48]

Cerebellar hemisphere, S 294 <0.001 8.76 0.029 [-48; -70; -30]

Fusiform gyrus, S 6.87 0.212 [-36; -72; -22]

Frontal insular cortex, D 171 <0.001 8.03 0.062 [42; 4; 6]

Operculum, D 5.68 0.64 [46; -2; 18]

Supramarginal gyrus, D 1465 <0.001 8.02 0.063 [56; -32; 52]

Postcentral gyrus, S 263 <0.001 7.97 0.067 [-42; -34; 66]

Supramarginal gyrus, S 6.04 0.482 [-48; -40; 54]

Middle frontal gyrus, D 201 <0.001 7.95 0.068 [38; 18; 48]

Supramarginal gyrus, S 588 <0.001 7.56 0.103 [-48; -28; 28]

ПPostcentral gyrus, S

Operculum, S

Cerebellar vermis 259 <0.001 6.86 0.214 [8; -60; -8]

Lingual gyrus, bilateral

Cerebellar hemisphere, S 97 0.006 6.60 0.282 [-2; -42; 0]

Note: MNI — Montreal Neurological Institute, FWE — family-wise error, T — T-statistics, D (dexter) — right, S (sinister) — left.

prefrontal cortex and the cerebellar hemispheres (bilaterally). In addition, activation of the parietal cortex was more prominent, and the activation region in the visual cortex was relatively small. During the maintenance period we detected activation regions in the parietal and prefrontal cortex bilaterally, while in the study described above the activation was observed mainly in the left hemisphere. At the retrieval period the regions of activation in the parietal cortex on the left and bilateral regions

of the operculum are consistent. Compared to [9], additional regions of activation were found in the right parietal cortex and the temporal lobe, but no activation of the left prefrontal cortex was observed. This difference could be due to the different task complexity (the researchers [9] used sequences of 4-6 letters, unlike 7 letters in our study) or due to the different duration of the encoding, maintenance and retrieval periods (2.16; 6.48 and 2.16 seconds respectively [9]). The differences could

Fig. 4. Temporal dynamics of BOLD signal intensity change in the clusters of activation in the left middle frontal gyrus, left inferior parietal lobule and left middle occipital gyrus, averaged over all participants and task blocks. All regions are shown on the surface of the left hemisphere. Solid line — periodically interpolated simulated signal, dashed line — filtered measured signal corrected for non-sphericity and head motion. The BOLD signal intensity is measured in arbitrary units. The origin corresponds to the rest state at the beginning of the block. Interpretation of these data in terms of neuronal events involved the subtraction of 5 seconds from the onsets of BOLD events to account for the hemodynamic delay

also result from the different task performance strategies: the activation pattern in the maintenance period of our experiment was similar to the activation detected during the working memory task performance with articulatory suppression [26], which could indicate the use of non-articulatory maintenance mechanisms by participants.

Activation of clusters in the functionally different regions during the encoding and maintenance periods can be explained by the combination of several types of activity: visual information processing, subvocal speaking, choosing the strategy of maintenance, as well as attention and motor activity. Thus, during the encoding period, significant activation of the occipital cortex regions can be observed, associated with the visual information processing [9]. The activation of left inferior frontal gyrus region at encoding period might be associated with the verbal information retention [12]. Bilateral activation of parietal cortex at the same time with activation of the inferior frontal gyrus is consistent with the preferential activation of dorsal attention network (DAN) during the encoding period [27, 28]. Activation of a region in the cerebellar hemisphere has also been shown in a meta-analysis and is consistent with the hypothesis that the right cerebellar hemisphere is involved in the verbal working memory [11].

The activation of FPCN components was revealed during the maintenance period [24], which is consistent with the data of a meta-analysis demonstrating a gradual increase in FPCN activity and a decrease in DAN activity from the encoding period to the maintenance period in delayed response tasks [27]. At the same time, a decrease in the total number of activation regions during the maintenance period compared to encoding and retrieval may be associated with fewer types of activity involved (only maintenance of information).

During the retrieval period, activation was detected in the parietal lobes bilaterally with the predominance of the right one, which might be due to attention processes [13, 29]. Activation of the right middle temporal gyrus may be due to attention switching [29]. As in the encoding period, activation of the occipital cortex might be associated with visual perception. Observed activation in the right anterior insula is consistent with the evidence of their control function [28]. As opposed to the meta-analysis results [27], in our study, the activation of both FPCN and DAN components has been confirmed in the retrieval period. Activation of the right prefrontal cortex during the retrieval period with presentation of a verbal stimulus emphasizes the ambiguity of the idea about the lateralization of the regions with a functional role in working memory.

According to the results of the BOLD signal changes analysis with a correction for the hemodynamic response delay, activation of the left dorsolateral prefrontal cortex achieves its

maximum in the middle of the encoding period. Then, during the maintenance period it experiences a slight decrease remaining relatively high, after that it increases slightly at the end of the maintenance period and decreases during the retrieval period (neural events in the range of 12-15 seconds from the beginning of the task). This is mostly consistent with the data of other researchers [9, 18, 20] and the hypothesis of the importance of persistent neural activity in working memory [4]. Compared to dorsolateral prefrontal cortex, the parietal cortex activation decreases more strongly in the middle of the maintenance period. Then it experiences a slight increase with a maximum at the beginning of the retrieval period, which is consistent with the results of another study and is especially characteristic for the high-load Sternberg task (6-9 letters) [18]. On the contrary, the occipital cortex activation achieves its maximum during the encoding period and then decreases, which suggests its role in the processes of visual stimulus perception, but not the working memory itself [9, 13, 18].

When interpreting the results, it is necessary to take into account the limitations of the method used to determine the task stages, during which every brain region is active. Since the task periods follow each other in a fixed order, it is impossible to separate the effects of each state on the BOLD signal completely using a model based on the canonical hemodynamic response function, in the absence of detailed information about this function for each subject and brain region. A number of approaches has been proposed in the literature to reduce the impact of this problem, but has not completely solved it [20, 30]. Other limitations include different individual task performance strategies, which could affect the group analysis. The short duration of the encoding and retrieval periods compared with the rest period should also be noted.

CONCLUSION

The study revealed the activation of prefrontal and posterior parietal cortex (parts of the frontoparietal control network, FPCN) in healthy volunteers during the verbal working memory delayed response task. The results are consistent with the ideas about the key role of prefrontal cortex in the working memory mechanisms of information maintenance, while during the encoding and retrieval periods predominant activation of the parietal cortex was found. Despite the limitations of the method, the study results are consistent with the results of other research using the Sternberg task for the investigation of visual and verbal working memory by fMRI, which speaks in favor of the reproducibility and possible use the results to study the distinct working memory mechanisms.

References

1. Archer JA, Lee A, Qiu A, Chen SA. Working memory, age and education: A lifespan fMRI study. PLoS One. 2018; 13 (3): e0194878.

2. Chai WJ, Abd Hamid AI, Abdullah JM. Working Memory From the Psychological and Neurosciences Perspectives: A Review. Front Psychol. 2018; 9: 401.

3. D'Esposito M, Postle BR. The cognitive neuroscience of working memory. Annu Rev Psychol. 2015; 66: 115-42.

4. Fuster JM. Unit activity in prefrontal cortex during delayed-response performance: neuronal correlates of transient memory. J Neurophysiol. 1973; 36 (1): 61-78.

5. Fuster JM, Bauer RH, Jervey JP. Functional interactions between inferotemporal and prefrontal cortex in a cognitive task. Brain Res.

1985; 330 (2): 299-307.

6. Riley MR, Constantinidis C. Role of Prefrontal Persistent Activity in Working Memory. Front Syst Neurosci. 2016; 9: 181.

7. Egli T, Coynel D, Spalek K, Fastenrath M, Freytag V, Heck A et al. Identification of Two Distinct Working Memory-Related Brain Networks in Healthy Young Adults. eNeuro. 2018; 5 (1): pii: ENEURO.0222-17.2018.

8. Cohen JR, Sreenivasan KK, D'Esposito M. Correspondence between stimulus encoding- and maintenance-related neural processes underlies successful working memory. Cereb Cortex. 2014; 24 (3): 593-9.

9. Narayanan NS, Prabhakaran V, Bunge SA, Christoff K, Fine EM,

Gabrieli JD. The role of the prefrontal cortex in the maintenance of verbal working memory: an event-related FMRI analysis. Neuropsychology. 2005; 19 (2): 223-32.

10. Beynel L, Davis SW, Crowell CA, Hilbig SA, Lim W, Nguyen D, et al. Online repetitive transcranial magnetic stimulation during working memory in younger and older adults: A randomized within-subject comparison. PLoS One. 2019; 14 (3): e0213707.

11. Emch M, von Bastian CC, Koch K. Neural correlates of verbal working memory: An fMRI meta-analysis. Front Hum Neurosci. 2019; 13: 180.

12. Rottschy C, Langner R, Dogan I, Reetz K, Laird AR, Schulz JB, et al. Modelling neural correlates of working memory: a coordinate-based meta-analysis. Neuroimage. 2012; 60 (1): 830-46.

13. Xu Y. Reevaluating the Sensory Account of Visual Working Memory Storage. Trends Cogn Sci. 2017; 21 (10): 794-815.

14. Sternberg, S. High-speed scanning in human memory. Science. 1966; 153 (3736): 652-4.

15. Brunoni AR, Vanderhasselt MA. Working memory improvement with non-invasive brain stimulation of the dorsolateral prefrontal cortex: a systematic review and meta-analysis. Brain Cogn. 2014; 86: 1-9.

16. Jiang Y, Guo Z, Xing G, He L, Peng H, Du F, et al. Effects of High-Frequency Transcranial Magnetic Stimulation for Cognitive Deficit in Schizophrenia: A Meta-Analysis. Front Psychiatry. 2019; 10: 135.

17. Sack AT, Kadosh RC, Schuhmann T, Moerel M, Walsh V Goebel R. Optimizing functional accuracy of TMS in cognitive studies: A comparison of methods. J Cogn Neurosci. 2009; 21 (2): 207-21.

18. Wen X, Wang H, Liu Z, Liu C, Li K, Ding M, et al. Dynamic Top-down Configuration by the Core Control System During Working Memory. Neuroscience. 2018; 391: 13-24.

19. Woodward TS, Feredoes E, Metzak PD, Takane Y, Manoach DS. Epoch-specific functional networks involved in working memory. Neuroimage. 2013; 65: 529-39.

20. Motes MA, Rypma B. Working memory component processes:

Isolating BOLD signal changes. Neuroimage. 2010; 49 (2): 1933-41.

21. Hamidi M, Tononi G, Postle BR. Evaluating the role of prefrontal and parietal cortices in memory-guided response with repetitive transcranial magnetic stimulation. Neuropsychologia. 2009; 47 (2): 295-302.

22. Oldfield RC. The assessment and analysis of handedness: The Edinburgh inventory. Neuropsychologia. 1971; 9 (1): 97-113.

23. Kanal E, Barkovich AJ, Bell C, Borgstede JP, Bradley Jr. WG, Froelich JW, et al. ACR guidance document on MR safe practices: 2013. J Magn Reson Imaging. 2013; 37 (3): 501-30.

24. Yeo BT, Krienen FM, Sepulcre J, Sabuncu MR, Lashkari D, Hollinshead M, et al. The organization of the human cerebral cortex estimated by intrinsic functional connectivity. J Neurophysiol. 2011; 106 (3): 1125-65. Available from: http://www.freesurfer. net/fswiki/CorticalParcellation_Yeo2011.

25. xjView. Version 9.6 [software]. Available from: http://www. alivelearn.net/xjview.

26. Trapp S, Mueller K, Lepsien J, Kraemer B, Gruber O. Different neural capacity limitations for articulatory and non-articulatory maintenance of verbal information. Exp Brain Res. 2014; 232 (2): 619-28.

27. Kim H. Neural activity during working memory encoding, maintenance, and retrieval: A network-based model and meta-analysis. Hum Brain Mapp. 2019; 40 (17): 4912-33.

28. Soreq E, Leech R, Hampshire A. Dynamic network coding of working-memory domains and working-memory processes. Nat Commun. 2019; 10 (1): 936.

29. Nee DE, Brown JW, Askren MK, Berman MG, Demiralp E, Krawitz A, et al. A meta-analysis of executive components of working memory. Cereb Cortex. 2013; 23 (2): 264-82.

30. Rozovskaya RI, Pechenkova EV, Mershina EA, Machinskaya RI. fMRI Study of Retention of Images with Different Emotional Valence in the Working Memory. Psychology. Journal of the Higher School of Economics. 2014; 11 (1): 27-48. Russian.

Литература

1. Archer JA, Lee A, Qiu A, Chen SA. Working memory, age and education: A lifespan fMRI study. PLoS One. 2018; 13 (3): e0194878.

2. Chai WJ, Abd Hamid AI, Abdullah JM. Working Memory From the Psychological and Neurosciences Perspectives: A Review. Front Psychol. 2018; 9: 401.

3. D'Esposito M, Postle BR. The cognitive neuroscience of working memory. Annu Rev Psychol. 2015; 66: 115-42.

4. Fuster JM. Unit activity in prefrontal cortex during delayed-response performance: neuronal correlates of transient memory. J Neurophysiol. 1973; 36 (1): 61-78.

5. Fuster JM, Bauer RH, Jervey JP. Functional interactions between inferotemporal and prefrontal cortex in a cognitive task. Brain Res. 1985; 330 (2): 299-307.

6. Riley MR, Constantinidis C. Role of Prefrontal Persistent Activity in Working Memory. Front Syst Neurosci. 2016; 9: 181.

7. Egli T, Coynel D, Spalek K, Fastenrath M, Freytag V, Heck A et al. Identification of Two Distinct Working Memory-Related Brain Networks in Healthy Young Adults. eNeuro. 2018; 5 (1): pii: ENEUR0.0222-17.2018.

8. Cohen JR, Sreenivasan KK, D'Esposito M. Correspondence between stimulus encoding- and maintenance-related neural processes underlies successful working memory. Cereb Cortex. 2014; 24 (3): 593-9.

9. Narayanan NS, Prabhakaran V, Bunge SA, Christoff K, Fine EM, Gabrieli JD. The role of the prefrontal cortex in the maintenance of verbal working memory: an event-related FMRI analysis. Neuropsychology. 2005; 19 (2): 223-32.