Takatoshi Satake
Ai Taki1,3
Kazumi Kasahara
Daisuke Yoshimaru
Tomokazu Tsurugizawa- 1Human Informatics and Interaction Research Institute, National Institute of Advanced Industrial Science and Technology (AIST), Ibaraki, Japan
- 2Graduate School of Comprehensive Human Science, University of Tsukuba, Tsukuba, Japan
- 3Faculty of Engineering, Information and Systems, University of Tsukuba, Tsukuba, Japan
- 4Division of Regenerative Medicine, The Jikei University School of Medicine, Tokyo, Japan
The N-back task is widely used to investigate working memory. Previous functional magnetic resonance imaging (fMRI) studies have shown that local brain activation depends on the difficulty of the N-back task. Recently, changes in functional connectivity and local activation during a task, such as a single-hand movement task, have been reported to give the distinct information. However, previous studies have not investigated functional connectivity changes in the entire brain during N-back tasks. In this study, we compared alterations in functional connectivity and local activation related to the difficulty of the N-back task. Because structural connectivity has been reported to be associated with local activation, we also investigated the relationship between structural connectivity and accuracy in a N-back task using diffusion tensor imaging (DTI). Changes in functional connectivity depend on the difficulty of the N-back task in a manner different from local activation, and the 2-back task is the best method for investigating working memory. This indicates that local activation and functional connectivity reflect different neuronal events during the N-back task. The top 10 structural connectivities associated with accuracy in the 2-back task were locally activated during the 2-back task. Therefore, structural connectivity as well as fMRI will be useful for predicting the accuracy of the 2-back task.
Introduction
Working memory is essential for the cognitive process that keeps small amounts of information active for later recall and manipulates limited amounts of information for current computation (Daneman, 1980; Baddeley, 1992). Working memory capacity is the amount of information that can be retained in mind for a short period. The amount of accessible information is determined by working memory capacity, which depends on the individual (Wilhelm et al., 2013). The n-back task, a continuous-recognition measure that uses sequences of stimuli, such as letters, is a simple tool for monitoring working memory processes (Kirchner, 1958). For this task, the participants were required to monitor a series of stimuli and respond whenever a stimulus was presented, the same as in the N trials. The difficulty of the n-back task varies with N (i.e., the 3-back task is more complicated than the 1-back task) (Picchioni et al., 2007; Lamichhane et al., 2020). Previous studies using functional magnetic resonance imaging (fMRI), which detects blood oxygenation level-dependent (BOLD) responses coupled with neuronal activation, have revealed the locations of local activation for N-back tasks. In general, the BOLD response is calculated using a general linear model (GLM), which is a model of the temporal convolution of the hemodynamic response function to the task period. Owen et al. reported local activation in frontoparietal regions, including the ventrolateral, dorsolateral, and frontopolar prefrontal cortices, in addition to the dorsal cingulate and premotor cortices, during an N-back task (Owen et al., 2005). Another meta-analysis revealed the largest clusters in the prefrontal and parietal cortices of the left hemisphere, including the middle frontal gyrus (MFG) and inferior parietal lobule, in young adults during an N-back task. Other areas include the medial frontal gyrus, insula, nuclei of the basal ganglia, and cerebellum (Yaple et al., 2019). Several studies have evaluated the effect of the N-back task difficulty by changing the number of N in the N-back task on the BOLD responses in the prefrontal and posterior parietal cortices. The amplitudes of the evoked responses in the frontal and parietal cortices increased as the number increased from 1- to 3-back (Braver et al., 1997; Jaeggi et al., 2003). Another study explored the N-back task, with scores ranging from N = 1 to N = 6. The BOLD response of the lateral prefrontal cortex was found to increase from 1-back to 3-back and then plateau until 6-back (Lamichhane et al., 2020). These results indicate that the difficulty of N-back task affects the degree of working memory processing in the brain.
These previous studies investigated local activation during the N-back task but did not assess functional connectivity. The functional connectivity is derived from the synchronization of neuronal oscillations between anatomically separated regions. Importantly, functional connectivity is related to cognitive function (Power et al., 2011). The generalized psychophysiological interaction (gPPI) is used to investigate functional connectivity during cognitive tasks (McLaren et al., 2012). GLM analysis indicates local neuronal activation, whereas gPPI showed a change in the functional connectivity between the affected brain regions during the task. This approach enables us to investigate functional connectivity, which is more sensitive and specific to the task phase (Tsurugizawa et al., 2023).
In addition to fMRI, the white matter microstructure, which can be estimated using diffusion tensor imaging (DTI), is essential for cognitive processing and memory. Several DTI studies on working memory tasks have revealed positive correlations between fractional anisotropy in the frontoparietal white matter and working memory task performance in children (Olesen et al., 2003; Klingberg, 2006; Darki and Klingberg, 2015) and adults (Schulze et al., 2011). The anatomical connections of the white matter between anatomically separated regions, called structural connectivity, can be calculated using DTI. Previous studies have suggested its use in predicting working memory (McKenna et al., 2015). However, to date, no study has compared the structural connectivity, functional connectivity, and local activation in the N-back task.
In the present study, we compared the functional connectivity, structural connectivity, and local activation during an N-back task. We used GLM and gPPI analyses to investigate local activation and functional connectivity depending on the difficulty of the N-back task. Because we found that the 2-back task was better for investigating working memory from behavior, the functional and structural connectivity associated with accuracy were also investigated.
Materials and methods
Participants
Twenty-five healthy young adults (14 males and 11 females; age 24.3 ± 7.32 years) participated in the fMRI experiment. All the participants were right-handed and had no history of alcohol abuse, neurological diseases, or learning disabilities. All experimental procedures and protocols were approved by the Institutional Review Board of the National Institute of Advanced Industrial Science and Technology.
Behavioral analysis
R Statistical Software (v4.3.1, R Core Team, 2023), RStudio (2023.06.1, Rstudio Team, 2023), Coin [v1.4.2, (Hothorn et al., 2006)] and Tidyverse [v2.0.0, (Wickham et al., 2019)] were used for statistical behavioral analyses.
fMRI N-back task and acquisition parameters
Participants were scanned with a 32-channel phased array receiving head coil (3.0 T scanner, Philips, Netherlands). The participants wore earplugs to reduce scanner noise, and foam padding was used to suppress head motion. T1-weighted 3D magnetization prepared-rapid Gradient Echo (MPRAGE) images were obtained using the following parameters: repetition time (TR) = 11 ms, echo time (TE) = 5.1 ms, flip angle = 8°, matrix = 368 × 315 × 257, and resolution = 0.7 × 0.76 × 0.7 mm3/voxel. The fMRI data were acquired using gradient echo Echo-Planar Imaging (EPI) with the following conditions: TR = 1,500 ms, TE = 30 ms, flip angle = 80°, multiband = 2, matrix = 76 × 76 × 44, resolution = 2.5 × 2.5 × 2.5 mm3/voxel. The fMRI data for each task were acquired for 8 min 45 s (350 volumes), and the resting-state fMRI data were acquired for 10 min 30 s (420 volumes).
Following the fMRI experiment, DTI data were acquired with diffusion-weighted spin echo EPI and the following parameters: TR = 7,263 ms, TE = 95 ms, matrix = 112 × 110 × 55, resolution = 2.0 × 2.0 × 2.0 mm3/voxel with 1 b = 0 s/mm2 image and with 64 diffusion directions defined evenly across the sphere with a diffusion weighting of b = 1,000 s/mm2.
N-back task with fMRI
The N-back task was performed with five letters: “a,” “b,” “c,” “d,” and “e” (Figure 1). The letters were randomly presented on a computer screen located outside the MRI bore and the participants viewed the letters through a tilted mirror. The tasks were performed in the order 1- of 3-back task. Each of the letters was presented for 2.5 s, twenty times in each stimulation block (in total 50 s), and the “+” mark was presented for the 40-s resting period (Figure 1). Five stimulation blocks were used in each experiment. When participants recognized that the current letter was the same as the previous letter, they pushed the button. The timing of the alphabet presentation and data collection was determined using PsychoPy.1 fMRI scanning with the “+” mark throughout MRI scan was performed as the negative control (Kumari et al., 2006; Yaple et al., 2021).
Figure 1. Schematic of the N-back task and resting-state fMRI. Schematic representation of the (A) N-back task and (B) resting-state fMRI. Five letters (“a,” “b,” “c,” “d,” and “e”) were presented randomly during the stimulation block. Letter presentation was repeated 20 times for each block. The “+” mark was presented during the resting period.
Preprocessing of the fMRI data
Statistical parametric mapping SPM12 software (Welcome Trust Center for Neuroimaging, United Kingdom) was used in the preprocessing steps, including slice timing correction, motion correction by realignment, normalization, and smoothing with a Gaussian filter (8 × 8 × 8 mm2/voxel of half-width at half-maximum). Following the preprocessing with SPM12, fMRI data were detrended, and slow periodic fluctuations were extracted using a bandpass filter (0.008–0.09 Hz) with CONN toolbox.2
General linear model (GLM) analysis
For the first level (fixed-effect analysis), an onset regressor defined the onset of the N-back task, and the block length was set to 50 s for each session. The hemodynamic response was modeled using the canonical hemodynamic response function in SPM12. Six head-movement parameters estimated during the realignment process were used as regressors to determine the effect of head motion. The resulting contrast images were calculated by setting the targeting regressor (regressor of the N-back task) to 1 and the others to zero in N-back tasks and no N-back task. The calculated contrasts were then used for a second-level analysis (random effect analysis). The contrasts of 1, 2, and 3-back were compared with no N-back contrast. A high-pass filter of 1/128 Hz was used to remove this detrending. The contrasts of the 1-, 2-, and 3-back tasks were compared with the no N-back contrast. Model estimation was performed using a paired t-test to compare each task with the negative control. The significant voxels were thresholded at p < 0.05, and family-wise error (FWE) was corrected with a cluster size. The activated regions were determined by counting the voxels in each region of interest (ROI). To investigate the correlation between fMRI beta values and task performance, the beta values of individuals within ROI were extracted from the contrast image that were calculated in first level analysis and were averaged.
ROI-To-ROI gPPI analysis
The ROI-to-ROI gPPI was performed using the CONN toolbox to assess the brain regions that interact in a task-dependent manner. A total of 132 ROIs from the automated anatomical labeling atlas were used for ROI analysis. ROI-to-ROI gPPI analysis was performed for every possible pairwise combination of the selected ROIs for each participant. The predictors included the estimated time course of the task (psychological term), time-series regressor of BOLD signal changes in each ROI (physiological term), and interactions between psychological and physiological terms (PPI term), which consisted of the product of the psychological term multiplied by the physiological term. For the first-level analysis, a PPI regressor (PPI term) was generated for each condition (resting and task) as the product of the ROI time series multiplied by the task effect, and the beta weight was calculated for all ROIs. A random effects analysis was used across participants at the group level, and a one-sample t-test was conducted to compare ROI-based connectivity at rest and during the tasks.
DTI analysis
DTI data were preprocessed by denoising, eddy current correction, and motion correction using MRtrix3 (Tournier et al., 2019). The DTI data were denoised using the dwidenoise command. Eddy currents and motion corrections were performed using the dwifslprepcroc command (Ning et al., 2021). A b0 image was registered on a T1-weighted image. T1-weighted images were then registered onto the standard MNI space with a nonlinear transformation using FNIRT in FSL (University of Oxford, United Kingdom) (Park et al., 2015). Nonlinear registration of the native diffusion space to the MNI standard space was performed using these parameters. Structural connectivity was computed from preprocessed DTI data using the DSI Studio software (Tsurugizawa et al., 2020). The preprocessed DTI data were then used for fiber tracking. First, the differential tractogram was obtained by placing 1,000,000 seeding points in the white matter using the following parameters: QA threshold = 0.1, angular threshold = 45°, step size = 0.5 mm, minimum length = 1 mm, and maximum length = 300 mm (Bauer et al., 2017). The same ROIs as gPPI were used for brain parcellation, and the connectivity matrix was calculated using the connecting tracks. The connectivity matrix was calculated by calculating the number of tracts passing through the two ROIs. The number of fibers is proportional to the fiber length from a tracking algorithm when fibers are calculated between two ROIs. The number of tracts between two ROIs should be scaled by the fiber length to account for this effect (Latora and Marchiori, 2001; Cheng et al., 2012). The number of tracts between two ROIs should be scaled by the fiber length to account for this effect. This number was then normalized by multiplying it with the sum of the inverses of the lengths.
Statistical tests
We performed simple regression analysis for task accuracy acquired during fMRI scanning as behavior analyses. Friedman non-parametric tests were used to compare task accuracy and response time between tasks. As a post hoc test, Wilcoxon signed-rank tests were conducted for each of the three combinations: 1-back vs. 2-back, 2-back vs. 3-back, and 1-back vs. 3-back.
We performed a paired t-test with Bonferroni correction (alpha significance of 0.05/3 = 0.017) on the accuracy and response time of the 1-, 2-, and 3-back tasks. For ROI-to-ROI gPPI, the statistical significance of ROI-ROI connectivity was assessed using a paired t-test or two-sample t-test with a corrected threshold of p < 0.05 and FDR-corrected using the CONN toolbox. The least-squares regression analysis between beta values, functional connectivity, structural connectivity, and accuracy of each task was performed.
Results
Task performance for the N-back task
Overall, accuracy negatively correlated with N-back task difficulty, whereas response time positively correlated with N-back task difficulty. The correlation between accuracy and response time was calculated using the least-squares method (R2 was 1.18 × 10−2 for 1-back, 0.18 for 2-back, and − 1.93 × 10−2 for 3-back) (Figure 2). The accuracies for 1-back, 2-back, and 3-back were 95.70 ± 4.98%, 75.12 ± 9.69%, and 52.75 ± 15.18%, respectively. The response times for 1-back, 2-back, and 3-back were 0.59 ± 0.09 s, 0.71 ± 0.17 s, and 0.90 ± 0.20 s, respectively. The 1-back task led to higher accuracy and shorter response time than the 2-back and 3-back tasks. Significant differences between groups were found for both response accuracy and reaction time (p-value = 1.39 × 10−11, chi-squared = 50 for accuracy, value of p = 7.12 × 10−9, chi-squared = 37.52 for response time, based on the Friedman’s test). Post-hoc analysis using the Wilcoxon test revealed significant differences for all 1-backs vs. 2-back (accuracy: p-value = 5.96 × 10−8, response time: p-value = 5.39 × 10−5), 2-back vs. 3-back (p-value = 5.96 × 10−8 for accuracy, p-value = 6.56 × 10−6 for response time), and 1-back vs. 3-back (p-value = 5.96 × 10−8 for accuracy, p-value = 1.29 × 10−7 for response time). These results were consistent with those of a previous study (Lamichhane et al., 2020).
Figure 2. Accuracy and response time for the 1-, 2-, and 3-back tasks. Correlation between accuracy and response time (ms) for tasks 1-back (red), 2-back (green), and 3-back (blue). The solid line indicates fitting using the least-squares method.
Altered functional connectivity during the N-back task
Significant changes in functional connectivity were observed in the 1-, 2-, and 3-back tasks (Figure 3). No significant changes in the functional connectivity were observed during the 1-back task. In contrast, functional connectivity in several brain regions changed differently in the 2-back and 3-back tasks. During the 2-back task, functional connectivity increased among several brain regions, including the left anterior supramaginal gyrus (aSMG), right paracingulate gyrus (PaCiG), bilateral frontal pole (FP), left superior lateral occipital cortex (sLOC), posterior cingulate gyrus (postCG), and bilateral temporal occipital fusiform cortex (ToFusC) (Figure 3A and Table 1). Functional connectivity decreased in several brain regions, including the right anterior inferior temporal gyrus (aITG), left temporal pole (TP), and bilateral occipital pole (OP) during the 2-back task compared with the rest period.
Figure 3. Significant changes in functional connectivity during the 2-and 3-back task. Circle plot of significant changes in functional connectivity for the (A) 2-back task and (B) 3-back task compared to the no-N-back task (p < 0.05, FDR-corrected). Positive color indicates an increased functional connectivity T-value. Negative colors indicate T-values of decreased functional connectivity. aTFusC, temporal fusiform cortex, anterior division; aITG, inferior temporal gyrus, anterior division; aMTG, middle temporal gyrus, anterior division; aSMG, supramarginal gyrus, anterior division; Cereb2, cerebelum crus2; CO, central opercular cortex; FP, frontal pole; HG, Heschl’s gyrus; IC, insular cortex; ICC, intracalcarine cortex; LG, lingual gyrus; MidFG, middle frontal gyrus; OFusG, occipital fusiform gyrus; OP, occipital pole; PaCiG, paracingulate gyrus; PC, cingulate gyrus, posterior division; pMTG, middle temporal gyrus, posterior division; PO, parietal operculum cortex; PostCG, postcentral gyrus; pPaHC, parahippocampal gyrus, posterior division; PreCG, precentral gyrus; pTFusC, temporal fusiform cortex, posterior division; SFG, superior frontal gyrus; sLOC, lateral occipital cortex, superior division; SMA, supplementary motor area; SubCalC, subcallosal cortex; TOFusC, temporal occipital fusiform cortex; toITG, inferior temporal gyrus, temporooccipital part; TP, temporal pole.
Table 1. The significant changes (p < 0.05, FDR-corrected) in functional connectivity for 2-back task.
The 3-back task evoked significant changes in functional connectivity in distinct brain regions compared with the 2-back task. Functional connectivity increased significantly in the right central/parietal opercular cortex (CO/PO) and right insular cortex (IC) (Figure 3B and Table 2). In particular, the right IC and bilateral opercula showed increased functional connectivity with the right frontal lobe and cingulate cortex. Functional connectivity between the right anterior medial temporal gyrus (aMTG) and right IC, between the right posterior middle temporal gyrus (pMTG) and left PO, and between the right CO and subcallosal cortex (SubCalC) increased during the 3-back task. Regions related to motor regulation (the left precentral gyrus (preCG) and supplementary motor area (SMA)) showed increased functional connectivity with the left accumbens and right occipital fusiform gyrus (OFusG). The basal ganglia also exhibited decreased functional connectivity in several regions. Functional connectivity between the bilateral superior frontal gyrus (SFG) and left caudate and between the left putamen and bilateral PaCiG was significantly decreased. Functional connectivity was also significantly decreased between the visual cortex and occipital fusiform cortex, including the bilateral lingual gyrus (LG) and bilateral intracranial compliance (ICC).
Table 2. The significant changes (p < 0.05, FDR-corrected) in functional connectivity for 3-back task.
Bold response during the N-back task
BOLD responses were observed in the bilateral MFG and bilateral SFG in all n-back tasks (Figure 4 and Tables 3–7). Positive BOLD responses were observed in the bilateral superior parietal lobule (SPL), bilateral supramarginal gyrus, posterior division, bilateral sLOC, bilateral frontal operculum cortex, bilateral vermis 4 5, right inferior frontal gyrus, pars opercularis, and right angular gyrus in the 2-back and 3-back tasks, but not in the 1-back task. Positive BOLD responses were observed only in the right thalamus and bilateral caudates in the 3-back task. A negative BOLD response in the bilateral FP, bilateral precuneus cortex, and left paracingulate gyrus (PaCiG) was observed in the 2-back and 3-back tasks; however, BOLD signals in no regions were negatively changed in the 1-back task. Negative BOLD responses were only observed in the left parahippocampal gyrus, posterior division, and hippocampus during the 3-back task.
Figure 4. BOLD signal increase during the N-back task. Significant BOLD signal changes in the (A) 1-back, (B) 2-back, and (C) 3-back tasks are shown. The color bar indicates the t-value. aMTG, middle temporal gyrus, anterior division; FO, frontal operculum cortex; FP, frontal pole; IFG, inferior frontal gyrus; MFG, middle frontal gyrus; PaCiG, paracingulate gyrus; SFG, superior frontal gyrus; SMG, supramarginal gyrus; SPL, superior parietal lobule; Ver45, vermis 4 5.
Table 3. The significant BOLD increase during 1-back task.
Table 4. The significant BOLD increase during 2-back task.
Table 5. The significant BOLD decrease during 2-back task.
Table 6. The significant BOLD increase during 3-back task.
Table 7. The significant BOLD decrease during 3-back task.
We compared the areas of significant change in the BOLD responses for the 1-, 2-, and 3-back tasks. The areas of BOLD signal increase and decrease were the largest for the 3-back task (Table 8).
Table 8. Cluster size of the significant change in each task.
The correlation between the fMRI beta value of the bilateral MFG and SFG, which were activated in all N-back tasks, and task accuracy was investigated. In 1-back and 3-back task had small R2 in these ROIs. The R2 in 1-back was 3.42 × 10−2 for right SFG, 4.66 × 10−5 for left SFG, 8.56× 10−5 for right MFG, and 5.27 × 10−3 for left MFG. The R2 in 3-back was 3.61 × 10−2 for right SFG, 6.13 × 10−2 for left SFG, 4.96 × 10−2 for right MFG, and 6.04 × 10−2 for left MFG. The 2-back task showed weak correlation in left MFG (R2 = 1.93 × 10−1), left SFG (R2 = 2.89 × 10−1) and left MFG (R2 = 2.64 × 10−1) and very weak correlation in right SFG (R2 = 9.01 × 10−3). These results indicate that there is no significant correlation between BOLD response and task performance (with an R2 > 0.3).
Correlation between structural/functional connectivity and accuracy during the N-back task
From the behavioral analysis, we used the accuracy of each N-back task to investigate the correlation between behavior and structural connectivity (Table 9). For 1-back task, a connection between cerebellum 3 and right supplementary motor area showed higher correlation (R2 > 0.4), and no connection had R2 > 0.5. There is highest correlation for 2-back task in n-back task (10 connections had R2 > 0.4, and 2 connections had R2 > 0.5). In contrast, 3-back task showed no connection that had R2 > 0.4. Notably, several structural connections had the correlation with accuracy in 1-back and 2-back tasks, i.e., left cerebellum 3 and right SMA, Vermis 3 and right SFG, right SPL and right MFG, and right frontal orbital cortex and right inferior frontal gyrus, pars triangularis. However, it should be noted that R2 for the accuracy of 1-back task was lower than that of 2-back task. Furthermore, a correlation was found between structural connectivity and accuracy in several brain regions activated during the 2-back task, such as the bilateral SPL, right SMA, right angular gyrus, right SFG, right MFG, and left thalamus.
Table 9. The top 10 structural connectivity in each task.
The simple regression analysis between functional connectivity and 2-back and 3-back task accuracies were also performed, but there was no functional connectivity with an R2 > 0.3.
Discussion
In the present study, we compared the local activation, functional connectivity, and structural connectivity associated with the N-back task in the same participants. Functional connectivity during the N-back task was altered in a manner that differed according to the difficulty of the N-back task. In contrast, GLM analysis revealed that the area of significant BOLD response during the N-back task increased in the 3-back task compared to that in the 1- and 2-back tasks, but the location of local activation was similar between the 2-back and 3-back tasks, except for the thalamus and caudate. The differing results between GLM and gPPI analyses can be explained by information on BOLD signals. GLM analysis detected the local BOLD response related to neuronal activation during the N-back task compared with the resting period. In contrast, gPPI analysis detects functional connectivity, which is the synchronization of neuronal fluctuations influenced by the N-back task. A previous study revealed that local activation is not necessarily related to connectivity changes in large-scale brain networks, which is consistent with the results of the present study (Gerchen and Kirsch, 2017). The structural connectivity associated with the accuracy of the 2-back task involved the regions activated during the 2-back task (R2 > 0.3). However, increased functional connectivity using gPPI analysis and local activation did not show a significant correlation with accuracy. Previous study showed the weak correlation between local activation and accuracy (R2 < 0.3), corresponding to current study (Lamichhane et al., 2020). The structural connectivity refers to the white matter pathways between regions, while functional connectivity and local activation refer to the node-to-node interactions and local activity between neurophysiologically active regions. These results indicate that task performance is related to the white matter pathways rather than the neuronal activity. Another study also shows that the training changes the white matter structure, indicating the relationship between white matter structure and task performance (Scholz et al., 2009).
Functional connectivity significantly differed between the 2-back and 3-back tasks, whereas it did not change significantly in the 1-back task (Figure 3A and Table 1). In the 2-back task, functional connectivity increased in the regions related to word perception and working memory. The inferior temporal cortex, including the ITG, temporo-occipital region, and aSMG, is essential for language comprehension and cognition (Kourtidou et al., 2022). The temporal cortex, including the left TP and right aMTG, is associated with word meaning (Mesulam et al., 2013). The posterior cingulate cortex is essential for working memory (Leech and Sharp, 2014). The visual cortex, which is related to word processing, was negatively correlated with several regions of the temporal cortex (left TP, right aMTG, left parahippocampal gyrus, and posterior division). Notably, most brain regions that showed increased functional connectivity were not locally activated (Figure 4), indicating a different manner of neuronal activity between the GLM and gPPI analyses. In contrast, the 3-back task increased the functional connectivity in regions mainly related to verbal processing and motor regulation (Figure 3B and Table 2). Functional connectivity was increased between the insular cortex/operculum, motor cortex, hippocampus, and basal ganglia. Although the operculum and the insular cortex are essential for the N-back task (Namkung et al., 2017), other regions are mainly involved in motor regulation (Lanciego et al., 2012). Functional connectivity between the occipital lobe, which is related to verbal processing, and the hippocampus was decreased in the 3-back task. Furthermore, functional connectivity decreased between the frontal lobe and basal ganglia, which are related to perception and motion control, in the 3-back task. This distinct functional connectivity in the N-back task may be due to the difficulty of the working memory task. Despite the high accuracy for all volunteers, the response time depended on the individual for the 1-back task, indicating that it was easy to perform for all participants. The response time and accuracy varied widely for the 3-back task, indicating that the 3-back task exceeded the working memory capacity of some volunteers. A weak correlation was found between accuracy and response time for the 2-back task. Therefore, we decided that the 2-back task is better to investigate working memory than the 1-, and 3-back tasks. Importantly, the increased functional connectivity for the 3-back task did not include the regions essential for working memory and word meaning, indicating that the 3-back task exceeded the capacity of working memory in some volunteers. The present study also shows the correlation between accuracy and local activations was higher in 2-back task than in 1- and 3-back tasks. This implies that the brain uses primitive brain networks related to intuitive decisions in complex working memory tests.
As described above, the local activation and functional connectivity are triggered by neuronal activation (Drew, 2019) and neuronal oscillation (Takeuchi et al., 2010) influenced by the task. A previous study revealed that local activation is not necessarily similar to connectivity changes in brain networks, which supports the results of the present study (Gerchen and Kirsch, 2017). Functional connectivity is related to the synchronization of neuronal oscillations. A previous study has revealed that local brain activation and task-dependent connectivity changes can be modulated in different directions during motor tasks. Local activation and functional connectivity are two aspects of the laterality of single-hand movements in different manner (Tsurugizawa et al., 2023). These results indicate that a combination of GLM (local activation) and gPPI (functional connectivity) analyses is better suited for investigating the neuronal mechanisms of working memory.
The top 10 structural connectivity associated with accuracy in the 2-back task were locally activated during the 2-back task. Structural connectivity with the right MFG, inferior frontal gyrus, and SPL correlated with accuracy in the 2-back task, and these regions were activated during the 1- and 2-back task. However, the R2 in 1-back task is smaller than 2-back. Because the SFG and MFG showed increased functional connectivity, as well as local activation and structural connectivity, these regions may be key regions for the N-back task. A previous study showed that the clustering coefficients of the right anterior cingulate gyrus and right inferior ventrolateral prefrontal cortex significantly increased from pre-training to post-training in a learning memory task (Caeyenberghs et al., 2016). The structural connectivity that identifies white matter bundles and connects the nodes of functional networks reflects the physical connections between anatomically separated regions. This physical network contributes to the local activation during working memory tasks. A previous study has revealed that structural connectivity is affected by training and working memory (Rogers et al., 2007). Structural connectivity predicts local activation, suggesting the increased sensitivity of this technique in identifying individual differences in task-based fMRI activation (Ekstrand et al., 2020). These results indicate that structural connectivity, rather than functional connectivity, is useful for predicting local activation. Notably, the observed functional connectivity for the 2-back task did not correspond to structural connectivity. Instead, previous study shows that resting-state functional connectivity can be predicted from structural connectivity (Chen and Wang, 2018). However, altered functional connectivity during working memory tasks is not linked to structural connectivity. This indicates that functional connectivity changes during the task provide distinct information compared with structural connectivity and BOLD responses. Previous studies have reported that resting state functional connectivity in several brain regions increases following training for working memory (Cole et al., 2021; Pongpipat et al., 2021; Yeung and Han, 2023) and structural connectivity (Roman et al., 2017). The default mode network also decreases when the working memory load increases (Ginestet and Simmons, 2011). As recent study indicates, the combination of structural and functional connectivity can be better predictor of working memory accuracy (Kim et al., 2021). As the current study investigated functional connectivity in naïve participants, future studies should investigate dynamic changes in structural and functional connectivity following working memory training to assess as the predictor of working memory accuracy.
In conclusion, both GLM (local activation) and gPPI (functional connectivity) analyses revealed a new aspect of working memory. Structural connectivity is useful for predicting working memory accuracy, although functional connectivity and local activation shows weak correlation with it. Therefore, using a comprehensive approach of fMRI and DTI is better for gaining an understanding of the neuronal mechanisms of working memory.
Limitation of the study
Large sample sizes have recently been required to obtain reliable results in fMRI studies, and meta-analyses using this database have been performed for N-back tasks (Owen et al., 2005; Yaple et al., 2019). Although the brain region of local activation observed during the N-back task is consistent with these studies, and the difficulty-dependent increase in the area of local activation has also been reported in previous studies (Lamichhane et al., 2020), the results of the current study should be investigated with a larger sample size in the future.
In the present study, 1, 2, 3-back tasks and no N-back task were conducted in different sessions. Therefore, it was not possible to compare the fMRI response in each n-back task directly. In the future, this analysis should be assessed with the data sets that includes the fMRI data of no-back, 1, 2, and 3-back tasks in the same run for the direct comparison.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving humans were approved by Institutional Review Board of the National Institute of Advanced Industrial Science and Technology. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
TS: Investigation, Validation, Visualization, Writing – original draft, Writing – review & editing. AT: Investigation, Writing – review & editing. KK: Methodology, Writing – review & editing. DY: Methodology, Writing – review & editing. TT: Conceptualization, Funding acquisition, Investigation, Methodology, Project administration, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. The assessment conducted in the laboratory of TT was supported by Grant-in-Aid for Challenging Research (Exploratory) (Grant number 21 K19464).
Acknowledgments
We thank Ms. Natsu Nozaki for assisting with the MRI experiments.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Footnotes
References
Bauer, C. M., Zajac, L. E., Koo, B. B., Killiany, R. J., and Merabet, L. B. (2017). Age-related changes in structural connectivity are improved using subject-specific thresholding. J. Neurosci. Methods 288, 45–56. doi: 10.1016/j.jneumeth.2017.06.010
Braver, T. S., Cohen, J. D., Nystrom, L. E., Jonides, J., Smith, E. E., and Noll, D. C. (1997). A parametric study of prefrontal cortex involvement in human working memory. NeuroImage 5, 49–62. doi: 10.1006/nimg.1996.0247
Caeyenberghs, K., Metzler-Baddeley, C., Foley, S., and Jones, D. K. (2016). Dynamics of the human structural connectome underlying working memory training. J. Neurosci. 36, 4056–4066. doi: 10.1523/JNEUROSCI.1973-15.2016
Chen, X., and Wang, Y. J. (2018). Predicting resting-state functional connectivity with efficient structural connectivity. IEEE/CAA J. Autom. Sin. 5, 1079–1088. doi: 10.1109/Jas.2017.7510880
Cheng, H., Wang, Y., Sheng, J., Kronenberger, W. G., Mathews, V. P., Hummer, T. A., et al. (2012). Characteristics and variability of structural networks derived from diffusion tensor imaging. NeuroImage 61, 1153–1164. doi: 10.1016/j.neuroimage.2012.03.036
Cole, M. W., Ito, T., Cocuzza, C., and Sanchez-Romero, R. (2021). The functional relevance of task-state functional connectivity. J. Neurosci. 41, 2684–2702. doi: 10.1523/JNEUROSCI.1713-20.2021
Daneman, M. (1980). Individual differences in working memory and reading. J. Verbal Learn. Verbal Behav. 19, 450–466. doi: 10.1016/S0022-5371(80)90312-6
Darki, F., and Klingberg, T. (2015). The role of fronto-parietal and fronto-striatal networks in the development of working memory: a longitudinal study. Cereb. Cortex 25, 1587–1595. doi: 10.1093/cercor/bht352
Drew, P. J. (2019). Vascular and neural basis of the BOLD signal. Curr. Opin. Neurobiol. 58, 61–69. doi: 10.1016/j.conb.2019.06.004
Ekstrand, C., Neudorf, J., Kress, S., and Borowsky, R. (2020). Structural connectivity predicts functional activation during lexical and sublexical reading. NeuroImage 218:117008. doi: 10.1016/j.neuroimage.2020.117008
Gerchen, M. F., and Kirsch, P. (2017). Combining task-related activation and connectivity analysis of fMRI data reveals complex modulation of brain networks. Hum. Brain Mapp. 38, 5726–5739. doi: 10.1002/hbm.23762
Ginestet, C. E., and Simmons, A. (2011). Statistical parametric network analysis of functional connectivity dynamics during a working memory task. NeuroImage 55, 688–704. doi: 10.1016/j.neuroimage.2010.11.030
Hothorn, T., Hornik, K., Van De Wiel, M. A., and Zeileis, A. (2006). A Lego system for conditional inference. Am. Stat. 60, 257–263. doi: 10.1198/000313006x118430
Jaeggi, S. M., Seewer, R., Nirkko, A. C., Eckstein, D., Schroth, G., Groner, R., et al. (2003). Does excessive memory load attenuate activation in the prefrontal cortex? Load-dependent processing in single and dual tasks: functional magnetic resonance imaging study. NeuroImage 19, 210–225. doi: 10.1016/s1053-8119(03)00098-3
Kim, M., Bao, J., Liu, K., Park, B. Y., Park, H., Baik, J. Y., et al. (2021). A structural enriched functional network: an application to predict brain cognitive performance. Med. Image Anal. 71:102026. doi: 10.1016/j.media.2021.102026
Kirchner, W. K. (1958). Age differences in short-term retention of rapidly changing information. J. Exp. Psychol. 55, 352–358. doi: 10.1037/h0043688
Klingberg, T. (2006). Development of a superior frontal-intraparietal network for visuo-spatial working memory. Neuropsychologia 44, 2171–2177. doi: 10.1016/j.neuropsychologia.2005.11.019
Kourtidou, E., Kasselimis, D., Angelopoulou, G., Karavasilis, E., Velonakis, G., Kelekis, N., et al. (2022). Specific disruption of the ventral anterior temporo-frontal network reveals key implications for language comprehension and cognition. Commun. Biol. 5:1077. doi: 10.1038/s42003-022-03983-9
Kumari, V., Aasen, I., Taylor, P., Ffytche, D. H., Das, M., Barkataki, I., et al. (2006). Neural dysfunction and violence in schizophrenia: an fMRI investigation. Schizophr. Res. 84, 144–164. doi: 10.1016/j.schres.2006.02.017
Lamichhane, B., Westbrook, A., Cole, M. W., and Braver, T. S. (2020). Exploring brain-behavior relationships in the N-back task. NeuroImage 212:116683. doi: 10.1016/j.neuroimage.2020.116683
Lanciego, J. L., Luquin, N., and Obeso, J. A. (2012). Functional neuroanatomy of the basal ganglia. Cold Spring Harb. Perspect. Med. 2:a009621. doi: 10.1101/cshperspect.a009621
Latora, V., and Marchiori, M. (2001). Efficient behavior of small-world networks. Phys. Rev. Lett. 87:198701. doi: 10.1103/PhysRevLett.87.198701
Leech, R., and Sharp, D. J. (2014). The role of the posterior cingulate cortex in cognition and disease. Brain 137, 12–32. doi: 10.1093/brain/awt162
McKenna, B. S., Theilmann, R. J., Sutherland, A. N., and Eyler, L. T. (2015). Fusing functional MRI and diffusion tensor imaging measures of brain function and structure to predict working memory and processing speed performance among inter-episode bipolar patients. J. Int. Neuropsychol. Soc. 21, 330–341. doi: 10.1017/S1355617715000314
McLaren, D. G., Ries, M. L., Xu, G., and Johnson, S. C. (2012). A generalized form of context-dependent psychophysiological interactions (gPPI): a comparison to standard approaches. NeuroImage 61, 1277–1286. doi: 10.1016/j.neuroimage.2012.03.068
Mesulam, M. M., Wieneke, C., Hurley, R., Rademaker, A., Thompson, C. K., Weintraub, S., et al. (2013). Words and objects at the tip of the left temporal lobe in primary progressive aphasia. Brain 136, 601–618. doi: 10.1093/brain/aws336
Namkung, H., Kim, S. H., and Sawa, A. (2017). The insula: an underestimated brain area in clinical neuroscience, psychiatry, and neurology. Trends Neurosci. 40, 200–207. doi: 10.1016/j.tins.2017.02.002
Ning, Y., Fang, M., Zhang, Y., Feng, S., Feng, Z., Liu, X., et al. (2021). Attention performance correlated with white matter structural brain networks in shift work disorder. Front. Psych. 12:802830. doi: 10.3389/fpsyt.2021.802830
Olesen, P. J., Nagy, Z., Westerberg, H., and Klingberg, T. (2003). Combined analysis of DTI and fMRI data reveals a joint maturation of white and grey matter in a fronto-parietal network. Brain Res. Cogn. Brain Res. 18, 48–57. doi: 10.1016/j.cogbrainres.2003.09.003
Owen, A. M., McMillan, K. M., Laird, A. R., and Bullmore, E. (2005). N-back working memory paradigm: a meta-analysis of normative functional neuroimaging studies. Hum. Brain Mapp. 25, 46–59. doi: 10.1002/hbm.20131
Park, B. Y., Seo, J., Yi, J., and Park, H. (2015). Structural and functional brain connectivity of people with obesity and prediction of body mass index using connectivity. PLoS One 10:e0141376. doi: 10.1371/journal.pone.0141376
Picchioni, M., Matthiasson, P., Broome, M., Giampietro, V., Brammer, M., Mathes, B., et al. (2007). Medial temporal lobe activity at recognition increases with the duration of mnemonic delay during an object working memory task. Hum. Brain Mapp. 28, 1235–1250. doi: 10.1002/hbm.20357
Pongpipat, E. E., Kennedy, K. M., Foster, C. M., Boylan, M. A., and Rodrigue, K. M. (2021). Functional connectivity within and between n-Back modulated regions: an adult lifespan psychophysiological interaction investigation. Brain Connect. 11, 103–118. doi: 10.1089/brain.2020.0791
Power, J. D., Cohen, A. L., Nelson, S. M., Wig, G. S., Barnes, K. A., Church, J. A., et al. (2011). Functional network organization of the human brain. Neuron 72, 665–678. doi: 10.1016/j.neuron.2011.09.006
Rogers, B. P., Morgan, V. L., Newton, A. T., and Gore, J. C. (2007). Assessing functional connectivity in the human brain by fMRI. Magn. Reson. Imaging 25, 1347–1357. doi: 10.1016/j.mri.2007.03.007
Roman, F. J., Iturria-Medina, Y., Martinez, K., Karama, S., Burgaleta, M., Evans, A. C., et al. (2017). Enhanced structural connectivity within a brain sub-network supporting working memory and engagement processes after cognitive training. Neurobiol. Learn. Mem. 141, 33–43. doi: 10.1016/j.nlm.2017.03.010
Scholz, J., Klein, M. C., Behrens, T. E., and Johansen-Berg, H. (2009). Training induces changes in white-matter architecture. Nat. Neurosci. 12, 1370–1371. doi: 10.1038/nn.2412
Schulze, E. T., Geary, E. K., Susmaras, T. M., Paliga, J. T., Maki, P. M., and Little, D. M. (2011). Anatomical correlates of age-related working memory declines. J. Aging Res. 2011:606871. doi: 10.4061/2011/606871
Takeuchi, H., Sekiguchi, A., Taki, Y., Yokoyama, S., Yomogida, Y., Komuro, N., et al. (2010). Training of working memory impacts structural connectivity. J. Neurosci. 30, 3297–3303. doi: 10.1523/JNEUROSCI.4611-09.2010
Tournier, J. D., Smith, R., Raffelt, D., Tabbara, R., Dhollander, T., Pietsch, M., et al. (2019). MRtrix3: a fast, flexible and open software framework for medical image processing and visualisation. NeuroImage 202:116137. doi: 10.1016/j.neuroimage.2019.116137
Tsurugizawa, T., Taki, A., Zalesky, A., and Kasahara, K. (2023). Increased interhemispheric functional connectivity during non-dominant hand movement in right-handed subjects. iScience 26:107592. doi: 10.1016/j.isci.2023.107592
Tsurugizawa, T., Tamada, K., Ono, N., Karakawa, S., Kodama, Y., Debacker, C., et al. (2020). Awake functional MRI detects neural circuit dysfunction in a mouse model of autism. Sci. Adv. 6:eaav4520. doi: 10.1126/sciadv.aav4520
Wickham, H., Averick, M., Bryan, J., Chang, W., McGowan, L., François, R., et al. (2019). Welcome to the Tidyverse. J. Open Source Softw. 4:1686. doi: 10.21105/joss.01686
Wilhelm, O., Hildebrandt, A., and Oberauer, K. (2013). What is working memory capacity, and how can we measure it? Front. Psychol. 4:433. doi: 10.3389/fpsyg.2013.00433
Yaple, Z. A., Stevens, W. D., and Arsalidou, M. (2019). Meta-analyses of the n-back working memory task: fMRI evidence of age-related changes in prefrontal cortex involvement across the adult lifespan. NeuroImage 196, 16–31. doi: 10.1016/j.neuroimage.2019.03.074
Yaple, Z. A., Tolomeo, S., and Yu, R. (2021). Mapping working memory-specific dysfunction using a transdiagnostic approach. Neuroimage Clin. 31:102747. doi: 10.1016/j.nicl.2021.102747
Keywords: functional MRI, diffusion tensor imaging, N-back task, structural connectivity, local activation
Citation: Satake T, Taki A, Kasahara K, Yoshimaru D and Tsurugizawa T (2024) Comparison of local activation, functional connectivity, and structural connectivity in the N-back task. Front. Neurosci. 18:1337976. doi: 10.3389/fnins.2024.1337976
Edited by:
Mitsunari Abe, National Center of Neurology and Psychiatry (Japan), JapanReviewed by:
Masamichi J. Hayashi, Center for Information and Neural Networks (CiNet), JapanYuchuan Zhuang, AbbVie, United States
Copyright © 2024 Satake, Taki, Kasahara, Yoshimaru and Tsurugizawa. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Tomokazu Tsurugizawa, dC10c3VydWdpemF3YUBhaXN0LmdvLmpw