Europe PMC

This website requires cookies, and the limited processing of your personal data in order to function. By using the site you are agreeing to this as outlined in our privacy notice and cookie policy.

Abstract 


Recent neuroimaging studies have shown that working memory (WM) task difficulty can be decoded from patterns of brain activation in the WM network during preparation to perform those tasks. The inter-regional connectivity among the WM regions during task preparation has not yet been investigated. We examined this question using the graph modeling methods IMaGES and LOFS, applied to the previously published fMRI data of Manelis and Reder (2013). In that study, subjects performed 1-, 2-, and 3-back tasks. Each block of n-back was preceded by a preparation period and followed by a rest period. The analyses of task-related brain activity identified a network of 18 regions that increased in activation from 1- to 3-back (Increase network) and a network of 17 regions that decreased in activation from 1- to 3-back (Decrease network). The graph analyses revealed two types of connectivity sub-networks within the Increase and Decrease networks: "default" and "preparation-related." The "default" connectivity was present not only during task performance, but also during task preparation and during rest. We propose that this sub-network may serve as a core system that allows one to quickly activate cognitive, perceptual and motor systems in response to the relevant stimuli. The "preparation-related" connectivity was present during task preparation and task performance, but not at rest, and depended on the n-back condition. The role of this sub-network may be to pre-activate a connectivity "road map" in order to establish a top-down and bottom-up regulation of attention prior to performance on WM tasks.

Free full text 


Logo of fronthumneuroLink to Publisher's site
Front Hum Neurosci. 2014; 8: 593.
Published online 2014 Aug 5. https://doi.org/10.3389/fnhum.2014.00593
PMCID: PMC4122182
PMID: 25140143

Effective connectivity among the working memory regions during preparation for and during performance of the n-back task

Abstract

Recent neuroimaging studies have shown that working memory (WM) task difficulty can be decoded from patterns of brain activation in the WM network during preparation to perform those tasks. The inter-regional connectivity among the WM regions during task preparation has not yet been investigated. We examined this question using the graph modeling methods IMaGES and LOFS, applied to the previously published fMRI data of Manelis and Reder (2013). In that study, subjects performed 1-, 2-, and 3-back tasks. Each block of n-back was preceded by a preparation period and followed by a rest period. The analyses of task-related brain activity identified a network of 18 regions that increased in activation from 1- to 3-back (Increase network) and a network of 17 regions that decreased in activation from 1- to 3-back (Decrease network). The graph analyses revealed two types of connectivity sub-networks within the Increase and Decrease networks: “default” and “preparation-related.” The “default” connectivity was present not only during task performance, but also during task preparation and during rest. We propose that this sub-network may serve as a core system that allows one to quickly activate cognitive, perceptual and motor systems in response to the relevant stimuli. The “preparation-related” connectivity was present during task preparation and task performance, but not at rest, and depended on the n-back condition. The role of this sub-network may be to pre-activate a connectivity “road map” in order to establish a top-down and bottom-up regulation of attention prior to performance on WM tasks.

Keywords: working memory, task preparation, fMRI, working memory network, connectivity, graph modeling, IMaGES, LOFS

Introduction

Working memory (WM) is a system involved in on-line maintenance and manipulation of information (Baddeley and Hitch, 1974; Baddeley, 2010). It has a limited capacity (Cowan, 2001; Baddeley, 2010) and is critically important for learning, reasoning and decision-making (Mishkin and Manning, 1978; Petrides and Milner, 1982; Curtis and D'Esposito, 2003; Müller and Knight, 2006). Extensive neuroimaging research demonstrated that when a WM task becomes more difficult, frontal, parietal and striatal regions increase in activation (see Owen et al., 2005; Rottschy et al., 2012, for reviews), while frontal medial and posterior cingulate cortices decrease in activation (e.g., McKiernan et al., 2003; Esposito et al., 2006; Mayer et al., 2010; Manelis and Reder, 2013). Several studies have shown that inter-regional connectivity among the regions in the WM network also changes as a function of WM load (Honey et al., 2002; Wendelken et al., 2008; Kim et al., 2012; Ma et al., 2012; Dima et al., 2014). For example, increases in WM load were characterized by increases in connectivity among the frontal regions (right and left middle frontal gyrus (MFG), and left inferior frontal gyrus (IFG) and supplementary motor area; Honey et al., 2002), between the left parietal and frontal cortices (Ma et al., 2012), and increased right hemisphere dominance for verbal n-back (Dima et al., 2014), but increased left hemisphere dominance in face-matching n-back (Kim et al., 2012). Recent neuroimaging studies revealed that the modulation of brain activation occurs not only during performance on WM tasks of varying difficulty, but also during preparation to perform these tasks when no maintenance or manipulation requirements were yet imposed (Altamura et al., 2010; Manelis and Reder, 2013). For example, activation in the left inferior frontal gyrus (LIFG), anterior cingulate/paracingulate cortex and the left intraparietal sulcus (LIPS) linearly increased during 1- vs. 2- vs. 3-back task performance, but linearly decreased during preparation to perform the corresponding task (Manelis and Reder, 2013). Using a linear SVM [support vector machine (Vapnik, 1995)] classifier, Manelis and Reder demonstrated that the upcoming n-back condition (1- vs. 2- vs. 3-back) could be accurately decoded from the patterns of brain activation recorded in the WM network during preparation periods preceding 1-, 2-, or 3-back tasks. Not only did subjects with higher classification accuracies have more distinct neural representations for each upcoming difficulty level (a reason for why their classification accuracies were higher) they also displayed faster response times (RT) overall and smaller differences between 1-back and 3-back conditions. These results suggest that formation of more distinct neural representations during preparation periods helps subjects to be more efficient during task performance.

If the different preparation periods can be distinguished based on the activation patterns within the WM network, it is also possible that they can be distinguished based on the connectivity within that network. The inter-regional connectivity among the WM regions during task preparation has not yet been investigated in WM studies. In this study, we tested two alternative hypotheses: (1) Given that no information maintenance or manipulation is required during any condition's preparation periods, the connectivity among the WM regions during these preparation periods will resemble the connectivity at rest; vs. (2) The connectivity among the WM regions during preparation periods will resemble the connectivity during task performance. Formation of task-related connectivity before task performance (that is, during task preparation) may be beneficial because such advanced formation may (1) free some neural resources during task performance, and (2) help to activate and transfer the task-related rules and intentions from task preparation to the task itself.

The goal of this paper is to examine effective connectivity within the WM network during WM task preparation, WM task performance and during the rest periods that separate the blocks of the WM tasks. Many effective connectivity algorithms have been used to identify connectivity among regions within a network. Many of these approaches use a confirmatory method that require a priori model specification (e.g., DCM), which makes an exhaustive model search across large network of regions almost impossible (Hanson et al., 2007). Another limitation of many effective connectivity approaches is that they were not specifically designed for a multi-subject fMRI data processing and may produce false statistical dependencies by directly combining time series across subjects. A recently developed graphical analysis using the Independent Multiple sample Greedy Equivalence Search (IMaGES) algorithm (Ramsey et al., 2010) used in combination with the Linear non-gaussian Orientation, Fixed Structure (LOFS) algorithm (Ramsey et al., 2011) overcomes these limitations. Because both algorithms were specifically designed to estimate causal relationships based on simultaneous processing of multiple time series of multiple subjects, they do not produce artifacts related to concatenating multiple time series for the analyses (Ramsey et al., 2011). In simulation studies, the combined IMaGES + LOFS method showed very accurate performance on simulated data (Ramsey et al., 2011), while other effective connectivity algorithms had trouble with identification of both inter-regional connections and directions of the identified connections (Smith et al., 2011). For example, on simulated data that consisted of 50 variables used in Smith et al. (2011), IMaGES discovered over 95% of the connections (see Section Effective connectivity analyses for detailed description of how connections were discovered) and LOFS correctly oriented over 80% of these connections (Ramsey et al., 2011).

Based on the fact that using the combination of IMaGES and LOFS produced the most accurate performance among all effective connectivity algorithms in simulation studies, that IMaGES and LOFS were specifically designed to deal with the multi-subject fMRI data and that these methods were successfully used to examine effective connectivity within the language (Boukrina and Graves, 2013; Boukrina et al., 2014) and social brain (Hanson et al., 2013) networks, we decided to use these algorithms in our study. In this paper, we applied IMaGES (to search for connections) and LOFS (to orient the connections discovered by IMaGES) algorithms on Manelis and Reder's (2013) fMRI data to examine effective connectivity within the WM network of regions that increased in activation from 1- to 3-back (Increase network) and within the WM network of regions that decreased in activation from 1- to 3-back (Decrease network). The connectivity among 18 regions comprising the Increase network and, separately, among 17 regions comprising the Decrease network was examined during the 1-, 2-, and 3-back tasks, during preparation for these tasks and during rest periods separating the n-back blocks. We focused on identifying two distinct groups of connections—those that were (1) common for all preparation periods, all n-back blocks and rest periods (“default connections”), and those that were (2) common for task preparation and the corresponding n-back conditions, but were not observed during rest periods (“preparation-related connections”).

Materials and methods

Subjects and task design

The details pertaining to the study participants, the n-back task and the fMRI methods were published in Manelis and Reder (2013). In short, the dataset consisted of 16 subjects (mean age = 24, 11 females), all of whom were right-handed, native speakers of English, with normal or corrected-to-normal vision. They were treated in accordance with the Carnegie Mellon University Institutional Review Board guidelines.

The subjects were scanned while doing the n-back task at the three levels of difficulty (1-, 2-, and 3-back) with words selected from the MRC Psycholinguistic Database as stimuli. The words were between 4 and 7 letters long and were repeated within, but not between, the blocks. The words were separated with a jittered interval of 2–8 s. The 12-trial blocks of 1-, 2-, and 3-back (10 blocks in each condition) were presented in random order. The duration of a trial was limited to 4 s. Each block was preceded by an 8-s instruction screen, which informed subjects about whether the upcoming block would be 1-, 2-, or 3-back by displaying “1-back,” “2-back,” or “3-back.” A 10–12-s rest period followed each block.

fMRI acquisition, preprocessing, and GLM analyses

As described in Manelis and Reder (2013), the fMRI data were collected using a Siemens 3T Verio MR system. We acquired a high-resolution structural image (0.8 × 0.8 × 0.8 mm) using MPRAGE (TR = 1800 ms, TE = 2.22 ms, FOV = 205, FA = 9°, number of slices = 256), functional data using a gradient-echo echo-planar sequence (TR = 2000 ms, TE = 30 ms, FOV = 205, FA = 79°, 36 slices, 3.2 × 3.2 × 3.2 mm), and field maps with the same resolution as the BOLD images using a gradient-echo sequence (TR = 394 ms, FA = 60°, TE = 5.1 and 7.56 ms).

The fMRI data were preprocessed and analyzed using FSL 4.1.7 (www.fmrib.ox.ac.uk/fsl). Preprocessing included non-linear noise reduction performed using SUSAN (http://fsl.fmrib.ox.ac.uk/fsl/fslwiki/SUSAN), motion correction with MCFLIRT (Jenkinson et al., 2002), fieldmap-based EPI unwarping using PRELUDE+FUGUE (Jenkinson, 2003), non-brain removal using BET (Smith, 2002), spatial smoothing using a Gaussian kernel of FWHM 6 mm, grand-mean intensity normalization of the entire 4D dataset by a single multiplicative factor, high-pass temporal filtering (Gaussian-weighted least-squares straight line fitting, with sigma = 50.0 s). The Probabilistic Independent Component Analysis (ICA; Beckmann and Smith, 2004), implemented using FSL's Multivariate Exploratory Linear Decomposition into Independent Components (http://www.fmrib.ox.ac.uk/fsl/melodic/index.html), served to identify “noise” components (Tohka et al., 2008; Kelly et al., 2010) that were then removed using the fsl_regfilt script.

The preprocessed data were used in the GLM analysis with 6 regressors (1-, 2-, and 3-back instruction periods, and 1-, 2-, and 3-back task performance blocks) that examined two linear trends in brain activation: 1-back<2-back<3-back and 1-back>2-back>3-back. Co-registration of BOLD images with the MNI152_T1_2 mm template was carried out using FLIRT (Jenkinson and Smith, 2001; Jenkinson et al., 2002). As described in Manelis and Reder (2013), a group analysis was conducted using the Randomise v2.1 tool (http://www.fmrib.ox.ac.uk/fsl/randomise/index.html) with the whole brain as a mask, 5 mm smoothing, 5000 permutations and correction for multiple comparisons at the voxel-wise FWE-controlled threshold p < 0.05). This analysis revealed a network of 18 regions that linearly increased activation from 1- to 3-back (Increase network; see Table Table1)1) and a network of 17 regions that linearly decreased activation from 1- to 3-back (Decrease network; see Table Table1)1) that were used in the effective connectivity analyses described below.

Table 1

Parametric changes in brain activation as a function of working memory load during task performance.

RegionAbbreviationn-voxelsZ-maxxyz
PARAMETRIC INCREASES IN ACTIVATION (1-back<2-back<3-back) DURING TASK PERFORMANCE
LAnterior cingulate c./paracingulate g.ACC8189.99−22242
RInferior parietal s.RIPS5038.2940−4644
LOrbitofrontal c.LOFc4479.66−3224−6
LInferior parietal s.LIPS4078.0−34−5240
LSuperior frontal g.LSFG3887.4−28862
RFrontal poleRFP3708.22424822
RPrecuneusRprecun3147.992−6248
LFrontal poleLFP3137.58−385610
LInferior frontal g.LIFG3008.21−40826
LMiddle frontal g.LMFG2777.66−483622
RSuperior frontal g.RSFG2626.87281260
RInsular c.Rins1778.883424−2
LBasal gangliaLbas1607.2−16014
RLateral occipital c., superiorRLOCs977.132−6646
LThalamusLthal526.49−8−188
RCerebellumRcereb275.8336−70−28
LLateral occipital c., superiorLLOCs266.42−16−7452
RMiddle frontal g.RMFG105.92443438
PARAMETRIC DECREASES IN ACTIVATION (3-back<2-back<1-back) DURING TASK PERFORMANCE
LPlanum polareLPlanPol159910.4−42−16−8
RParietal operculumRParOperc13058.6748−3022
LAnterior cingulate/Medial frontal c.LFmed9688.22−254−4
BJuxtapositional lobule c.BJuxt9169.02−848
LPosterior cingulate g., posteriorLPCCp5918.29−10−5028
RRight pre−central g.RPrec3127.136−1848
RPlanum polareRPlanPol2287.14420−18
LPosterior cingulate g., anteriorLPCCa2008.47−14−3038
LPost−central/pre−central g.LPostPrec648.2−36−1840
RLateral occipital c., inferiorRLOCi606.5254−708
RPost−central g., anteriorRPosta316.4156−1648
RTemporal poleRTP176.554024−26
BSubcallosal c.BSub135.9016−12
RFrontal poleRFP135.853636−16
RPost−central g., middle partRPostm12630−3270
LOccipital poleLOccipP115.69−26−9424
RPost−central g., posteriorRPostp105.7126−3656

The images were thresholded at voxel-wise FWE-corrected p < 0.05. The abbreviations for all regions are in the third column. L, left; R, right; B, bilateral; g., gyrus; s., sulcus; c., cortex.

Effective connectivity analyses

The effective connectivity analyses were conducted using graphical causal modeling using IMaGES (the Independent Multiple sample Greedy Equivalence Search) and LOFS (Linear non-gaussian Orientation, Fixed Structure) algorithms (Ramsey et al., 2010, 2011, 2014; Mumford and Ramsey, 2014) implemented using the TETRAD IV (version 5.0.0-1; http://www.phil.cmu.edu/projects/tetrad) software. Separate connectivity analyses were conducted for the regions within the Increase network and the regions within the Decrease network. First, we extracted the time series from each of the ROIs identified in the Manelis and Reder (2013) study for each subject. All those ROIs have already been corrected for multiple comparisons in the original GLM that identified those ROIs. Each of these ROIs is a node in the network whose connectivity we examine in the analyses described below. Each of the three n-back conditions was associated with 10 preparation periods lasting for 4 TRs each, which provided us with a total of 40 data points per preparation condition per subject per ROI. The length of each block of n-back was at least 21 TRs (but could be longer given the self-paced nature of the task), which provided us with at least 210 data points per n-back condition per subject per ROI. There were 30 rest periods in the experiment. The length of each rest periods was 5–6 TRs. Only the last 3–4 TRs were included in the data analyses, which provided us with at least 90 data points per subject per ROI (3 TRs × 30 rest periods).

Second, the extracted time series combined across subjects and ROIs were used as input to the IMaGES algorithm with increasing penalty discounts (Ramsey et al., 2010). IMaGES is a Bayesian search algorithm that produces a Markov equivalence class of models that have the same structure (the same connections between the nodes without considering the direction of those connections). For each set of regions, IMaGES starts with an empty graph. It tests all possible models with one connection and selects one model with the highest Bayesian Information Criterion (BIC) score averaged across several datasets. The algorithm continues to add connections between nodes until further connections do not improve the BIC score. At that point, the process is reversed, removing connections from the model, one at a time until the BIC score can no longer be improved. In some cases, Greedy Equivalence Search can produce the graphs where the three nodes (e.g., three brain regions) are connected to each other (i.e., “triangulation”). Given that we are trying to estimate causal relationships among latent variables, such “triangulation” can lead to false conclusions (Ramsey et al., 2010). “Triangulation” and the inherent possibility of spurious causal connections can be avoided by increasing the penalty function in the BIC score (Ramsey et al., 2010). In this study, we used the option “find first non-triangular” to search for a graph that does not contain “triangulation” [directed acyclic graph (DAG)].

Third, after IMaGES identified a DAG for the set of regions, the DAG is fed to the LOFS algorithm (Ramsey et al., 2011). LOFS determines the orientation (direction) of each connection by exploiting the fact that the residuals of any incorrect linear model will be more Gaussian than the residuals of the correct model with independent non-Gaussian sources of error (Ramsey et al., 2011). In our study, the degree of non-Gaussianity was estimated using the Anderson–Darling score (Anderson and Darling, 1952). After LOFS oriented the connection, each graph consisted of nodes and arrows (or edges, or connections) that connect some of those nodes thus depicting causal relationships between them. For example, the arrow in the LMFG→LIFG expression shows that the changes in LMFG activation influence the changes in LIFG activation. Finally, after the connections were detected and oriented, the outcome of the LOFS algorithm is submitted to a SEM (structural equation modeling) estimator that estimates model goodness of fit to each set of data by estimating the values of parameters for a SEM parametric model with a regression optimizer. Accurate regression estimates presuppose that the input parametric model be a DAG, and its associated statistics are based on a linear, Gaussian model.

The effective connectivity within Increase and Decrease networks was assessed for each of the three n-back conditions, for each of the three corresponding preparation periods and also at rest. The rest periods were analyzed starting 4 s after the end of a preceding n-back block. The effective connectivity for the seven conditions (rest, 1, 2, 3-back preparation periods and 1, 2, 3-back task performance periods) were compared to identify the connections that were common for all these conditions and those that were common for the corresponding preparation and n-back conditions but were absent at rest.

Results

Behavioral and neuroimaging activation analyses

The details of the behavioral analysis and the univariate analysis of neuroimaging data are described in Manelis and Reder (2013). Consistent with previous work (e.g., Braver et al., 1997; Nystrom et al., 2000), they found that, as the n-back task difficulty increased, performance slowed [F(2, 30) = 46.4, p < 0.001] and became less accurate [F(2, 30) = 22.5, p < 0.001]. Table Table11 and Figure Figure11 report regions (and their corresponding abbreviations) that increased in activation as the task difficulty increased, and also those regions that decreased in activation as the task difficulty increased.

An external file that holds a picture, illustration, etc.
Object name is fnhum-08-00593-g0001.jpg

Increase (1-back<2-back<3-back, shown in red) and Decrease (1-back>2-back>3-back, shown in blue) networks identified in the GLM analysis. The abbreviations for the brain regions are described in Table Table1.1. Left hemisphere is on the right and right hemisphere is on the left.

Effective connectivity analysis

The results of IMaGES and LOFS analyses are presented in Tables Tables22 (for Increase network) and and33 (for Decrease network). The tables report mean and standard deviations of SEM coefficients across subjects. The connections between pairs of regions are in the left column. The directions of the connections are shown with the directed arrows “→.” For example, the sentence A→B means that changes in activation in the region A cause changes in activation in the region B. Sometimes, IMaGES would find that a certain connection existed in several conditions; however, LOFS determined that the direction of this connectivity varied among the specific conditions. In cases such as these (i.e., when the direction of the connection was different from that shown in the left column of Tables Tables2,2, ,3),3), we placed the arrow “←” left from the SEM coefficient. Although Tables Tables2,2, ,33 report all connections and their directions identified by IMaGES and LOFS, the discussion will focus only on those connections that were common for all seven conditions and on the connections that were common for the corresponding task preparation and task performance periods (e.g., 1-back preparation and 1-back performance) but were absent at rest. Figure Figure22 (Increase network) and Figure Figure33 (Decrease network) also focus only on these two sets of the findings.

Table 2

Inter-regional connectivity within the Increase network.

ConnectionsIntro11-backIntro22-backIntro33-backRest
RFP→RMFG0.72 (0.3)0.94 (0.3)← 0.59 (0.2)0.96 (0.3)← 0.22 (0.2)0.97 (0.3)← 0.21 (0.1)
LIPS→LIFG← 0.28 (0.2)0.42 (0.2)← 0.45 (0.2)0.52 (0.2)0.7 (0.3)0.34 (0.2)0.32 (0.3)
LMFG→LIFG0.45 (0.2)0.4 (0.2)0.69 (0.2)← 0.51 (0.3)← 1.08 (0.4)0.31 (0.1)0.34 (0.1)
LFP→RFP0.43 (0.3)← 0.82 (0.2)← 0.81 (0.2)← 0.82 (0.2)0.34 (0.2)0.49 (0.1)0.36 (0.2)
LFP→LMFG← 0.93 (0.2)0.69 (0.2)0.33 (0.2)0.26 (0.2)← 0.89 (0.2)0.38 (0.2)0.41 (0.2)
ACC→LOFc← 0.89 (0.2)0.76 (0.2)0.51 (0.2)0.77 (0.2)0.78 (0.2)← 0.86 (0.2)0.43 (0.2)
LOFc→Rins0.86 (0.2)0.74 (0.1)← 0.41 (0.2)0.76 (0.1)0.66 (0.2)0.78 (0.1)← 0.47 (0.1)
Lbas→Lthal1.07 (0.3)1.08 (0.3)← 0.33 (0.2)0.97 (0.2)1.03 (0.4)← 0.29 (0.1)← 0.48 (0.2)
Rprecun→LLOCs← 0.59 (0.3)0.68 (0.2)0.75 (0.3)0.71 (0.2)0.5 (0.3)← 0.61 (0.2)0.70 (0.2)
RIPS→LIPS0.29 (0.2)← 0.61 (0.2)0.51 (0.1)0.82 (0.2)0.78 (0.2)0.65 (0.1)0.81 (0.2)
RSFG→LSFG0.74 (0.2)0.4 (0.2)← 0.91 (0.2)0.5 (0.2)0.67 (0.2)← 0.61 (0.1)← 0.85 (0.1)
LIFG→ACC0.33 (0.3)0.61 (0.2)0.39 (0.2)0.39 (0.2)0.24 (0.1)0.19 (0.2)
LSFG→LMFG0.99 (0.4)← 0.27 (0.1)0.66 (0.3)0.34 (0.2)0.62 (0.1)0.6 (0.2)
RIPS→RLOCs0.94 (0.2)0.96 (0.3)0.96 (0.3)0.78 (0.3)0.92 (0.2)0.73 (0.3)
RIPS→RFP0.28 (0.3)← 0.24 (0.1)0.39 (0.2)← 0.64 (0.2)0.22 (0.2)
LIPS-LSFG← 0.17 (0.2)0.17 (0.1)0.39 (0.2)0.74 (0.2)
RSFG→RFP0.32 (0.3)← 0.33 (0.2)0.21 (0.2)
LLOCs→RLOCs0.25 (0.2)← 0.39 (0.3)0.27 (0.2)
RSFG→ACC0.75 (0.2)
RFP→ACC← 0.00 (0.3)0.34 (0.1)0.28 (0.1)← 0.5 (0.2)
RSFG→RIPS0.18 (0.2)0.69 (0.2)← 0.85 (0.2)0.63 (0.4)
ACC→Lbas0.66 (0.1)0.37 (0.2)0.59 (0.2)0.66 (0.3)0.38 (0.2)
LSFG→ACC0.9 (0.3)0.35 (0.2)← 0.76 (0.2)
LIPS→LLOCs← 0.46 (0.1)0.55 (0.5)0.82 (0.4)
RIPS→Rprecun1.19 (0.3)← 0.17 (0.2)0.48 (0.3)
Rcereb→Rprecun0.21 (0.3)0.13 (0.4)
RLOCs→Rprecun0.56 (0.4)0.86 (0.4)
LFP→LIPS0.16 (0.1)
LFP→Rprecun0.82 (0.3)
LMFG→RMFG0.34 (0.3)
Rcereb→LIPS0.13 (0.1)
RIPS→Rins0.19 (0.2)
RLOCs→LIPS0.16 (0.1)
RMFG→LSFG0.07 (0.2)
Rprecun→Rcereb0.74 (0.3)
Rprecun→RFP0.08 (0.2)

The arrow “←” shows the connections whose direction was opposite to one described in the left column titled “Connections.” Standard deviations of mean SEM coefficients across the subjects are in parentheses. The connections formed during the preparation periods and carried over to the corresponding task performance periods, but absent during rest, are shown in bold.

Table 3

Inter-regional connectivity within the Decrease network.

ConnectionsIntro11-backIntro22-backIntro33-backRest
LPCCa→LPCCp0.59 (0.4)1.1 (0.3← 0.17 (0.1)0.55 (0.3)0.83 (0.5)1.13 (0.2)1.04 (0.3)
LPCCp→LFmed← 0.48 (0.4)0.64 (0.2)0.52 (0.3)← 0.54 (0.2)0.5 (0.3)0.73 (0.2)0.73 (0.1)
RFP→LFmed← 0.59 (0.2)0.25 (0.2)0.35 (0.3)0.73 (0.2)0.42 (0.4)← 0.58 (0.2)← 0.48 (0.2)
BJuxt→LPCCa← 1.3 (0.4)← 0.82 (0.2)0.52 (0.1)0.63 (0.1)0.61 (0.1)0.63 (0.1)0.55 (0.1)
BJuxt→LPostPrec0.61 (0.2)0.6 (0.2)0.61 (0.2)0.59 (0.2)← 1.06 (0.3)0.62 (0.2)← 0.55 (0.3)
BJuxt→RPrec0.73 (0.2)← 0.39 (0.1)0.74 (0.3)0.7 (0.2)0.26 (0.2)0.69 (0.2)0.37 (0.2)
RLOCi→LOccipP0.69 (0.2)0.71 (0.3)0.7 (0.3)0.69 (0.2)0.64 (0.4)← 0.31 (0.2)0.67 (0.2)
RParOperc→LPlanPol← 0.72 (0.2)0.63 (0.2)← 1.04 (0.2)← 0.67 (0.1)0.47 (0.2)0.51 (0.1)0.74 (0.1)
RPrec→RPosta0.74 (0.3)0.77 (0.2)0.78 (0.2)0.76 (0.2)← 0.35 (0.2)0.76 (0.2)0.71 (0.2)
RPrec→RPostp0.84 (0.2)0.46 (0.3)0.81 (0.2)0.72 (0.2)← 0.38 (0.2)0.76 (0.2)← 0.46 (0.1)
BJuxt→RParOperc0.31 (0.2)0.76 (0.1)0.32 (0.1)← 0.52 (0.1)
RPostp→RPostm1.17 (0.3)← 0.31 (0.2)1.1 (0.4)1.09 (0.4)
RParOperc→RPlanPol0.86 (0.1)0.86 (0.2)← 0.55 (0.2)
RTP→RPlanPol0.18 (0.2)0.34 (0.2)← 0.72 (0.5)
LPCCa→RLOCi0.95 (0.4)0.80 (0.2)
RPostp→LPCCa0.1 (0.09)
RParOperc→RLOCi0.64 (0.3)0.58 (0.4)0.64 (0.2)0.48 (0.2)
RPlanPol→LPlanPola← 0.9 (0.2)0.14 (0.1)0.29 (0.2)0.15 (0.1)
BJuxt→LPlanPol0.47 (0.1)0.24 (0.1)0.29 (0.1)
RPrec→RPostm1.12 (0.3)1.02 (0.4)1.05 (0.4)
RTP→LFmed0.11 (0.1)0.17 (0.2)
RFP→RTP0.85 (0.6)
BJuxt→RTP0.7 (0.5)
RLOCi→BJuxt0.06 (0.1)

The arrow “←” shows the connections whose direction was opposite to one described in the left column titled “Connections.” Standard deviations of mean SEM coefficients across the subjects are in parentheses. The connections formed during the preparation periods and carried over to the corresponding task performance periods, but absent during rest, are shown in bold.

An external file that holds a picture, illustration, etc.
Object name is fnhum-08-00593-g0002.jpg

Connectivity within the Increase network. Black lines depict the connections common for 1-, 2-, 3-back preparation, 1-, 2-, 3-back task performance and also rest periods. The red line represents the connection common for 1-back preparation and task, but not rest periods. Blue lines represent the connections common for 2-back preparation and task, but not rest periods. Green lines represent the connections common for 3-back preparation and task, but not rest periods. Please note that the direction of each arrow represents the directionality for the majority of connections (not all of them). The lines without arrows represent the connections whose directions differed for preparation and task performance conditions.

An external file that holds a picture, illustration, etc.
Object name is fnhum-08-00593-g0003.jpg

Connectivity within the Decrease network. Black lines depict the connections common for 1-, 2-, 3-back preparation, 1-, 2-, 3-back task performance and also rest periods. The red line represents the connection common for 1-back preparation and task, but not rest, periods. Blue lines represent the connections common for 2-back preparation and task, but not rest, periods. The green line represents the connection common for 3-back preparation and task, but not rest, periods. Please note that the direction of each arrow represents the directionality for the majority of connections (not all of them). The lines without arrows represent the connections whose directions differed for preparation and task performance conditions.

Connections common across all conditions in the experiment including rest: Increase network

Within the Increase network, 11 connections were present during rest, all three preparatory periods and all three task performance conditions. While the direction of those connections varied depending on the condition, in general, those connections formed a network linking frontal and parietal regions as well as the right and left hemispheres: RMFG-RFP, RFP-LFP, LFP-LMFG, LMFG-LIFG, LIFG-LIPS, LIPS-RIPS, Rins-LOFc, LOFc-ACC, Lbas-Lthal, LLOCs-Rprecun, RSFG-LSFG (black arrows in Figure Figure2).2). Given the verbal nature of the task, it is not surprising that there were more left-sided connections than right-sided ones.

Connections that existed during preparation periods and corresponding task conditions, but were absent at rest: Increase network

The graph modeling analysis revealed that some connections were present during corresponding preparation and task periods, but not at rest (colored arrows and lines in Figure Figure2).2). The number of common connections increased for more difficult tasks. One common connection (RFP-ACC) was found for 1-back preparation and 1-back task performance. Two common connections (RSFG-RIPS and ACC-Lbas) were found for 2-back preparation and 2-back task performance. Five common connections (RFP-ACC, ACC-Lbas, ACC-LSFG, LIPS-LLOCs, and RIPS-Rprecun) were found for 3-back preparation and 3-back task performance periods.

Connections common across all conditions in the experiment including rest: Decrease network

Within the Decrease network, 10 connections (BJuxt-RPrec, Rprec-RPosta, Rprec-Rpostp, BJuxt-LPostPrec, BJuxt-LPCCa, LPCCa-LPCCp, LPCCp-LFmed, RFP-LFmed, RParOperc-LPlanPol, RLOCi-LOccip) existed throughout the experiment during rest, all preparatory and all task performance conditions (Table (Table3,3, Figure Figure3).3). Those connections linked the right and left hemispheres (e.g., RParOperc-LPlanPol, RLOCi-LOccip) as well as the frontal, parietal and occipital regions, often through the structures located along the midline (e.g., BJuxt-RPrec, BJuxt-LPostPrec, BJuxt-LPCCa, LPCCa-LPCCp, LPCCp-LFmed). BJuxt was a point of origin for several connections as the changes in activation in this region often caused changes in activation in LpostPrec, LPCCa, and RPrec. RPrec, in turn, often caused changes in activation in RPosta and RPostp. In contrast, the LFmed was a point of convergence as its activation was often caused by changes in activation in RFP and LPPCp.

Connections that existed during preparation periods and corresponding task conditions, but were absent at rest: Decrease network

In the Decrease network, there was a common LPlanPolar-RPlanPolar connection for 1-back preparation and 1-back task performance. In the 2-back condition, there were common RPrec→RPostm and RparOperc→RLOCi connections. In the 3-back condition, there was a common BJuxt→LPlanPolar connection.

Connections observed during preparation periods and at rest, but absent during task performance

In the Increase network, there were no connections that were common for the preparation periods and rest that were not also present during task performance. In the Decrease network, RTP-RPlanPol connectivity was found during preparation for the 1- and 3-back tasks and also at rest, but not during n-back task performance. In the preparation conditions, activation in RTP caused activation in RPlanPol, while at rest, activation in RPlanPol caused activation in RTP. LPCCa→RLOCi connectivity was observed during preparation for the 3-back task and at rest but not during other conditions.

Discussion

This is the first study to examine effective connectivity within the WM network during preparation to perform the n-back task, during performance on these tasks, and during rest periods that separate the n-back blocks. We applied recently developed graph modeling methods IMaGES and LOFS to Manelis and Reder's (2013) fMRI data in order to examine connectivity within the Increase and Decrease networks. The regions within the Increase network increased in activation with an increase in WM load and were consistent with those previously reported in other WM studies (see Owen et al., 2005, for a review). The regions within the Decrease network decreased in activation when WM load or attentional demands increased (e.g., McKiernan et al., 2003; Esposito et al., 2006; Mayer et al., 2010). These included the regions associated with the default mode network (e.g., Raichle et al., 2001; Greicius et al., 2003; Buckner et al., 2008). The results of the graph analyses revealed two types of connectivity sub-networks within the Increase and Decrease networks: “default” and “preparation-related.”

“Default” sub-network within the increase and decrease networks

We called one sub-network “default” connectivity because those connections existed regardless of whether a subject was preparing for the task, was performing the task or was at rest. “Default” connections were found in both Increase (RMFG-RFP, RFP-LFP, LFP-LMFG, LMFG-LIFG, LIFG-LIPS, LIPS-RIPS, Rins-LOFc, LOFc-ACC, Lbas-Lthal, LLOCs-Rprecun, RSFG-LSFG) and Decrease (BJuxt-RPrec, Rprec-RPosta, Rprec-Rpostp, BJuxt-LPostPrec, BJuxt-LPCCa, LPCCa-LPCCp, LPCCp-LFmed, RFP-LFmed, RParOperc-LPlanPol, RLOCi-LOccip) networks. Those connections were present across all conditions despite the fact that activation in these regions depended on WM load in the task. While it may seem surprising that the patterns of connectivity did not resemble the patterns of activation and did not change in a linear manner, it is possible that stable connectivity is beneficial in those situations in which activity in brain regions is determined by external requirements that unpredictably change during the experiment.

The “default” connections linked 17 of the 18 regions in the Increase network and 14 of the 17 regions in the Decrease network, thus connecting most of the regions within the WM network to allow integration and quick propagation of incoming information when the task requirements changed. Analogous to the default mode network (e.g., Raichle et al., 2001; Greicius et al., 2003; Buckner et al., 2008), “default” connectivity may play a fundamental role in monitoring the environment, cognitive requirements, and motor responses. PFC and IPS belong to the attentional network (Corbetta, 1998) and are thought to also be involved in task preparation (Brass and von Cramon, 2004). Connectivity among these regions modulates attentional control (Wang et al., 2010), while disconnection among these regions results in impaired attention (Neufang et al., 2011). “Default” connectivity between PFC and IPS regions may help to maintain alertness throughout the experiment thereby enabling responses to stimuli in a timely manner. “Default” connectivity across primary, sensory and motor regions may help individuals to respond “as quickly as possible” when prompted by a task.

Honey et al. (2002) found connectivity between PFC (LIFG) and posterior parietal cortex during the 1- and 2-back conditions. In our study, we found LIFG-LIPS connectivity not only for the 1, 2, and 3-back conditions, but also for all the preparatory and rest conditions. Given the verbal nature of our task and the fact that this connection may be important for mediating articulatory rehearsal (Honey et al., 2002), the presence of this connection in all conditions may have been necessary to facilitate the input of new verbal information during the task. Ma et al. (2012) reported that WM load modulates connectivity within the fronto-parietal network. In contrast to their finding that WM load modulates left parietal→LIFG and right parietal→left parietal connectivity, we found that these connections were present even at rest, and were not modulated by either actual or expected task difficulty. Such differences can be explain by several factors that include differences in the tasks, differences between DCM and IMaGES in modeling effective connectivity and the fact that we examined connectivity not only at task, but also during preparatory and rest periods.

The Decrease network includes several regions that belong to the default mode network (e.g., PCC, Fmed). The function of the default mode network is still debated. For example, one group of researchers argues that its activation is related to task-unrelated thoughts (e.g., McKiernan et al., 2006), another that it reflects fundamental functional organization (e.g., Raichle and Snyder, 2007; Vincent et al., 2007), and still others argue that it reflects recent experiences (e.g., Albert et al., 2009; Hasson et al., 2009; Tambini et al., 2010). Connectivity among the regions that decrease in activation during task performance was often examined at rest rather than during task performance; however, some studies have shown that these regions might be equally important for cognitive functioning as those regions that increase in activation during a task (e.g., Sambataro et al., 2010; Yakushev et al., 2013). Our study supports this latter idea by showing that both Increase and Decrease networks contain “default” connections whose role may be to enable an immediate response to changes in the environment.

“Preparation-related” sub-network within the increase and decrease networks

Connectivity in the “preparation-related” sub-network was formed during task preparation and carried over to task performance, but not rest periods. The number of such connections depended on task difficulty. Two connections (one in the Increase network [RFP-ACC] and one in the Decrease network [RPlanPol-LPlanPol]) comprised the “preparation-related” connectivity sub-network during the 1-back task. Four connections (two in the Increase network [RSFG-RIPS and ACC-Lbas] and two in the Decrease network [RPrec→RPostm and RparOperc→RLOCi]) comprised the “preparation-related” sub-network during the 2-back task. Six connections (five in the Increase network [RFP-ACC, ACC-Lbas, ACC-LSFG, LIPS-LLOCs, and RIPS-Rprecun] and one in the Decrease network [BJuxt→LPlanPolar]) comprised “preparation-related” sub-network during 3-back.

If the “preparation-related” connectivity within the WM load networks reflect only general task preparation and rule activation, then the number of connections formed during preparation should not depend on the anticipated task difficulty. However, in this study, the number of “preparation-related” connections increased with an increase in the level of expected task difficulty, suggesting that preparation is specific to the level of cognitive demand in the anticipated task. Because different n-back conditions likely require the use of different cognitive strategies, different preparatory conditions were associated with the pre-formation of different connections. Many of those connections involved the ACC, which plays a role in conflict monitoring (e.g., Botvinick et al., 2001, 2004) and in anticipation of conflict monitoring (Sohn et al., 2007). The 1-back condition is the least demanding. It was associated with only a weak (based on the SEM coefficient) ACC→RFP connectivity during preparation and RFP→ACC connectivity at task. Considering that the RFP is associated with time-based prospective memory tasks for both words and pictures (Volle et al., 2011) and with visuospatial prospective memory (Costa et al., 2013), the ACC→RFP connectivity may reflect the regulation processes that change future intentions based on the current state of conflict (e.g., response error) detected by the ACC. The RFP→ACC connectivity, in contrast, may reflect the top-down process that regulates the perception of conflict based on the adjusted prospective goals and intentions.

The 2- and 3-back conditions are more difficult than the 1-back condition and likely involve different strategies to ensure optimal performance. Preparatory processes presumably involve inhibiting the n-back rules specific to the previous n-back block, activating the rules specific for the upcoming block, and establishing some level of cognitive control before a given n-back block starts. Both preparation and task performance in the 2-back and 3-back conditions involved the ACC→Lbas connectivity: the increase in the ACC activation caused the increase in the Lbas activation. The ACC is involved in conflict monitoring (Botvinick et al., 2004). The conflict may arise from the need to use new rules for the upcoming block, from the emotional reaction to the objective or subjective task difficulty and from the perception of response errors. These information may be task-irrelevant and may interfere with task performance. After the ACC detects such information, the activation is spread to the Lbas that inhibits irrelevant information (Yehene et al., 2008). Thus, the ACC→Lbas connectivity may serve for detection and inhibition of irrelevant to the current task information to ensure optimal task performance. During the 3-back task performance, activation in the ACC causes not only a change in Lbas activation but also changes in LSFG and RFP activation, thus spreading activation to frontal regions involved in on-line monitoring and manipulation of information [LSFG (du Boisgueheneuc et al., 2006)] as well as in prospective memory (RFP; e.g., Volle et al., 2011; Costa et al., 2013). In contrast, during preparation periods preceding 3-back trials, activation in the LSFG and RFP caused changes in ACC activation, thus preparing this region for optimal functioning during the demanding 3-back task.

“Preparation-related” connectivity within the parietal cortex (RIPS-Rprecun) and between the parietal and occipital cortices (LIPS-LLOCs) was unique for 3-back preparation and task performance. Together with the “default” RIPS-LIPS and RPrecun-LLOCs connections, they formed a fully connected parieto-occipital network whose function may have been to prepare to integrate and to integrate information about the spatial position of each stimulus in the n-back sequence. This information could then be transferred to the LIFG through the “default” LIPS→LIFG connection for further processing (i.e., evaluation, updating and manipulation) in the frontal cortex. Taken together, the role of the “preparation-related” connectivity sub-network is to establish a top-down and bottom-up regulation of attention prior to performance on a difficult WM task and to pre-activate a connectivity “road map” for subsequent task performance. This early formation of connectivity may be an efficient way to cope with the high processing demands during a task by decreasing the number of connections that have to be formed and allowing more resources to be allocated to the formation of other connections during the task.

Does the connectivity among the WM regions during preparation periods resembles the connectivity at rest?

Neither the preparation nor the rest periods require maintenance or manipulation of information on-line. As such, it is reasonable to propose that the preparation-related connectivity might resemble the connectivity at rest. To test this hypothesis we examined the connections that were common for task preparation and rest, but were absent during task performance. Surprisingly, we found no such connections in the Increase network. In the Decrease network, we identified two connections common for preparation and rest, but not task performance: One for the 3-back preparation (RTP-RPlanPol and LPCCa→RLOCi) and one such connection for the 1-back preparation (RTP-RPlanPol). The fact that there were more common connections for task preparation and task performance (but not rest) than for task preparation and rest (but not task performance) suggests that the connectivity during task preparation is more similar to that during task performance than to the connectivity at rest. These results also suggest that, despite the fact that the task preparation periods do not involve any active information processing, task preparation is an active state whose role is to “bridge” resting state and information processing phases in the experiment by providing timely connection and disconnection within the WM networks.

Limitations

One limitation of the current study is that the task performance blocks contain more data points than rest and preparation periods (>200 vs. >90 vs. 40 data points per condition). One might wonder whether the fewer data points in the analyses of the preparation periods might result in reduced power. However, based on a simulation study by Ramsey et al. (2011), we doubt that this is an issue for our study. First, Ramsey et al. point out that one of the most important factors for the identification and orientation of the connections is the number of subjects in the data set. They report a monotonic decrease in error rates with an increase in the number of subjects and indicate that, with as few as 10 subjects in a data set, identification and orientation of connections are accurate. Given that our data set included data from 16 subjects our calculations should be accurate. Second, the simulation in Ramsey et al. included longer sessions (close to the length of our n-back task blocks) and shorter sessions (close in length to our task preparation blocks). Given that the data from 10 (rather than one) simulated subjects were analyzed using the combined IMaGES and LOFS method, the accuracy of edge detection/discovery was very high (close to 100%) and quite similar for longer and shorter sessions. Orientation precision and recall were somewhat lower for shorter than for longer sessions, but the accuracy for shorter sessions was still very high (above 80%). This latter finding suggests the possibility that the accuracy of edge orientation in our study was higher for task performance than for the rest periods and task preparation.

The main hypotheses of our study were related to the correspondence of task preparation connectivity with the task performance connectivity and the connectivity at rest. Therefore, even though we identified several connections that were unique to a preparatory or a task performance condition (e.g., LFP→LIPS for the 3-back task condition, or Rprecun→RFP for the 1-back preparation condition), or that were present for some conditions, but not the others (e.g., RSFG→RIPS connection was observed for the 1-back task, 2-back preparation and task, and 3-back preparation conditions, but not for the 1-back preparation or 3-back task conditions), we did not discuss those connections. The presence of unique connections suggests that each task preparation or task performance condition may require using unique strategies of information processing. Unfortunately, the current design does not allow us to examine these possibilities. Our study also does not allow us to examine individual difference in connectivity within the Increase and Decrease networks. Each connectivity model was derived from multiple subjects. The model fit varies across subjects and may, or may not, be related to a given subject's behavioral performance. While it is possible that faster and/or more accurate performance is associated with formation of some specific connections, our study does not have enough power to test this possibility.

Conclusion

In summary, the present study provides novel findings about the relationship between connectivity during task preparation, task performance and rest periods by demonstrating that the connectivity among the regions within the WM network is not limited to task performance. Even though preparation periods did not require any active information maintenance or manipulation, effective connectivity during task preparation did not resemble connectivity during rest periods. Instead, two connectivity sub-networks were identified. The “default” connectivity sub-network was present in all conditions including rest. We proposed that this sub-network might serve as a core system that allows one to quickly activate cognitive, perceptual and motor systems in response to the relevant stimuli. “Preparation-related” connectivity exists during task preparation and task performance, but not at rest. It is specific for each difficulty level and likely “pre-activates” cognitive resources important for performance on each specific task. One role of such “pre-activation” may be to free neural resources during task performance and to help activate and carry over the task-related rules and intentions from task preparation to task performance. Future research using graph modeling should determine whether “default” and “preparation-related” connectivity is a phenomenon generalizable to all cognitive tasks or specific to just the n-back task.

Conflict of interest statement

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.

Acknowledgments

This work was supported by a National Institute of Mental Health training grant T32MH019983. The authors declare no competing financial interests.

References

  • Albert N. B., Robertson E. M., Mehta P., Miall R. C. (2009). Resting state networks and memory consolidation. Commun. Integr. Biol. 2, 530–532 10.4161/cib.2.6.9612 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
  • Altamura M., Goldberg T. E., Elvevåg B., Holroyd T., Carver F. W., Weinberger D. R., et al. (2010). Prefrontal cortex modulation during anticipation of working memory demands as revealed by magnetoencephalography. Int. J. Biomed. Imaging 2010:840416 10.1155/2010/840416 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
  • Anderson T. W., Darling D. A. (1952). Asymptotic theory of certain“goodness of fit” criteria based on stochastic processes. Ann. Math. Stat. 23, 193–212 10.1214/aoms/1177729437 [CrossRef] [Google Scholar]
  • Baddeley A. (2010). Working memory. Curr. Biol. 20, 136–140 10.1016/j.cub.2009.12.014 [Abstract] [CrossRef] [Google Scholar]
  • Baddeley A. D., Hitch G. J. (1974). Working memory, in Recent Advances in Learning and Motivation, Vol. 8. ed Bower G. A. (New York, NY: Academic Press; ), 47–90 [Google Scholar]
  • Beckmann C. F., Smith S. M. (2004). Probabilistic independent component analysis for functional magnetic resonance imaging. IEEE Trans. Med. Imaging 23, 137–152 10.1109/TMI.2003.822821 [Abstract] [CrossRef] [Google Scholar]
  • Botvinick M. M., Braver T. S., Barch D. M., Carter C. S., Cohen J. D. (2001). Conflict monitoring and cognitive control. Psychol. Rev. 108, 624–652 10.1037/0033-295X.108.3.624 [Abstract] [CrossRef] [Google Scholar]
  • Botvinick M. M., Cohen J. D., Carter C. S. (2004). Conflict monitoring and anterior cingulate cortex: an update. Trends Cogn. Sci. 8, 539–546 10.1016/j.tics.2004.10.003 [Abstract] [CrossRef] [Google Scholar]
  • Boukrina O., Graves W. W. (2013). Neural networks underlying contributions from semantics in reading aloud. Front. Hum. Neurosci. 7:518 10.3389/fnhum.2013.00518 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
  • Boukrina O., Hanson S. J., Hanson C. (2014). Modeling activation and effective connectivity of VWFA in same script bilinguals. Hum. Brain Mapp. 35, 2543–2560 10.1002/hbm.22348 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
  • Brass M., von Cramon D. Y. (2004). Decomposing components of task preparation with functional magnetic resonance imaging. J. Cogn. Neurosci. 16, 609–620 10.1162/089892904323057335 [Abstract] [CrossRef] [Google Scholar]
  • Braver T. S., Cohen J. D., Nystrom L. E., Jonides J., Smith E. E., Noll D. C. (1997). A parametric study of frontal cortex involvement in human working memory. Neuroimage 5, 49–62 10.1006/nimg.1996.0247 [Abstract] [CrossRef] [Google Scholar]
  • Buckner R. L., Andrews-Hanna J. R., Schacter D. L. (2008). The brain's default network: anatomy, function, and relevance to disease. Ann. N.Y. Acad. Sci. 1124, 1–38 10.1196/annals.1440.011 [Abstract] [CrossRef] [Google Scholar]
  • Corbetta M. (1998). Frontoparietal cortical networks for directing attention and the eye to visual locations: identical, independent, or overlapping neural systems? Proc. Natl. Acad. Sci. U.S.A. 95, 831–838 [Europe PMC free article] [Abstract] [Google Scholar]
  • Costa A., Oliveri M., Barban F., Bonnì S., Koch G., Caltagirone C., et al. (2013). The right frontopolar cortex is involved in visual-spatial prospective memory. PLoS ONE 8:e56039 10.1371/journal.pone.0056039 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
  • Cowan N. (2001). The magical number 4 in short-term memory: a reconsideration of mental storage capacity. Behav. Brain. Sci. 24, 87–114 10.1017/S0140525X01003922 [Abstract] [CrossRef] [Google Scholar]
  • Curtis C. E., D'Esposito M. (2003). Persistent activity in the prefrontal cortex during working memory. Trends Cogn. Sci. 7, 415–423 10.1016/S1364-6613(03)00197-9 [Abstract] [CrossRef] [Google Scholar]
  • Dima D., Jogia J., Frangou S. (2014). Dynamic causal modeling of load-dependent modulation of effective connectivity within the verbal working memory network. Hum. Brain Mapp. 35, 3025–3035 10.1002/hbm.22382 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
  • du Boisgueheneuc F., Levy R., Volle E., Seassau M., Duffau H., Kinkingnehun S., et al. (2006). Functions of the left superior frontal gyrus in humans: a lesion study. Brain 129, 3315–3328 10.1093/brain/awl244 [Abstract] [CrossRef] [Google Scholar]
  • Esposito F., Bertolino A., Scarabino T., Latorre V., Blasi G., Popolizio T., et al. (2006). Independent component model of the default-mode brain function: assessing the impact of active thinking. Brain Res. Bull. 70, 263–269 10.1016/j.brainresbull.2006.06.012 [Abstract] [CrossRef] [Google Scholar]
  • Greicius M. D., Krasnow B., Reiss A. L., Menon V. (2003). Functional connectivity in the resting brain: a network analysis of the default mode hypothesis. Proc. Natl. Acad. Sci. U.S.A. 100, 253–258 10.1073/pnas.0135058100 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
  • Hanson C., Hanson S. J., Ramsey J., Glymour C. (2013). Atypical effective connectivity of social brain networks in individuals with autism. Brain Connect. 3, 578–589 10.1089/brain.2013.0161 [Abstract] [CrossRef] [Google Scholar]
  • Hanson S. J., Hanson C., Halchenko Y., Matsuka T., Zaimi A. (2007). Bottom-up and top-down brain functional connectivity underlying comprehension of everyday visual action. Brain Struct. Funct. 212, 231–244 10.1007/s00429-007-0160-2 [Abstract] [CrossRef] [Google Scholar]
  • Hasson U., Nusbaum H. C., Small S. L. (2009). Task-dependent organization of brain regions active during rest. Proc. Natl. Acad. Sci. U.S.A. 106, 10841–10846 10.1073/pnas.0903253106 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
  • Honey G. D., Fu C. H. Y., Kim J., Brammer M. J., Croudace T. J., Suckling J., et al. (2002). Effects of verbal working memory load on corticocortical connectivity modeled by path analysis of functional magnetic resonance imaging data. Neuroimage 17, 573–582 10.1006/nimg.2002.1193 [Abstract] [CrossRef] [Google Scholar]
  • Jenkinson M. (2003). A fast, automated, n-dimensional phase unwrapping algorithm. Magn. Reson. Med. 49, 193–197 10.1002/mrm.10354 [Abstract] [CrossRef] [Google Scholar]
  • Jenkinson M., Bannister P., Brady M., Smith S. (2002). Improved optimisation for the robust and accurate linear registration and motion correction of brain images. Neuroimage 17, 825–841 10.1006/nimg.2002.1132 [Abstract] [CrossRef] [Google Scholar]
  • Jenkinson M., Smith S. (2001). A global optimisation method for robust affine registration of brain images. Med. Image Anal. 5, 143–156 10.1016/S1361-8415(01)00036-6 [Abstract] [CrossRef] [Google Scholar]
  • Kelly R. E., Jr., Alexopoulos G. S., Wang Z., Gunning F. M., Murphy C. F., Morimoto S. S., et al. (2010). Visual inspection of independent components: defining a procedure for artifact removal from fMRI data. J. Neurosci. Methods 189, 233–245 10.1016/j.jneumeth.2010.03.028 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
  • Kim J. S., Jung W. H., Kang D.-H., Park J.-Y., Jang J. H., Choi J.-S., et al. (2012). Changes in effective connectivity according to working memory load: an FMRI study of face and location working memory tasks. Psychiatry Invest. 9, 283–292 10.4306/pi.2012.9.3.283 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
  • Ma L., Steinberg J. L., Hasan K. M., Narayana P. A., Kramer L. A., Moeller F. G. (2012). Working memory load modulation of parieto-frontal connections: evidence from dynamic causal modeling. Hum. Brain Mapp. 33, 1850–1867 10.1002/hbm.21329 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
  • Manelis A., Reder L. M. (2013). He Who is well prepared has half won the battle: an fMRI Study of task preparation. Cereb. Cortex. [Epub ahead of print]. 10.1093/cercor/bht262 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
  • Mayer J. S., Roebroeck A., Maurer K., Linden D. E. J. (2010). Specialization in the default mode: task-induced brain deactivations dissociate between visual working memory and attention. Hum. Brain Mapp. 31, 126–139 10.1002/hbm.20850 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
  • McKiernan K. A., D'Angelo B. R., Kaufman J. N., Binder J. R. (2006). Interrupting the “stream of consciousness:” an fMRI investigation. Neuroimage 29, 1185–1191 10.1016/j.neuroimage.2005.09.030 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
  • McKiernan K. A., Kaufman J. N., Kucera-Thompson J., Binder J. R. (2003). A parametric manipulation of factors affecting task-induced deactivation in functional neuroimaging. J. Cogn. Neurosci. 15, 394–408 10.1162/089892903321593117 [Abstract] [CrossRef] [Google Scholar]
  • Mishkin M., Manning F. J. (1978). Non-spatial memory after selective prefrontal lesions in monkeys. Brain Res. 143, 313–323 [Abstract] [Google Scholar]
  • Müller N. G., Knight R. T. (2006). The functional neuroanatomy of working memory: contributions of human brain lesion studies. Neuroscience 139, 51–58 10.1016/j.neuroscience.2005.09.018 [Abstract] [CrossRef] [Google Scholar]
  • Mumford J. A., Ramsey J. D. (2014). Bayesian networks for fMRI: a primer. Neuroimage 86, 573–582 10.1016/j.neuroimage.2013.10.020 [Abstract] [CrossRef] [Google Scholar]
  • Neufang S., Akhrif A., Riedl V., Förstl H., Kurz A., Zimmer C., et al. (2011). Disconnection of frontal and parietal areas contributes to impaired attention in very early Alzheimer's disease. J. Alzheimers Dis. 25, 309–321 10.3233/JAD-2011-102154 [Abstract] [CrossRef] [Google Scholar]
  • Nystrom L. E., Braver T. S., Sabb F. W., Delgado M. R., Noll D. C., Cohen J. D. (2000). Working memory for letters, shapes, and locations: fMRI evidence against stimulus-based regional organization in human prefrontal cortex. Neuroimage 11, 424–446 10.1006/nimg.2000.0572 [Abstract] [CrossRef] [Google Scholar]
  • Owen A. M., McMillan K. M., Laird A. R., Bullmore E. (2005). N-back working memory paradigm: a meta-analysis of normative functional neuroimaging studies. Hum. Brain Mapp. 25, 46–59 10.1002/hbm.20131 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
  • Petrides M., Milner B. (1982). Deficits on subject-ordered tasks after frontal- and temporal-lobe lesions in man. Neuropsychologia 20, 249–262 [Abstract] [Google Scholar]
  • Raichle M. E., MacLeod A. M., Snyder A. Z., Powers W. J., Gusnard D. A., Shulman G. L. (2001). A default mode of brain function. Proc. Natl. Acad. Sci. U.S.A. 98, 676–682 10.1073/pnas.98.2.676 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
  • Raichle M. E., Snyder A. Z. (2007). A default mode of brain function: a brief history of an evolving idea. Neuroimage 37, 1083–1090 10.1016/j.neuroimage.2007.02.041 [Abstract] [CrossRef] [Google Scholar]
  • Ramsey J. D., Hanson S. J., Hanson C., Halchenko Y. O., Poldrack R. A., Glymour C. (2010). Six problems for causal inference from fMRI. Neuroimage 49, 1545–1558 10.1016/j.neuroimage.2009.08.065 [Abstract] [CrossRef] [Google Scholar]
  • Ramsey J. D., Hanson S. J., Glymour C. (2011). Multi-subject search correctly identifies causal connections and most causal directions in the DCM models of the Smith et al. simulation study. Neuroimage 58, 838–848 10.1016/j.neuroimage.2011.06.068 [Abstract] [CrossRef] [Google Scholar]
  • Ramsey J. D., Sanchez-Romero R., Glymour C. (2014). Non-Gaussian methods and high-pass filters in the estimation of effective connections. Neuroimage 84, 986–1006 10.1016/j.neuroimage.2013.09.062 [Abstract] [CrossRef] [Google Scholar]
  • Rottschy C., Langner R., Dogan I., Reetz K., Laird A. R., Schulz J. B., et al. (2012). Modelling neural correlates of working memory: a coordinate-based meta-analysis. Neuroimage 60, 830–846 10.1016/j.neuroimage.2011.11.050 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
  • Sambataro F., Murty V. P., Callicott J. H., Tan H.-Y., Das S., Weinberger D. R., et al. (2010). Age-related alterations in default mode network: impact on working memory performance. Neurobiol. Aging 31, 839–852 10.1016/j.neurobiolaging.2008.05.022 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
  • Smith S. (2002). Fast robust automated brain extraction. Hum. Brain Mapp. 17, 143–155 10.1002/hbm.10062 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
  • Smith S. M., Miller K. L., Salimi- Khorshidi G., Webster M., Beckmann C., Nichols T., et al. (2011). Network modeling methods for fMRI. Neuroimage 54, 875–891 10.1016/j.neuroimage.2010.08.063 [Abstract] [CrossRef] [Google Scholar]
  • Sohn M.-H., Albert M. V., Jung K., Carter C. S., Anderson J. R. (2007). Anticipation of conflict monitoring in the anterior cingulate cortex and the prefrontal cortex. Proc. Natl. Acad. Sci. U.S.A. 104, 10330–10334 10.1073/pnas.0703225104 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
  • Tambini A., Ketz N., Davachi L. (2010). Enhanced brain correlations during rest are related to memory for recent experiences. Neuron 65, 280–290 10.1016/j.neuron.2010.01.001 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
  • Tohka J., Foerde K., Aron A. R., Tom S. M., Toga A. W., Poldrack R. A. (2008). Automatic independent component labeling for artifact removal in fMRI. Neuroimage 39, 1227–1245 10.1016/j.neuroimage.2007.10.013 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
  • Vapnik V. (1995). The Nature of Statistical Learning Theory. New York, NY: Springer [Google Scholar]
  • Vincent J. L., Patel G. H., Fox M. D., Snyder A. Z., Baker J. T., Van Essen D. C., et al. (2007). Intrinsic functional architecture in the anaesthetized monkey brain. Nature 447, 83–86 10.1038/nature05758 [Abstract] [CrossRef] [Google Scholar]
  • Volle E., Gonen-Yaacovi G., Costello A. d. L., Gilbert S. J., Burgess P. W. (2011). The role of rostral prefrontal cortex in prospective memory: a voxel-based lesion study. Neuropsychologia 49, 2185–2198 10.1016/j.neuropsychologia.2011.02.045 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
  • Wang L., Liu X., Guise K. G., Knight R. T., Ghajar J., Fan J. (2010). Effective connectivity of the fronto-parietal network during attentional control. J. Cogn. Neurosci. 22, 543–553 10.1162/jocn.2009.21210 [Abstract] [CrossRef] [Google Scholar]
  • Wendelken C., Bunge S. A., Carter C. S. (2008). Maintaining structured information: an investigation into functions of parietal and lateral prefrontal cortices. Neuropsychologia 46, 665–678 10.1016/j.neuropsychologia.2007.09.015 [Abstract] [CrossRef] [Google Scholar]
  • Yakushev I., Chételat G., Fischer F. U., Landeau B., Bastin C., Scheurich A., et al. (2013). Metabolic and structural connectivity within the default mode network relates to working memory performance in young healthy adults. Neuroimage 79, 184–190 10.1016/j.neuroimage.2013.04.069 [Abstract] [CrossRef] [Google Scholar]
  • Yehene E., Meiran N., Soroker N. (2008). Basal ganglia play a unique role in task switching within the frontal-subcortical circuits: evidence from patients with focal lesions. J. Cogn. Neurosci. 20, 1079–1093 10.1162/jocn.2008.20077 [Abstract] [CrossRef] [Google Scholar]

Articles from Frontiers in Human Neuroscience are provided here courtesy of Frontiers Media SA

Citations & impact 


Impact metrics

Jump to Citations

Citations of article over time

Alternative metrics

Altmetric item for https://www.altmetric.com/details/2513795
Altmetric
Discover the attention surrounding your research
https://www.altmetric.com/details/2513795

Article citations


Go to all (9) article citations

Similar Articles 


To arrive at the top five similar articles we use a word-weighted algorithm to compare words from the Title and Abstract of each citation.


Funding 


Funders who supported this work.

NIMH NIH HHS (1)